The population history of northeastern Siberia since the Pleistocene.

Northeastern Siberia has been inhabited by humans for more than 40,000 years but its deep population history remains poorly understood. Here we investigate the late Pleistocene population history of northeastern Siberia through analyses of 34 newly recovered ancient genomes that date to between 31,000 and 600 years ago. We document complex population dynamics during this period, including at least three major migration events: an initial peopling by a previously unknown Palaeolithic population of 'Ancient North Siberians' who are distantly related to early West Eurasian hunter-gatherers; the arrival of East Asian-related peoples, which gave rise to 'Ancient Palaeo-Siberians' who are closely related to contemporary communities from far-northeastern Siberia (such as the Koryaks), as well as Native Americans; and a Holocene migration of other East Asian-related peoples, who we name 'Neo-Siberians', and from whom many contemporary Siberians are descended. Each of these population expansions largely replaced the earlier inhabitants, and ultimately generated the mosaic genetic make-up of contemporary peoples who inhabit a vast area across northern Eurasia and the Americas.

A genomic history of Aboriginal Australia.

The population history of Aboriginal Australians remains largely uncharacterized. Here we generate high-coverage genomes for 83 Aboriginal Australians (speakers of Pama-Nyungan languages) and 25 Papuans from the New Guinea Highlands. We find that Papuan and Aboriginal Australian ancestors diversified 25-40 thousand years ago (kya), suggesting pre-Holocene population structure in the ancient continent of Sahul (Australia, New Guinea and Tasmania). However, all of the studied Aboriginal Australians descend from a single founding population that differentiated ~10-32 kya. We infer a population expansion in northeast Australia during the Holocene epoch (past 10,000 years) associated with limited gene flow from this region to the rest of Australia, consistent with the spread of the Pama-Nyungan languages. We estimate that Aboriginal Australians and Papuans diverged from Eurasians 51-72 kya, following a single out-of-Africa dispersal, and subsequently admixed with archaic populations. Finally, we report evidence of selection in Aboriginal Australians potentially associated with living in the desert.

Distance from sub-Saharan Africa predicts mutational load in diverse human genomes.

The Out-of-Africa (OOA) dispersal ∼ 50,000 y ago is characterized by a series of founder events as modern humans expanded into multiple continents. Population genetics theory predicts an increase of mutational load in populations undergoing serial founder effects during range expansions. To test this hypothesis, we have sequenced full genomes and high-coverage exomes from seven geographically divergent human populations from Namibia, Congo, Algeria, Pakistan, Cambodia, Siberia, and Mexico. We find that individual genomes vary modestly in the overall number of predicted deleterious alleles. We show via spatially explicit simulations that the observed distribution of deleterious allele frequencies is consistent with the OOA dispersal, particularly under a model where deleterious mutations are recessive. We conclude that there is a strong signal of purifying selection at conserved genomic positions within Africa, but that many predicted deleterious mutations have evolved as if they were neutral during the expansion out of Africa. Under a model where selection is inversely related to dominance, we show that OOA populations are likely to have a higher mutation load due to increased allele frequencies of nearly neutral variants that are recessive or partially recessive.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication.

The domestication of the horse ∼ 5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski's horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the "cost of domestication" hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.

Robust demographic inference from genomic and SNP data.

We introduce a flexible and robust simulation-based framework to infer demographic parameters from the site frequency spectrum (SFS) computed on large genomic datasets. We show that our composite-likelihood approach allows one to study evolutionary models of arbitrary complexity, which cannot be tackled by other current likelihood-based methods. For simple scenarios, our approach compares favorably in terms of accuracy and speed with ∂a∂i, the current reference in the field, while showing better convergence properties for complex models. We first apply our methodology to non-coding genomic SNP data from four human populations. To infer their demographic history, we compare neutral evolutionary models of increasing complexity, including unsampled populations. We further show the versatility of our framework by extending it to the inference of demographic parameters from SNP chips with known ascertainment, such as that recently released by Affymetrix to study human origins. Whereas previous ways of handling ascertained SNPs were either restricted to a single population or only allowed the inference of divergence time between a pair of populations, our framework can correctly infer parameters of more complex models including the divergence of several populations, bottlenecks and migration. We apply this approach to the reconstruction of African demography using two distinct ascertained human SNP panels studied under two evolutionary models. The two SNP panels lead to globally very similar estimates and confidence intervals, and suggest an ancient divergence (>110 Ky) between Yoruba and San populations. Our methodology appears well suited to the study of complex scenarios from large genomic data sets.

