Ancient DNA from the first European farmers in 7500-year-old Neolithic sites

Additive partitioning of multispecies distributional aggregation of local assemblages

The distribution of species is not random in space. At the finest-resolution spatial scale, that is, field sampling locations, distributional aggregation level of different species would be determined by various factors, for example spatial autocorrelation or environmental filtering. However, few studies have quantitatively measured the importance of these factors. In this study, inspired by the statistical properties of a Markov transition model, we propose a novel additive framework to partition local multispecies distributional aggregation levels for sequential sampling-derived field biodiversity data. The framework partitions the spatial distributional aggregation of different species into two independent components: regional abundance variability and the local spatial inertia effect. Empirical studies from field amphibian surveys through line-transect sampling in southwestern China (Minya Konka) and central-southern Vietnam showed that local spatial inertia was always the dominant mechanism structuring the local occurrence and distributional aggregation of amphibians in the two regions with a latitudinal gradient from 1200 to nearly 4000 m. However, regional abundance variability is still nonnegligible in highly diverse tropical regions (i.e. Vietnam) where the altitude is not higher than 2000 m. In summary, we propose a novel framework that shows that the multispecies distributional aggregation level can be structured by two additive components. The two partitioned components could be theoretically independent. These findings are expected to deepen our understanding of the local community structure from the perspective of both spatial distribution and regional diversity patterns. The partitioning framework might have potential applications in field ecology and macroecology research.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity

STUDY QUESTION

Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins?

SUMMARY ANSWER

We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3.

WHAT IS KNOWN ALREADY

The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures.

STUDY DESIGN, SIZE, DURATION

We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation.

MAIN RESULTS AND THE ROLE OF CHANCE

This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle.

LARGE SCALE DATA

The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/).

LIMITATIONS, REASONS FOR CAUTION

Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility.

WIDER IMPLICATIONS OF THE FINDINGS

About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility.

STUDY FUNDING/COMPETING INTEREST(S)

Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Skłodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests.

TRIAL REGISTRATION NUMBER

N/A.

Reconstructing the post-glacial spread of the sand fly Phlebotomus mascittii Grassi, 1908 (Diptera: Psychodidae) in Europe

Phlebotomine sand flies (Diptera: Phlebotominae) are the principal vectors of Leishmania spp. (Kinetoplastida: Trypanosomatidae). In Central Europe, Phlebotomus mascittii is the predominant species, but largely understudied. To better understand factors driving its current distribution, we infer patterns of genetic diversity by testing for signals of population expansion based on two mitochondrial genes and model current and past climate and habitat suitability for seven post-glacial maximum periods, taking 19 climatic variables into account. Consequently, we elucidate their connections by environmental-geographical network analysis. Most analyzed populations share a main haplotype tracing back to a single glacial maximum refuge area on the Mediterranean coasts of South France, which is supported by network analysis. The rapid range expansion of Ph. mascittii likely started in the early mid-Holocene epoch until today and its spread possibly followed two routes. The first one was through northern France to Germany and then Belgium, and the second across the Ligurian coast through present-day Slovenia to Austria, toward the northern Balkans. Here we present a combined approach to reveal glacial refugia and post-glacial spread of Ph. mascittii and observed discrepancies between the modelled and the current known distribution might reveal yet overlooked populations and potential further spread.

Refining the Global Phylogeny of Mitochondrial N1a, X, and HV2 Haplogroups Based on Rare Mitogenomes from Croatian Isolates

Mitochondrial DNA (mtDNA) has been used for decades as a predominant tool in population genetics and as a valuable addition to forensic genetic research, owing to its unique maternal inheritance pattern that enables the tracing of individuals along the maternal lineage across numerous generations. The dynamic interplay between evolutionary forces, primarily genetic drift, bottlenecks, and the founder effect, can exert significant influence on genetic profiles. Consequently, the Adriatic islands have accumulated a subset of lineages that exhibits remarkable absence or rarity within other European populations. This distinctive genetic composition underscores the islands' potential as a significant resource in phylogenetic research, with implications reaching beyond regional boundaries to contribute to a global understanding. In the initial attempt to expand the mitochondrial forensic database of the Croatian population with haplotypes from small isolated communities, we sequenced mitogenomes of rare haplogroups from different Croatian island and mainland populations using next-generation sequencing (NGS). In the next step and based on the obtained results, we refined the global phylogeny of haplogroup N1a, HV2, and X by analyzing rare haplotypes, which are absent from the current phylogenetic tree. The trees were based on 16 novel and 52 previously published samples, revealing completely novel branches in the X and HV2 haplogroups and a new European cluster in the ancestral N1a variant, previously believed to be an exclusively African-Asian haplogroup. The research emphasizes the importance of investigating geographically isolated populations and their unique characteristics within a global context.

Methodological Changes in the Field of Paleogenetics

Paleogenetics has significantly changed since its inception almost forty years ago. Initially, molecular techniques available to the researchers offered minimal possibilities for ancient DNA analysis. The subsequent expansion of the scientific tool cabinet allowed for more remarkable achievements, combined has with the newfound popularity of this budding field of science. Finally, a breakthrough was made with the development of next-generation sequencing (NGS) technologies and the update of DNA isolation protocols, through which even very fragmented aDNA samples could be used to sequence whole genomes. In this paper, we review the achievements made thus far and compare the methodologies utilized in this field of science, discussing their benefits and challenges.

Environmental DNA as an innovative technique to identify the origins of falsified antimalarial tablets-a pilot study of the pharmabiome

Falsified medicines are a major threat to global health. Antimalarial drugs have been particularly targeted by criminals. As DNA analysis has revolutionized forensic criminology, we hypothesized that these techniques could also be used to investigate the origins of falsified medicines. Medicines may contain diverse adventitious biological contamination, and the sealed nature of blister-packages may capture and preserve genetic signals from the manufacturing processes allowing identification of production source(s). We conducted a blinded pilot study to determine if such environmental DNA (eDNA) could be detected in eleven samples of falsified and genuine artesunate antimalarial tablets, collected in SE Asia, which could be indicative of origin. Massively Parallel Sequencing (MPS) was used to characterize microbial and eukaryote diversity. Two mitochondrial DNA analysis approaches were explored to detect the presence of human DNA. Trace eDNA from these low biomass samples demonstrated sample specific signals using two target markers. Significant differences in bacterial and eukaryote DNA community structures were observed between genuine and falsified tablets and between different packaging types of falsified artesunate. Human DNA, which was indicative of likely east Asian ancestry, was found in falsified tablets. This pilot study of the 'pharmabiome' shows the potential of environmental DNA as a powerful forensic tool to assist with the identification of the environments, and hence location and timing, of the source and manufacture of falsified medicines, establish links between seizures and complement existing tools to build a more complete picture of criminal trade routes. The finding of human DNA in tablets raises important ethical issues that need to be addressed.

