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.

Perceptions and beliefs of community gatekeepers about genomic risk information in African cleft research.

BACKGROUND: A fundamental ethical issue in African genomics research is how socio-cultural factors impact perspectives, acceptance, and utility of genomic information, especially in stigmatizing conditions like orofacial clefts (OFCs). Previous research has shown that gatekeepers (e.g., religious, political, family or community leaders) wield considerable influence on the decision-making capabilities of their members, including health issues. Thus, their perspectives can inform the design of engagement strategies and increase exposure to the benefits of genomics testing/research. This is especially important for Africans underrepresented in genomic research. Our study aims to investigate the perspectives of gatekeepers concerning genomic risk information (GRI) in the presence of OFCs in a sub-Saharan African cohort. METHODS: Twenty-five focus group discussions (FGDs) consisting of 214 gatekeepers (religious, community, ethnic leaders, and traditional birth attendants) in Lagos, Nigeria, explored the opinions of participants on genomic risk information (GRI), OFC experience, and the possibility of involvement in collaborative decision-making in Lagos, Nigeria. Transcripts generated from audio recordings were coded and analyzed in NVivo using thematic analysis. RESULTS: Three main themes-knowledge, beliefs, and willingness to act-emerged from exploring the perspective of gatekeepers about GRI in this group. We observed mixed opinions regarding the acceptance of GRI. Many participants believed their role is to guide and support members when they receive results; this is based on the level of trust their members have in them. However, participants felt they would need to be trained by medical experts to do this. Also, religious and cultural beliefs were crucial to determining participants' understanding of OFCs and the acceptance and utilization of GRI. CONCLUSIONS: Incorporating cultural sensitivity into public engagement could help develop appropriate strategies to manage conflicting ideologies surrounding genomic information in African communities. This will allow for more widespread access to the advances in genomics research in underrepresented populations. We also recommend a synergistic relationship between community health specialists/scientists, and community leaders, including spiritual providers to better understand and utilize GRI.

Parents and Provider Perspectives on the Return of Genomic Findings for Cleft Families in Africa.

BACKGROUND: Inadequate knowledge among health care providers (HCPs) and parents of affected children limits the understanding and utility of secondary genetic findings (SFs) in under-represented populations in genomics research. SFs arise from deep DNA sequencing done for research or diagnostic purposes and may burden patients and their families despite their potential health importance. This study aims to evaluate the perspective of both groups regarding SFs and their choices in the return of results from genetic testing in the context of orofacial clefts. METHODS: Using an online survey, we evaluated the experiences of 252 HCPs and 197 parents across participating cleft clinics in Ghana and Nigeria toward the return of SFs across several domains. RESULTS: Only 1.6% of the HCPs felt they had an expert understanding of when and how to incorporate genomic medicine into practice, while 50.0% agreed that all SFs should be returned to patients. About 95.4% of parents were willing to receive all the information from genetic testing (including SFs), while the majority cited physicians as their primary information source (64%). CONCLUSIONS: Overall, parents and providers were aware that genetic testing could help in the clinical management of diseases. However, they cited a lack of knowledge about genomic medicine, uncertain clinical utility, and lack of available learning resources as barriers. The knowledge gained from this study will assist with developing guidelines and policies to guide providers on the return of SFs in sub-Saharan Africa and across the continent.

Sickle cell allele HBB-rs334(T) is associated with decreased risk of childhood Burkitt lymphoma in East Africa.

