Reproducibility of ancient DNA sequences from extinct Pleistocene fauna

New populations of the black-flanked rock-wallaby (Petrogale lateralis) from the Little Sandy Desert and Murchison, Western Australia

During two independent fauna surveys, rock-wallaby (Petrogale) scats were recorded from caves located outside the current known Petrogale distribution. Scats collected from Desert Queen Baths (Little Sandy Desert, Western Australia, 2012), and the Barr Smith Range (Murchison, Western Australia, 2015) were genetically analysed and a follow-up motion camera survey confirmed an extant rock-wallaby population at Desert Queen Baths. The combination of sampling techniques overcame the detection difficulties associated with rare and cryptic taxa, and together were important in establishing the presence of Petrogale lateralis from regions where the species has been poorly documented. At both locations, P. lateralis scats were recorded from deep caves situated close to permanent water, reflecting the species’ physiological constraints in the arid zone. These records represent significant range extensions of a highly threatened macropod.

Metabarcoding of Environmental DNA Samples to Explore the Use of Uranium Mine Containment Ponds as a Water Source for Wildlife

Understanding how anthropogenic impacts on the landscape affect wildlife requires a knowledge of community assemblages. Species surveys are the first step in assessing community structure, and recent molecular applications such as metabarcoding and environmental DNA analyses have been proposed as an additional and complementary wildlife survey method. Here, we test eDNA metabarcoding as a survey tool to examine the potential use of uranium mine containment ponds as water sources by wildlife. We tested samples from surface water near mines and from one mine containment pond using two markers, 12S and 16S rRNA gene amplicons, to survey for vertebrate species. We recovered large numbers of sequence reads from taxa expected to be in the area and from less common or hard to observe taxa such as the tiger salamander and gray fox. Detection of these two species is of note because they were not observed in a previous species assessment, and tiger salamander DNA was found in the mine containment pond sample. We also found that sample concentration by centrifugation was a more efficient and more feasible method than filtration in these highly turbid surface waters. Ultimately, the use of eDNA metabarcoding could allow for a better understanding of the area’s overall biodiversity and community composition as well as aid current ecotoxicological risk assessment work.

DNA metabarcoding for diet analysis and biodiversity: A case study using the endangered Australian sea lion (Neophoca cinerea)

The analysis of apex predator diet has the ability to deliver valuable insights into ecosystem health, and the potential impacts a predator might have on commercially relevant species. The Australian sea lion (Neophoca cinerea) is an endemic apex predator and one of the world's most endangered pinnipeds. Given that prey availability is vital to the survival of top predators, this study set out to understand what dietary information DNA metabarcoding could yield from 36 sea lion scats collected across 1,500 km of its distribution in southwest Western Australia. A combination of PCR assays were designed to target a variety of potential sea lion prey, including mammals, fish, crustaceans, cephalopods, and birds. Over 1.2 million metabarcodes identified six classes from three phyla, together representing over 80 taxa. The results confirm that the Australian sea lion is a wide-ranging opportunistic predator that consumes an array of mainly demersal fauna. Further, the important commercial species Sepioteuthis australis (southern calamari squid) and Panulirus cygnus (western rock lobster) were detected, but were present in <25% of samples. Some of the taxa identified, such as fish, sharks and rays, clarify previous knowledge of sea lion prey, and some, such as eel taxa and two gastropod species, represent new dietary insights. Even with modest sample sizes, a spatial analysis of taxa and operational taxonomic units found within the scat shows significant differences in diet between many of the sample locations and identifies the primary taxa that are driving this variance. This study provides new insights into the diet of this endangered predator and confirms the efficacy of DNA metabarcoding of scat as a noninvasive tool to more broadly define regional biodiversity.

Assessing the feasibility of fly based surveillance of wildlife infectious diseases

Monitoring wildlife infectious agents requires acquiring samples suitable for analyses, which is often logistically demanding. A possible alternative to invasive or non-invasive sampling of wild-living vertebrates is the use of vertebrate material contained in invertebrates feeding on them, their feces, or their remains. Carrion flies have been shown to contain vertebrate DNA; here we investigate whether they might also be suitable for wildlife pathogen detection. We collected 498 flies in Taï National Park, Côte d’Ivoire, a tropical rainforest and examined them for adenoviruses (family Adenoviridae), whose DNA is frequently shed in feces of local mammals. Adenoviral DNA was detected in 6/142 mammal-positive flies. Phylogenetic analyses revealed that five of these sequences were closely related to sequences obtained from local non-human primates, while the sixth sequence was closely related to a murine adenovirus. Next-generation sequencing-based DNA-profiling of the meals of the respective flies identified putative hosts that were a good fit to those suggested by adenoviral sequence affinities. We conclude that, while characterizing the genetic diversity of wildlife infectious agents through fly-based monitoring may not be cost-efficient, this method could probably be used to detect the genetic material of wildlife infectious agents causing wildlife mass mortality in pristine areas.

High diversity of picornaviruses in rats from different continents revealed by deep sequencing

Outbreaks of zoonotic diseases in humans and livestock are not uncommon, and an important component in containment of such emerging viral diseases is rapid and reliable diagnostics. Such methods are often PCR-based and hence require the availability of sequence data from the pathogen. Rattus norvegicus (R. norvegicus) is a known reservoir for important zoonotic pathogens. Transmission may be direct via contact with the animal, for example, through exposure to its faecal matter, or indirectly mediated by arthropod vectors. Here we investigated the viral content in rat faecal matter (n=29) collected from two continents by analyzing 2.2 billion next-generation sequencing reads derived from both DNA and RNA. Among other virus families, we found sequences from members of the Picornaviridae to be abundant in the microbiome of all the samples. Here we describe the diversity of the picornavirus-like contigs including near-full-length genomes closely related to the Boone cardiovirus and Theiler's encephalomyelitis virus. From this study, we conclude that picornaviruses within R. norvegicus are more diverse than previously recognized. The virome of R. norvegicus should be investigated further to assess the full potential for zoonotic virus transmission.

