Sequence capture phylogenomics of historical ethanol‐preserved museum specimens: Unlocking the rest of the vault

Natural history collections play a crucial role in biodiversity research, and museum specimens are increasingly being incorporated into modern genetics‐based studies. Sequence capture methods have proven incredibly useful for phylogenomics, providing the additional ability to sequence historical museum specimens with highly degraded DNA, which until recently have been deemed less valuable for genetic work. The successful sequencing of ultraconserved elements (UCEs) from historical museum specimens has been demonstrated on multiple tissue types including dried bird skins, formalin‐fixed squamates and pinned insects. However, no study has thoroughly demonstrated this approach for historical ethanol‐preserved museum specimens. Alongside sequencing of “fresh” specimens preserved in >95% ethanol and stored at ?80°C, we used extraction techniques specifically designed for degraded DNA coupled with sequence capture protocols to sequence UCEs from historical museum specimens preserved in 70%–80% ethanol and stored at room temperature, the standard for such ethanol‐preserved museum collections. Across 35 fresh and 15 historical museum samples of the arachnid order Opiliones, an average of 345 UCE loci were included in phylogenomic matrices, with museum samples ranging from six to 495 loci. We successfully demonstrate the inclusion of historical ethanol‐preserved museum specimens in modern sequence capture phylogenomic studies, show a high frequency of variant bases at the species and population levels, and from off‐target reads successfully recover multiple loci traditionally sequenced in multilocus studies including mitochondrial loci and nuclear rRNA loci. The methods detailed in this study will allow researchers to potentially acquire genetic data from millions of ethanol‐preserved museum specimens held in collections worldwide.     

Amplification success rate of DNA from museum skin collections: a case study of stoats from 18 museums

There are large museum collections of mammalian skins and we wished to determine their usefulness for DNA‐based evolutionary and conservation studies. Methods derived from the ancient DNA approach were used to process samples from skins of the stoat (Mustela erminea) from 18 museums in 11 countries. Successful polymerase chain reaction (PCR) amplification was achieved in 56%, 46% and 40% of all samples for 0.65‐, 0.70‐ and 0.78‐kb PCR products of mitochondrial DNA, respectively. Some of the best‐preserved skin samples were those in tight‐fitting containers in a dry and cold environment. With care in their preservation, mammalian skin collections may be a good source of DNA.     

DNA content and distribution in ancient feathers and potential to reconstruct the plumage of extinct avian taxa

Feathers are known to contain amplifiable DNA at their base (calamus) and have provided an important genetic source from museum specimens. However, feathers in subfossil deposits generally only preserve the upper shaft and feather ‘vane’ which are thought to be unsuitable for DNA analysis. We analyse subfossil moa feathers from Holocene New Zealand rockshelter sites and demonstrate that both ancient DNA and plumage information can be recovered from their upper portion, allowing species identification and a means to reconstruct the appearance of extinct taxa. These ancient DNA sequences indicate that the distal portions of feathers are an untapped resource for studies of museum, palaeontological and modern specimens. We investigate the potential to reconstruct the plumage of pre-historically extinct avian taxa using subfossil remains, rather than assuming morphological uniformity with closely related extant taxa. To test the notion of colour persistence in subfossil feathers, we perform digital comparisons of feathers of the red-crowned parakeet (Cyanoramphus novaezelandiae novaezelandiae) excavated from the same horizons as the moa feathers, with modern samples. The results suggest that the coloration of the moa feathers is authentic, and computer software is used to perform plumage reconstructions of moa based on subfossil remains.     

Molecular and Immunohistochemical Characterization of Historical Long-Term Preserved Fixed Tissues from Different Human Organs

