Sequencing the nuclear genome of the extinct woolly mammoth

Broader incorporation of bioinformatics in education: opportunities and challenges

The major opportunities for broader incorporation of bioinformatics in education can be placed into three general categories: general applicability of bioinformatics in life science and related curricula; inherent fit of bioinformatics for promoting student learning in most biology programs; and the general experience and associated comfort students have with computers and technology. Conversely, the major challenges for broader incorporation of bioinformatics in education can be placed into three general categories: required infrastructure and logistics; instructor knowledge of bioinformatics and continuing education; and the breadth of bioinformatics, and the diversity of students and educational objectives. Broader incorporation of bioinformatics at all education levels requires overcoming the challenges to using transformative computer-requiring learning activities, assisting faculty in collecting assessment data on mastery of student learning outcomes, as well as creating more faculty development opportunities that span diverse skill levels, with an emphasis placed on providing resource materials that are kept up-to-date as the field and tools change.

Neanderthal genomics and the evolution of modern humans

Humans possess unique physical and cognitive characteristics relative to other primates. Comparative analyses of the human and chimpanzee genomes are beginning to reveal sequence changes on the human lineage that may have contributed to the evolution of human traits. However, these studies cannot identify the genetic differences that distinguish modern humans from archaic human species. Here, I will discuss efforts to obtain genomic sequence from Neanderthal, the closest known relative of modern humans. Recent studies in this nascent field have focused on developing methods to recover nuclear DNA from Neanderthal remains. The success of these early studies has inspired a Neanderthal genome project, which promises to produce a reference Neanderthal genome sequence in the near future. Technical issues, such as the level of Neanderthal sequence coverage that can realistically be obtained from a single specimen and the presence of modern human contaminating sequences, reduce the detection of authentic human-Neanderthal sequence differences but may be remedied by methodological improvements. More critical for the utility of a Neanderthal genome sequence is the evolutionary relationship of humans and Neanderthals. Current evidence suggests that the modern human and Neanderthal lineages diverged before the emergence of contemporary humans. A fraction of biologically relevant human-chimpanzee sequence differences are thus likely to have arisen and become fixed exclusively on the modern human lineage. A reconstructed Neanderthal genome sequence could be integrated into human-primate genome comparisons to help reveal the evolutionary genetic events that produced modern humans.

Next Generation Sequencing of Ancient DNA: Requirements, Strategies and Perspectives

The invention of next-generation-sequencing has revolutionized almost all fields of genetics, but few have profited from it as much as the field of ancient DNA research. From its beginnings as an interesting but rather marginal discipline, ancient DNA research is now on its way into the centre of evolutionary biology. In less than a year from its invention next-generation-sequencing had increased the amount of DNA sequence data available from extinct organisms by several orders of magnitude. Ancient DNA research is now not only adding a temporal aspect to evolutionary studies and allowing for the observation of evolution in real time, it also provides important data to help understand the origins of our own species. Here we review progress that has been made in next-generation-sequencing of ancient DNA over the past five years and evaluate sequencing strategies and future directions.

Analysis of ancient human genomes: using next generation sequencing, 20-fold coverage of the genome of a 4,000-year-old human from Greenland has been obtained

High-capacity sequencing technologies have dramatically reduced both the cost and time required to generate complete human genome sequences. Besides expanding our knowledge about existing diversity, the nature of these technologies makes it possible to extend knowledge in yet another dimension: time. Recently, the complete genome sequence of a 4,000-year-old human from the Saqqaq culture of Greenland was determined to 20-fold coverage. These data make it possible to investigate the population affinities of this enigmatic culture and, by identifying several phenotypic traits of this individual, provide a limited glimpse into how these people may have looked. While undoubtedly a milestone in ancient DNA research, the cost to generate an ancient genome, even from such an exceptionally preserved specimen, remains out of reach for most. Nonetheless, recently developed DNA capture methods, already applied to Neanderthal and fossil human mitochondrial DNA, may soon make large-scale genome-wide analysis of ancient human diversity a reality, providing a fresh look at human population history.

