Katerina Guschanski

Temporal patterns of nucleotide misincorporations and DNA fragmentation in ancient DNA

DNA that survives in museum specimens, bones and other tissues recovered by archaeologists is invariably fragmented and chemically modified. The extent to which such modifications accumulate over time is largely unknown but could potentially be used to differentiate between endogenous old DNA and present-day DNA contaminating specimens and experiments. Here we examine mitochondrial DNA sequences from tissue remains that vary in age between 18 and 60,000 years with respect to three molecular features: fragment length, base composition at strand breaks, and apparent C to T substitutions. We find that fragment length does not decrease consistently over time and that strand breaks occur preferentially before purine residues by what may be at least two different molecular mechanisms that are not yet understood. In contrast, the frequency of apparent C to T substitutions towards the 5′-ends of molecules tends to increase over time. These nucleotide misincorporations are thus a useful tool to distinguish recent from ancient DNA sources in specimens that have not been subjected to unusual or harsh treatments.

Birth and expression evolution of mammalian microRNA genes

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression, yet their origins and functional evolution in mammals remain little understood due to the lack of appropriate comparative data. Using RNA sequencing, we have generated extensive and comparable miRNA data for five organs in six species that represent all main mammalian lineages and birds (the evolutionary outgroup) with the aim to unravel the evolution of mammalian miRNAs. Our analyses reveal an overall expansion of miRNA repertoires in mammals, with threefold accelerated birth rates of miRNA families in placentals and marsupials, facilitated by the de novo emergence of miRNAs in host gene introns. Generally, our analyses suggest a high rate of miRNA family turnover in mammals with many newly emerged miRNA families being lost soon after their formation. Selectively preserved mammalian miRNA families gradually evolved higher expression levels, as well as altered mature sequences and target gene repertoires, and were apparently mainly recruited to exert regulatory functions in nervous tissues. However, miRNAs that originated on the X chromosome evolved high expression levels and potentially diverse functions during spermatogenesis, including meiosis, through selectively driven duplication-divergence processes. Overall, our study thus provides detailed insights into the birth and evolution of mammalian miRNA genes and the associated selective forces.

Two‐step multiplex polymerase chain reaction improves the speed and accuracy of genotyping using DNA from noninvasive and museum samples

Many studies in molecular ecology rely upon the genotyping of large numbers of low‐quantity DNA extracts derived from noninvasive or museum specimens. To overcome low amplification success rates and avoid genotyping errors such as allelic dropout and false alleles, multiple polymerase chain reaction (PCR) replicates for each sample are typically used. Recently, two‐step multiplex procedures have been introduced which drastically increase the success rate and efficiency of genotyping. However, controversy still exists concerning the amount of replication needed for suitable control of error. Here we describe the use of a two‐step multiplex PCR procedure that allows rapid genotyping using at least 19 different microsatellite loci. We applied this approach to quantified amounts of noninvasive DNAs from western chimpanzee, western gorilla, mountain gorilla and black and white colobus faecal samples, as well as to DNA from ~100‐year‐old gorilla teeth from museums. Analysis of over 45 000 PCRs revealed average success rates of > 90% using faecal DNAs and 74% using museum specimen DNAs. Average allelic dropout rates were substantially reduced compared to those obtained using conventional singleplex PCR protocols, and reliable genotyping using low (< 25 pg) amounts of template DNA was possible. However, four to five replicates of apparently homozygous results are needed to avoid allelic dropout when using the lowest concentration DNAs (< 50 pg/reaction), suggesting that use of protocols allowing routine acceptance of homozygous genotypes after as few as three replicates may lead to unanticipated errors when applied to low‐concentration DNAs.

Next-Generation Museomics Disentangles One of the Largest Primate Radiations

