A multi-cohort genome-wide association study in African ancestry individuals reveals risk loci for primary open-angle glaucoma.

Primary open-angle glaucoma (POAG), the leading cause of irreversible blindness worldwide, disproportionately affects individuals of African ancestry. We conducted a genome-wide association study (GWAS) for POAG in 11,275 individuals of African ancestry (6,003 cases; 5,272 controls). We detected 46 risk loci associated with POAG at genome-wide significance. Replication and post-GWAS analyses, including functionally informed fine-mapping, multiple trait co-localization, and in silico validation, implicated two previously undescribed variants (rs1666698 mapping to DBF4P2; rs34957764 mapping to ROCK1P1) and one previously associated variant (rs11824032 mapping to ARHGEF12) as likely causal. For individuals of African ancestry, a polygenic risk score (PRS) for POAG from our mega-analysis (African ancestry individuals) outperformed a PRS from summary statistics of a much larger GWAS derived from European ancestry individuals. This study quantifies the genetic architecture similarities and differences between African and non-African ancestry populations for this blinding disease.

The Missing Diversity in Human Genetic Studies.

The majority of studies of genetic association with disease have been performed in Europeans. This European bias has important implications for risk prediction of diseases across global populations. In this commentary, we justify the need to study more diverse populations using both empirical examples and theoretical reasoning.

The African Society of Human Genetics successfully launches global data science workshops.

To accelerate the impact of African genomics on human health, data science skills and awareness of Africa's rich genetic diversity must be strengthened globally. We describe the first African genomics data science workshop, implemented by the African Society of Human Genetics (AfSHG) and international partners, providing a framework for future workshops.

Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets.

