Recharacterization of ancient DNA miscoding lesions: insights in the era of sequencing-by-synthesis

Compacted hair in broken teeth reveals dietary prey of historic lions

With recent advances, nuclear genome data for phylogenomic analyses can now be sequenced from minuscule quantities of DNA1 and from specimens that are more than a million years old.2 DNA analysis from hair is a well-established approach3 widely used in forensic science4 and wildlife conservation.5 Hair samples can be effectively decontaminated6 and can be used to identify the mammalian species from which the hair was shed.7,8 We aimed to use advances optimized for degraded DNA to systematically identify dietary prey species from hair compacted in the teeth of two Tsavo lions that lived during the 1890s in Kenya (see description of samples in the STAR Methods and Patterson9 and Kerbis Peterhans and Gnoske10 for background on the Tsavo "man-eaters"). Analysis of hair DNA identified giraffe, human, oryx, waterbuck, wildebeest, and zebra as prey and also identified hair that originated from lion. DNA preservation allowed for analyses of complete mitogenome profiles of zebra, giraffe, and lion. Giraffe mitogenomes are phylogeographically partitioned, and we found that the lions ate at least two individuals that belong to a subspecies of Masai giraffe (Giraffa tippelskirchi tippelskirchi) typically found in southeast Kenya. The lion mitogenome from a hair sample was identical to the Tsavo lion endogenous mitogenome and most closely matched other East African lions from Kenya and Tanzania. Our approach enables a better understanding of the hunting behaviors, diets, and ecology of historical individuals, populations, and species and holds promise for extinct populations and species.

Nucleic Acids Persistence-Benefits and Limitations in Forensic Genetics

The analysis of genetic material may be the only way to identify an unknown person or solve a criminal case. Often, the conditions in which the genetic material was found determine the choice of the analytical method. Hence, it is extremely important to understand the influence of various factors, both external and internal, on genetic material. The review presents information on DNA and RNA persistence, depending on the chemical and physical factors affecting the genetic material integrity. One of the factors taken into account is the time elapsing to genetic material recovery. Temperature can both preserve the genetic material or lead to its rapid degradation. Radiation, aquatic environments, and various types of chemical and physical factors also affect the genetic material quality. The substances used during the forensic process, i.e., for biological trace visualization or maceration, are also discussed. Proper analysis of genetic material degradation can help determine the post-mortem interval (PMI) or time since deposition (TsD), which may play a key role in criminal cases.

Ancient DNA analysis from epoxy resin Biodur®-embedded bones

For microscopic investigation, archaeological bone samples are often embedded in Biodur^® epoxy resin. This study wants to test whether it is possible to extract DNA suitable for PCR amplification from this sample type. For eight individuals a set of samples - each consisting of a Biodur-embedded femur sample, a native femur sample and a control sample of different anatomical origin - were submitted to organic DNA extraction. The extraction success was tested by autosomal short tandem repeat amplification. Seven out of eight Biodur-embedded femur samples revealed successful amplification results. If Biodur-embedded bone material exists from earlier microscopic investigations, our results encourage the use of this sample type as a source for genetic research.

Using in silico predicted ancestral genomes to improve the efficiency of paleogenome reconstruction

Paleogenomics is the nascent discipline concerned with sequencing and analysis of genome-scale information from historic, ancient, and even extinct samples. While once inconceivable due to the challenges of DNA damage, contamination, and the technical limitations of PCR-based Sanger sequencing, following the dawn of the second-generation sequencing revolution, it has rapidly become a reality. However, a significant challenge facing ancient DNA studies on extinct species is the lack of closely related reference genomes against which to map the sequencing reads from ancient samples. Although bioinformatic efforts to improve the assemblies have focused mainly in mapping algorithms, in this article we explore the potential of an alternative approach, namely using reconstructed ancestral genome as reference for mapping DNA sequences of ancient samples. Specifically, we present a preliminary proof of concept for a general framework and demonstrate how under certain evolutionary divergence thresholds, considerable mapping improvements can be easily obtained.

Were domestic camelids present on the prehispanic South American agricultural frontier? An ancient DNA study

The southern boundary of prehispanic farming in South America occurs in central Mendoza Province, Argentina at approximately 34 degrees south latitude. Archaeological evidence of farming includes the recovery of macrobotanical remains of cultigens and isotopic chemistry of human bone. Since the 1990s, archaeologists have also hypothesized that the llama (Lama glama), a domesticated South American camelid, was also herded near the southern boundary of prehispanic farming. The remains of a wild congeneric camelid, the guanaco (Lama guanicoe), however, are common in archaeological sites throughout Mendoza Province. It is difficult to distinguish bones of the domestic llama from wild guanaco in terms of osteological morphology, and therefore, claims that llama were in geographic areas where guanaco were also present based on osteometric analysis alone remain equivocal. A recent study, for example, claimed that twenty-five percent of the camelid remains from the high elevation Andes site of Laguna del Diamante S4 were identified based on osteometric evidence as domestic llama, but guanaco are also a likely candidate since the two species overlap in size. We test the hypothesis that domesticated camelids occurred in prehispanic, southern Mendoza through analysis of ancient DNA. We generated whole mitochondrial genome datasets from 41 samples from southern Mendoza late Holocene archaeological sites, located between 450 and 3400 meters above sea level (masl). All camelid samples from those sites were identified as guanaco; thus, we have no evidence to support the hypothesis that the domestic llama occurred in prehispanic southern Mendoza.

Single-strand DNA library preparation improves sequencing of formalin-fixed and paraffin-embedded (FFPE) cancer DNA

DNA derived from formalin-fixed and paraffin-embedded (FFPE) tissue has been a challenge to large-scale genomic sequencing, due to its low quality and quantities. Improved techniques enabling the genome-wide analysis of FFPE material would be of great value, both from a research and clinical perspective.

Comparing a single-strand DNA library preparation method originally developed for ancient DNA to conventional protocols using double-stranded DNA derived from FFPE material we obtain on average 900-fold more library molecules and improved sequence complexity from as little as 5 ng input DNA. FFPE DNA is highly fragmented, usually below 100bp, and up to 60% of reads start after or end prior to adenine residues, suggesting that crosslinks predominate at adenine residues. Similar to ancient DNA, C > T substitutions are slightly increased with maximum rates up to 3% at the ends of molecules. In whole exome sequencing of single-strand libraries from lung, breast, colorectal, prostate and skin cancers we identify known cancer mutations. In summary, we show that single-strand library preparation enables genomic sequencing, even from low amounts of degraded FFPE DNA. This method provides a clear advantage both in research and clinical settings, where FFPE material (e.g. from biopsies) often is the only source of DNA available. Improving the genetic characterization that can be performed on conventional archived FFPE tissue, the single-strand library preparation allows scarce samples to be used in personalized medicine and enables larger sample sizes in future sequencing studies.