Variation in the substitution rates among the human mitochondrial haplogroup U sublineages

Resolving the absolute timescale of phylogenetic trees stipulates reliable estimates for the rate of DNA sequence evolution. For this end, various calibration methods have been developed and studied intensively. Intraspecific rate variation among distinct genetic lineages, however, has gained less attention. Here, we have assessed lineage-specific molecular rates of human mitochondrial DNA (mtDNA) by performing tip-calibrated Bayesian phylogenetic analyses. Tip-calibration, as opposed to traditional nodal time stamps from dated fossil evidence or geological events, is based on sample ages and becoming ever more feasible as ancient DNA data from radiocarbon-dated samples accumulate. We focus on subhaplogroups U2, U4, U5a, and U5b, the data including ancient mtDNA genomes from ¹⁴C-dated samples (n = 234), contemporary genomes (n = 301), and two outgroup sequences from haplogroup R. The obtained molecular rates depended on the data sets (with or without contemporary sequences), suggesting time-dependency. More notable was the rate variation between haplogroups: U4 and U5a stand out having a substantially higher rate than U5b. This is also reflected in the divergence times obtained (U5a: 17,700 years and U5b: 29,700 years), a disparity not reported previously. After ruling out various alternative causes (e.g., selection, sampling, and sequence quality), we propose that the substitution rates have been influenced by demographic histories, widely different among populations where U4/U5a or U5b are frequent. As with the Y-chromosomal subhaplogroup R1b, the mitochondrial U4 and U5a have been associated with remarkable range extensions of the Yamnaya culture in the Bronze Age.

Biodiversity underpins fisheries resilience to exploitation in the Amazon river basin

Inland fisheries feed greater than 150 million people globally, yet their status is rarely assessed due to their socio-ecological complexity and pervasive lack of data. Here, we leverage an unprecedented landings time series from the Amazon, Earth's largest river basin, together with theoretical food web models to examine (i) taxonomic and trait-based signatures of exploitation in inland fish landings and (ii) implications of changing biodiversity for fisheries resilience. In both landings time series and theory, we find that multi-species exploitation of diverse inland fisheries results in a hump-shaped landings evenness curve. Along this trajectory, abundant and large species are sequentially replaced with faster growing and smaller species. Further theoretical analysis indicates that harvests can be maintained for a period of time but that continued biodiversity depletion reduces the pool of compensating species and consequently diminishes fisheries resilience. Critically, higher fisheries biodiversity can delay fishery collapse. Although existing landings data provide an incomplete snapshot of long-term dynamics, our results suggest that multi-species exploitation is affecting freshwater biodiversity and eroding fisheries resilience in the Amazon. More broadly, we conclude that trends in landings evenness could characterize multi-species fisheries development and aid in assessing their sustainability.

Evolution and dispersal of mitochondrial DNA haplogroup U5 in Northern Europe: insights from an unsupervised learning approach to phylogeography

Background

We combined an unsupervised learning methodology for analyzing mitogenome sequences with maximum likelihood (ML) phylogenetics to make detailed inferences about the evolution and diversification of mitochondrial DNA (mtDNA) haplogroup U5, which appears at high frequencies in northern Europe.

Methods

Haplogroup U5 mitogenome sequences were gathered from GenBank. The hierarchal Bayesian Analysis of Population Structure (hierBAPS) method was used to generate groups of sequences that were then projected onto a rooted maximum likelihood (ML) phylogenetic tree to visualize the pattern of clustering. The haplogroup statuses of the individual sequences were assessed using Haplogrep2.

Results

A total of 23 hierBAPS groups were identified, all of which corresponded to subclades defined in Phylotree, v.17. The hierBAPS groups projected onto the ML phylogeny accurately clustered all haplotypes belonging to a specific haplogroup in accordance with Haplogrep2. By incorporating the geographic source of each sequence and subclade age estimates into this framework, inferences about the diversification of U5 mtDNAs were made. Haplogroup U5 has been present in northern Europe since the Mesolithic, and spread in both eastern and western directions, undergoing significant diversification within Scandinavia. A review of historical and archeological evidence attests to some of the population interactions contributing to this pattern.

Conclusions

The hierBAPS algorithm accurately grouped mitogenome sequences into subclades in a phylogenetically robust manner. This analysis provided new insights into the phylogeographic structure of haplogroup U5 diversity in northern Europe, revealing a detailed perspective on the diversity of subclades in this region and their distribution in Scandinavian populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08572-y.

Genomic and dietary discontinuities during the Mesolithic and Neolithic in Sicily

Sicily is a key region for understanding the agricultural transition in the Mediterranean because of its central position. Here, we present genomic and stable isotopic data for 19 prehistoric Sicilians covering the Mesolithic to Bronze Age periods (10,700–4,100 yBP). We find that Early Mesolithic hunter-gatherers (HGs) from Sicily are a highly drifted lineage of the Early Holocene western European HGs, whereas Late Mesolithic HGs carry ∼20% ancestry related to northern and (south) eastern European HGs, indicating substantial gene flow. Early Neolithic farmers are genetically most similar to farmers from the Balkans and Greece, with only ∼7% of ancestry from local Mesolithic HGs. The genetic discontinuities during the Mesolithic and Early Neolithic match the changes in material culture and diet. Three outlying individuals dated to ∼8,000 yBP; however, suggest that hunter-gatherers interacted with incoming farmers at Grotta dell’Uzzo, resulting in a mixed economy and diet for a brief interlude at the Mesolithic-Neolithic transition.

•

Genetic transition between Early Mesolithic and Late Mesolithic hunter-gatherers

•

A near-complete genetic turnover during the Mesolithic-Neolithic transition

•

Exchange of subsistence practices between hunter-gatherers and early farmers

Mitochondrial DNA Consensus Calling and Quality Filtering for Constructing Ancient Human Mitogenomes: Comparison of Two Widely Applied Methods