Leveraging linkage disequilibrium (LD) patterns as representative of population substructure enables the discovery of additive association signals in genome-wide association studies (GWASs). Standard GWASs are well-powered to interrogate additive models; however, new approaches are required for invesigating other modes of inheritance such as dominance and epistasis. Epistasis, or non-additive interaction between genes, exists across the genome but often goes undetected because of a lack of statistical power. Furthermore, the adoption of LD pruning as customary in standard GWASs excludes detection of sites that are in LD but might underlie the genetic architecture of complex traits. We hypothesize that uncovering long-range interactions between loci with strong LD due to epistatic selection can elucidate genetic mechanisms underlying common diseases. To investigate this hypothesis, we tested for associations between 23 common diseases and 5,625,845 epistatic SNP-SNP pairs (determined by Ohta's D statistics) in long-range LD (>0.25 cM). Across five disease phenotypes, we identified one significant and four near-significant associations that replicated in two large genotype-phenotype datasets (UK Biobank and eMERGE). The genes that were most likely involved in the replicated associations were (1) members of highly conserved gene families with complex roles in multiple pathways, (2) essential genes, and/or (3) genes that were associated in the literature with complex traits that display variable expressivity. These results support the highly pleiotropic and conserved nature of variants in long-range LD under epistatic selection. Our work supports the hypothesis that epistatic interactions regulate diverse clinical mechanisms and might especially be driving factors in conditions with a wide range of phenotypic outcomes.

Tuberculosis severity associates with variants and eQTLs related to vascular biology and infection-induced inflammation.

BACKGROUND: Tuberculosis (TB) remains a major public health problem globally, even compared to COVID-19. Genome-wide studies have failed to discover genes that explain a large proportion of genetic risk for adult pulmonary TB, and even fewer have examined genetic factors underlying TB severity, an intermediate trait impacting disease experience, quality of life, and risk of mortality. No prior severity analyses used a genome-wide approach. METHODS AND FINDINGS: As part of our ongoing household contact study in Kampala, Uganda, we conducted a genome-wide association study (GWAS) of TB severity measured by TBScore, in two independent cohorts of culture-confirmed adult TB cases (n = 149 and n = 179). We identified 3 SNPs (P<1.0 x 10-7) including one on chromosome 5, rs1848553, that was GWAS significant (meta-analysis p = 2.97x10-8). All three SNPs are in introns of RGS7BP and have effect sizes corresponding to clinically meaningful reductions in disease severity. RGS7BP is highly expressed in blood vessels and plays a role in infectious disease pathogenesis. Other genes with suggestive associations defined gene sets involved in platelet homeostasis and transport of organic anions. To explore functional implications of the TB severity-associated variants, we conducted eQTL analyses using expression data from Mtb-stimulated monocyte-derived macrophages. A single variant (rs2976562) associated with monocyte SLA expression (p = 0.03) and subsequent analyses indicated that SLA downregulation following MTB stimulation associated with increased TB severity. Src Like Adaptor (SLAP-1), encoded by SLA, is highly expressed in immune cells and negatively regulates T cell receptor signaling, providing a potential mechanistic link to TB severity. CONCLUSIONS: These analyses reveal new insights into the genetics of TB severity with regulation of platelet homeostasis and vascular biology being central to consequences for active TB patients. This analysis also reveals genes that regulate inflammation can lead to differences in severity. Our findings provide an important step in improving TB patient outcomes.

Impact of natural selection on global patterns of genetic variation and association with clinical phenotypes at genes involved in SARS-CoV-2 infection.

Human genomic diversity has been shaped by both ancient and ongoing challenges from viruses. The current coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a devastating impact on population health. However, genetic diversity and evolutionary forces impacting host genes related to SARS-CoV-2 infection are not well understood. We investigated global patterns of genetic variation and signatures of natural selection at host genes relevant to SARS-CoV-2 infection (angiotensin converting enzyme 2 [ACE2], transmembrane protease serine 2 [TMPRSS2], dipeptidyl peptidase 4 [DPP4], and lymphocyte antigen 6 complex locus E [LY6E]). We analyzed data from 2,012 ethnically diverse Africans and 15,977 individuals of European and African ancestry with electronic health records and integrated with global data from the 1000 Genomes Project. At ACE2, we identified 41 nonsynonymous variants that were rare in most populations, several of which impact protein function. However, three nonsynonymous variants (rs138390800, rs147311723, and rs145437639) were common among central African hunter-gatherers from Cameroon (minor allele frequency 0.083 to 0.164) and are on haplotypes that exhibit signatures of positive selection. We identify signatures of selection impacting variation at regulatory regions influencing ACE2 expression in multiple African populations. At TMPRSS2, we identified 13 amino acid changes that are adaptive and specific to the human lineage compared with the chimpanzee genome. Genetic variants that are targets of natural selection are associated with clinical phenotypes common in patients with COVID-19. Our study provides insights into global variation at host genes related to SARS-CoV-2 infection, which have been shaped by natural selection in some populations, possibly due to prior viral infections.

Hemizygosity can reveal variant pathogenicity on the X-chromosome.

Pathogenic variants on the X-chromosome can have more severe consequences for hemizygous males, while heterozygote females can avoid severe consequences due to diploidy and the capacity for nonrandom expression. Thus, when an allele is more common in females this could indicate that it increases the probability of early death in the male hemizygous state, which can be considered a measure of pathogenicity. Importantly, large-scale genomic data now makes it possible to compare allele proportions between the sexes. To discover pathogenic variants on the X-chromosome, we analyzed exome data from 125,748 ancestrally diverse participants in the Genome Aggregation Database (gnomAD). After filtering out duplicates and extremely rare variants, 44,606 of the original 348,221 remained for analysis. We divided the proportion of variant alleles in females by the proportion in males for all variant sites, and then placed each variant into one of three a priori categories: (1) Reference (Primarily synonymous and intronic), (2) Unlikely-to-be-tolerated (Primarily missense), and (3) Least-likely-to-be-tolerated (Primarily frameshift). To assess the impact of ploidy, we compared the distribution of these ratios between pseudoautosomal and non-pseudoautosomal regions. In the non-pseudoautosomal regions, mean female-to-male ratios were lowest among Reference (2.40), greater for Unlikely-to-be-tolerated (2.77) and highest for Least-likely-to-be-tolerated (3.28) variants. Corresponding ratios were lower in the pseudoautosomal regions (1.52, 1.57, and 1.68, respectively), with the most extreme ratio being just below 11. Because pathogenic effects in the pseudoautosomal regions should not drive ratio increases, this maximum ratio provides an upper bound for baseline noise. In the non-pseudoautosomal regions, 319 variants had a ratio over 11. In sum, we identified a measure with a dataset specific threshold for identifying pathogenicity in non-pseudoautosomal X-chromosome variants: the female-to-male allele proportion ratio.

Late-onset neonatal sepsis: genetic differences by sex and involvement of the NOTCH pathway.

BACKGROUND: Late-Onset Neonatal Sepsis (LOS) is a rare condition, involving widespread infection, immune disruption, organ dysfunction, and often death. Because exposure to pathogens is not completely preventable, identifying susceptibility factors is critical to characterizing the pathophysiology and developing interventions. Prior studies demonstrated both genetics and infant sex influence susceptibility. Our study was designed to identify LOS associated genetic variants. METHODS: We performed an exploratory genome wide association study (GWAS) with 224 LOS cases and 273 controls from six European countries. LOS was defined as sepsis presenting from 3 to 90 days of age; diagnosis was established by clinical criteria consensus guidelines. We tested for association with both autosomal and X-chromosome variants in the total sample and in sex-stratified analyses. RESULTS: In total, 71 SNPs associated with neonatal sepsis at p < 1 × 10-4 in at least one analysis. Most importantly, sex-stratified analyses revealed associations with multiple SNPs (28 in males and 16 in females), but no variants from single-sex analyses associated with sepsis in the other sex. Pathway analyses showed NOTCH signaling is over-represented among genes linked to these SNPS. CONCLUSION: Our results indicate genetic susceptibility to LOS is sexually dimorphic and corroborate that NOTCH signaling plays a role in determining risk. IMPACT: Genes associate with late onset neonatal sepsis. Notch pathway genes are overrepresented in associations with sepsis. Genes associating with sepsis do not overlap between males and females. Sexual dimorphism can lead to sex specific treatment of sepsis.

Interaction between host genes and Mycobacterium tuberculosis lineage can affect tuberculosis severity: Evidence for coevolution?

Genetic studies of both the human host and Mycobacterium tuberculosis (MTB) demonstrate independent association with tuberculosis (TB) risk. However, neither explains a large portion of disease risk or severity. Based on studies in other infectious diseases and animal models of TB, we hypothesized that the genomes of the two interact to modulate risk of developing active TB or increasing the severity of disease, when present. We examined this hypothesis in our TB household contact study in Kampala, Uganda, in which there were 3 MTB lineages of which L4-Ugandan (L4.6) is the most recent. TB severity, measured using the Bandim TBscore, was modeled as a function of host SNP genotype, MTB lineage, and their interaction, within two independent cohorts of TB cases, N = 113 and 121. No association was found between lineage and severity, but association between multiple polymorphisms in IL12B and TBscore was replicated in two independent cohorts (most significant rs3212227, combined p = 0.0006), supporting previous associations of IL12B with TB susceptibility. We also observed significant interaction between a single nucleotide polymorphism (SNP) in SLC11A1 and the L4-Ugandan lineage in both cohorts (rs17235409, meta p = 0.0002). Interestingly, the presence of the L4-Uganda lineage in the presence of the ancestral human allele associated with more severe disease. These findings demonstrate that IL12B is associated with severity of TB in addition to susceptibility, and that the association between TB severity and human genetics can be due to an interaction between genes in the two species, consistent with host-pathogen coevolution in TB.

Genetics and geography of leukocyte telomere length in sub-Saharan Africans.

Leukocyte telomere length (LTL) might be causal in cardiovascular disease and major cancers. To elucidate the roles of genetics and geography in LTL variability across humans, we compared LTL measured in 1295 sub-Saharan Africans (SSAs) with 559 African-Americans (AAms) and 2464 European-Americans (EAms). LTL differed significantly across SSAs (P = 0.003), with the San from Botswana (with the oldest genomic ancestry) having the longest LTL and populations from Ethiopia having the shortest LTL. SSAs had significantly longer LTL than AAms [P = 6.5(e-16)] whose LTL was significantly longer than EAms [P = 2.5(e-7)]. Genetic variation in SSAs explained 52% of LTL variance versus 27% in AAms and 34% in EAms. Adjustment for genetic variation removed the LTL differences among SSAs. LTL genetic variation among SSAs, with the longest LTL in the San, supports the hypothesis that longer LTL was ancestral in humans. Identifying factors driving LTL variation in Africa may have important ramifications for LTL-associated diseases.

Estimating prevalence of human traits among populations from polygenic risk scores.

The genetic basis of phenotypic variation across populations has not been well explained for most traits. Several factors may cause disparities, from variation in environments to divergent population genetic structure. We hypothesized that a population-level polygenic risk score (PRS) can explain phenotypic variation among geographic populations based solely on risk allele frequencies. We applied a population-specific PRS (psPRS) to 26 populations from the 1000 Genomes to four phenotypes: lactase persistence (LP), melanoma, multiple sclerosis (MS) and height. Our models assumed additive genetic architecture among the polymorphisms in the psPRSs, as is convention. Linear psPRSs explained a significant proportion of trait variance ranging from 0.32 for height in men to 0.88 for melanoma. The best models for LP and height were linear, while those for melanoma and MS were nonlinear. As not all variants in a PRS may confer similar, or even any, risk among diverse populations, we also filtered out SNPs to assess whether variance explained was improved using psPRSs with fewer SNPs. Variance explained usually improved with fewer SNPs in the psPRS and was as high as 0.99 for height in men using only 548 of the initial 4208 SNPs. That reducing SNPs improves psPRSs performance may indicate that missing heritability is partially due to complex architecture that does not mandate additivity, undiscovered variants or spurious associations in the databases. We demonstrated that PRS-based analyses can be used across diverse populations and phenotypes for population prediction and that these comparisons can identify the universal risk variants.

Methylome-wide Analysis Reveals Epigenetic Marks Associated With Resistance to Tuberculosis in Human Immunodeficiency Virus-Infected Individuals From East Africa.

BACKGROUND: Tuberculosis (TB) is the most deadly infectious disease globally and is highly prevalent in the developing world. For individuals infected with both Mycobacterium tuberculosis (Mtb) and human immunodeficiency virus (HIV), the risk of active TB is 10% or more annually. Previously, we identified in a genome-wide association study (GWAS) a region on chromosome 5 associated with resistance to TB, which included epigenetic marks that could influence gene regulation. We hypothesized that HIV-infected individuals exposed to Mtb who remain disease free carry epigenetic changes that strongly protect them from active TB. METHODS: We conducted a methylome-wide study in HIV-infected, TB-exposed cohorts from Uganda and Tanzania and integrated data from our GWAS. RESULTS: We identified 3 regions of interest that included markers that were differentially methylated between TB cases and controls with latent TB infection: chromosome 1 (RNF220, P = 4 × 10-5), chromosome 2 (between COPS8 and COL6A3, P = 2.7 × 10-5), and chromosome 5 (CEP72, P = 1.3 × 10-5). These methylation results co-localized with associated single-nucleotide polymorphisms (SNPs), methylation QTLs, and methylation × SNP interaction effects. These markers were in regions with regulatory markers for cells involved in TB immunity and/or lung. CONCLUSIONS: Epigenetic regulation is a potential biologic factor underlying resistance to TB in immunocompromised individuals that can act in conjunction with genetic variants.

Helicobacter pylori infection causes both protective and deleterious effects in human health and disease.

Infection with Helicobacter pylori (H. pylori) is necessary but not sufficient for the development of gastric cancer, the third leading cause of cancer death globally. H. pylori infection affects over half of people globally; however, it does not affect populations uniformly. H. pylori infection rates are declining in western industrialized countries but are plateauing in developing and newly industrialized countries where gastric cancer is most prevalent. Despite H. pylori infection being the primary causative agent for gastric cancer, H. pylori infection can also cause other effects, detrimental or beneficial, throughout an individual's life, with the beneficial effects often being seen in childhood and the deleterious effects in adulthood. H. pylori is an ancient bacterium and its likelihood of affecting disease or health is dependent on both human and bacterial genetics that have co-evolved over millennia. In this review, we focus on the impact of infection and its genetic bases in different populations and diseases throughout an individual's lifespan, highlighting the benefits of individualized treatment and argue that universal eradication of H. pylori in its host may cause more harm than good for those infected with H. pylori.

Correction: Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

[This corrects the article DOI: 10.1371/journal.pgen.1007394.].

Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

Preterm birth is a leading cause of morbidity and mortality in infants. Genetic and environmental factors play a role in the susceptibility to preterm birth, but despite many investigations, the genetic basis for preterm birth remain largely unknown. Our objective was to identify rare, possibly damaging, nucleotide variants in mothers from families with recurrent spontaneous preterm births (SPTB). DNA samples from 17 Finnish mothers who delivered at least one infant preterm were subjected to whole exome sequencing. All mothers were of northern Finnish origin and were from seven multiplex families. Additional replication samples of European origin consisted of 93 Danish sister pairs (and two sister triads), all with a history of a preterm delivery. Rare exonic variants (frequency <1%) were analyzed to identify genes and pathways likely to affect SPTB susceptibility. We identified rare, possibly damaging, variants in genes that were common to multiple affected individuals. The glucocorticoid receptor signaling pathway was the most significant (p<1.7e-8) with genes containing these variants in a subgroup of ten Finnish mothers, each having had 2-4 SPTBs. This pathway was replicated among the Danish sister pairs. A gene in this pathway, heat shock protein family A (Hsp70) member 1 like (HSPA1L), contains two likely damaging missense alleles that were found in four different Finnish families. One of the variants (rs34620296) had a higher frequency in cases compared to controls (0.0025 vs. 0.0010, p = 0.002) in a large preterm birth genome-wide association study (GWAS) consisting of mothers of general European ancestry. Sister pairs in replication samples also shared rare, likely damaging HSPA1L variants. Furthermore, in silico analysis predicted an additional phosphorylation site generated by rs34620296 that could potentially affect chaperone activity or HSPA1L protein stability. Finally, in vitro functional experiment showed a link between HSPA1L activity and decidualization. In conclusion, rare, likely damaging, variants in HSPA1L were observed in multiple families with recurrent SPTB.

A comparison of two workflows for regulome and transcriptome-based prioritization of genetic variants associated with myocardial mass.

A typical task arising from main effect analyses in a Genome Wide Association Study (GWAS) is to identify single nucleotide polymorphisms (SNPs), in linkage disequilibrium with the observed signals, that are likely causal variants and the affected genes. The affected genes may not be those closest to associating SNPs. Functional genomics data from relevant tissues are believed to be helpful in selecting likely causal SNPs and interpreting implicated biological mechanisms, ultimately facilitating prevention and treatment in the case of a disease trait. These data are typically used post GWAS analyses to fine-map the statistically significant signals identified agnostically by testing all SNPs and applying a multiple testing correction. The number of tested SNPs is typically in the millions, so the multiple testing burden is high. Motivated by this, in this study we investigated an alternative workflow, which consists in utilizing the available functional genomics data as a first step to reduce the number of SNPs tested for association. We analyzed GWAS on electrocardiographic QRS duration using these two workflows. The alternative workflow identified more SNPs, including some residing in loci not discovered with the typical workflow. Moreover, the latter are corroborated by other reports on QRS duration. This indicates the potential value of incorporating functional genomics information at the onset in GWAS analyses.

Dissecting maternal and fetal genetic effects underlying the associations between maternal phenotypes, birth outcomes, and adult phenotypes: A mendelian-randomization and haplotype-based genetic score analysis in 10,734 mother-infant pairs.

BACKGROUND: Many maternal traits are associated with a neonate's gestational duration, birth weight, and birth length. These birth outcomes are subsequently associated with late-onset health conditions. The causal mechanisms and the relative contributions of maternal and fetal genetic effects behind these observed associations are unresolved. METHODS AND FINDINGS: Based on 10,734 mother-infant duos of European ancestry from the UK, Northern Europe, Australia, and North America, we constructed haplotype genetic scores using single-nucleotide polymorphisms (SNPs) known to be associated with adult height, body mass index (BMI), blood pressure (BP), fasting plasma glucose (FPG), and type 2 diabetes (T2D). Using these scores as genetic instruments, we estimated the maternal and fetal genetic effects underlying the observed associations between maternal phenotypes and pregnancy outcomes. We also used infant-specific birth weight genetic scores as instrument and examined the effects of fetal growth on pregnancy outcomes, maternal BP, and glucose levels during pregnancy. The maternal nontransmitted haplotype score for height was significantly associated with gestational duration (p = 2.2 × 10-4). Both maternal and paternal transmitted height haplotype scores were highly significantly associated with birth weight and length (p < 1 × 10-17). The maternal transmitted BMI scores were associated with birth weight with a significant maternal effect (p = 1.6 × 10-4). Both maternal and paternal transmitted BP scores were negatively associated with birth weight with a significant fetal effect (p = 9.4 × 10-3), whereas BP alleles were significantly associated with gestational duration and preterm birth through maternal effects (p = 3.3 × 10-2 and p = 4.5 × 10-3, respectively). The nontransmitted haplotype score for FPG was strongly associated with birth weight (p = 4.7 × 10-6); however, the glucose-increasing alleles in the fetus were associated with reduced birth weight through a fetal effect (p = 2.2 × 10-3). The haplotype scores for T2D were associated with birth weight in a similar way but with a weaker maternal effect (p = 6.4 × 10-3) and a stronger fetal effect (p = 1.3 × 10-5). The paternal transmitted birth weight score was significantly associated with reduced gestational duration (p = 1.8 × 10-4) and increased maternal systolic BP during pregnancy (p = 2.2 × 10-2). The major limitations of the study include missing and heterogenous phenotype data in some data sets and different instrumental strength of genetic scores for different phenotypic traits. CONCLUSIONS: We found that both maternal height and fetal growth are important factors in shaping the duration of gestation: genetically elevated maternal height is associated with longer gestational duration, whereas alleles that increase fetal growth are associated with shorter gestational duration. Fetal growth is influenced by both maternal and fetal effects and can reciprocally influence maternal phenotypes: taller maternal stature, higher maternal BMI, and higher maternal blood glucose are associated with larger birth size through maternal effects; in the fetus, the height- and metabolic-risk-increasing alleles are associated with increased and decreased birth size, respectively; alleles raising birth weight in the fetus are associated with shorter gestational duration and higher maternal BP. These maternal and fetal genetic effects may explain the observed associations between the studied maternal phenotypes and birth outcomes, as well as the life-course associations between these birth outcomes and adult phenotypes.

Association of preeclampsia with infant APOL1 genotype in African Americans.

BACKGROUND: Black women in the United States and Africa are at an increased risk for preeclampsia. Allelic variants in the gene for apolipoprotein LI, APOL1, are found only in populations of African ancestry, and have been shown to contribute significant risk for kidney disease. Recent studies suggest these APOL1 variants also may contribute risk for preeclampsia. METHODS: The association of preeclampsia with carriage of APOL1 risk alleles was evaluated in a case-control study of deliveries from black women at a single center in Cleveland, Ohio that included gross and histopathologic evaluations of placental tissues (395 cases and 282 controls). Using logistic regression models, associations between fetal APOL1 genotype and preeclampsia were evaluated using several case definitions based on prematurity and severity of preeclampsia, with uncomplicated term pregnancies as controls. Associations between APOL1 genotype and pathological features were also examined. RESULTS: The infant APOL1 genotype was significantly associated with preeclampsia in a dominant inheritance pattern with odds ratio of 1.41 (P=0.029, 95% CI 1.037, 1.926). Stratifying preeclampsia cases by preterm birth, significant associations were detected for both recessive (O.R.=1.70, P=0.038) and additive (O.R.=1.33, P=0.028) inheritance patterns. APOL1 genotype, however, was not significantly associated with pathological changes or other perinatal observations. CONCLUSIONS: Preeclampsia appears to be another disease associated with APOL1 variants, however, further studies are needed to increase confidence in the mode of inheritance. By understanding the association of APOL1 variants with preeclampsia, genetic screening tests for APOL1 may be useful to predict at-risk pregnancies and targeted interventions may be developed to improve pregnancy outcomes.