RADseq as a valuable tool for plants with large genomes—A case study in cycads

Full genome sequencing of organisms with large and complex genomes is intractable and cost ineffective under most research budgets. Cycads (Cycadales) represent one of the oldest lineages of the extant seed plants and, partly due to their age, have incredibly large genomes up to ~60 Gbp. Restriction site‐associated DNA sequencing (RADseq) offers an approach to find genome‐wide informative markers and has proven to be effective with both model and nonmodel organisms. We tested the application of RADseq using ezRAD across all 10 genera of the Cycadales including an example data set of Cycas calcicola representing 72 samples from natural populations. Using previously available plastid and mitochondrial genomes as references, reads were mapped recovering plastid and mitochondrial genome regions and nuclear markers for all of the genera. De novo assembly generated up to 138,407 high‐depth clusters and up to 1,705 phylogenetically informative loci for the genera, and 4,421 loci for the example assembly of C. calcicola. The number of loci recovered by de novo assembly was lower than previous RADseq studies, yet still sufficient for downstream analysis. However, the number of markers could be increased by relaxing our assembly parameters, especially for the C. calcicola data set. Our results demonstrate the successful application of RADseq across the Cycadales to generate a large number of markers for all genomic compartments, despite the large number of plastids present in a typical plant cell. Our modified protocol was adapted to be applied to cycads and other organisms with large genomes to yield many informative genome‐wide markers.

• Publication year

  2019

• Venue

  Molecular Ecology Resources

• Publication date

  2019-09-30

• Fields of study

  Biology, Medicine, Environmental Science

• Identifiers
  DOI 10.1111/1755-0998.13085PMID 31484214
• External record

  Open on Semantic Scholar

• Source metadata

  Semantic Scholar, PubMed

• No claims are published for this paper.

• No concepts are published for this paper.

• Oryza sativa

  2021influential reference
• ipyrad: Interactive assembly and analysis of RADseq datasets

  2020cited by this paper
• Population genomic analysis of mango (Mangifera indica) suggests a complex history of domestication.

  2019influential reference
• Tangled banks: A landscape genomic evaluation of Wallace's Riverine barrier hypothesis for three Amazon plant species

  2019influential reference
• Genome Size Diversity and Its Impact on the Evolution of Land Plants

  2018cited by this paper
• PSMC (pairwise sequentially Markovian coalescent) analysis of RAD (restriction site associated DNA) sequencing data

  2017cited by this paper
• Evidence for an ancient whole genome duplication in the cycad lineage

  2017cited by this paper
• Unbroken: RADseq remains a powerful tool for understanding the genetics of adaptation in natural populations

  2017cited by this paper
• Phylogenomic inferences from reference-mapped and de novo assembled short-read sequence data using RADseq sequencing of California white oaks (Quercus section Quercus).

  2017cited by this paper
• Genomic diversity guides conservation strategies among rare terrestrial orchid species when taxonomy remains uncertain

  2017cited by this paper
• Comparative Analysis of the Complete Chloroplast Genomes of Five Quercus Species

  2016cited by this paper
• Misconceptions on Missing Data in RAD‐seq Phylogenetics with a Deep‐scale Example from Flowering Plants

  2016influential reference
• The complete plastid genome sequence of Abies koreana (Pinaceae: Abietoideae)

  2016cited by this paper
• Ginkgo and Welwitschia Mitogenomes Reveal Extreme Contrasts in Gymnosperm Mitochondrial Evolution.

  2016influential reference
• Extensive Natural Variation in Arabidopsis Seed Mucilage Structure

  2016cited by this paper
• Draft genome of the living fossil Ginkgo biloba

  2016cited by this paper
• Genome sizes for 71 species of Zamia (Cycadales: Zamiaceae) correspond with three different biogeographic regions

  2016cited by this paper
• Utilizing RADseq data for phylogenetic analysis of challenging taxonomic groups: A case study in Carex sect. Racemosae.

  2016influential reference
• Characterization of the complete chloroplast genome sequence of Primula veris (Ericales: Primulaceae)

  2016cited by this paper
• Harnessing the power of RADseq for ecological and evolutionary genomics

  2016influential reference
• Adaptive divergence in the monkey flower Mimulus guttatus is maintained by a chromosomal inversion

  2015influential reference
• On the Robustness to Gene Tree Estimation Error (or lack thereof) of Coalescent-Based Species Tree Methods.

  2015cited by this paper
• Processes Driving the Adaptive Radiation of a Tropical Tree ( Diospyros , Ebenaceae) in New Caledonia, a Biodiversity Hotspot

  2015influential reference
• Whole genome duplication in coast redwood (Sequoia sempervirens) and its implications for explaining the rarity of polyploidy in conifers

  2015cited by this paper
• ASTRAL-II: coalescent-based species tree estimation with many hundreds of taxa and thousands of genes

  2015cited by this paper
