Genomic treasure troves: complete genome sequencing of herbarium and insect museum specimens

Ostracod (Ostracoda, Crustacea) genomics - Promises and challenges

Ostracods are well-suited model organisms for evolutionary research. Classic genetic techniques have mostly been used for phylogenetic studies on Ostracoda and were somewhat affected by the lack of large numbers of suitable markers. Genomic methods with their huge potential have so far rarely been applied to this group of crustaceans. We provide relevant examples of genomic studies on other organisms to propose future avenues of genomic ostracod research. At the same time, we suggest solutions to the potential problems in ostracods that the application of genomic techniques might present.

Successful Recovery of Nuclear Protein-Coding Genes from Small Insects in Museums Using Illumina Sequencing

In this paper we explore high-throughput Illumina sequencing of nuclear protein-coding, ribosomal, and mitochondrial genes in small, dried insects stored in natural history collections. We sequenced one tenebrionid beetle and 12 carabid beetles ranging in size from 3.7 to 9.7 mm in length that have been stored in various museums for 4 to 84 years. Although we chose a number of old, small specimens for which we expected low sequence recovery, we successfully recovered at least some low-copy nuclear protein-coding genes from all specimens. For example, in one 56-year-old beetle, 4.4 mm in length, our de novo assembly recovered about 63% of approximately 41,900 nucleotides in a target suite of 67 nuclear protein-coding gene fragments, and 70% using a reference-based assembly. Even in the least successfully sequenced carabid specimen, reference-based assembly yielded fragments that were at least 50% of the target length for 34 of 67 nuclear protein-coding gene fragments. Exploration of alternative references for reference-based assembly revealed few signs of bias created by the reference. For all specimens we recovered almost complete copies of ribosomal and mitochondrial genes. We verified the general accuracy of the sequences through comparisons with sequences obtained from PCR and Sanger sequencing, including of conspecific, fresh specimens, and through phylogenetic analysis that tested the placement of sequences in predicted regions. A few possible inaccuracies in the sequences were detected, but these rarely affected the phylogenetic placement of the samples. Although our sample sizes are low, an exploratory regression study suggests that the dominant factor in predicting success at recovering nuclear protein-coding genes is a high number of Illumina reads, with success at PCR of COI and killing by immersion in ethanol being secondary factors; in analyses of only high-read samples, the primary significant explanatory variable was body length, with small beetles being more successfully sequenced.

RAD-seq data point to a northern origin of the arctic-alpine genus Cassiope (Ericaceae)

Many arctic-alpine plants display a highly disjunct distribution between the Arctic/Boreal regions and the southern Asian mountains. Two main hypotheses have been proposed to explain the origin of this biogeographic pattern: (1) south-to-north migration in the late Pliocene/early Pleistocene, and (2) north-to-south migration during the Miocene. The genus Cassiope is disjunctly distributed between the Arctic/Boreal regions and the Himalayan-Hengduan Mountains (HHM) and was selected to test these hypotheses. We constructed a fossil-calibrated phylogeny of Ericaceae using two plastid regions to estimate the crown group age of Cassiope, and used sequence data from thousands of loci produced by restriction site associated DNA sequencing (RAD-seq) to reconstruct the phylogeny of Cassiope. We also performed Bayesian divergence time analysis and biogeographic analysis. The Cassiope crown group was estimated to have originated in the Miocene, which predates the onset of Northern hemisphere glaciation. All HHM species formed a clade together with one eastern Siberian species, and this clade was sister to all other Arctic/Boreal species. This topology implies a northern origin of Cassiope, which is confirmed by our biogeographic analysis. Our results thus suggest that the ancient north-to-south migration hypothesis is most consistent with the origin of Cassiope.

Sequence capture of ultraconserved elements from bird museum specimens

New DNA sequencing technologies are allowing researchers to explore the genomes of the millions of natural history specimens collected prior to the molecular era. Yet, we know little about how well specific next-generation sequencing (NGS) techniques work with the degraded DNA typically extracted from museum specimens. Here, we use one type of NGS approach, sequence capture of ultraconserved elements (UCEs), to collect data from bird museum specimens as old as 120 years. We targeted 5060 UCE loci in 27 western scrub-jays (Aphelocoma californica) representing three evolutionary lineages that could be species, and we collected an average of 3749 UCE loci containing 4460 single nucleotide polymorphisms (SNPs). Despite older specimens producing fewer and shorter loci in general, we collected thousands of markers from even the oldest specimens. More sequencing reads per individual helped to boost the number of UCE loci we recovered from older specimens, but more sequencing was not as successful at increasing the length of loci. We detected contamination in some samples and determined that contamination was more prevalent in older samples that were subject to less sequencing. For the phylogeny generated from concatenated UCE loci, contamination led to incorrect placement of some individuals. In contrast, a species tree constructed from SNPs called within UCE loci correctly placed individuals into three monophyletic groups, perhaps because of the stricter analytical procedures used for SNP calling. This study and other recent studies on the genomics of museum specimens have profound implications for natural history collections, where millions of older specimens should now be considered genomic resources.