Retrieving high-quality endogenous ancient DNA (aDNA) poses several challenges, including low molecular copy number, high rates of fragmentation, damage at read termini, and potential presence of exogenous contaminant DNA. All these factors complicate a reliable reconstruction of consensus aDNA sequences in reads from high-throughput sequencing platforms. Here, we report findings from a thorough evaluation of two alternative tools (ANGSD and schmutzi) aimed at overcoming these issues and constructing high-quality ancient mitogenomes. Raw genomic data (BAM/FASTQ) from a total of 17 previously published whole ancient human genomes ranging from the 14th to the 7th millennium BCE were retrieved and mitochondrial consensus sequences were reconstructed using different quality filters, with their accuracy measured and compared. Moreover, the influence of different sequence parameters (number of reads, sequenced bases, mean coverage, and rate of deamination and contamination) as predictors of derived sequence quality was evaluated. Complete mitogenomes were successfully reconstructed for all ancient samples, and for the majority of them, filtering substantially improved mtDNA consensus calling and haplogroup prediction. Overall, the schmutzi pipeline, which estimates and takes into consideration exogenous contamination, appeared to have the edge over the much faster and user-friendly alternative method (ANGSD) in moderate to high coverage samples (>1,000,000 reads). ANGSD, however, through its read termini trimming filter, showed better capabilities in calling the consensus sequence from low-quality samples. Among all the predictors of overall sample quality examined, the strongest correlation was found for the available number of sequence reads and bases. In the process, we report a previously unassigned haplogroup (U3b) for an Early Chalcolithic individual from Southern Anatolia/Northern Levant.

Simulated patterns of mitochondrial diversity are consistent with partial population turnover in Bronze Age Central Europe

OBJECTIVES

The analysis of ancient mitochondrial DNA from osteological remains has challenged previous conclusions drawn from the analysis of mitochondrial DNA from present populations, notably by revealing an absence of genetic continuity between the Neolithic and modern populations in Central Europe. Our study investigates how to reconcile these contradictions at the mitochondrial level using a modeling approach.

MATERIALS AND METHODS

We used a spatially explicit computational framework to simulate ancient and modern DNA sequences under various evolutionary scenarios of post Neolithic demographic events and compared the genetic diversity of the simulated and observed mitochondrial sequences. We investigated which-if any-scenarios were able to reproduce statistics of genetic diversity similar to those observed, with a focus on the haplogroup N1a, associated with the spread of early Neolithic farmers.

RESULTS

Demographic fluctuations during the Neolithic transition or subsequent demographic collapses after this period, that is, due to epidemics such as plague, are not sufficient to explain the signal of population discontinuity detected on the mitochondrial DNA in Central Europe. Only a scenario involving a substantial genetic input due to the arrival of migrants after the Neolithic transition, possibly during the Bronze Age, is compatible with observed patterns of genetic diversity.

DISCUSSION

Our results corroborate paleogenomic studies, since out of the alternative hypotheses tested, the best one that was able to recover observed patterns of mitochondrial diversity in modern and ancient Central European populations was one were immigration of populations from the Pontic steppes during the Bronze Age was explicitly simulated.

Bioarchaeological insights into the last plague of Imola (1630-1632)

The plague of 1630-1632 was one of the deadliest plague epidemics to ever hit Northern Italy, and for many of the affected regions, it was also the last. While accounts on plague during the early 1630s in Florence and Milan are frequent, much less is known about the city of Imola. We analyzed the full skeletal assemblage of four mass graves (n = 133 individuals) at the Lazaretto dell'Osservanza, which date back to the outbreak of 1630-1632 in Imola and evaluated our results by integrating new archival sources. The skeletons showed little evidence of physical trauma and were covered by multiple layers of lime, which is characteristic for epidemic mass mortality sites. We screened 15 teeth for Yersinia pestis aDNA and were able to confirm the presence of plague in Imola via metagenomic analysis. Additionally, we studied a contemporaneous register, in which a friar recorded patient outcomes at the lazaretto during the last year of the epidemic. Our multidisciplinary approach combining historical, osteological and genomic data provided a unique opportunity to reconstruct an in-depth picture of the last plague of Imola through the city's main lazaretto.

Observed Ecological Communities Are Formed by Species Combinations That Are among the Most Likely to Persist under Changing Environments

AbstractDespite the rich biodiversity found in nature, it is unclear to what extent some combinations of interacting species, while conceivable in a given place and time, may never be realized. Yet solving this problem is important for understanding the role of randomness and predictability in the assembly of ecological communities. Here we show that the specific combinations of interacting species that emerge from the ecological dynamics within regional species pools are not all equally likely to be seen; rather, they are among the most likely to persist under changing environments. First, we use niche-based competition matrices and Lotka-Volterra models to demonstrate that realized combinations of interacting species are more likely to persist under random parameter perturbations than the majority of potential combinations with the same number of species that could have been formed from the regional pool. We then corroborate our theoretical results using a 10-year observational study, recording 88 plant-herbivore communities across three different forest successional stages. By inferring and validating plant-mediated communities of competing herbivore species, we find that observed combinations of herbivores have an expected probability of species persistence higher than half of all potential combinations. Our findings open up the opportunity to establish a formal probabilistic and predictive understanding of the composition of ecological communities.

Genetic ancestry changes in Stone to Bronze Age transition in the East European plain

The transition from Stone to Bronze Age in Central and Western Europe was a period of major population movements originating from the Ponto-Caspian Steppe. Here, we report new genome-wide sequence data from 30 individuals north of this area, from the understudied western part of present-day Russia, including 3 Stone Age hunter-gatherers (10,800 to 4250 cal BCE) and 26 Bronze Age farmers from the Corded Ware complex Fatyanovo Culture (2900 to 2050 cal BCE). We show that Eastern hunter-gatherer ancestry was present in northwestern Russia already from around 10,000 BCE. Furthermore, we see a change in ancestry with the arrival of farming-Fatyanovo Culture individuals were genetically similar to other Corded Ware cultures, carrying a mixture of Steppe and European early farmer ancestry. Thus, they likely originate from a fast migration toward the northeast from somewhere near modern-day Ukraine-the closest area where these ancestries coexisted from around 3000 BCE.

A genomic and historical synthesis of plague in 18th century Eurasia

Plague continued to afflict Europe for more than five centuries after the Black Death. Yet, by the 17th century, the dynamics of plague had changed, leading to its slow decline in Western Europe over the subsequent 200 y, a period for which only one genome was previously available. Using a multidisciplinary approach, combining genomic and historical data, we assembled Y. pestis genomes from nine individuals covering four Eurasian sites and placed them into an historical context within the established phylogeny. CHE1 (Chechnya, Russia, 18th century) is now the latest Second Plague Pandemic genome and the first non-European sample in the post-Black Death lineage. Its placement in the phylogeny and our synthesis point toward the existence of an extra-European reservoir feeding plague into Western Europe in multiple waves. By considering socioeconomic, ecological, and climatic factors we highlight the importance of a noneurocentric approach for the discussion on Second Plague Pandemic dynamics in Europe.