Road blocks on paleogenomes--polymerase extension profiling reveals the frequency of blocking lesions in ancient DNA

Although the last few years have seen great progress in DNA sequence retrieval from fossil specimens, some of the characteristics of ancient DNA remain poorly understood. This is particularly true for blocking lesions, i.e. chemical alterations that cannot be bypassed by DNA polymerases and thus prevent amplification and subsequent sequencing of affected molecules. Some studies have concluded that the vast majority of ancient DNA molecules carry blocking lesions, suggesting that the removal, repair or bypass of blocking lesions might dramatically increase both the time depth and geographical range of specimens available for ancient DNA analysis. However, previous studies used very indirect detection methods that did not provide conclusive estimates on the frequency of blocking lesions in endogenous ancient DNA. We developed a new method, polymerase extension profiling (PEP), that directly reveals occurrences of polymerase stalling on DNA templates. By sequencing thousands of single primer extension products using PEP methodology, we have for the first time directly identified blocking lesions in ancient DNA on a single molecule level. Although we found clear evidence for blocking lesions in three out of four ancient samples, no more than 40% of the molecules were affected in any of the samples, indicating that such modifications are far less frequent in ancient DNA than previously thought.

Targeted investigation of the Neandertal genome by array-based sequence capture

It is now possible to perform whole-genome shotgun sequencing as well as capture of specific genomic regions for extinct organisms. However, targeted resequencing of large parts of nuclear genomes has yet to be demonstrated for ancient DNA. Here we show that hybridization capture on microarrays can successfully recover more than a megabase of target regions from Neandertal DNA even in the presence of approximately 99.8% microbial DNA. Using this approach, we have sequenced approximately 14,000 protein-coding positions inferred to have changed on the human lineage since the last common ancestor shared with chimpanzees. By generating the sequence of one Neandertal and 50 present-day humans at these positions, we have identified 88 amino acid substitutions that have become fixed in humans since our divergence from the Neandertals.

A draft sequence of the Neandertal genome

Neandertals, the closest evolutionary relatives of present-day humans, lived in large parts of Europe and western Asia before disappearing 30,000 years ago. We present a draft sequence of the Neandertal genome composed of more than 4 billion nucleotides from three individuals. Comparisons of the Neandertal genome to the genomes of five present-day humans from different parts of the world identify a number of genomic regions that may have been affected by positive selection in ancestral modern humans, including genes involved in metabolism and in cognitive and skeletal development. We show that Neandertals shared more genetic variants with present-day humans in Eurasia than with present-day humans in sub-Saharan Africa, suggesting that gene flow from Neandertals into the ancestors of non-Africans occurred before the divergence of Eurasian groups from each other.

Computational challenges in the analysis of ancient DNA

A new method of next-generation sequencing analysis is presented which takes into account the biases characteristic of ancient, including Neandertal, DNA samples.

High-throughput sequencing technologies have opened up a new avenue for studying extinct organisms. Here we identify and quantify biases introduced by particular characteristics of ancient DNA samples. These analyses demonstrate the importance of closely related genomic sequence for correctly identifying and classifying bona fide endogenous DNA fragments. We show that more accurate genome divergence estimates from ancient DNA sequence can be attained using at least two outgroup genomes and appropriate filtering.

Discovery and evaluation of candidate sex-determining genes and xenobiotics in the gonads of lake sturgeon (Acipenser fulvescens)

