Rates of evolution in ancient DNA from Adélie penguins

Analysis of the sequence variations in the Mhc DRB1-like gene of the endangered Humboldt penguin (Spheniscus humboldti)

The Major Histocompatibility Complex (Mhc) genomic region of many vertebrates is known to contain at least one highly polymorphic class II gene that is homologous in sequence to one or other of the human Mhc DRB1 class II genes. The diversity of the avian Mhc class II gene sequences have been extensively studied in chickens, quails, and some songbirds, but have been largely ignored in the oceanic birds, including the flightless penguins. We have previously reported that several penguin species have a high degree of polymorphism on exon 2 of the Mhc class II DRB1-like gene. In this study, we present for the first time the complete nucleotide sequences of exon 2, intron 2, and exon 3 of the DRB1-like gene of 20 Humboldt penguins, a species that is presently vulnerable to the dangers of extinction. The Humboldt DRB1-like nucleotide and amino acid sequences reveal at least eight unique alleles. Phylogenetic analysis of all the available avian DRB-like sequences showed that, of five penguin species and nine other bird species, the sequences of the Humboldt penguins grouped most closely to the Little penguin and the mallard, respectively. The present analysis confirms that the sequence variations of the Mhc class II gene, DRB1, are useful for discriminating among individuals within the same penguin population as well those within different penguin population groups and species.

Ancient DNA as a multidisciplinary experience

Investigation into DNA from archeological remains offers an inestimable tool for unraveling the history of humankind. However, a series of basic and technical difficulties renders the analysis of ancient DNA (aDNA) molecules troublesome, depending either on their own peculiar characteristics or on the complexity of processes affecting the bone matrix over time, all compromising the preservation of ancient DNA. This review underlines the contribution of many different disciplines, in particular molecular biology and genetics, to overcome these obstacles. The role of each expertise is illustrated to appropriately address the questions arising in aDNA investigations.

Studying the effect of environmental change on biotic evolution: past genetic contributions, current work and future directions

Evolutionary geneticists currently face a major scientific opportunity when integrating across the rapidly increasing amount of genetic data and existing biological scenarios based on ecology, fossils or climate models. Although genetic data acquisition and analysis have improved tremendously, several limitations remain. Here, we discuss the feedback between history and genetic variation in the face of environmental change with increasing taxonomic and temporal scale, as well as the major challenges that lie ahead. In particular, we focus on recent developments in two promising genetic methods, those of 'phylochronology' and 'molecular clocks'. With the advent of ancient DNA techniques, we can now directly sample the recent past. We illustrate this amazing and largely untapped utility of ancient DNA extracted from accurately dated localities with documented environmental changes. Innovative statistical analyses of these genetic data expose the direct effect of recent environmental change on genetic endurance, or maintenance of genetic variation. The 'molecular clock' (assumption of a linear relationship between genetic distance and evolutionary time) has been used extensively in phylogenetic studies to infer time and correlation between lineage divergence time and concurrent environmental change. Several studies at both population and species scale support a persuasive relationship between particular perturbation events and time of biotic divergence. However, we are still a way from gleaning an overall pattern to this relationship, which is a prerequisite to ultimately understanding the mechanisms by which past environments have shaped the evolutionary trajectory. Current obstacles include as-yet undecided reasons behind the frequent discrepancy between molecular and fossil time estimates, and the frequent lack of consideration of extensive confidence intervals around time estimates. We suggest that use and interpretation of both ancient DNA and molecular clocks is most effective when results are synthesized with palaeontological (fossil) and ecological (life history) information.

Genetic Analyses from Ancient DNA

