Characterizing DNA preservation in degraded specimens of Amara alpina (Carabidae: Coleoptera)

Developing the Protocol Infrastructure for DNA Sequencing Natural History Collections

Intentionally preserved biological material in natural history collections represents a vast repository of biodiversity. Advances in laboratory and sequencing technologies have made these specimens increasingly accessible for genomic analyses, offering a window into the genetic past of species and often permitting access to information that can no longer be sampled in the wild. Due to their age, preparation and storage conditions, DNA retrieved from museum and herbarium specimens is often poor in yield, heavily fragmented and biochemically modified. This not only poses methodological challenges in recovering nucleotide sequences, but also makes such investigations susceptible to environmental and laboratory contamination. In this paper, we review the practical challenges associated with making the recovery of DNA sequence data from museum collections more routine. We first review key operational principles and issues to address, to guide the decision-making process and dialogue between researchers and curators about when and how to sample museum specimens for genomic analyses. We then outline the range of steps that can be taken to reduce the likelihood of contamination including laboratory set-ups, workflows and working practices. We finish by presenting a series of case studies, each focusing on protocol practicalities for the application of different mainstream methodologies to museum specimens including: (i) shotgun sequencing of insect mitogenomes, (ii) whole genome sequencing of insects, (iii) genome skimming to recover plant plastid genomes from herbarium specimens, (iv) target capture of multi-locus nuclear sequences from herbarium specimens, (v) RAD-sequencing of bird specimens and (vi) shotgun sequencing of ancient bovid bone samples.

Complete mitochondrial genomes from museum specimens clarify millipede evolution in the Eastern Arc Mountains

The Eastern Arc Mountains in Tanzania represent a hotspot for biological diversity of global importance. The level of endemism is high, and Eastern Arc biodiversity has been studied extensively in vertebrates and invertebrates, including millipedes. However, millipede evolution is vastly understudied at the molecular level. Therefore, we used next-generation ‘shotgun’ sequencing to obtain mitochondrial genome sequences of 26 museum specimens, representing six genera and 12 millipede species found across the Eastern Arc Mountains. Bayesian and maximum likelihood methods yielded consistent topologies with high node support, confirming a high level of congruence between molecular and morphological analyses. The only exception was a Tropostreptus sigmatospinus individual from Zanzibar, which was placed outside an otherwise monophyletic group consisting of mainland individuals of the same assumed species. For two species with a distribution across several mountain blocks (Tropostreptus sigmatospinus and Tropostreptus hamatus), each mountain population represents a distinct monophyletic lineage. In contrast, we also observe that distinct species exist sympatrically in the same montane forests, indicative of older speciation events that are not defined by current forest distribution. Our results are important for understanding speciation mechanisms in montane rain forests and highlight that ethanol-preserved invertebrates exhibit a tremendous potential for genomic analyses.

Old specimens for old branches: Assessing effects of sample age in resolving a rapid Neotropical radiation of squirrels

Ultraconserved Elements (UCEs) have been useful to resolve challenging phylogenies of non-model clades, unpuzzling long-conflicted relationships in key branches of the Tree of Life at both deep and shallow levels. UCEs are often reliably recovered from historical samples, unlocking a vast number of preserved natural history specimens for analysis. However, the extent to which sample age and preservation method impact UCE recovery as well as downstream inferences remains unclear. Furthermore, there is an ongoing debate on how to curate, filter, and properly analyze UCE data when locus recovery is uneven across sample age and quality. In the present study we address these questions with an empirical dataset composed of over 3800 UCE loci from 219 historical and modern samples of Sciuridae, a globally distributed and ecologically important family of rodents. We provide a genome-scale phylogeny of two squirrel subfamilies (Sciurillinae and Sciurinae: Sciurini) and investigate their placement within Sciuridae. For historical specimens, recovery of UCE loci and mean length per locus were inversely related to sample age; deeper sequencing improved the number of UCE loci recovered but not locus length. Most of our phylogenetic inferences-performed on six datasets with alternative data-filtering strategies, and using three distinct optimality criteria-resulted in distinct topologies. Datasets containing more loci (40% and 50% taxa representativeness matrices) yielded more concordant topologies and higher support values than strictly filtered datasets (60% matrices) particularly with IQ-Tree and SVDquartets, while filtering based on information content provided better topological resolution for inferences with the coalescent gene-tree based approach in ASTRAL-III. We resolved deep relationships in Sciuridae (including among the five currently recognized subfamilies) and relationships among the deepest branches of Sciurini, but conflicting relationships remain at both genus- and species-levels for the rapid Neotropical tree squirrel radiation. Our results suggest that phylogenomic consensus can be difficult and heavily influenced by the age of available samples and the filtering steps used to optimize dataset properties.