2000-year-old pathogen genomes reconstructed from metagenomic analysis of Egyptian mummified individuals

BACKGROUND

Recent advances in sequencing have facilitated large-scale analyses of the metagenomic composition of different samples, including the environmental microbiome of air, water, and soil, as well as the microbiome of living humans and other animals. Analyses of the microbiome of ancient human samples may provide insights into human health and disease, as well as pathogen evolution, but the field is still in its very early stages and considered highly challenging.

RESULTS

The metagenomic and pathogen content of Egyptian mummified individuals from different time periods was investigated via genetic analysis of the microbial composition of various tissues. The analysis of the dental calculus' microbiome identified Red Complex bacteria, which are correlated with periodontal diseases. From bone and soft tissue, genomes of two ancient pathogens, a 2200-year-old Mycobacterium leprae strain and a 2000-year-old human hepatitis B virus, were successfully reconstructed.

CONCLUSIONS

The results show the reliability of metagenomic studies on Egyptian mummified individuals and the potential to use them as a source for the extraction of ancient pathogen DNA.

Species multidimensional effects explain idiosyncratic responses of communities to environmental change

Environmental change can alter species' abundances within communities consistently; for example, increasing all abundances by the same percentage, or more idiosyncratically. Here, we show how comparing effects of temperature on species grown in isolation and when grown together helps our understanding of how ecological communities more generally respond to environmental change. In particular, we find that the shape of the feasibility domain (the parameter space of carrying capacities compatible with positive species' abundances) helps to explain the composition of experimental microbial communities under changing environmental conditions. First, we introduce a measure to quantify the asymmetry of a community's feasibility domain using the column vectors of the corresponding interaction matrix. These column vectors describe the effects each species has on all other species in the community (hereafter referred to as species' multidimensional effects). We show that as the asymmetry of the feasibility domain increases the relationship between species' abundance when grown together and when grown in isolation weakens. We then show that microbial communities experiencing different temperature environments exhibit patterns consistent with this theory. Specifically, communities at warmer temperatures show relatively more asymmetry; thus, the idiosyncrasy of responses is higher compared with that in communities at cooler temperatures. These results suggest that while species' interactions are typically defined at the pairwise level, multispecies dynamics can be better understood by focusing on the effects of these interactions at the community level.

Ancient genome-wide DNA from France highlights the complexity of interactions between Mesolithic hunter-gatherers and Neolithic farmers

Starting from 12,000 years ago in the Middle East, the Neolithic lifestyle spread across Europe via separate continental and Mediterranean routes. Genomes from early European farmers have shown a clear Near Eastern/Anatolian genetic affinity with limited contribution from hunter-gatherers. However, no genomic data are available from modern-day France, where both routes converged, as evidenced by a mosaic cultural pattern. Here, we present genome-wide data from 101 individuals from 12 sites covering today's France and Germany from the Mesolithic (N = 3) to the Neolithic (N = 98) (7000-3000 BCE). Using the genetic substructure observed in European hunter-gatherers, we characterize diverse patterns of admixture in different regions, consistent with both routes of expansion. Early western European farmers show a higher proportion of distinctly western hunter-gatherer ancestry compared to central/southeastern farmers. Our data highlight the complexity of the biological interactions during the Neolithic expansion by revealing major regional variations.

Interaction Between Root Exudates of the Poisonous Plant Stellera chamaejasme L. and Arbuscular Mycorrhizal Fungi on the Growth of Leymus chinensis (Trin.) Tzvel

The growth of a large number of poisonous plants is an indicator of grassland degradation. Releasing allelochemicals through root exudates is one of the strategies with which poisonous plants affect neighboring plants in nature. Arbuscular mycorrhizal fungi (AMF) can form a mutualistic symbiosis with most of the higher plants. However, the manner of interaction between root exudates of poisonous plants and AMF on neighboring herbage in grasslands remains poorly understood. Stellera chamaejasme L., a common poisonous plant with approved allelopathy, is widely distributed with the dominant grass of Leymus chinensis in the degradeds of Northern China. In this study, we investigated the addition of S. chamaejasme root exudates (SRE), the inoculation of AMF, and their interaction on the growth and tissue nitrogen contents of L. chinensis, the characteristics of rhizosphere AMF, and soil physicochemical properties. Results showed that SRE had significant effects on ramet number, aboveground biomass, and total nitrogen of L. chinensis in a concentration dependent manner. Additionally, SRE had a significant negative effect on the rate of mycorrhiza infection and spore density of the AMF. Meanwhile, the addition of SRE significantly affected soil pH, electrical conductivity, available nitrogen (AN), available phosphorus (AP), total nitrogen (TN), and total carbon (TC) contents; while neither inoculation of AMF itself nor the interaction of AMF with SRE significantly affected the growth of L. chinensis. The interaction between AMF and SRE dramatically changed the pH, AP, and TC of rhizosphere soil. Therefore, we suggested SRE of S. chamaejasme affected the growth of L. chinensis by altering soil pH and nutrient availability. AMF could change the effect of SRE on soil nutrients and have the potential to regulate the allelopathic effects of S. chamaejasme and the interspecific interaction between the two plant species. We have provided new evidence for the allelopathic mechanism of S. chamaejasme and the regulation effects of AMF on the interspecific relationship between poisonous plants and neighboring plants. Our findings reveal the complex interplay between the root exudates of poisonous plants and rhizosphere AMF in regulating population growth and dynamics of neighboring plants in degraded grassland ecosystems.

The genetic landscape of Arab Population, Chechens and Circassians subpopulations from Jordan through HV1 and HV2 regions of mtDNA

Mitochondrial DNA (mtDNA) is widely used in several fields including medical genetics, forensic science, genetic genealogy, and evolutionary anthropology. In this study, mtDNA haplotype diversity was determined for 293 unrelated subjects from Jordanian population (Circassians, Chechens, and the original inhabitants of Jordan). A total of 102 haplotypes were identified and analyzed among the populations to describe the maternal lineage landscape. Our results revealed that the distribution of mtDNA haplotype frequencies among the three populations showed disparity and significant differences when compared to each other. We also constructed mitochondrial haplotype classification trees for the three populations to determine the phylogenetic relationship of mtDNA haplotype variants, and we observed clear differences in the distribution of maternal genetic ancestries, especially between Arab and the minority ethnic populations. To our knowledge, this study is the first, to date, to characterize mitochondrial haplotypes and haplotype distributions in a population-based sample from the Jordanian population. It provides a powerful reference for future studies investigating the contribution of mtDNA variation to human health and disease and studying population history and evolution by comparing the mtDNA haplotypes to other populations.

