Search

Abstract Background Epigenetic modifications, particularly DNA methylation (DNAm) in cord blood, are an important biological marker of how external exposures during gestation can influence the in-utero environment and subsequent offspring development. Despite the recognized importance of DNAm during gestation, comparative studies to determine the consistency of these epigenetic signals across different ethnic groups are largely absent. To address this gap, we first performed epigenome-wide association studies (EWAS) of gestational age (GA) using newborn cord blood DNAm comparatively in a white European (n = 342) and a South Asian (n = 490) birth cohort living in Canada. Then, we capitalized on established cord blood epigenetic GA clocks to examine the associations between maternal exposures, offspring characteristics and epigenetic GA, as well as GA acceleration, defined as the residual difference between epigenetic and chronological GA at birth. Results Individual EWASs confirmed 1,211 and 1,543 differentially methylated CpGs previously reported to be associated with GA, in white European and South Asian cohorts, respectively, with a similar distribution of effects. We confirmed that Bohlin’s cord blood GA clock was robustly correlated with GA in white Europeans (r = 0.71; p = 6.0 × 10–54) and South Asians (r = 0.66; p = 6.9 × 10–64). In both cohorts, Bohlin’s clock was positively associated with newborn weight and length and negatively associated with parity, newborn female sex, and gestational diabetes. Exclusive to South Asians, the GA clock was positively associated with the newborn ponderal index, while pre-pregnancy weight and gestational weight gain were strongly predictive of increased epigenetic GA in white Europeans. Important predictors of GA acceleration included gestational diabetes mellitus, newborn sex, and parity in both cohorts. Conclusions These results demonstrate the consistent DNAm signatures of GA and the utility of Bohlin’s GA clock across the two populations. Although the overall pattern of DNAm is similar, its connections with the mother's environment and the baby's anthropometrics can differ between the two groups. Further research is needed to understand these unique relationships.

Abstract Background Polygenic prediction studies in continental Africans are scarce. Africa’s genetic and environmental diversity pose a challenge that limits the generalizability of polygenic risk scores (PRS) for body mass index (BMI) within the continent. Studies to understand the factors that affect PRS variability within Africa are required. Methods Using the first multi-ancestry genome-wide association study (GWAS) meta-analysis for BMI involving continental Africans, we derived a multi-ancestry PRS and compared its performance to a European ancestry-specific PRS in continental Africans (AWI-Gen study) and a European cohort (Estonian Biobank). We then evaluated the factors affecting the performance of the PRS in Africans which included fine-mapping resolution, allele frequencies, linkage disequilibrium patterns, and PRS-environment interactions. Results Polygenic prediction of BMI in continental Africans is poor compared to that in European ancestry individuals. However, we show that the multi-ancestry PRS is more predictive than the European ancestry-specific PRS due to its improved fine-mapping resolution. We noted regional variation in polygenic prediction across Africa’s East, South, and West regions, which was driven by a complex interplay of the PRS with environmental factors, such as physical activity, smoking, alcohol intake, and socioeconomic status. Conclusions Our findings highlight the role of gene-environment interactions in PRS prediction variability in Africa. PRS methods that correct for these interactions, coupled with the increased representation of Africans in GWAS, may improve PRS prediction in Africa.

Uromodulin (UMOD) is a major risk gene for monogenic and complex forms of kidney disease. The encoded kidney-specific protein uromodulin is highly abundant in urine and related to chronic kidney disease, hypertension, and pathogen defense. To gain insights into potential systemic roles, we performed genome-wide screens of circulating uromodulin using complementary antibody-based and aptamer-based assays. We detected 3 and 10 distinct significant loci, respectively. Integration of antibody-based results at the UMOD locus with functional genomics data (RNA-Seq, ATAC-Seq, Hi-C) of primary human kidney tissue highlighted an upstream variant with differential accessibility and transcription in uromodulin-synthesizing kidney cells as underlying the observed cis effect. Shared association patterns with complex traits, including chronic kidney disease and blood pressure, placed the PRKAG2 locus in the same pathway as UMOD. Experimental validation of the third antibody-based locus, B4GALNT2, showed that the p.Cys466Arg variant of the encoded N-acetylgalactosaminyltransferase had a loss-of-function effect leading to higher serum uromodulin levels. Aptamer-based results pointed to enzymes writing glycan marks present on uromodulin and to their receptors in the circulation, suggesting that this assay permits investigating uromodulin’s complex glycosylation rather than its quantitative levels. Overall, our study provides insights into circulating uromodulin and its emerging functions.