• Phylogeny of Morella rubra and Its Relatives (Myricaceae) and Genetic Resources of Chinese Bayberry Using RAD Sequencing

  2015influential reference
• Historical introgression among the American live oaks and the comparative nature of tests for introgression

  2015cited by this paper
• Genomic exploration and molecular marker development in a large and complex conifer genome using RADseq and mRNAseq

  2015cited by this paper
• Thousands of RAD-seq Loci Fully Resolve the Phylogeny of the Highly Disjunct Arctic-Alpine Genus Diapensia (Diapensiaceae)

  2015cited by this paper
• Does complete plastid genome sequencing improve species discrimination and phylogenetic resolution in Araucaria?

  2015cited by this paper
• Evolutionary Stasis in Cycad Plastomes and the First Case of Plastome GC-Biased Gene Conversion

  2015cited by this paper
• Restriction site‐associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference

  2015cited by this paper
• Trimmomatic: a flexible trimmer for Illumina sequence data

  2014cited by this paper
• An improved 2b-RAD approach (I2b-RAD) offering genotyping tested by a rice (Oryza sativa L.) F2 population

  2014cited by this paper
• ASTRAL: genome-scale coalescent-based species tree estimation

  2014cited by this paper
• Disk covering methods improve phylogenomic analyses

  2014cited by this paper
• Sequencing and Assembly of the 22-Gb Loblolly Pine Genome

  2014cited by this paper
• Unique Features of the Loblolly Pine (Pinus taeda L.) Megagenome Revealed Through Sequence Annotation

  2014cited by this paper
• Efficient de novo assembly of large and complex genomes by massively parallel sequencing of Fosmid pools

  2014cited by this paper
• Trade‐offs and utility of alternative RADseq methods: Reply to Puritz et al.

  2014cited by this paper
• Quartet Inference from SNP Data Under the Coalescent Model

  2014cited by this paper
• Conserved genetic regions across angiosperms as tools to develop single‐copy nuclear markers in gymnosperms: an example using cycads

  2014cited by this paper
• Recent novel approaches for population genomics data analysis

  2014cited by this paper
• Can a Botanic Garden Cycad Collection Capture the Genetic Diversity in a Wild Population?

  2014cited by this paper
• Phylogenetic analysis at deep timescales: unreliable gene trees, bypassed hidden support, and the coalescence/concatalescence conundrum.

  2014cited by this paper
• Genomic evidence for the parallel evolution of coastal forms in the Senecio lautus complex

  2013influential reference
• Genome sequence of dwarf birch (Betula nana) and cross‐species RAD markers

  2013cited by this paper
• Read and assembly metrics inconsequential for clinical utility of whole-genome sequencing in mapping outbreaks

  2013cited by this paper
• RADseq underestimates diversity and introduces genealogical biases due to nonrandom haplotype sampling

  2013cited by this paper
• Genotyping‐by‐sequencing in ecological and conservation genomics

  2013cited by this paper
• SNP discovery using Paired‐End RAD‐tag sequencing on pooled genomic DNA of Sisymbrium austriacum (Brassicaceae)

  2013influential reference
• ezRAD: a simplified method for genomic genotyping in non-model organisms

  2013influential reference
• Inferring Phylogeny and Introgression using RADseq Data: An Example from Flowering Plants (Pedicularis: Orobanchaceae)

  2013cited by this paper
• Phylogeny of the cycads based on multiple single-copy nuclear genes: congruence of concatenated parsimony, likelihood and species tree inference methods.

  2013cited by this paper
• PyRAD: assembly of de novo RADseq loci for phylogenetic analyses

  2013cited by this paper
• Zamia (Cycadales: Zamiaceae) on Puerto Rico: asymmetric genetic differentiation and the hypothesis of multiple introductions.

  2012cited by this paper
• Conifer genome sizes of 172 species, covering 64 of 67 genera, range from 8 to 72 picogram

  2012cited by this paper
• Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species

  2012cited by this paper
• Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data

  2012cited by this paper
• Recent Synchronous Radiation of a Living Fossil

  2011cited by this paper
• Complete plastid genome sequences of three Rosids (Castanea, Prunus, Theobroma): evidence for at least two independent transfers of rpl22 to the nucleus.

  2011cited by this paper
• RADSeq : next-generation population genetics 417

  2011cited by this paper
• Quality control and preprocessing of metagenomic datasets

  2011influential reference
• Genome sizes for all genera of Cycadales.

  2011influential reference
• When North and South don’t mix: genetic connectivity of a recently endangered oceanic cycad, Cycas micronesica, in Guam using EST‐microsatellites

  2010cited by this paper
• The largest eukaryotic genome of them all

  2010cited by this paper
• RADSeq: next-generation population genetics.

  2010cited by this paper
• IUCN Red List of Threatened Species

  2010cited by this paper
• A maximum pseudo-likelihood approach for estimating species trees under the coalescent model

  2010cited by this paper
• Coalescent methods for estimating phylogenetic trees.

  2009cited by this paper