Guenons (tribe Cercopithecini) are one of the most diverse groups of primates. They occupy all of sub-Saharan Africa and show great variation in ecology, behavior, and morphology. This variation led to the description of over 60 species and subspecies. Here, using next-generation DNA sequencing (NGS) in combination with targeted DNA capture, we sequenced 92 mitochondrial genomes from museum-preserved specimens as old as 117 years. We infer evolutionary relationships and estimate divergence times of almost all guenon taxa based on mitochondrial genome sequences. Using this phylogenetic framework, we infer divergence dates and reconstruct ancestral geographic ranges. We conclude that the extraordinary radiation of guenons has been a complex process driven by, among other factors, localized fluctuations of African forest cover. We find incongruences between phylogenetic trees reconstructed from mitochondrial and nuclear DNA sequences, which can be explained by either incomplete lineage sorting or hybridization. Furthermore, having produced the largest mitochondrial DNA data set from museum specimens, we document how NGS technologies can “unlock” museum collections, thereby helping to unravel the tree-of-life. [Museum collection; next-generation DNA sequencing; primate radiation; speciation; target capture.]

Historical sampling reveals dramatic demographic changes in western gorilla populations

BackgroundToday many large mammals live in small, fragmented populations, but it is often unclear whether this subdivision is the result of long-term or recent events. Demographic modeling using genetic data can estimate changes in long-term population sizes while temporal sampling provides a way to compare genetic variation present today with that sampled in the past. In order to better understand the dynamics associated with the divergences of great ape populations, these analytical approaches were applied to western gorillas (Gorilla gorilla) and in particular to the isolated and Critically Endangered Cross River gorilla subspecies (G. g. diehli).ResultsWe used microsatellite genotypes from museum specimens and contemporary samples of Cross River gorillas to infer both the long-term and recent population history. We find that Cross River gorillas diverged from the ancestral western gorilla population ~17,800 years ago (95% HDI: 760, 63,245 years). However, gene flow ceased only ~420 years ago (95% HDI: 200, 16,256 years), followed by a bottleneck beginning ~320 years ago (95% HDI: 200, 2,825 years) that caused a 60-fold decrease in the effective population size of Cross River gorillas. Direct comparison of heterozygosity estimates from museum and contemporary samples suggests a loss of genetic variation over the last 100 years.ConclusionsThe composite history of western gorillas could plausibly be explained by climatic oscillations inducing environmental changes in western equatorial Africa that would have allowed gorilla populations to expand over time but ultimately isolate the Cross River gorillas, which thereafter exhibited a dramatic population size reduction. The recent decrease in the Cross River population is accordingly most likely attributable to increasing anthropogenic pressure over the last several hundred years. Isolation of diverging populations with prolonged concomitant gene flow, but not secondary admixture, appears to be a typical characteristic of the population histories of African great apes, including gorillas, chimpanzees and bonobos.

Females Shape the Genetic Structure of a Gorilla Population

Dispersal, one of the key life-history features of a species, influences gene flow and, consequently, the genetic structuring of populations. Landscape characteristics such as rivers, mountains, or habitat fragmentation affect dispersal and result in broad-scale genetic structuring of various mammalian species [1-5]. However, less attention has been paid to studying how dispersal is influenced by finer-scale microgeographic variation in a continuous habitat. Here we investigate the genetic structure of a closed population of approximately 300 endangered mountain gorillas living in multiple groups in a small (331 km(2)) forest in southwestern Uganda. In a species in which both sexes routinely disperse, population genetic structure in females was influenced by distance, altitude, and plant community composition, whereas males were not geographically structured. The effect of distance fits the observed tendency of females to transfer to neighboring groups, whereas the effects of altitude and vegetation reflect the changing species composition of locally available food resources. These results suggest that individual dietary preferences are important even in a highly mobile species living amid abundant food, and we propose that preference for natal habitats will influence dispersal decisions in many other vertebrate taxa.

Using genetics to understand the dynamics of wild primate populations