Abstract Background With the legalization of cannabis in Canada, there is an increase trend in use. Cannabis has been known to have several health implications, one of which is the development of cannabis use disorder (CUD). CUD is more common in males than females, as well as in certain ethnic groups such as Native Americans. Additionally, both environmental and genetic risk factors have been found for cannabis use. The objective of this systematic review will be to summarize the genetic variants associated with cannabis use which have reached borderline genome-wide significance. Methods This systematic review will incorporate articles that have performed a genome-wide association study (GWAS) investigating cannabis use. MEDLINE, Web of Science, EMBASE, GWAS Catalog, GWAS Central, and NIH Database of Genotype and Phenotype will be searched using a comprehensive search strategy. The quality of genetic association studies (Q-Genie) tool will be utilized to assess the quality of the included studies. All screening and data extraction will occur independently by two authors. If feasible, a random-effects meta-analysis will be conducted on pooled odds ratios of single nucleotide polymorphisms reaching borderline genome-wide significance. Discussion This systematic review will synthesize available GWAS on cannabis use. Results from this review will inform and direct further investigation of genetic variants associated with cannabis use. Systematic review registration PROSPERO CRD42020176016

Abstract Background Hearing impairment (HI) genes are poorly studied in African populations. Methods We used whole exome sequencing (WES) to investigate pathogenic and likely pathogenic (PLP) variants in 10 individuals with HI, from four multiplex families from Cameroon, two of which were previously unresolved with a targeted gene enrichment (TGE) panel of 116 genes. In silico protein modelling, western blotting and live imaging of transfected HEK293 cells were performed to study protein structure and functions. Results All PLP variants previously identified with TGE were replicated. In one previously unresolved family, we found a homozygous frameshift PLP variant in GRXCR2 (OMIM: 615762), NM_001080516.1(GRXCR2):c.251delC p.(Ile85SerfsTer33), in two affected siblings; and additionally, in 1/80 unrelated individuals affected with non‐syndromic hearing impairment (NSHI). The GRXCR2‐c.251delC variant introduced a premature stop codon, leading to truncation and loss of a zinc‐finger domain. Fluorescence confocal microscopy tracked the wild‐type GRXCR2 protein to the cellular membrane, unlike the mutated GRXCR2 protein. Conclusion This study confirms GRXCR2 as a HI‐associated gene. GRXCR2 should be included to the currently available TGE panels for HI diagnosis.

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10(-8)-1.2×10(-43)). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p