Burkitt lymphoma (BL) is an aggressive B-cell lymphoma that significantly contributes to childhood cancer burden in sub-Saharan Africa. Plasmodium falciparum, which causes malaria, is geographically associated with BL, but the evidence remains insufficient for causal inference. Inference could be strengthened by demonstrating that mendelian genes known to protect against malaria-such as the sickle cell trait variant, HBB-rs334(T)-also protect against BL. We investigated this hypothesis among 800 BL cases and 3845 controls in four East African countries using genome-scan data to detect polymorphisms in 22 genes known to affect malaria risk. We fit generalized linear mixed models to estimate odds ratios (OR) and 95% confidence intervals (95% CI), controlling for age, sex, country, and ancestry. The ORs of the loci with BL and P. falciparum infection among controls were correlated (Spearman's ρ = 0.37, p = .039). HBB-rs334(T) was associated with lower P. falciparum infection risk among controls (OR = 0.752, 95% CI 0.628-0.9; p = .00189) and BL risk (OR = 0.687, 95% CI 0.533-0.885; p = .0037). ABO-rs8176703(T) was associated with decreased risk of BL (OR = 0.591, 95% CI 0.379-0.992; p = .00271), but not of P. falciparum infection. Our results increase support for the etiological correlation between P. falciparum and BL risk.

Mendelian randomization analyses suggest a causal role for circulating GIP and IL-1RA levels in homeostatic model assessment-derived measures of ß-cell function and insulin sensitivity in Africans without type 2 diabetes.

BACKGROUND: In vitro and in vivo studies have shown that certain cytokines and hormones may play a role in the development and progression of type 2 diabetes (T2D). However, studies on their role in T2D in humans are scarce. We evaluated associations between 11 circulating cytokines and hormones with T2D among a population of sub-Saharan Africans and tested for causal relationships using Mendelian randomization (MR) analyses. METHODS: We used logistic regression analysis adjusted for age, sex, body mass index, and recruitment country to regress levels of 11 cytokines and hormones (adipsin, leptin, visfatin, PAI-1, GIP, GLP-1, ghrelin, resistin, IL-6, IL-10, IL-1RA) on T2D among Ghanaians, Nigerians, and Kenyans from the Africa America Diabetes Mellitus study including 2276 individuals with T2D and 2790 non-T2D individuals. Similar linear regression models were fitted with homeostatic modelling assessments of insulin sensitivity (HOMA-S) and ß-cell function (HOMA-B) as dependent variables among non-T2D individuals (n = 2790). We used 35 genetic variants previously associated with at least one of these 11 cytokines and hormones among non-T2D individuals as instrumental variables in univariable and multivariable MR analyses. Statistical significance was set at 0.0045 (0.05/11 cytokines and hormones). RESULTS: Circulating GIP and IL-1RA levels were associated with T2D. Nine of the 11 cytokines and hormones (exceptions GLP-1 and IL-6) were associated with HOMA-S, HOMA-B, or both among non-T2D individuals. Two-stage least squares MR analysis provided evidence for a causal effect of GIP and IL-RA on HOMA-S and HOMA-B in multivariable analyses (GIP ~ HOMA-S ß = - 0.67, P-value = 1.88 × 10-6 and HOMA-B ß = 0.59, P-value = 1.88 × 10-5; IL-1RA ~ HOMA-S ß = - 0.51, P-value = 8.49 × 10-5 and HOMA-B ß = 0.48, P-value = 5.71 × 10-4). IL-RA was partly mediated via BMI (30-34%), but GIP was not. Inverse variance weighted MR analysis provided evidence for a causal effect of adipsin on T2D (multivariable OR = 1.83, P-value = 9.79 × 10-6), though these associations were not consistent in all sensitivity analyses. CONCLUSIONS: The findings of this comprehensive MR analysis indicate that circulating GIP and IL-1RA levels are causal for reduced insulin sensitivity and increased ß-cell function. GIP's effect being independent of BMI suggests that circulating levels of GIP could be a promising early biomarker for T2D risk. Our MR analyses do not provide conclusive evidence for a causal role of other circulating cytokines in T2D among sub-Saharan Africans.

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.

Clinically actionable secondary findings in 130 triads from sub-Saharan African families with non-syndromic orofacial clefts.

INTRODUCTION: The frequency and implications of secondary findings (SFs) from genomic testing data have been extensively researched. However, little is known about the frequency or reporting of SFs in Africans, who are underrepresented in large-scale population genomic studies. The availability of data from the first whole-genome sequencing for orofacial clefts in an African population motivated this investigation. METHODS: In total, 130 case-parent trios were analyzed for SFs within the ACMG SFv.3.0 list genes. Additionally, we filtered for four more genes (HBB, HSD32B, G6PD and ACADM). RESULTS: We identified 246 unique variants in 55 genes; five variants in four genes were classified as pathogenic or likely pathogenic (P/LP). The P/LP variants were seen in 2.3% (9/390) of the subjects, a frequency higher than ~1% reported for diverse ethnicities. On the ACMG list, pathogenic variants were observed in PRKAG (p. Glu183Lys). Variants in the PALB2 (p. Glu159Ter), RYR1 (p. Arg2163Leu) and LDLR (p. Asn564Ser) genes were predicted to be LP. CONCLUSION: This study provides information on the frequency and pathogenicity of SFs in an African cohort. Early risk detection will help reduce disease burden and contribute to efforts to increase knowledge of the distribution and impact of actionable genomic variants in diverse populations.

Ambiguous genitalia, giant congenital melanocytic nevus and subpulmonic outlet ventricular septal defect in an African child with Neurofibromatosis 1.

Neurofibromatosis type 1 is an autosomal dominant multisystemic disease caused by mutation of the neurofibromin (NF1) gene located on chromosome 17q11. We report a case of Neurofibromatosis 1 with ambiguous genitalia, giant congenital melanocytic nevus, and associated subpulmonic outlet ventricular septal defect, hitherto unreported in sub-Saharan Africa. In addition, a literature review of congenital heart diseases associated with Neurofibromatosis 1 is presented.

Ancestral and environmental patterns in the association between triglycerides and other cardiometabolic risk factors.

BACKGROUND: West Africans and African Americans with substantial (∼80%) West African ancestry are characterized by low levels of triglycerides (TG) compared to East Africans and Europeans. The impact of these varying TG levels on other cardiometabolic risk factors is unclear. We compared the strength of association between TG with hypertension, blood pressure, BMI, waist circumference, type 2 diabetes (T2D), and fasting glucose across West African (WA), East African (EA), and European (EU) ancestry populations residing in three vastly different environmental settings: sub-Saharan Africa, United States, and Europe. METHODS: We analysed data from four cross-sectional studies that included WA in sub-Saharan Africa (n = 7201), the U.S. (n = 4390), and Europe (n = 6436), EA in sub-Saharan Africa (n = 781), and EU in the U.S. (n = 8670) and Europe (n = 4541). Linear regression analyses were used to test the association between TG and cardiometabolic risk factors. FINDINGS: Higher adjusted regression coefficients were observed in EU compared with WA ancestry for TG on hypertension (EU ß [95% CI]: 0.179 [0.156, 0.203], WA ß [95% CI]: 0.102 [0.086, 0.118]), BMI (EU ß [95% CI]: 0.028 [0.027, 0.030], WA ß [95% CI]: 0.015 [0.014, 0.016]), and waist circumference (EU ß [95% CI]: 0.013 [0.013, 0.014], WA ß [95% CI]: 0.009 [0.008, 0.009) (all ancestry × trait interaction P-values <0.05), irrespective of environmental differences within ancestry groups. Less consistency was observed among EA. Associations of TG with T2D did not follow ancestry patterns, with substantial variation observed between environments. INTERPRETATION: TG may not be an equally strong associated with other established cardiometabolic risk factors in West and East Africans in contrast to European ancestry populations. The value of TG for identifying individuals at high risk for developing metabolic disorders needs to be re-evaluated for African ancestry populations. FUNDING: National Institutes of Health, European Commission, Dutch Heart Foundation, Netherlands Organization for Health Research and Development, Centers for Disease Control and Prevention.

Cardiovascular Risk Factor Burden and Association With CKD in Ghana and Nigeria.

Introduction: Cardiovascular disease is the leading cause of morbidity and mortality in patients with chronic kidney disease (CKD); however, the burden of cardiovascular risk factors in patients with CKD in Africa is not well characterized. We determined the prevalence of selected cardiovascular risk factors, and association with CKD in the Human Heredity for Health in Africa Kidney Disease Research Network study. Methods: We recruited patients with and without CKD in Ghana and Nigeria. CKD was defined as estimated glomerular filtration rate of <60 ml/min per 1.73 m2 and/or albuminuria as albumin-to-creatinine ratio <3.0 mg/mmol (<30 mg/g) for ≥3 months. We assessed self-reported (physician-diagnosis and/or use of medication) hypertension, diabetes, and elevated cholesterol; and self-reported smoking as cardiovascular risk factors. Association between the risk factors and CKD was determined by multivariate logistic regression. Results: We enrolled 8396 participants (cases with CKD, 3956), with 56% females. The mean age (45.5 ± 15.1 years) did not differ between patients and control group. The prevalence of hypertension (59%), diabetes (20%), and elevated cholesterol (9.9%), was higher in CKD patients than in the control participants (P < 0.001). Prevalence of risk factors was higher in Ghana than in Nigeria. Hypertension (adjusted odds ratio [aOR] = 1.69 [1.43-2.01, P < 0.001]), elevated cholesterol (aOR = 2.0 [1.39-2.86, P < 0.001]), age >50 years, and body mass index (BMI) <18.5 kg/m2 were independently associated with CKD. The association of diabetes and smoking with CKD was modified by other risk factors. Conclusion: Cardiovascular risk factors are prevalent in middle-aged adult patients with CKD in Ghana and Nigeria, with higher proportions in Ghana than in Nigeria. Hypertension, elevated cholesterol, and underweight were independently associated with CKD.

Hajdu-Cheney syndrome with atypical cardiovascular abnormalities.

Hajdu-Cheney syndrome is an ultra-rare autosomal dominant disorder caused by a heterozygous variant in NOTCH2 gene. Characteristic features include osteolysis, distinct facial appearance, skull deformity, joint laxity, osteoporosis, and short stature. Associated abnormalities are congenital heart disease, congenital defects of the kidney, and neurological problems. Here, we present the first reported case of an African child with a variant in NOTCH2 gene and features of Hajdu-Cheney syndrome in whom we detected a congenital heart defect that has not been previously reported in association with the syndrome. To appropriately characterize this disease and document correct proportion of cardiovascular malformation associations, echocardiography is recommended for all cases of Hajdu Cheney syndrome.

Damaging Mutations in AFDN Contribute to Risk of Nonsyndromic Cleft Lip With or Without Cleft Palate.

Novel or rare damaging mutations have been implicated in the developmental pathogenesis of nonsyndromic cleft lip with or without cleft palate (nsCL ± P). Thus, we investigated the human genome for high-impact mutations that could explain the risk of nsCL ± P in our cohorts.We conducted next-generation sequencing (NGS) analysis of 130 nsCL ± P case-parent African trios to identify pathogenic variants that contribute to the risk of clefting. We replicated this analysis using whole-exome sequence data from a Brazilian nsCL ± P cohort. Computational analyses were then used to predict the mechanism by which these variants could result in increased risks for nsCL ± P.We discovered damaging mutations within the AFDN gene, a cell adhesion molecule (CAMs) that was previously shown to contribute to cleft palate in mice. These mutations include p.Met1164Ile, p.Thr453Asn, p.Pro1638Ala, p.Arg669Gln, p.Ala1717Val, and p.Arg1596His. We also discovered a novel splicing p.Leu1588Leu mutation in this protein. Computational analysis suggests that these amino acid changes affect the interactions with other cleft-associated genes including nectins (PVRL1, PVRL2, PVRL3, and PVRL4) CDH1, CTNNA1, and CTNND1.This is the first report on the contribution of AFDN to the risk for nsCL ± P in humans. AFDN encodes AFADIN, an important CAM that forms calcium-independent complexes with nectins 1 and 4 (encoded by the genes PVRL1 and PVRL4). This discovery shows the power of NGS analysis of multiethnic cleft samples in combination with a computational approach in the understanding of the pathogenesis of nsCL ± P.

Mendelian randomization study reveals a causal relationship between adiponectin and LDL cholesterol in Africans.

Adiponectin has been associated with cardiometabolic traits in observational studies across populations, yet it is unclear if these associations are causal. We performed Mendelian randomization (MR) analysis to assess the relationship between adiponectin and cardiometabolic traits in sub-Saharan Africans. We constructed a polygenic risk score (PRS) for adiponectin levels across 3354 unrelated sub-Saharan Africans. The PRS was used as the instrumental variable in two-stage least-squares MR analysis to assess its association with insulin resistance, HDL, LDL, total cholesterol, triglycerides, blood pressure, Type 2 Diabetes (T2D), and hypertension. The adiponectin PRS was causally related with LDL (ß = 0.55, 95%CI 0.07-1.04, P-value = 0.024) but not the other traits. This association was observed in both overweight/obese and normal weight individuals, but only reached statistical significance among overweight/obese individuals (ß = 0.55, 95%CI 0.01-1.08, P-value = 0.045). In normal weight individuals, the adiponectin PRS was associated with T2D (OR = 0.13, 95%CI 0.02-0.73, P-value = 0.021), and in men with HDL (ß = 1.03, 95%CI 0.14-1.92, P-value = 0.023). The findings of this first MR study in sub-Saharan Africans support a causal relationship of adiponectin with LDL, with T2D in normal weight individuals only, and with HDL in men only. These observations add to the small but growing literature on adiponectin MR studies.

Whole-genome sequencing reveals de-novo mutations associated with nonsyndromic cleft lip/palate.

The majority (85%) of nonsyndromic cleft lip with or without cleft palate (nsCL/P) cases occur sporadically, suggesting a role for de novo mutations (DNMs) in the etiology of nsCL/P. To identify high impact protein-altering DNMs that contribute to the risk of nsCL/P, we conducted whole-genome sequencing (WGS) analyses in 130 African case-parent trios (affected probands and unaffected parents). We identified 162 high confidence protein-altering DNMs some of which are based on available evidence, contribute to the risk of nsCL/P. These include novel protein-truncating DNMs in the ACTL6A, ARHGAP10, MINK1, TMEM5 and TTN genes; as well as missense variants in ACAN, DHRS3, DLX6, EPHB2, FKBP10, KMT2D, RECQL4, SEMA3C, SEMA4D, SHH, TP63, and TULP4. Many of these protein-altering DNMs were predicted to be pathogenic. Analysis using mouse transcriptomics data showed that some of these genes are expressed during the development of primary and secondary palate. Gene-set enrichment analysis of the protein-altering DNMs identified palatal development and neural crest migration among the few processes that were significantly enriched. These processes are directly involved in the etiopathogenesis of clefting. The analysis of the coding sequence in the WGS data provides more evidence of the opportunity for novel findings in the African genome.

An epigenome-wide association study of insulin resistance in African Americans.

BACKGROUND: African Americans have a high risk for type 2 diabetes (T2D) and insulin resistance. Studies among other population groups have identified DNA methylation loci associated with insulin resistance, but data in African Americans are lacking. Using DNA methylation profiles of blood samples obtained from the Illumina Infinium® HumanMethylation450 BeadChip, we performed an epigenome-wide association study to identify DNA methylation loci associated with insulin resistance among 136 non-diabetic, unrelated African American men (mean age 41.6 years) from the Howard University Family Study. RESULTS: We identified three differentially methylated positions (DMPs) for homeostatic model assessment of insulin resistance (HOMA-IR) at 5% FDR. One DMP (cg14013695, HOXA5) is a known locus among Mexican Americans, while the other two DMPs are novel-cg00456326 (OSR1; beta = 0.027) and cg20259981 (ST18; beta = 0.010). Although the cg00456326 DMP is novel, the OSR1 gene has previously been found associated with both insulin resistance and T2D in Europeans. The genes HOXA5 and ST18 have been implicated in biological processes relevant to insulin resistance. Differential methylation at the significant HOXA5 and OSR1 DMPs is associated with differences in gene expression in the iMETHYL database. Analysis of differentially methylated regions (DMRs) did not identify any epigenome-wide DMRs for HOMA-IR. We tested transferability of HOMA-IR associated DMPs from five previous EWAS in Mexican Americans, Indian Asians, Europeans, and European ancestry Americans. Out of the 730 previously reported HOMA-IR DMPs, 47 (6.4%) were associated with HOMA-IR in this cohort of African Americans. CONCLUSIONS: The findings from our study suggest substantial differences in DNA methylation patterns associated with insulin resistance across populations. Two of the DMPs we identified in African Americans have not been reported in other populations, and we found low transferability of HOMA-IR DMPs reported in other populations in African Americans. More work in African-ancestry populations is needed to confirm our findings as well as functional analyses to understand how such DNA methylation alterations contribute to T2D pathology.

Genome-Wide Scan for Parent-of-Origin Effects in a sub-Saharan African Cohort With Nonsyndromic Cleft Lip and/or Cleft Palate (CL/P).

OBJECTIVE: Nonsyndromic cleft lip and/or cleft palate (NSCL/P) have multifactorial etiology where genetic factors, gene-environment interactions, stochastic factors, gene-gene interactions, and parent-of-origin effects (POEs) play cardinal roles. POEs arise when the parental origin of alleles differentially impacts the phenotype of the offspring. The aim of this study was to identify POEs that can increase risk for NSCL/P in humans using a genome-wide dataset. METHODS: The samples (174 case-parent trios from Ghana, Ethiopia, and Nigeria) included in this study were from the African only genome wide association studies (GWAS) that was published in 2019. Genotyping of individual DNA using over 2 million multiethnic and African ancestry-specific single-nucleotide polymorphisms from the Illumina Multi-Ethnic Genotyping Array v2 15070954 A2 (genome build GRCh37/hg19) was done at the Center for Inherited Diseases Research. After quality control checks, PLINK was employed to carry out POE analysis employing the pooled subphenotypes of NSCL/P. RESULTS: We observed possible hints of POEs at a cluster of genes at a 1 mega base pair window at the major histocompatibility complex class 1 locus on chromosome 6, as well as at other loci encompassing candidate genes such as ASB18, ANKEF1, AGAP1, GABRD, HHAT, CCT7, DNMT3A, EPHA7, FOXO3, lncRNAs, microRNA, antisense RNAs, ZNRD1, ZFAT, and ZBTB16. CONCLUSION: Findings from our study suggest that some loci may increase the risk for NSCL/P through POEs. Additional studies are required to confirm these suggestive loci in NSCL/P etiology.

Variant analyses of candidate genes in orofacial clefts in multi-ethnic populations.

OBJECTIVES: Cleft lip with/without cleft palate and cleft palate only is congenital birth defects where the upper lip and/or palate fail to fuse properly during embryonic facial development. Affecting ~1.2/1000 live births worldwide, these orofacial clefts impose significant social and financial burdens on affected individuals and their families. Orofacial clefts have a complex etiology resulting from genetic variants combined with environmental covariates. Recent genome-wide association studies and whole-exome sequencing for orofacial clefts identified significant genetic associations and variants in several genes. Of these, we investigated the role of common/rare variants in SHH, RORA, MRPL53, ACVR1, and GDF11. MATERIALS AND METHODS: We sequenced these five genes in 1255 multi-ethnic cleft lip with/without palate and cleft palate only samples in order to find variants that may provide potential explanations for the missing heritability of orofacial clefts. Rare and novel variants were further analyzed using in silico predictive tools. RESULTS: Ninteen total variants of interest were found, with variant types including stop-gain, missense, synonymous, intronic, and splice-site variants. Of these, 3 novel missense variants were found, one in SHH, one in RORA, and one in GDF11. CONCLUSION: This study provides evidence that variants in SHH, RORA, MRPL53, ACVR1, and GDF11 may contribute to risk of orofacial clefts in various populations.

Epigenome-wide association study of serum urate reveals insights into urate co-regulation and the SLC2A9 locus.

Elevated serum urate levels, a complex trait and major risk factor for incident gout, are correlated with cardiometabolic traits via incompletely understood mechanisms. DNA methylation in whole blood captures genetic and environmental influences and is assessed in transethnic meta-analysis of epigenome-wide association studies (EWAS) of serum urate (discovery, n = 12,474, replication, n = 5522). The 100 replicated, epigenome-wide significant (p < 1.1E-7) CpGs explain 11.6% of the serum urate variance. At SLC2A9, the serum urate locus with the largest effect in genome-wide association studies (GWAS), five CpGs are associated with SLC2A9 gene expression. Four CpGs at SLC2A9 have significant causal effects on serum urate levels and/or gout, and two of these partly mediate the effects of urate-associated GWAS variants. In other genes, including SLC7A11 and PHGDH, 17 urate-associated CpGs are associated with conditions defining metabolic syndrome, suggesting that these CpGs may represent a blood DNA methylation signature of cardiometabolic risk factors. This study demonstrates that EWAS can provide new insights into GWAS loci and the correlation of serum urate with other complex traits.

Meta-analyses identify DNA methylation associated with kidney function and damage.

Chronic kidney disease is a major public health burden. Elevated urinary albumin-to-creatinine ratio is a measure of kidney damage, and used to diagnose and stage chronic kidney disease. To extend the knowledge on regulatory mechanisms related to kidney function and disease, we conducted a blood-based epigenome-wide association study for estimated glomerular filtration rate (n = 33,605) and urinary albumin-to-creatinine ratio (n = 15,068) and detected 69 and seven CpG sites where DNA methylation was associated with the respective trait. The majority of these findings showed directionally consistent associations with the respective clinical outcomes chronic kidney disease and moderately increased albuminuria. Associations of DNA methylation with kidney function, such as CpGs at JAZF1, PELI1 and CHD2 were validated in kidney tissue. Methylation at PHRF1, LDB2, CSRNP1 and IRF5 indicated causal effects on kidney function. Enrichment analyses revealed pathways related to hemostasis and blood cell migration for estimated glomerular filtration rate, and immune cell activation and response for urinary albumin-to-creatinineratio-associated CpGs.

Genome-wide analyses of multiple obesity-related cytokines and hormones informs biology of cardiometabolic traits.

BACKGROUND: A complex set of perturbations occur in cytokines and hormones in the etiopathogenesis of obesity and related cardiometabolic conditions such as type 2 diabetes (T2D). Evidence for the genetic regulation of these cytokines and hormones is limited, particularly in African-ancestry populations. In order to improve our understanding of the biology of cardiometabolic traits, we investigated the genetic architecture of a large panel of obesity- related cytokines and hormones among Africans with replication analyses in African Americans. METHODS: We performed genome-wide association studies (GWAS) in 4432 continental Africans, enrolled from Ghana, Kenya, and Nigeria as part of the Africa America Diabetes Mellitus (AADM) study, for 13 obesity-related cytokines and hormones, including adipsin, glucose-dependent insulinotropic peptide (GIP), glucagon-like peptide-1 (GLP-1), interleukin-1 receptor antagonist (IL1-RA), interleukin-6 (IL-6), interleukin-10 (IL-10), leptin, plasminogen activator inhibitor-1 (PAI-1), resistin, visfatin, insulin, glucagon, and ghrelin. Exact and local replication analyses were conducted in African Americans (n = 7990). The effects of sex, body mass index (BMI), and T2D on results were investigated through stratified analyses. RESULTS: GWAS identified 39 significant (P value < 5 × 10-8) loci across all 13 traits. Notably, 14 loci were African-ancestry specific. In this first GWAS for adipsin and ghrelin, we detected 13 and 4 genome-wide significant loci respectively. Stratified analyses by sex, BMI, and T2D showed a strong effect of these variables on detected loci. Eight novel loci were successfully replicated: adipsin (3), GIP (1), GLP-1 (1), and insulin (3). Annotation of these loci revealed promising links between these adipocytokines and cardiometabolic outcomes as illustrated by rs201751833 for adipsin and blood pressure and locus rs759790 for insulin level and T2D in lean individuals. CONCLUSIONS: Our study identified genetic variants underlying variation in multiple adipocytokines, including the first loci for adipsin and ghrelin. We identified population differences in variants associated with adipocytokines and highlight the importance of stratification for discovery of loci. The high number of African-specific loci detected emphasizes the need for GWAS in African-ancestry populations, as these loci could not have been detected in other populations. Overall, our work contributes to the understanding of the biology linking adipocytokines to cardiometabolic traits.