While much can be learned about primates by means of observation, the slow life history of many primates means that even decades of dedicated effort cannot illuminate long-term evolutionary processes. For example, while the size of a contemporary population can be estimated from field censuses, it is often desirable to know whether a population has been constant or changing in size over a time frame of hundreds or thousands of years. Even the nature of “a population” is open to question, and the extent to which individuals successfully disperse among defined populations is also difficult to estimate by using observational methods alone. Researchers have thus turned to genetic methods to examine the size, structure, and evolutionary histories of primate populations. Many results have been gained by study of sequence variation of maternally inherited mitochondrial DNA, but in recent years researchers have been increasingly focusing on analysis of short, highly variable microsatellite segments in the autosomal genome for a high-resolution view of evolutionary processes involving both sexes. In this review we describe some of the insights thus gained, and discuss the likely impact on this field of new technologies such as high-throughput DNA sequencing.

Recent divergences and size decreases of eastern gorilla populations

Compared with other African apes, eastern gorillas (Gorilla beringei) have been little studied genetically. We used analysis of autosomal DNA genotypes obtained from non-invasively collected faecal samples to estimate the evolutionary histories of the two extant mountain gorilla populations and the closely related eastern lowland gorillas. Our results suggest that eastern lowland gorillas and mountain gorillas split beginning some 10 000 years ago, followed 5000 years ago by the split of the two mountain gorilla populations of Bwindi Impenetrable National Park and the Virungas Massif. All three populations have decreased in effective population size, with particularly substantial 10-fold decreases for the mountain gorillas. These dynamics probably reflect responses to habitat changes resulting from climate fluctuations over the past 20 000 years as well as increasing human influence in this densely populated region in the last several thousand years.

The evolution of duplicate gene expression in mammalian organs

Gene duplications generate genomic raw material that allows the emergence of novel functions, likely facilitating adaptive evolutionary innovations. However, global assessments of the functional and evolutionary relevance of duplicate genes in mammals were until recently limited by the lack of appropriate comparative data. Here, we report a large-scale study of the expression evolution of DNA-based functional gene duplicates in three major mammalian lineages (placental mammals, marsupials, egg-laying monotremes) and birds, on the basis of RNA sequencing (RNA-seq) data from nine species and eight organs. We observe dynamic changes in tissue expression preference of paralogs with different duplication ages, suggesting differential contribution of paralogs to specific organ functions during vertebrate evolution. Specifically, we show that paralogs that emerged in the common ancestor of bony vertebrates are enriched for genes with brain-specific expression and provide evidence for differential forces underlying the preferential emergence of young testis- and liver-specific expressed genes. Further analyses uncovered that the overall spatial expression profiles of gene families tend to be conserved, with several exceptions of pronounced tissue specificity shifts among lineage-specific gene family expansions. Finally, we trace new lineage-specific genes that may have contributed to the specific biology of mammalian organs, including the little-studied placenta. Overall, our study provides novel and taxonomically broad evidence for the differential contribution of duplicate genes to tissue-specific transcriptomes and for their importance for the phenotypic evolution of vertebrates.

Whole mitochondrial genome capture from faecal samples and museum-preserved specimens

Population‐scale molecular studies of endangered and cryptic species are often limited by access to high‐quality samples. The use of noninvasively collected samples or museum‐preserved specimens reduces the pressure on modern populations by removing the need to capture and handle live animals. However, endogenous DNA content in such samples is low, making shotgun sequencing a financially prohibitive approach. Here, we apply a target enrichment method to retrieve mitochondrial genomes from 65 museum specimens and 56 noninvasively collected faecal samples of two endangered great ape species, Grauers gorilla and the eastern chimpanzee. We show that the applied method is suitable for a wide range of sample types that differ in endogenous DNA content, increasing the proportion of target reads to over 300‐fold. By systematically evaluating biases introduced during target enrichment of pooled museum samples, we show that capture is less efficient for fragments shorter or longer than the baits, that the proportion of human contaminating reads increases postcapture although capture efficiency is lower for human compared to gorilla fragments with a gorilla‐generated bait, and that the rate of jumping PCR is considerable, but can be controlled for with a double‐barcoding approach. We succeed in capturing complete mitochondrial genomes from faecal samples, but observe reduced capture efficiency as sequence divergence increases between the bait and target species. As previously shown for museum specimens, we demonstrate here that mitochondrial genome capture from field‐collected faecal samples is a robust and reliable approach for population‐wide studies of nonmodel organisms.