Summary: Chronic kidney disease is a leading cause of death and disability globally and impacts individuals of African ancestry (AFR) or with ancestry in the Americas (AMS) who are under-represented in genome-wide association studies (GWASs) of kidney function. To address this bias, we conducted a large meta-analysis of GWASs of estimated glomerular filtration rate (eGFR) in 145,732 AFR and AMS individuals. We identified 41 loci at genome-wide significance (p

Abstract Background Environmental exposures in utero which modify DNA methylation may have a long-lasting impact on health and disease in offspring. We aimed to identify and replicate previously published genomic loci where DNA methylation changes are attributable to in utero exposures in the NutriGen birth cohort studies Alliance. Methods We reviewed the literature to identify differentially methylated sites of newborn DNA which are associated with the following five traits of interest maternal diabetes, pre-pregnancy body mass index (BMI), diet during pregnancy, smoking, and gestational age. We then attempted to replicate these published associations in the Canadian Healthy Infant Longitudinal Development (CHILD) and the South Asian birth cohort (START) cord blood epigenome-wide data. Results We screened 68 full-text articles and identified a total of 17 cord blood epigenome-wide association studies (EWAS) of the traits of interest. Out of the 290 CpG sites reported, 19 were identified in more than one study; all of them associated with maternal smoking. In CHILD and START EWAS, thousands of sites associated with gestational age were identified and maintained significance after correction for multiple testing. In CHILD, there was differential methylation observed for 8 of the published maternal smoking sites. No other traits tested (i.e., folate levels, gestational diabetes, birthweight) replicated in the CHILD or START cohorts. Conclusions Maternal smoking during pregnancy and gestational age are strongly associated with differential methylation in offspring cord blood, as assessed in the EWAS literature and our birth cohorts. There are a limited number of reported methylation sites associated in more than two independent studies related to pregnancy. Additional large studies of diverse populations with fine phenotyping are needed to produce robust epigenome-wide data in order to further elucidate the effect of intrauterine exposures on the infants’ methylome.

Abstract Aims The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure. Methods and results The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome‐wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow‐up following heart failure diagnosis ranged from 2 to 116 months. Forty‐nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34–90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of ≥1.10 for common variants (allele frequency ≥ 0.05) and ≥1.20 for low‐frequency variants (allele frequency 0.01–0.05) at P

Abstract Background Claes-Jensen syndrome is an X-linked inherited intellectual disability caused by mutations in the KDM5C gene. Kdm5c is a histone lysine demethylase involved in histone modifications and chromatin remodeling. Males with hemizygous mutations in KDM5C present with intellectual disability and facial dysmorphism, while most heterozygous female carriers are asymptomatic. We hypothesized that loss of Kdm5c function may influence other components of the epigenomic machinery including DNA methylation in affected patients. Results Genome-wide DNA methylation analysis of 7 male patients affected with Claes-Jensen syndrome and 56 age- and sex-matched controls identified a specific DNA methylation defect (epi-signature) in the peripheral blood of these patients, including 1769 individual CpGs and 9 genomic regions. Six healthy female carriers showed less pronounced but distinctive changes in the same regions enabling their differentiation from both patients and controls. Highly specific computational model using the most significant methylation changes demonstrated 100% accuracy in differentiating patients, carriers, and controls in the training cohort, which was confirmed on a separate cohort of patients and carriers. The 100% specificity of this unique epi-signature was further confirmed on additional 500 unaffected controls and 600 patients with intellectual disability and developmental delay, including other patient cohorts with previously described epi-signatures. Conclusion Peripheral blood epi-signature in Claes-Jensen syndrome can be used for molecular diagnosis and carrier identification and assist with interpretation of genetic variants of unknown clinical significance in the KDM5C gene.

Kabuki syndrome (KS) is caused by mutations in KMT2D, which is a histone methyltransferase involved in methylation of H3K4, a histone marker associated with DNA methylation. Analysis of >450,000 CpGs in 24 KS patients with pathogenic mutations in KMT2D and 216 controls, identified 24 genomic regions, along with 1,504 CpG sites with significant DNA methylation changes including a number of Hox genes and the MYO1F gene. Using the most differentiating and significant probes and regions we developed a “methylation variant pathogenicity (MVP) score,” which enables 100% sensitive and specific identification of individuals with KS, which was confirmed using multiple public and internal patient DNA methylation databases. We also demonstrated the ability of the MVP score to accurately reclassify variants of unknown significance in subjects with apparent clinical features of KS, enabling its potential use in molecular diagnostics. These findings provide novel insights into the molecular etiology of KS and illustrate that DNA methylation patterns can be interpreted as ‘epigenetic echoes’ in certain clinical disorders.

Phenotypic variance heterogeneity across genotypes at a single nucleotide polymorphism (SNP) may reflect underlying gene-environment (G×E) or gene-gene interactions. We modeled variance heterogeneity for blood lipids and BMI in up to 44,211 participants and investigated relationships between variance effects (Pv), G×E interaction effects (with smoking and physical activity), and marginal genetic effects (Pm). Correlations between Pv and Pm were stronger for SNPs with established marginal effects (Spearman's ρ = 0.401 for triglycerides, and ρ = 0.236 for BMI) compared to all SNPs. When Pv and Pm were compared for all pruned SNPs, only BMI was statistically significant (Spearman's ρ = 0.010). Overall, SNPs with established marginal effects were overrepresented in the nominally significant part of the Pv distribution (Pbinomial

A polygenic model of inheritance, whereby hundreds or thousands of weakly associated variants contribute to a trait's heritability, has been proposed to underlie the genetic architecture of complex traits. However, relatively few genetic variants have been positively identified so far and they collectively explain only a small fraction of the predicted heritability. We hypothesized that joint association of multiple weakly associated variants over large chromosomal regions contributes to complex traits variance. Confirmation of such regional associations can help identify new loci and lead to a better understanding of known ones. To test this hypothesis, we first characterized the ability of commonly used genetic association models to identify large region joint associations. Through theoretical derivation and simulation, we showed that multivariate linear models where multiple SNPs are included as independent predictors have the most favorable association profile. Based on these results, we tested for large region association with height in 3,740 European participants from the Health and Retirement Study (HRS) study. Adjusting for SNPs with known association with height, we demonstrated clustering of weak associations (p = 2x10-4) in regions extending up to 433.0 Kb from known height loci. The contribution of regional associations to phenotypic variance was estimated at 0.172 (95% CI 0.063-0.279; p < 0.001), which compared favorably to 0.129 explained by known height variants. Conversely, we showed that suggestively associated regions are enriched for known height loci. To extend our findings to other traits, we also tested BMI, HDLc and CRP for large region associations, with consistent results for CRP. Our results demonstrate the presence of large region joint associations and suggest these can be used to pinpoint weakly associated SNPs.

Numerous obesity loci have been identified using genome-wide association studies. A UK study indicated that physical activity may attenuate the cumulative effect of 12 of these loci, but replication studies are lacking. Therefore, we tested whether the aggregate effect of these loci is diminished in adults of European ancestry reporting high levels of physical activity. Twelve obesity-susceptibility loci were genotyped or imputed in 111,421 participants. A genetic risk score (GRS) was calculated by summing the BMI-associated alleles of each genetic variant. Physical activity was assessed using self-administered questionnaires. Multiplicative interactions between the GRS and physical activity on BMI were tested in linear and logistic regression models in each cohort, with adjustment for age, age(2), sex, study center (for multicenter studies), and the marginal terms for physical activity and the GRS. These results were combined using meta-analysis weighted by cohort sample size. The meta-analysis yielded a statistically significant GRS × physical activity interaction effect estimate (Pinteraction = 0.015). However, a statistically significant interaction effect was only apparent in North American cohorts (n = 39,810, Pinteraction = 0.014 vs. n = 71,611, Pinteraction = 0.275 for Europeans). In secondary analyses, both the FTO rs1121980 (Pinteraction = 0.003) and the SEC16B rs10913469 (Pinteraction = 0.025) variants showed evidence of SNP × physical activity interactions. This meta-analysis of 111,421 individuals provides further support for an interaction between physical activity and a GRS in obesity disposition, although these findings hinge on the inclusion of cohorts from North America, indicating that these results are either population-specific or non-causal.

Soluble ICAM-1 (sICAM-1) is an endothelium-derived inflammatory marker that has been associated with diverse conditions such as myocardial infarction, diabetes, stroke, and malaria. Despite evidence for a heritable component to sICAM-1 levels, few genetic loci have been identified so far. To comprehensively address this issue, we performed a genome-wide association analysis of sICAM-1 concentration in 22,435 apparently healthy women from the Women's Genome Health Study. While our results confirm the previously reported associations at the ABO and ICAM1 loci, four novel associations were identified in the vicinity of NFKBIK (rs3136642, P = 5.4 × 10(-9)), PNPLA3 (rs738409, P = 5.8 × 10(-9)), RELA (rs1049728, P = 2.7 × 10(-16)), and SH2B3 (rs3184504, P = 2.9 × 10(-17)). Two loci, NFKBIB and RELA, are involved in NFKB signaling pathway; PNPLA3 is known for its association with fatty liver disease; and SH3B2 has been associated with a multitude of traits and disease including myocardial infarction. These associations provide insights into the genetic regulation of sICAM-1 levels and implicate these loci in the regulation of endothelial function.

Testing for genetic effects on mean values of a quantitative trait has been a very successful strategy. However, most studies to date have not explored genetic effects on the variance of quantitative traits as a relevant consequence of genetic variation. In this report, we demonstrate that, under plausible scenarios of genetic interaction, the variance of a quantitative trait is expected to differ among the three possible genotypes of a biallelic SNP. Leveraging this observation with Levene's test of equality of variance, we propose a novel method to prioritize SNPs for subsequent gene-gene and gene-environment testing. This method has the advantageous characteristic that the interacting covariate need not be known or measured for a SNP to be prioritized. Using simulations, we show that this method has increased power over exhaustive search under certain conditions. We further investigate the utility of variance per genotype by examining data from the Women's Genome Health Study. Using this dataset, we identify new interactions between the LEPR SNP rs12753193 and body mass index in the prediction of C-reactive protein levels, between the ICAM1 SNP rs1799969 and smoking in the prediction of soluble ICAM-1 levels, and between the PNPLA3 SNP rs738409 and body mass index in the prediction of soluble ICAM-1 levels. These results demonstrate the utility of our approach and provide novel genetic insight into the relationship among obesity, smoking, and inflammation.

While conventional LDL-C, HDL-C, and triglyceride measurements reflect aggregate properties of plasma lipoprotein fractions, NMR-based measurements more accurately reflect lipoprotein particle concentrations according to class (LDL, HDL, and VLDL) and particle size (small, medium, and large). The concentrations of these lipoprotein sub-fractions may be related to risk of cardiovascular disease and related metabolic disorders. We performed a genome-wide association study of 17 lipoprotein measures determined by NMR together with LDL-C, HDL-C, triglycerides, ApoA1, and ApoB in 17,296 women from the Women's Genome Health Study (WGHS). Among 36 loci with genome-wide significance (P

Type 2 diabetes is a leading cause of morbidity and mortality. While genetic variants have been found to influence the risk of type 2 diabetes mellitus, relatively few studies have focused on genes associated with glycated hemoglobin, an index of the mean blood glucose concentration of the preceding 8-12 weeks. Epidemiologic studies and randomized clinical trials have documented the relationship between glycated hemoglobin levels and the development of long-term complications in diabetes; moreover, higher glycated hemoglobin levels in the subdiabetic range have been shown to predict type 2 diabetes risk and cardiovascular disease. To examine the common genetic determinants of glycated hemoglobin levels, we performed a genome-wide association study that evaluated 337,343 SNPs in 14,618 apparently healthy Caucasian women. The results show that glycated hemoglobin levels are associated with genetic variation at the GCK (rs730497; P = 2.8 x 10(-12)), SLC30A8 (rs13266634; P = 9.8 x 10(-8)), G6PC2 (rs1402837; P = 6.8 x 10(-10)), and HK1 (rs7072268; P = 6.4 x 10(-9)) loci. While associations at the GCK, SLC30A8, and G6PC2 loci are confirmatory, the findings at HK1 are novel. We were able to replicate this novel association in an independent validation sample of 455 additional non-diabetic men and women. HK1 encodes the enzyme hexokinase, the first step in glycolysis and a likely candidate for the control of glucose metabolism. This observed genetic association between glycated hemoglobin levels and HK1 polymorphisms paves the way for further studies of the role of HK1 in hemoglobin glycation, glucose metabolism, and diabetes.

While circulating levels of soluble Intercellular Adhesion Molecule 1 (sICAM-1) have been associated with diverse conditions including myocardial infarction, stroke, malaria, and diabetes, comprehensive analysis of the common genetic determinants of sICAM-1 is not available. In a genome-wide association study conducted among 6,578 participants in the Women's Genome Health Study, we find that three SNPs at the ICAM1 (19p13.2) locus (rs1799969, rs5498 and rs281437) are non-redundantly associated with plasma sICAM-1 concentrations at a genome-wide significance level (P