Francesca Luca

Maternal traces of deep common ancestry and asymmetric gene flow between Pygmy hunter–gatherers and Bantu-speaking farmers

Two groups of populations with completely different lifestyles—the Pygmy hunter–gatherers and the Bantu-speaking farmers—coexist in Central Africa. We investigated the origins of these two groups and the interactions between them, by analyzing mtDNA variation in 1,404 individuals from 20 farming populations and 9 Pygmy populations from Central Africa, with the aim of shedding light on one of the most fascinating cultural transitions in human evolution (the transition from hunting and gathering to agriculture). Our data indicate that this region was colonized gradually, with an initial L1c-rich ancestral population ultimately giving rise to current-day farmers, who display various L1c clades, and to Pygmies, in whom L1c1a is the only surviving clade. Detailed phylogenetic analysis of complete mtDNA sequences for L1c1a showed this clade to be autochthonous to Central Africa, with its most recent branches shared between farmers and Pygmies. Coalescence analyses revealed that these two groups arose through a complex evolutionary process characterized by (i) initial divergence of the ancestors of contemporary Pygmies from an ancestral Central African population no more than ≈70,000 years ago, (ii) a period of isolation between the two groups, accounting for their phenotypic differences, (iii) long-standing asymmetric maternal gene flow from Pygmies to the ancestors of the farming populations, beginning no more than ≈40,000 years ago and persisting until a few thousand years ago, and (iv) enrichment of the maternal gene pool of the ancestors of the farming populations by the arrival and/or subsequent demographic expansion of L0a, L2, and L3 carriers.

Y chromosomal haplogroup J as a signature of the post-neolithic colonization of Europe

In order to attain a finer reconstruction of the peopling of southern and central-eastern Europe from the Levant, we determined the frequencies of eight lineages internal to the Y chromosomal haplogroup J, defined by biallelic markers, in 22 population samples obtained with a fine-grained sampling scheme. Our results partially resolve a major multifurcation of lineages within the haplogroup. Analyses of molecular variance show that the area covered by haplogroup J dispersal is characterized by a significant degree of molecular radiation for unique event polymorphisms within the haplogroup, with a higher incidence of the most derived sub-haplogroups on the northern Mediterranean coast, from Turkey westward; here, J diversity is not simply a subset of that present in the area in which this haplogroup first originated. Dating estimates, based on simple tandem repeat loci (STR) diversity within each lineage, confirmed the presence of a major population structuring at the time of spread of haplogroup J in Europe and a punctuation in the peopling of this continent in the post-Neolithic, compatible with the expansion of the Greek world. We also present here, for the first time, a novel method for comparative dating of lineages, free of assumptions of STR mutation rates.

The Timing of Selection at the Human FOXP2 Gene

Krause J, Lalueza-Fox C, Orlando L, et al. recently examined patterns of genetic variation at FOXP2 in 2 Neanderthals. This gene is of particular interest because it is involved in speech and language and was previously shown to harbor the signature of recent positive selection. The authors found the same 2 amino acid substitutions in Neanderthals as in modern humans. Assuming that these sites were the targets of selection and no interbreeding between the 2 groups, they concluded that selection at FOXP2 occurred before the populations split, over 300 thousand years ago. Here, we show that the data are unlikely under this scenario but may instead be consistent with low rates of gene flow between modern humans and Neanderthals. We also collect additional data and introduce a modeling framework to estimate levels of modern human contamination of the Neanderthal samples. We find that, depending on the assumptions, additional control experiments may be needed to rule out contamination at FOXP2.

Evolutionary Adaptations to Dietary Changes

Through cultural innovation and changes in habitat and ecology, there have been a number of major dietary shifts in human evolution, including meat eating, cooking, and those associated with plant and animal domestication. The identification of signatures of adaptations to such dietary changes in the genome of extant primates (including humans) may shed light not only on the evolutionary history of our species, but also on the mechanisms that underlie common metabolic diseases in modern human populations. In this review, we provide a brief overview of the major dietary shifts that occurred during hominin evolution, and we discuss the methods and approaches used to identify signals of natural selection in patterns of sequence variation. We then review the results of studies aimed at detecting the genetic loci that played a major role in dietary adaptations and conclude by outlining the potential of future studies in this area.

Interactions between Glucocorticoid Treatment and Cis-Regulatory Polymorphisms Contribute to Cellular Response Phenotypes

Glucocorticoids (GCs) mediate physiological responses to environmental stress and are commonly used as pharmaceuticals. GCs act primarily through the GC receptor (GR, a transcription factor). Despite their clear biomedical importance, little is known about the genetic architecture of variation in GC response. Here we provide an initial assessment of variability in the cellular response to GC treatment by profiling gene expression and protein secretion in 114 EBV-transformed B lymphocytes of African and European ancestry. We found that genetic variation affects the response of nearby genes and exhibits distinctive patterns of genotype-treatment interactions, with genotypic effects evident in either only GC-treated or only control-treated conditions. Using a novel statistical framework, we identified interactions that influence the expression of 26 genes known to play central roles in GC-related pathways (e.g. NQO1, AIRE, and SGK1) and that influence the secretion of IL6.

Clinal patterns of human Y chromosomal diversity in continental Italy and Greece are dominated by drift and founder effects.