In silico assessment of primers for eDNA studies using PrimerTree and application to characterize the biodiversity surrounding the Cuyahoga River

Analysis of environmental DNA (eDNA) enables the detection of species of interest from water and soil samples, typically using species-specific PCR. Here, we describe a method to characterize the biodiversity of a given environment by amplifying eDNA using primer pairs targeting a wide range of taxa and high-throughput sequencing for species identification. We tested this approach on 91 water samples of 40 mL collected along the Cuyahoga River (Ohio, USA). We amplified eDNA using 12 primer pairs targeting mammals, fish, amphibians, birds, bryophytes, arthropods, copepods, plants and several microorganism taxa and sequenced all PCR products simultaneously by high-throughput sequencing. Overall, we identified DNA sequences from 15 species of fish, 17 species of mammals, 8 species of birds, 15 species of arthropods, one turtle and one salamander. Interestingly, in addition to aquatic and semi-aquatic animals, we identified DNA from terrestrial species that live near the Cuyahoga River. We also identified DNA from one Asian carp species invasive to the Great Lakes but that had not been previously reported in the Cuyahoga River. Our study shows that analysis of eDNA extracted from small water samples using wide-range PCR amplification combined with high-throughput sequencing can provide a broad perspective on biological diversity.

Unbiased Characterization of Anopheles Mosquito Blood Meals by Targeted High-Throughput Sequencing

Understanding mosquito host choice is important for assessing vector competence or identifying disease reservoirs. Unfortunately, the availability of an unbiased method for comprehensively evaluating the composition of insect blood meals is very limited, as most current molecular assays only test for the presence of a few pre-selected species. These approaches also have limited ability to identify the presence of multiple mammalian hosts in a single blood meal. Here, we describe a novel high-throughput sequencing method that enables analysis of 96 mosquitoes simultaneously and provides a comprehensive and quantitative perspective on the composition of each blood meal. We validated in silico that universal primers targeting the mammalian mitochondrial 16S ribosomal RNA genes (16S rRNA) should amplify more than 95% of the mammalian 16S rRNA sequences present in the NCBI nucleotide database. We applied this method to 442 female Anopheles punctulatus s. l. mosquitoes collected in Papua New Guinea (PNG). While human (52.9%), dog (15.8%) and pig (29.2%) were the most common hosts identified in our study, we also detected DNA from mice, one marsupial species and two bat species. Our analyses also revealed that 16.3% of the mosquitoes fed on more than one host. Analysis of the human mitochondrial hypervariable region I in 102 human blood meals showed that 5 (4.9%) of the mosquitoes unambiguously fed on more than one person. Overall, analysis of PNG mosquitoes illustrates the potential of this approach to identify unsuspected hosts and characterize mixed blood meals, and shows how this approach can be adapted to evaluate inter-individual variations among human blood meals. Furthermore, this approach can be applied to any disease-transmitting arthropod and can be easily customized to investigate non-mammalian host sources.

Combined DNA, toxicological and heavy metal analyses provides an auditing toolkit to improve pharmacovigilance of traditional Chinese medicine (TCM)

Globally, there has been an increase in the use of herbal remedies including traditional Chinese medicine (TCM). There is a perception that products are natural, safe and effectively regulated, however, regulatory agencies are hampered by a lack of a toolkit to audit ingredient lists, adulterants and constituent active compounds. Here, for the first time, a multidisciplinary approach to assessing the molecular content of 26 TCMs is described. Next generation DNA sequencing is combined with toxicological and heavy metal screening by separation techniques and mass spectrometry (MS) to provide a comprehensive audit. Genetic analysis revealed that 50% of samples contained DNA of undeclared plant or animal taxa, including an endangered species of Panthera (snow leopard). In 50% of the TCMs, an undeclared pharmaceutical agent was detected including warfarin, dexamethasone, diclofenac, cyproheptadine and paracetamol. Mass spectrometry revealed heavy metals including arsenic, lead and cadmium, one with a level of arsenic >10 times the acceptable limit. The study showed 92% of the TCMs examined were found to have some form of contamination and/or substitution. This study demonstrates that a combination of molecular methodologies can provide an effective means by which to audit complementary and alternative medicines.

Investigating the zoonotic origin of the West African Ebola epidemic

The severe Ebola virus disease epidemic occurring in West Africa stems from a single zoonotic transmission event to a 2-year-old boy in Meliandou, Guinea. We investigated the zoonotic origins of the epidemic using wildlife surveys, interviews, and molecular analyses of bat and environmental samples. We found no evidence for a concurrent outbreak in larger wildlife. Exposure to fruit bats is common in the region, but the index case may have been infected by playing in a hollow tree housing a colony of insectivorous free-tailed bats (Mops condylurus). Bats in this family have previously been discussed as potential sources for Ebola virus outbreaks, and experimental data have shown that this species can survive experimental infection. These analyses expand the range of possible Ebola virus sources to include insectivorous bats and reiterate the importance of broader sampling efforts for understanding Ebola virus ecology.

A cautionary note on fecal sampling and molecular epidemiology in predatory wild great apes