Human mitochondrial DNA lineages in Iron-Age Fennoscandia suggest incipient admixture and eastern introduction of farming-related maternal ancestry

Human ancient DNA studies have revealed high mobility in Europe’s past, and have helped to decode the human history on the Eurasian continent. Northeastern Europe, especially north of the Baltic Sea, however, remains less well understood largely due to the lack of preserved human remains. Finland, with a divergent population history from most of Europe, offers a unique perspective to hunter-gatherer way of life, but thus far genetic information on prehistoric human groups in Finland is nearly absent. Here we report 103 complete ancient mitochondrial genomes from human remains dated to AD 300–1800, and explore mtDNA diversity associated with hunter-gatherers and Neolithic farmers. The results indicate largely unadmixed mtDNA pools of differing ancestries from Iron-Age on, suggesting a rather late genetic shift from hunter-gatherers towards farmers in North-East Europe. Furthermore, the data suggest eastern introduction of farmer-related haplogroups into Finland, contradicting contemporary genetic patterns in Finns.

On the Consequences of the Interdependence of Stabilizing and Equalizing Mechanisms

We present an overlooked but important property of modern coexistence theory (MCT), along with two key new results and their consequences. The overlooked property is that stabilizing mechanisms (increasing species' niche differences) and equalizing mechanisms (reducing species' fitness differences) have two distinct sets of meanings within MCT: one in a two-species context and another in a general multispecies context. We demonstrate that the two-species framework is not a special case of the multispecies one, and therefore these two parallel frameworks must be studied independently. Our first result is that, using the two-species framework and mechanistic consumer-resource models, stabilizing and equalizing mechanisms exhibit complex interdependence, such that changing one will simultaneously change the other. Furthermore, the nature and direction of this simultaneous change sensitively depend on model parameters. The second result states that while MCT is often seen as bridging niche and neutral modes of coexistence by building a niche-neutrality continuum, the interdependence between stabilizing and equalizing mechanisms acts to break this continuum under almost any biologically relevant circumstance. We conclude that the complex entanglement of stabilizing and equalizing terms makes their impact on coexistence difficult to understand, but by seeing them as aggregated effects (rather than underlying causes) of coexistence, we may increase our understanding of ecological dynamics.

Disentangling ecologically equivalent from neutral species: The mechanisms of population regulation matter

The neutral theory of biodiversity explored the structure of a community of ecologically equivalent species. Such species are expected to display community drift dynamics analogous to neutral alleles undergoing genetic drift. While entire communities of species are not ecologically equivalent, recent field experiments have documented the existence of guilds of such neutral species embedded in real food webs. What demographic outcomes of the interactions within and between species in these guilds are expected to produce ecological drift versus coexistence remains unclear. To address this issue, and guide empirical testing, we consider models of a guild of ecologically equivalent competitors feeding on a single resource to explore when community drift should manifest. We show that community drift dynamics only emerge when the density-dependent effects of each species on itself are identical to its density-dependent effects on every other guild member. In contrast, if each guild member directly limits itself more than it limits the abundance of other guild members, all species in the guild are coexisting, even though they all are ecologically equivalent with respect to their interactions with species outside the guild (i.e. resources, predators, mutualists). Hence, considering only interspecific ecological differences generating density dependence, and not fully accounting for the preponderance of mechanisms causing intraspecific density dependence, will provide an incomplete picture for segregating between neutrality and coexistence. We also identify critical experiments necessary to disentangle guilds of ecologically equivalent species from those experiencing ecological drift, as well as provide an overview of ways of incorporating a mechanistic basis into studies of species coexistence and neutrality. Identifying these characteristics, and the mechanistic basis underlying community structure, is not merely an exercise in clarifying the semantics of coexistence and neutral theories, but rather reflects key differences that must exist among community members in order to determine how and why communities are structured.

Broadly inflicted stressors can cause ecosystem thinning

Many anthropogenic stressors broadly inflict mortality or reduce fecundity, including habitat destruction, pollution, climate change, invasive species, and multispecies harvesting. Here, we show-in four analytical models of interspecies competition-that broadly inflicted stressors disproportionately cause competitive exclusions within groups of ecologically similar species. As a result, we predict that ecosystems become progressively thinner-that is, they have progressively less functional redundancy-as broadly inflicted stressors become progressively more intense. This may negatively affect the temporal stability of ecosystem functions, but it also buffers ecosystem productivity against stress by favoring species less sensitive to the stressors. Our main result follows from the weak limiting similarity principle: species with more similar ecological niches compete more strongly, and their coexistence can be upset by smaller perturbations. We show that stressors can cause indirect competitive exclusions at much lower stressor intensity than needed to directly cause species extinction, consistent with the finding of empirical studies that species interactions are often the proximal drivers of local extinctions. The excluded species are more sensitive to the stressor relative to their ecologically similar competitors. Moreover, broadly inflicted stressors may cause hydra effects-where higher stressor intensity results in higher abundance for a species with lower sensitivity to the stressor than its competitors. Correlations between stressor impacts and ecological niches reduce the potential for indirect competitive exclusions, but they consequently also reduce the buffering effect of ecosystem thinning on ecosystem productivity. Our findings suggest that ecosystems experiencing stress may continue to provision ecosystem services but lose functional redundancy and stability.

Integrative approach using Yersinia pestis genomes to revisit the historical landscape of plague during the Medieval Period

Over the last few years, genomic studies on Yersinia pestis, the causative agent of all known plague epidemics, have considerably increased in numbers, spanning a period of about 5,000 y. Nonetheless, questions concerning historical reservoirs and routes of transmission remain open. Here, we present and describe five genomes from the second half of the 14th century and reconstruct the evolutionary history of Y. pestis by reanalyzing previously published genomes and by building a comprehensive phylogeny focused on strains attributed to the Second Plague Pandemic (14th to 18th century). Corroborated by historical and ecological evidence, the presented phylogeny, which includes our Y. pestis genomes, could support the hypothesis of an entry of plague into Western European ports through distinct waves of introduction during the Medieval Period, possibly by means of fur trade routes, as well as the recirculation of plague within the human population via trade routes and human movement.

Emergence and Spread of Basal Lineages of Yersinia pestis during the Neolithic Decline

