Improving ancient DNA read mapping against modern reference genomes

ABSTRACT: 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.     

The history of the plant embryo. Terminology and visualization from ancient until modern times

Since ancient times comparisons between embryonic forms of humans, animals, and plants are known. In deciphering a plant embryo and its development, one applied a specific zoomorphic terminology. Until the 17th century naturalists who studied plants were inspired by the concepts of ancient natural philosophy. Since then plant embryos are visualized by drawings and diagrammatic sketches. In the 18th century the embryo became an important issue in debates concerning theories of generation and the analogy between animal egg and vegetable seed was emphasized. Due to the cell theory and refined microscopic techniques around 1850 botanists described the 'plant embryo' as an aggregate of cells. The 20th century profited of further technical improvements, e.g. microphotography, SEM, and plant tissue culture. The spell of the word embryo integrated morphological and morphogenetic data in botany and served as a code to decipher botanic forms. The paper presents a short history of the concept of 'plant embryo' and focuses on its distinct meanings over centuries, supported by different images.     

Reanalysis of Eurasian population history: ancient DNA evidence of population affinities

Mitochondrial hypervariable region I genetic data from ancient populations at two sites in Asia-Linzi in Shandong (northern China) and Egyin Gol in Mongolia-were reanalyzed to detect population affinities. Data from 51 modern populations were used to generate distance measures (Fst's) to the two ancient populations. The tests first analyzed relationships at the regional level and then compiled the top regional matches for an overall comparison to the two probe populations. The reanalysis showed that the Egyin Gol and Linzi populations have clear distinctions in genetic affinity. The Egyin Gol population as a whole appears to bear close affinities with modern populations of northern East Asia. The Linzi population seems to have some genetic affinities with the West, as suggested by the original analysis, although the original attribution of "European-like" seems to be misleading. We suggest that the Linzi individuals are potentially related to early Iranians, who are thought to have been widespread in parts of Central Eurasia and the steppe regions in the first millennium B.C., although some significant admixture between a number of populations of varying origin cannot be ruled out. We also examine the effect of sequence length on this type of genetic data analysis and discuss the results of previous studies on the Linzi sample.     

Pathocoenoses ancient and modern

Problems concerning the concept of biocoenosis in ecology (the antecedent of the pathocoenosis concept) are discussed first of all. Six main problems are identified: the problem of emergent properties of ecological communities; the problem of ambiguity; the problem of heterogeneity; the boundary problem; the problem of retrospective differential diagnosis; and the problem of explaining change over time. The rest of the paper gives illustrations of these problems in relation to human pathogens drawn mainly from the interactions of malaria with other diseases, particularly but not exclusively in the Mediterranean world, from antiquity through to modern times.     

Reconstructing the evolutionary history of China: a caveat about inferences drawn from ancient DNA

The decipherment of the meager information provided by short fragments of ancient mitochondrial DNA (mtDNA) is notoriously difficult but is regarded as a most promising way toward reconstructing the past from the genetic perspective. By haplogroup-specific hypervariable segment (HVS) motif search and matching or near-matching with available modern data sets, most of the ancient mtDNAs can be tentatively assigned to haplogroups, which are often subcontinent specific. Further typing for mtDNA haplogroup-diagnostic coding region polymorphisms, however, is indispensable for establishing the geographic/genetic affinities of ancient samples with less ambiguity. In the present study, we sequenced a fragment (approximately 982 bp) of the mtDNA control region in 76 Han individuals from Taian, Shandong, China, and we combined these data with previously reported samples from Zibo and Qingdao, Shandong. The reanalysis of two previously published ancient mtDNA population data sets from Linzi (same province) then indicates that the ancient populations had features in common with the modern populations from south China rather than any specific affinity to the European mtDNA pool. Our results highlight that ancient mtDNA data obtained under different sampling schemes and subject to potential contamination can easily create the impression of drastic spatiotemporal changes in the genetic structure of a regional population during the past few thousand years if inappropriate methods of data analysis are employed.     

Classical catalase: ancient and modern

This review describes the historical difficulties in devising a kinetically satisfactory mechanism for the classical catalase after its identification as a unique catalytic entity in 1902 and prior to the breakthrough 1947 analysis by Chance and co-workers which led to the identification of peroxide compounds I and II. The role of protons in the formation of these two ferryl complexes is discussed and current problems of inhibitory ligand and hydrogen donor binding at the active site are outlined, especially the multiple roles involving formate or formic acid. A previous mechanism of NADPH-dependent catalase protection against substrate inhibition is defended. A revised model linking the catalytic ('catalatic') action and the one-electron side reactions involving compound II is suggested. And it is concluded that, contrary to an idea proposed in 1963 that eukaryotic catalase might be a 'fossil enzyme', current thinking gives it a central role in the redox protective processes of long term importance for human and other eukaryotic and prokaryotic life.     

Mitochondrial genomes from major lizard families suggest their phylogenetic relationships and ancient radiations

In placental mammals and birds, molecular data generally support a view that they diverged into their ordinal groups in good response to mid-Cretaceous continental fragmentations. However, such divergence patterns have rarely been studied for reptiles for which phylogenetic relationships among their major groups have not yet been established molecularly. Here, I determined complete or nearly complete mitochondrial DNA sequences from seven lizard families and reconstructed phylogenetic relationships between major lizard families. When snakes were included, maximum likelihood analysis did not support a morphological view of the snakes-varanoids affinity, although several other competing hypotheses on the position of snakes still cannot be discriminated presumably due to extremely long branches of the snake lineages. I also conducted clock-free Bayesian analyses to show that divergence times between major lizard families were centered in Triassic-Jurassic times. Thus, lizards include much deeper divergences than the mammals and birds and they appear to have already radiated into various families prior to the mid-Cretaceous major continental fragmentation.     

Science of breeding and heredity from ancient Persia to modern Iran

About 1700 years BC, the prophet Zoroaster declared equal right for women and men to choose their "own ways." There is much evidence that ancient Persians believed in the equal contribution of women and men toward producing a child, and all its hereditary characteristics.Even more surprising are the phrases in Vandidad book, which were gathered by Mobedans in the Mad dynasty about egg extraction (gametes) from animal reproductive organs (gonads) and their storage for future conception.Centuries later, Western philosopher beliefs in regard to reproduction were contrary to Persian knowledge. The Greek philosophers believed that man's water (semen) contains all human characteristics, and the female uterus is only responsible for nurturing and development of fetus. After detection of the ovum (de Graaf 2(nd) half 17 century) Malpigy proposed the preformation theory (ovist) which means there is a miniature human inside ovum, that grows after Semen has entered the uterus and grow into a well-developed fetus. This hypothesis was later delegated to spermatozoa. These contradictory and inappropriate beliefs were subject to discussions and dispute, until C.E. Wolf demonstrated that the embryo is a product of the fertilization of ovum by spermatozoa.800 years prior this the sage Ferdowsi "The Great Iranian Poet" explains nicely the equal participation of man and woman in the production of the fetus and transmission of characters.After the renaissance and especially in recent years, tremendous achievements have been made in unraveling biological secrets of reproduction. There was no work o n genetics in Iran until 1936, when a genetic course was added to the biology curriculum in related colleges and universities; Iranian Genetics Society was founded in 1966, initiating a steady movement in this field.Although there was an inevitable gap during the revolution and war in our country, now there is great effort by researchers to eliminate the gap and bring us into the mainstream of world science, and development in biomedical sciences in the third millennium.     

Deciphering the ancient and complex evolutionary history of human arylamine N-acetyltransferase genes

The human N-acetyltransferase genes NAT1 and NAT2 encode two phase-II enzymes that metabolize various drugs and carcinogens. Functional variability at these genes has been associated with adverse drug reactions and cancer susceptibility. Mutations in NAT2 leading to the so-called slow-acetylation phenotype reach high frequencies worldwide, which questions the significance of altered acetylation in human adaptation. To investigate the role of population history and natural selection in shaping NATs variation, we characterized genetic diversity through the resequencing and genotyping of NAT1, NAT2, and the pseudogene NATP in a collection of 13 different populations with distinct ethnic backgrounds and demographic pasts. This combined study design allowed us to define a detailed map of linkage disequilibrium of the NATs region as well as to perform a number of sequence-based neutrality tests and the long-range haplotype (LRH) test. Our data revealed distinctive patterns of variability for the two genes: the reduced diversity observed at NAT1 is consistent with the action of purifying selection, whereas NAT2 functional variation contributes to high levels of diversity. In addition, the LRH test identified a particular NAT2 haplotype (NAT2*5B) under recent positive selection in western/central Eurasians. This haplotype harbors the mutation 341T-->C and encodes the "slowest-acetylator" NAT2 enzyme, suggesting a general selective advantage for the slow-acetylator phenotype. Interestingly, the NAT2*5B haplotype, which seems to have conferred a selective advantage during the past approximately 6,500 years, exhibits today the strongest association with susceptibility to bladder cancer and adverse drug reactions. On the whole, the patterns observed for NAT2 well illustrate how geographically and temporally fluctuating xenobiotic environments may have influenced not only our genome variability but also our present-day susceptibility to disease.     

Environment, population, and biology: a short history of modern epidemiology

Since its origin in the 19th century, epidemiology has faced an internal tension between an approach oriented toward biology and the study of mechanisms, and an approach oriented toward populations and their interactions with the environment. Initially, this tension took the form of an opposition between microbiology and statistics. We describe the early roots of the quantitative approach to health and disease and several historical examples of the above tension. The search for the causes of pellagra exemplifies our thesis. In Italy, where pellagra was endemic, contrasting opinions coexisted between the hypothesis of contaminated maize, supported by Cesare Lombroso, and the hypothesis of a prevailing role of poverty and poor nutrition. In the United States, Joseph Goldberger found no evidence for the hypothesis of contaminated maize or for a microbiological agent, but recognized the central role of nutrition. The "cure" Goldberger proposed was land reform, but he continued studying the disease from a mechanistic point of view; shortly after his death, niacin deficiency was identified as the cause of pellagra. The tension between mechanistic and population-based studies is still present within epidemiology and is in fact essential for the success of the discipline.     

Evolutionary landscape of amphibians emerging from ancient freshwater fish inferred from complete mitochondrial genomes

It is very interesting that the only extant marine amphibian is the marine frog, Fejervarya cancrivora. This study investigated the reasons for this apparent rarity by conducting a phylogenetic tree analysis of the complete mitochondrial genomes from 14 amphibians, 67 freshwater fishes, four migratory fishes, 35 saltwater fishes, and one hemichordate. The results showed that amphibians, living fossil fishes, and the common ancestors of modern fishes are phylogenetically separated. In general, amphibians, living fossil fishes, saltwater fishes, and freshwater fishes are clustered in different clades. This suggests that the ancestor of living amphibians arose from a type of primordial freshwater fish, rather than the coelacanth, lungfish, or modern saltwater fish. Modern freshwater fish and modern saltwater fish were probably separated from a common ancestor by a single event, caused by crustal movement.     

Molecular phylogenetics employing modern and ancient DNA

Comparative studies of DNA in recent populations and characterisation of ancient hereditary material have contributed very interesting facts to our understanding of evolution of modern mankind. Analysis of DNA homology in related species, assessment of mutations and polymorphisms in various populations and new DNA sequence data from prehistoric finds allowed - via sophisticated DNA extraction techniques, PCR, sequencing and digitalised processing of genetic information - insights into possible roots of Homo sapiens and related species, migration patterns and ancient cultural habits, thus enrhing the palaeoanthropological discipline. However, a presentation of this development would not be complete without pointing towards the methodological limitations and manifold presentations burdened with artifacts, data misinterpretation and unjustified conclusions. Presently, this modern field of research is in its consolidation phase and new parameters for quality control and authentication are being implemented to avoid spectacular but unfounded reports. It is expected that most of the problems connected to old biomolecules may be closely related to fossilisation parameters. The future challenge will be the full understanding of the complex and multi-faceted processes underlying diagenesis, including the elucidation of nucleic acid postmortem damage".     

Inferring population history from genealogical trees

Inference about population history from DNA sequence data has become increasingly popular. For human populations, questions about whether a population has been expanding and when expansion began are often the focus of attention. For viral populations, questions about the epidemiological history of a virus, e.g., HIV-1 and Hepatitis C, are often of interest. In this paper I address the following question: Can population history be accurately inferred from single locus DNA data? An idealised world is considered in which the tree relating a sample of n non-recombining and selectively neutral DNA sequences is observed, rather than just the sequences themselves. This approach provides an upper limit to the information that possibly can be extracted from a sample. It is shown, based on Kingman's (1982a) coalescent process, that consistent estimation of parameters describing population history (e.g., a growth rate) cannot be achieved for increasing sample size, n. This is worse than often found for estimators of genetic parameters, e.g., the mutation rate typically converges at rate \(\) under the assumption that all historical mutations can be observed in the sample. In addition, various results for the distribution of maximum likelihood estimators are presented.     

Human pelvis and long bones reveal differential preservation of ancient population history and migration out of Africa

One of the main events in the history of our species has been our expansion out of Africa. A clear signature of this expansion has been found on global patterns of neutral genetic variation, whereby a serial founder effect accompanied the colonization of new regions, in turn creating a wilhin-pupulation decrease in neutral genetic diversity with increasing distance from Africa. This same distinctive pattern has also been described for cranial and dental morphological variation in human populations distributed across the globe. Here, we used a data set of postcranial linear measurements for 30 globally distributed human populations, and a climatic data set of minimum annual temperature, maximum annual temperature, and precipitation in order to separate for the first time the relative effect of neutral demographic processes and climatic selection on four long (limb) bones (femur, tibia, radius, and humerus) versus the pelvic bones of the human appendicular skeleton. We implemented a stepwise regression procedure in which phenotypic variance is assumed to be affected by the iterative founder events that accompanied human expansion from Africa, as well as by climate. This model included, as independent factors, geographic distance from central Africa, the three climatic variables, and all possible interactions between the three climatic variables. We excluded all nonsignificant factors by backward stepwise elimination with the aim of identifying the minimal model significantly explaining variation in the phenotypic data. Our results indicate a sharp difference in the way the pelvis and the limb bones reflect the neutral signature of the out-of-Africa expansion. Consistent with previous analyses of the cranium and dentition, pelvic shape variation shows a significant within-population decrease with increasing distance from Africa. However, no such pattern could be found in the long bones. Rather, in the case of both the tibia and the femur, a significant relationship between population-level variance and minimum temperature was demonstrated. Hence, in the case of these limb bones, it is probable that the effects of climatic selection have obliterated the demographic signature of human dispersal from Africa. Our finding mat pelvic variation exhibits the neutral effects of demographic history suggests that consideration of this skeletal element might be used to shed light on factors of human population history, just as the cranium has done.