About 20 years ago, DNA sequences were separately described from the quagga (a type of zebra) and an ancient Egyptian individual. What made these DNA sequences exceptional was that they were derived from 140- and 2400-year-old specimens. However, ancient DNA research, defined broadly as the retrieval of DNA sequences from museum specimens, archaeological finds, fossil remains, and other unusual sources of DNA, only really became feasible with the advent of techniques for the enzymatic amplification of specific DNA sequences. Today, reports of analyses of specimens hundreds, thousands, and even millions of years old are almost commonplace. But can all these results be believed? In this paper, we critically assess the state of ancient DNA research. In particular, we discuss the precautions and criteria necessary to ascertain to the greatest extent possible that results represent authentic ancient DNA sequences. We also highlight some significant results and areas of promising future research.

Rise and Fall of the Beringian Steppe Bison

The widespread extinctions of large mammals at the end of the Pleistocene epoch have often been attributed to the depredations of humans; here we present genetic evidence that questions this assumption. We used ancient DNA and Bayesian techniques to reconstruct a detailed genetic history of bison throughout the late Pleistocene and Holocene epochs. Our analyses depict a large diverse population living throughout Beringia until around 37,000 years before the present, when the population's genetic diversity began to decline dramatically. The timing of this decline correlates with environmental changes associated with the onset of the last glacial cycle, whereas archaeological evidence does not support the presence of large populations of humans in Eastern Beringia until more than 15,000 years later.

Serial SimCoal: A population genetics model for data from multiple populations and points in time

We present Serial SimCoal, a program that models population genetic data from multiple time points, as with ancient DNA data. An extension of SIMCOAL, it also allows simultaneous modeling of complex demographic histories, and migration between multiple populations. Further, we incorporate a statistical package to calculate relevant summary statistics, which, for the first time allows users to investigate the statistical power provided by, conduct hypothesis-testing with, and explore sample size limitations of ancient DNA data.

AVAILABILITY

Source code and Windows/Mac executables at http://www.stanford.edu/group/hadlylab/ssc.html

CONTACT

senka@stanford.edu.

2000 Year-old ancient equids: an ancient-DNA lesson from pompeii remains

Ancient DNA extracted from 2000 year-old equine bones was examined in order to amplify mitochondrial and nuclear DNA fragments. A specific equine satellite-type sequence representing 3.7%-11% of the entire equine genome, proved to be a suitable target to address the question of the presence of aDNA in ancient bones. The PCR strategy designed to investigate this specific target also allowed us to calculate the molecular weight of amplifiable DNA fragments. Sequencing of a 370 bp DNA fragment of mitochondrial control region allowed the comparison of ancient DNA sequences with those of modern horses to assess their genetic relationship. The 16S rRNA mitochondrial gene was also examined to unravel the post-mortem base modification feature and to test the status of Pompeian equids taxon on the basis of a Mae III restriction site polymorphism.

The incomplete natural history of mitochondria

Mitochondrial DNA (mtDNA) has been used to study molecular ecology and phylogeography for 25 years. Much important information has been gained in this way, but it is time to reflect on the biology of the mitochondrion itself and consider opportunities for evolutionary studies of the organelle itself and its ecology, biochemistry and physiology. This review has four sections. First, we review aspects of the natural history of mitochondria and their DNA to show that it is a unique molecule with specific characteristics that differ from nuclear DNA. We do not attempt to cover the plethora of differences between mitochondrial and nuclear DNA; rather we spotlight differences that can cause significant bias when inferring demographic properties of populations and/or the evolutionary history of species. We focus on recombination, effective population size and mutation rate. Second, we explore some of the difficulties in interpreting phylogeographical data from mtDNA data alone and suggest a broader use of multiple nuclear markers. We argue that mtDNA is not a sufficient marker for phylogeographical studies if the focus of the investigation is the species and not the organelle. We focus on the potential bias caused by introgression. Third, we show that it is not safe to assume a priori that mtDNA evolves as a strictly neutral marker because both direct and indirect selection influence mitochondria. We outline some of the statistical tests of neutrality that can, and should, be applied to mtDNA sequence data prior to making any global statements concerning the history of the organism. We conclude with a critical examination of the neglected biology of mitochondria and point out several surprising gaps in the state of our knowledge about this important organelle. Here we limelight mitochondrial ecology, sexually antagonistic selection, life-history evolution including ageing and disease, and the evolution of mitochondrial inheritance.

