Population History and Gene Divergence in Native Mexicans Inferred from 76 Human Exomes.

Native American genetic variation remains underrepresented in most catalogs of human genome sequencing data. Previous genotyping efforts have revealed that Mexico's Indigenous population is highly differentiated and substructured, thus potentially harboring higher proportions of private genetic variants of functional and biomedical relevance. Here we have targeted the coding fraction of the genome and characterized its full site frequency spectrum by sequencing 76 exomes from five Indigenous populations across Mexico. Using diffusion approximations, we modeled the demographic history of Indigenous populations from Mexico with northern and southern ethnic groups splitting 7.2 KYA and subsequently diverging locally 6.5 and 5.7 KYA, respectively. Selection scans for positive selection revealed BCL2L13 and KBTBD8 genes as potential candidates for adaptive evolution in Rarámuris and Triquis, respectively. BCL2L13 is highly expressed in skeletal muscle and could be related to physical endurance, a well-known phenotype of the northern Mexico Rarámuri. The KBTBD8 gene has been associated with idiopathic short stature and we found it to be highly differentiated in Triqui, a southern Indigenous group from Oaxaca whose height is extremely low compared to other Native populations

Kidney function in tenofovir disoproxil fumarate-based oral pre-exposure prophylaxis users: a systematic review and meta-analysis of published literature and a multi-country meta-analysis of individual participant data

BACKGROUND: Previous WHO guidance on tenofovir disoproxil fumarate-based oral pre-exposure prophylaxis (PrEP) suggests measuring creatinine levels at PrEP initiation and regularly afterwards, which might represent barriers to PrEP implementation and uptake. We aimed to systematically review published literature on kidney toxicity among tenofovir disoproxil fumarate-based oral PrEP users and conducted an individual participant data meta-analysis (IPDMA) on kidney function among PrEP users in a global implementation project dataset. METHODS: In this systematic review and meta-analysis we searched PubMed up to June 30, 2021, for randomised controlled trials (RCTs) or cohort studies that reported on graded kidney-related adverse events among oral PrEP users (tenofovir disoproxil fumarate-based PrEP alone or in combination with emtricitabine or lamivudine). We extracted summary data and conducted meta-analyses with random-effects models to estimate relative risks of grade 1 and higher and grade 2 and higher kidney-related adverse events, measured by elevated serum creatinine or decline in estimated creatinine clearance or estimated glomerular filtration rate. The IPDMA included (largely unpublished) individual participant data from 17 PrEP implementation projects and two RCTs. Estimated baseline creatinine clearance and creatinine clearance change after initiation were described by age, gender, and comorbidities. We used random-effects regressions to estimate the risk in decline of creatinine clearance to less than 60 mL/min. FINDINGS: We identified 62 unique records and included 17 articles reporting on 11 RCTs with 13 523 participants in meta-analyses. PrEP use was associated with increased risk of grade 1 and higher kidney adverse events (pooled odds ratio [OR] 1·49, 95% CI 1·22-1·81; I2=25%) and grade 2 and higher events (OR 1·75, 0·68-4·49; I2=0%), although the grade 2 and higher association was not statistically significant and events were rare (13 out of 6764 in the intervention group vs six out of 6782 in the control group). The IPDMA included 18 676 individuals from 15 countries (1453 [7·8%] from RCTs) and 79 (0·42%) had a baseline estimated creatinine clearance of less than 60 mL/min (increasing proportions with increasing age). Longitudinal analyses included 14 368 PrEP users and 349 (2·43%) individuals had a decline to less than 60 mL/min creatinine clearance, with higher risks associated with increasing age and baseline creatinine clearance of 60·00-89·99 mL/min (adjusted hazard ratio [aHR] 8·49, 95% CI 6·44-11·20) and less than 60 mL/min (aHR 20·83, 12·83-33·82). INTERPRETATION: RCTs suggest that risks of kidney-related adverse events among tenofovir disoproxil fumarate-based oral PrEP users are increased but generally mild and small. Our global PrEP user analysis found varying risks by age and baseline creatinine clearance. Kidney function screening and monitoring might focus on older individuals, those with baseline creatinine clearance of less than 90 mL/min, and those with kidney-related comorbidities. Less frequent or optional screening among younger individuals without kidney-related comorbidities may reduce barriers to PrEP implementation and use. FUNDING: Unitaid, Bill & Melinda Gates Foundation, WHO

Mechanisms Underlying Interferon-γ-Induced Priming of Microglial Reactive Oxygen Species Production.

