Incorporating historical museum specimens into molecular systematic and conservation genetics research

Uncovering the genomic and metagenomic research potential in old ethanol-preserved snakes

Natural history museum collections worldwide represent a tremendous resource of information on past and present biodiversity. Fish, reptiles, amphibians and many invertebrate collections have often been preserved in ethanol for decades or centuries and our knowledge on the genomic and metagenomic research potential of such material is limited. Here, we use ancient DNA protocols, combined with shotgun sequencing to test the molecular preservation in liver, skin and bone tissue from five old (1842 to 1964) museum specimens of the common garter snake (Thamnophis sirtalis). When mapping reads to a T. sirtalis reference genome, we find that the DNA molecules are highly damaged with short average sequence lengths (38-64 bp) and high C-T deamination, ranging from 9% to 21% at the first position. Despite this, the samples displayed relatively high endogenous DNA content, ranging from 26% to 56%, revealing that genome-scale analyses are indeed possible from all specimens and tissues included here. Of the three tested types of tissue, bone shows marginally but significantly higher DNA quality in these metrics. Though at least one of the snakes had been exposed to formalin, neither the concentration nor the quality of the obtained DNA was affected. Lastly, we demonstrate that these specimens display a diverse and tissue-specific microbial genetic profile, thus offering authentic metagenomic data despite being submerged in ethanol for many years. Our results emphasize that historical museum collections continue to offer an invaluable source of information in the era of genomics.

Species Separation within, and Preliminary Phylogeny for, the Leafhopper Genus Anoscopus with Particular Reference to the Putative British Endemic Anoscopus duffieldi (Hemiptera: Cicadellidae)

The subfamily Aphrodinae (Hemiptera: Cicadellidae) contains ~33 species in Europe within four genera. Species in two genera in particular, Aphrodes and Anoscopus, have proved to be difficult to distinguish morphologically. Our aim was to determine the status of the putative species Anoscopus duffieldi, found only on the RSPB Nature Reserve at Dungeness, Kent, a possible rare UK endemic. DNA from samples of all seven UK Anoscopus species (plus Anoscopusalpinus from the Czech Republic) were sequenced using parts of the mitochondrial cytochrome oxidase I and 16S rRNA genes. Bayesian inference phylogenies were created. Specimens of each species clustered into monophyletic groups, except for Anoscopusalbifrons, A. duffieldi and Anoscopuslimicola. Two A. albifrons specimens grouped with A. duffieldi repeatedly with strong support, and the remaining A. albifrons clustered within A. limicola. Genetic distances suggest that A. albifrons and A. limicola are a single interbreeding population (0% divergence), while A. albifrons and A. duffieldi diverged by only 0.28%. Shared haplotypes between A. albifrons, A. limicola and A. duffieldi strongly suggest interbreeding, although misidentification may also explain these topologies. However, all A. duffieldi clustered together in the trees. A conservative approach might be to treat A. duffieldi, until other evidence is forthcoming, as a possible endemic subspecies.

Dispersal versus vicariance in the Aegean: combining molecular and morphological phylogenies of eastern Mediterranean Dendarus (Coleoptera: Tenebrionidae) sheds new light on the phylogeography of the Aegean area

The Aegean archipelago, as an ‘evolutionary laboratory of nature’, is an ideal model for research in phylogeography. In this area, the darkling beetles of the genus Dendarus (distributed from Morocco to the Caucasus) exhibit a high level of diversity with 36 species, 27 of which are island endemics. However, their taxonomy is complex and unstable, having undergone continuous revision to address extensive morphological and ecological plasticity. Here, we examine the phylogenetic relationships of 23 species from Greece and Turkey, using mitochondrial and nuclear DNA sequences and 61 morphological characters, to unveil their phylogeny in the Aegean. This helps to clarify phylogeographic scenarios and historical processes that shaped the observed patterns. The analyses reveal 13 distinct lineages with several para- and polyphyletic cases that correspond to three major phylogroups [south/south-east Aegean (D. foraminosus complex, D. rhodius, D. sporadicus, D. wettsteini); central to north Aegean, Turkey and mainland Greece (D. crenulatus, D. moesiacus group, D. sinuatus complex, D. stygius) and mainland Greece (D. messenius, D. paganettii)], indicating the need for further taxonomic re-evaluation. Lineage topology and phylogeography suggest a spatial and temporal sequence of geographic isolation, following either a vicariant or a dispersal model coincident with major palaeogeographic separations in the Aegean.

Consideration of sample source for establishing reliable genetic microsatellite data from mammalian carnivore specimens held in natural history collections

