Andrea Benazzo

ABC as a flexible framework to estimate demography over space and time: some cons, many pros

The analysis of genetic variation to estimate demographic and historical parameters and to quantitatively compare alternative scenarios recently gained a powerful and flexible approach: the Approximate Bayesian Computation (ABC). The likelihood functions does not need to be theoretically specified, but posterior distributions can be approximated by simulation even assuming very complex population models including both natural and human‐induced processes. Prior information can be easily incorporated and the quality of the results can be analysed with rather limited additional effort. ABC is not a statistical analysis per se, but rather a statistical framework and any specific application is a sort of hybrid between a simulation and a data‐analysis study. Complete software packages performing the necessary steps under a set of models and for specific genetic markers are already available, but the flexibility of the method is better exploited combining different programs. Many questions relevant in ecology can be addressed using ABC, but adequate amount of time should be dedicated to decide among alternative options and to evaluate the results. In this paper we will describe and critically comment on the different steps of an ABC analysis, analyse some of the published applications of ABC and provide user guidelines.

Impact of type 2 diabetes susceptibility variants on quantitative glycemic traits reveals mechanistic heterogeneity

Patients with established type 2 diabetes display both β-cell dysfunction and insulin resistance. To define fundamental processes leading to the diabetic state, we examined the relationship between type 2 diabetes risk variants at 37 established susceptibility loci, and indices of proinsulin processing, insulin secretion, and insulin sensitivity. We included data from up to 58,614 nondiabetic subjects with basal measures and 17,327 with dynamic measures. We used additive genetic models with adjustment for sex, age, and BMI, followed by fixed-effects, inverse-variance meta-analyses. Cluster analyses grouped risk loci into five major categories based on their relationship to these continuous glycemic phenotypes. The first cluster (PPARG, KLF14, IRS1, GCKR) was characterized by primary effects on insulin sensitivity. The second cluster (MTNR1B, GCK) featured risk alleles associated with reduced insulin secretion and fasting hyperglycemia. ARAP1 constituted a third cluster characterized by defects in insulin processing. A fourth cluster (TCF7L2, SLC30A8, HHEX/IDE, CDKAL1, CDKN2A/2B) was defined by loci influencing insulin processing and secretion without a detectable change in fasting glucose levels. The final group contained 20 risk loci with no clear-cut associations to continuous glycemic traits. By assembling extensive data on continuous glycemic traits, we have exposed the diverse mechanisms whereby type 2 diabetes risk variants impact disease predisposition.

Inferring Genealogical Processes from Patterns of Bronze-Age and Modern DNA Variation in Sardinia

The ancient inhabitants of a region are often regarded as ancestral, and hence genetically related, to the modern dwellers (for instance, in studies of admixture), but so far, this assumption has not been tested empirically using ancient DNA data. We studied mitochondrial DNA (mtDNA) variation in Sardinia, across a time span of 2,500 years, comparing 23 Bronze-Age (nuragic) mtDNA sequences with those of 254 modern individuals from two regions, Ogliastra (a likely genetic isolate) and Gallura, and considering the possible impact of gene flow from mainland Italy. To understand the genealogical relationships between past and present populations, we developed seven explicit demographic models; we tested whether these models can account for the levels and patterns of genetic diversity in the data and which one does it best. Extensive simulation based on a serial coalescent algorithm allowed us to compare the posterior probability of each model and estimate the relevant evolutionary (mutation and migration rates) and demographic (effective population sizes, times since population splits) parameters, by approximate Bayesian computations. We then validated the analyses by investigating how well parameters estimated from the simulated data can reproduce the observed data set. We show that a direct genealogical continuity between Bronze-Age Sardinians and the current people of Ogliastra, but not Gallura, has a much higher probability than any alternative scenarios and that genetic diversity in Gallura evolved largely independently, owing in part to gene flow from the mainland.

Large-scale recent expansion of European patrilineages shown by population resequencing