Genetic response to climatic change: insights from ancient DNA and phylochronology

Understanding how climatic change impacts biological diversity is critical to conservation. Yet despite demonstrated effects of climatic perturbation on geographic ranges and population persistence, surprisingly little is known of the genetic response of species. Even less is known over ecologically long time scales pertinent to understanding the interplay between microevolution and environmental change. Here, we present a study of population variation by directly tracking genetic change and population size in two geographically widespread mammal species (Microtus montanus and Thomomys talpoides) during late-Holocene climatic change. We use ancient DNA to compare two independent estimates of population size (ecological and genetic) and corroborate our results with gene diversity and serial coalescent simulations. Our data and analyses indicate that, with population size decreasing at times of climatic change, some species will exhibit declining gene diversity as expected from simple population genetic models, whereas others will not. While our results could be consistent with selection, independent lines of evidence implicate differences in gene flow, which depends on the life history strategy of species.

Phylogeography of a high-dispersal New Zealand sea-star: does upwelling block gene-flow?

New Zealand's (NZ) geographical isolation, extensive coastline and well-characterized oceanography offer a valuable system for marine biogeographical research. Here we use mtDNA control region sequences in the abundant endemic sea-star Patiriella regularis to test the following literature-based predictions: that coastal upwelling disrupts north-south gene flow and promotes population differentiation (hypothesis 1); and that an invasive Tasmanian population of the species was introduced anthropogenically from southern New Zealand (hypothesis 2). We sequenced 114 samples from 22 geographical locations, including nine sites from North Island, nine from South Island, one from Stewart Island and three from Tasmania. Our analysis of these sequences revealed an abundance of shallow phylogenetic lineages within P. regularis (68 haplotypes, mean divergence 0.9%). We detected significant genetic heterogeneity between pooled samples from northern vs. southern New Zealand (FST = 0.072; P = 0.0002), consistent with the hypothesis that upwelling disrupts gene flow between these regions (hypothesis 1). However, we are currently unable to rule out the alternative hypothesis that Cook Strait represents a barrier to dispersal (North Island vs. South Island; FST = 0.031; P = 0.0467). The detection of significant spatial structure in NZ samples is consistent with restricted gene flow, and the strong structure evident in northern NZ may be facilitated by distinct ocean current systems. Four shared haplotypes and nonsignificant differentiation (FST = 0.025; P = 0.2525) between southern New Zealand and Tasmanian samples is consistent with an anthropogenic origin for the latter population (hypothesis 2).

Temporal changes in allele frequencies and low effective population size in greater prairie-chickens

The number of greater prairie-chickens in Wisconsin has decreased by 91% since 1932. The current population of approximately 1500 birds exists primarily in four isolated management areas. In previous studies of the Wisconsin populations we documented low levels of genetic variation at microsatellite loci and the mitochondrial DNA control region. Here we investigate changes in genetic structure between the four management areas in Wisconsin over the last 50 years. We estimated the harmonic mean effective population size (Ne) over the last 50 years by comparing allele frequencies from the early 1950s with those from contemporary samples. Using a pseudo-likelihood approach that accounted for migration, estimates of Ne (15-32 prairie-chickens within each management area) were 10 times lower than census numbers from booming-ground counts. These low estimates of Ne are consistent with increased habitat fragmentation and an increase in genetic isolation between management areas over the last 50 years. The reduction of gene flow between areas has reduced Ne, increased genetic drift and, consequently, reduced genetic variation. These results have immediate consequences for the conservation of the prairie-chicken, and highlight the importance of how mating systems and limited dispersal may exacerbate the loss of genetic variation in fragmented populations.

Continuing the debate on the role of Quaternary environmental change for macroevolution

