Investigating the role of demography and selection in genome scale patterns of common and rare variant diversity in humans

In the last decade, an unprecedented increase in the availability of whole genome sequence (WGS) data has reshaped the field of human evolutionary genomics. However, many earlier sequencing projects like the HapMap and 1000 Genomes panels focussed on a limited set of populations. Therefore, more research has been required to better characterise genetic diversity in understudied regions, such as Island Southeast Asia and Siberia. This thesis contributes to this ongoing effort in the form of three partially related subprojects. Firstly, population structure and local adaptations in Southeast Asia were investigated using novel autosomal 730,000 SNP data from 146 individuals in the context of a larger worldwide panel of 1,825 humans. The Kankanaey Igorot from the highlands of the Philippine Mountain Province were highlighted as the closest living representatives of the source population that may have given rise to the Austronesian expansion. Furthermore, consistent with archaeological, cultural and linguistic evidence of Indian influence in Southeast Asia starting from 2.5 kya South Asian admixture in the region was estimated to date back to the last couple of thousand years. To provide an unbiased high-resolution picture of the patterns of functional and rare variants worldwide high coverage WGS data from 483 individuals (including 379 novel genomes) were analysed. Ingenuity Variant Analysis and the Ensembl Variant Effect Predictor were applied to a subset of these genomes (n = 382) to create a repository of functional and deleterious variants. Evidence for purifying selection in genes involved in pigmentation and immune defence against viruses was detected in African populations. The most differentiated sites across continental groups were integrated with haplotype-based selection tests and annotations from functional databases to pinpoint disease and metabolism-related candidate loci. A subset of the WGS dataset, designed to maximise coverage of diverse ethnic groups (n = 447), was screened for variants occurring exclusively in two individuals in a heterozygous state (f2 variants). It was shown that f2 sharing correlates well with the results of CHROMOPAINTER, a state-of-the-art method to detect recent gene flow, and, allows for the detection of cryptic relatedness among distant populations. This was demonstrated by an example of a previously undetected low-scale African ancestry component in the South American Calchaquíes putatively related to the transatlantic slave trade

Evidence of Early-Stage Selection on EPAS1 and GPR126 Genes in Andean High Altitude Populations.

The aim of this study is to identify genetic variants that harbour signatures of recent positive selection and may facilitate physiological adaptations to hypobaric hypoxia. To achieve this, we conducted whole genome sequencing and lung function tests in 19 Argentinean highlanders (>3500 m) comparing them to 16 Native American lowlanders. We developed a new statistical procedure using a combination of population branch statistics (PBS) and number of segregating sites by length (nSL) to detect beneficial alleles that arose since the settlement of the Andes and are currently present in 15-50% of the population. We identified two missense variants as significant targets of selection. One of these variants, located within the GPR126 gene, has been previously associated with the forced expiratory volume/forced vital capacity ratio. The other novel missense variant mapped to the EPAS1 gene encoding the hypoxia inducible factor 2α. EPAS1 is known to be the major selection candidate gene in Tibetans. The derived allele of GPR126 is associated with lung function in our sample of highlanders (p < 0.05). These variants may contribute to the physiological adaptations to hypobaric hypoxia, possibly by altering lung function. The new statistical approach might be a useful tool to detect selected variants in population studies

Genomic analyses inform on migration events during the peopling of Eurasia.