Imputation approach for deducing a complete mitogenome sequence from low-depth-coverage next-generation sequencing data: application to ancient remains from the Moon Pyramid, Mexico

It is considered that more than 15 depths of coverage are necessary for next-generation sequencing (NGS) data to obtain reliable complete nucleotide sequences of the mitogenome. However, it is difficult to satisfy this requirement for all nucleotide positions because of problems obtaining a uniform depth of coverage for poorly preserved materials. Thus, we propose an imputation approach that allows a complete mitogenome sequence to be deduced from low-depth-coverage NGS data. We used different types of mitogenome data files as panels for imputation: a worldwide panel comprising all the major haplogroups, a worldwide panel comprising sequences belonging to the estimated haplogroup alone, a panel comprising sequences from the population most closely related to an individual under investigation, and a panel comprising sequences belonging to the estimated haplogroup from the population most closely related to an individual under investigation. The number of missing nucleotides was drastically reduced in all the panels, but the contents obtained by imputation were quite different among the panels. The efficiency of the imputation method differed according to the panels used. The missing nucleotides were most credibly imputed using sequences of the estimated haplogroup from the population most closely related to the individual under investigation as a panel.

Analysis of Allelic Imbalance in Rice Hybrids Under Water Stress and Association of Asymmetrically Expressed Genes with Drought-Response QTLs

BACKGROUND

Information on the effect of stress on the allele-specific expression (ASE) profile of rice hybrids is limited. More so, the association of allelically imbalanced genes to important traits is yet to be understood. Here we assessed allelic imbalance (AI) in the heterozygote state of rice under non- and water-stress treatments and determined association of asymmetrically expressed genes with grain yield (GY) under drought stress by in-silico co-localization analysis and selective genotyping. The genotypes IR64, Apo and their F1 hybrid (IR64 × Apo) were grown under normal and water-limiting conditions. We sequenced the total RNA transcripts for all genotypes then reconstructed the two chromosomes in the heterozygote.

RESULTS

We are able to estimate the transcript abundance of and the differential expression (DE) between the two parent-specific alleles in the rice hybrids. The magnitude and direction of AI are classified into two categories: (1) symmetrical or biallelic and (2) asymmetrical. The latter can be further classified as either IR64- or Apo-favoring gene. Analysis showed that in the hybrids grown under non-stress conditions, 179 and 183 favor Apo- and IR64-specific alleles, respectively. Hence, the number of IR64- and Apo-favoring genes is relatively equal. Under water-stress conditions, 179 and 255 favor Apo- and IR64-specific alleles, respectively, indicating that the number of allelically imbalanced genes is skewed towards IR64. This is nearly 40-60 % preference for Apo and IR64 alleles, respectively, to the hybrid transcriptome. We also observed genes which exhibit allele preference switching when exposed to water-stress conditions. Results of in-silico co-localization procedure and selective genotyping of Apo/IR64 F3:5 progenies revealed significant association of several asymmetrically expressed genes with GY under drought stress conditions.

CONCLUSION

Our data suggest that water stress skews AI on a genome-wide scale towards the IR64 allele, the cross-specific maternal allele. Several asymmetrically expressed genes are strongly associated with GY under drought stress which may shed hints that genes associated with important traits are allelically imbalanced. Our approach of integrating hybrid expression analysis and QTL mapping analysis may be an efficient strategy for shortlisting candidate genes for gene discovery.

Considering DNA damage when interpreting mtDNA heteroplasmy in deep sequencing data

Resolution of mitochondrial (mt) DNA heteroplasmy is now possible when applying a massively parallel sequencing (MPS) approach, including minor components down to 1%. However, reporting thresholds and interpretation criteria will need to be established for calling heteroplasmic variants that address a number of important topics, one of which is DNA damage. We assessed the impact of increasing amounts of DNA damage on the interpretation of minor component sequence variants in the mtDNA control region, including low-level mixed sites. A passive approach was used to evaluate the impact of storage conditions, and an active approach was employed to accelerate the process of hydrolytic damage (for example, replication errors associated with depurination events). The patterns of damage were compared and assessed in relation to damage typically encountered in poor quality samples. As expected, the number of miscoding lesions increased as conditions worsened. Single nucleotide polymorphisms (SNPs) associated with miscoding lesions were indistinguishable from innate heteroplasmy and were most often observed as 1-2% of the total sequencing reads. Numerous examples of miscoding lesions above 2% were identified, including two complete changes in the nucleotide sequence, presenting a challenge when assessing the placement of reporting thresholds for heteroplasmy. To mitigate the impact, replication of miscoding lesions was not observed in stored samples, and was rarely seen in data associated with accelerated hydrolysis. In addition, a significant decrease in the expected transition:transversion ratio was observed, providing a useful tool for predicting the presence of damage-induced lesions. The results of this study directly impact MPS analysis of minor sequence variants from poorly preserved DNA extracts, and when biological samples have been exposed to agents that induce DNA damage. These findings are particularly relevant to clinical and forensic investigations.

Detecting hybridization using ancient DNA

It is well established that related species hybridize and that this can have varied but significant effects on speciation and environmental adaptation. It should therefore come as no surprise that hybridization is not limited to species that are alive today. In the last several decades, advances in technologies for recovering and sequencing DNA from fossil remains have enabled the assembly of high-coverage genome sequences for a growing diversity of organisms, including many that are extinct. Thanks to the development of new statistical approaches for detecting and quantifying admixture from genomic data, genomes from extinct populations have proven useful both in revealing previously unknown hybridization events and informing the study of hybridization between living organisms. Here, we review some of the key recent statistical innovations for detecting ancient hybridization using genomewide sequence data and discuss how these innovations have revised our understanding of human evolutionary history.

Molecular analysis of ancient caries

An 84 base pair sequence of the Streptococcus mutans virulence factor, known as dextranase, has been obtained from 10 individuals from the Bronze Age to the Modern Era in Europe and from before and after the colonization in America. Modern samples show four polymorphic sites that have not been found in the ancient samples studied so far. The nucleotide and haplotype diversity of this region have increased over time, which could be reflecting the footprint of a population expansion. While this segment has apparently evolved according to neutral evolution, we have been able to detect one site that is under positive selection pressure both in present and past populations. This study is a first step to study the evolution of this microorganism, analysed using direct evidence obtained from ancient remains.