The Quaternary has been a period of dramatic environmental change for the past 1.8 Myr, with major shifts in distributions and abundances of terrestrial and marine organisms. The evolutionary consequences of this have been debated since the nineteenth century. However, the lack of accurate relative and absolute time-scales for evolutions and environmental change inhibited progress. We do now have an understanding of time-scales. Palaeoecology has demonstrated the individualistic nature of species' response to environmental change, but lacks a means of determining ancestry. DNA characterization of modern populations in relation to their distributions nicely complements palaeoecological results by contributing ancestry. The chance to understand how species originate and the causal factors of speciation (environmental change or otherwise) may be within reach.

Genetic characterization of Pompeii and Herculaneum Equidae buried by Vesuvius in 79 AD

DNA extracted from the skeletons of five equids discovered in a Pompeii stable and of a horse found in Herculaneum was investigated. Amino acid racemization level was consistent with the presence of DNA. Post-mortem base modifications were excluded by sequencing a 146 bp fragment of the 16S rRNA mitochondrial gene. Sequencing of a 370 bp fragment of mitochondrial (mt)DNA control region allowed the construction of a phylogenetic tree that, along with sequencing of nuclear genes (epsilon globin, gamma interferon, and p53) fragments, gave us the possibility to address some questions puzzling archaeologists. What animals-donkeys, horses, or crossbreeds-were they? And, given they had been evidently assigned to one specific job, were they all akin or were they animals with different mitochondrial haplotypes? The conclusions provided by molecular analysis show that the Pompeii remains are those of horses and mules. Furthermore one of the equids (CAV5) seems to belong to a haplotype, which is either not yet documented in the GenBank or has since disappeared. As its characteristics closely recall those of donkeys, which is the out group chosen to construct the tree, that appears to have evolved within the Equidae family much earlier than horses, this assumption seems to be nearer the truth.

Nuclear DNA sequences detect species limits in ancient moa

Ancient DNA studies have typically used multi-copy mitochondrial DNA sequences. This is largely because single-locus nuclear genes have been difficult to recover from sub-fossil material, restricting the scope of ancient DNA research. Here, we have isolated single-locus nuclear DNA markers to assign the sex of 115 extinct moa and, in combination with a mitochondrial DNA phylogeny, tested competing hypotheses about the specific status of moa taxa. Moa were large ratite birds that showed extreme size variation both within and among species. For some taxa, this large variation was hypothesized to represent sexual dimorphism, while for others it was argued to reflect the existence of different species. Our results show that moa were characterized by extreme reverse sexual dimorphism and as a result we have been able to clarify the number of moa species. For example, we show that the three recognized 'species' of Dinornis comprised only two monophyletic groups and that two of these 'species' comprised individuals of one sex only. This study also illustrates that single-locus nuclear DNA sequences can be consistently recovered from ancient material.

Epidemiology, genetic diversity, and evolution of endemic feline immunodeficiency virus in a population of wild cougars

Within the large body of research on retroviruses, the distribution and evolution of endemic retroviruses in natural host populations have so far received little attention. In this study, the epidemiology, genetic diversity, and molecular evolution of feline immunodeficiency virus specific to cougars (FIVpco) was examined using blood samples collected over several years from a free-ranging cougar population in the western United States. The virus prevalence was 58% in this population (n = 52) and increased significantly with host age. Based on phylogenetic analysis of fragments of envelope (env) and polymerase (pol) genes, two genetically distinct lineages of FIVpco were found to cooccur in the population but not in the same individuals. Within each of the virus lineages, geographically nearby isolates formed monophyletic clusters of closely related viruses. Sequence diversity for env within a host rarely exceeded 1%, and the evolution of this gene was dominated by purifying selection. For both pol and env, our data indicate mean rates of molecular evolution of 1 to 3% per 10 years. These results support the premise that FIVpco is well adapted to its cougar host and provide a basis for comparing lentivirus evolution in endemic and epidemic infections in natural hosts.