Modern pyrosequencing has the potential to uncover many interesting aspects of genome evolution, even in lineages where genomic resources are scarce. In particular, 454 pyrosequencing of nonmodel species has been used to characterize expressed sequence tags, xenobiotics, gene ontologies, and relative levels of gene expression. Herein, we use pyrosequencing to study the evolution of genes expressed in the gonads of a polyploid fish, the lake sturgeon (Acipenser fulvescens). Using 454 pyrosequencing of transcribed genes, we produced more than 125 MB of sequence data from 473,577 high-quality sequencing reads. Sequences that passed stringent quality control thresholds were assembled into 12,791 male contigs and 32,629 female contigs. Average depth of coverage was 4.2 x for the male assembly and 5.5x for the female assembly. Analytical rarefaction indicates that our assemblies include most of the genes expressed in lake sturgeon gonads. Over 86,700 sequencing reads were assigned gene ontologies, many to general housekeeping genes like protein, RNA, and ion binding genes. We searched specifically for sex determining genes and documented significant sex differences in the expression of two genes involved in animal sex determination, DMRT1 and TRA-1. DMRT1 is the master sex determining gene in birds and in medaka (Oryzias latipes) whereas TRA-1 helps direct sexual differentiation in nematodes. We also searched the lake sturgeon assembly for evidence of xenobiotic organisms that may exist as endosymbionts. Our results suggest that exogenous parasites (trematodes) and pathogens (protozoans) apparently have infected lake sturgeon gonads, and the trematodes have horizontally transferred some genes to the lake sturgeon genome.

Ancient human genome sequence of an extinct Palaeo-Eskimo

We report here the genome sequence of an ancient human. Obtained from approximately 4,000-year-old permafrost-preserved hair, the genome represents a male individual from the first known culture to settle in Greenland. Sequenced to an average depth of 20x, we recover 79% of the diploid genome, an amount close to the practical limit of current sequencing technologies. We identify 353,151 high-confidence single-nucleotide polymorphisms (SNPs), of which 6.8% have not been reported previously. We estimate raw read contamination to be no higher than 0.8%. We use functional SNP assessment to assign possible phenotypic characteristics of the individual that belonged to a culture whose location has yielded only trace human remains. We compare the high-confidence SNPs to those of contemporary populations to find the populations most closely related to the individual. This provides evidence for a migration from Siberia into the New World some 5,500 years ago, independent of that giving rise to the modern Native Americans and Inuit.

Fossil avian eggshell preserves ancient DNA

Owing to exceptional biomolecule preservation, fossil avian eggshell has been used extensively in geochronology and palaeodietary studies. Here, we show, to our knowledge, for the first time that fossil eggshell is a previously unrecognized source of ancient DNA (aDNA). We describe the successful isolation and amplification of DNA from fossil eggshell up to 19 ka old. aDNA was successfully characterized from eggshell obtained from New Zealand (extinct moa and ducks), Madagascar (extinct elephant birds) and Australia (emu and owl). Our data demonstrate excellent preservation of the nucleic acids, evidenced by retrieval of both mitochondrial and nuclear DNA from many of the samples. Using confocal microscopy and quantitative PCR, this study critically evaluates approaches to maximize DNA recovery from powdered eggshell. Our quantitative PCR experiments also demonstrate that moa eggshell has approximately 125 times lower bacterial load than bone, making it a highly suitable substrate for high-throughput sequencing approaches. Importantly, the preservation of DNA in Pleistocene eggshell from Australia and Holocene deposits from Madagascar indicates that eggshell is an excellent substrate for the long-term preservation of DNA in warmer climates. The successful recovery of DNA from this substrate has implications in a number of scientific disciplines; most notably archaeology and palaeontology, where genotypes and/or DNA-based species identifications can add significantly to our understanding of diets, environments, past biodiversity and evolutionary processes.

A complete mitochondrial genome sequence from a mesolithic wild aurochs (Bos primigenius)

BACKGROUND

The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample.

METHODOLOGY

DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738+/-68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences.

CONCLUSIONS

For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.

A rapid column-based ancient DNA extraction method for increased sample throughput