Gourds and squashes (Cucurbita spp.) adapted to megafaunal extinction and ecological anachronism through domestication

The genus Cucurbita (squashes, pumpkins, gourds) contains numerous domesticated lineages with ancient New World origins. It was broadly distributed in the past but has declined to the point that several of the crops' progenitor species are scarce or unknown in the wild. We hypothesize that Holocene ecological shifts and megafaunal extinctions severely impacted wild Cucurbita, whereas their domestic counterparts adapted to changing conditions via symbiosis with human cultivators. First, we used high-throughput sequencing to analyze complete plastid genomes of 91 total Cucurbita samples, comprising ancient (n = 19), modern wild (n = 30), and modern domestic (n = 42) taxa. This analysis demonstrates independent domestication in eastern North America, evidence of a previously unknown pathway to domestication in northeastern Mexico, and broad archaeological distributions of taxa currently unknown in the wild. Further, sequence similarity between distant wild populations suggests recent fragmentation. Collectively, these results point to wild-type declines coinciding with widespread domestication. Second, we hypothesize that the disappearance of large herbivores struck a critical ecological blow against wild Cucurbita, and we take initial steps to consider this hypothesis through cross-mammal analyses of bitter taste receptor gene repertoires. Directly, megafauna consumed Cucurbita fruits and dispersed their seeds; wild Cucurbita were likely left without mutualistic dispersal partners in the Holocene because they are unpalatable to smaller surviving mammals with more bitter taste receptor genes. Indirectly, megafauna maintained mosaic-like landscapes ideal for Cucurbita, and vegetative changes following the megafaunal extinctions likely crowded out their disturbed-ground niche. Thus, anthropogenic landscapes provided favorable growth habitats and willing dispersal partners in the wake of ecological upheaval.

Examining the uncertain origin and management role of martens on Prince of Wales Island, Alaska

Conservation biologists are generally united in efforts to curtail the spread of non-native species globally. However, the colonization history of a species is not always certain, and whether a species is considered non-native or native depends on the conservation benchmark. Such ambiguities have led to inconsistent management. Within the Tongass National Forest of Alaska, the status of American marten (Martes americana) on the largest, most biologically diverse and deforested island, Prince of Wales (POW), is unclear. Ten martens were released to POW in the early 1930s, and it was generally believed to be the founding event, although this has been questioned. The uncertainty surrounding when and how martens colonized POW complicates management, especially because martens were selected as a design species for the Tongass. To explore the history of martens of POW we reviewed other plausible routes of colonization; genetically and isotopically analyzed putative marten fossils deposited in the late Pleistocene and early Holocene to verify marten occupancy of POW; and used contemporary genetic data from martens on POW and the mainland in coalescent simulations to identify the probable source of the present-day marten population on POW. We found evidence for multiple routes of colonization by forest-associated mammals beginning in the Holocene, which were likely used by American martens to naturally colonize POW. Although we cannot rule out human-assisted movement of martens by Alaskan Natives or fur trappers, we suggest that martens be managed for persistence on POW. More generally, our findings illustrate the difficulty of labeling species as non-native or native, even when genetic and paleo-ecological data are available, and support the notion that community resilience or species invasiveness should be prioritized when making management decisions rather than more subjective and less certain conservation benchmarks.

Reconstructing ancient genomes and epigenomes

Research involving ancient DNA (aDNA) has experienced a true technological revolution in recent years through advances in the recovery of aDNA and, particularly, through applications of high-throughput sequencing. Formerly restricted to the analysis of only limited amounts of genetic information, aDNA studies have now progressed to whole-genome sequencing for an increasing number of ancient individuals and extinct species, as well as to epigenomic characterization. Such advances have enabled the sequencing of specimens of up to 1 million years old, which, owing to their extensive DNA damage and contamination, were previously not amenable to genetic analyses. In this Review, we discuss these varied technical challenges and solutions for sequencing ancient genomes and epigenomes.

Ancient pathogen genomics: insights into timing and adaptation

Disease is a major cause of natural selection affecting human evolution, whether through a sudden pandemic or persistent morbidity and mortality. Recent contributions in the field of ancient pathogen genomics have advanced our understanding of the antiquity and nature of human-pathogen interactions through time. Technical advancements have facilitated the recovery, enrichment, and high-throughput sequencing of pathogen and parasite DNA from archived and archaeological remains. These time-stamped genomes are crucial for calibrating molecular clocks to infer the timing of evolutionary events, while providing finer-grain resolution to phylogenetic reconstructions and complex biogeographical patterns. Additionally, genome scale data allow better identification of substitutions linked to adaptations of the pathogen to their human hosts. As methodology continues to improve, ancient genomes of humans and their diverse microbiomes from a range of eras and archaeological contexts will enable population-level ancient analyses in the near future and a better understanding of their co-evolutionary history.

Reconstructing the colonization history of lost wolf lineages by the analysis of the mitochondrial genome

The grey wolves (Canis lupus) originally inhabited major parts of the Northern hemisphere, but many local populations became extinct. Two lineages of wolves in Japan, namely, Japanese or Honshu (C. l. hodophilax) and Ezo or Hokkaido (C. l. hattai) wolves, rapidly went extinct between 100 and 120years ago. Here we analyse the complete mitochondrial genome sequences from ancient specimens and reconstruct the colonization history of the two extinct subspecies. We show a unique status of Japanese wolves in wolf phylogeny, suggesting their long time separation from other grey wolf populations. Japanese wolves appeared to have colonized the Japanese archipelago in the Late Pleistocene (ca. 25,000-125,000years ago). By contrast, Ezo wolves, which are clearly separated from Japanese wolves in phylogeny, are likely to have arrived at Japan relatively recently (<14,000years ago). Interestingly, their colonization history to Japan tallies well with the dynamics of wolf populations in Europe and America during the last several millennia. Our analyses suggest that at least several thousands of wolves once inhabited in the Japanese archipelago. Our analyses also show that an enigmatic clade of domestic dogs is likely to have originated from rare admixture events between male dogs and female Japanese wolves.

Genotyping of ancient Mycobacterium tuberculosis strains reveals historic genetic diversity