Microglial priming and enhanced reactivity to secondary insults cause substantial neuronal damage and are hallmarks of brain aging, traumatic brain injury and neurodegenerative diseases. It is, thus, of particular interest to identify mechanisms involved in microglial priming. Here, we demonstrate that priming of microglia with interferon-γ (IFN γ) substantially enhanced production of reactive oxygen species (ROS) following stimulation of microglia with ATP. Priming of microglial ROS production was substantially reduced by inhibition of p38 MAPK activity with SB203580, by increases in intracellular glutathione levels with N-Acetyl-L-cysteine, by blockade of NADPH oxidase subunit NOX2 activity with gp91ds-tat or by inhibition of nitric oxide production with L-NAME. Together, our data indicate that priming of microglial ROS production involves reduction of intracellular glutathione levels, upregulation of NADPH oxidase subunit NOX2 and increases in nitric oxide production, and suggest that these simultaneously occurring processes result in enhanced production of neurotoxic peroxynitrite. Furthermore, IFNγ-induced priming of microglial ROS production was reduced upon blockade of Kir2.1 inward rectifier K+ channels with ML133. Inhibitory effects of ML133 on microglial priming were mediated via regulation of intracellular glutathione levels and nitric oxide production. These data suggest that microglial Kir2.1 channels may represent novel therapeutic targets to inhibit excessive ROS production by primed microglia in brain pathology

The time scale of recombination rate evolution in great apes

We present three linkage-disequilibrium (LD)-based recombination maps generated using whole-genome sequence data from 10 Nigerian chimpanzees, 13 bonobos, and 15 western gorillas, collected as part of the Great Ape Genome Project (Prado-Martinez J, et al. 2013. Great ape genetic diversity and population history. Nature 499:471-475). We also identified species-specific recombination hotspots in each group using a modified LDhot framework, which greatly improves statistical power to detect hotspots at varying strengths. We show that fewer hotspots are shared among chimpanzee subspecies than within human populations, further narrowing the time scale of complete hotspot turnover. Further, using species-specific PRDM9 sequences to predict potential binding sites (PBS), we show higher predicted PRDM9 binding in recombination hotspots as compared to matched cold spot regions in multiple great ape species, including at least one chimpanzee subspecies. We found that correlations between broad-scale recombination rates decline more rapidly than nucleotide divergence between species. We also compared the skew of recombination rates at centromeres and telomeres between species and show a skew from chromosome means extending as far as 10-15Mb from chromosome ends. Further, we examined broad-scale recombination rate changes near a translocation in gorillas and found minimal differences as compared to other great ape species perhaps because the coordinates relative to the chromosome ends were unaffected. Finally, on the basis of multiple linear regression analysis, we found that various correlates of recombination rate persist throughout the African great apes including repeats, diversity, and divergence. Our study is the first to analyze within- And between-species genome-wide recombination rate variation in several close relatives

Identification of novel risk loci, causal insights, and heritable risk for Parkinson's disease: a meta-analysis of genome-wide association studies

