Untargeted metabolomic profiling reveals molecular signatures associated with type 2 diabetes in Nigerians.

BACKGROUND: Type 2 diabetes (T2D) has reached epidemic proportions globally, including in Africa. However, molecular studies to understand the pathophysiology of T2D remain scarce outside Europe and North America. The aims of this study are to use an untargeted metabolomics approach to identify: (a) metabolites that are differentially expressed between individuals with and without T2D and (b) a metabolic signature associated with T2D in a population of Sub-Saharan Africa (SSA). METHODS: A total of 580 adult Nigerians from the Africa America Diabetes Mellitus (AADM) study were studied. The discovery study included 310 individuals (210 without T2D, 100 with T2D). Metabolites in plasma were assessed by reverse phase, ultra-performance liquid chromatography and mass spectrometry (RP)/UPLC-MS/MS methods on the Metabolon Platform. Welch's two-sample t-test was used to identify differentially expressed metabolites (DEMs), followed by the construction of a biomarker panel using a random forest (RF) algorithm. The biomarker panel was evaluated in a replication sample of 270 individuals (110 without T2D and 160 with T2D) from the same study. RESULTS: Untargeted metabolomic analyses revealed 280 DEMs between individuals with and without T2D. The DEMs predominantly belonged to the lipid (51%, 142/280), amino acid (21%, 59/280), xenobiotics (13%, 35/280), carbohydrate (4%, 10/280) and nucleotide (4%, 10/280) super pathways. At the sub-pathway level, glycolysis, free fatty acid, bile metabolism, and branched chain amino acid catabolism were altered in T2D individuals. A 10-metabolite biomarker panel including glucose, gluconate, mannose, mannonate, 1,5-anhydroglucitol, fructose, fructosyl-lysine, 1-carboxylethylleucine, metformin, and methyl-glucopyranoside predicted T2D with an area under the curve (AUC) of 0.924 (95% CI: 0.845-0.966) and a predicted accuracy of 89.3%. The panel was validated with a similar AUC (0.935, 95% CI 0.906-0.958) in the replication cohort. The 10 metabolites in the biomarker panel correlated significantly with several T2D-related glycemic indices, including Hba1C, insulin resistance (HOMA-IR), and diabetes duration. CONCLUSIONS: We demonstrate that metabolomic dysregulation associated with T2D in Nigerians affects multiple processes, including glycolysis, free fatty acid and bile metabolism, and branched chain amino acid catabolism. Our study replicated previous findings in other populations and identified a metabolic signature that could be used as a biomarker panel of T2D risk and glycemic control thus enhancing our knowledge of molecular pathophysiologic changes in T2D. The metabolomics dataset generated in this study represents an invaluable addition to publicly available multi-omics data on understudied African ancestry populations.

Random forest classifiers trained on simulated data enable accurate short read-based genotyping of structural variants in the alpha globin region at Chr16p13.3.

In regions where reads don't align well to a reference, it is generally difficult to characterize structural variation using short read sequencing. Here, we utilize machine learning classifiers and short sequence reads to genotype structural variants in the alpha globin locus on chromosome 16, a medically-relevant region that is challenging to genotype in individuals. Using models trained only with simulated data, we accurately genotype two hard-to-distinguish deletions in two separate human cohorts. Furthermore, population allele frequencies produced by our methods across a wide set of ancestries agree more closely with previously-determined frequencies than those obtained using currently available genotyping software.

Unappreciated subcontinental admixture in Europeans and European Americans and implications for genetic epidemiology studies.

European-ancestry populations are recognized as stratified but not as admixed, implying that residual confounding by locus-specific ancestry can affect studies of association, polygenic adaptation, and polygenic risk scores. We integrate individual-level genome-wide data from ~19,000 European-ancestry individuals across 79 European populations and five European American cohorts. We generate a new reference panel that captures ancestral diversity missed by both the 1000 Genomes and Human Genome Diversity Projects. Both Europeans and European Americans are admixed at the subcontinental level, with admixture dates differing among subgroups of European Americans. After adjustment for both genome-wide and locus-specific ancestry, associations between a highly differentiated variant in LCT (rs4988235) and height or LDL-cholesterol were confirmed to be false positives whereas the association between LCT and body mass index was genuine. We provide formal evidence of subcontinental admixture in individuals with European ancestry, which, if not properly accounted for, can produce spurious results in genetic epidemiology studies.

Universal genome-wide association studies: Powerful joint ancestry and association testing.