The evolutionary history of the Mycobacterium tuberculosis complex (MTBC) has previously been studied by analysis of sequence diversity in extant strains, but not addressed by direct examination of strain genotypes in archaeological remains. Here, we use ancient DNA sequencing to type 11 single nucleotide polymorphisms and two large sequence polymorphisms in the MTBC strains present in 10 archaeological samples from skeletons from Britain and Europe dating to the second–nineteenth centuries AD. The results enable us to assign the strains to groupings and lineages recognized in the extant MTBC. We show that at least during the eighteenth–nineteenth centuries AD, strains of M. tuberculosis belonging to different genetic groups were present in Britain at the same time, possibly even at a single location, and we present evidence for a mixed infection in at least one individual. Our study shows that ancient DNA typing applied to multiple samples can provide sufficiently detailed information to contribute to both archaeological and evolutionary knowledge of the history of tuberculosis.

Emulsion PCR-coupled target enrichment: an effective fishing method for high-throughput sequencing of poorly preserved ancient DNA

Due to the difficulties in deep sequencing, high-throughput sequencing of ancient DNA has been limited to exceptionally well-preserved ancient materials. The primary factor is microbial attack popularly observed in the buried materials, and it causes drastic increase in relative ratio of microbial DNA in the extracted DNA. We present a unified strategy in which emulsion PCR is coupled with target enrichment followed by next-generation sequencing. The method made it possible to obtain efficiently non-duplicated reads mapped to target sequences of interest, and this can achieve deep and reliable sequencing of ancient DNA from typical materials, even though poorly preserved.

Further evaluation of the efficacy of contamination removal from bone surfaces

Studies of low copy number (LCN) and degraded DNA are prone to contamination from exogenous DNA sources that in some cases out-compete endogenous DNA in PCR amplification, thus leading to false positives and/or aberrant results. Particularly problematic is contamination that is inadvertently deposited on the surfaces of bones through direct handling. Whereas some previous studies have shown that contamination removal is possible by subjecting samples to sodium hypochlorite prior to DNA extraction, others caution that such treatment can destroy a majority of the molecules endogenous to the sample. To further explore this topic, we experimentally contaminated ancient northern fur seal (Callorhinus ursinus) ribs with human DNA and treated them with sodium hypochlorite to remove that contamination. Our findings are consistent with previous studies that found sodium hypochlorite to be highly efficient (~81-99%) at contamination removal; however, there emerged no treatment capable of removing 100% of the contamination across all of the experiments. Moreover, the ability to estimate the degree of damage to endogenous northern fur seal molecules was compromised due to the inherent variability of preserved mtDNA across the bones, and the presence of co-extracted PCR inhibitors.

Ancient DNA damage

Under favorable conditions DNA can survive for thousands of years in the remains of dead organisms. The DNA extracted from such remains is invariably degraded to a small average size by processes that at least partly involve depurination. It also contains large amounts of deaminated cytosine residues that are accumulated toward the ends of the molecules, as well as several other lesions that are less well characterized.

Evaluating the efficacy of various thermo-stable polymerases against co-extracted PCR inhibitors in ancient DNA samples

DNA from ancient and forensic specimens is often co-extracted with unknown amounts of unknown PCR inhibitors, which can lead to underestimated DNA concentrations, allelic drop-out, and/or false-negative results. It is not surprising, in this case, that numerous methods have been developed to remove PCR inhibitors or subdue their effects. One simple and cost effective approach could be the adoption of a polymerase that overcomes or is less affected by PCR inhibitors. In this study, nine different polymerases were evaluated for their efficacy against PCR inhibitors co-extracted with DNA from 63 ancient salmon vertebrae. These samples were excavated from two archeological sites located at the Dionisio Point locality on the northern end of Galiano Island in coastal southwestern British Columbia, Canada and date to 700-1000 and 1300-1500 years before present. Previously, DNA extracts from samples studied from this locality were determined to be largely inhibited to PCR amplification. In the present study, Omni Klentaq LA (DNA Polymerase Technology, Inc.) outperformed the other 8 polymerases in two measures: (1) its success in genetic species identification of these vertebrae, and (2) its ability to amplify an ancient DNA positive control when spiked with a volume of potentially inhibited extract from the vertebrae.

Analysis of 454 sequencing error rate, error sources, and artifact recombination for detection of Low-frequency drug resistance mutations in HIV-1 DNA

Background

454 sequencing technology is a promising approach for characterizing HIV-1 populations and for identifying low frequency mutations. The utility of 454 technology for determining allele frequencies and linkage associations in HIV infected individuals has not been extensively investigated. We evaluated the performance of 454 sequencing for characterizing HIV populations with defined allele frequencies.

Results

We constructed two HIV-1 RT clones. Clone A was a wild type sequence. Clone B was identical to clone A except it contained 13 introduced drug resistant mutations. The clones were mixed at ratios ranging from 1% to 50% and were amplified by standard PCR conditions and by PCR conditions aimed at reducing PCR-based recombination. The products were sequenced using 454 pyrosequencing. Sequence analysis from standard PCR amplification revealed that 14% of all sequencing reads from a sample with a 50:50 mixture of wild type and mutant DNA were recombinants. The majority of the recombinants were the result of a single crossover event which can happen during PCR when the DNA polymerase terminates synthesis prematurely. The incompletely extended template then competes for primer sites in subsequent rounds of PCR. Although less often, a spectrum of other distinct crossover patterns was also detected. In addition, we observed point mutation errors ranging from 0.01% to 1.0% per base as well as indel (insertion and deletion) errors ranging from 0.02% to nearly 50%. The point errors (single nucleotide substitution errors) were mainly introduced during PCR while indels were the result of pyrosequencing. We then used new PCR conditions designed to reduce PCR-based recombination. Using these new conditions, the frequency of recombination was reduced 27-fold. The new conditions had no effect on point mutation errors. We found that 454 pyrosequencing was capable of identifying minority HIV-1 mutations at frequencies down to 0.1% at some nucleotide positions.

Conclusion

Standard PCR amplification results in a high frequency of PCR-introduced recombination precluding its use for linkage analysis of HIV populations using 454 pyrosequencing. We designed a new PCR protocol that resulted in a much lower recombination frequency and provided a powerful technique for linkage analysis and haplotype determination in HIV-1 populations. Our analyses of 454 sequencing results also demonstrated that at some specific HIV-1 drug resistant sites, mutations can reliably be detected at frequencies down to 0.1%.

PyroTRF-ID: a novel bioinformatics methodology for the affiliation of terminal-restriction fragments using 16S rRNA gene pyrosequencing data

Background