High-coverage whole-genome sequence studies have so far focused on a limited number of geographically restricted populations, or been targeted at specific diseases, such as cancer. Nevertheless, the availability of high-resolution genomic data has led to the development of new methodologies for inferring population history and refuelled the debate on the mutation rate in humans. Here we present the Estonian Biocentre Human Genome Diversity Panel (EGDP), a dataset of 483 high-coverage human genomes from 148 populations worldwide, including 379 new genomes from 125 populations, which we group into diversity and selection sets. We analyse this dataset to refine estimates of continent-wide patterns of heterozygosity, long- and short-distance gene flow, archaic admixture, and changes in effective population size through time as well as for signals of positive or balancing selection. We find a genetic signature in present-day Papuans that suggests that at least 2% of their genome originates from an early and largely extinct expansion of anatomically modern humans (AMHs) out of Africa. Together with evidence from the western Asian fossil record, and admixture between AMHs and Neanderthals predating the main Eurasian expansion, our results contribute to the mounting evidence for the presence of AMHs out of Africa earlier than 75,000 years ago.Support was provided by: Estonian Research Infrastructure Roadmap grant no 3.2.0304.11-0312; Australian Research Council Discovery grants (DP110102635 and DP140101405) (D.M.L., M.W. and E.W.); Danish National Research Foundation; the Lundbeck Foundation and KU2016 (E.W.); ERC Starting Investigator grant (FP7 - 261213) (T.K.); Estonian Research Council grant PUT766 (G.C. and M.K.); EU European Regional Development Fund through the Centre of Excellence in Genomics to Estonian Biocentre (R.V.; M.Me. and A.Me.), and Centre of Excellence for Genomics and Translational Medicine Project No. 2014-2020.4.01.15-0012 to EGC of UT (A.Me.) and EBC (M.Me.); Estonian Institutional Research grant IUT24-1 (L.S., M.J., A.K., B.Y., K.T., C.B.M., Le.S., H.Sa., S.L., D.M.B., E.M., R.V., G.H., M.K., G.C., T.K. and M.Me.) and IUT20-60 (A.Me.); French Ministry of Foreign and European Affairs and French ANR grant number ANR-14-CE31-0013-01 (F.-X.R.); Gates Cambridge Trust Funding (E.J.); ICG SB RAS (No. VI.58.1.1) (D.V.L.); Leverhulme Programme grant no. RP2011-R-045 (A.B.M., P.G. and M.G.T.); Ministry of Education and Science of Russia; Project 6.656.2014/K (S.A.F.); NEFREX grant funded by the European Union (People Marie Curie Actions; International Research Staff Exchange Scheme; call FP7-PEOPLE-2012-IRSES-number 318979) (M.Me., G.H. and M.K.); NIH grants 5DP1ES022577 05, 1R01DK104339-01, and 1R01GM113657-01 (S.Tis.); Russian Foundation for Basic Research (grant N 14-06-00180a) (M.G.); Russian Foundation for Basic Research; grant 16-04-00890 (O.B. and E.B); Russian Science Foundation grant 14-14-00827 (O.B.); The Russian Foundation for Basic Research (14-04-00725-a), The Russian Humanitarian Scientific Foundation (13-11-02014) and the Program of the Basic Research of the RAS Presidium “Biological diversity” (E.K.K.); Wellcome Trust and Royal Society grant WT104125AIA & the Bristol Advanced Computing Research Centre (http://www.bris.ac.uk/acrc/) (D.J.L.); Wellcome Trust grant 098051 (Q.A.; C.T.-S. and Y.X.); Wellcome Trust Senior Research Fellowship grant 100719/Z/12/Z (M.G.T.); Young Explorers Grant from the National Geographic Society (8900-11) (C.A.E.); ERC Consolidator Grant 647787 ‘LocalAdaptatio’ (A.Ma.); Program of the RAS Presidium “Basic research for the development of the Russian Arctic” (B.M.); Russian Foundation for Basic Research grant 16-06-00303 (E.B.); a Rutherford Fellowship (RDF-10-MAU-001) from the Royal Society of New Zealand (M.P.C.)

Same old in Middle Neolithic diets? A stable isotope study of bone collagen from the burial community of Jechtingen,southwest Germany

Dietary patterns and subsistence strategies of the early farmers of the central European Linearbandkeramik culture (LBK) have received much research attention, but only little is known about subsequent Neolithic time periods. In order to address continuity or alterations through time, this study explores bone collagen carbon and nitrogen isotope data from the Middle Neolithic cemetery of Jechtingen in southwest Germany, attributed to the ceramic phases of Großgartach, Planig-Friedberg and Rössen (4850/4800–4500 BC). The dataset comprises 77 individuals making it the most extensive from this period to date. The average values of the adults are −20.6 ± 0.3‰ for δ13C and 9.4 ± 0.8‰ for δ15N (n = 62) indicating a mixed diet of animal and plant proteins based on a terrestrial C3 ecosystem. There is no indication of dietary differences due to age or sex. The δ13C values follow regional trends emerging from previously investigated LBK and Middle Neolithic sites. They appear to be primarily driven by environmental factors, such as annual precipitation rates. Concerning δ15N, LBK and Middle Neolithic datasets overlap widely, and Jechtingen falls well into the data ranges of other previously published sites. The stable isotope data do not hint at remarkable changes of the dietary habits between the food-producing communities of the LBK and the subsequent Middle Neolithic

Response to Wyckelsma et&nbsp;al.: Loss of α-actinin-3 during human evolution provides superior cold resilience and muscle heat generation

The common loss-of-function mutation R577X in the structural muscle protein ACTN3 emerged as a potential target of positive selection from early studies and has been the focus of insightful physiological work suggesting a significant impact on muscle metabolism. Adaptation to cold climates has been proposed as a key adaptive mechanism explaining its global allele frequency patterns. Here, we re-examine this hypothesis analyzing modern (n = 3,626) and ancient (n = 1,651) genomic data by using allele-frequency as well as haplotype homozygosity-based methods. The presented results are more consistent with genetic drift rather than selection in cold climates as the main driver of the ACTN3 R577X frequency distribution in human populations across the world. This Matters Arising paper is in response to Wyckelsma et al. (2021),(1) published in The American Journal of Human Genetics. See also the response by Wyckelsma et al. (2022),(2) published in this issue

The Common H202D Variant in GDF-15 Does Not Affect Its Bioactivity but Can Significantly Interfere with Measurement of Its Circulating Levels.

BACKGROUND: There is growing interest in the measurement of growth differentiation factor 15 (GDF-15) in a range of disorders associated with cachexia. We undertook studies to determine whether a common histidine (H) to aspartate (D) variant at position 202 in the pro-peptide (position 6 in the mature peptide) interfered with its detection by 3 of the most commonly used immunoassays. METHODS: Three synthetic GDF-15-forms (HH homo-, HD hetero-, and DD-homodimers) were measured after serial dilution using Roche Elecsys®, R&amp;D QuantikineTM ELISA, and MSD R&amp;D DuoSet® immunoassays. GDF-15 concentrations were measured by the Roche and the MSD R&amp;D immunoassays in 173 genotyped participants (61 HH homozygotes, 59 HD heterozygotes, and 53 DD homozygotes). For the comparative statistical analyses of the GDF-15 concentrations, we used non-parametric tests, in particular Bland-Altman difference (bias) plots and Passing-Bablok regression. The bioactivity of the 2 different homodimers was compared in a cell-based assay in HEK293S-SRF-RET/GFRAL cells. RESULTS: The Roche assay detected H- and D-containing peptides similarly but the R&amp;D reagents (Quantikine and DuoSet) consistently underreported GDF-15 concentrations in the presence of the D variant. DD dimers had recoveries of approximately 45% while HD dimers recoveries were 62% to 78%. In human serum samples, the GDF-15 concentrations reported by the R&amp;D assay were a median of 4% lower for HH, a median of 36% lower for HD, and a median of 61% lower for DD compared to the Roche assay. The bioactivities of the HH and DD peptides were indistinguishable. CONCLUSIONS: The D variant of GDF-15 substantially affects its measurement by a commonly used immunoassay, a finding that has clear implications for its interpretation in research and clinical settings

Loss-of-function mutations in the melanocortin 4 receptor in a UK birth cohort.

Mutations in the melanocortin 4 receptor gene (MC4R) are associated with obesity but little is known about the prevalence and impact of such mutations throughout human growth and development. We examined the MC4R coding sequence in 5,724 participants from the Avon Longitudinal Study of Parents and Children, functionally characterized all nonsynonymous MC4R variants and examined their association with anthropometric phenotypes from childhood to early adulthood. The frequency of heterozygous loss-of-function (LoF) mutations in MC4R was ~1 in 337 (0.30%), considerably higher than previous estimates. At age 18 years, mean differences in body weight, body mass index and fat mass between carriers and noncarriers of LoF mutations were 17.76 kg (95% CI 9.41, 26.10), 4.84 kg m-2 (95% CI 2.19, 7.49) and 14.78 kg (95% CI 8.56, 20.99), respectively. MC4R LoF mutations may be more common than previously reported and carriers of such variants may enter adult life with a substantial burden of excess adiposity

Penetrance of pathogenic genetic variants associated with premature ovarian insufficiency.

Premature ovarian insufficiency (POI) affects 1% of women and is a leading cause of infertility. It is often considered to be a monogenic disorder, with pathogenic variants in ~100 genes described in the literature. We sought to systematically evaluate the penetrance of variants in these genes using exome sequence data in 104,733 women from the UK Biobank, 2,231 (1.14%) of whom reported at natural menopause under the age of 40 years. We found limited evidence to support any previously reported autosomal dominant effect. For nearly all heterozygous effects on previously reported POI genes, we ruled out even modest penetrance, with 99.9% (13,699 out of 13,708) of all protein-truncating variants found in reproductively healthy women. We found evidence of haploinsufficiency effects in several genes, including TWNK (1.54 years earlier menopause, P = 1.59 × 10-6) and SOHLH2 (3.48 years earlier menopause, P = 1.03 × 10-4). Collectively, our results suggest that, for the vast majority of women, POI is not caused by autosomal dominant variants either in genes previously reported or currently evaluated in clinical diagnostic panels. Our findings, plus previous studies, suggest that most POI cases are likely oligogenic or polygenic in nature, which has important implications for future clinical genetic studies, and genetic counseling for families affected by POI

Ancient human parallel lineages within North America contributed to a coastal expansion.

Little is known regarding the first people to enter the Americas and their genetic legacy. Genomic analysis of the oldest human remains from the Americas showed a direct relationship between a Clovis-related ancestral population and all modern Central and South Americans as well as a deep split separating them from North Americans in Canada. We present 91 ancient human genomes from California and Southwestern Ontario and demonstrate the existence of two distinct ancestries in North America, which possibly split south of the ice sheets. A contribution from both of these ancestral populations is found in all modern Central and South Americans. The proportions of these two ancestries in ancient and modern populations are consistent with a coastal dispersal and multiple admixture events

Ancient human parallel lineages within North America contributed to a coastal expansion.

Little is known regarding the first people to enter the Americas and their genetic legacy. Genomic analysis of the oldest human remains from the Americas showed a direct relationship between a Clovis-related ancestral population and all modern Central and South Americans as well as a deep split separating them from North Americans in Canada. We present 91 ancient human genomes from California and Southwestern Ontario and demonstrate the existence of two distinct ancestries in North America, which possibly split south of the ice sheets. A contribution from both of these ancestral populations is found in all modern Central and South Americans. The proportions of these two ancestries in ancient and modern populations are consistent with a coastal dispersal and multiple admixture events