Ancient DNA from lakes. Insights into plant response to past climate changes

Considerable advances during the last decade in DNA sequencing technology and bioinformatics have enabled in-depth study of the genome from a larger number of ancient samples derived from a larger variety of substrates and environments. In particular, extraction and analysis of ancient DNA from lake sediments has provided important clues regarding the evolution of past environments in relation to climate change. This information is essential for understanding the importance of climate and plants' dispersal ability in determining their distributional ranges, and for predicting the response of plants to ongoing climate change

Genetics and evolution of the Mediterranean Abies species

This thesis summarizes and discusses results o f five separate studies in which molecular techniques have been used to study the genetic variability and evolution o f the Abies taxa occurring in the Mediterranean region. In particular, the investigation focused on the rare species Abies nebrodensis (Lojac.) Mattei, endemic to the island of Sicily, and the three neighbouring species A. alba (Mill.), A. cephalonica (Laud.) and A. numidica (De Lann.). The main aim of the studies was to determine the amount and distribution of the genetic variability within and among Mediterranean taxa of Abies, at both the nuclear and chloroplast levels, in order to elucidate their origin and evolution and to shed light on the taxonomic position o f A. nebrodensis. In studies I, II and V allozyme markers were used to provide information on the level and distribution o f genetic variation among and within natural populations of A. alba, A. cephalonica, A. nebrodensis and A. numidica and to estimate the outcrossing rate within A. alba. In studies III and IV, DNA markers from the chloroplast genome were developed and employed at the intra- and inter-specific levels to estimate the degree of cpDNA variation in the genus and to derive inferences concerning species relationships. Two different approaches were used: the first involved a comparative restriction-site analysis of ten different amplified chloroplast DNA fragments and the second involved the analysis of six chloroplast hypervariable repetitive simple-sequence repeats (cpSSRs or microsatellites). Analyses of both allozyme and cpDNA data indicated that A. nebrodensis differs from the neighbouring Abies species and justified its classification as a separate species. Based on results from studies III, IV and V we suggested that A. nebrodensis originated during the Messinian crises of the latest part o f the Miocene, through hybridisation between A. alba and A. numidica. During the warmer period of the Holocene A. nebrodensis became more isolated from both A. numidica and A. alba, promoting further divergence. Results from studies III, IV and V agreed well with the hypothesis that the genus Abies had an Asiatic origin and then differentiated along an east-west axis (beginning in the Miocene) in the Balkan and Middle East regions. Taxa like A. alba, A. bornmuelleriana, A. nordmanniana, A. cephalonica and A. cilicica were the first to differentiate from the common Abies ancestor. Abies pinsapo and A. numidica appear to have differentiated much later, and to have remained isolated for a long time after their settlement in Spain and North Africa. The lower levels o f polymorphism observed in these two species appear to be a consequence of genetic drift due to isolation and small population size

Phylogeography and population genetics of Acanthophyllum squarrosum complex (Caryophyllaceae) in the Irano-Turanian region

Acanthophyllum squarrosum and two closely related species, A. heratense and A. laxiusculum (Caryophyllaceae), form a complex that covers parts of subalpine steppes of the Irano-Turanian (IT) region. In this study, we explored the genetic structure and phylogeography of this complex based on partial sequences of two chloroplasts (psbA–trnH and rpl32–trnL (UAG)) and two nuclear (EST24 and nrITS) DNA regions. We analysed 80 individuals from eight populations and detected 12 chloroplast haplotypes, 16 and eight nuclear alleles in EST24 and nrITS sequences, respectively. Phylogenetic trees and haplotype networks did not show distinct genetic groups in the complex and this could be explained by incomplete lineage sorting or introgression between species. Divergence time analysis revealed a Quaternary origin for A. squarrosum complex at approximately 1.8 million years ago (Mya) and the neutrality test results indicated that this complex experienced a recent population expansion. AMOVA analysis of the chloroplast regions showed a significant genetic differentiation among populations and low genetic differentiation within populations, but opposite results were found with nuclear markers, implying introgression between A. squarrosum complex populations

Population dynamics and genetic changes of Picea abies in the South Carpathians revealed by pollen and ancient DNA analyses