We explored the spatial distribution of human Y chromosomal diversity on a microgeographic scale, by typing 30 population samples from closely spaced locations in Italy and Greece for 9 haplogroups and their internal microsatellite variation. We confirm a significant difference in the composition of the Y chromosomal gene pools of the two countries. However, within each country, heterogeneity is not organized along the lines of clinal variation deduced from studies on larger spatial scales. Microsatellite data indicate that local increases of haplogroup frequencies can be often explained by a limited number of founders. We conclude that local founder or drift effects are the main determinants in shaping the microgeographic Y chromosomal diversity.

Bacterial infection remodels the DNA methylation landscape of human dendritic cells

DNA methylation is an epigenetic mark thought to be robust to environmental perturbations on a short time scale. Here, we challenge that view by demonstrating that the infection of human dendritic cells (DCs) with a live pathogenic bacteria is associated with rapid and active demethylation at thousands of loci, independent of cell division. We performed an integrated analysis of data on genome-wide DNA methylation, histone mark patterns, chromatin accessibility, and gene expression, before and after infection. We found that infection-induced demethylation rarely occurs at promoter regions and instead localizes to distal enhancer elements, including those that regulate the activation of key immune transcription factors. Active demethylation is associated with extensive epigenetic remodeling, including the gain of histone activation marks and increased chromatin accessibility, and is strongly predictive of changes in the expression levels of nearby genes. Collectively, our observations show that active, rapid changes in DNA methylation in enhancers play a previously unappreciated role in regulating the transcriptional response to infection, even in nonproliferating cells.

Multiple Advantageous Amino Acid Variants in the NAT2 Gene in Human Populations

Background Genetic variation at NAT2 has been long recognized as the cause of differential ability to metabolize a wide variety of drugs of therapeutic use. Here, we explore the pattern of genetic variation in 12 human populations that significantly extend the geographic range and resolution of previous surveys, to test the hypothesis that different dietary regimens and lifestyles may explain inter-population differences in NAT2 variation. Methodology/Principal Findings The entire coding region was resequenced in 98 subjects and six polymorphic positions were genotyped in 150 additional subjects. A single previously undescribed variant was found (34T>C; 12Y>H). Several aspects of the data do not fit the expectations of a neutral model, as assessed by coalescent simulations. Tajimas D is positive in all populations, indicating an excess of intermediate alleles. The level of between-population differentiation is low, and is mainly accounted for by the proportion of fast vs. slow acetylators. However, haplotype frequencies significantly differ across groups of populations with different subsistence. Conclusions/Significance Data on the structure of haplotypes and their frequencies are compatible with a model in which slow-causing variants were present in widely dispersed populations before major shifts to pastoralism and/or agriculture. In this model, slow-causing mutations gained a selective advantage in populations shifting from hunting-gathering to pastoralism/agriculture. We suggest the diminished dietary availability of folates resulting from the nutritional shift, as the possible cause of the fitness increase associated to haplotypes carrying mutations that reduce enzymatic activity.

A reduced representation approach to population genetic analyses and applications to human evolution

Second-generation sequencing technologies allow surveys of sequence variation on an unprecedented scale. However, despite the rapid decrease in sequencing costs, collecting whole-genome sequence data on a population scale is still prohibitive for many laboratories. We have implemented an inexpensive, reduced representation protocol for preparing resequencing targets, and we have developed the analytical tools necessary for making population genetic inferences. This approach can be applied to any species for which a draft or complete reference genome sequence is available. The new tools we have developed include methods for aligning reads, calling genotypes, and incorporating sample-specific sequencing error rates in the estimate of evolutionary parameters. When applied to 19 individuals from a total of 18 human populations, our approach allowed sampling regions that are largely overlapping across individuals and that are representative of the entire genome. The resequencing data were used to test the serial founder model of human dispersal and to estimate the time of the Out of Africa migration. Our results also represent the first attempt to provide a time frame for the colonization of Australia based on large-scale resequencing data.

Adaptive Variation Regulates the Expression of the Human SGK1 Gene in Response to Stress

The Serum and Glucocorticoid-regulated Kinase1 (SGK1) gene is a target of the glucocorticoid receptor (GR) and is central to the stress response in many human tissues. Because environmental stress varies across habitats, we hypothesized that natural selection shaped the geographic distribution of genetic variants regulating the level of SGK1 expression following GR activation. By combining population genetics and molecular biology methods, we identified a variant (rs9493857) with marked allele frequency differences between populations of African and European ancestry and with a strong correlation between allele frequency and latitude in worldwide population samples. This SNP is located in a GR-binding region upstream of SGK1 that was identified using a GR ChIP-chip. SNP rs9493857 also lies within a predicted binding site for Oct1, a transcription factor known to cooperate with the GR in the transactivation of target genes. Using ChIP assays, we show that both GR and Oct1 bind to this region and that the ancestral allele at rs9493857 binds the GR-Oct1 complex more efficiently than the derived allele. Finally, using a reporter gene assay, we demonstrate that the ancestral allele is associated with increased glucocorticoid-dependent gene expression when compared to the derived allele. Our results suggest a novel paradigm in which hormonal responsiveness is modulated by sequence variation in the regulatory regions of nuclear receptor target genes. Identifying such functional variants may shed light on the mechanisms underlying inter-individual variation in response to environmental stressors and to hormonal therapy, as well as in the susceptibility to hormone-dependent diseases.