In molecular microbial ecology, massive sequencing is gradually replacing classical fingerprinting techniques such as terminal-restriction fragment length polymorphism (T-RFLP) combined with cloning-sequencing for the characterization of microbiomes. Here, a bioinformatics methodology for pyrosequencing-based T-RF identification (PyroTRF-ID) was developed to combine pyrosequencing and T-RFLP approaches for the description of microbial communities. The strength of this methodology relies on the identification of T-RFs by comparison of experimental and digital T-RFLP profiles obtained from the same samples. DNA extracts were subjected to amplification of the 16S rRNA gene pool, T-RFLP with the HaeIII restriction enzyme, 454 tag encoded FLX amplicon pyrosequencing, and PyroTRF-ID analysis. Digital T-RFLP profiles were generated from the denoised full pyrosequencing datasets, and the sequences contributing to each digital T-RF were classified to taxonomic bins using the Greengenes reference database. The method was tested both on bacterial communities found in chloroethene-contaminated groundwater samples and in aerobic granular sludge biofilms originating from wastewater treatment systems.

Results

PyroTRF-ID was efficient for high-throughput mapping and digital T-RFLP profiling of pyrosequencing datasets. After denoising, a dataset comprising ca. 10′000 reads of 300 to 500 bp was typically processed within ca. 20 minutes on a high-performance computing cluster, running on a Linux-related CentOS 5.5 operating system, enabling parallel processing of multiple samples. Both digital and experimental T-RFLP profiles were aligned with maximum cross-correlation coefficients of 0.71 and 0.92 for high- and low-complexity environments, respectively. On average, 63±18% of all experimental T-RFs (30 to 93 peaks per sample) were affiliated to phylotypes.

Conclusions

PyroTRF-ID profits from complementary advantages of pyrosequencing and T-RFLP and is particularly adapted for optimizing laboratory and computational efforts to describe microbial communities and their dynamics in any biological system. The high resolution of the microbial community composition is provided by pyrosequencing, which can be performed on a restricted set of selected samples, whereas T-RFLP enables simultaneous fingerprinting of numerous samples at relatively low cost and is especially adapted for routine analysis and follow-up of microbial communities on the long run.

High potential for using DNA from ancient herring bones to inform modern fisheries management and conservation

Pacific herring (Clupea pallasi) are an abundant and important component of the coastal ecosystems for the west coast of North America. Current Canadian federal herring management assumes five regional herring populations in British Columbia with a high degree of exchange between units, and few distinct local populations within them. Indigenous traditional knowledge and historic sources, however, suggest that locally adapted, distinct regional herring populations may have been more prevalent in the past. Within the last century, the combined effects of commercial fishing and other anthropogenic factors have resulted in severe declines of herring populations, with contemporary populations potentially reflecting only the remnants of a previously more abundant and genetically diverse metapopulation. Through the analysis of 85 archaeological herring bones, this study attempted to reconstruct the genetic diversity and population structure of ancient herring populations using three different marker systems (mitochondrial DNA (mtDNA), microsatellites and SNPs). A high success rate (91%) of DNA recovery was obtained from the extremely small herring bone samples (often <10 mg). The ancient herring mtDNA revealed high haplotype diversity comparable to modern populations, although population discrimination was not possible due to the limited power of the mtDNA marker. Ancient microsatellite diversity was also similar to modern samples, but the data quality was compromised by large allele drop-out and stuttering. In contrast, SNPs were found to have low error rates with no evidence for deviations from Hardy-Weinberg equilibrium, and simulations indicated high power to detect genetic differentiation if loci under selection are used. This study demonstrates that SNPs may be the most effective and feasible approach to survey genetic population structure in ancient remains, and further efforts should be made to screen for high differentiation markers.This study provides the much needed foundation for wider scale studies on temporal genetic variation in herring, with important implications for herring fisheries management, Aboriginal title rights and herring conservation.

Reconstructing the history of a fragmented and heavily exploited red deer population using ancient and contemporary DNA

Background

Red deer (Cervus elaphus) have been an important human resource for millennia, experiencing intensive human influence through habitat alterations, hunting and translocation of animals. In this study we investigate a time series of ancient and contemporary DNA from Norwegian red deer spanning about 7,000 years. Our main aim was to investigate how increasing agricultural land use, hunting pressure and possibly human mediated translocation of animals have affected the genetic diversity on a long-term scale.

Results

We obtained mtDNA (D-loop) sequences from 73 ancient specimens. These show higher genetic diversity in ancient compared to extant samples, with the highest diversity preceding the onset of agricultural intensification in the Early Iron Age. Using standard diversity indices, Bayesian skyline plot and approximate Bayesian computation, we detected a population reduction which was more prolonged than, but not as severe as, historic documents indicate. There are signs of substantial changes in haplotype frequencies primarily due to loss of haplotypes through genetic drift. There is no indication of human mediated translocations into the Norwegian population. All the Norwegian sequences show a western European origin, from which the Norwegian lineage diverged approximately 15,000 years ago.

Conclusions

Our results provide direct insight into the effects of increasing habitat fragmentation and human hunting pressure on genetic diversity and structure of red deer populations. They also shed light on the northward post-glacial colonisation process of red deer in Europe and suggest increased precision in inferring past demographic events when including both ancient and contemporary DNA.

Heat degradation of eukaryotic and bacterial DNA: an experimental model for paleomicrobiology

Background

Theoretical models suggest that DNA degradation would sharply limit the PCR-based detection of both eukaryotic and prokaryotic DNA within ancient specimens. However, the relative extent of decay of eukaryote and prokaryote DNA over time is a matter of debate. In this study, the murine macrophage cell line J774, alone or infected with Mycobacterium smegmatis bacteria, were killed after exposure to 90°C dry heat for intervals ranging from 1 to 48 h in order to compare eukaryotic cells, extracellular bacteria and intracellular bacteria. The sizes of the resulting mycobacterial rpoB and murine rpb2 homologous gene fragments were then determined by real-time PCR and fluorescent probing.

Findings

The cycle threshold (Ct) values of PCR-amplified DNA fragments from J774 cells and the M. smegmatis negative controls (without heat exposure) varied from 26–33 for the J774 rpb2 gene fragments and from 24–29 for M. smegmatis rpoB fragments. After 90°C dry heat incubation for up to 48 h, the Ct values of test samples increased relative to those of the controls for each amplicon size. For each dry heat exposure time, the Ct values of the 146-149-bp fragments were lower than those of 746-747-bp fragments. During the 4- to 24-h dry heat incubation, the non-infected J774 cell DNA was degraded into 597-bp rpb2 fragments. After 48 h, however, only 450-bp rpb2 fragments of both non-infected and infected J774 cells could be amplified. In contrast, the 746-bp rpoB fragments of M. smegmatis DNA could be amplified after the 48-h dry heat exposure in all experiments. Infected and non-infected J774 cell DNA was degraded more rapidly than M. smegmatis DNA after dry heat exposure (ANOVA test, p < 0.05).