Genetic analyses using museum specimens and ancient DNA from fossil samples are becoming increasingly important in phylogenetic and especially population genetic studies. Recent progress in ancient DNA sequencing technologies has substantially increased DNA sequence yields and, in combination with barcoding methods, has enabled large-scale studies using any type of DNA. Moreover, more and more studies now use nuclear DNA sequences in addition to mitochondrial ones. Unfortunately, nuclear DNA is, due to its much lower copy number in living cells compared to mitochondrial DNA, much more difficult to obtain from low-quality samples. Therefore, a DNA extraction method that optimizes DNA yields from low-quality samples and at the same time allows processing many samples within a short time frame is immediately required. In fact, the major bottleneck in the analysis process using samples containing low amounts of degraded DNA now lies in the extraction of samples, as column-based methods using commercial kits are fast but have proven to give very low yields, while more efficient methods are generally very time-consuming. Here, we present a method that combines the high DNA yield of batch-based silica extraction with the time-efficiency of column-based methods. Our results on Pleistocene cave bear samples show that DNA yields are quantitatively comparable, and in fact even slightly better than with silica batch extraction, while at the same time the number of samples that can conveniently be processed in parallel increases and both bench time and costs decrease using this method. Thus, this method is suited for harvesting the power of high-throughput sequencing using the DNA preserved in the millions of paleontological and museums specimens.

Molecular evolution of dentin phosphoprotein among toothed and toothless animals

Background

Dentin sialophosphoprotein (DSPP) is the largest member of the SIBLING family and is the most abundant noncollagenous protein in dentin. DSPP is also expressed in non-mineralized tissues including metabolically active ductal epithelia and some cancers. Its function, however, is poorly defined. The carboxy-terminal fragment, dentin phosphoprotein (DPP) is encoded predominantly by a large repetitive domain that requires separate cloning/sequencing reactions and is, therefore, often incomplete in genomic databases. Comparison of DPP sequences from at least one member of each major branch in the mammalian evolutionary tree (including some "toothless" mammals) as well as one reptile and bird may help delineate its possible functions in both dentin and ductal epithelia.

Results

The BMP1-cleavage and translation-termination domains were sufficiently conserved to permit amplification/cloning/sequencing of most species' DPP. While the integrin-binding domain, RGD, was present in about half of species, only vestigial remnants of this tripeptide were identified in the others. The number of tandem repeats of the nominal SerSerAsp phosphorylation motif in toothed mammals (including baleen whale and platypus which lack teeth as adults), ranged from ~75 (elephant) to >230 (human). These repeats were not perfect, however, and patterns of intervening sequences highlight the rapidity of changes among even closely related species. Two toothless anteater species have evolved different sets of nonsense mutations shortly after their BMP1 motifs suggesting that while cleavage may be important for DSPP processing in other tissues, the DPP domain itself may be required only in dentin. The lizard DSPP had an intact BMP1 site, a remnant RGD motif, as well as a distinctly different Ser/Asp-rich domain compared to mammals.

Conclusions

The DPP domain of DSPP was found to change dramatically within mammals and was lost in two truly toothless animals. The defining aspect of DPP, the long repeating phosphorylation domain, apparently undergoes frequent slip replication and recombination events that rapidly change specific patterns but not its overall biochemical character in toothed animals. Species may have to co-evolve protein processing mechanisms, however, to handle increased lengths of DSP repeats. While the RGD domain is lost in many species, some evolutionary pressure to maintain integrin binding can be observed.

Species identification by analysis of bone collagen using matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry

Species identification of fragmentary bone, such as in rendered meat and bone meal or from archaeological sites, is often difficult in the absence of clear morphological markers. Here we present a robust method of analysing genus-specific collagen peptides by mass spectrometry simply by using solid-phase extraction (a C18 ZipTip) for peptide purification, rather than liquid chromatography/mass spectrometry (LC/MS). Analysis of the collagen from 32 different mammal species identified a total of 92 peptide markers that could be used for species identification, for example, in processed food and animal feed. A set of ancient (>100 ka@10 degrees C) bone samples was also analysed to show that the proposed method has applications to archaeological bone identification.