Using targeted enrichment of nuclear genes to increase phylogenetic resolution in the neotropical rain forest genus Inga (Leguminosae: Mimosoideae)

Evolutionary radiations are prominent and pervasive across many plant lineages in diverse geographical and ecological settings; in neotropical rainforests there is growing evidence suggesting that a significant fraction of species richness is the result of recent radiations. Understanding the evolutionary trajectories and mechanisms underlying these radiations demands much greater phylogenetic resolution than is currently available for these groups. The neotropical tree genus Inga (Leguminosae) is a good example, with ~300 extant species and a crown age of 2-10 MY, yet over 6 kb of plastid and nuclear DNA sequence data gives only poor phylogenetic resolution among species. Here we explore the use of larger-scale nuclear gene data obtained though targeted enrichment to increase phylogenetic resolution within Inga. Transcriptome data from three Inga species were used to select 264 nuclear loci for targeted enrichment and sequencing. Following quality control to remove probable paralogs from these sequence data, the final dataset comprised 259,313 bases from 194 loci for 24 accessions representing 22 Inga species and an outgroup (Zygia). Bayesian phylogenies reconstructed using either all loci concatenated or a gene-tree/species-tree approach yielded highly resolved phylogenies. We used coalescent approaches to show that the same targeted enrichment data also have significant power to discriminate among alternative within-species population histories within the widespread species I. umbellifera. In either application, targeted enrichment simplifies the informatics challenge of identifying orthologous loci associated with de novo genome sequencing. We conclude that targeted enrichment provides the large volumes of phylogenetically-informative sequence data required to resolve relationships within recent plant species radiations, both at the species level and for within-species phylogeographic studies.

Ancient population genomics and the study of evolution

Recently, the study of ancient DNA (aDNA) has been greatly enhanced by the development of second-generation DNA sequencing technologies and targeted enrichment strategies. These developments have allowed the recovery of several complete ancient genomes, a result that would have been considered virtually impossible only a decade ago. Prior to these developments, aDNA research was largely focused on the recovery of short DNA sequences and their use in the study of phylogenetic relationships, molecular rates, species identification and population structure. However, it is now possible to sequence a large number of modern and ancient complete genomes from a single species and thereby study the genomic patterns of evolutionary change over time. Such a study would herald the beginnings of ancient population genomics and its use in the study of evolution. Species that are amenable to such large-scale studies warrant increased research effort. We report here progress on a population genomic study of the Adélie penguin (Pygoscelis adeliae). This species is ideally suited to ancient population genomic research because both modern and ancient samples are abundant in the permafrost conditions of Antarctica. This species will enable us to directly address many of the fundamental questions in ecology and evolution.

The Highly Reduced Plastome of Mycoheterotrophic Sciaphila (Triuridaceae) Is Colinear with Its Green Relatives and Is under Strong Purifying Selection

The enigmatic monocot family Triuridaceae provides a potentially useful model system for studying the effects of an ancient loss of photosynthesis on the plant plastid genome, as all of its members are mycoheterotrophic and achlorophyllous. However, few studies have placed the family in a comparative context, and its phylogenetic placement is only partly resolved. It was also unclear whether any taxa in this family have retained a plastid genome. Here, we used genome survey sequencing to retrieve plastid genome data for Sciaphila densiflora (Triuridaceae) and ten autotrophic relatives in the orders Dioscoreales and Pandanales. We recovered a highly reduced plastome for Sciaphila that is nearly colinear with Carludovica palmata, a photosynthetic relative that belongs to its sister group in Pandanales, Cyclanthaceae-Pandanaceae. This phylogenetic placement is well supported and robust to a broad range of analytical assumptions in maximum-likelihood inference, and is congruent with recent findings based on nuclear and mitochondrial evidence. The 28 genes retained in the S. densiflora plastid genome are involved in translation and other nonphotosynthetic functions, and we demonstrate that nearly all of the 18 protein-coding genes are under strong purifying selection. Our study confirms the utility of whole plastid genome data in phylogenetic studies of highly modified heterotrophic plants, even when they have substantially elevated rates of substitution.