New Light on Historical Specimens Reveals a New Species of Ladybird (Coleoptera: Coccinellidae): Morphological, Museomic, and Phylogenetic Analyses

Natural history collections house an important source of genetic data from yet unexplored biological diversity. Molecular data from museum specimens remain underexploited, which is mainly due to the degradation of DNA from specimens over time. However, Next-Generation Sequencing (NGS) technology can now be used to sequence "old" specimens. Indeed, many of these specimens are unique samples of nomenclatural types and can be crucial for resolving systematic or biogeographic scientific questions. Two ladybird beetle specimens from Patagonia corresponding to a new species of the genus Eriopis Mulsant were found in the collections of the Muséum national d'Histoire naturelle (MNHN), Paris. Here, we describe Eriopis patagonia Salazar, sp. nov. Total DNA of one of the two specimens was sequenced by NGS using a paired-end Illumina approach. We reconstruct and characterize the mitochondrial genome of this species (16,194 bp). Then, the protein-coding genes (PCGs) and ribosomal RNAs (rRNAs) were used to infer by maximum likelihood and Bayesian Inference the phylogenetic position of E. patagonia among 27 representatives of Coccinellidae. Phylogenetic analysis confirmed the position of Eriopis as sister group to Cycloneda Crotch. Hence, we highlight the high potential of sequencing technology for extracting molecular information from old specimens, which are used here for the systematic study of a genus, while demonstrating the importance of preserving biological collections.

DNA from resin-embedded organisms: Past, present and future

Past claims have been made for fossil DNA recovery from various organisms (bacteria, plants, insects and mammals, including humans) dating back in time from thousands to several million years BP. However, many of these recoveries, especially those described from million-year-old amber (fossil resin), have faced criticism as being the result of modern environmental contamination and for lack of reproducibility. Using modern genomic techniques, DNA can be obtained with confidence from a variety of substrates (e.g. bones, teeth, gum, museum specimens and fossil insects) of different ages, albeit always less than one million years BP, and results can also be obtained from much older materials using palaeoproteomics. Nevertheless, new attempts to determine if ancient DNA (aDNA) is present in insects preserved in 40 000-year old sub-fossilised resin, the precursor of amber, have been unsuccessful or not well documented. Resin-embedded specimens are therefore regarded as unsuitable for genetic studies. However, we demonstrate here, for the first time, that although a labile molecule, DNA is still present in platypodine beetles (Coleoptera: Curculionidae) embedded in six-year-old and two-year-old resin pieces from Hymenaea verrucosa (Angiospermae: Fabaceae) collected in Madagascar. We describe an optimised method which meets all the requirements and precautions for aDNA experiments for our purpose: to explore the DNA preservation limits in resin. Our objective is far from starting an uncontrolled search for aDNA in amber as it was in the past, but to start resolving basic aspects from the DNA preservation in resin and search from the most modern samples to the ancient ones, step by step. We conclude that it is therefore possible to study genomics from resin-embedded organisms, although the time limits remain to be determined.

Declining genetic diversity of European honeybees along the twentieth century

The European honeybee (Apis mellifera) is a key pollinator and has in the last decades suffered significant population decline. A combination of factors, including decrease in genetic diversity and introduction of Varroa mites, have been suggested to be responsible for these losses, but no definitive cause has yet been appointed. In Europe not only have wild colonies been severely affected, but managed hives have had a massive decline in numbers. To test the hypothesis that honeybees’ genetic diversity has decreased in the recent past, we used reduced representation genome sequencing of 40 historical honeybee specimens collected in Natural History collections across Europe and compared them to genomic data from 40 individuals from extant populations (collected post 2006). Our results are consistent with the existence of five evolutionary lineages as previously described, and show a decrease in genetic diversity between historical and extant individuals of the same lineage, as well as high levels of admixture in historical specimens. Our data confirm that a loss of genetic diversity has occurred during the last century, potentially increasing honeybees’ vulnerability to contemporary ecological and anthropogenic stressors.