Evidence for polygenic adaptation to pathogens in the human genome.

Most approaches aiming at finding genes involved in adaptive events have focused on the detection of outlier loci, which resulted in the discovery of individually "significant" genes with strong effects. However, a collection of small effect mutations could have a large effect on a given biological pathway that includes many genes, and such a polygenic mode of adaptation has not been systematically investigated in humans. We propose here to evidence polygenic selection by detecting signals of adaptation at the pathway or gene set level instead of analyzing single independent genes. Using a gene-set enrichment test to identify genome-wide signals of adaptation among human populations, we find that most pathways globally enriched for signals of positive selection are either directly or indirectly involved in immune response. We also find evidence for long-distance genotypic linkage disequilibrium, suggesting functional epistatic interactions between members of the same pathway. Our results show that past interactions with pathogens have elicited widespread and coordinated genomic responses, and suggest that adaptation to pathogens can be considered as a primary example of polygenic selection.

Haplotype Inference Using Long-Read Nanopore Sequencing: Application to GSTA1 Promoter.

Recovering true haplotypes can have important clinical consequences. The laboratory process is difficult and is, therefore, most often done through inference. In this paper, we show that when using the Oxford nanopore sequencing technology, we could recover the true haplotypes of the GSTA1 promoter region. Eight LCL cell lines with potentially ambiguous haplotypes were used to characterize the efficacy of Oxford nanopore sequencing to phase the correct GSTA1 promoter haplotypes. The results were compared to Sanger sequencing and inferred haplotypes in the 1000 genomes project. The average read length was 813 bp out of a total PCR length of 1336 bp. The best coverage of sequencing was in the middle of the PCR product and decreased to 50% at the PCR ends. SNPs separated by less than 200 bp showed > 90% of correct haplotypes, while at the distance of 1089 bp, this proportion still exceeded 58%. The number of cycles influences the generation of hybrid haplotypes but not extension or annealing time. The results demonstrate that this long sequencing reads methodology, can accurately determine the haplotypes without the need for inference. The technology proved to be robust but the success of phasing nonetheless depends on the distances and frequencies of SNPs.

17q Gain in Neuroblastoma: A Review of Clinical and Biological Implications.

Neuroblastoma (NB) is the most frequent extracranial solid childhood tumor. Despite advances in the understanding and treatment of this disease, the prognosis in cases of high-risk NB is still poor. 17q gain has been shown to be the most frequent genomic alteration in NB. However, the significance of this remains unclear because of its high frequency and association with other genetic modifications, particularly segmental chromosomal aberrations, 1p and 11q deletions, and MYCN amplification, all of which are also associated with a poor clinical prognosis. This work reviewed the evidence on the clinical and biological significance of 17q gain. It strongly supports the significance of 17q gain in the development of NB and its importance as a clinically relevant marker. However, it is crucial to distinguish between whole and partial chromosome 17q gains. The most important breakpoints appear to be at 17q12 and 17q21. The former distinguishes between whole and partial chromosome 17q gain; the latter is a site of IGF2BP1 and NME1 genes that appear to be the main oncogenes responsible for the functional effects of 17q gain.

A novel integrative multi-omics approach to unravel the genetic determinants of rare diseases with application in sinusoidal obstruction syndrome.

