An approach to identify gene-environment interactions and reveal new biological insight in complex traits.

There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the Gene-Lifestyle Interactions Working Group within the Cohorts for Heart and Aging Research in Genetic Epidemiology Consortium has been spearheading efforts to investigate G × E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identify and confirm 5 loci (6 independent signals) interacted with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrate that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated heritability is significant (P < 0.02) for low-density lipoprotein cholesterol and triglycerides in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.

A new Approach to Identify Gene-Environment Interactions and Reveal New Biological Insight in Complex traits.

There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the CHARGE Gene-Lifestyle Interactions Working Group has been spearheading efforts to investigate G×E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identified and confirmed 5 loci (6 independent signals) interacting with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrated that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated contribution ranges from 1.76% to 14.05% of SNP heritability of serum lipids in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.

Genetic susceptibility to cognitive decline following craniospinal irradiation for pediatric central nervous system tumors.

BACKGROUND: Survivors of pediatric central nervous system (CNS) tumors treated with craniospinal irradiation (CSI) exhibit long-term cognitive difficulties. Goals of this study were to evaluate longitudinal effects of candidate and novel genetic variants on cognitive decline following CSI. METHODS: Intelligence quotient (IQ), working memory (WM), and processing speed (PS) were longitudinally collected from patients treated with CSI (n = 241). Genotype-by-time interactions were evaluated using mixed-effects linear regression to identify common variants (minor allele frequency > 1%) associated with cognitive performance change. Novel variants associated with cognitive decline (P < 5 × 10-5) in individuals of European ancestry (n = 163) were considered replicated if they demonstrated consistent genotype-by-time interactions (P < .05) in individuals of non-European ancestries (n = 78) and achieved genome-wide statistical significance (P < 5 × 10-8) in a meta-analysis across ancestry groups. RESULTS: Participants were mostly males (65%) diagnosed with embryonal tumors (98%) at a median age of 8.3 years. Overall, 1150 neurocognitive evaluations were obtained (median = 5, range: 2-10 per participant). One of the five loci previously associated with cognitive outcomes in pediatric CNS tumors survivors demonstrated significant time-dependent IQ declines (PPARA rs6008197, P = .004). Two variants associated with IQ in the general population were associated with declines in IQ after Bonferroni correction (rs9348721, P = 1.7 × 10-5; rs31771, P = 7.8 × 10-4). In genome-wide analyses, we identified novel loci associated with accelerated declines in IQ (rs116595313, meta-P = 9.4 × 10-9), WM (rs17774009, meta-P = 4.2 × 10-9), and PS (rs77467524, meta-P = 1.5 × 10-8; rs17630683, meta-P = 2.0 × 10-8; rs73249323, meta-P = 3.1 × 10-8). CONCLUSIONS: Inherited genetic variants involved in baseline cognitive functioning and novel susceptibility loci jointly influence the degree of treatment-associated cognitive decline in pediatric CNS tumor survivors.

Novel genetic regulators of fibrinogen synthesis identified by an in vitro experimental platform.