The vast majority of human populations and individuals have mixed ancestry. Consequently, adjustment for locus-specific ancestry is essential for genetic association studies. To empower association studies for all populations, it is necessary to integrate effects of locus-specific ancestry and genotype. We developed a joint test of ancestry and association that can be performed with summary statistics, is independent of study design, can take advantage of locus-specific ancestry effects to boost power in association testing, and can utilize association effects to fine map admixture peaks. We illustrate the test using the association between serum triglycerides and LPL. By combining data from African Americans, European Americans, and West Africans, we identify three conditionally independent variants with varying amounts of ancestrally differentiated allele frequencies. Using out-of-sample data, we demonstrate improved prediction achievable by accounting for multiple causal variants and locus-specific ancestry effects at a single locus.

Overview of Admixture Mapping.

Admixture mapping is a powerful method of gene mapping for diseases or traits that show differential risk by ancestry. Admixture mapping has been applied most often to Americans who trace ancestry to various combinations of Native Americans, Europeans, and West Africans. Recent developments in admixture mapping include improvements in the methods and reference data needed to make inferences about ancestry as well as extensions of the mapping approach in the framework of linear mixed models. This overview outlines the key concepts of admixture mapping. It describes approaches for inferring local ancestry and provides strategies for performing admixture mapping depending on the study design. Finally, linkage analysis, association analysis, and admixture mapping are compared, with an emphasis on integrating admixture mapping and association testing. Published 2023. This article is a U.S. Government work and is in the public domain in the USA.

Discovery of a rapidly evolving yeast defense factor, KTD1, against the secreted killer toxin K28.

Secreted protein toxins are widely used weapons in conflicts between organisms. Elucidating how organisms genetically adapt to defend themselves against these toxins is fundamental to understanding the coevolutionary dynamics of competing organisms. Within yeast communities, "killer" toxins are secreted to kill nearby sensitive yeast, providing a fitness advantage in competitive growth environments. Natural yeast isolates vary in their sensitivity to these toxins, but to date, no polymorphic genetic factors contributing to defense have been identified. We investigated the variation in resistance to the killer toxin K28 across diverse natural isolates of the Saccharomyces cerevisiae population. Using large-scale linkage mapping, we discovered a novel defense factor, which we named KTD1. We identified many KTD1 alleles, which provided different levels of K28 resistance. KTD1 is a member of the DUP240 gene family of unknown function, which is rapidly evolving in a region spanning its two encoded transmembrane helices. We found that this domain is critical to KTD1's protective ability. Our findings implicate KTD1 as a key polymorphic factor in the defense against K28 toxin.

Additive genetic effect of GCKR, G6PC2, and SLC30A8 variants on fasting glucose levels and risk of type 2 diabetes.

Impaired glucose tolerance is a major risk factor for type 2 diabetes (T2D) and several cardiometabolic disorders. To identify genetic loci underlying fasting glucose levels, we conducted an analysis of 9,232 individuals of European ancestry who at enrollment were either nondiabetic or had untreated type 2 diabetes. Multivariable linear mixed models were used to test for associations between fasting glucose and 7.9 million SNPs, with adjustment for age, body mass index (BMI), sex, significant principal components of the genotypes, and cryptic relatedness. Three previously discovered loci were genome-wide significant, with the lead SNPs being rs1260326, a missense variant in GCKR (p = 1.06×10-8); rs560887, an intronic variant in G6PC2 (p = 3.39×10-11); and rs13266634, a missense variant in SLC30A8 (p = 4.28×10-10). Fine mapping, genome-wide conditional analysis, and functional annotation indicated that the three loci were independently associated with fasting glucose. Each copy of an alternate allele at any of these three SNPs was associated with a reduction of 0.012 mmol/L in fasting glucose levels (p = 8.0×10-28), and this association was replicated in trans-ethnic analysis of 14,303 individuals (p = 2.2×10-16). The three SNPs were jointly associated with significantly reduced T2D risk, with an odds ratio (95% CI) of 0.93 (0.88, 0.98) per protective allele. Our findings implicate additive effects across pathophysiological pathways involved in type 2 diabetes, including glycolysis, gluconeogenesis, and insulin secretion. Since none of the individuals homozygous for the alternate alleles at all three loci has T2D, it might be possible to use a genetic predictor of fasting glucose levels to identify individuals at low vs. high risk of developing type 2 diabetes.

Evolutionary genetics and acclimatization in nephrology.