BACKGROUND: Genotype-phenotype analyses of rare diseases often suffer from a lack of power, due to small sample size, which makes identifying significant associations difficult. Sinusoidal obstruction syndrome (SOS) of the liver is a rare but life-threatening complication of hematopoietic stem cell transplantation (HSCT). The alkylating agent busulfan is commonly used in HSCT and known to trigger SOS. We developed a novel pipeline to identify genetic determinants in rare diseases by combining in vitro information with clinical whole-exome sequencing (WES) data and applied it in SOS patients and controls. METHODS: First, we analysed differential gene expression in six lymphoblastoid cell lines (LCLs) before and after incubation with busulfan. Second, we used WES data from 87 HSCT patients and estimated the association with SOS at the SNP and the gene levels. We then combined the results of the expression and the association analyses into an association statistic at the gene level. We used an over-representation analysis to functionally characterize the genes that were associated with a significant combined test statistic. RESULTS: After treatment of LCLs with busulfan, 1708 genes were significantly up-, and 1385 down-regulated. The combination of the expression experiment and the association analysis of WES data into a single test statistic revealed 35 genes associated with the outcome. These genes are involved in various biological functions and processes, such as "Cell growth and death", "Signalling molecules and interaction", "Cancer", and "Infectious disease". CONCLUSIONS: This novel data analysis pipeline integrates two independent omics datasets and increases statistical power for identifying genotype-phenotype associations. The analysis of the transcriptomics profile of cell lines treated with busulfan and WES data from HSCT patients allowed us to identify potential genetic contributors to SOS. Our pipeline could be useful for identifying genetic contributors to other rare diseases where limited power renders genome-wide analyses unpromising. TRIAL REGISTRATION: For the clinical dataset: Clinicaltrials.gov: NCT01257854. https://clinicaltrials.gov/ct2/history/NCT01257854.

Differentiation of IL-26+ TH17 intermediates into IL-17A producers via epithelial crosstalk in psoriasis.

Interleukin (IL)-26 is a TH17 cytokine with known antimicrobial and pro-inflammatory functions. However, the precise role of IL-26 in the context of pathogenic TH17 responses is unknown. Here we identify a population of blood TH17 intermediates that produce high levels of IL-26 and differentiate into IL-17A-producing TH17 cells upon TGF-ß1 exposure. By combining single cell RNA sequencing, TCR sequencing and spatial transcriptomics we show that this process occurs in psoriatic skin. In fact, IL-26+ TH17 intermediates infiltrating psoriatic skin induce TGF-ß1 expression in basal keratinocytes and thereby promote their own differentiation into IL-17A-producing cells. Thus, our study identifies IL-26-producing cells as an early differentiation stage of TH17 cells that infiltrates psoriatic skin and controls its own maturation into IL17A-producing TH17 cells, via epithelial crosstalk involving paracrine production of TGF-ß1.

Apelin-driven endothelial cell migration sustains intestinal progenitor cells and tumor growth.

Stem and progenitor cells residing in the intestinal crypts drive the majority of colorectal cancers (CRCs), yet vascular contribution to this niche remains largely unexplored. VEGFA is a key driver of physiological and tumor angiogenesis. Accordingly, current anti-angiogenic cancer therapies target the VEGFA pathway. Here we report that in CRC expansion of the stem/progenitor pool in intestinal crypts requires VEGFA-independent growth and remodeling of blood vessels. Epithelial transformation induced expression of the endothelial peptide apelin, directs migration of distant venous endothelial cells towards progenitor niche vessels ensuring optimal perfusion. In the absence of apelin, loss of injury-inducible PROX1+ epithelial progenitors inhibited both incipient and advanced intestinal tumor growth. Our results establish fundamental principles for the reciprocal communication between vasculature and the intestinal progenitor niche and provide a mechanism for resistance to VEGFA-targeting drugs in CRCs.

The analysis of GSTA1 promoter genetic and functional diversity of human populations.