Next-generation sampling: Pairing genomics with herbarium specimens provides species-level signal in Solidago (Asteraceae)

PREMISE OF THE STUDY

The ability to conduct species delimitation and phylogeny reconstruction with genomic data sets obtained exclusively from herbarium specimens would rapidly enhance our knowledge of large, taxonomically contentious plant genera. In this study, the utility of genotyping by sequencing is assessed in the notoriously difficult genus Solidago (Asteraceae) by attempting to obtain an informative single-nucleotide polymorphism data set from a set of specimens collected between 1970 and 2010.

METHODS

Reduced representation libraries were prepared and Illumina-sequenced from 95 Solidago herbarium specimen DNAs, and resulting reads were processed with the nonreference Universal Network-Enabled Analysis Kit (UNEAK) pipeline. Multidimensional clustering was used to assess the correspondence between genetic groups and morphologically defined species.

RESULTS

Library construction and sequencing were successful in 93 of 95 samples. The UNEAK pipeline identified 8470 single-nucleotide polymorphisms, and a filtered data set was analyzed for each of three Solidago subsections. Although results varied, clustering identified genomic groups that often corresponded to currently recognized species or groups of closely related species.

DISCUSSION

These results suggest that genotyping by sequencing is broadly applicable to DNAs obtained from herbarium specimens. The data obtained and their biological signal suggest that pairing genomics with large-scale herbarium sampling is a promising strategy in species-rich plant groups.

Next-generation sampling: Pairing genomics with herbarium specimens provides species-level signal in Solidago (Asteraceae)

PREMISE OF THE STUDY

The ability to conduct species delimitation and phylogeny reconstruction with genomic data sets obtained exclusively from herbarium specimens would rapidly enhance our knowledge of large, taxonomically contentious plant genera. In this study, the utility of genotyping by sequencing is assessed in the notoriously difficult genus Solidago (Asteraceae) by attempting to obtain an informative single-nucleotide polymorphism data set from a set of specimens collected between 1970 and 2010.

METHODS

Reduced representation libraries were prepared and Illumina-sequenced from 95 Solidago herbarium specimen DNAs, and resulting reads were processed with the nonreference Universal Network-Enabled Analysis Kit (UNEAK) pipeline. Multidimensional clustering was used to assess the correspondence between genetic groups and morphologically defined species.

RESULTS

Library construction and sequencing were successful in 93 of 95 samples. The UNEAK pipeline identified 8470 single-nucleotide polymorphisms, and a filtered data set was analyzed for each of three Solidago subsections. Although results varied, clustering identified genomic groups that often corresponded to currently recognized species or groups of closely related species.

DISCUSSION

These results suggest that genotyping by sequencing is broadly applicable to DNAs obtained from herbarium specimens. The data obtained and their biological signal suggest that pairing genomics with large-scale herbarium sampling is a promising strategy in species-rich plant groups.

The mitochondrial genome of the peach fruit fly, Bactrocera zonata (Saunders) (Diptera: Tephritidae): Complete DNA sequence, genome organization, and phylogenetic analysis with other tephritids using next generation DNA sequencing

Mitochondrial genome can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. The complete 15,935 bp mitochondrial genome of Bactrocera zonata (Diptera: Tephritidae), is assembled from Illumina MiSeq read data. The mitogenome information for B. zonata was compared to the homologous sequences of other tephritids. Annotation indicated that the structure and orientation of 13 protein coding genes (PCGs), 22 tRNA and 2 rRNA sequences were typical of, and similar to, the ten closely related tephritid species. The nucleotide composition shows heavily biased toward As and Ts accounting 73.34% and exhibits a slightly positive AT skew, which is similar to other known tephritid species. All PCGs are initiated by ATN codons, except for cox1 with TCG and atp8 with GTG. Nine PCGs use a common stop codon of TAA or TAG, whereas the remaining four use an incomplete termination codon T or TA likely to be completed by adenylation. All tRNAs have the typical clover-leaf structure, with an exception for trnS((AGN)). Four short intergenic spacers showed high degree of conservation among B. zonata and other ten tephritids. A poly(T) stretch at the 5' end followed by [TA(A)]n-like stretch and a tandem repeats of 39 bp has been observed in CR. The analysis of gene evolutionary rate revealed that the cox1 and atp6 exhibits lowest and highest gene substitution rates, respectively than other genes. The phylogenetic relationships based on Maximum Likelihood method using all protein-coding genes and two ribosomal RNA genes confirmed that B. zonata is closely related to Bactrocera correcta, Bactrocera carambolae, Bactrocera papayae, and Bactrocera philippinensis and Bactrocera dorsalis belonging to B. dorsalis species complex forms a monophyletic clade, which is in accordance with the traditional morphological classification and recent molecular works.