Conclusion

In this study, mycobacterial DNA was more resistant to dry-heat stress than eukaryotic DNA. Therefore, the detection of large, experimental, ancient mycobacterial DNA fragments is a suitable approach for paleomicrobiological studies.

Earliest Mexican Turkeys (Meleagris gallopavo) in the Maya Region: implications for pre-Hispanic animal trade and the timing of turkey domestication

Late Preclassic (300 BC-AD 100) turkey remains identified at the archaeological site of El Mirador (Petén, Guatemala) represent the earliest evidence of the Mexican turkey (Meleagris gallopavo) in the ancient Maya world. Archaeological, zooarchaeological, and ancient DNA evidence combine to confirm the identification and context. The natural pre-Hispanic range of the Mexican turkey does not extend south of central Mexico, making the species non-local to the Maya area where another species, the ocellated turkey (Meleagris ocellata), is indigenous. Prior to this discovery, the earliest evidence of M. gallopavo in the Maya area dated to approximately one thousand years later. The El Mirador specimens therefore represent previously unrecorded Preclassic exchange of animals from northern Mesoamerica to the Maya cultural region. As the earliest evidence of M. gallopavo found outside its natural geographic range, the El Mirador turkeys also represent the earliest indirect evidence for Mesoamerican turkey rearing or domestication. The presence of male, female and sub-adult turkeys, and reduced flight morphology further suggests that the El Mirador turkeys were raised in captivity. This supports an argument for the origins of turkey husbandry or at least captive rearing in the Preclassic.

Genome-wide identification of Brassica napus microRNAs and their targets in response to cadmium

MicroRNAs (miRNAs) are a distinct class of small RNAs in plants that not only regulate biological processes but also regulate response to environmental stresses. The toxic heavy metal cadmium (Cd) induces expression of several miRNAs in rapeseed (Brassica napus), but it is not known on a genome-wide scale how the expression of miRNAs and their target genes, is regulated by Cd. In this study, four small RNA libraries and four degradome libraries were constructed from Cd-treated and non-Cd-treated roots and shoots of B. napus seedlings. Using high-throughput sequencing, the study identified 84 conserved and non-conserved miRNAs (belonging to 37 miRNA families) from Cd-treated and non-treated B. napus, including 19 miRNA members that were not identified before. Some of the miRNAs were validated by RNA gel blotting. Most of the identified miRNAs were found to be differentially expressed in roots/shoots or regulated by Cd exposure. The study simultaneously identified 802 targets for the 37 (24 conserved and 13 non-conserved) miRNA families, from which there are 200, 537, and 65 targets, belonging to categories I, II, and III, respectively. In category I alone, many novel targets for miRNAs were identified and shown to be involved in plant response to Cd.

Regulation of leaf morphology by microRNA394 and its target LEAF CURLING RESPONSIVENESS

The present study identified Arabidopsis miR394 and its target, an F-box (SKP1-Cullin/CDC53-F-box) gene At1g27340 (here referred to as LEAF CURLING RESPONSIVENESS, LCR), for regulation of leaf curling-related morphology. The loss-of-function lcr mutants exhibit pleiotropic defects with semi-dwarfism, altered leaf shape and a shorter stem. Overexpression of an miR394-resistant version of LCR under the 35S promoter (35S:m5LCR) and target mimicry MIM394 resulted in a curled-down leaf defect. Conversely, transgenic plants overexpressing 35S:MIR394a/b display a curled-up leaf phenotype. Detailed analyses show that there is a certain level of LCR that is optimal for leaf morphology, but lower or higher levels lead to abnormal leaf development, indicating that expression of miR394 in the leaf lamina is necessary for proper leaf morphology. Because the phytohormone auxin plays a crucial role in leaf morphogenesis and patterning, the DR5-GUS reporter gene was used to monitor the auxin response. We show that DR5 expression patterns in lcr and 35S::m5LCR plants were significantly different from those in the wild type. Also, overexpression of LCR in 35S::m5LCR plants drastically decreased the expression of the auxin-responsive genes IAA3, AXR3 and IAMT1, whereas increased expression of the genes was found in 35S::MIR394a plants. These results indicate that miR394 and its target LCR are involved in the regulation of leaf development.

Improving ancient DNA read mapping against modern reference genomes

Background

Next-Generation Sequencing has revolutionized our approach to ancient DNA (aDNA) research, by providing complete genomic sequences of ancient individuals and extinct species. However, the recovery of genetic material from long-dead organisms is still complicated by a number of issues, including post-mortem DNA damage and high levels of environmental contamination. Together with error profiles specific to the type of sequencing platforms used, these specificities could limit our ability to map sequencing reads against modern reference genomes and therefore limit our ability to identify endogenous ancient reads, reducing the efficiency of shotgun sequencing aDNA.

Results

In this study, we compare different computational methods for improving the accuracy and sensitivity of aDNA sequence identification, based on shotgun sequencing reads recovered from Pleistocene horse extracts using Illumina GAIIx and Helicos Heliscope platforms. We show that the performance of the Burrows Wheeler Aligner (BWA), that has been developed for mapping of undamaged sequencing reads using platforms with low rates of indel-types of sequencing errors, can be employed at acceptable run-times by modifying default parameters in a platform-specific manner. We also examine if trimming likely damaged positions at read ends can increase the recovery of genuine aDNA fragments and if accurate identification of human contamination can be achieved using a strategy previously suggested based on best hit filtering. We show that combining our different mapping and filtering approaches can increase the number of high-quality endogenous hits recovered by up to 33%.

Conclusions

We have shown that Illumina and Helicos sequences recovered from aDNA extracts could not be aligned to modern reference genomes with the same efficiency unless mapping parameters are optimized for the specific types of errors generated by these platforms and by post-mortem DNA damage. Our findings have important implications for future aDNA research, as we define mapping guidelines that improve our ability to identify genuine aDNA sequences, which in turn could improve the genotyping accuracy of ancient specimens. Our framework provides a significant improvement to the standard procedures used for characterizing ancient genomes, which is challenged by contamination and often low amounts of DNA material.

Ancient human DNA