University and museum collections are very important sources of biological samples that can be used to asses the past and present genetic diversity of many species. Modern genetic and immunohistochemical techniques can be used on long-term preserved fixed tissues from museum specimens to answer epidemiological questions. A proof of principle was established to apply modern molecular genetics and immunohistochemical methods to these old specimens and to verify the original diagnosis. We analysed 19 specimens from our university collection including human organs that had been in fixative for more than 80 years. The tissues originated from lung, colon, brain, heart, adrenal gland, uterus and skin. We isolated amplifiable DNA from these wet preparations and performed mutational analysis of BRAF, KRAS and EGFR. The tissues were also embedded in paraffin and used for modern histology and immunohistochemistry. Our data show that amplifiable DNA is extractable and ranged from 0.25 to 22.77 μg of total DNA. In three specimens BRAF^V600E or KRAS^G12D mutations were found. Additionally, expression of different proteins like vimentin and GFAP was detected immunohistochemical in six investigated specimens. On the basis of our results the original diagnosis was altered in three specimens. Our work showed that it is possible to extract amplifiable DNA suitable for sequence analysis from long-term fixed tissue. Furthermore, histology and immunohistochemistry is feasible in specimens fixed long time ago. We conclude that these old preparations are suitable for further epidemiological research and that our methods open up new opportunities for future studies.     

Note on using nuclear 28S rDNA for sequencing ancient and strongly degraded insect DNA

A protocol using insect specimens or parts thereof allows for sequencing of sections of nuclear 28S rDNA. In the present note it is demonstrated that this protocol can readily be applied to strongly degraded DNA (ancient, fixed or contaminated). Primers that are specifically designed to discriminate against human DNA but also other non‐arthropod species are tested on a range of species covering all insect groups (59 insect species from all 33 orders). Additionally, the samples represent a selection of various, mostly DNA‐degrading, preservation methods, including the most common fixatives used for morphological investigations and for long‐term storage in collections. Successful amplification was possible for all tested samples including ca. 200 year‐old dried museum specimens as well as for over 4000 year‐old fossil insects embedded in copal. When the NCBI database contained information on the tested species an unambiguous taxonomic discrimination was possible. This approach is based on a standardized protocol that guarantees easy application. This note presents primer pairs for 28S rDNA that can be a useful tool for ancient DNA (aDNA) research.     

Complete mitochondrial genomes of eleven extinct or possibly extinct bird species

Natural history museum collections represent a vast source of ancient and historical DNA samples from extinct taxa that can be utilized by high‐throughput sequencing tools to reveal novel genetic and phylogenetic information about them. Here, we report on the successful sequencing of complete mitochondrial genome sequences (mitogenomes) from eleven extinct bird species, using de novo assembly of short sequences derived from toepad samples of degraded DNA from museum specimens. For two species (the Passenger Pigeon Ectopistes migratorius and the South Island Piopio Turnagra capensis), whole mitogenomes were already available from recent studies, whereas for five others (the Great Auk Pinguinis impennis, the Imperial Woodpecker Campehilus imperialis, the Huia Heteralocha acutirostris, the Kauai Oo Moho braccathus and the South Island Kokako Callaeas cinereus), there were partial mitochondrial sequences available for comparison. For all seven species, we found sequence similarities of >98%. For the remaining four species (the Kamao Myadestes myadestinus, the Paradise Parrot Psephotellus pulcherrimus, the Ou Psittirostra psittacea and the Lesser Akialoa Akialoa obscura), there was no sequence information available for comparison, so we conducted blast searches and phylogenetic analyses to determine their phylogenetic positions and identify their closest extant relatives. These mitogenomes will be valuable for future analyses of avian phylogenetics and illustrate the importance of museum collections as repositories for genomics resources.     

Genomics and museum specimens

Nearly 25 years ago, Allan Wilson and colleagues isolated DNA sequences from museum specimens of kangaroo rats (Dipodomys panamintinus) and compared these sequences with those from freshly collected animals (Thomas et al. 1990 ). The museum specimens had been collected up to 78 years earlier, so the two samples provided a direct temporal comparison of patterns of genetic variation. This was not the first time DNA sequences had been isolated from preserved material, but it was the first time it had been carried out with a population sample. Population geneticists often try to make inferences about the influence of historical processes such as selection, drift, mutation and migration on patterns of genetic variation in the present. The work of Wilson and colleagues was important in part because it suggested a way in which population geneticists could actually study genetic change in natural populations through time, much the same way that experimentalists can do with artificial populations in the laboratory. Indeed, the work of Thomas et al. ( 1990 ) spawned dozens of studies in which museum specimens were used to compare historical and present‐day genetic diversity (reviewed in Wandeler et al. 2007 ). All of these studies, however, were limited by the same fundamental problem: old DNA is degraded into short fragments. As a consequence, these studies mostly involved PCR amplification of short templates, usually short stretches of mitochondrial DNA or microsatellites. In this issue, Bi et al. ( 2013 ) report a breakthrough that should open the door to studies of genomic variation in museum specimens. They used target enrichment (exon capture) and next‐generation (Illumina) sequencing to compare patterns of genetic variation in historic and present‐day population samples of alpine chipmunks (Tamias alpinus) (Fig. 1). The historic samples came from specimens collected in 1915, so the temporal span of this comparison is nearly 100 years.     