The proportion of Europeans descending from Neolithic farmers ∼10 thousand years ago (KYA) or Palaeolithic hunter-gatherers has been much debated. The male-specific region of the Y chromosome (MSY) has been widely applied to this question, but unbiased estimates of diversity and time depth have been lacking. Here we show that European patrilineages underwent a recent continent-wide expansion. Resequencing of 3.7 Mb of MSY DNA in 334 males, comprising 17 European and Middle Eastern populations, defines a phylogeny containing 5,996 single-nucleotide polymorphisms. Dating indicates that three major lineages (I1, R1a and R1b), accounting for 64% of our sample, have very recent coalescent times, ranging between 3.5 and 7.3 KYA. A continuous swathe of 13/17 populations share similar histories featuring a demographic expansion starting ∼2.1–4.2 KYA. Our results are compatible with ancient MSY DNA data, and contrast with data on mitochondrial DNA, indicating a widespread male-specific phenomenon that focuses interest on the social structure of Bronze Age Europe.

Chimpanzee genomic diversity reveals ancient admixture with bonobos

Of chimpanzees and bonobos Modern non-African human genomes contain genomic remnants that suggest that there was interbreeding between ancient humans and archaic hominoid lineages. Now, de Manuel et al. show similar ancestral interbreeding between the ancestors of todays chimpanzees and bonobos (see the Perspective by Hoelzel). The study also provides population-specific genetic markers that may be valuable for conservation efforts. Science, this issue p. 477; see also p. 414 Genome sequences reveal ancient interbreeding between chimpanzees and bonobos. Our closest living relatives, chimpanzees and bonobos, have a complex demographic history. We analyzed the high-coverage whole genomes of 75 wild-born chimpanzees and bonobos from 10 countries in Africa. We found that chimpanzee population substructure makes genetic information a good predictor of geographic origin at country and regional scales. Multiple lines of evidence suggest that gene flow occurred from bonobos into the ancestors of central and eastern chimpanzees between 200,000 and 550,000 years ago, probably with subsequent spread into Nigeria-Cameroon chimpanzees. Together with another, possibly more recent contact (after 200,000 years ago), bonobos contributed less than 1% to the central chimpanzee genomes. Admixture thus appears to have been widespread during hominid evolution.

Decreased Nucleotide and Expression Diversity and Modified Coexpression Patterns Characterize Domestication in the Common Bean.

About 60% of the nucleotide diversity was lost during domestication of the common bean. The whole pattern of gene expression has also been affected, with changes in the patterns of coexpression among genes and 18% reduction in the overall diversity of gene expression. About 9% of the genes were selected during domestication, which is associated with further reduced diversity of expression. Using RNA sequencing technology and de novo transcriptome assembly, we compared representative sets of wild and domesticated accessions of common bean (Phaseolus vulgaris) from Mesoamerica. RNA was extracted at the first true-leaf stage, and de novo assembly was used to develop a reference transcriptome; the final data set consists of ∼190,000 single nucleotide polymorphisms from 27,243 contigs in expressed genomic regions. A drastic reduction in nucleotide diversity (∼60%) is evident for the domesticated form, compared with the wild form, and almost 50% of the contigs that are polymorphic were brought to fixation by domestication. In parallel, the effects of domestication decreased the diversity of gene expression (18%). While the coexpression networks for the wild and domesticated accessions demonstrate similar seminal network properties, they show distinct community structures that are enriched for different molecular functions. After simulating the demographic dynamics during domestication, we found that 9% of the genes were actively selected during domestication. We also show that selection induced a further reduction in the diversity of gene expression (26%) and was associated with 5-fold enrichment of differentially expressed genes. While there is substantial evidence of positive selection associated with domestication, in a few cases, this selection has increased the nucleotide diversity in the domesticated pool at target loci associated with abiotic stress responses, flowering time, and morphology.

Comparing models on the genealogical relationships among Neandertal, Cro-Magnoid and modern Europeans by serial coalescent simulations

Populations of anatomically archaic (Neandertal) and early modern (Cro-Magnoid) humans are jointly documented in the European fossil record, in the period between 40 000 and 25 000 years BP, but the large differences between their cultures, morphologies and DNAs suggest that the two groups were not close relatives. However, it is still unclear whether any genealogical continuity between them can be ruled out. Here, we simulated a broad range of demographic scenarios by means of a serial coalescence algorithm in which Neandertals, Cro-Magnoids and modern Europeans were either part of the same mitochondrial genealogy or of two separate genealogies. Mutation rates, population sizes, population structure and demographic growth rates varied across simulations. All models in which anatomically modern (that is, Cro-Magnoid and current) Europeans belong to a distinct genealogy performed better than any model in which the three groups were assigned to the same mitochondrial genealogy. The maximum admissible level of gene flow between Neandertals and the ancestors of current Europeans is 0.001% per generation, one order of magnitude lower than estimated in previous studies not considering genetic data on Cro-Magnoid people.