The use of museum specimens with high-throughput DNA sequencers

Natural history collections have long been used by morphologists, anatomists, and taxonomists to probe the evolutionary process and describe biological diversity. These biological archives also offer great opportunities for genetic research in taxonomy, conservation, systematics, and population biology. They allow assays of past populations, including those of extinct species, giving context to present patterns of genetic variation and direct measures of evolutionary processes. Despite this potential, museum specimens are difficult to work with because natural postmortem processes and preservation methods fragment and damage DNA. These problems have restricted geneticists' ability to use natural history collections primarily by limiting how much of the genome can be surveyed. Recent advances in DNA sequencing technology, however, have radically changed this, making truly genomic studies from museum specimens possible. We review the opportunities and drawbacks of the use of museum specimens, and suggest how to best execute projects when incorporating such samples. Several high-throughput (HT) sequencing methodologies, including whole genome shotgun sequencing, sequence capture, and restriction digests (demonstrated here), can be used with archived biomaterials.

Computational analyses of ancient pathogen DNA from herbarium samples: challenges and prospects

The application of DNA sequencing technology to the study of ancient DNA has enabled the reconstruction of past epidemics from genomes of historically important plant-associated microbes. Recently, the genome sequences of the potato late blight pathogen Phytophthora infestans were analyzed from 19th century herbarium specimens. These herbarium samples originated from infected potatoes collected during and after the Irish potato famine. Herbaria have therefore great potential to help elucidate past epidemics of crops, date the emergence of pathogens, and inform about past pathogen population dynamics. DNA preservation in herbarium samples was unexpectedly good, raising the possibility of a whole new research area in plant and microbial genomics. However, the recovered DNA can be extremely fragmented resulting in specific challenges in reconstructing genome sequences. Here we review some of the challenges in computational analyses of ancient DNA from herbarium samples. We also applied the recently developed linkage method to haplotype reconstruction of diploid or polyploid genomes from fragmented ancient DNA.

Using targeted enrichment of nuclear genes to increase phylogenetic resolution in the neotropical rain forest genus Inga (Leguminosae: Mimosoideae)

Evolutionary radiations are prominent and pervasive across many plant lineages in diverse geographical and ecological settings; in neotropical rainforests there is growing evidence suggesting that a significant fraction of species richness is the result of recent radiations. Understanding the evolutionary trajectories and mechanisms underlying these radiations demands much greater phylogenetic resolution than is currently available for these groups. The neotropical tree genus Inga (Leguminosae) is a good example, with ~300 extant species and a crown age of 2-10 MY, yet over 6 kb of plastid and nuclear DNA sequence data gives only poor phylogenetic resolution among species. Here we explore the use of larger-scale nuclear gene data obtained though targeted enrichment to increase phylogenetic resolution within Inga. Transcriptome data from three Inga species were used to select 264 nuclear loci for targeted enrichment and sequencing. Following quality control to remove probable paralogs from these sequence data, the final dataset comprised 259,313 bases from 194 loci for 24 accessions representing 22 Inga species and an outgroup (Zygia). Bayesian phylogenies reconstructed using either all loci concatenated or a gene-tree/species-tree approach yielded highly resolved phylogenies. We used coalescent approaches to show that the same targeted enrichment data also have significant power to discriminate among alternative within-species population histories within the widespread species I. umbellifera. In either application, targeted enrichment simplifies the informatics challenge of identifying orthologous loci associated with de novo genome sequencing. We conclude that targeted enrichment provides the large volumes of phylogenetically-informative sequence data required to resolve relationships within recent plant species radiations, both at the species level and for within-species phylogeographic studies.

From museums to genomics: old herbarium specimens shed light on a C3 to C4 transition

Collections of specimens held by natural history museums are invaluable material for biodiversity inventory and evolutionary studies, with specimens accumulated over 300 years readily available for sampling. Unfortunately, most museum specimens yield low-quality DNA. Recent advances in sequencing technologies, so called next-generation sequencing, are revolutionizing phylogenetic investigations at a deep level. Here, the Illumina technology (HiSeq) was used on herbarium specimens of Sartidia (subfamily Aristidoideae, Poaceae), a small African-Malagasy grass lineage (six species) characteristic of wooded savannas, which is the C3 sister group of Stipagrostis, an important C4 genus from Africa and SW Asia. Complete chloroplast and nuclear ribosomal sequences were assembled for two Sartidia species, one of which (S. perrieri) is only known from a single specimen collected in Madagascar 100 years ago. Partial sequences of a few single-copy genes encoding phosphoenolpyruvate carboxylases (ppc) and malic enzymes (nadpme) were also assembled. Based on these data, the phylogenetic position of Malagasy Sartidia in the subfamily Aristidoideae was investigated and the biogeographical history of this genus was analysed with full species sampling. The evolutionary history of two genes for C4 photosynthesis (ppc-aL1b and nadpme-IV) in the group was also investigated. The gene encoding the C4 phosphoenolpyruvate caroxylase of Stipagrostis is absent from S. dewinteri suggesting that it is not essential in C3 members of the group, which might have favoured its recruitment into a new metabolic pathway. Altogether, the inclusion of historical museum specimens in phylogenomic analyses of biodiversity opens new avenues for evolutionary studies.