Background: Studies on allele length polymorphism designate several glacial refugia for Norway spruce (Picea abies) in the South Carpathian Mountains, but infer only limited expansion from these refugia after the last glaciation. To better understand the genetic dynamics of a South Carpathian spruce lineage, we compared ancient DNA from 10,700 and 11,000-year-old spruce pollen and macrofossils retrieved from Holocene lake sediment in the Retezat Mountains with DNA extracted from extant material from the same site. We used eight primer pairs that amplified short and variable regions of the spruce cpDNA. In addition, from the same lake sediment we obtained a 15,000-years-long pollen accumulation rate (PAR) record for spruce that helped us to infer changes in population size at this site. Results: We obtained successful amplifications for Norway spruce from 17 out of 462 pollen grains tested, while the macrofossil material provided 22 DNA sequences. Two fossil sequences were found to be unique to the ancient material. Population genetic statistics showed higher genetic diversity in the ancient individuals compared to the extant ones. Similarly, statistically significant Ks and Kst values showed a considerable level of differentiation between extant and ancient populations at the same loci. Lateglacial and Holocene PAR values suggested that population size of the ancient population was small, in the range of 1/10 or 1/5 of the extant population. PAR analysis also detected two periods of rapid population growths (from ca. 11,100 and 3900 calibrated years before present (cal yr BP)) and three bottlenecks (around 9180, 7200 and 2200 cal yr BP), likely triggered by climatic change and human impact. Conclusion: Our results suggest that the paternal lineages observed today in the Retezat Mountains persisted at this site at least since the early Holocene. Combination of the results from the genetic and the PAR analyses furthermore suggests that the higher level of genetic variation found in the ancient populations and the loss of ancient allele types detected in the extant individuals were likely due to the repeated bottlenecks during the Holocene. This study demonstrates how past population size changes inferred from PAR records can be efficiently used in combination with ancient DNA studies. The joint application of palaeoecological and population genetic analyses proved to be a powerful tool to understand the influence of past population demographic changes on the haplotype diversity and genetic composition of forest tree species

Population dynamics and genetic changes of Picea abies in the South Carpathians revealed by pollen and ancient DNA analyses

Background: Studies on allele length polymorphism designate several glacial refugia for Norway spruce (Picea abies) in the South Carpathian Mountains, but infer only limited expansion from these refugia after the last glaciation. To better understand the genetic dynamics of a South Carpathian spruce lineage, we compared ancient DNA from 10,700 and 11,000-year-old spruce pollen and macrofossils retrieved from Holocene lake sediment in the Retezat Mountains with DNA extracted from extant material from the same site. We used eight primer pairs that amplified short and variable regions of the spruce cpDNA. In addition, from the same lake sediment we obtained a 15,000-years-long pollen accumulation rate (PAR) record for spruce that helped us to infer changes in population size at this site.Results: We obtained successful amplifications for Norway spruce from 17 out of 462 pollen grains tested, while the macrofossil material provided 22 DNA sequences. Two fossil sequences were found to be unique to the ancient material. Population genetic statistics showed higher genetic diversity in the ancient individuals compared to the extant ones. Similarly, statistically significant Ks and Kst values showed a considerable level of differentiation between extant and ancient populations at the same loci.Lateglacial and Holocene PAR values suggested that population size of the ancient population was small, in the range of 1/10 or 1/5 of the extant population. PAR analysis also detected two periods of rapid population growths (from ca. 11,100 and 3900 calibrated years before present (cal yr BP)) and three bottlenecks (around 9180, 7200 and 2200 cal yr BP), likely triggered by climatic change and human impact.Conclusion: Our results suggest that the paternal lineages observed today in the Retezat Mountains persisted at this site at least since the early Holocene. Combination of the results from the genetic and the PAR analyses furthermore suggests that the higher level of genetic variation found in the ancient populations and the loss of ancient allele types detected in the extant individuals were likely due to the repeated bottlenecks during the Holocene; however our limited sample size did not allow us to exclude sampling effect.This study demonstrates how past population size changes inferred from PAR records can be efficiently used in combination with ancient DNA studies. The joint application of palaeoecological and population genetics analyses proved to be a powerful tool to understand the influence of past population demographic changes on the haplotype diversity and genetic composition of forest tree species

Genetic diversity and divergence at the Arbutus unedo L. (Ericaceae) westernmost distribution limit