Shotgun metagenomics of biological stains using ultra-deep DNA sequencing

A detailed molecular analysis of blood or other biological stains at a crime scene is often hampered by the low quantity and quality of the extractable DNA. However, the determination of the origin and composition of a stain is in most cases a prerequisite for the final elucidation of a criminal case. Standard methodologies, e.g. amplification of DNA followed by microsatellite typing or mitochondrial DNA sequencing, are often not sensitive enough to result in sufficient and conclusive data. We have applied ultra-deep DNA sequencing using the 454 pyrosequencing technology on a whole genome amplified (WGA) environmental biological stain, which was analysed unsuccessfully with standard methodologies following WGA. With the combination of WGA and 454 pyrosequencing, however, we were able to generate 7242 single sequences with an average length of 195bp. A total of 1,441,971bp DNA sequences were generated and compared with public DNA sequence databases. Using RepeatMasker and basic logical alignment search tool (BLAST) searches against known microbial and mammalian genomes it was possible to determine the metagenomic composition of the stain, i.e. 4.2% bacterial DNA, 0.3% viral DNA, 2.7% fungal DNA, 10.3% mammalian repetitive DNA, 0.9% porcine DNA, 0.13% human DNA and 81.5% DNA of unknown origin. Our data demonstrate that 454 pyrosequencing has the potential to become a powerful tool not only in basic research but also in the metagenomic analysis of biological trace materials for forensic genetics.

Reading between the LINEs to see into the past

Transposable elements (TEs) are an important source of genome diversity and play a crucial role in genome evolution. A recent study by Zhao et al. describes novel patterns of TE diversification in the genome of the extinct mammoth Mammuthus primigenius. Analysis of Mammuthus has provided a unique genome landscape, a pivotal species for understanding TEs and genome evolution and hints at the diversity we verge on discovering by expanding our taxonomic sampling among genomes. Strategies based on this work might also revolutionize investigations of the interface between TE dynamics and genome diversity.

Processed pseudogenes: the 'fossilized footprints' of past gene expression

Although our knowledge of the genes and genomes of extinct organisms is improving as a result of progress in sequencing ancient DNA, the transcriptomes of extinct organisms remain inaccessible, owing to the rapid degradation of messenger RNA after death. We provide empirical evidence that gene expression levels in the reproductive tissues of mice and during early mouse development correlate highly with the rate of inherited retroposition: the source of processed pseudogenes in the genome. Thus, processed pseudogenes might serve as fossilized footprints of the expression of their parent genes, shedding light on ancient transcriptomes that could provide significant insights into the evolution of gene expression.

Ancient mitogenomics

The mitochondrial genome has been the traditional focus of most research into ancient DNA, owing to its high copy number and population-level variability. Despite this long-standing interest in mitochondrial DNA, it was only in 2001 that the first complete ancient mitogenomic sequences were obtained. As a result of various methodological developments, including the introduction of high-throughput sequencing techniques, the total number of ancient mitogenome sequences has increased rapidly over the past few years. In this review, we present a brief history of ancient mitogenomics and describe the technical challenges that face researchers in the field. We catalogue the diverse sequencing methods and source materials used to obtain ancient mitogenomic sequences, summarise the associated genetic and phylogenetic studies that have been conducted, and evaluate the future prospects of the field.

The 'stem cell' concept: is it holding us back?

Developmental biology, regenerative medicine and cancer biology are increasingly occupied with the molecular characterization of stem cells. Yet recent work adds to a growing body of literature suggesting that 'stemness' cannot be reduced to the molecular features of cell types, and is instead an emergent property of cell lineages under feedback control.

COMPASS: a program for generating serial samples under an infinite sites model

SUMMARY

The program COMPASS can generate samples that have been collected at various points in time from a population that is evolving according to a Wright-Fisher model. The samples are generated using coalescence simulations permitting various demographic scenarios and the program uses an infinite sites model to generate polymorphism data for the samples. By generating serially sampled population-genetic data, COMPASS allows investigating properties of polymorphism data that has been collected at different time points, and aid in making inference from ancient polymorphism data.