Unlocking the vault: next‐generation museum population genomics

Natural history museum collections provide unique resources for understanding how species respond to environmental change, including the abrupt, anthropogenic climate change of the past century. Ideally, researchers would conduct genome‐scale screening of museum specimens to explore the evolutionary consequences of environmental changes, but to date such analyses have been severely limited by the numerous challenges of working with the highly degraded DNA typical of historic samples. Here, we circumvent these challenges by using custom, multiplexed, exon capture to enrich and sequence ~11 000 exons (~4 Mb) from early 20th‐century museum skins. We used this approach to test for changes in genomic diversity accompanying a climate‐related range retraction in the alpine chipmunks (Tamias alpinus) in the high Sierra Nevada area of California, USA. We developed robust bioinformatic pipelines that rigorously detect and filter out base misincorporations in DNA derived from skins, most of which likely resulted from postmortem damage. Furthermore, to accommodate genotyping uncertainties associated with low‐medium coverage data, we applied a recently developed probabilistic method to call single‐nucleotide polymorphisms and estimate allele frequencies and the joint site frequency spectrum. Our results show increased genetic subdivision following range retraction, but no change in overall genetic diversity at either nonsynonymous or synonymous sites. This case study showcases the advantages of integrating emerging genomic and statistical tools in museum collection‐based population genomic applications. Such technical advances greatly enhance the value of museum collections, even where a pre‐existing reference is lacking and points to a broad range of potential applications in evolutionary and conservation biology.     

Museum specimens provide reliable SNP data for population genomic analysis of a widely distributed but threatened cockatoo species

Natural history museums harbour a plethora of biological specimens which are of potential use in population and conservation genetic studies. Although technical advancements in museum genomics have enabled genome‐wide markers to be generated from aged museum specimens, the suitability of these data for robust biological inference is not well characterized. The aim of this study was to test the utility of museum specimens in population and conservation genomics by assessing the biological and technical validity of single nucleotide polymorphism (SNP) data derived from such samples. To achieve this, we generated thousands of SNPs from 47 red‐tailed black cockatoo (Calyptorhychus banksii) traditional museum samples (i.e. samples that were not collected with the primary intent of DNA analysis) and 113 fresh tissue samples (cryopreserved liver/muscle) using a restriction site‐associated DNA marker approach (DArTseq^?). Thousands of SNPs were successfully generated from most of the traditional museum samples (with a mean age of 44 years, ranging from 5 to 123 years), although 38% did not provide useful data. These SNPs exhibited higher error rates and contained significantly more missing data compared with SNPs from fresh tissue samples, likely due to considerable DNA fragmentation. However, based on simulation results, the level of genotyping error had a negligible effect on inference of population structure in this species. We did identify a bias towards low diversity SNPs in older samples that appears to compromise temporal inferences of genetic diversity. This study demonstrates the utility of a RADseq‐based method to produce reliable genome‐wide SNP data from traditional museum specimens.     

Extracting DNA from ‘jaws’: high yield and quality from archived tiger shark (Galeocerdo cuvier) skeletal material