BACKGROUND: Fibrinogen has an established, essential role in both coagulation and inflammatory pathways, and these processes are deeply intertwined in the development of thrombotic and atherosclerotic diseases. Previous studies aimed to better understand the (patho) physiological actions of fibrinogen by characterizing the genomic contribution to circulating fibrinogen levels. OBJECTIVES: Establish an in vitro approach to define functional roles between genes within these loci and fibrinogen synthesis. METHODS: Candidate genes were selected on the basis of their proximity to genetic variants associated with fibrinogen levels and expression in hepatocytes and HepG2 cells. HepG2 cells were transfected with small interfering RNAs targeting candidate genes and cultured in the absence or presence of the proinflammatory cytokine interleukin-6. Effects on fibrinogen protein production, gene expression, and cell growth were assessed by immunoblotting, real-time polymerase chain reaction, and cell counts, respectively. RESULTS: HepG2 cells secreted fibrinogen, and stimulation with interleukin-6 increased fibrinogen production by 3.4 ± 1.2 fold. In the absence of interleukin-6, small interfering RNA knockdown of FGA, IL6R, or EEPD1 decreased fibrinogen production, and knockdown of LEPR, PDIA5, PLEC, SHANK3, or CPS1 increased production. In the presence of interleukin-6, knockdown of FGA, IL6R, or ATXN2L decreased fibrinogen production. Knockdown of FGA, IL6R, EEPD1, LEPR, PDIA5, PLEC, or CPS1 altered transcription of one or more fibrinogen genes. Knocking down ATXN2L suppressed inducible but not basal fibrinogen production via a post-transcriptional mechanism. CONCLUSIONS: We established an in vitro platform to define the impact of select gene products on fibrinogen production. Genes identified in our screen may reveal cellular mechanisms that drive fibrinogen production as well as fibrin(ogen)-mediated (patho)physiological mechanisms.

Multi-omics insights into the biological mechanisms underlying statistical gene-by-lifestyle interactions with smoking and alcohol consumption.

Though both genetic and lifestyle factors are known to influence cardiometabolic outcomes, less attention has been given to whether lifestyle exposures can alter the association between a genetic variant and these outcomes. The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium's Gene-Lifestyle Interactions Working Group has recently published investigations of genome-wide gene-environment interactions in large multi-ancestry meta-analyses with a focus on cigarette smoking and alcohol consumption as lifestyle factors and blood pressure and serum lipids as outcomes. Further description of the biological mechanisms underlying these statistical interactions would represent a significant advance in our understanding of gene-environment interactions, yet accessing and harmonizing individual-level genetic and 'omics data is challenging. Here, we demonstrate the coordinated use of summary-level data for gene-lifestyle interaction associations on up to 600,000 individuals, differential methylation data, and gene expression data for the characterization and prioritization of loci for future follow-up analyses. Using this approach, we identify 48 genes for which there are multiple sources of functional support for the identified gene-lifestyle interaction. We also identified five genes for which differential expression was observed by the same lifestyle factor for which a gene-lifestyle interaction was found. For instance, in gene-lifestyle interaction analysis, the T allele of rs6490056 (ALDH2) was associated with higher systolic blood pressure, and a larger effect was observed in smokers compared to non-smokers. In gene expression studies, this allele is associated with decreased expression of ALDH2, which is part of a major oxidative pathway. Other results show increased expression of ALDH2 among smokers. Oxidative stress is known to contribute to worsening blood pressure. Together these data support the hypothesis that rs6490056 reduces expression of ALDH2, which raises oxidative stress, leading to an increase in blood pressure, with a stronger effect among smokers, in whom the burden of oxidative stress is greater. Other genes for which the aggregation of data types suggest a potential mechanism include: GCNT4×current smoking (HDL), PTPRZ1×ever-smoking (HDL), SYN2×current smoking (pulse pressure), and TMEM116×ever-smoking (mean arterial pressure). This work demonstrates the utility of careful curation of summary-level data from a variety of sources to prioritize gene-lifestyle interaction loci for follow-up analyses.

Whole genome sequencing identifies structural variants contributing to hematologic traits in the NHLBI TOPMed program.