Fecal samples are an important source of information on parasites (viruses, prokaryotes, or eukaryotes) infecting wild great apes. Molecular analysis of fecal samples has already been used for deciphering the origins of major human pathogens such as HIV-1 or Plasmodium falciparum. However, for apes that hunt (chimpanzees and bonobos), detection of parasite nucleic acids may reflect either true infection of the host of interest or ingestion of an infected prey, for example, another non-human primate. To determine the potential magnitude of this issue, we estimated the prevalence of prey DNA in fecal samples obtained from two wild chimpanzee communities. We observed values >15%, which are higher than or close to the fecal detection rates of many great ape parasites. Contamination of fecal samples with parasite DNA from dietary origin may therefore occasionally impact non-invasive epidemiological studies. This problem can be addressed (at least partially) by monitoring the presence of prey DNA.

From benchtop to desktop: important considerations when designing amplicon sequencing workflows

Amplicon sequencing has been the method of choice in many high-throughput DNA sequencing (HTS) applications. To date there has been a heavy focus on the means by which to analyse the burgeoning amount of data afforded by HTS. In contrast, there has been a distinct lack of attention paid to considerations surrounding the importance of sample preparation and the fidelity of library generation. No amount of high-end bioinformatics can compensate for poorly prepared samples and it is therefore imperative that careful attention is given to sample preparation and library generation within workflows, especially those involving multiple PCR steps. This paper redresses this imbalance by focusing on aspects pertaining to the benchtop within typical amplicon workflows: sample screening, the target region, and library generation. Empirical data is provided to illustrate the scope of the problem. Lastly, the impact of various data analysis parameters is also investigated in the context of how the data was initially generated. It is hoped this paper may serve to highlight the importance of pre-analysis workflows in achieving meaningful, future-proof data that can be analysed appropriately. As amplicon sequencing gains traction in a variety of diagnostic applications from forensics to environmental DNA (eDNA) it is paramount workflows and analytics are both fit for purpose.

Tag jumps illuminated--reducing sequence-to-sample misidentifications in metabarcoding studies

Metabarcoding of environmental samples on second-generation sequencing platforms has rapidly become a valuable tool for ecological studies. A fundamental assumption of this approach is the reliance on being able to track tagged amplicons back to the samples from which they originated. In this study, we address the problem of sequences in metabarcoding sequencing outputs with false combinations of used tags (tag jumps). Unless these sequences can be identified and excluded from downstream analyses, tag jumps creating sequences with false, but already used tag combinations, can cause incorrect assignment of sequences to samples and artificially inflate diversity. In this study, we document and investigate tag jumping in metabarcoding studies on Illumina sequencing platforms by amplifying mixed-template extracts obtained from bat droppings and leech gut contents with tagged generic arthropod and mammal primers, respectively. We found that an average of 2.6% and 2.1% of sequences had tag combinations, which could be explained by tag jumping in the leech and bat diet study, respectively. We suggest that tag jumping can happen during blunt-ending of pools of tagged amplicons during library build and as a consequence of chimera formation during bulk amplification of tagged amplicons during library index PCR. We argue that tag jumping and contamination between libraries represents a considerable challenge for Illumina-based metabarcoding studies, and suggest measures to avoid false assignment of tag jumping-derived sequences to samples.

The microbiome of New World vultures

Vultures are scavengers that fill a key ecosystem niche, in which they have evolved a remarkable tolerance to bacterial toxins in decaying meat. Here we report the first deep metagenomic analysis of the vulture microbiome. Through face and gut comparisons of 50 vultures representing two species, we demonstrate a remarkably conserved low diversity of gut microbial flora. The gut samples contained an average of 76 operational taxonomic units (OTUs) per specimen, compared with 528 OTUs on the facial skin. Clostridia and Fusobacteria, widely pathogenic to other vertebrates, dominate the vulture's gut microbiota. We reveal a likely faecal-oral-gut route for their origin. DNA of prey species detectable on facial swabs was completely degraded in the gut samples from most vultures, suggesting that the gastrointestinal tracts of vultures are extremely selective. Our findings show a strong adaption of vultures and their bacteria to their food source, exemplifying a specialized host-microbial alliance.

Targeted detection of mammalian species using carrion fly-derived DNA

DNA analysis from carrion flies (iDNA analysis) has recently been promoted as a powerful tool for cost- and time-efficient monitoring of wildlife. While originally applied to identify any mammalian species present in an area, it should also allow for targeted detection of species and individuals. Using carrion flies captured in the Taï National Park, Côte d'Ivoire, we assessed this possibility by (i) screening carrion fly DNA extracts with nonspecific and species-specific PCR systems, respectively, targeting mitochondrial DNA (mtDNA) fragments of any mammal or of Jentink's duiker (Cephalophus jentinki), three colobine monkeys (subfamily Colobinae) and sooty mangabey (Cercocebus atys); and (ii) genotyping carrion fly extracts containing sooty mangabey mtDNA. In comparison with the nonspecific PCR assay, the use of specific PCRs increased the frequency of detection of target species up to threefold. Detection rates partially reflected relative abundances of target species in the area. Amplification of seven microsatellite loci from carrion flies positive for sooty mangabey mtDNA yielded an average PCR success of 46%, showing that the identification of individuals is, to some extent, possible. Regression analysis of microsatellite PCR success and mtDNA concentration revealed that, among all carrion flies analysed for this study, 1% contained amounts of mammal mtDNA sufficient to attempt genotyping with potentially high success. We conclude that carrion fly-derived DNA analysis represents a promising tool for targeted monitoring of mammals in their natural habitat.

Metabarcoding avian diets at airports: implications for birdstrike hazard management planning

BACKGROUND