Archived specimens are highly valuable sources of DNA for retrospective genetic/genomic analysis. However, often limited effort has been made to evaluate and optimize extraction methods, which may be crucial for downstream applications. Here, we assessed and optimized the usefulness of abundant archived skeletal material from sharks as a source of DNA for temporal genomic studies. Six different methods for DNA extraction, encompassing two different commercial kits and three different protocols, were applied to material, so‐called bio‐swarf, from contemporary and archived jaws and vertebrae of tiger sharks (Galeocerdo cuvier). Protocols were compared for DNA yield and quality using a qPCR approach. For jaw swarf, all methods provided relatively high DNA yield and quality, while large differences in yield between protocols were observed for vertebrae. Similar results were obtained from samples of white shark (Carcharodon carcharias). Application of the optimized methods to 38 museum and private angler trophy specimens dating back to 1912 yielded sufficient DNA for downstream genomic analysis for 68% of the samples. No clear relationships between age of samples, DNA quality and quantity were observed, likely reflecting different preparation and storage methods for the trophies. Trial sequencing of DNA capture genomic libraries using 20 000 baits revealed that a significant proportion of captured sequences were derived from tiger sharks. This study demonstrates that archived shark jaws and vertebrae are potential high‐yield sources of DNA for genomic‐scale analysis. It also highlights that even for similar tissue types, a careful evaluation of extraction protocols can vastly improve DNA yield.     

How old can we go? Evaluating the age limit for effective DNA recovery from historical insect specimens

Historical museum specimens are valuable for exploring population genetics and evolutionary questions because they can provide snapshots of morphological and genetic characteristics from populations over space and time. Unfortunately, DNA found in older museum specimens is frequently degraded, so obtaining genotypes from many individual samples necessary for rigorous molecular population genetic studies is challenging. Previous studies have varied greatly in their success at obtaining genotypes from older preserved insect material. Many well-intentioned collection curators have used research results showing poor preservation of DNA preserved in museum specimens to inform curatorial best practices, in some cases choosing not to allow DNA extraction by destructive sampling because, in their estimation, the likelihood of success would be low. Recent methodological advances in DNA extraction, amplification, and genotyping have allowed some researchers to include mid-19th century samples in molecular genetic analyses. Here we present a robust, high-throughput, and low-cost DNA extraction and genotyping protocol for historical insect specimens employing restriction digests of PCR products followed by high sensitivity electrophoresis. Using this technique, we obtained mitochondrial haplotypes for 100% of 48 New World Junonia butterfly specimens (Nymphalidae) ranging in age from pre-1813 to 1909 and show that the haplotype frequencies obtained are statistically indistinguishable from 20th-century and contemporary reference populations of Junonia (1632 specimens) matched by geographic region. As most extant insect specimens were collected after 1813, based on our findings we would expect that many or even most pinned specimens preserved in museum collections contain usable DNA for mitochondrial haplotyping.     

Micro‐CT X‐rays do not fragment DNA in preserved bird skins

Most zoological systematics studies are currently based on morphological features, molecular traits or a combination of both to reconstruct animals’ phylogenetic history. Increasingly, morphological studies of museum specimens are using X‐ray computed tomography to visualize internal morphology, because of its ‘non‐destructive’ nature. However, it is not known whether CT can fragment the size of DNA extracted from museum specimens, as has been demonstrated to occur in living cells. This question is of paramount importance for collections based research because X‐rays may reduce the amount of data obtainable from specimens. In our study, we tested whether exposure of museum bird skins to typical CT X‐ray energies (for visualization of the skeleton) increased DNA strand fragmentation, a key factor for the success of downstream molecular applications. For the present study, we extracted DNA from shavings of 24 prepared and dried bird skins (100⁺ years) footpads before and after CT scanning. The pre‐ and post‐CT fragmentation profiles were assessed using a capillary electrophoresis high‐precision instrument (Agilent Bioanalyzer). Comparison of the most common strand length in each DNA sample (relative mass) revealed no significant difference unexposed and exposed tissue (paired t‐test p = 0.463). In conclusion, we found no further quantifiable degradation of DNA strand length under standard X‐ray exposure obtained from our bird skins sample. Differences in museum preservation techniques probably had a greater effect on variation of pre‐CT DNA fragmentation.     

Phylogenomics using formalin‐fixed and 100+ year‐old intractable natural history specimens

Museum specimens provide a wealth of information to biologists, but obtaining genetic data from formalin‐fixed and fluid‐preserved specimens remains challenging. While DNA sequences have been recovered from such specimens, most approaches are time‐consuming and produce low data quality and quantity. Here, we use a modified DNA extraction protocol combined with high‐throughput sequencing to recover DNA from formalin‐fixed and fluid‐preserved snakes that were collected over a century ago and for which little or no modern genetic materials exist in public collections. We successfully extracted DNA and sequenced ultraconserved elements ( = 2318 loci) from 10 fluid‐preserved snakes and included them in a phylogeny with modern samples. This phylogeny demonstrates the general use of such specimens in phylogenomic studies and provides evidence for the placement of enigmatic snakes, such as the rare and never‐before sequenced Indian Xylophis stenorhynchus. Our study emphasizes the relevance of museum collections in modern research and simultaneously provides a protocol that may prove useful for specimens that have been previously intractable for DNA sequencing.     