Phylogeography of the rock partridge (Alectoris graeca)

We used mitochondrial DNA control-region and microsatellite data to infer the evolutionary history and past demographic changes in 332 rock partridges (Alectoris graeca) sampled from throughout the species' distribution range, with the exception of the central Balkans region. Maternal and biparental DNA markers indicated concordantly that rock partridge populations are structured geographically (mtDNA phiST = 0.86, microsatellite FST = 0.35; RST = 0.31; P < 0.001). Phylogenetic analyses of 22 mtDNA haplotypes identified two major phylogroups (supported by bootstrap values = 93%), splitting partridges from Sicily vs. all the other sampled populations at an average Tamura-Nei genetic distance of 0.035, which corresponds to 65% of the average distance between closely related species of Alectoris. Coalescent estimates of divergence times suggested that rock partridges in Sicily were isolated for more than 200000 years. This deep subdivision was confirmed by multivariate, Bayesian clustering and population assignment analyses of microsatellite genotypes, which supported also a subdivision of partridges from the Alps vs. populations in the Apennines, Albania and Greece. Partridges in the Apennines and Albania-Greece were probably connected by gene flow since recently through a late Pleistocene Adriatic landbridge. Deglaciated Alps were probably colonized by distinct and, perhaps, not yet sampled source populations. Bottleneck and mismatch analyses indicate that rock partridges have lost variability through past population declines, and did not expand recently. Deglaciated areas could have been recolonized without any strong demographic expansion. Genetic data partially supported subspecies subdivisions, and allowed delimiting distinct conservation units. Rock partridges in Sicily, formally recognized as A. g. whitakeri, met the criteria for a distinct evolutionary significant unit.

Estimating change rates of genetic markers using serial samples: applications to the transposon IS6110 in Mycobacterium tuberculosis

In infectious disease epidemiology, it is useful to know how quickly genetic markers of pathogenic agents evolve while inside hosts. We propose a modular framework with which these genotype change rates can be estimated. The estimation scheme requires a model of the underlying process of genetic change, a detection scheme that filters this process into observable quantities, and a monitoring scheme that describes the timing of observations. We study a linear "birth-shift-death" model for change in transposable element genotypes, obtaining maximum-likelihood estimators for various detection and monitoring schemes. The method is applied to serial genotypes of the transposon IS6110 in Mycobacterium tuberculosis. The estimated birth rate of 0.0161 (events per copy of the transposon per year) and death rate of 0.0108 are both significantly larger than the estimated shift rate of 0.0018. The sum of these estimates, which corresponds to a "half-life" of 2.4 years for a typical strain that has 10 copies of the element, substantially exceeds a previous estimate of 0.0135 total changes per copy per year. We consider experimental design issues that enable the precision of estimates to be improved. We also discuss extensions to other markers and implications for molecular epidemiology.

The pedigree rate of sequence divergence in the human mitochondrial genome: there is a difference between phylogenetic and pedigree rates

We have extended our previous analysis of the pedigree rate of control-region divergence in the human mitochondrial genome. One new germline mutation in the mitochondrial DNA (mtDNA) control region was detected among 185 transmission events (generations) from five Leber hereditary optic neuropathy (LHON) pedigrees. Pooling the LHON pedigree analyses yields a control-region divergence rate of 1.0 mutation/bp/10(6) years (Myr). When the results from eight published studies that used a similar approach were pooled with the LHON pedigree studies, totaling >2,600 transmission events, a pedigree divergence rate of 0.95 mutations/bp/Myr for the control region was obtained with a 99.5% confidence interval of 0.53-1.57. Taken together, the cumulative results support the original conclusion that the pedigree divergence rate for the control region is approximately 10-fold higher than that obtained with phylogenetic analyses. There is no evidence that any one factor explains this discrepancy, and the possible roles of mutational hotspots (rate heterogeneity), selection, and random genetic drift and the limitations of phylogenetic approaches to deal with high levels of homoplasy are discussed. In addition, we have extended our pedigree analysis of divergence in the mtDNA coding region. Finally, divergence of complete mtDNA sequences was analyzed in two tissues, white blood cells and skeletal muscle, from each of 17 individuals. In three of these individuals, there were four instances in which an mtDNA mutation was found in one tissue but not in the other. These results are discussed in terms of the occurrence of somatic mtDNA mutations.

