Hybridization ddRAD‐sequencing for population genomics of nonmodel plants using highly degraded historical specimen DNA

• Published in: Molecular ecology resources 2020-09, Vol.20 (5), p.1228-1247
• Main Authors: Lang, Patricia L. M., Weiß, Clemens L., Kersten, Sonja, Latorre, Sergio M., Nagel, Sarah, Nickel, Birgit, Meyer, Matthias, Burbano, Hernán A.
• Format: Article
• Language: English
• Subjects: ancient DNA; Anthropogenic factors; Baits; capture; Deoxyribonucleic acid; DNA; DNA damage; DNA fragmentation; DNA sequencing; Evolution; Genetic diversity; Genomes; herbarium; Hybridization; hybridization double‐digest RADseq; Libraries; nonmodel species; Nucleotide sequence; Representations; Species
• ISSN: 1755-098X; 1755-0998
• DOI: 10.1111/1755-0998.13168

Species’ responses at the genetic level are key to understanding the long‐term consequences of anthropogenic global change. Herbaria document such responses, and, with contemporary sampling, provide high‐resolution time‐series of plant evolutionary change. Characterizing genetic diversity is straightforward for model species with small genomes and a reference sequence. For nonmodel species—with small or large genomes—diversity is traditionally assessed using restriction‐enzyme‐based sequencing. However, age‐related DNA damage and fragmentation preclude the use of this approach for ancient herbarium DNA. Here, we combine reduced‐representation sequencing and hybridization‐capture to overcome this challenge and efficiently compare contemporary and historical specimens. Specifically, we describe how homemade DNA baits can be produced from reduced‐representation libraries of fresh samples, and used to efficiently enrich historical libraries for the same fraction of the genome to produce compatible sets of sequence data from both types of material. Applying this approach to both Arabidopsis thaliana and the nonmodel plant Cardamine bulbifera, we discovered polymorphisms de novo in an unbiased, reference‐free manner. We show that the recovered genetic variation recapitulates known genetic diversity in A. thaliana, and recovers geographical origin in both species and over time, independent of bait diversity. Hence, our method enables fast, cost‐efficient, large‐scale integration of contemporary and historical specimens for assessment of genome‐wide genetic trends over time, independent of genome size and presence of a reference genome.