Between 5,000 and 6,000 years ago, many Neolithic societies declined throughout western Eurasia due to a combination of factors that are still largely debated. Here, we report the discovery and genome reconstruction of Yersinia pestis, the etiological agent of plague, in Neolithic farmers in Sweden, pre-dating and basal to all modern and ancient known strains of this pathogen. We investigated the history of this strain by combining phylogenetic and molecular clock analyses of the bacterial genome, detailed archaeological information, and genomic analyses from infected individuals and hundreds of ancient human samples across Eurasia. These analyses revealed that multiple and independent lineages of Y. pestis branched and expanded across Eurasia during the Neolithic decline, spreading most likely through early trade networks rather than massive human migrations. Our results are consistent with the existence of a prehistoric plague pandemic that likely contributed to the decay of Neolithic populations in Europe.

Bayesian estimation of partial population continuity using ancient DNA and spatially explicit simulations

The retrieval of ancient DNA from osteological material provides direct evidence of human genetic diversity in the past. Ancient DNA samples are often used to investigate whether there was population continuity in the settlement history of an area. Methods based on the serial coalescent algorithm have been developed to test whether the population continuity hypothesis can be statistically rejected by analysing DNA samples from the same region but of different ages. Rejection of this hypothesis is indicative of a large genetic shift, possibly due to immigration occurring between two sampling times. However, this approach is only able to reject a model of full continuity model (a total absence of genetic input from outside), but admixture between local and immigrant populations may lead to partial continuity. We have recently developed a method to test for population continuity that explicitly considers the spatial and temporal dynamics of populations. Here, we extended this approach to estimate the proportion of genetic continuity between two populations, using ancient genetic samples. We applied our original approach to the question of the Neolithic transition in Central Europe. Our results confirmed the rejection of full continuity, but our approach represents an important step forward by estimating the relative contribution of immigrant farmers and of local hunter-gatherers to the final Central European Neolithic genetic pool. Furthermore, we show that a substantial proportion of genes brought by the farmers in this region were assimilated from other hunter-gatherer populations along the way from Anatolia, which was not detectable by previous continuity tests. Our approach is also able to jointly estimate demographic parameters, as we show here by finding both low density and low migration rate for pre-Neolithic hunter-gatherers. It provides a useful tool for the analysis of the numerous ancient DNA data sets that are currently being produced for many different species.

Maternal genetic features of the Iron Age Tagar population from Southern Siberia (1st millennium BC)

Early nomads in the Eurasian steppes since the beginning of the 1st millennium BC played a key role in the formation of the cultural and genetic landscape of populations of a significant part of Eurasia, from Eastern Europe to Eastern Central Asia. Numerous archaeological cultures associated with early nomads have been discovered throughout the Eurasian steppe belt. The Tagar archaeological culture existed in the Minusinsk basin (Sayan Mountains, Southern Siberia, Russia) in the northeastern periphery of the Eurasian steppe belt from the 8th to 1st century BC during the pre-Scythian, Scythian, and Early Xiongnu-Sarmatian periods. In this study, we evaluated mtDNA diversity in the Tagar population based on representative series (N = 79) belonging to all chronological stages of the culture. The Tagar population had a mixed mtDNA pool dominated by Western Eurasian haplogroups and subgroups (H, HV6, HV*, I, K, T, U2e, U4, U5a, and U*) and, to a lesser degree, Eastern Eurasian haplogroups (A*, A8, C*, C5, D, G2a, and F1b). The Tagar population showed a similar mtDNA pool structure to those of other Iron Age populations representing the "Scythian World." We observed particularly high similarity between the Tagar and Classic Scythians from the North Pontic region. Our results support the assumption that genetic components introduced by Bronze Age migrants from Western Eurasia contributed to the formation of the genetic composition of Scythian period populations in Southern Siberia. Another important component of the Tagar mtDNA pool was autochthonous East Eurasian lineages, some of which (A8 and C4a2a) are potential markers of the westward genetic influence of the eastern populations of the Scythian period. Our results suggest a genetic continuity (at least partial) between the Early, Middle, and Late Tagar populations.

Whole mitochondrial genome diversity in two Hungarian populations

Complete mitochondrial genomics is an effective tool for studying the demographic history of human populations, but there is still a deficit of mitogenomic data in European populations. In this paper, we present results of study of variability of 80 complete mitochondrial genomes in two Hungarian populations from eastern part of Hungary (Szeged and Debrecen areas). The genetic diversity of Hungarian mitogenomes is remarkably high, reaching 99.9% in a combined sample. According to the analysis of molecular variance (AMOVA), European populations showed a low, but statistically significant level of between-population differentiation (Fst = 0.61%, p = 0), and two Hungarian populations demonstrate lack of between-population differences. Phylogeographic analysis allowed us to identify 71 different mtDNA sub-clades in Hungarians, sixteen of which are novel. Analysis of ancestry-informative mtDNA sub-clades revealed a complex genetic structure associated with the genetic impact of populations from different parts of Eurasia, though the contribution from European populations is the most pronounced. At least 8% of ancestry-informative haplotypes found in Hungarians demonstrate similarity with East and West Slavic populations (sub-clades H1c23a, H2a1c1, J2b1a6, T2b25a1, U4a2e, K1c1j, and I1a1c), while the influence of Siberian populations is not so noticeable (sub-clades A12a, C4a1a, and probably U4b1a4).

Structural stability as a consistent predictor of phenological events

The timing of the first and last seasonal appearance of a species in a community typically follows a pattern that is governed by temporal factors. While it has been shown that changes in the environment are linked to phenological changes, the direction of this link appears elusive and context-dependent. Thus, finding consistent predictors of phenological events is of central importance for a better assessment of expected changes in the temporal dynamics of ecological communities. Here we introduce a measure of structural stability derived from species interaction networks as an estimator of the expected range of environmental conditions compatible with the existence of a community. We test this measure as a predictor of changes in species richness recorded on a daily basis in a high-arctic plant-pollinator community during two spring seasons. We find that our measure of structural stability is the only consistent predictor of changes in species richness among different ecological and environmental variables. Our findings suggest that measures based on the notion of structural stability can synthesize the expected variation of environmental conditions tolerated by a community, and explain more consistently the phenological changes observed in ecological communities.

Investigating Holocene human population history in North Asia using ancient mitogenomes