Background Genome-wide association studies (GWAS) in Parkinson's disease have increased the scope of biological knowledge about the disease over the past decade. We aimed to use the largest aggregate of GWAS data to identify novel risk loci and gain further insight into the causes of Parkinson's disease. Methods We did a meta-analysis of 17 datasets from Parkinson's disease GWAS available from European ancestry samples to nominate novel loci for disease risk. These datasets incorporated all available data. We then used these data to estimate heritable risk and develop predictive models of this heritability. We also used large gene expression and methylation resources to examine possible functional consequences as well as tissue, cell type, and biological pathway enrichments for the identified risk factors. Additionally, we examined shared genetic risk between Parkinson's disease and other phenotypes of interest via genetic correlations followed by Mendelian randomisation. Findings Between Oct 1, 2017, and Aug 9, 2018, we analysed 7·8 million single nucleotide polymorphisms in 37 688 cases, 18 618 UK Biobank proxy-cases (ie, individuals who do not have Parkinson's disease but have a first degree relative that does), and 1·4 million controls. We identified 90 independent genome-wide significant risk signals across 78 genomic regions, including 38 novel independent risk signals in 37 loci. These 90 variants explained 16–36% of the heritable risk of Parkinson's disease depending on prevalence. Integrating methylation and expression data within a Mendelian randomisation framework identified putatively associated genes at 70 risk signals underlying GWAS loci for follow-up functional studies. Tissue-specific expression enrichment analyses suggested Parkinson's disease loci were heavily brain-enriched, with specific neuronal cell types being implicated from single cell data. We found significant genetic correlations with brain volumes (false discovery rate-adjusted p=0·0035 for intracranial volume, p=0·024 for putamen volume), smoking status (p=0·024), and educational attainment (p=0·038). Mendelian randomisation between cognitive performance and Parkinson's disease risk showed a robust association (p=8·00 × 10−7). Interpretation These data provide the most comprehensive survey of genetic risk within Parkinson's disease to date, to the best of our knowledge, by revealing many additional Parkinson's disease risk loci, providing a biological context for these risk factors, and showing that a considerable genetic component of this disease remains unidentified. These associations derived from European ancestry datasets will need to be followed-up with more diverse data. Funding The National Institute on Aging at the National Institutes of Health (USA), The Michael J Fox Foundation, and The Parkinson's Foundation (see appendix for full list of funding sources)

Genome-wide analysis of 53,400 people with irritable bowel syndrome highlights shared genetic pathways with mood and anxiety disorders

Funder: Kennedy Trust Rheumatology Research Prize StudentshipFunder: DFG Cluster of Excellence “Precision Medicine in Chronic In-flammation” (PMI; ID: EXC2167)Funder: EC | EC Seventh Framework Programm | FP7 Ideas: European Research Council (FP7-IDEAS-ERC - Specific Programme: “Ideas” Implementing the Seventh Framework Programme of the European Community for Research, Technological Development and Demonstration Activities (2007 to 2013)); doi: https://doi.org/10.13039/100011199; Grant(s): 715772Funder: NWO-VIDI grant 016.178.056, the Netherlands Heart Foundation CVON grant 2018-27, and NWO Gravitation grant ExposomeNLFunder: Li Ka Shing Foundation (Li Ka Shing Foundation Limited); doi: https://doi.org/10.13039/100007421Abstract: Irritable bowel syndrome (IBS) results from disordered brain–gut interactions. Identifying susceptibility genes could highlight the underlying pathophysiological mechanisms. We designed a digestive health questionnaire for UK Biobank and combined identified cases with IBS with independent cohorts. We conducted a genome-wide association study with 53,400 cases and 433,201 controls and replicated significant associations in a 23andMe panel (205,252 cases and 1,384,055 controls). Our study identified and confirmed six genetic susceptibility loci for IBS. Implicated genes included NCAM1, CADM2, PHF2/FAM120A, DOCK9, CKAP2/TPTE2P3 and BAG6. The first four are associated with mood and anxiety disorders, expressed in the nervous system, or both. Mirroring this, we also found strong genome-wide correlation between the risk of IBS and anxiety, neuroticism and depression (rg > 0.5). Additional analyses suggested this arises due to shared pathogenic pathways rather than, for example, anxiety causing abdominal symptoms. Implicated mechanisms require further exploration to help understand the altered brain–gut interactions underlying IBS

Population genetic analysis of the DARC locus (Duffy) reveals adaptation from standing variation associated with malaria resistance in humans.

The human DARC (Duffy antigen receptor for chemokines) gene encodes a membrane-bound chemokine receptor crucial for the infection of red blood cells by Plasmodium vivax, a major causative agent of malaria. Of the three major allelic classes segregating in human populations, the FY*O allele has been shown to protect against P. vivax infection and is at near fixation in sub-Saharan Africa, while FY*B and FY*A are common in Europe and Asia, respectively. Due to the combination of strong geographic differentiation and association with malaria resistance, DARC is considered a canonical example of positive selection in humans. Despite this, details of the timing and mode of selection at DARC remain poorly understood. Here, we use sequencing data from over 1,000 individuals in twenty-one human populations, as well as ancient human genomes, to perform a fine-scale investigation of the evolutionary history of DARC. We estimate the time to most recent common ancestor (TMRCA) of the most common FY*O haplotype to be 42 kya (95% CI: 34-49 kya). We infer the FY*O null mutation swept to fixation in Africa from standing variation with very low initial frequency (0.1%) and a selection coefficient of 0.043 (95% CI:0.011-0.18), which is among the strongest estimated in the human genome. We estimate the TMRCA of the FY*A mutation in non-Africans to be 57 kya (95% CI: 48-65 kya) and infer that, prior to the sweep of FY*O, all three alleles were segregating in Africa, as highly diverged populations from Asia and ≠Khomani San hunter-gatherers share the same FY*A haplotypes. We test multiple models of admixture that may account for this observation and reject recent Asian or European admixture as the cause

Population genetic analysis of the DARC locus (Duffy) reveals adaptation from standing variation associated with malaria resistance in humans

<div><p>The human DARC (Duffy antigen receptor for chemokines) gene encodes a membrane-bound chemokine receptor crucial for the infection of red blood cells by <i>Plasmodium vivax</i>, a major causative agent of malaria. Of the three major allelic classes segregating in human populations, the FY*O allele has been shown to protect against <i>P. vivax</i> infection and is at near fixation in sub-Saharan Africa, while FY*B and FY*A are common in Europe and Asia, respectively. Due to the combination of strong geographic differentiation and association with malaria resistance, DARC is considered a canonical example of positive selection in humans. Despite this, details of the timing and mode of selection at DARC remain poorly understood. Here, we use sequencing data from over 1,000 individuals in twenty-one human populations, as well as ancient human genomes, to perform a fine-scale investigation of the evolutionary history of DARC. We estimate the time to most recent common ancestor (T</p><p>_MRCA</p>) of the most common FY*O haplotype to be 42 kya (95% CI: 34–49 kya). We infer the FY*O null mutation swept to fixation in Africa from standing variation with very low initial frequency (0.1%) and a selection coefficient of 0.043 (95% CI:0.011–0.18), which is among the strongest estimated in the human genome. We estimate the T<p>_MRCA</p> of the FY*A mutation in non-Africans to be 57 kya (95% CI: 48–65 kya) and infer that, prior to the sweep of FY*O, all three alleles were segregating in Africa, as highly diverged populations from Asia and ≠Khomani San hunter-gatherers share the same FY*A haplotypes. We test multiple models of admixture that may account for this observation and reject recent Asian or European admixture as the cause.<p></p></div