Specimens from natural history collections (NHCs) are increasingly being used for genetic studies and can provide information on extinct populations, facilitate comparisons of historical and contemporary populations, produce baseline data before environmental changes, and elucidate patterns of change. Destructive sampling for DNA may be in disagreement with NHC goals of long-term care and maintenance. Differentiating quality among sample sources can direct destructive sampling to the source predicted to yield the highest quality DNA and most reliable data, potentially reducing damage to specimens, laboratory costs, and genotyping errors. We used the kit fox (Vulpes macrotis) as a model species and evaluated the quality and reliability of genetic data obtained from carnivoran specimens via three different sample sources: cranial bones, nasal bones, and toepads. We quantified variation in microsatellite amplification success and genotyping error rates and assessed the reliability of source-specific genic data. Toepads had the highest amplification success rates and lowest genotyping error rates. Shorter loci had higher amplification success and lower allelic dropout rates than longer loci. There were substantial differences in the reliability of resulting multilocus genotypes. Toepads produced the most reliable data, required the fewest replicates, and therefore, had the lowest costs to achieve reliable data. Our results demonstrate that the quality of DNA obtained from specimens varies by sample source and can inform NHCs when evaluating requests for destructive sampling. Our results suggest that prior to large-scale specimen sampling, researchers should conduct pilot studies to differentiate among source-specific data reliability, identify high performing loci, reduce costs of analyses, and minimize destructive sampling.

Evidence of introgressive hybridization between Stenella coeruleoalba and Delphinus delphis in the Greek Seas

Natural interspecific hybridization might be more important for the evolutionary history and speciation of animals than previously thought, considering several demographic and life history traits as well as habitat disturbance as factors that promote it. In this aspect, cetaceans comprise an interesting case in which the occurrence of sympatric species in mixed associations provides excellent opportunities for interspecific sexual interaction and the potential for hybridization. Here, we present evidence of natural hybridization for two cetacean species commonly occurring in the Greek Seas (Stenella coeruleoalba and Delphinus delphis), which naturally overlap in the Gulf of Corinth by analyzing highly resolving microsatellite DNA markers and mitochondrial DNA sequences in skin samples from 45 individuals of S. coeruleoalba, 12 D. delphis and three intermediate morphs. Employing several phylogenetic and population genetic approaches, we found 15 individuals that are potential hybrids including the three intermediate morphs, verifying the occurrence of natural hybridization between species of different genera. Their hybrids are fertile and able to reproduce not only with the other hybrids but also with each of the two-parental species. However, current evidence does not allow firm conclusions whether hybridization might constitute a step towards the generation of a new species and/or the swan song of an already existing species (i.e., D. delphis). Given that the focal species form mixed pods in several areas of Mediterranean, this study is an excellent opportunity to understand the mechanisms leading to hybridization in the context of gene flow and urges for the evaluation of the genetic status of common dolphins in the Mediterranean.

Assessment of alternatives to environmental toxic formalin for DNA conservation in biological specimens

One essential step of museum and clinical specimen preservation is immersion in a fixative fluid to prevent degradation. Formalin is the most largely used fixative, but its benefit is balanced with its toxic and carcinogenic status. Moreover, because formalin-fixation impairs nucleic acids recovery and quality, current museum wet collections and formalin-fixed, paraffin-embedded clinical samples do not represent optimal tanks of molecular information. Our study has been developed to compare formalin to two alternative fixatives (RCL2® and ethanol) in a context of molecular exploitation. Based on a unique protocol, we created mammalian fixed collections, simulated the impact of time on preservation using an artificial ageing treatment and followed the evolution of specimens' DNA quality. DNA extraction yield, purity, visual integrity and qualitative and quantitative ability to amplify the Cox1 gene were assessed. Our results show that both RCL2 and ethanol exhibit better performances than formalin. They do not impair DNA extraction yield, and more importantly, DNA alteration is delayed over the preservation step. The use of RCL2 or ethanol as fixative in biological collections may insure a better exploitation of the genetic resources they propose.

A reinvestigation of phylogeny and divergence times of the Ablepharus kitaibelii species complex (Sauria, Scincidae) based on mtDNA and nuDNA genes

Morphological and DNA data support that the East Mediterranean snake-eyed skink Ablepharus kitaibelii represents a species complex that includes four species A. kitaibelii, A. budaki, A. chernovi, and A. rueppellii, highlighting the need of its taxonomic reevaluation. Here, we used Bayesian and Maximum Likelihood methods to estimate the phylogenetic relationships of all members of the complex based on two mitochondrial (cyt b, 16S rRNA) and two nuclear markers (MC1R, and NKTR) and using Chalcides, Eumeces, and Eutropis as outgroups. The biogeographic history of the complex was also investigated through the application of several phylogeographic (BEAST) and biogeographic (BBM) analyses. Paleogeographic and paleoclimatic data were used to support the inferred phylogeographic patterns. The A. kitaibelli species complex exhibits high genetic diversity, revealing cases of hidden diversity and cases of non-monophyletic species such as A. kitaibelii and A. budaki. Our results indicate that A. pannonicus branches off first and a group that comprises specimens of A. kitaibelli and A. budaki from Kastelorizo Island group (southeast Greece) and southwest Turkey, respectively is differentiated from the rest A. kitaibelli and A. budaki populations and may represent a new species. The estimated divergence times place the origin of the complex in the Middle Miocene (∼16Mya) and the divergence of most currently recognized species in the Late Miocene. The inferred ancestral distribution suggests that the complex originated in Anatolia, supposing that several vicariance and dispersal events that are related with the formation of the Mid-Aegean Trench, the Anatolian Diagonal and the orogenesis of the mountain chains in southern and eastern Anatolia have led to current distribution pattern of A. kitaibelii species complex in the Balkans and Middle East.

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.