GSTA1 encodes a member of a family of enzymes that function to add glutathione to target electrophilic compounds, including carcinogens, therapeutic drugs, environmental toxins, and products of oxidative stress. GSTA1 has several functional SNPs within its promoter region that are responsible for a change in its expression by altering promoter function. This study aims to investigate distributions of GSTA1 promoter haplotypes across different human populations and to assess their impact on the expression of GSTA1. PHASE 2.1.1 was used to infer haplotypes and diplotypes of six GSTA1 promoter SNPs on 2501 individuals from 26 populations classified by the 1000 Genomes Project into five super-populations that included Africa (N = 660), America (N = 347), East Asia (N = 504), Europe (N = 502), and South Asia (N = 488). We used pairwise FST analysis to compare sub-populations and luciferase reporter assay (LRA) to evaluate the impact of each SNP on activation of transcription and interaction with other SNPs. The distributions of GSTA1 promoter haplotypes and diplotypes were significantly different among the different human populations. Three new promoter haplotypes were found in the African super-population. LRA demonstrated that SNPs at -52 and -69 has the most impact on GSTA1 expression, however other SNPs have a significant impact on transcriptional activity. Based on LRA, a new model of cis-elements interaction is presented. Due to the significant differences in GSTA1 diplotype population frequencies, future pharmacogenomics or disease-related studies would benefit from the inclusion of the complete GSTA1 promoter haplotype based on the newly proposed metabolic grouping derived from the LRA results.

Classifying dementia progression using microbial profiling of saliva.

INTRODUCTION: There is increasing evidence linking periodontal infections to Alzheimer's disease (AD). Saliva sampling can reveal information about the host and pathogen interactions that can inform about physiological and pathological brain states. METHODS: A cross-sectional cohort of age-matched participants (78) was segmented according to their chemosensory (University of Pennsylvania Smell Identification Test; UPSIT) and cognitive scores (Mini-Mental State Exam; MMSE and clinical dementia rating; CDR). Mid-morning saliva was sampled from each participant and processed for microbiome composition and cytokine analysis. Linear discriminant analysis (LDA) was used to unravel specific changes in microbial and immunological signatures and logistic regression analysis (LRA) was employed to identify taxa that varied in abundance among patient groups. RESULTS: Using olfaction we distinguish in the cognitively normal population a segment with high chemosensory scores (CNh, 27) and another segment with chemosensory scores (CNr, 16) as low as mild cognitive impairment (MCI, 21) but higher than the AD group (17). We could identify stage-specific microbial signatures changes but no clearly distinct cytokine profiles. Periodontal pathogen species as Filifactor villosus decline with the increasing severity of AD, whereas opportunistic oral bacteria such as Leptotrichia wadei show a significant enrichment in MCI. CONCLUSIONS: The salivary microbiome indicates stage-dependent changes in oral bacteria favoring opportunistic species at the expense of periodontal bacteria, whereas the inflammatory profiles remain mainly unchanged in the sampled population.

Emergence of young human genes after a burst of retroposition in primates.

The origin of new genes through gene duplication is fundamental to the evolution of lineage- or species-specific phenotypic traits. In this report, we estimate the number of functional retrogenes on the lineage leading to humans generated by the high rate of retroposition (retroduplication) in primates. Extensive comparative sequencing and expression studies coupled with evolutionary analyses and simulations suggest that a significant proportion of recent retrocopies represent bona fide human genes. We estimate that at least one new retrogene per million years emerged on the human lineage during the past approximately 63 million years of primate evolution. Detailed analysis of a subset of the data shows that the majority of retrogenes are specifically expressed in testis, whereas their parental genes show broad expression patterns. Consistently, most retrogenes evolved functional roles in spermatogenesis. Proteins encoded by X chromosome-derived retrogenes were strongly preserved by purifying selection following the duplication event, supporting the view that they may act as functional autosomal substitutes during X-inactivation of late spermatogenesis genes. Also, some retrogenes acquired a new or more adapted function driven by positive selection. We conclude that retroduplication significantly contributed to the formation of recent human genes and that most new retrogenes were progressively recruited during primate evolution by natural and/or sexual selection to enhance male germline function.