Wildlife collisions with aircraft cost the airline industry billions of dollars per annum and represent a public safety risk. Clearly, adapting aerodrome habitats to become less attractive to hazardous wildlife will reduce the incidence of collisions. Formulating effective habitat management strategies relies on accurate species identification of high-risk species. This can be successfully achieved for all strikes either through morphology and/or DNA-based identifications. Beyond species identification, dietary analysis of birdstrike gut contents can provide valuable intelligence for airport hazard management practices in regards to what food is attracting which species to aerodromes. Here, we present birdstrike identification and dietary data from Perth Airport, Western Australia, an aerodrome that saw approximately 140,000 aircraft movements in 2012. Next-generation high throughput DNA sequencing was employed to investigate 77 carcasses from 16 bird species collected over a 12-month period. Five DNA markers, which broadly characterize vertebrates, invertebrates and plants, were used to target three animal mitochondrial genes (12S rRNA, 16S rRNA, and COI) and a plastid gene (trnL) from DNA extracted from birdstrike carcass gastrointestinal tracts.

RESULTS

Over 151,000 DNA sequences were generated, filtered and analyzed by a fusion-tag amplicon sequencing approach. Across the 77 carcasses, the most commonly identified vertebrate was Mus musculus (house mouse). Acrididae (grasshoppers) was the most common invertebrate family identified, and Poaceae (grasses) the most commonly identified plant family. The DNA-based dietary data has the potential to provide some key insights into feeding ecologies within and around the aerodrome.

CONCLUSIONS

The data generated here, together with the methodological approach, will greatly assist in the development of hazard management plans and, in combination with existing observational studies, provide an improved way to monitor the effectiveness of mitigation strategies (for example, netting of water, grass type, insecticides and so on) at aerodromes. It is hoped that with the insights provided by dietary data, airports will be able to allocate financial resources to the areas that will achieve the best outcomes for birdstrike reduction.

Scrapheap challenge: a novel bulk-bone metabarcoding method to investigate ancient DNA in faunal assemblages

Highly fragmented and morphologically indistinct fossil bone is common in archaeological and paleontological deposits but unfortunately it is of little use in compiling faunal assemblages. The development of a cost-effective methodology to taxonomically identify bulk bone is therefore a key challenge. Here, an ancient DNA methodology using high-throughput sequencing is developed to survey and analyse thousands of archaeological bones from southwest Australia. Fossils were collectively ground together depending on which of fifteen stratigraphical layers they were excavated from. By generating fifteen synthetic blends of bulk bone powder, each corresponding to a chronologically distinct layer, samples could be collectively analysed in an efficient manner. A diverse range of taxa, including endemic, extirpated and hitherto unrecorded taxa, dating back to c.46,000 years BP was characterized. The method is a novel, cost-effective use for unidentifiable bone fragments and a powerful molecular tool for surveying fossils that otherwise end up on the taxonomic “scrapheap”.

Deep sequencing of plant and animal DNA contained within traditional Chinese medicines reveals legality issues and health safety concerns

Traditional Chinese medicine (TCM) has been practiced for thousands of years, but only within the last few decades has its use become more widespread outside of Asia. Concerns continue to be raised about the efficacy, legality, and safety of many popular complementary alternative medicines, including TCMs. Ingredients of some TCMs are known to include derivatives of endangered, trade-restricted species of plants and animals, and therefore contravene the Convention on International Trade in Endangered Species (CITES) legislation. Chromatographic studies have detected the presence of heavy metals and plant toxins within some TCMs, and there are numerous cases of adverse reactions. It is in the interests of both biodiversity conservation and public safety that techniques are developed to screen medicinals like TCMs. Targeting both the p-loop region of the plastid trnL gene and the mitochondrial 16S ribosomal RNA gene, over 49,000 amplicon sequence reads were generated from 15 TCM samples presented in the form of powders, tablets, capsules, bile flakes, and herbal teas. Here we show that second-generation, high-throughput sequencing (HTS) of DNA represents an effective means to genetically audit organic ingredients within complex TCMs. Comparison of DNA sequence data to reference databases revealed the presence of 68 different plant families and included genera, such as Ephedra and Asarum, that are potentially toxic. Similarly, animal families were identified that include genera that are classified as vulnerable, endangered, or critically endangered, including Asiatic black bear (Ursus thibetanus) and Saiga antelope (Saiga tatarica). Bovidae, Cervidae, and Bufonidae DNA were also detected in many of the TCM samples and were rarely declared on the product packaging. This study demonstrates that deep sequencing via HTS is an efficient and cost-effective way to audit highly processed TCM products and will assist in monitoring their legality and safety especially when plant reference databases become better established.

Pre-Clovis mastodon hunting 13,800 years ago at the Manis site, Washington

The tip of a projectile point made of mastodon bone is embedded in a rib of a single disarticulated mastodon at the Manis site in the state of Washington. Radiocarbon dating and DNA analysis show that the rib is associated with the other remains and dates to 13,800 years ago. Thus, osseous projectile points, common to the Beringian Upper Paleolithic and Clovis, were made and used during pre-Clovis times in North America. The Manis site, combined with evidence of mammoth hunting at sites in Wisconsin, provides evidence that people were hunting proboscideans at least two millennia before Clovis.

Blocking human contaminant DNA during PCR allows amplification of rare mammal species from sedimentary ancient DNA