The role of integrative taxonomy in the conservation management of cryptic species: the taxonomic status of endangered earless dragons (Agamidae: Tympanocryptis) in the grasslands of Queensland, Australia

Molecular phylogenetics is increasingly highlighting the prevalence of cryptic species, where morphologically similar organisms have long independent evolutionary histories. When such cryptic species are known to be declining in numbers and are at risk of extinction due to a range of threatening processes, the disjunction between molecular systematics research and conservation policy becomes a significant problem. We investigate the taxonomic status of Tympanocryptis populations in Queensland, which have previously been assigned to T. tetraporophora, using three species delimitation approaches. The taxonomic uncertainties in this species-group are of particular importance in the Darling Downs Earless Dragon (T. cf. tetraporophora), which is ranked as an endangered 'species' of high priority for conservation by the Queensland Department of Environment and Heritage Protection. We undertook a morphological study, integrated with a comprehensive genetic study and species delimitation analyses, to investigate the species status of populations in the region. Phylogenetic analyses of two gene regions (mtDNA: ND2; nuclear: RAG1) revealed high levels of genetic divergence between populations, indicating isolation over long evolutionary time frames, and strongly supporting two independent evolutionary lineages in southeastern Queensland, from the Darling Downs, and a third in the Gulf Region of northern Queensland. Of the three species delimitation protocols used, we found integrative taxonomy the most applicable to this cryptic species complex. Our study demonstrates the utility of integrative taxonomy as a species delimitation approach in cryptic complexes of species with conservation significance, where limited numbers of specimens are available.

Museum genomics: low-cost and high-accuracy genetic data from historical specimens

Natural history collections are unparalleled repositories of geographical and temporal variation in faunal conditions. Molecular studies offer an opportunity to uncover much of this variation; however, genetic studies of historical museum specimens typically rely on extracting highly degraded and chemically modified DNA samples from skins, skulls or other dried samples. Despite this limitation, obtaining short fragments of DNA sequences using traditional PCR amplification of DNA has been the primary method for genetic study of historical specimens. Few laboratories have succeeded in obtaining genome-scale sequences from historical specimens and then only with considerable effort and cost. Here, we describe a low-cost approach using high-throughput next-generation sequencing to obtain reliable genome-scale sequence data from a traditionally preserved mammal skin and skull using a simple extraction protocol. We show that single-nucleotide polymorphisms (SNPs) from the genome sequences obtained independently from the skin and from the skull are highly repeatable compared to a reference genome.

A century of genetic change and metapopulation dynamics in the Galápagos warbler finches (Certhidea)

Populations that are connected by immigrants play an important role in evolutionary and conservation biology, yet we have little direct evidence of how such metapopulations change genetically over evolutionary time. We compared historic (1894-1906) to modern (1988-2006) genetic variation in 11 populations of warbler finches at 14 microsatellite loci. Although several lines of evidence suggest that Darwin's finches may be in decline, we found that the genetic diversity of warbler finches has not generally declined, and broad-scale patterns of variation remained similar over time. Contrary to expectations, inferred population sizes have generally increased over time (6-8%) as have immigration rates (8-16%), which may reflect a recent increase in the frequency and intensity of El Niño events. Individual island populations showed significant declines (18-19%) and also substantial gains (18-20%) in allelic richness over time. Changes in genetic diversity were correlated with changes in immigration rates, but did not correspond to population size or human disturbance. These results reflect the expected stabilizing properties of whole metapopulations over time. However, the dramatic and unpredictable changes observed in individual populations during this short time interval suggests that care should be taken when monitoring individual population fragments with snapshots of genetic variation.

Primers for identification of type and other archived specimens of Aphrodes leafhoppers (Hemiptera, Cicadellidae)

Primers were developed for leafhoppers of the genus Aphrodes amplifying 84-244 bp fragments of the mitochondrial cytochrome oxidase subunit I gene. DNA was extracted from legs of over 100-year-old archived museum specimens, amplified and sequenced. The fragments contain sufficient variation to unequivocally identify the different species. The majority of the analysed museum specimens, including three specimens of the syntype series for the UK endemic species A. aestuarina (Edwards 1908), were found to have been assigned to the wrong species. This work clearly underlines the need to validate museum specimens using molecular methods where identity is in doubt, based on reliable standards for species discrimination.