A method to generate multilocus barcodes of pinned insect specimens using MiSeq

For molecular insect identification, amplicon sequencing methods are recommended because they offer a cost-effective approach for targeting small sets of informative genes from multiple samples. In this context, high-throughput multilocus amplicon sequencing has been achieved using the MiSeq Illumina sequencing platform. However, this approach generates short gene fragments of <500 bp, which then have to be overlapped using bioinformatics to achieve longer sequence lengths. This increases the risk of generating chimeric sequences or leads to the formation of incomplete loci. Here, we propose a modified nested amplicon sequencing method for targeting multiple loci from pinned insect specimens using the MiSeq Illumina platform. The modification exists in using a three-step nested PCR approach targeting near full-length loci in the initial PCR and subsequently amplifying short fragments of between 300 and 350 bp for high-throughput sequencing using Illumina chemistry. Using this method, we generated 407 sequences of three loci from 86% of all the specimens sequenced. Out of 103 pinned bee specimens of replicated species, 71% passed the 95% sequence similarity threshold between species replicates. This method worked best for pinned specimens aged between 0 and 5 years, with a limit of 10 years for pinned and 14 years for ethanol-preserved specimens. Hence, our method overcomes some of the challenges of amplicon sequencing using short read next generation sequencing and improves the possibility of creating high-quality multilocus barcodes from insect collections.

Evaluation of DNA degradation and establishment of a degradation analysis model for Lepidoptera specimens

Millions of museum specimens are integral to biodiversity studies; however, DNA degradation may limit the ability to obtain DNA sequences. In this study, a degradation analysis model for Lepidoptera specimens was established. Based on this model, we revealed the characteristics of DNA fragment distribution caused by external DNA damage factors during specimen preservation. We found that the degree of DNA degradation increased over time; DNA degradation of spread and dried adult specimens was significantly higher than that in the folded and formalin-fixed larval specimens. However, the effects of folding wings on DNA degradation and the effects of the preservation method/stage (formalin-fixed larval vs air-dried adult specimens) were different for different species.

Utilizing field collected insects for next generation sequencing: Effects of sampling, storage, and DNA extraction methods

DNA sequencing technologies continue to advance the biological sciences, expanding opportunities for genomic studies of non-model organisms for basic and applied questions. Despite these opportunities, many next generation sequencing protocols have been developed assuming a substantial quantity of high molecular weight DNA (>100 ng), which can be difficult to obtain for many study systems. In particular, the ability to sequence field-collected specimens that exhibit varying levels of DNA degradation remains largely unexplored. In this study we investigate the influence of five traditional insect capture and curation methods on Double-Digest Restriction Enzyme Associated DNA (ddRAD) sequencing success for three wild bee species. We sequenced a total of 105 specimens (between 7-13 specimens per species and treatment). We additionally investigated how different DNA quality metrics (including pre-sequence concentration and contamination) predicted downstream sequencing success, and also compared two DNA extraction methods. We report successful library preparation for all specimens, with all treatments and extraction methods producing enough highly reliable loci for population genetic analyses. Although results varied between species, we found that specimens collected by net sampling directly into 100% EtOH, or by passive trapping followed by 100% EtOH storage before pinning tended to produce higher quality ddRAD assemblies, likely as a result of rapid specimen desiccation. Surprisingly, we found that specimens preserved in propylene glycol during field sampling exhibited lower-quality assemblies. We provide recommendations for each treatment, extraction method, and DNA quality assessment, and further encourage researchers to consider utilizing a wider variety of specimens for genomic analyses.

Genome skimming is a low-cost and robust strategy to assemble complete mitochondrial genomes from ethanol preserved specimens in biodiversity studies