Relaxed Selection During a Recent Human Expansion.

Humans have colonized the planet through a series of range expansions, which deeply impacted genetic diversity in newly settled areas and potentially increased the frequency of deleterious mutations on expanding wave fronts. To test this prediction, we studied the genomic diversity of French Canadians who colonized Quebec in the 17th century. We used historical information and records from ∼4000 ascending genealogies to select individuals whose ancestors lived mostly on the colonizing wave front and individuals whose ancestors remained in the core of the settlement. Comparison of exomic diversity reveals that: (i) both new and low-frequency variants are significantly more deleterious in front than in core individuals, (ii) equally deleterious mutations are at higher frequencies in front individuals, and (iii) front individuals are two times more likely to be homozygous for rare very deleterious mutations present in Europeans. These differences have emerged in the past six to nine generations and cannot be explained by differential inbreeding, but are consistent with relaxed selection mainly due to higher rates of genetic drift on the wave front. Demographic inference and modeling of the evolution of rare variants suggest lower effective size on the front, and lead to an estimation of selection coefficients that increase with conservation scores. Even though range expansions have had a relatively limited impact on the overall fitness of French Canadians, they could explain the higher prevalence of recessive genetic diseases in recently settled regions of Quebec.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

Accumulation of Deleterious Mutations During Bacterial Range Expansions.

Recent theory predicts that the fitness of pioneer populations can decline when species expand their range, due to high rates of genetic drift on wave fronts making selection less efficient at purging deleterious variants. To test these predictions, we studied the fate of mutator bacteria expanding their range for 1650 generations on agar plates. In agreement with theory, we find that growth abilities of strains with a high mutation rate (HMR lines) decreased significantly over time, unlike strains with a lower mutation rate (LMR lines) that present three to four times fewer mutations. Estimation of the distribution of fitness effect under a spatially explicit model reveals a mean negative effect for new mutations (-0.38%), but it suggests that both advantageous and deleterious mutations have accumulated during the experiment. Furthermore, the fitness of HMR lines measured in different environments has decreased relative to the ancestor strain, whereas that of LMR lines remained unchanged. Contrastingly, strains with a HMR evolving in a well-mixed environment accumulated less mutations than agar-evolved strains and showed an increased fitness relative to the ancestor. Our results suggest that spatially expanding species are affected by deleterious mutations, leading to a drastic impairment of their evolutionary potential.

Ancient genomes show social and reproductive behavior of early Upper Paleolithic foragers.

Present-day hunter-gatherers (HGs) live in multilevel social groups essential to sustain a population structure characterized by limited levels of within-band relatedness and inbreeding. When these wider social networks evolved among HGs is unknown. To investigate whether the contemporary HG strategy was already present in the Upper Paleolithic, we used complete genome sequences from Sunghir, a site dated to ~34,000 years before the present, containing multiple anatomically modern human individuals. We show that individuals at Sunghir derive from a population of small effective size, with limited kinship and levels of inbreeding similar to HG populations. Our findings suggest that Upper Paleolithic social organization was similar to that of living HGs, with limited relatedness within residential groups embedded in a larger mating network.

Genetic surfing in human populations: from genes to genomes.

Genetic surfing describes the spatial spread and increase in frequency of variants that are not lost by genetic drift and serial migrant sampling during a range expansion. Genetic surfing does not modify the total number of derived alleles in a population or in an individual genome, but it leads to a loss of heterozygosity along the expansion axis, implying that derived alleles are more often in homozygous state. Genetic surfing also affects selected variants on the wave front, making them behave almost like neutral variants during the expansion. In agreement with theoretical predictions, human genomic data reveals an increase in recessive mutation load with distance from Africa, an expansion load likely to have developed during the expansions of human populations out of Africa.

Chimpanzee genomic diversity reveals ancient admixture with 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.