Archaeogenomic studies have largely elucidated human population history in West Eurasia during the Stone Age. However, despite being a broad geographical region of significant cultural and linguistic diversity, little is known about the population history in North Asia. We present complete mitochondrial genome sequences together with stable isotope data for 41 serially sampled ancient individuals from North Asia, dated between c.13,790 BP and c.1,380 BP extending from the Palaeolithic to the Iron Age. Analyses of mitochondrial DNA sequences and haplogroup data of these individuals revealed the highest genetic affinity to present-day North Asian populations of the same geographical region suggesting a possible long-term maternal genetic continuity in the region. We observed a decrease in genetic diversity over time and a reduction of maternal effective population size (N_e) approximately seven thousand years before present. Coalescent simulations were consistent with genetic continuity between present day individuals and individuals dating to 7,000 BP, 4,800 BP or 3,000 BP. Meanwhile, genetic differences observed between 7,000 BP and 3,000 BP as well as between 4,800 BP and 3,000 BP were inconsistent with genetic drift alone, suggesting gene flow into the region from distant gene pools or structure within the population. These results indicate that despite some level of continuity between ancient groups and present-day populations, the region exhibits a complex demographic history during the Holocene.

Dispersals of the Siberian Y-chromosome haplogroup Q in Eurasia

The human Y-chromosome has proven to be a powerful tool for tracing the paternal history of human populations and genealogical ancestors. The human Y-chromosome haplogroup Q is the most frequent haplogroup in the Americas. Previous studies have traced the origin of haplogroup Q to the region around Central Asia and Southern Siberia. Although the diversity of haplogroup Q in the Americas has been studied in detail, investigations on the diffusion of haplogroup Q in Eurasia and Africa are still limited. In this study, we collected 39 samples from China and Russia, investigated 432 samples from previous studies of haplogroup Q, and analyzed the single nucleotide polymorphism (SNP) subclades Q1a1a1-M120, Q1a2a1-L54, Q1a1b-M25, Q1a2-M346, Q1a2a1a2-L804, Q1a2b2-F1161, Q1b1a-M378, and Q1b1a1-L245. Through NETWORK and BATWING analyses, we found that the subclades of haplogroup Q continued to disperse from Central Asia and Southern Siberia during the past 10,000 years. Apart from its migration through the Beringia to the Americas, haplogroup Q also moved from Asia to the south and to the west during the Neolithic period, and subsequently to the whole of Eurasia and part of Africa.

Reconciling evidence from ancient and contemporary genomes: a major source for the European Neolithic within Mediterranean Europe

Important gaps remain in our understanding of the spread of farming into Europe, due partly to apparent contradictions between studies of contemporary genetic variation and ancient DNA. It seems clear that farming was introduced into central, northern, and eastern Europe from the south by pioneer colonization. It is often argued that these dispersals originated in the Near East, where the potential source genetic pool resembles that of the early European farmers, but clear ancient DNA evidence from Mediterranean Europe is lacking, and there are suggestions that Mediterranean Europe may have resembled the Near East more than the rest of Europe in the Mesolithic. Here, we test this proposal by dating mitogenome founder lineages from the Near East in different regions of Europe. We find that whereas the lineages date mainly to the Neolithic in central Europe and Iberia, they largely date to the Late Glacial period in central/eastern Mediterranean Europe. This supports a scenario in which the genetic pool of Mediterranean Europe was partly a result of Late Glacial expansions from a Near Eastern refuge, and that this formed an important source pool for subsequent Neolithic expansions into the rest of Europe.

Investigating population continuity with ancient DNA under a spatially explicit simulation framework

BACKGROUND

Recent advances in sequencing technologies have allowed for the retrieval of ancient DNA data (aDNA) from skeletal remains, providing direct genetic snapshots from diverse periods of human prehistory. Comparing samples taken in the same region but at different times, hereafter called "serial samples", may indicate whether there is continuity in the peopling history of that area or whether an immigration of a genetically different population has occurred between the two sampling times. However, the exploration of genetic relationships between serial samples generally ignores their geographical locations and the spatiotemporal dynamics of populations. Here, we present a new coalescent-based, spatially explicit modelling approach to investigate population continuity using aDNA, which includes two fundamental elements neglected in previous methods: population structure and migration. The approach also considers the extensive temporal and geographical variance that is commonly found in aDNA population samples.

RESULTS

We first showed that our spatially explicit approach is more conservative than the previous (panmictic) approach and should be preferred to test for population continuity, especially when small and isolated populations are considered. We then applied our method to two mitochondrial datasets from Germany and France, both including modern and ancient lineages dating from the early Neolithic. The results clearly reject population continuity for the maternal line over the last 7500 years for the German dataset but not for the French dataset, suggesting regional heterogeneity in post-Neolithic migratory processes.

CONCLUSIONS

Here, we demonstrate the benefits of using a spatially explicit method when investigating population continuity with aDNA. It constitutes an improvement over panmictic methods by considering the spatiotemporal dynamics of genetic lineages and the precise location of ancient samples. The method can be used to investigate population continuity between any pair of serial samples (ancient-ancient or ancient-modern) and to investigate more complex evolutionary scenarios. Although we based our study on mitochondrial DNA sequences, diploid molecular markers of different types (DNA, SNP, STR) can also be simulated with our approach. It thus constitutes a promising tool for the analysis of the numerous aDNA datasets being produced, including genome wide data, in humans but also in many other species.

The maternal genetic make-up of the Iberian Peninsula between the Neolithic and the Early Bronze Age

Agriculture first reached the Iberian Peninsula around 5700 BCE. However, little is known about the genetic structure and changes of prehistoric populations in different geographic areas of Iberia. In our study, we focus on the maternal genetic makeup of the Neolithic (~ 5500–3000 BCE), Chalcolithic (~ 3000–2200 BCE) and Early Bronze Age (~ 2200–1500 BCE). We report ancient mitochondrial DNA results of 213 individuals (151 HVS-I sequences) from the northeast, central, southeast and southwest regions and thus on the largest archaeogenetic dataset from the Peninsula to date. Similar to other parts of Europe, we observe a discontinuity between hunter-gatherers and the first farmers of the Neolithic. During the subsequent periods, we detect regional continuity of Early Neolithic lineages across Iberia, however the genetic contribution of hunter-gatherers is generally higher than in other parts of Europe and varies regionally. In contrast to ancient DNA findings from Central Europe, we do not observe a major turnover in the mtDNA record of the Iberian Late Chalcolithic and Early Bronze Age, suggesting that the population history of the Iberian Peninsula is distinct in character.

Beyond pairwise mechanisms of species coexistence in complex communities