Global loss of biodiversity is an ongoing process that concerns both local and global authorities. Studies of biodiversity mainly involve traditional methods using morphological characters and molecular protocols. However, conventional methods are a time consuming and resource demanding task. The development of high-throughput sequencing (HTS) techniques has reshaped the way we explore biodiversity and opened a path to new questions and novel empirical approaches. With the emergence of HTS, sequencing the complete mitochondrial genome became more accessible, and the number of genome sequences published has increased exponentially during the last decades. Despite the current state of knowledge about the potential of mitogenomics in phylogenetics, this is still a relatively under-explored area for a multitude of taxonomic groups, especially for those without commercial relevance, non-models organisms and with preserved DNA. Here we take the first step to assemble and annotate the genomes from HTS data using a new protocol of genome skimming which will offer an opportunity to extend the field of mitogenomics to under-studied organisms. We extracted genomic DNA from specimens preserved in ethanol. We used Nextera XT DNA to prepare indexed paired-end libraries since it is a powerful tool for working with diverse samples, requiring a low amount of input DNA. We sequenced the samples in two different Illumina platform (MiSeq or NextSeq 550). We trimmed raw reads, filtered and had their quality tested accordingly. We performed the assembly using a baiting and iterative mapping strategy, and the annotated the putative mitochondrion through a semi-automatic procedure. We applied the contiguity index to access the completeness of each new mitogenome. Our results reveal the efficiency of the proposed method to recover the whole mitogenomes of preserved DNA from non-model organisms even if there are gene rearrangement in the specimens. Our findings suggest the potential of combining the adequate platform and library to the genome skimming as an innovative approach, which opens a new range of possibilities of its use to obtain molecular data from organisms with different levels of preservation.

Biodiversity seen through the perspective of insects: 10 simple rules on methodological choices and experimental design for genomic studies

Massively parallel DNA sequencing opens up opportunities for bridging multiple temporal and spatial dimensions in biodiversity research, thanks to its efficiency to recover millions of nucleotide polymorphisms. Here, we identify the current status, discuss the main challenges, and look into future perspectives on biodiversity genomics focusing on insects, which arguably constitute the most diverse and ecologically important group among all animals. We suggest 10 simple rules that provide a succinct step-by-step guide and best-practices to anyone interested in biodiversity research through the study of insect genomics. To this end, we review relevant literature on biodiversity and evolutionary research in the field of entomology. Our compilation is targeted at researchers and students who may not yet be specialists in entomology or molecular biology. We foresee that the genomic revolution and its application to the study of non-model insect lineages will represent a major leap to our understanding of insect diversity.

Using host species traits to understand the Wolbachia infection distribution across terrestrial beetles

Knowledge of Wolbachia prevalence with respect to its hosts is restricted mainly to taxonomic/phylogenetic context. In contrast, relations between infection and most host’s ecological and biological traits are poorly understood. This study aimed to elaborate on relations between bacteria and its beetle hosts in taxonomic and the ecological contexts. In particular, the goal is to verify which ecological and biological traits of beetles could cause them to be prone to be infected. Verification of Wolbachia infection status across 297 beetle taxa showed that approximately 27% of taxa are infected by supergroups A and B. Only minor support for coevolution between bacteria and its beetle hosts was observed in some genera of beetles, but in general coevolution between beetles and Wolbachia was rejected. Some traits of beetles were found to be unrelated to Wolbachia prevalence (type of range and thermal preferences); some traits were related with ambiguous effects (habitats, distribution, mobility and body size); some were substantially related (reproduction mode and trophy). The aforementioned summary does not show obvious patterns of Wolbachia prevalence and diversity in relation to host taxonomy, biology, and ecology. As both Wolbachia and Coleoptera are diverse groups, this lack of clear patterns is probably a reflection of nature, which is characterised by highly diversified and probably unstable relations.

Current state of knowledge on Wolbachia infection among Coleoptera: a systematic review

Background

Despite great progress in studies on Wolbachia infection in insects, the knowledge about its relations with beetle species, populations and individuals, and the effects of bacteria on these hosts, is still unsatisfactory. In this review we summarize the current state of knowledge about Wolbachia occurrence and interactions with Coleopteran hosts.

Methods

An intensive search of the available literature resulted in the selection of 86 publications that describe the relevant details about Wolbachia presence among beetles. These publications were then examined with respect to the distribution and taxonomy of infected hosts and diversity of Wolbachia found in beetles. Sequences of Wolbachia genes (16S rDNA, ftsZ) were used for the phylogenetic analyses.