Accelerated protein evolution and origins of human-specific features: Foxp2 as an example

Genes responsible for human-specific phenotypes may have been under altered selective pressures in human evolution and thus exhibit changes in substitution rate and pattern at the protein sequence level. Using comparative analysis of human, chimpanzee, and mouse protein sequences, we identified two genes (PRM2 and FOXP2) with significantly enhanced evolutionary rates in the hominid lineage. PRM2 is a histone-like protein essential to spermatogenesis and was previously reported to be a likely target of sexual selection in humans and chimpanzees. FOXP2 is a transcription factor involved in speech and language development. Human FOXP2 experienced a >60-fold increase in substitution rate and incorporated two fixed amino acid changes in a broadly defined transcription suppression domain. A survey of a diverse group of placental mammals reveals the uniqueness of the human FOXP2 sequence and a population genetic analysis indicates possible adaptive selection behind the accelerated evolution. Taken together, our results suggest an important role that FOXP2 may have played in the origin of human speech and demonstrate a strategy for identifying candidate genes underlying the emergences of human-specific features.

The role of ecologic diversification in sibling speciation of Empidonax flycatchers (Tyrannidae): multigene evidence from mtDNA

Avian genera characterized by sibling species with distinctive habitat preferences present an evolutionary enigma in view of the more commonplace occurrence of syntopic congeners that differ strikingly in colour and pattern. No existing theory has explained the evolutionary background that led to these differences. Here we propose that great phenotypic similarity among some groups of sibling species limits their coexistence and that clues to their radiation can be seen in patterns of geographical occurrence. To illustrate our thesis we focused on the New World flycatcher genus Empidonax, a group of 15 species notorious for their great phenotypic similarity. Using 3069 base pairs of mitochondrial DNA from four genes, we produced a complete molecular phylogeny that identified four clades, three of which represent close relatives. The fourth clade includes only E. virescens, which apparently has no close living relatives. The majority of species, including many distant relatives, are completely (58.1%) or essentially (6.7%) allopatric in breeding distribution and exhibit striking ecological segregation into distinctive climate-vegetation zones. Even where ranges overlap, occupancy of the same habitat by different species is rare. Phylogenetic and distributional patterns in Empidonax suggest a peripatric model of stepwise colonization and then range expansion of small groups of pioneers during glacial periods into initially enlarging, distinctive habitats destined to be widespread during interglacials. Vicariance is not indicated in the absence of barriers of appropriate age and geographical position. Rapoport's rule that northern species have larger ranges than southern species is strongly supported.

Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data

Molecular sequences obtained at different sampling times from populations of rapidly evolving pathogens and from ancient subfossil and fossil sources are increasingly available with modern sequencing technology. Here, we present a Bayesian statistical inference approach to the joint estimation of mutation rate and population size that incorporates the uncertainty in the genealogy of such temporally spaced sequences by using Markov chain Monte Carlo (MCMC) integration. The Kingman coalescent model is used to describe the time structure of the ancestral tree. We recover information about the unknown true ancestral coalescent tree, population size, and the overall mutation rate from temporally spaced data, that is, from nucleotide sequences gathered at different times, from different individuals, in an evolving haploid population. We briefly discuss the methodological implications and show what can be inferred, in various practically relevant states of prior knowledge. We develop extensions for exponentially growing population size and joint estimation of substitution model parameters. We illustrate some of the important features of this approach on a genealogy of HIV-1 envelope (env) partial sequences.