Ancient DNA analysis identifies marine mollusc shells as new metagenomic archives of the past

Marine mollusc shells enclose a wealth of information on coastal organisms and their environment. Their life history traits as well as (palaeo‐) environmental conditions, including temperature, food availability, salinity and pollution, can be traced through the analysis of their shell (micro‐) structure and biogeochemical composition. Adding to this list, the DNA entrapped in shell carbonate biominerals potentially offers a novel and complementary proxy both for reconstructing palaeoenvironments and tracking mollusc evolutionary trajectories. Here, we assess this potential by applying DNA extraction, high‐throughput shotgun DNA sequencing and metagenomic analyses to marine mollusc shells spanning the last ~7,000 years. We report successful DNA extraction from shells, including a variety of ancient specimens, and find that DNA recovery is highly dependent on their biomineral structure, carbonate layer preservation and disease state. We demonstrate positive taxonomic identification of mollusc species using a combination of mitochondrial DNA genomes, barcodes, genome‐scale data and metagenomic approaches. We also find shell biominerals to contain a diversity of microbial DNA from the marine environment. Finally, we reconstruct genomic sequences of organisms closely related to the Vibrio tapetis bacteria from Manila clam shells previously diagnosed with Brown Ring Disease. Our results reveal marine mollusc shells as novel genetic archives of the past, which opens new perspectives in ancient DNA research, with the potential to reconstruct the evolutionary history of molluscs, microbial communities and pathogens in the face of environmental changes. Other future applications include conservation of endangered mollusc species and aquaculture management.     

Extraction of nuclear DNA from bone of skeletonized and fluid‐preserved museum specimens

Obtaining DNA sequences, particularly nuclear DNA, from museum specimens is challenging. We sequenced nuclear DNA from small bone fragments of skeletonized and fluid‐fixed museum specimens of squamate reptiles by using a forensic protocol developed for isolating DNA from human bones. The method yielded high quality nuclear DNA sequences from bones taken from 11 of 21 (52.4%) skeletonized or desiccated specimens, the oldest of which dated back to 1938, and 1 of 9 (11.1%) fluid‐preserved specimens, which was collected in 1957.     

A new cost‐effective and fast direct PCR protocol for insects based on PBS buffer

Insect DNA barcoding is a species identification technique used in biodiversity assessment and ecological studies. However, DNA extraction can result in the loss of up to 70% of DNA. Recent research has reported that direct PCR can overcome this issue. However, the success rates could still be improved, and tissues used for direct PCR could not be reused for further genetic studies. Here, we developed a direct PCR workflow that incorporates a 2‐min sample preparation in PBS‐buffer step for fast and effective universal insect species identification. The developed protocol achieved 100% success rates for amplification in six orders: Mantodea, Phasmatodea, Neuroptera, Odonata, Blattodea and Orthoptera. High and moderate success rates were obtained for five other species: Lepidoptera (97.3%), Coleoptera (93.8%), Diptera (90.5%), Hemiptera (81.8%) and Hymenoptera (75.0%). High‐quality sequencing data were also obtained from these amplifiable products, allowing confidence in species identification. The method was sensitive down to 1/4th of a 1‐mm fragment of leg or body and its success rates with oven‐dried, ethanol‐preserved, food, bat guano and museum specimens were 100%, 98.6%, 90.0%, 84.0% and 30.0%, respectively. In addition, the pre‐PCR solution (PBS with insect tissues) could be used for further DNA extraction if needed. The workflow will be beneficial in the fields of insect taxonomy and ecological studies due to its low cost, simplicity and applicability to highly degraded specimens.     

Shotgun microbial profiling of fossil remains