Genome-wide association studies have identified thousands of single nucleotide variants and small indels that contribute to variation in hematologic traits. While structural variants are known to cause rare blood or hematopoietic disorders, the genome-wide contribution of structural variants to quantitative blood cell trait variation is unknown. Here we utilized whole genome sequencing data in ancestrally diverse participants of the NHLBI Trans Omics for Precision Medicine program (N = 50,675) to detect structural variants associated with hematologic traits. Using single variant tests, we assessed the association of common and rare structural variants with red cell-, white cell-, and platelet-related quantitative traits and observed 21 independent signals (12 common and 9 rare) reaching genome-wide significance. The majority of these associations (N = 18) replicated in independent datasets. In genome-editing experiments, we provide evidence that a deletion associated with lower monocyte counts leads to disruption of an S1PR3 monocyte enhancer and decreased S1PR3 expression.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program.

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV1] and its ratio to forced vital capacity [FEV1/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV1 and FEV1/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Gene-lifestyle interactions in the genomics of human complex traits.

The role and biological significance of gene-environment interactions in human traits and diseases remain poorly understood. To address these questions, the CHARGE Gene-Lifestyle Interactions Working Group conducted series of genome-wide interaction studies (GWIS) involving up to 610,475 individuals across four ancestries for three lipids and four blood pressure traits, while accounting for interaction effects with drinking and smoking exposures. Here we used GWIS summary statistics from these studies to decipher potential differences in genetic associations and G×E interactions across phenotype-exposure-ancestry combinations, and to derive insights on the potential mechanistic underlying G×E through in-silico functional analyses. Our analyses show first that interaction effects likely contribute to the commonly reported ancestry-specific genetic effect in complex traits, and second, that some phenotype-exposures pairs are more likely to benefit from a greater detection power when accounting for interactions. It also highlighted modest correlation between marginal and interaction effects, providing material for future methodological development and biological discussions. We also estimated contributions to phenotypic variance, including in particular the genetic heritability conditional on the exposure, and heritability partitioned across a range of functional annotations and cell types. In these analyses, we found multiple instances of potential heterogeneity of functional partitions between exposed and unexposed individuals, providing new evidence for likely exposure-specific genetic pathways. Finally, along this work, we identified potential biases in methods used to jointly meta-analyze genetic and interaction effects. We performed simulations to characterize these limitations and to provide the community with guidelines for future G×E studies.

Rare and low-frequency exonic variants and gene-by-smoking interactions in pulmonary function.

Genome-wide association studies have identified numerous common genetic variants associated with spirometric measures of pulmonary function, including forced expiratory volume in one second (FEV1), forced vital capacity, and their ratio. However, variants with lower minor allele frequencies are less explored. We conducted a large-scale gene-smoking interaction meta-analysis on exonic rare and low-frequency variants involving 44,429 individuals of European ancestry in the discovery stage and sought replication in the UK BiLEVE study with 45,133 European ancestry samples and UK Biobank study with 59,478 samples. We leveraged data on cigarette smoking, the major environmental risk factor for reduced lung function, by testing gene-by-smoking interaction effects only and simultaneously testing the genetic main effects and interaction effects. The most statistically significant signal that replicated was a previously reported low-frequency signal in GPR126, distinct from common variant associations in this gene. Although only nominal replication was obtained for a top rare variant signal rs142935352 in one of the two studies, interaction and joint tests for current smoking and PDE3B were significantly associated with FEV1. This study investigates the utility of assessing gene-by-smoking interactions and underscores their effects on potential pulmonary function.

Supplemental Association of Clonal Hematopoiesis With Incident Heart Failure.

BACKGROUND: Age-related clonal hematopoiesis of indeterminate potential (CHIP), defined as clonally expanded leukemogenic sequence variations (particularly in DNMT3A, TET2, ASXL1, and JAK2) in asymptomatic individuals, is associated with cardiovascular events, including recurrent heart failure (HF). OBJECTIVES: This study sought to evaluate whether CHIP is associated with incident HF. METHODS: CHIP status was obtained from whole exome or genome sequencing of blood DNA in participants without prevalent HF or hematological malignancy from 5 cohorts. Cox proportional hazards models were performed within each cohort, adjusting for demographic and clinical risk factors, followed by fixed-effect meta-analyses. Large CHIP clones (defined as variant allele frequency >10%), HF with or without baseline coronary heart disease, and left ventricular ejection fraction were evaluated in secondary analyses. RESULTS: Of 56,597 individuals (59% women, mean age 58 years at baseline), 3,406 (6%) had CHIP, and 4,694 developed HF (8.3%) over up to 20 years of follow-up. CHIP was prospectively associated with a 25% increased risk of HF in meta-analysis (hazard ratio: 1.25; 95% confidence interval: 1.13-1.38) with consistent associations across cohorts. ASXL1, TET2, and JAK2 sequence variations were each associated with an increased risk of HF, whereas DNMT3A sequence variations were not associated with HF. Secondary analyses suggested large CHIP was associated with a greater risk of HF (hazard ratio: 1.29; 95% confidence interval: 1.15-1.44), and the associations for CHIP on HF with and without prior coronary heart disease were homogenous. ASXL1 sequence variations were associated with reduced left ventricular ejection fraction. CONCLUSIONS: CHIP, particularly sequence variations in ASXL1, TET2, and JAK2, represents a new risk factor for HF.

Leveraging a health information exchange for analyses of COVID-19 outcomes including an example application using smoking history and mortality.

Understanding sociodemographic, behavioral, clinical, and laboratory risk factors in patients diagnosed with COVID-19 is critically important, and requires building large and diverse COVID-19 cohorts with both retrospective information and prospective follow-up. A large Health Information Exchange (HIE) in Southeast Texas, which assembles and shares electronic health information among providers to facilitate patient care, was leveraged to identify COVID-19 patients, create a cohort, and identify risk factors for both favorable and unfavorable outcomes. The initial sample consists of 8,874 COVID-19 patients ascertained from the pandemic's onset to June 12th, 2020 and was created for the analyses shown here. We gathered demographic, lifestyle, laboratory, and clinical data from patient's encounters across the healthcare system. Tobacco use history was examined as a potential risk factor for COVID-19 fatality along with age, gender, race/ethnicity, body mass index (BMI), and number of comorbidities. Of the 8,874 patients included in the cohort, 475 died from COVID-19. Of the 5,356 patients who had information on history of tobacco use, over 26% were current or former tobacco users. Multivariable logistic regression showed that the odds of COVID-19 fatality increased among those who were older (odds ratio = 1.07, 95% CI 1.06, 1.08), male (1.91, 95% CI 1.58, 2.31), and had a history of tobacco use (2.45, 95% CI 1.93, 3.11). History of tobacco use remained significantly associated (1.65, 95% CI 1.27, 2.13) with COVID-19 fatality after adjusting for age, gender, and race/ethnicity. This effort demonstrates the impact of having an HIE to rapidly identify a cohort, aggregate sociodemographic, behavioral, clinical and laboratory data across disparate healthcare providers electronic health record (HER) systems, and follow the cohort over time. These HIE capabilities enable clinical specialists and epidemiologists to conduct outcomes analyses during the current COVID-19 pandemic and beyond. Tobacco use appears to be an important risk factor for COVID-19 related death.

Maternal Hypertension-Related Genotypes and Congenital Heart Defects.

BACKGROUND: Maternal hypertension has been associated with congenital heart defect occurrence in several studies. We assessed whether maternal genotypes associated with this condition were also associated with congenital heart defect occurrence. METHODS: We used data from the National Birth Defects Prevention Study to identify non-Hispanic white (NHW) and Hispanic women with (cases) and without (controls) a pregnancy in which a select simple, isolated heart defect was present between 1999 and 2011. We genotyped 29 hypertension-related single nucleotide polymorphisms (SNPs). We conducted logistic regression analyses separately by race/ethnicity to assess the relationship between the presence of any congenital heart defect and each SNP and an overall blood pressure genetic risk score (GRS). All analyses were then repeated to assess 4 separate congenital heart defect subtypes. RESULTS: Four hypertension-related variants were associated with congenital heart defects among NHW women (N = 1,568 with affected pregnancies). For example, 1 intronic variant in ARHGAP2, rs633185, was associated with conotruncal defects (odds ratio [OR]: 1.3, 95% confidence interval [CI]: 1.1-1.6). Additionally, 2 variants were associated with congenital heart defects among Hispanic women (N = 489 with affected pregnancies). The GRS had a significant association with septal defects (OR: 2.1, 95% CI: 1.2-3.5) among NHW women. CONCLUSIONS: We replicated a previously reported association between rs633185 and conotruncal defects. Although additional hypertension-related SNPs were also associated with congenital heart defects, more work is needed to better understand the relationship between genetic risk for maternal hypertension and congenital heart defects occurrence.

Identifying blood pressure loci whose effects are modulated by multiple lifestyle exposures.

Although multiple lifestyle exposures simultaneously impact blood pressure (BP) and cardiovascular health, most analysis so far has considered each single lifestyle exposure (e.g., smoking) at a time. Here, we exploit gene-multiple lifestyle exposure interactions to find novel BP loci. For each of 6,254 Framingham Heart Study participants, we computed lifestyle risk score (LRS) value by aggregating the risk of four lifestyle exposures (smoking, alcohol, education, and physical activity) on BP. Using the LRS, we performed genome-wide gene-environment interaction analysis in systolic and diastolic BP using the joint 2 degree of freedom (DF) and 1 DF interaction tests. We identified one genome-wide significant (p < 5 × 10-8 ) and 11 suggestive (p < 1 × 10-6 ) loci. Gene-environment analysis using single lifestyle exposures identified only one of the 12 loci. Nine of the 12 BP loci detected were novel. Loci detected by the LRS were located within or nearby genes with biologically plausible roles in the pathophysiology of hypertension, including KALRN, VIPR2, SNX1, and DAPK2. Our results suggest that simultaneous consideration of multiple lifestyle exposures in gene-environment interaction analysis can identify additional loci missed by single lifestyle approaches.

A genome-wide association study identifies new loci for factor VII and implicates factor VII in ischemic stroke etiology.

Factor VII (FVII) is an important component of the coagulation cascade. Few genetic loci regulating FVII activity and/or levels have been discovered to date. We conducted a meta-analysis of 9 genome-wide association studies of plasma FVII levels (7 FVII activity and 2 FVII antigen) among 27 495 participants of European and African ancestry. Each study performed ancestry-specific association analyses. Inverse variance weighted meta-analysis was performed within each ancestry group and then combined for a trans-ancestry meta-analysis. Our primary analysis included the 7 studies that measured FVII activity, and a secondary analysis included all 9 studies. We provided functional genomic validation for newly identified significant loci by silencing candidate genes in a human liver cell line (HuH7) using small-interfering RNA and then measuring F7 messenger RNA and FVII protein expression. Lastly, we used meta-analysis results to perform Mendelian randomization analysis to estimate the causal effect of FVII activity on coronary artery disease, ischemic stroke (IS), and venous thromboembolism. We identified 2 novel (REEP3 and JAZF1-AS1) and 6 known loci associated with FVII activity, explaining 19.0% of the phenotypic variance. Adding FVII antigen data to the meta-analysis did not result in the discovery of further loci. Silencing REEP3 in HuH7 cells upregulated FVII, whereas silencing JAZF1 downregulated FVII. Mendelian randomization analyses suggest that FVII activity has a positive causal effect on the risk of IS. Variants at REEP3 and JAZF1 contribute to FVII activity by regulating F7 expression levels. FVII activity appears to contribute to the etiology of IS in the general population.

RNA sequence analysis of inducible pluripotent stem cell-derived cardiomyocytes reveals altered expression of DNA damage and cell cycle genes in response to doxorubicin.

Anthracyclines, such as doxorubicin, are highly effective chemotherapeutic agents, yet are associated with increased risk of cardiotoxicity. The genes and pathways involved in the development of heart damage following doxorubicin exposure in humans remain elusive. Our objective was to explore time- and dose-dependent changes in gene expression via RNA sequence (RNAseq) that mediate doxorubicin response in human iPSC-cardiomyocytes following 50, 150, or 450 nM exposure for 2, 7, or 12 days. Clustering and differential expression analyses were conducted to identify genes with altered expression. Samples clustered in dose and time-dependent manners, and MCM5, PRC1, NUSAP1, CENPF, CCNB1, MELK, AURKB, and RACGAP1 were consistently significantly differentially expressed between untreated and treated conditions. These genes were also significantly downregulated in pairwise analyses, which was validated by reverse transcription polymerase chain reaction (RT-PCR). Pathway analysis identified the top canonical pathways involved in response, implicating DNA damage repair response and the cell cycle as having roles in the development of doxorubicin-induced cardiotoxicity in the human cardiomyocyte.

CD163+ macrophages promote angiogenesis and vascular permeability accompanied by inflammation in atherosclerosis.

Intake of hemoglobin by the hemoglobin-haptoglobin receptor CD163 leads to a distinct alternative non-foam cell antiinflammatory macrophage phenotype that was previously considered atheroprotective. Here, we reveal an unexpected but important pathogenic role for these macrophages in atherosclerosis. Using human atherosclerotic samples, cultured cells, and a mouse model of advanced atherosclerosis, we investigated the role of intraplaque hemorrhage on macrophage function with respect to angiogenesis, vascular permeability, inflammation, and plaque progression. In human atherosclerotic lesions, CD163+ macrophages were associated with plaque progression, microvascularity, and a high level of HIF1α and VEGF-A expression. We observed irregular vascular endothelial cadherin in intraplaque microvessels surrounded by CD163+ macrophages. Within these cells, activation of HIF1α via inhibition of prolyl hydroxylases promoted VEGF-mediated increases in intraplaque angiogenesis, vascular permeability, and inflammatory cell recruitment. CD163+ macrophages increased intraplaque endothelial VCAM expression and plaque inflammation. Subjects with homozygous minor alleles of the SNP rs7136716 had elevated microvessel density, increased expression of CD163 in ruptured coronary plaques, and a higher risk of myocardial infarction and coronary heart disease in population cohorts. Thus, our findings highlight a nonlipid-driven mechanism by which alternative macrophages promote plaque angiogenesis, leakiness, inflammation, and progression via the CD163/HIF1α/VEGF-A pathway.

GWAS for male-pattern baldness identifies 71 susceptibility loci explaining 38% of the risk.

Male pattern baldness (MPB) or androgenetic alopecia is one of the most common conditions affecting men, reaching a prevalence of ~50% by the age of 50; however, the known genes explain little of the heritability. Here, we present the results of a genome-wide association study including more than 70,000 men, identifying 71 independently replicated loci, of which 30 are novel. These loci explain 38% of the risk, suggesting that MPB is less genetically complex than other complex traits. We show that many of these loci contain genes that are relevant to the pathology and highlight pathways and functions underlying baldness. Finally, despite only showing genome-wide genetic correlation with height, pathway-specific genetic correlations are significant for traits including lifespan and cancer. Our study not only greatly increases the number of MPB loci, illuminating the genetic architecture, but also provides a new approach to disentangling the shared biological pathways underlying complex diseases.

Mutations in folate transporter genes and risk for human myelomeningocele.

The molecular mechanisms linking folate deficiency and neural tube defect (NTD) risk in offspring remain unclear. Folate transporters (SLC19A1, SLC46A1, SLC25A32, and FOLH1) and folate receptors (FOLR1, FOLR2, and FOLR3) are suggested to play essential roles in transporting folate from maternal intestinal lumen to the developing embryo. Loss of function variants in these genes may affect folate availability and contribute to NTD risk. This study examines whether variants within the folate transporter and receptor genes are associated with an increased risk for myelomeningocele (MM). Exons and their flanking intron sequences of 348 MM subjects were sequenced using the Sanger sequencing method and/or next generation sequencing to identify variants. Frequencies of alleles of single nucleotide polymorphisms (SNPs) in MM subjects were compared to those from ethnically matched reference populations to evaluate alleles' associated risk for MM. We identified eight novel variants in SLC19A1 and twelve novel variants in FOLR1, FOLR2, and FOLR3. Pathogenic variants include c.1265delG in SLC19A1 resulting in an early stop codon, four large insertion deletion variants in FOLR3, and a stop_gain variant in FOLR3. No new variants were identified in SLC46A1, SLC25A32, or FOLH1. In SLC19A1, c.80A>G (rs1051266) was not associated with our MM cohort; we did observe a variant allele G frequency of 61.7%, higher than previously reported in other NTD populations. In conclusion, we discovered novel loss of function variants in genes involved in folate transport in MM subjects. Our results support the growing evidence of associations between genes involved in folate transport and susceptibility to NTDs.

Whole-genome sequencing study of serum peptide levels: the Atherosclerosis Risk in Communities study.

Oligopeptides are important markers of protein metabolism, as they are cleaved from larger polypeptides and proteins. Genetic association studies may help elucidate their origin and function. In 1,552 European Americans and 1,872 African Americans of the Atherosclerosis Risk in Communities study, we performed whole-genome and whole-exome sequencing and measured serum levels of 25 peptides. Common variants (minor allele frequency > 5%) were analysed individually. We grouped low-frequency variants (minor allele frequency ≤ 5%) by a genome-wide sliding window using region-based aggregate tests. Furthermore, low-frequency regulatory variants were grouped by gene, as were functional coding variants. All analyses were performed separately in each ancestry group and then meta-analysed. We identified 22 common variant associations with peptide levels (P-value < 4.2 × 10-10), including 16 novel gene-peptide pairs. Notably, variants in kinin-kallikrein genes KNG1, F12, KLKB1, and ACE were associated with several different peptides. Variants in KLKB1 and ACE were associated with a fragment of complement component 3f. Both common variants and low-frequency coding variants in CPN1 were associated with a fibrinogen cleavage peptide. Four sliding windows were significantly associated with peptide levels (P-value < 4.2 × 10-10). Our results highlight the importance of the kinin-kallikrein system in the regulation of serum peptide levels, strengthen the evidence for a broad link between the kinin-kallikrein and complement systems, and suggest a role of CPN1 in the conversion of fibrinogen to fibrin.

Hypertension Susceptibility Loci are Associated with Anthracycline-related Cardiotoxicity in Long-term Childhood Cancer Survivors.

Anthracycline-based chemotherapy is associated with dose-dependent, irreversible damage to the heart. Childhood cancer survivors with hypertension after anthracycline exposure are at increased risk of cardiotoxicity, leading to the hypothesis that genetic susceptibility loci for hypertension may serve as predictors for development of late cardiotoxicity. Therefore, we determined the association between 12 GWAS-identified hypertension-susceptibility loci and cardiotoxicity in a cohort of long-term childhood cancer survivors (N = 108) who received anthracyclines and were screened for cardiac function via echocardiograms. Hypertension-susceptibility alleles of PLCE1:rs9327264 and ATP2B1:rs17249754 were significantly associated with cardiotoxicity risk conferring a protective effect with a 64% (95% CI: 0.18-0.76, P = 0.0068) and 74% (95% CI: 0.07-0.96, P = 0.040) reduction in risk, respectively. In RNAseq experiments of human induced pluripotent stem cell (iPSC) derived cardiomyocytes treated with doxorubicin, both PLCE1 and ATP2B1 displayed anthracycline-dependent gene expression profiles. In silico functional assessment further supported this relationship - rs9327264 in PLCE1 (P = 0.0080) and ATP2B1 expression (P = 0.0079) were both significantly associated with daunorubicin IC50 values in a panel of lymphoblastoid cell lines. Our findings demonstrate that the hypertension-susceptibility variants in PLCE1 and ATP2B1 confer a protective effect on risk of developing anthracycline-related cardiotoxicity, and functional analyses suggest that these genes are influenced by exposure to anthracyclines.