The contribution of palaeogenetic data to the study of various aspects of hominin biology and evolution has been significant, and has the potential to increase substantially with the widespread implementation of next generation sequencing techniques. Here we discuss the present state-of-the-art of ancient human DNA analysis and the characteristics of hominin aDNA that make sequence validation particularly complex. A brief overview of the development of anthropological palaeogenetic analysis is given to illustrate the technical challenges motivating recent technological advancements.

Case study: ancient sloth DNA recovered from hairs preserved in paleofeces

Ancient hair, which has proved to be an excellent source of well-preserved ancient DNA, is often preserved in paleofeces. Here, we separate and wash hair shafts preserved in a paleofecal specimen believed to be from a Darwin's ground sloth, Mylodon darwinii. After extracting DNA from the recovered and cleaned hair using a protocol optimized for DNA extraction from keratinous substrates, we amplify 12S and 16S rDNA sequences from the DNA extract. As expected, the recovered sequences most closely match previously published sequences of M. darwinii. Our results demonstrate that hair preserved in paleofeces, even from temperate cave environments, is an effective source of ancient DNA.

Quantitative real-time PCR in aDNA research

Quantitative real-time PCR (qPCR) is a technique that is widely used in the field of ancient DNA (aDNA). Quantitative PCR can be used to optimize aDNA extraction methodologies, to detect PCR inhibition, and to quantify aDNA libraries for use in high-throughput sequencing. In this chapter, we outline factors that need to be considered when developing efficient SYBR Green qPCR assays. We describe how to setup qPCR standards of known copy number and provide some useful tips regarding interpretation of qPCR data generated from aDNA templates.

Capillary electrophoresis of human mtDNA control region sequences from highly degraded samples using short mtDNA amplicons

The forensic applications of mtDNA sequencing center primarily on samples that are either highly degraded or contain little or no nuclear DNA, since the testing of these sample types is often unsuccessful with more widely used nuclear STR profiling assays. In these cases, sequence data from the noncoding mtDNA control region are targeted due to its high variability. However, the ease of authentic DNA recovery and the strategy used for recovery depend strictly on the quality of the sample. In this chapter, we will cover mitochondrial DNA sequencing procedures for short mtDNA amplicons which range in size from 100 to 350 bp. Generally speaking, amplicons of this size are required only for the most degraded specimens, and the protocols described here have been specifically developed for recalcitrant human skeletal remains encountered during the course of a large-scale missing persons' identification effort. DNA templates from these types of specimens tend to exhibit various forms of intrastrand damage that, in turn, manifest as artifacts in the sequence data. Because these artifacts are not generally observed among sequence data from pristine templates, we address the particular data idiosyncrasies that warrant additional scrutiny. Although this chapter will primarily highlight this particular application, the basic experimental parameters and data considerations should easily extend to other applications and/or sample types. The protocols described here have been deliberately designed to produce raw sequence electropherograms and final mtDNA profiles that adhere to the strictest forensic guidelines in terms of overall data quality.

Ancient DNA from an Early Neolithic Iberian population supports a pioneer colonization by first farmers

The Neolithic transition has been widely debated particularly regarding the extent to which this revolution implied a demographic expansion from the Near East. We attempted to shed some light on this process in northeastern Iberia by combining ancient DNA (aDNA) data from Early Neolithic settlers and published DNA data from Middle Neolithic and modern samples from the same region. We successfully extracted and amplified mitochondrial DNA from 13 human specimens, found at three archaeological sites dated back to the Cardial culture in the Early Neolithic (Can Sadurní and Chaves) and to the Late Early Neolithic (Sant Pau del Camp). We found that haplogroups with a low frequency in modern populations-N* and X1-are found at higher frequencies in our Early Neolithic population (∼31%). Genetic differentiation between Early and Middle Neolithic populations was significant (F(ST) ∼0.13, P<10(-5)), suggesting that genetic drift played an important role at this time. To improve our understanding of the Neolithic demographic processes, we used a Bayesian coalescence-based simulation approach to identify the most likely of three demographic scenarios that might explain the genetic data. The three scenarios were chosen to reflect archaeological knowledge and previous genetic studies using similar inferential approaches. We found that models that ignore population structure, as previously used in aDNA studies, are unlikely to explain the data. Our results are compatible with a pioneer colonization of northeastern Iberia at the Early Neolithic characterized by the arrival of small genetically distinctive groups, showing cultural and genetic connections with the Near East.

Evaluating the impact of post-mortem damage in ancient DNA: a theoretical approach

The growth of ancient DNA research has offered exceptional opportunities and raised great expectations, but has also presented some considerable challenges. One of the ongoing issues is the impact of post-mortem damage in DNA molecules. Nucleotide alterations and DNA strand breakages lead to a significant decrease in the quantity of DNA molecules of useful length in a sample and to errors in the final DNA sequences obtained. We present a model of age-dependent DNA damage and quantify the influence of that damage on subsequent steps in the sequencing process, including the polymerase chain reaction and cloning. Calculations using our model show that deposition conditions, rather than the age of a sample, have the greatest influence on the level of DNA damage. In turn, this affects the probability of interpreting an erroneous (possessing damage-derived mutations) sequence as being authentic. We also evaluated the effect of post-mortem damage on real data sets using a Bayesian phylogenetic approach. According to our study, damage-derived sequence alterations appear to have little impact on the final DNA sequences. This indicates the effectiveness of current methods for sequence authentication and validation.

MACSE: Multiple Alignment of Coding SEquences accounting for frameshifts and stop codons

Until now the most efficient solution to align nucleotide sequences containing open reading frames was to use indirect procedures that align amino acid translation before reporting the inferred gap positions at the codon level. There are two important pitfalls with this approach. Firstly, any premature stop codon impedes using such a strategy. Secondly, each sequence is translated with the same reading frame from beginning to end, so that the presence of a single additional nucleotide leads to both aberrant translation and alignment.

We present an algorithm that has the same space and time complexity as the classical Needleman-Wunsch algorithm while accommodating sequencing errors and other biological deviations from the coding frame. The resulting pairwise coding sequence alignment method was extended to a multiple sequence alignment (MSA) algorithm implemented in a program called MACSE (Multiple Alignment of Coding SEquences accounting for frameshifts and stop codons). MACSE is the first automatic solution to align protein-coding gene datasets containing non-functional sequences (pseudogenes) without disrupting the underlying codon structure. It has also proved useful in detecting undocumented frameshifts in public database sequences and in aligning next-generation sequencing reads/contigs against a reference coding sequence.