Phylogenetic properties of 50 nuclear loci in Medicago (Leguminosae) generated using multiplexed sequence capture and next-generation sequencing

Next-generation sequencing technology has increased the capacity to generate molecular data for plant biological research, including phylogenetics, and can potentially contribute to resolving complex phylogenetic problems. The evolutionary history of Medicago L. (Leguminosae: Trifoliae) remains unresolved due to incongruence between published phylogenies. Identification of the processes causing this genealogical incongruence is essential for the inference of a correct species phylogeny of the genus and requires that more molecular data, preferably from low-copy nuclear genes, are obtained across different species. Here we report the development of 50 novel LCN markers in Medicago and assess the phylogenetic properties of each marker. We used the genomic resources available for Medicago truncatula Gaertn., hybridisation-based gene enrichment (sequence capture) techniques and Next-Generation Sequencing to generate sequences. This alternative proves to be a cost-effective approach to amplicon sequencing in phylogenetic studies at the genus or tribe level and allows for an increase in number and size of targeted loci. Substitution rate estimates for each of the 50 loci are provided, and an overview of the variation in substitution rates among a large number of low-copy nuclear genes in plants is presented for the first time. Aligned sequences of major species lineages of Medicago and its sister genus are made available and can be used in further probe development for sequence-capture of the same markers.

Population genomics of marine fishes: next-generation prospects and challenges

Over the past few years, technological advances have facilitated giant leaps forward in our ability to generate genome-wide molecular data, offering exciting opportunities for gaining new insights into the ecology and evolution of species where genomic information is still limited. Marine fishes are valuable organisms for advancing our understanding of evolution on historical and contemporary time scales, and here we highlight areas in which research on these species is likely to be particularly important in the near future. These include possibilities for gaining insights into processes on ecological time scales, identifying genomic signatures associated with population divergence under gene flow, and determining the genetic basis of phenotypic traits. We also consider future challenges pertaining to the implementation of genome-wide coverage through next-generation sequencing and genotyping methods in marine fishes. Complications associated with fast decay of linkage disequilibrium, as expected for species with large effective population sizes, and the possibility that adaptation is associated with both soft selective sweeps and polygenic selection, leaving complex genomic signatures in natural populations, are likely to challenge future studies. However, the combination of high genome coverage and new statistical developments offers promising solutions. Thus, the next generation of studies is likely to truly facilitate the transition from population genetics to population genomics in marine fishes. This transition will advance our understanding of basic evolutionary processes and will offer new possibilities for conservation and management of valuable marine resources.

The more the merrier: recent hybridization and polyploidy in cardamine

This article describes the use of cytogenomic and molecular approaches to explore the origin and evolution of Cardamine schulzii, a textbook example of a recent allopolyploid, in its ~110-year history of human-induced hybridization and allopolyploidy in the Swiss Alps. Triploids are typically viewed as bridges between diploids and tetraploids but rarely as parental genomes of high-level hybrids and polyploids. The genome of the triploid semifertile hybrid Cardamine × insueta (2n = 24, RRA) was shown to combine the parental genomes of two diploid (2n = 2x = 16) species, Cardamine amara (AA) and Cardamine rivularis (RR). These parental genomes have remained structurally stable within the triploid genome over the >100 years since its origin. Furthermore, we provide compelling evidence that the alleged recent polyploid C. schulzii is not an autohexaploid derivative of C. × insueta. Instead, at least two hybridization events involving C. × insueta and the hypotetraploid Cardamine pratensis (PPPP, 2n = 4x-2 = 30) have resulted in the origin of the trigenomic hypopentaploid (2n = 5x-2 = 38, PPRRA) and hypohexaploid (2n = 6x-2 = 46, PPPPRA). These data show that the semifertile triploid hybrid can promote a merger of three different genomes and demonstrate how important it is to reexamine the routinely repeated textbook examples using modern techniques.