The tremendous diversity of species in ecological communities has motivated a century of research into the mechanisms that maintain biodiversity. However, much of this work examines the coexistence of just pairs of competitors. This approach ignores those mechanisms of coexistence that emerge only in diverse competitive networks. Despite the potential for these mechanisms to create conditions under which the loss of one competitor triggers the loss of others, we lack the knowledge needed to judge their importance for coexistence in nature. Progress requires borrowing insight from the study of multitrophic interaction networks, and coupling empirical data to models of competition.

Extensive Farming in Estonia Started through a Sex-Biased Migration from the Steppe

The transition from hunting and gathering to farming in Europe was brought upon by arrival of new people carrying novel material culture and genetic ancestry. The exact nature and scale of the transition-both material and genetic-varied in different parts of Europe [1-7]. Farming-based economies appear relatively late in Northeast Europe, and the extent to which they involve change in genetic ancestry is not fully understood due to the lack of relevant ancient DNA data. Here we present the results from new low-coverage whole-genome shotgun sequence data from five hunter-gatherers and five first farmers of Estonia whose remains date to 4,500 to 6,300 years before present. We find evidence of significant differences between the two groups in the composition of autosomal as well as mtDNA, X chromosome, and Y chromosome ancestries. We find that Estonian hunter-gatherers of Comb Ceramic culture are closest to Eastern hunter-gatherers, which is in contrast to earlier hunter-gatherers from the Baltics, who are close to Western hunter-gatherers [8, 9]. The Estonian first farmers of Corded Ware culture show high similarity in their autosomes with European hunter-gatherers, Steppe Eneolithic and Bronze Age populations, and European Late Neolithic/Bronze Age populations, while their X chromosomes are in addition equally closely related to European and Anatolian and Levantine early farmers. These findings suggest that the shift to intensive cultivation and animal husbandry in Estonia was triggered by the arrival of new people with predominantly Steppe ancestry but whose ancestors had undergone sex-specific admixture with early farmers with Anatolian ancestry.

Carriers of mitochondrial DNA macrohaplogroup R colonized Eurasia and Australasia from a southeast Asia core area

BACKGROUND

The colonization of Eurasia and Australasia by African modern humans has been explained, nearly unanimously, as the result of a quick southern coastal dispersal route through the Arabian Peninsula, the Indian subcontinent, and the Indochinese Peninsula, to reach Australia around 50 kya. The phylogeny and phylogeography of the major mitochondrial DNA Eurasian haplogroups M and N have played the main role in giving molecular genetics support to that scenario. However, using the same molecular tools, a northern route across central Asia has been invoked as an alternative that is more conciliatory with the fossil record of East Asia. Here, we assess as the Eurasian macrohaplogroup R fits in the northern path.

RESULTS

Haplogroup U, with a founder age around 50 kya, is one of the oldest clades of macrohaplogroup R in western Asia. The main branches of U expanded in successive waves across West, Central and South Asia before the Last Glacial Maximum. All these dispersions had rather overlapping ranges. Some of them, as those of U6 and U3, reached North Africa. At the other end of Asia, in Wallacea, another branch of macrohaplogroup R, haplogroup P, also independently expanded in the area around 52 kya, in this case as isolated bursts geographically well structured, with autochthonous branches in Australia, New Guinea, and the Philippines.

CONCLUSIONS

Coeval independently dispersals around 50 kya of the West Asia haplogroup U and the Wallacea haplogroup P, points to a halfway core area in southeast Asia as the most probable centre of expansion of macrohaplogroup R, what fits in the phylogeographic pattern of its ancestor, macrohaplogroup N, for which a northern route and a southeast Asian origin has been already proposed.

Origin and spread of human mitochondrial DNA haplogroup U7

Human mitochondrial DNA haplogroup U is among the initial maternal founders in Southwest Asia and Europe and one that best indicates matrilineal genetic continuity between late Pleistocene hunter-gatherer groups and present-day populations of Europe. While most haplogroup U subclades are older than 30 thousand years, the comparatively recent coalescence time of the extant variation of haplogroup U7 (~16–19 thousand years ago) suggests that its current distribution is the consequence of more recent dispersal events, despite its wide geographical range across Europe, the Near East and South Asia. Here we report 267 new U7 mitogenomes that – analysed alongside 100 published ones – enable us to discern at least two distinct temporal phases of dispersal, both of which most likely emanated from the Near East. The earlier one began prior to the Holocene (~11.5 thousand years ago) towards South Asia, while the later dispersal took place more recently towards Mediterranean Europe during the Neolithic (~8 thousand years ago). These findings imply that the carriers of haplogroup U7 spread to South Asia and Europe before the suggested Bronze Age expansion of Indo-European languages from the Pontic-Caspian Steppe region.

A genetic chronology for the Indian Subcontinent points to heavily sex-biased dispersals

BACKGROUND

India is a patchwork of tribal and non-tribal populations that speak many different languages from various language families. Indo-European, spoken across northern and central India, and also in Pakistan and Bangladesh, has been frequently connected to the so-called "Indo-Aryan invasions" from Central Asia ~3.5 ka and the establishment of the caste system, but the extent of immigration at this time remains extremely controversial. South India, on the other hand, is dominated by Dravidian languages. India displays a high level of endogamy due to its strict social boundaries, and high genetic drift as a result of long-term isolation which, together with a very complex history, makes the genetic study of Indian populations challenging.

RESULTS

We have combined a detailed, high-resolution mitogenome analysis with summaries of autosomal data and Y-chromosome lineages to establish a settlement chronology for the Indian Subcontinent. Maternal lineages document the earliest settlement ~55-65 ka (thousand years ago), and major population shifts in the later Pleistocene that explain previous dating discrepancies and neutrality violation. Whilst current genome-wide analyses conflate all dispersals from Southwest and Central Asia, we were able to tease out from the mitogenome data distinct dispersal episodes dating from between the Last Glacial Maximum to the Bronze Age. Moreover, we found an extremely marked sex bias by comparing the different genetic systems.

CONCLUSIONS

Maternal lineages primarily reflect earlier, pre-Holocene processes, and paternal lineages predominantly episodes within the last 10 ka. In particular, genetic influx from Central Asia in the Bronze Age was strongly male-driven, consistent with the patriarchal, patrilocal and patrilineal social structure attributed to the inferred pastoralist early Indo-European society. This was part of a much wider process of Indo-European expansion, with an ultimate source in the Pontic-Caspian region, which carried closely related Y-chromosome lineages, a smaller fraction of autosomal genome-wide variation and an even smaller fraction of mitogenomes across a vast swathe of Eurasia between 5 and 3.5 ka.