MACSE is distributed as an open-source java file executable with freely available source code and can be used via a web interface at: http://mbb.univ-montp2.fr/macse.

True single-molecule DNA sequencing of a pleistocene horse bone

Second-generation sequencing platforms have revolutionized the field of ancient DNA, opening access to complete genomes of past individuals and extinct species. However, these platforms are dependent on library construction and amplification steps that may result in sequences that do not reflect the original DNA template composition. This is particularly true for ancient DNA, where templates have undergone extensive damage post-mortem. Here, we report the results of the first "true single molecule sequencing" of ancient DNA. We generated 115.9 Mb and 76.9 Mb of DNA sequences from a permafrost-preserved Pleistocene horse bone using the Helicos HeliScope and Illumina GAIIx platforms, respectively. We find that the percentage of endogenous DNA sequences derived from the horse is higher among the Helicos data than Illumina data. This result indicates that the molecular biology tools used to generate sequencing libraries of ancient DNA molecules, as required for second-generation sequencing, introduce biases into the data that reduce the efficiency of the sequencing process and limit our ability to fully explore the molecular complexity of ancient DNA extracts. We demonstrate that simple modifications to the standard Helicos DNA template preparation protocol further increase the proportion of horse DNA for this sample by threefold. Comparison of Helicos-specific biases and sequence errors in modern DNA with those in ancient DNA also reveals extensive cytosine deamination damage at the 3' ends of ancient templates, indicating the presence of 3'-sequence overhangs. Our results suggest that paleogenomes could be sequenced in an unprecedented manner by combining current second- and third-generation sequencing approaches.

Molecular identification of bacteria by total sequence screening: determining the cause of death in ancient human subjects

Research of ancient pathogens in ancient human skeletons has been mainly carried out on the basis of one essential historical or archaeological observation, permitting specific pathogens to be targeted. Detection of ancient human pathogens without such evidence is more difficult, since the quantity and quality of ancient DNA, as well as the environmental bacteria potentially present in the sample, limit the analyses possible. Using human lung tissue and/or teeth samples from burials in eastern Siberia, dating from the end of 17(th) to the 19(th) century, we propose a methodology that includes the: 1) amplification of all 16S rDNA gene sequences present in each sample; 2) identification of all bacterial DNA sequences with a degree of identity ≥ 95%, according to quality criteria; 3) identification and confirmation of bacterial pathogens by the amplification of the rpoB gene; and 4) establishment of authenticity criteria for ancient DNA. This study demonstrates that from teeth samples originating from ancient human subjects, we can realise: 1) the correct identification of bacterial molecular sequence signatures by quality criteria; 2) the separation of environmental and pathogenic bacterial 16S rDNA sequences; 3) the distribution of bacterial species for each subject and for each burial; and 4) the characterisation of bacteria specific to the permafrost. Moreover, we identified three pathogens in different teeth samples by 16S rDNA sequence amplification: Bordetella sp., Streptococcus pneumoniae and Shigella dysenteriae. We tested for the presence of these pathogens by amplifying the rpoB gene. For the first time, we confirmed sequences from Bordetella pertussis in the lungs of an ancient male Siberian subject, whose grave dated from the end of the 17(th) century to the early 18(th) century.

An algorithm for detecting high frequency copy number polymorphisms using SNP arrays

We present a general algorithm for the detection of genomic variants using the Illumina iSelect platform. The Illumina iSelect platform is designed to detect SNPs, but our algorithm allows for the detections of more general forms of variations, including copy number polymorphisms and microsatellites. The algorithm does not rely on a priori information of the type of polymorphism being studied and is designed to genotype call a large number of individuals simultaneously. The algorithm proceeds by initially normalizing intensity and correcting for batch effects. Then each marker is clustered using a modified Gaussian mixture model where we account for variances in the expression of an individuals and the variance measured in bead level intensities of a probe/marker pair. Finally, these clusters are used to determine genotypes. The algorithm was then run on a dataset of 35,000 Icelandic individuals.

Silencing of Her2, CCNB1 and PKC Genes by siRNA Results in Prolonged Retardation of Neuroblastoma Cell Division

Deregulation of the expression of the genes that are involved in the control of the cell cycle impairs cellular differentiation and leads to cell death. This process can result in uncontrollable cell proliferation and, subsequently, cancer development. In this study, we examined the effect of the silencing of cancer-related genes by small interfering RNAs (siRNA) targeted at mRNAof Her2, cyclin B1 (CCNB1), and protein kinase C(PKC) on the proliferation of human cancer cells of different origins. Maximum silencing ofCCNB1,Her2(in KB-3-1, SK-N-MC, MCF-7 cells), andPKC(in MCF-7 cells) was achieved 72 h after transfection of the corresponding siRNAs, and 12 days after the transfection, the initial levels of the target mRNAs were fully recovered. Silencing ofHer2,CCNB1,andPKCdifferently effected the proliferation of the cell lines under study. The most pronounced antiproliferative action of the investigated siRNAs was observed in neuroblastoma SK-N-MC cells (3 - 10-fold reduction in the proliferation rate) even after the recovery of the initial levels of expression ofthe Her2,CCNB1, andPKС genes. The obtained data indicate that theCCNB1 andPKCgenes can be used as targets in the development of drugs for neuroblastoma treatment.

To clone or not to clone: method analysis for retrieving consensus sequences in ancient DNA samples

The challenges associated with the retrieval and authentication of ancient DNA (aDNA) evidence are principally due to post-mortem damage which makes ancient samples particularly prone to contamination from "modern" DNA sources. The necessity for authentication of results has led many aDNA researchers to adopt methods considered to be "gold standards" in the field, including cloning aDNA amplicons as opposed to directly sequencing them. However, no standardized protocol has emerged regarding the necessary number of clones to sequence, how a consensus sequence is most appropriately derived, or how results should be reported in the literature. In addition, there has been no systematic demonstration of the degree to which direct sequences are affected by damage or whether direct sequencing would provide disparate results from a consensus of clones.To address this issue, a comparative study was designed to examine both cloned and direct sequences amplified from ∼3,500 year-old ancient northern fur seal DNA extracts. Majority rules and the Consensus Confidence Program were used to generate consensus sequences for each individual from the cloned sequences, which exhibited damage at 31 of 139 base pairs across all clones. In no instance did the consensus of clones differ from the direct sequence. This study demonstrates that, when appropriate, cloning need not be the default method, but instead, should be used as a measure of authentication on a case-by-case basis, especially when this practice adds time and cost to studies where it may be superfluous.