Evolutionary processes, including mutation, migration and natural selection, have influenced the prevalence and distribution of various disorders in humans. However, despite a few well-known examples, such as the APOL1 variants - which have undergone positive genetic selection for their ability to confer resistance to Trypanosoma brucei infection but confer a higher risk of chronic kidney disease - little is known about the effects of evolutionary processes that have shaped genetic variation on kidney disease. An understanding of basic concepts in evolutionary genetics provides an opportunity to consider how findings from ancient and archaic genomes could inform our knowledge of evolution and provide insights into how population migration and genetic admixture have shaped the current distribution and landscape of human kidney-associated diseases. Differences in exposures to infectious agents, environmental toxins, dietary components and climate also have the potential to influence the evolutionary genetics of kidneys. Of note, selective pressure on loci associated with kidney disease is often from non-kidney diseases, and thus it is important to understand how the link between genome-wide selected loci and kidney disease occurs in relation to secondary nephropathies.

GWAS in Africans identifies novel lipids loci and demonstrates heterogenous association within Africa.

Serum lipids are biomarkers of cardiometabolic disease risk, and understanding genomic factors contributing to their distribution is of interest. Studies of lipids in Africans are rare, though it is expected that such studies could identify novel loci. We conducted a GWAS of 4317 Africans enrolled from Nigeria, Ghana and Kenya. We evaluated linear mixed models of high-density lipoprotein cholesterol (HDLC), low-density lipoprotein cholesterol (LDLC), total cholesterol (CHOL), triglycerides (TG) and TG/HDLC. Replication was attempted in 9542 African Americans (AA). In our main analysis, we identified 28 novel associations in Africans. Of the 18 of these that could be tested in AA, three associations replicated (GPNMB-TG, ENPP1-TG and SMARCA4-LDLC). Five additional novel loci were discovered upon meta-analysis with AA (rs138282551-TG, PGBD5-HDLC, CD80-TG/HDLC, SLC44A1-CHOL and TLL2-CHOL). Analyses considering only those with predominantly West African ancestry (Nigeria, Ghana and AA) yielded new insights: ORC5-LDLC and chr20:60973327-CHOL. Among our novel findings are some loci with known connections to lipids pathways. For instance, rs147706369 (TLL2) alters a regulatory motif for sterol regulatory element-binding proteins, a family of transcription factors that control the expression of a range of enzymes involved in cholesterol, fatty acid and TG synthesis, and rs115749422 (SMARCA4), an independent association near the known LDLR locus that is rare or absent in populations without African ancestry. These findings demonstrate the utility of conducting genomic analyses in Africans for discovering novel loci and provide some preliminary evidence for caution against treating 'African ancestry' as a monolithic category.

A UGT1A1 variant is associated with serum total bilirubin levels, which are causal for hypertension in African-ancestry individuals.

Serum bilirubin is associated with several clinical outcomes, including hypertension, type 2 diabetes (T2D), and drug metabolism. Here, we describe findings from our genome-wide association studies (GWAS) of serum (TBIL) using a generalized linear mixed model in West Africans (n = 1127), with adjustment for age, sex, body mass index, T2D, significant principal components of population structure, and cryptic relatedness. Genome-wide conditional analysis and CAVIARBF were used to fine map significant loci. The causal effect of TBIL on hypertension was assessed by Mendelian randomization (MR) using the GWAS findings as instrumental variables (IVs) in African Americans (n = 3,067). The SNP rs887829 (UGT1A1) was significantly associated with TBIL levels (effect allele (T) frequency = 0.49, ß (SE) = 0.59 (0.04), p = 9.13 × 10-54). Genome-wide conditional analysis and regional fine mapping pointed to rs887829 as a possible causal variant with a posterior inclusion probability of 0.99. The T allele of rs887829 is associated with lower hepatic expression of UGT1A1. Using rs887829 as an IV, two-stage least-squares MR showed a causal effect of bilirubin on hypertension (ß = -0.76, 95% CI [-1.52, -0.01], p = 0.0459). Our finding confirms that UGT1A1 influences bilirubin levels. Notably, lower TBIL is causally associated with the increased risk of hypertension.

Genetic risk scores for cardiometabolic traits in sub-Saharan African populations.