Genealogical discontinuities among Etruscan, Medieval and contemporary Tuscans

The available mitochondrial DNA (mtDNA) data do not point to clear genetic relationships between current Tuscans and the Bronze-Age inhabitants of Tuscany, the Etruscans. To understand how and when such a genetic discontinuity may have arisen, we extracted and typed the mtDNAs of 27 medieval Tuscans from an initial sample of 61, spanning a period between the 10th and 15th century AD. We then tested by serial coalescent simulation various models describing the genealogical relationships among past and current inhabitants of Tuscany, the latter including three samples (from Murlo, Volterra, and Casentino) that were recently claimed to be of Etruscan descent. Etruscans and medieval Tuscans share three mitochondrial haplotypes but fall in distinct branches of the mitochondrial genealogy in the only model that proved compatible with the data. Under that model, contemporary people of Tuscany show clear genetic relationships with Medieval people, but not with the Etruscans, along the female lines. No evidence of excess mutation was found in the Etruscan DNAs by a Bayesian test, and so there is no reason to suspect that these results are biased by systematic contamination of the ancient sequences or laboratory artefacts. Extensive demographic changes before AD 1000 are thus the simplest explanation for the differences between the contemporary and the Bronze-Age mtDNAs of Tuscany. Accordingly, genealogical continuity between ancient and modern populations of the same area does not seem a safe general assumption, but rather a hypothesis that, when possible, should be tested using ancient DNA analysis.

No evidence of Neandertal admixture in the mitochondrial genomes of early European modern humans and contemporary Europeans.

Neandertals, the archaic human form documented in Eurasia until 29,000 years ago, share no mitochondrial haplotype with modern Europeans. Whether this means that the two groups were reproductively isolated is controversial, and indeed nuclear data have been interpreted as suggesting that they admixed. We explored the range of demographic parameters that may have generated the observed mitochondrial diversity, simulating 3.0 million genealogies under six models differing as for the relationships among contemporary Europeans, Neandertals, and Upper Palaeolithic European early modern humans (EEMH), who coexisted with Neandertals for millennia. We compared by Approximate Bayesian Computations the simulation results with mitochondrial diversity in 7 Neandertals, 3 EEMH, and 150 opportunely chosen modern Europeans. A model of genealogical continuity between EEMH and contemporary Europeans, with no Neandertal contribution, received overwhelming support from the analyses. The maximum degree of Neandertal admixture, under the model of gene flow supported by nuclear data, was estimated at 1.5%, but this model proved 20-32 times less likely than a model without any gene flow. Nuclear and mitochondrial evidence might be reconciled if smaller population sizes led to faster lineage sorting for mitochondrial DNA, and Neandertals shared a longer period of common ancestry with the non-Africans than with the Africans ancestors.

Across language families: Genome diversity mirrors linguistic variation within Europe.

ABSTRACT Objectives: The notion that patterns of linguistic and biological variation may cast light on each other and on population histories dates back to Darwins times; yet, turning this intuition into a proper research program has met with serious methodological difficulties, especially affecting language comparisons. This article takes advantage of two new tools of comparative linguistics: a refined list of Indo‐European cognate words, and a novel method of language comparison estimating linguistic diversity from a universal inventory of grammatical polymorphisms, and hence enabling comparison even across different families. We corroborated the method and used it to compare patterns of linguistic and genomic variation in Europe. Materials and Methods: Two sets of linguistic distances, lexical and syntactic, were inferred from these data and compared with measures of geographic and genomic distance through a series of matrix correlation tests. Linguistic and genomic trees were also estimated and compared. A method (Treemix) was used to infer migration episodes after the main population splits. Results: We observed significant correlations between genomic and linguistic diversity, the latter inferred from data on both Indo‐European and non‐Indo‐European languages. Contrary to previous observations, on the European scale, language proved a better predictor of genomic differences than geography. Inferred episodes of genetic admixture following the main population splits found convincing correlates also in the linguistic realm. Discussion: These results pave the ground for previously unfeasible cross‐disciplinary analyses at the worldwide scale, encompassing populations of distant language families. Am J Phys Anthropol 157:630–640, 2015.