Analyses of degraded DNA are typically hampered by contamination, especially when employing universal primers such as commonly used in environmental DNA studies. In addition to false-positive results, the amplification of contaminant DNA may cause false-negative results because of competition, or bias, during the PCR. In this study, we test the utility of human-specific blocking primers in mammal diversity analyses of ancient permafrost samples from Siberia. Using quantitative PCR (qPCR) on human and mammoth DNA, we first optimized the design and concentration of blocking primer in the PCR. Subsequently, 454 pyrosequencing of ancient permafrost samples amplified with and without the addition of blocking primer revealed that DNA sequences from a diversity of mammalian representatives of the Beringian megafauna were retrieved only when the blocking primer was added to the PCR. Notably, we observe the first retrieval of woolly rhinoceros (Coelodonta antiquitatis) DNA from ancient permafrost cores. In contrast, reactions without blocking primer resulted in complete dominance by human DNA sequences. These results demonstrate that in ancient environmental analyses, the PCR can be biased towards the amplification of contaminant sequences to such an extent that retrieval of the endogenous DNA is severely restricted. The application of blocking primers is a promising tool to avoid this bias and can greatly enhance the quantity and the diversity of the endogenous DNA sequences that are amplified.

Meta-barcoding of 'dirt' DNA from soil reflects vertebrate biodiversity

DNA molecules originating from animals and plants can be retrieved directly from sediments and have been used for reconstructing both contemporary and past ecosystems. However, the extent to which such 'dirt' DNA reflects taxonomic richness and structural diversity remains contentious. Here, we couple second generation high-throughput sequencing with 16S mitochondrial DNA (mtDNA) meta-barcoding, to explore the accuracy and sensitivity of 'dirt' DNA as an indicator of vertebrate diversity, from soil sampled at safari parks, zoological gardens and farms with known species compositions. PCR amplification was successful in the full pH range of the investigated soils (6.2 ± 0.2 to 8.3 ± 0.2), but inhibition was detected in extracts from soil of high organic content. DNA movement (leaching) through strata was evident in some sporadic cases and is influenced by soil texture and structure. We find that DNA from the soil surface reflects overall taxonomic richness and relative biomass of individual species. However, one species that was recently introduced was not detected. Furthermore, animal behaviour was shown to influence DNA deposition rates. The approach potentially provides a quick methodological alternative to classical ecological surveys of biodiversity, and most reliable results are obtained with spatial sample replicates, while relative amounts of soil processed per site is of less importance.

An in silico approach for the evaluation of DNA barcodes

BACKGROUND

DNA barcoding is a key tool for assessing biodiversity in both taxonomic and environmental studies. Essential features of barcodes include their applicability to a wide spectrum of taxa and their ability to identify even closely related species. Several DNA regions have been proposed as barcodes and the region selected strongly influences the output of a study. However, formal comparisons between barcodes remained limited until now. Here we present a standard method for evaluating barcode quality, based on the use of a new bioinformatic tool that performs in silico PCR over large databases. We illustrate this approach by comparing the taxonomic coverage and the resolution of several DNA regions already proposed for the barcoding of vertebrates. To assess the relationship between in silico and in vitro PCR, we also developed specific primers amplifying different species of Felidae, and we tested them using both kinds of PCR RESULTS: Tests on specific primers confirmed the correspondence between in silico and in vitro PCR. Nevertheless, results of in silico and in vitro PCRs can be somehow different, also because tuning PCR conditions can increase the performance of primers with limited taxonomic coverage. The in silico evaluation of DNA barcodes showed a strong variation of taxonomic coverage (i.e., universality): barcodes based on highly degenerated primers and those corresponding to the conserved region of the Cyt-b showed the highest coverage. As expected, longer barcodes had a better resolution than shorter ones, which are however more convenient for ecological studies analysing environmental samples.

CONCLUSIONS

In silico PCR could be used to improve the performance of a study, by allowing the preliminary comparison of several DNA regions in order to identify the most appropriate barcode depending on the study aims.

The use of coded PCR primers enables high-throughput sequencing of multiple homolog amplification products by 454 parallel sequencing

Background

The invention of the Genome Sequence 20™ DNA Sequencing System (454 parallel sequencing platform) has enabled the rapid and high-volume production of sequence data. Until now, however, individual emulsion PCR (emPCR) reactions and subsequent sequencing runs have been unable to combine template DNA from multiple individuals, as homologous sequences cannot be subsequently assigned to their original sources.

Methodology

We use conventional PCR with 5′-nucleotide tagged primers to generate homologous DNA amplification products from multiple specimens, followed by sequencing through the high-throughput Genome Sequence 20™ DNA Sequencing System (GS20, Roche/454 Life Sciences). Each DNA sequence is subsequently traced back to its individual source through 5′tag-analysis.

Conclusions

We demonstrate that this new approach enables the assignment of virtually all the generated DNA sequences to the correct source once sequencing anomalies are accounted for (miss-assignment rate<0.4%). Therefore, the method enables accurate sequencing and assignment of homologous DNA sequences from multiple sources in single high-throughput GS20 run. We observe a bias in the distribution of the differently tagged primers that is dependent on the 5′ nucleotide of the tag. In particular, primers 5′ labelled with a cytosine are heavily overrepresented among the final sequences, while those 5′ labelled with a thymine are strongly underrepresented. A weaker bias also exists with regards to the distribution of the sequences as sorted by the second nucleotide of the dinucleotide tags. As the results are based on a single GS20 run, the general applicability of the approach requires confirmation. However, our experiments demonstrate that 5′primer tagging is a useful method in which the sequencing power of the GS20 can be applied to PCR-based assays of multiple homologous PCR products. The new approach will be of value to a broad range of research areas, such as those of comparative genomics, complete mitochondrial analyses, population genetics, and phylogenetics.

Higher-order phylogeny of modern birds (Theropoda, Aves: Neornithes) based on comparative anatomy. II. Analysis and discussion