Mediterranean forests are fragile ecosystems vulnerable to recent global warming and reduction of precipitation, and a long-term negative effect is expected on vegetation with increasing drought and in areas burnt by fires. We investigated the spatial distribution of genetic variation of Arbutus unedo in the western Iberia Peninsula, using plastid markers with conservation and provenance regions design purposes. This species is currently undergoing an intense domestication process in the region, and, like other species, is increasingly under the threat from climate change, habitat fragmentation and wildfires. We sampled 451 trees from 15 natural populations from different ecological conditions spanning the whole species' distribution range in the region. We applied Bayesian analysis and identified four clusters (north, centre, south, and a single-population cluster). Hierarchical AMOVA showed higher differentiation among clusters than among populations within clusters. The relatively low within-clusters differentiation can be explained by a common postglacial history of nearby populations. The genetic structure found, supported by the few available palaeobotanical records, cannot exclude the hypothesis of two independent A. unedo refugia in western Iberia Peninsula during the Last Glacial Maximum. Based on the results we recommend a conservation strategy by selecting populations for conservation based on their allelic richness and diversity and careful seed transfer consistent with current species' genetic structure

Norway spruce postglacial recolonization of Fennoscandia

Contrasting theories exist regarding how Norway spruce recolonized Fennoscandia after the last glaciation. Here, the authors provide evidences from sedimentary ancient DNA and modern population genomics to support that Norway spruce was present in southern Fennoscandia shortly after deglaciation and the early Holocene migration from the east. Contrasting theories exist regarding how Norway spruce (Picea abies) recolonized Fennoscandia after the last glaciation and both early Holocene establishments from western microrefugia and late Holocene colonization from the east have been postulated. Here, we show that Norway spruce was present in southern Fennoscandia as early as 14.7 +/- 0.1 cal. kyr BP and that the millennia-old clonal spruce trees present today in central Sweden likely arrived with an early Holocene migration from the east. Our findings are based on ancient sedimentary DNA from multiple European sites (N = 15) combined with nuclear and mitochondrial DNA analysis of ancient clonal (N = 135) and contemporary spruce forest trees (N = 129) from central Sweden. Our other findings imply that Norway spruce was present shortly after deglaciation at the margins of the Scandinavian Ice Sheet, and support previously disputed finds of pollen in southern Sweden claiming spruce establishment during the Lateglacial.Peer reviewe

Isolated Grauer's gorilla populations differ in diet and gut microbiome

The animal gut microbiome has been implicated in a number of key biological processes, ranging from digestion to behavior, and has also been suggested to facilitate local adaptation. Yet studies in wild animals rarely compare multiple populations that differ ecologically, which is the level at which local adaptation may occur. Further, few studies simultaneously characterize diet and gut microbiome from the same sample, despite their likely interdependence. Here, we investigate the interplay between diet and gut microbiome in three geographically isolated populations of the critically endangered Grauer's gorilla, which we show to be genetically differentiated. We find population- and social group-specific dietary and gut microbial profiles and co-variation between diet and gut microbiome, despite the presence of core microbial taxa. There was no detectable effect of age, and only marginal effects of sex and genetic relatedness on the microbiome. Diet differed considerably across populations, with the high-altitude population consuming a lower diversity of plants compared to low-altitude populations, consistent with plant availability constraining diet. The observed pattern of covariation between diet and gut microbiome is likely a result of long-term social and environmental factors. Our study suggests that the gut microbiome is sufficiently plastic to support flexible food selection and hence contribute to local adaptation

Strengthening global-change science by integrating aeDNA with paleoecoinformatics

Ancient environmental DNA (aeDNA) data are close to enabling insights into past global-scale biodiversity dynamics at unprecedented taxonomic extent and resolution. However, achieving this potential requires solutions that bridge bioinformatics and paleoecoinformatics. Essential needs include support for dynamic taxonomic inferences, dynamic age inferences, and precise stratigraphic depth. Moreover, aeDNA data are complex and heterogeneous, generated by dispersed researcher networks, with methods advancing rapidly. Hence, expert community governance and curation are essential to building high-value data resources. Immediate recommendations include uploading metabarcoding-based taxonomic inventories into paleoecoinformatic resources, building linkages among open bioinformatic and paleoecoinformatic data resources, harmonizing aeDNA processing workflows, and expanding community data governance. These advances will enable transformative insights into global-scale biodiversity dynamics during large environmental and anthropogenic changes