Millions to billions of DNA sequences can now be generated from ancient skeletal remains thanks to the massive throughput of next‐generation sequencing platforms. Except in cases of exceptional endogenous DNA preservation, most of the sequences isolated from fossil material do not originate from the specimen of interest, but instead reflect environmental organisms that colonized the specimen after death. Here, we characterize the microbial diversity recovered from seven c. 200‐ to 13 000‐year‐old horse bones collected from northern Siberia. We use a robust, taxonomy‐based assignment approach to identify the microorganisms present in ancient DNA extracts and quantify their relative abundance. Our results suggest that molecular preservation niches exist within ancient samples that can potentially be used to characterize the environments from which the remains are recovered. In addition, microbial community profiling of the seven specimens revealed site‐specific environmental signatures. These microbial communities appear to comprise mainly organisms that colonized the fossils recently. Our approach significantly extends the amount of useful data that can be recovered from ancient specimens using a shotgun sequencing approach. In future, it may be possible to correlate, for example, the accumulation of postmortem DNA damage with the presence and/or abundance of particular microbes.     

Index hopping on the Illumina HiseqX platform and its consequences for ancient DNA studies

The high‐throughput capacities of the Illumina sequencing platforms and the possibility to label samples individually have encouraged wide use of sample multiplexing. However, this practice results in read misassignment (usually <1%) across samples sequenced on the same lane. Alarmingly high rates of read misassignment of up to 10% were reported for lllumina sequencing machines with exclusion amplification chemistry. This may make use of these platforms prohibitive, particularly in studies that rely on low‐quantity and low‐quality samples, such as historical and archaeological specimens. Here, we use barcodes, short sequences that are ligated to both ends of the DNA insert, to directly quantify the rate of index hopping in 100‐year old museum‐preserved gorilla (Gorilla beringei) samples. Correcting for multiple sources of noise, we identify on average 0.470% of reads containing a hopped index. We show that sample‐specific quantity of misassigned reads depends on the number of reads that any given sample contributes to the total sequencing pool, so that samples with few sequenced reads receive the greatest proportion of misassigned reads. This particularly affects ancient DNA samples, as these frequently differ in their DNA quantity and endogenous content. Through simulations we show that even low rates of index hopping, as reported here, can lead to biases in ancient DNA studies when multiplexing samples with vastly different quantities of endogenous material.     

Nuclear DNA from a 180‐year‐old study skin reveals the phylogenetic position of the Kinglet Calyptura Calyptura cristata (Passeriformes: Tyrannides)

The Kinglet Calyptura Calyptura cristata is one of the most enigmatic bird species in South America, known only from specimens collected in the 19th century and a few recent observations. Knowledge of its biology is scanty and its systematic position is obscure. Traditionally, Calyptura was placed in the Cotingidae, but associated with genera that are now known to fall outside the Cotingidae. In an attempt to clarify its phylogenetic position, sequence data from four nuclear markers were obtained from a 180‐year‐old museum study skin of Calyptura, and incorporated into a comprehensive dataset of tyrant flycatchers, cotingas, manakins and allies. Our analyses demonstrate that Calyptura is most closely related to Platyrinchus and Neopipo and that these three genera constitute a deep branch in the clade containing the Rhynchocyclidae (tody‐tyrants and flatbills) and Tyrannidae (typical tyrant flycatchers). The Calyptura specimen is one of the oldest avian museum specimens from which a substantial amount of nuclear DNA sequence data have been obtained, and highlights the immense value of museum collections for DNA‐based phylogenetic studies.     

古DNA实时荧光定量PCR实验中标准品的制备

实时荧光定量PCR技术通过对PCR每一循环扩增产物的实时检测,可对模板的精确拷贝数进行绝对定量,从而用于古DNA实验中提取和扩增条件的比较和优化.本研究采用异硫氰酸胍碱裂解-SiO2吸附的方法,从采自黑龙江省的晚更新世斑鬣狗化石材料中提取得到了斑鬣狗线粒体基因组古DNA.经常规PCR扩增后,将纯化的扩增产物克隆到微生物体内使其大量复制,再用M13通用引物扩增出含少量外源DNA的古DNA目标片段,从而建立了适用于古DNA荧光定量PCR扩增的标准品的制备方法.经检测分析,运用该方法制备的标准品性质稳定,能够准确地指示反应体系中较为精确的古DNA模板拷贝数,从而反映古DNA的提取和扩增效率,用于比较并优化古DNA提取和扩增条件.