In recent years, avian systematics has been characterized by a diminished reliance on morphological cladistics of modern taxa, intensive palaeornithogical research stimulated by new discoveries and an inundation by analyses based on DNA sequences. Unfortunately, in contrast to significant insights into basal origins, the broad picture of neornithine phylogeny remains largely unresolved. Morphological studies have emphasized characters of use in palaeontological contexts. Molecular studies, following disillusionment with the pioneering, but non-cladistic, work of Sibley and Ahlquist, have differed markedly from each other and from morphological works in both methods and findings. Consequently, at the turn of the millennium, points of robust agreement among schools concerning higher-order neornithine phylogeny have been limited to the two basalmost and several mid-level, primary groups. This paper describes a phylogenetic (cladistic) analysis of 150 taxa of Neornithes, including exemplars from all non-passeriform families, and subordinal representatives of Passeriformes. Thirty-five outgroup taxa encompassing Crocodylia, predominately theropod Dinosauria, and selected Mesozoic birds were used to root the trees. Based on study of specimens and the literature, 2954 morphological characters were defined; these characters have been described in a companion work, approximately one-third of which were multistate (i.e. comprised at least three states), and states within more than one-half of these multistate characters were ordered for analysis. Complete heuristic searches using 10 000 random-addition replicates recovered a total solution set of 97 well-resolved, most-parsimonious trees (MPTs). The set of MPTs was confirmed by an expanded heuristic search based on 10 000 random-addition replicates and a full ratchet-augmented exploration to ascertain global optima. A strict consensus tree of MPTs included only six trichotomies, i.e. nodes differing topologically among MPTs. Bootstrapping (based on 10 000 replicates) percentages and ratchet-minimized support (Bremer) indices indicated most nodes to be robust. Several fossil Neornithes (e.g. Dinornithiformes, Aepyornithiformes) were placed within the ingroup a posteriori either through unconstrained, heursitic searches based on the complete matrix augmented by these taxa separately or using backbone-constraints. Analysis confirmed the topology among outgroup Theropoda and achieved robust resolution at virtually all levels of the Neornithes. Findings included monophyly of the palaeognathous birds, comprising the sister taxa Tinamiformes and ratites, respectively, and the Anseriformes and Galliformes as monophyletic sister-groups, together forming the sister-group to other Neornithes exclusive of the Palaeognathae (Neoaves). Noteworthy inferences include: (i) the sister-group to remaining Neoaves comprises a diversity of marine and wading birds; (ii) Podicipedidae are the sister-group of Gaviidae, and not closely related to the Phoenicopteridae, as recently suggested; (iii) the traditional Pelecaniformes, including the shoebill (Balaeniceps rex) as sister-taxon to other members, are monophyletic; (iv) traditional Ciconiiformes are monophyletic; (v) Strigiformes and Falconiformes are sister-groups; (vi) Cathartidae is the sister-group of the remaining Falconiformes; (vii) Ralliformes (Rallidae and Heliornithidae) are the sister-group to the monophyletic Charadriiformes, with the traditionally composed Gruiformes and Turniciformes (Turnicidae and Mesitornithidae) sequentially paraphyletic to the entire foregoing clade; (viii) Opisthocomus hoazin is the sister-taxon to the Cuculiformes (including the Musophagidae); (ix) traditional Caprimulgiformes are monophyletic and the sister-group of the Apodiformes; (x) Trogoniformes are the sister-group of Coliiformes; (xi) Coraciiformes, Piciformes and Passeriformes are mutually monophyletic and closely related; and (xii) the Galbulae are retained within the Piciformes. Unresolved portions of the Neornithes (nodes having more than one most-parsimonious solution) comprised three parts of the tree: (a) several interfamilial nodes within the Charadriiformes; (b) a trichotomy comprising the (i) Psittaciformes, (ii) Columbiformes and (iii) Trogonomorphae (Trogoniformes, Coliiformes) + Passerimorphae (Coraciiformes, Piciformes, Passeriformes); and (c) a trichotomy comprising the Coraciiformes, Piciformes and Passeriformes. The remaining polytomies were among outgroups, although several of the highest-order nodes were only marginally supported; however, the majority of nodes were resolved and met or surpassed conventional standards of support. Quantitative comparisons with alternative hypotheses, examination of highly supportive and diagnostic characters for higher taxa, correspondences with prior studies, complementarity and philosophical differences with palaeontological phylogenetics, promises and challenges of palaeogeography and calibration of evolutionary rates of birds, and classes of promising evidence and future directions of study are reviewed. Homology, as applied to avian examples of apparent homologues, is considered in terms of recent theory, and a revised annotated classification of higher-order taxa of Neornithes and other closely related Theropoda is proposed. (c) 2007 The Linnean Society of London, Zoological Journal of the Linnean Society, 2007, 149, 1-95.

Ancient genes of Saccharomyces cerevisiae

Amber is a plant resin mainly produced by coniferous trees that, after entrapping a variety of living beings, was subjected to a process of fossilization until it turned into yellowish, translucent stones. It is also one of the best sources of ancient DNA on which to perform studies on evolution. Here a method for the sterilization of amber that allows reliable ancient DNA extraction with no actual DNA contamination is described. Working with insects taken from amber, it was possible to amplify theATP9,PGU1andrRNA18Sancient genes ofSaccharomyces cerevisiaecorresponding to samples from the Miocene and Oligocene. After comparison of the current genes with their ancient (up to 35–40 million years) counterparts it was concluded that essential genes such asrRNA18Sare highly conserved and that even normal ‘house-keeping’ genes, such asPGU1, are strikingly conserved along the millions of years thatS. cerevisiaehas evolved.

Mitochondrial cytochrome b of the Lyakhov mammoth (Proboscidea, Mammalia): new data and phylogenetic analyses of Elephantidae