AVAILABILITY

The program and the manual are available at: http://www.egs.uu.se/evbiol/Research/JakobssonLab/compass.html

CONTACT

mattias.jakobsson@ebc.uu.se.

Genetic variation at hair length candidate genes in elephants and the extinct woolly mammoth

Background

Like humans, the living elephants are unusual among mammals in being sparsely covered with hair. Relative to extant elephants, the extinct woolly mammoth, Mammuthus primigenius, had a dense hair cover and extremely long hair, which likely were adaptations to its subarctic habitat. The fibroblast growth factor 5 (FGF5) gene affects hair length in a diverse set of mammalian species. Mutations in FGF5 lead to recessive long hair phenotypes in mice, dogs, and cats; and the gene has been implicated in hair length variation in rabbits. Thus, FGF5 represents a leading candidate gene for the phenotypic differences in hair length notable between extant elephants and the woolly mammoth. We therefore sequenced the three exons (except for the 3' UTR) and a portion of the promoter of FGF5 from the living elephantid species (Asian, African savanna and African forest elephants) and, using protocols for ancient DNA, from a woolly mammoth.

Results

Between the extant elephants and the mammoth, two single base substitutions were observed in FGF5, neither of which alters the amino acid sequence. Modeling of the protein structure suggests that the elephantid proteins fold similarly to the human FGF5 protein. Bioinformatics analyses and DNA sequencing of another locus that has been implicated in hair cover in humans, type I hair keratin pseudogene (KRTHAP1), also yielded negative results. Interestingly, KRTHAP1 is a pseudogene in elephantids as in humans (although fully functional in non-human primates).

Conclusion

The data suggest that the coding sequence of the FGF5 gene is not the critical determinant of hair length differences among elephantids. The results are discussed in the context of hairlessness among mammals and in terms of the potential impact of large body size, subarctic conditions, and an aquatic ancestor on hair cover in the Proboscidea.

The Neandertal genome and ancient DNA authenticity

Recent advances in high-thoughput DNA sequencing have made genome-scale analyses of genomes of extinct organisms possible. With these new opportunities come new difficulties in assessing the authenticity of the DNA sequences retrieved. We discuss how these difficulties can be addressed, particularly with regard to analyses of the Neandertal genome. We argue that only direct assays of DNA sequence positions in which Neandertals differ from all contemporary humans can serve as a reliable means to estimate human contamination. Indirect measures, such as the extent of DNA fragmentation, nucleotide misincorporations, or comparison of derived allele frequencies in different fragment size classes, are unreliable. Fortunately, interim approaches based on mtDNA differences between Neandertals and current humans, detection of male contamination through Y chromosomal sequences, and repeated sequencing from the same fossil to detect autosomal contamination allow initial large-scale sequencing of Neandertal genomes. This will result in the discovery of fixed differences in the nuclear genome between Neandertals and current humans that can serve as future direct assays for contamination. For analyses of other fossil hominins, which may become possible in the future, we suggest a similar ‘boot-strap' approach in which interim approaches are applied until sufficient data for more definitive direct assays are acquired.

Direct multiplex sequencing (DMPS)--a novel method for targeted high-throughput sequencing of ancient and highly degraded DNA

Although the emergence of high-throughput sequencing technologies has enabled whole-genome sequencing from extinct organisms, little progress has been made in accelerating targeted sequencing from highly degraded DNA. Here, we present a novel and highly sensitive method for targeted sequencing of ancient and degraded DNA, which couples multiplex PCR directly with sample barcoding and high-throughput sequencing. Using this approach, we obtained a 96% complete mitochondrial genome data set from 31 cave bear (Ursus spelaeus) samples using only two 454 Life Sciences (Roche) GS FLX runs. In contrast to previous studies relying only on short sequence fragments, the overlapping portion of our data comprises almost 10 kb of replicated mitochondrial genome sequence, allowing for the unambiguous differentiation of three major cave bear clades. Our method opens up the opportunity to simultaneously generate many kilobases of overlapping sequence data from large sets of difficult samples, such as museum specimens, medical collections, or forensic samples. Embedded in our approach, we present a new protocol for the construction of barcoded sequencing libraries, which is compatible with all current high-throughput technologies and can be performed entirely in plate setup.

Tracking the past: interspersed repeats in an extinct Afrotherian mammal, Mammuthus primigenius

The woolly mammoth (Mammuthus primigenius) died out about several thousand years ago, yet recent paleogenomic studies have successfully recovered genetic information from both the mitochondrial and nuclear genomes of this extinct species. Mammoths belong to Afrotheria, a group of mammals exhibiting extreme morphological diversity and large genome sizes. In this study, we found that the mammoth genome contains a larger proportion of interspersed repeats than any other mammalian genome reported so far, in which the proliferation of the RTE family of retrotransposons (covering 12% of the genome) may be the main reason for an increased genome size. Phylogenetic analysis showed that RTEs in mammoth are closely related to the family BovB/RTE. The incongruence of the reconstructed RTE phylogeny indicates that RTEs in mammoth may be acquired through an ancient lateral gene transfer event. A recent proliferation of SINEs was also found in the probocidean lineage, whereas the Afrotherian-wide SINEs in mammoth have undergone a rather flat and stepwise expansion. Comparisons of the transposable elements (TEs) between mammoth and other mammals may shed light on the evolutionary history of TEs in various mammalian lineages.

Paleogenomics in a temperate environment: shotgun sequencing from an extinct Mediterranean caprine

Background

Numerous endemic mammals, including dwarf elephants, goats, hippos and deers, evolved in isolation in the Mediterranean islands during the Pliocene and Pleistocene. Most of them subsequently became extinct during the Holocene. Recently developed high-throughput sequencing technologies could provide a unique tool for retrieving genomic data from these extinct species, making it possible to study their evolutionary history and the genetic bases underlying their particular, sometimes unique, adaptations.

Methodology/Principals Findings

A DNA extraction of a ∼6,000 year-old bone sample from an extinct caprine (Myotragus balearicus) from the Balearic Islands in the Western Mediterranean, has been subjected to shotgun sequencing with the GS FLX 454 platform. Only 0.27% of the resulting sequences, identified from alignments with the cow genome and comprising 15,832 nucleotides, with an average length of 60 nucleotides, proved to be endogenous.

Conclusions

A phylogenetic tree generated with Myotragus sequences and those from other artiodactyls displays an identical topology to that generated from mitochondrial DNA data. Despite being in an unfavourable thermal environment, which explains the low yield of endogenous sequences, our study demonstrates that it is possible to obtain genomic data from extinct species from temperate regions.

Next-generation sequencing of vertebrate experimental organisms

Next-generation sequencing technologies are revolutionizing biology by allowing for genome-wide transcription factor binding-site profiling, transcriptome sequencing, and more recently, whole-genome resequencing. While it is currently not possible to generate complete de novo assemblies of higher-vertebrate genomes using next-generation sequencing, improvements in sequence read lengths and throughput, coupled with new assembly algorithms for large data sets, will soon make this a reality. These developments will in turn spawn a revolution in how genomic data are used to understand genetics and how model organisms are used for disease gene discovery. This review provides an overview of the current next-generation sequencing platforms and the newest computational tools for the analysis of next-generation sequencing data. We also describe how next-generation sequencing may be applied in the context of vertebrate model organism genetics.

Analysis of complete mitochondrial genomes from extinct and extant rhinoceroses reveals lack of phylogenetic resolution

Background

The scientific literature contains many examples where DNA sequence analyses have been used to provide definitive answers to phylogenetic problems that traditional (non-DNA based) approaches alone have failed to resolve. One notable example concerns the rhinoceroses, a group for which several contradictory phylogenies were proposed on the basis of morphology, then apparently resolved using mitochondrial DNA fragments.

Results

In this study we report the first complete mitochondrial genome sequences of the extinct ice-age woolly rhinoceros (Coelodonta antiquitatis), and the threatened Javan (Rhinoceros sondaicus), Sumatran (Dicerorhinus sumatrensis), and black (Diceros bicornis) rhinoceroses. In combination with the previously published mitochondrial genomes of the white (Ceratotherium simum) and Indian (Rhinoceros unicornis) rhinoceroses, this data set putatively enables reconstruction of the rhinoceros phylogeny. While the six species cluster into three strongly supported sister-pairings: (i) The black/white, (ii) the woolly/Sumatran, and (iii) the Javan/Indian, resolution of the higher-level relationships has no statistical support. The phylogenetic signal from individual genes is highly diffuse, with mixed topological support from different genes. Furthermore, the choice of outgroup (horse vs tapir) has considerable effect on reconstruction of the phylogeny. The lack of resolution is suggestive of a hard polytomy at the base of crown-group Rhinocerotidae, and this is supported by an investigation of the relative branch lengths.

Conclusion

Satisfactory resolution of the rhinoceros phylogeny may not be achievable without additional analyses of substantial amounts of nuclear DNA. This study provides a compelling demonstration that, in spite of substantial sequence length, there are significant limitations with single-locus phylogenetics. We expect further examples of this to appear as next-generation, large-scale sequencing of complete mitochondrial genomes becomes commonplace in evolutionary studies.

"The human factor in classification is nowhere more evident than in dealing with this superfamily (Rhinocerotoidea)." G. G. Simpson (1945)

The mitochondrial genome sequence of the Tasmanian tiger (Thylacinus cynocephalus)

We report the first two complete mitochondrial genome sequences of the thylacine (Thylacinus cynocephalus), or so-called Tasmanian tiger, extinct since 1936. The thylacine's phylogenetic position within australidelphian marsupials has long been debated, and here we provide strong support for the thylacine's basal position in Dasyuromorphia, aided by mitochondrial genome sequence that we generated from the extant numbat (Myrmecobius fasciatus). Surprisingly, both of our thylacine sequences differ by 11%-15% from putative thylacine mitochondrial genes in GenBank, with one of our samples originating from a direct offspring of the previously sequenced individual. Our data sample each mitochondrial nucleotide an average of 50 times, thereby providing the first high-fidelity reference sequence for thylacine population genetics. Our two sequences differ in only five nucleotides out of 15,452, hinting at a very low genetic diversity shortly before extinction. Despite the samples' heavy contamination with bacterial and human DNA and their temperate storage history, we estimate that as much as one-third of the total DNA in each sample is from the thylacine. The microbial content of the two thylacine samples was subjected to metagenomic analysis, and showed striking differences between a wild-captured individual and a born-in-captivity one. This study therefore adds to the growing evidence that extensive sequencing of museum collections is both feasible and desirable, and can yield complete genomes.

Resurrection of a bull by cloning from organs frozen without cryoprotectant in a -80 degrees c freezer for a decade

Frozen animal tissues without cryoprotectant have been thought to be inappropriate for use as a nuclear donor for somatic cell nuclear transfer (SCNT). We report the cloning of a bull using cells retrieved from testicles that had been taken from a dead animal and frozen without cryoprotectant in a −80°C freezer for 10 years. We obtained live cells from defrosted pieces of the spermatic cords of frozen testicles. The cells proliferated actively in culture and were apparently normal. We transferred 16 SCNT embryos from these cells into 16 synchronized recipient animals. We obtained five pregnancies and four cloned calves developed to term. Our results indicate that complete genome sets are maintained in mammalian organs even after long-term frozen-storage without cryoprotectant, and that live clones can be produced from the recovered cells.