• How well do we understand the overall backbone of cycad phylogeny? New insights from a large, multigene plastid data set.

  2008cited by this paper
• Chloroplast genome (cpDNA) of Cycas taitungensis and 56 cp protein-coding genes of Gnetum parvifolium: insights into cpDNA evolution and phylogeny of extant seed plants.

  2007influential reference
• Genetic Variation in the Endangered Endemic Species Cycas fairylakea (Cycadaceae) in China and Implications for Conservation

  2006cited by this paper
• Mapping and sequencing complex genomes: let's get physical!

  2004cited by this paper
• PAUP* [Phylogenetic Analysis Using Parsimony (and Other Methods)]

  2004cited by this paper
• ISSR variation in the endemic and endangered plant Cycas guizhouensis (Cycadaceae).

  2004cited by this paper
• Responses of cycads with different life histories to the impact of plant collecting: simulation models to determine important life history stages and population recovery times

  2003cited by this paper
• Cycads: Status Survey and Conservation Action Plan.

  2003cited by this paper
• PAUP*: Phylogenetic analysis using parsimony (*and other methods), Version 4.0b10

  2002cited by this paper
• PAUP* Phylogenetic analysis using parsimony (*and other methods) Version 4

  2000cited by this paper
• Monitoring endemic plant extinction in Veracruz, Mexico

  1998cited by this paper
• Chloroplast DNA of black pine retains a residual inverted repeat lacking rRNA genes: nucleotide sequences of trnQ, trnK, psbA, trnI and trnH and the absence of rps16

  1992cited by this paper

• A reduced-representation sequencing method for identifying sex-specific molecular markers in plants with large genomes

  2025cites this paper
• Integrative species delimitation of the Cycas sexseminifera complex (Cycadaceae) utilizing genomic and morphological data

  2025cites this paper
• Phylogenomic incongruence in Castanopsis (Fagaceae) is primarily caused by gene flow rather than incomplete lineage sorting.

  2025cites this paper
• A genomic perspective on adaptation to heterogeneous environments for Cycas segmentifida (Cycadaceae), an endemic cycad in southwest China.

  2025cites this paper
• Next‐generation sequencing sheds new light on a major lineage of Acanthaceae, clarifying previously uncertain relationships in Barlerieae

  2025cites this paper
• Evidence of a Historical Metapopulation: Conservation Genomics and Management Strategies of Cycas candida (Cycadaceae)

  2024influential citation
• Approach of Genetic Diversity of Lippia alba (Mill) and Petiveria alliacea L.: Medicinal Plants of Colombia

  2024cites this paper
• Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes).

  2023cites this paper
• Optimising ddRAD sequencing for population genomic studies with ddgRADer

  2023cites this paper
• Species relationships and phylogenetic diversity of the African genus Encephalartos Lehm. (Zamiaceae)

  2023cites this paper
• De Novo SNP Discovery and Genotyping of Masson Pine (Pinus massoniana Lamb.) via Genotyping-by-Sequencing

  2023cites this paper
• Smart breeding driven by advances in sequencing technology

  2023cites this paper
• Exploring the impact of read clustering thresholds on RADseq-based systematics: an empirical example from European amphibians

  2023cites this paper
• Optimizing ddRAD sequencing for population genomic studies with ddgRADer

  2023cites this paper
• Conservation genomics of Dioon holmgrenii (Zamiaceae) reveals a history of range expansion, fragmentation, and isolation of populations

  2023cites this paper
• Genomics and conservation: Guidance from training to analyses and applications

  2023cites this paper
• A Tribute to Nathalie Seenamah Nagalingum (1975–2022)

  2023cites this paper
• The RadOrgMiner pipeline: Automated genotyping of organellar loci from RADseq data

  2022cites this paper
• Conservation genomics of an Australian cycad Cycas calcicola, and the Absence of Key Genotypes in Botanic Gardens

  2022cites this paper
• Target sequence capture of Barnadesioideae (Compositae) demonstrates the utility of low coverage loci in phylogenomic analyses.

  2022cites this paper
• A RADseq Phylogeny of Barleria (Acanthaceae) Resolves Fine-Scale Relationships.

  2022cites this paper
• A turn in species conservation for hairpin banksias: demonstration of oversplitting leads to better management of diversity

  2022cites this paper
• Genome-Wide Evidence for Complex Hybridization and Demographic History in a Group of Cycas From China

  2021cites this paper
• The Endangered Species Dioon edule in the Sierra Madre Oriental in San Luis Potosí: Demography and Genetic Diversity

  2021cites this paper
• DNA degradation in fish: Practical solutions and guidelines to improve DNA preservation for genomic research

  2020cites this paper
• Effective double‐digest RAD sequencing and genotyping despite large genome size

  2020cites this paper
• Retracing origins of exceptional cycads in botanical collections to increase conservation value

  2020cites this paper
• ipyrad: Interactive assembly and analysis of RADseq datasets

  2020cites this paper