BACKGROUND: There is growing support for the use of genetic risk scores (GRS) in routine clinical settings. Due to the limited diversity of current genomic discovery samples, there are concerns that the predictive power of GRS will be limited in non-European ancestry populations. GRS for cardiometabolic traits were evaluated in sub-Saharan Africans in comparison with African Americans and European Americans. METHODS: We evaluated the predictive utility of GRS for 12 cardiometabolic traits in sub-Saharan Africans (AF; n = 5200), African Americans (AA; n = 9139) and European Americans (EUR; n = 9594). GRS were constructed as weighted sums of the number of risk alleles. Predictive utility was assessed using the additional phenotypic variance explained and the increase in discriminatory ability over traditional risk factors [age, sex and body mass index (BMI)], with adjustment for ancestry-derived principal components. RESULTS: Across all traits, GRS showed up to a 5-fold and 20-fold greater predictive utility in EUR relative to AA and AF, respectively. Predictive utility was most consistent for lipid traits, with percentage increase in explained variation attributable to GRS ranging from 10.6% to 127.1% among EUR, 26.6% to 65.8% among AA and 2.4% to 37.5% among AF. These differences were recapitulated in the discriminatory power, whereby the predictive utility of GRS was 4-fold greater in EUR relative to AA and up to 44-fold greater in EUR relative to AF. Obesity and blood pressure traits showed a similar pattern of greater predictive utility among EUR. CONCLUSIONS: This work demonstrates the poorer performance of GRS in AF and highlights the need to improve representation of multiple ethnic populations in genomic studies to ensure equitable clinical translation of GRS.

Trans-ethnic meta-analysis identifies new loci associated with longitudinal blood pressure traits.

Genome-wide association studies (GWAS) have identified thousands of genetic loci associated with cross-sectional blood pressure (BP) traits; however, GWAS based on longitudinal BP have been underexplored. We performed ethnic-specific and trans-ethnic GWAS meta-analysis using longitudinal and cross-sectional BP data of 33,720 individuals from five cohorts in the US and one in Brazil. In addition to identifying several known loci, we identified thirteen novel loci with nine based on longitudinal and four on cross-sectional BP traits. Most of the novel loci were ethnic- or study-specific, with the majority identified in African Americans (AA). Four of these discoveries showed additional evidence of association in independent datasets, including an intergenic variant (rs4060030, p = 7.3 × 10-9) with reported regulatory function. We observed a high correlation between the meta-analysis results for baseline and longitudinal average BP (rho = 0.48). BP trajectory results were more correlated with those of average BP (rho = 0.35) than baseline BP(rho = 0.18). Heritability estimates trended higher for longitudinal traits than for cross-sectional traits, providing evidence for different genetic architectures. Furthermore, the longitudinal data identified up to 20% more BP known associations than did cross-sectional data. Our analyses of longitudinal BP data in diverse ethnic groups identified novel BP loci associated with BP trajectory, indicating a need for further longitudinal GWAS on BP and other age-related traits.

Time-to-event modeling of hypertension reveals the nonexistence of true controls.

Given a lifetime risk of ~90% by the ninth decade of life, it is unknown if there are true controls for hypertension in epidemiological and genetic studies. Here, we compared Bayesian logistic and time-to-event approaches to modeling hypertension. The median age at hypertension was approximately a decade earlier in African Americans than in European Americans or Mexican Americans. The probability of being free of hypertension at 85 years of age in African Americans was less than half that in European Americans or Mexican Americans. In all groups, baseline hazard rates increased until nearly 60 years of age and then decreased but did not reach zero. Taken together, modeling of the baseline hazard function of hypertension suggests that there are no true controls and that controls in logistic regression are cases with a late age of onset.

High-depth African genomes inform human migration and health.

The African continent is regarded as the cradle of modern humans and African genomes contain more genetic variation than those from any other continent, yet only a fraction of the genetic diversity among African individuals has been surveyed1. Here we performed whole-genome sequencing analyses of 426 individuals-comprising 50 ethnolinguistic groups, including previously unsampled populations-to explore the breadth of genomic diversity across Africa. We uncovered more than 3 million previously undescribed variants, most of which were found among individuals from newly sampled ethnolinguistic groups, as well as 62 previously unreported loci that are under strong selection, which were predominantly found in genes that are involved in viral immunity, DNA repair and metabolism. We observed complex patterns of ancestral admixture and putative-damaging and novel variation, both within and between populations, alongside evidence that Zambia was a likely intermediate site along the routes of expansion of Bantu-speaking populations. Pathogenic variants in genes that are currently characterized as medically relevant were uncommon-but in other genes, variants denoted as 'likely pathogenic' in the ClinVar database were commonly observed. Collectively, these findings refine our current understanding of continental migration, identify gene flow and the response to human disease as strong drivers of genome-level population variation, and underscore the scientific imperative for a broader characterization of the genomic diversity of African individuals to understand human ancestry and improve health.

Genetic modifiers of long-term survival in sickle cell anemia.