The phylogenetic relationships between recent Elephantidae (Proboscidea, Mammalia), that is to say extant elephants (Asian and African) and extinct woolly mammoth, have remained unclear to date. The prevailing morphological scheme (mammoth grouped with Asian elephant) is either supported or questioned by the molecular results. Recently, the monophyly of woolly mammoths on mitochondrial grounds has been demonstrated (Thomas, et al., 2000), but it conflicts with previous studies (Barriel et al., 1999; Derenko et al., 1997). Here, we report the partial sequencing of two mitochondrial genes: 128 bp of 12S rDNA and 561 bp of cytochrome b for the Lyakhov mammoth, a 49,000-year-old Siberian individual. We use the most comprehensive sample of mammoth (11 sequences) to determine whether the sequences achieved by former studies were congruent or not. The monophyly of a major subset of mammoths sequences (including ours) is recovered. Such a result is assumed to be a good criterion for ascertaining the origin of ancient DNA. Our sequence is incongruent with that of Yang et al. (1996), though obtained for the same individual. As far as the latter sequence is concerned, a contamination by non-identified exogenous DNA is suspected. The robustness and reliability of the sister group relation between Mammuthus primigenius and Loxodonta africana are examined: down-weighting saturated substitutions has no impact on the topology; analyzing data partitions proves that the support of this clade can be assigned to the most conservative phylogenetic signal; insufficient taxonomic and/or characters sampling contributed to former discordant conclusions. We therefore assume the monophyly of "real mammoth sequences" and the (Mammuthus, Loxodonta) clade.

Evolution of endogenous retrovirus-like elements of the woolly mammoth (Mammuthus primigenius) and its relatives

Endogenous retrovirus-like elements characterizable by a leucine tRNA primer (ERV-Ls) are reiterated genomic sequences known to be widespread in mammals, including humans. They may have arisen from an ancestral foamy virus-like element by successful germ line infection followed by copy number expansion. However, among mammals, only primates and rodents have thus far exhibited high copy number amplification and sequence diversification. Conventionally, empirical studies of proviral amplification and diversification have been limited to extant species, but taxa having good Quaternary fossil records could potentially be investigated using the techniques of "ancient" DNA research. To examine evolutionary parameters of ERV-Ls across both time and taxa, we characterized this proviral class in the extinct woolly mammoth (Mammuthus primigenius) and living elephants, as well as extant members of the larger clade to which they belong (Uranotheria, a group containing proboscideans, sirenians, hyraxes, and their extinct relatives). Ungulates and carnivores previously analyzed demonstrated low copy numbers of ERV-L sequences, and thus it was expected that uranotheres should as well. Here, we show that all uranothere taxa exhibit unexpectedly numerous and diverse ERV-L sequence complements, indicating active expansion within this group of lineages. Selection is the most parsimonious explanation for observed differences in ERV-L distribution and frequency, with relative success being reflected in the persistence of certain elements over a variety of sampled time depths (as can be observed by comparing sequences from fossil and extant elephantid samples).

A molecular phylogeny of two extinct sloths

Xenarthra (Edentata) is an extremely diverse mammalian order whose modern representatives are the armadillos, anteaters, and sloths. The phylogeny of these groups is poorly resolved. This is particularly true for the sloths (phyllophagans), originally a large and diverse group now reduced to two genera in two different families. Both morphological analyses and molecular analyses of rDNA genes of living and extinct sloths have been used with limited success to elucidate their phylogeny. In an attempt to clarify relationships among the sloths, DNA was extracted and mitochondrial cytochrome b gene sequences were determined from representatives of two extinct groups of sloths (Mylodontidae and Megatheriidae), their two living relatives (two-toed sloths [Megalonychidae], three-toed sloths [Bradypodidae]), anteaters and armadillos. A consistent feature of the latter two species was the nuclear copies of cytochrome b gene sequences. Several methods of phylogenetic reconstruction were applied to the sequences determined, and the results were compared with 12S rDNA sequences obtained in previous studies. The cytochrome b gene exhibited a phylogenetic resolving power similar to that of the 12S rDNA sequences. When both data sets were combined, they tended to support the grouping of two-toed sloths with mylodontids and three-toed sloths with megatheriids. The results strengthen the view that the two families of living sloths adapted independently to an arboreal life-style.

Protein preservation and DNA retrieval from ancient tissues

The retrieval of DNA from fossils remains controversial. To substantiate claims of DNA recovery, one needs additional information on the preservation of other molecules within the same sample. Flash pyrolysis with GC and MS was used to assess the quality of protein preservation in 11 archaeological and paleontological remains, some of which have yielded ancient DNA sequences authenticated via a number of criteria and some of which have consistently failed to yield any meaningful DNA. Several samples, including the Neanderthal-type specimen from which DNA sequences were recently reported, yielded abundant pyrolysis products assigned to 2,5-diketopiperazines of proline-containing dipeptides. The relative amounts of these products provide a good index of the amount of peptide hydrolysis and DNA preservation. Of these samples, four stem from arctic or subarctic regions, emphasizing the importance of cooler temperatures for the preservation of macromolecules. Flash pyrolysis with GC and MS offers a rapid and effective method for assessing fossils for the possibility of DNA preservation.

What does palaeontology contribute to systematics in a molecular world?