Results

The collected publications revealed that Wolbachia has been confirmed in 204 beetle species and that the estimated average prevalence of this bacteria across beetle species is 38.3% and varies greatly across families and genera (0–88% infected members) and is much lower (c. 13%) in geographic studies. The majority of the examined and infected beetles were from Europe and East Asia. The most intensively studied have been two groups of herbivorous beetles: Curculionidae and Chrysomelidae. Coleoptera harbor Wolbachia belonging to three supergroups: F found in only three species, and A and B found in similar numbers of beetles (including some doubly infected); however the latter two were most prevalent in different families. A total of 59% of species with precise data were found to be totally infected. Single infections were found in 69% of species and others were doubly- or multiply-infected. Wolbachia caused numerous effects on its beetle hosts, including selective sweep with host mtDNA (found in 3% of species), cytoplasmic incompatibility (detected in c. 6% of beetles) and other effects related to reproduction or development (like male-killing, possible parthenogenesis or haplodiploidy induction, and egg development). Phylogenetic reconstructions for Wolbachia genes rejected cospeciation between these bacteria and Coleoptera, with minor exceptions found in some Hydraenidae, Curculionidae and Chrysomelidae. In contrast, horizontal transmission of bacteria has been suspected or proven in numerous cases (e.g., among beetles sharing habitats and/or host plants).

Discussion

The present knowledge about Wolbachia infection across beetle species and populations is very uneven. Even the basic data about infection status in species and frequency of infected species across genera and families is very superficial, as only c. 0.15% of all beetle species have been tested so far. Future studies on Wolbachia diversity in Coleoptera should still be based on the Multi-locus Sequence Typing system, and next-generation sequencing technologies will be important for uncovering Wolbachia relations with host evolution and ecology, as well as with other, co-occurring endosymbiotic bacteria.

Ancient DNA Resolves the History of Tetragnatha (Araneae, Tetragnathidae) Spiders on Rapa Nui

Rapa Nui is one of the most remote islands in the world. As a young island, its biota is a consequence of both natural dispersals over the last ~1 million years and recent human introductions. It therefore provides an opportunity to study a unique community assemblage. Here, we extract DNA from museum-preserved and newly field-collected spiders from the genus Tetragnatha to explore their history on Rapa Nui. Using an optimized protocol to recover ancient DNA from museum-preserved spiders, we sequence and assemble partial mitochondrial genomes from nine Tetragnatha species, two of which were found on Rapa Nui, and estimate the evolutionary relationships between these and other Tetragnatha species. Our phylogeny shows that the two Rapa Nui species are not closely related. One, the possibly extinct, T. paschae, is nested within a circumtropical species complex (T. nitens), and the other (Tetragnatha sp. Rapa Nui) appears to be a recent human introduction. Our results highlight the power of ancient DNA approaches in identifying cryptic and rare species, which can contribute to our understanding of the global distribution of biodiversity in all taxonomic lineages.

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.

Singing from the grave: DNA from a 180 year old type specimen confirms the identity of Chrysoperla carnea (Stephens)

Historically serving as repositories for morphologically-based taxonomic research, natural history collections are now increasingly being targeted in studies utilizing DNA data. The development of advanced molecular techniques has facilitated extraction of useable DNA from old specimens, including type material. Sequencing diagnostic molecular markers from type material enables accurate species designation, especially where modern taxonomic hypotheses confirm morphologically cryptic species complexes. One such example is Chrysoperla carnea (Stephens), which belongs to a complex of about 20 cryptic species, most of which can only be reliably distinguished by their pre-mating courtship songs or by DNA analysis. The subtle morphological variation in the group has led to disagreement over the previous designation of the lectotype for C. carnea, an issue that has been further compounded because Chrysoperla carnea is a highly valued biological control agent in arable crops. Archival DNA extraction and sequencing from the 180 year old lectotype specimen, combined with Bayesian and Likelihood based phylogenetic analyses of modern specimens from the entire complex, were used to establish unambiguously the true identity of Chrysoperla carnea.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication

The domestication of the horse ∼ 5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski's horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the "cost of domestication" hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.