BACKGROUND: Sickle cell anemia (SCA) is a clinically heterogeneous, monogenic disorder. Medical care has less-than-optimal impact on clinical outcomes in SCA in Africa due to several factors, including patient accessibility, poor access to resources, and non-availability of specific effective interventions for SCA. METHODS: Against this background, we investigated 192 African participants who underwent whole exome sequencing. Participants included 105 SCA patients spanning variable clinical expression: a "long survivor" group (age over 40 years), a "stroke" group (at least one episode of overt stroke), and a "random" group (patients younger than 40 years without overt cerebrovascular disease). Fifty-eight ethnically matched homozygous hemoglobin A controls were also studied. Findings were validated in an independently recruited sample of 29 SCA patients. Statistical significance of the mutational burden of deleterious and loss-of-function variants per gene against a null model was estimated for each group, and gene-set association tests were conducted to test differences between groups. RESULTS: In the "long survivor" group, deleterious/loss-of-function variants were enriched in genes including CLCN6 (a voltage-dependent chloride channel for which rare deleterious variants have been associated with lower blood pressure) and OGHDL (important in arginine metabolism, which is a therapeutic target in SCA). In the "stroke" group, significant genes implicated were associated with increased activity of the blood coagulation cascade and increased complement activation, for example, SERPINC1, which encodes antithrombin. Oxidative stress and glutamate biosynthesis pathways were enriched in "long survivors" group. Published transcriptomic evidence provides functional support for the role of the identified pathways. CONCLUSIONS: This study provides new gene sets that contribute to variability in clinical expression of SCA. Identified genes and pathways suggest new avenues for other interventions.

Ancient familial Mediterranean fever mutations in human pyrin and resistance to Yersinia pestis.

Familial Mediterranean fever (FMF) is an autoinflammatory disease caused by homozygous or compound heterozygous gain-of-function mutations in MEFV, which encodes pyrin, an inflammasome protein. Heterozygous carrier frequencies for multiple MEFV mutations are high in several Mediterranean populations, suggesting that they confer selective advantage. Among 2,313 Turkish people, we found extended haplotype homozygosity flanking FMF-associated mutations, indicating evolutionarily recent positive selection of FMF-associated mutations. Two pathogenic pyrin variants independently arose >1,800 years ago. Mutant pyrin interacts less avidly with Yersinia pestis virulence factor YopM than with wild-type human pyrin, thereby attenuating YopM-induced interleukin (IL)-1ß suppression. Relative to healthy controls, leukocytes from patients with FMF harboring homozygous or compound heterozygous mutations and from asymptomatic heterozygous carriers released heightened IL-1ß specifically in response to Y. pestis. Y. pestis-infected MefvM680I/M680I FMF knock-in mice exhibited IL-1-dependent increased survival relative to wild-type knock-in mice. Thus, FMF mutations that were positively selected in Mediterranean populations confer heightened resistance to Y. pestis.

Admixture mapping identifies genetic regions associated with blood pressure phenotypes in African Americans.

Hypertension occurs at a higher rate in African Americans than in European Americans. Based on the assumption that causal variants are more frequently found on DNA segments inherited from the ancestral population with higher disease risk, we employed admixture mapping to identify genetic loci with excess local African ancestry associated with blood pressure. Chromosomal regions 1q21.2-21.3, 4p15.1, 19q12 and 20p13 were significantly associated with diastolic blood pressure (ß = 5.28, -7.94, -6.82 and 5.89, P-value = 6.39E-04, 2.07E-04, 6.56E-05 and 5.04E-04, respectively); 1q21.2-21.3 and 19q12 were also significantly associated with mean arterial pressure (ß = 5.86 and -6.40, P-value = 5.32E-04 and 6.37E-04, respectively). We further selected SNPs that had large allele frequency differences within these regions and tested their association with blood pressure. SNP rs4815428 was significantly associated with diastolic blood pressure after Bonferroni correction (ß = -2.42, P-value = 9.57E-04), and it partially explained the admixture mapping signal at 20p13. SNPs rs771205 (ß = -1.99, P-value = 3.37E-03), rs3126067, rs2184953 and rs58001094 (the latter three exhibit strong linkage disequilibrium, ß = -2.3, P-value = 1.4E-03) were identified to be significantly associated with mean arterial pressure, and together they fully explained the admixture signal at 1q21.2-21.3. Although no SNP at 4p15.1 showed large ancestral allele frequency differences in our dataset, we detected association at low-frequency African-specific variants that mapped predominantly to the gene PCDH7, which is most highly expressed in aorta. Our results suggest that these regions may harbor genetic variants that contribute to the different prevalence of hypertension.