Palaeontology provides the only direct record for morphological and genetic change through time and uniquely contributes to systematics in two ways: by providing access to denser taxon sampling than is otherwise possible and by dating divergence times. Claims that ancient DNA has survived millions of years in certain fossils suggested the possibility that palaeontology could contribute directly to molecular systematic studies. Unfortunately, none of the supposed geologically ancient DNA records stands up to detailed scrutiny and fossils therefore contribute primarily through the morphological information they preserve. Denser taxon sampling can improve the accuracy of phylogenetic estimates primarily through allowing better discrimination of homoplasy from homology. This in turn leads to more accurate hypotheses of character transformation. Denser taxon sampling also offers the opportunity for more accurate rooting, since more characters can be polarized by reference to a stem-group taxon than to an extant sister-group taxon. Missing data can be a problem for fossils, but is not crippling. Finally the temporal order of clade appearances in the fossil record can provide ancillary evidence for selecting a working phylogeny from among a number of equally most parsimonious cladograms.

New uses for old DNA

Several years have elapsed since the last report of million-year-old DNA, coinciding with increased standards for experimental procedures in ancient DNA research. Whereas many earlier studies are now regarded as erroneous, the recent successful characterisation of Neanderthal DNA has set new standards for the field. Researchers continue to find new ways to exploit preserved genetic information in studies of more recent remains, widening the utility of ancient DNA.

Flightlessness and phylogeny amongst endemic rails (Aves:Rallidae) of the New Zealand region

The phylogenetic relationships of a number of flightless and volant rails have been investigated using mtDNA sequence data. The third domain of the small ribosomal subunit (12S) has been sequenced for 22 taxa, and part of the 5' end of the cytochrome-b gene has been sequenced for 12 taxa. Additional sequences were obtained from outgroup taxa, two species of jacana, sarus crane, spur-winged plover and kagu. Extinct rails were investigated using DNA extracted from subfossil bones, and in cases where fresh material could not be obtained from other extant taxa, feathers and museum skins were used as sources of DNA. Phylogenetic trees produced from these data have topologies that are, in general, consistent with data from DNA-DNA hybridization studies and recent interpretations based on morphology. Gallinula chloropus moorhen) groups basally with Fulica (coots), Amaurornis (= Megacrex) ineptus falls within the Gallirallus/Rallus group, and Gallinula (= Porphyrula) martinica is basal to Porphyrio (swamphens) and should probably be placed in that genus. Subspecies of Porphyrio porphyrio are paraphyletic with respect to Porphyrio mantelli (takahe). The Northern Hemisphere Rallus aquaticus is basal to the south-western Pacific Rallus (or Gallirallus) group. The flightless Rallus philippensis dieffenbachii is close to Rallus modestus and distinct from the volant Rallus philippensis, and is evidently a separate species. Porzana (crakes) appears to be more closely associated with Porphyrio than Rallus. Deep relationships among the rails remain poorly resolved. Rhynochetus jubatus (kagu) is closer to the cranes than the rails in this analysis. Genetic distances between flightless rails and their volant counterparts varied considerably with observed 12S sequence distances, ranging from 0.3% (Porphyrio porphyrio melanotus and P. mantelli mantelli) to 7.6% (Rallus modestus and Rallus philippensis). This may be taken as an indication of the rapidity with which flightlessness can evolve, and of the persistence of flightless taxa. Genetic data supported the notion that flightless taxa were independently derived, sometimes from similar colonizing ancestors. The morphology of flightless rails is apparently frequently dominated by evolutionary parallelism although similarity of external appearance is not an indication of the extent of genetic divergence. In some cases taxa that are genetically close are morphologically distinct from one another (e.g. Rallus (philippensis) dieffenbachii and R. modestus), whilst some morphologically similar taxa are evidently independently derived (e.g. Porphyio mantelli hochstetteri and P.m. mantelli).

Problems of reproducibility--does geologically ancient DNA survive in amber-preserved insects?

Apparently ancient DNA has been reported from amber-preserved insects many millions of years old. Rigorous attempts to reproduce these DNA sequences from amber- and copal-preserved bees and flies have failed to detect any authentic ancient insect DNA. Lack of reproducibility suggests that DNA does not survive over millions of years even in amber, the most promising of fossil environments.

A blind testing design for authenticating ancient DNA sequences

Reproducibility is a serious concern among researchers of ancient DNA. We designed a blind testing procedure to evaluate laboratory accuracy and authenticity of ancient DNA obtained from closely related extant and extinct species. Soft tissue and bones of fossil and contemporary museum proboscideans were collected and identified based on morphology by one researcher, and other researchers carried out DNA testing on the samples, which were assigned anonymous numbers. DNA extracted using three principal isolation methods served as template in PCR amplifications of a segment of the cytochrome b gene (mitochondrial genome), and the PCR product was directly sequenced and analyzed. The results show that such a blind testing design performed in one laboratory, when coupled with phylogenetic analysis, can nonarbitrarily test the consistency and reliability of ancient DNA results. Such reproducible results obtained from the blind testing can increase confidence in the authenticity of ancient sequences obtained from postmortem specimens and avoid bias in phylogenetic analysis. A blind testing design may be applicable as an alternative to confirm ancient DNA results in one laboratory when independent testing by two laboratories is not available.

DNA damage and DNA sequence retrieval from ancient tissues

Gas chromatography/mass spectrometry (GC/MS) was used to determine the amounts of eight oxidative base modifications in DNA extracted from 11 specimens of bones and soft tissues, ranging in age from 40 to >50 000 years. Among the compounds assayed hydantoin derivatives of pyrimidines were quantitatively dominant. From five of the specimens endogenous ancient DNA sequences could be amplified by PCR. The DNA from these specimens contained substantially lower amounts of hydantoins than the six specimens from which no DNA could be amplified. Other types of damage, e.g. oxidation products of purines, did not correlate with the inability to retrieve DNA sequences. Furthermore, all samples with low amounts of damage and from which DNA could be amplified stemmed from regions where low temperatures have prevailed throughout the burial period of the specimens.