Maternal arsenic exposure modifies associations between arsenic, folate and arsenic metabolism gene variants, and spina bifida risk: A case‒control study in Bangladesh.

BACKGROUND: Spina bifida is a type of neural tube defect (NTD); NTDs are developmental malformations of the spinal cord that result from failure of neural tube closure during embryogenesis and are likely caused by interactions between genetic and environmental factors. Arsenic induces NTDs in animal models, and studies demonstrate that mice with genetic defects related to folate metabolism are more susceptible to arsenic's effects. We sought to determine whether 25 single-nucleotide polymorphisms (SNPs) in genes involved in folate and arsenic metabolism modified the associations between maternal arsenic exposure and risk of spina bifida (a common NTD) among a hospital-based case-control study population in Bangladesh. METHODS: We used data from 262 mothers and 220 infants who participated in a case‒control study at the National Institutes of Neurosciences & Hospital and Dhaka Shishu Hospital in Dhaka, Bangladesh. Neurosurgeons assessed infants using physical examinations, review of imaging, and we collected histories using questionnaires. We assessed arsenic from mothers' toenails using inductively coupled plasma mass spectrometry (ICP-MS), and we genotyped participants using the Illumina Global Screening Array v1.0. We chose candidate genes and SNPs through a review of the literature. We assessed SNP-environment interactions using interaction terms and stratified models, and we assessed gene-environment interactions using interaction sequence/SNP-set kernel association tests (iSKAT). RESULTS: The median toenail arsenic concentration was 0.42 µg/g (interquartile range [IQR]: 0.27-0.86) among mothers of cases and 0.47 µg/g (IQR: 0.30-0.97) among mothers of controls. We found an two SNPs in the infants' AS3MT gene (rs11191454 and rs7085104) and one SNP in mothers' DNMT1 gene (rs2228611) were associated with increased odds of spina bifida in the setting of high arsenic exposure (rs11191454, OR 3.01, 95% CI: 1.28-7.09; rs7085104, OR 2.33, 95% CI: 1.20-4.and rs2228611, OR 2.11, 95% CI: 1.11-4.01), along with significant SNP-arsenic interactions. iSKAT analyses revealed significant interactions between mothers' toenail concentrations and infants' AS3MT and MTR genes (p = 0.02), and mothers' CBS gene (p = 0.05). CONCLUSIONS: Our results support the hypothesis that arsenic increases spina bifida risk via interactions with folate and arsenic metabolic pathways and suggests that individuals in the population who have certain genetic polymorphisms in genes involved with arsenic and folate metabolism may be more susceptible than others to the arsenic teratogenicity.

TOLLIP and MUC5B modulate the effect of ambient NO2 on respiratory symptoms in infancy.

BACKGROUND: Current knowledge suggests that the gene region containing MUC5B and TOLLIP plays a role in airway defence and airway inflammation, and hence respiratory disease. It is also known that exposure to air pollution increases susceptibility to respiratory disease. We aimed to study whether the effect of air pollutants on the immune response and respiratory symptoms in infants may be modified by polymorphisms in MUC5B and TOLLIP genes. METHODS: 359 healthy term infants from the prospective Basel-Bern Infant Lung Development (BILD) birth cohort were included in the study. The main outcome was the score of weekly assessed respiratory symptoms in the first year of life. Using the candidate gene approach, we selected 10 single nucleotide polymorphisms (SNPs) from the MUC5B and TOLLIP regions. Nitrogen dioxide (NO2) and particulate matter ≤10 µm in aerodynamic diameter (PM10) exposure was estimated on a weekly basis. We used generalised additive mixed models adjusted for known covariates. To validate our results in vitro, cells from a lung epithelial cell line were downregulated in TOLLIP expression and exposed to diesel particulate matter (DPM) and polyinosinic-polycytidylic acid. RESULTS: Significant interaction was observed between modelled air pollution (weekly NO2 exposure) and 5 SNPs within MUC5B and TOLLIP genes regarding respiratory symptoms as outcome: E.g., infants carrying minor alleles of rs5744034, rs3793965 and rs3750920 (all TOLLIP) had an increased risk of respiratory symptoms with increasing NO2 exposure. In vitro experiments showed that cells downregulated for TOLLIP react differently to environmental pollutant exposure with DPM and viral stimulation. CONCLUSION: Our findings suggest that the effect of air pollution on respiratory symptoms in infancy may be influenced by the genotype of specific SNPs from the MUC5B and TOLLIP regions. For validation of the findings, we provided in vitro evidence for the interaction of TOLLIP with air pollution.

Interaction of Th17 differentiations-related gene polymorphisms and environmental factors contributing to the disease classification, complications, and surgical risks of Crohn's disease in the Chinese Han population.

OBJECTIVES: Few studies have been conducted on gene-environment interactions in the Chinese population with Crohn's disease (CD). We aimed to investigate the association between single nucleotide polymorphisms (SNPs) on the T helper 17 (Th17) cell and CD susceptibility/performance in Chinese individuals. METHODS: We conducted a case-control and case-only study at the Peking Union Medical College Hospital. Four SNPs related to the Th17 cell pathway genes were prioritized, including rs2284553 (interferon gamma receptor 2), rs7517847 (interleukin 23 receptor), rs7773324 (interferon regulatory factor 4), and rs4263839 (tumor necrosis factor superfamily 15). SNP frequency was calculated, and gene-environment interaction was assessed by multifactor dimensionality reduction analysis. RESULTS: Altogether 159 CD patients and 316 healthy controls were included. All analyzed SNPs were found in Hardy-Weinberg equilibrium (P > 0.05). The frequency of rs2284553-A allele and rs4263839-A allele were lower in CD patients compared with controls (P < 0.05). While the rs4263839-A allele was more prevalent in ileocolonic CD patients than in those with isolated small intestinal or colonic disease (P = 0.035). Gene-environment interactions revealed associations between rs2284553 and breastfeeding, sunshine exposure, and fridge-stored food, affecting age at diagnosis, intestinal involvement, and intestinal stricture. Interaction of rs4263839 and breastfeeding influenced small intestinal lesions and intestinal stricture in CD. CONCLUSIONS: This study provided information on the genetic background in Chinese CD patients. Incorporating these SNPs into predictive models may improve risk assessment and outcome prediction. Gene-environment interaction contributes to the understanding of CD pathogenesis.

Association of genetic risk, lifestyle, and their interaction with obesity and obesity-related morbidities.

The extent to which modifiable lifestyle factors offset the determined genetic risk of obesity and obesity-related morbidities remains unknown. We explored how the interaction between genetic and lifestyle factors influences the risk of obesity and obesity-related morbidities. The polygenic score for body mass index was calculated to quantify inherited susceptibility to obesity in 338,645 UK Biobank European participants, and a composite lifestyle score was derived from five obesogenic factors (physical activity, diet, sedentary behavior, alcohol consumption, and sleep duration). We observed significant interaction between high genetic risk and poor lifestyles (pinteraction < 0.001). Absolute differences in obesity risk between those who adhere to healthy lifestyles and those who do not had gradually expanded with an increase in polygenic score. Despite a high genetic risk for obesity, individuals can prevent obesity-related morbidities by adhering to a healthy lifestyle and maintaining a normal body weight. Healthy lifestyles should be promoted irrespective of genetic background.

BHAFT: Bayesian heredity-constrained accelerated failure time models for detecting gene-environment interactions in survival analysis.

In addition to considering the main effects, understanding gene-environment (G × E) interactions is imperative for determining the etiology of diseases and the factors that affect their prognosis. In the existing statistical framework for censored survival outcomes, there are several challenges in detecting G × E interactions, such as handling high-dimensional omics data, diverse environmental factors, and algorithmic complications in survival analysis. The effect heredity principle has widely been used in studies involving interaction identification because it incorporates the dependence of the main and interaction effects. However, Bayesian survival models that incorporate the assumption of this principle have not been developed. Therefore, we propose Bayesian heredity-constrained accelerated failure time (BHAFT) models for identifying main and interaction (M-I) effects with novel spike-and-slab or regularized horseshoe priors to incorporate the assumption of effect heredity principle. The R package rstan was used to fit the proposed models. Extensive simulations demonstrated that BHAFT models had outperformed other existing models in terms of signal identification, coefficient estimation, and prognosis prediction. Biologically plausible G × E interactions associated with the prognosis of lung adenocarcinoma were identified using our proposed model. Notably, BHAFT models incorporating the effect heredity principle could identify both main and interaction effects, which are highly useful in exploring G × E interactions in high-dimensional survival analysis. The code and data used in our paper are available at https://github.com/SunNa-bayesian/BHAFT.

LDER-GE estimates phenotypic variance component of gene-environment interactions in human complex traits accurately with GE interaction summary statistics and full LD information.

Gene-environment (GE) interactions are essential in understanding human complex traits. Identifying these interactions is necessary for deciphering the biological basis of such traits. In this study, we review state-of-art methods for estimating the proportion of phenotypic variance explained by genome-wide GE interactions and introduce a novel statistical method Linkage-Disequilibrium Eigenvalue Regression for Gene-Environment interactions (LDER-GE). LDER-GE improves the accuracy of estimating the phenotypic variance component explained by genome-wide GE interactions using large-scale biobank association summary statistics. LDER-GE leverages the complete Linkage Disequilibrium (LD) matrix, as opposed to only the diagonal squared LD matrix utilized by LDSC (Linkage Disequilibrium Score)-based methods. Our extensive simulation studies demonstrate that LDER-GE performs better than LDSC-based approaches by enhancing statistical efficiency by ~23%. This improvement is equivalent to a sample size increase of around 51%. Additionally, LDER-GE effectively controls type-I error rate and produces unbiased results. We conducted an analysis using UK Biobank data, comprising 307 259 unrelated European-Ancestry subjects and 966 766 variants, across 217 environmental covariate-phenotype (E-Y) pairs. LDER-GE identified 34 significant E-Y pairs while LDSC-based method only identified 23 significant E-Y pairs with 22 overlapped with LDER-GE. Furthermore, we employed LDER-GE to estimate the aggregated variance component attributed to multiple GE interactions, leading to an increase in the explained phenotypic variance with GE interactions compared to considering main genetic effects only. Our results suggest the importance of impacts of GE interactions on human complex traits.

Parental genetic effects on the offspring`s phenotype without direct transmission of the parental gene itself - pathophysiology and clinical evidence.

Parental genes can influence the phenotype of their offspring through genomic-epigenomic interactions even without the direct inheritance of specific parental genotypes. Maternal genetic variations can affect the ovarian and intrauterine environments and potentially alter lactation behaviors, impacting offspring nutrition and health outcomes independently of the fetal genome. Similarly, paternal genetic changes can affect the endocrine system and vascular functions in the testes, influencing sperm quality and seminal fluid composition. These changes can initiate early epigenetic modifications in sperm, including alterations in microRNAs, tRNA-derived small RNAs, and DNA methylation patterns. These epigenetic modifications might induce further changes in target organs of the offspring, leading to modified gene expression and phenotypic outcomes without transmitting the original parental genetic alterations. This review presents clinical evidence supporting this hypothesis and discusses the potential underlying molecular mechanisms. Parental gene-offspring epigenome-offspring phenotype interactions have been observed in neurocognitive disorders as well as cardio-renal diseases.

A scalable and robust variance components method reveals insights into the architecture of gene-environment interactions underlying complex traits.

Understanding the contribution of gene-environment interactions (GxE) to complex trait variation can provide insights into disease mechanisms, explain sources of heritability, and improve genetic risk prediction. While large biobanks with genetic and deep phenotypic data hold promise for obtaining novel insights into GxE, our understanding of GxE architecture in complex traits remains limited. We introduce a method to estimate the proportion of trait variance explained by GxE (GxE heritability) and additive genetic effects (additive heritability) across the genome and within specific genomic annotations. We show that our method is accurate in simulations and computationally efficient for biobank-scale datasets. We applied our method to common array SNPs (MAF ≥1%), fifty quantitative traits, and four environmental variables (smoking, sex, age, and statin usage) in unrelated white British individuals in the UK Biobank. We found 68 trait-E pairs with significant genome-wide GxE heritability (p<0.05/200) with a ratio of GxE to additive heritability of ≈6.8% on average. Analyzing ≈8 million imputed SNPs (MAF ≥0.1%), we documented an approximate 28% increase in genome-wide GxE heritability compared to array SNPs. We partitioned GxE heritability across minor allele frequency (MAF) and local linkage disequilibrium (LD) values, revealing that, like additive allelic effects, GxE allelic effects tend to increase with decreasing MAF and LD. Analyzing GxE heritability near genes highly expressed in specific tissues, we find significant brain-specific enrichment for body mass index (BMI) and basal metabolic rate in the context of smoking and adipose-specific enrichment for waist-hip ratio (WHR) in the context of sex.

Molecular Profiling and the Interaction of Somatic Mutations with Transcriptomic Profiles in Non-Melanoma Skin Cancer (NMSC) in a Population Exposed to Arsenic.

Exposure to inorganic arsenic (As) is recognized as a risk factor for non-melanoma skin cancer (NMSC). We followed up with 7000 adults for 6 years who were exposed to As. During follow-up, 2.2% of the males and 1.3% of the females developed basal cell carcinoma (BCC), while 0.4% of the male and 0.2% of the female participants developed squamous cell carcinoma (SCC). Using a panel of more than 400 cancer-related genes, we detected somatic mutations (SMs) in the first 32 NMSC samples (BCC = 26 and SCC = 6) by comparing paired (tissue-blood) samples from the same individual and then comparing them to the SM in healthy skin tissue from 16 participants. We identified (a) a list of NMSC-associated SMs, (b) SMs present in both NMSC and healthy skin, and (c) SMs found only in healthy skin. We also demonstrate that the presence of non-synonymous SMs in the top mutated genes (like PTCH1, NOTCH1, SYNE1, PKHD1 in BCC and TP53 in SCC) significantly affects the magnitude of differential expressions of major genes and gene pathways (basal cell carcinoma pathways, NOTCH signaling, IL-17 signaling, p53 signaling, Wnt signaling pathway). These findings may help select groups of patients for targeted therapy, like hedgehog signaling inhibitors, IL17 inhibitors, etc., in the future.

Interaction between Risk Single-Nucleotide Polymorphisms of Developmental Dyslexia and Parental Education on Reading Ability: Evidence for Differential Susceptibility Theory.

While genetic and environmental factors have been shown as predictors of children's reading ability, the interaction effects of identified genetic risk susceptibility and the specified environment for reading ability have rarely been investigated. The current study assessed potential gene-environment (G×E) interactions on reading ability in 1477 school-aged children. The gene-environment interactions on character recognition were investigated by an exploratory analysis between the risk single-nucleotide polymorphisms (SNPs), which were discovered by previous genome-wide association studies of developmental dyslexia (DD), and parental education (PE). The re-parameterized regression analysis suggested that this G×E interaction conformed to the strong differential susceptibility model. The results showed that rs281238 exhibits a significant interaction with PE on character recognition. Children with the "T" genotype profited from high PE, whereas they performed worse in low PE environments, but "CC" genotype children were not malleable in different PE environments. This study provided initial evidence for how the significant SNPs in developmental dyslexia GWA studies affect children's reading performance by interacting with the environmental factor of parental education.

Associations between residential environments and late-onset schizophrenia in UK Biobank: Interaction with genetic risk factor.

BACKGROUND: Environment and genes both contribute to schizophrenia. However, the impact of different natural environments surrounding residential addresses on schizophrenia in urban settings remains unknown. This study aimed to investigate the association of urbanisation, measured by residential environments, with late-onset schizophrenia and explore whether genetic risk for schizophrenia modified the associations. METHODS: We examined the associations between residential environments and late-onset schizophrenia and its interaction with genetic risk factors in UK Biobank, followed from 2006 to 2010 (baseline) to Dec 2021. Residential environments, including greenspace, domestic garden, blue space, and total natural environment, were evaluated using land use coverage percentage. The polygenic risk score (PRS) of schizophrenia was derived using a Bayesian approach and adjusted it against ancestry. Cox proportional hazard regression model was used to assess the associations between per interquartile (IQR) increase of each type of residential environments and late-onset schizophrenia. Interactive effects of PRS and residential environments on late-onset schizophrenia were assessed on both additive and multiplicative scales. RESULTS: A total of 393,680 participants were included in the analysis, with 844 cases of late-onset schizophrenia being observed after 12.8 years of follow-up. Within 300 m buffer surrounding the residential addresses, per interquartile increase in greenspace (31.5 %) and total natural environment (34.4 %) were both associated with an 11 % (HR = 0.89, 95 % CI 0.80, 0.99) lower risk of late-onset schizophrenia. Domestic garden and blue space did not show significant protective effects on late-onset schizophrenia. A strong dose-response relationship between schizophrenia PRS and schizophrenia was found, while no additive or multiplicative interaction effects were present between residential environments and PRS on late-onset schizophrenia. CONCLUSION: Residential greenspace and total natural environment may protect against late-onset schizophrenia in older people regardless of genetic risk. These findings shed light on the prevention of schizophrenia and urban planning to optimise ecosystem benefits linked to schizophrenia.

Unveiling the Influence of Copy Number Variations on Genetic Diversity and Adaptive Evolution in China's Native Pig Breeds via Whole-Genome Resequencing.

Copy number variations (CNVs) critically influence individual genetic diversity and phenotypic traits. In this study, we employed whole-genome resequencing technology to conduct an in-depth analysis of 50 pigs from five local swine populations [Rongchang pig (RC), Wuzhishan pig (WZS), Tibetan pig (T), Yorkshire (YL) and Landrace (LR)], aiming to assess their genetic potential and explore their prospects in the field of animal model applications. We identified a total of 96,466 CNVs, which were subsequently integrated into 7112 non-redundant CNVRs, encompassing 1.3% of the swine genome. Functional enrichment analysis of the genes within these CNVRs revealed significant associations with sensory perception, energy metabolism, and neural-related pathways. Further selective scan analyses of the local pig breeds RC, T, WZS, along with YL and LR, uncovered that for the RC variety, the genes PLA2G10 and ABCA8 were found to be closely related to fat metabolism and cardiovascular health. In the T breed, the genes NCF2 and CSGALNACT1 were associated with immune response and connective tissue characteristics. As for the WZS breed, the genes PLIN4 and CPB2 were primarily linked to fat storage and anti-inflammatory responses. In summary, this research underscores the pivotal role of CNVs in fostering the diversity and adaptive evolution of pig breeds while also offering valuable insights for further exploration of the advantageous genetic traits inherent to China's local pig breeds. This facilitates the creation of experimental animal models tailored to the specific characteristics of these breeds, contributing to the advancement of livestock and biomedical research.

The common p.Ile291Val variant of ERLIN1 enhances TM6SF2 function and is associated with protection against MASLD.

BACKGROUND: The ERLIN1 p.Ile291Val single-nucleotide polymorphism (rs2862954) is associated with protection from steatotic liver disease (SLD), but effects of this variant on metabolic phenotypes remain uncertain. METHODS: Metabolic phenotypes and outcomes associated with ERLIN1 p.Ile291Val were analyzed by using a genome-first approach in the UK Biobank (UKB), Penn Medicine BioBank (PMBB), and All of Us cohort. FINDINGS: ERLIN1 p.Ile291Val carriers exhibited significantly lower serum levels of alanine aminotransferase and aspartate aminotransferase as well as higher levels of triglycerides, low-density lipoprotein cholesterol, Apolipoprotein B, high-density lipoprotein cholesterol, and Apolipoprotein A1 in UKB, and these values were affected by ERLIN1 p.Ile291Val in an allele-dose-dependent manner. Homozygous ERLIN1 p.Ile291Val carriers had a significantly reduced risk of developing metabolic dysfunction-associated SLD (MASLD, adjusted odds ratio [aOR] = 0.92, 95% confidence interval [CI], 0.88-0.96). The protective effect of this variant was enhanced in patients with alcoholic liver disease. Our results were replicated in PMBB and the All of Us cohort. Strikingly, the protective effects of ERLIN1 p.Ile291Val were not apparent in individuals carrying the TM6SF2 p.Glu167Lys variant associated with increased risk of SLD. We analyzed the effects of predicted loss-of-function ERLIN1 variants and found that they had opposite effects, namely reduced plasma lipids, suggesting that ERLIN1 p.Ile291Val may be a gain-of-function variant. CONCLUSION: Our study contributes to a better understanding of ERLIN1 by investigating a coding variant that has emerged as a potential gain-of-function mutation with protective effects against MASLD development.

Disadvantaged social status contributed to sleep disorders: An observational and genome-wide gene-environment interaction analysis.

BACKGROUND: Sleep is a natural and essential physiological need for individuals. Our study aimed to research the associations between accumulated social risks and sleep disorders. METHODS: In this study, we came up with a polysocial risk score (PsRS), which is a cumulative social risk index composed of 13 social determinants of health. This research includes 239,165 individuals with sleep disorders and social determinants of health data from the UK Biobank cohort. First, logistic regression models were performed to examine the associations of social determinants of health and sleep disorders, including chronotype, narcolepsy, insomnia, snoring, short and long sleep duration. Then, PsRS was calculated based on statistically significant social determinants of health for each sleep disorder. Third, a genome-wide gene-environment interaction study was conducted to explore the interactions between single-nucleotide polymorphisms and PsRS in relation to sleep disorders. RESULTS: Higher PsRS scores were associated with worse sleep status, with the adjusted odds ratio (OR) ranging from 1.10 (95% Confidence interval [CI]: 1.09-1.11) to 1.29 (95% CI: 1.27-1.30) for sleep disorders. Emotional stress (OR = 1.36, 95% CI: 1.28-1.43) and not in paid employment (OR = 2.62, 95% CI: 2.51-2.74) were found to have significant contributions for sleep disorders. Moreover, multiple single-nucleotide polymorphisms were discovered to have interactions with PsRS, such as FRAS1 (P = 2.57 × 10-14) and CACNA1A (P = 8.62 × 10-14) for narcolepsy, and ACKR3 (P = 1.24 × 10-8) for long sleep. CONCLUSIONS: Our findings suggested that cumulative social risks was associated with sleep disorders, while the interactions between genetic susceptibility and disadvantaged social status are risk factors for the development of sleep disorders.

Genetic propensity to depression and the role of partnership status.

Social relationships and genetic propensity are known to affect depression risk, but their joint effects are poorly understood. This study examined the association of a polygenic index for depression with time to antidepressant (AD) purchasing and the moderating role of partnership status. We analysed data from 30,192 Finnish individuals who participated in the FINRISK and Health 2000 and 2011 surveys and had register and medication data available. We measured genetic risk with a polygenic index (PGI) for depression. Depression was assessed through antidepressant purchases. We estimated an accelerated failure time model with partnership status as time-varying and different sets of confounder adjustments. The predicted cumulative hazard of antidepressant purchasing varied across PGI and partnership status. At follow-up year 10, being widowed was associated with the largest cumulative hazard of 0.34 (95%CI: 0.28-0.39) in the 80th and 0.20 (95%CI: 0.17-0.23) in the 20th PGI percentile, followed by divorced, single, married and cohabiting. Cohabiting was associated with a cumulative hazard of 0.19 (95%CI: 0.16-0.23) in the 80th and 0.11 (95%CI: 0.1-0.13) in the 20th PGI percentile. We found no evidence for an interaction between the PGI and partnership status. Results were robust to different model specifications, gender stratification, and the choice of PGI. Although antidepressant purchasing correlated with both PGI and partnership status, we found no evidence that partnership status could partially offset or amplify the association between the PGI for depression and antidepressant purchasing incidence.

The Association of Urinary Sodium Excretion with Glaucoma and Related Traits in a Large United Kingdom Population.

PURPOSE: Excessive dietary sodium intake has known adverse effects on intravascular fluid volume and systemic blood pressure, which may influence intraocular pressure (IOP) and glaucoma risk. This study aimed to assess the association of urinary sodium excretion, a biomarker of dietary intake, with glaucoma and related traits, and determine whether this relationship is modified by genetic susceptibility to disease. DESIGN: Cross-sectional observational and gene-environment interaction analyses in the population-based UK Biobank study. PARTICIPANTS: Up to 103 634 individuals (mean age: 57 years; 51% women) with complete urinary, ocular, and covariable data. METHODS: Urine sodium:creatinine ratio (UNa:Cr; mmol:mmol) was calculated from a midstream urine sample. Ocular parameters were measured as part of a comprehensive eye examination, and glaucoma case ascertainment was through a combination of self-report and linked national hospital records. Genetic susceptibility to glaucoma was calculated based on a glaucoma polygenic risk score comprising 2673 common genetic variants. Multivariable linear and logistic regression, adjusted for key sociodemographic, medical, anthropometric, and lifestyle factors, were used to model associations and gene-environment interactions. MAIN OUTCOME MEASURES: Corneal-compensated IOP, OCT derived macular retinal nerve fiber layer and ganglion cell-inner plexiform layer (GCIPL) thickness, and prevalent glaucoma. RESULTS: In maximally adjusted regression models, a 1 standard deviation increase in UNa:Cr was associated with higher IOP (0.14 mmHg; 95% confidence interval [CI], 0.12-0.17; P < 0.001) and greater prevalence of glaucoma (odds ratio, 1.11; 95% CI, 1.07-1.14; P < 0.001) but not macular retinal nerve fiber layer or ganglion cell-inner plexiform layer thickness. Compared with those with UNa:Cr in the lowest quintile, those in the highest quintile had significantly higher IOP (0.45 mmHg; 95% CI, 0.36-0.53, P < 0.001) and prevalence of glaucoma (odds ratio, 1.30; 95% CI, 1.17-1.45; P < 0.001). Stronger associations with glaucoma (P interaction = 0.001) were noted in participants with a higher glaucoma polygenic risk score. CONCLUSIONS: Urinary sodium excretion, a biomarker of dietary intake, may represent an important modifiable risk factor for glaucoma, especially in individuals at high underlying genetic risk. These findings warrant further investigation because they may have important clinical and public health implications. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Genome-wide investigation of exogenous female hormones, genetic variation, and venous thromboembolism risk.

BACKGROUND: Increased risk of venous thromboembolism (VTE) is a life-threatening side effect for users of oral contraceptives (OCs) or hormone therapy (HT). OBJECTIVES: To investigate the potential for genetic predisposition to VTE in OC or HT users, we conducted a gene-by-environment case-only meta-analysis of genome-wide association studies (GWAS). METHODS: Use or nonuse of OCs (7 studies) or HT (8 studies) at the time of the VTE event was determined by pharmacy records or self-report. A synergy index (SI) was modeled for each variant in each study and submultiplicative/supramultiplicative gene-by-environment interactions were estimated. The SI parameters were first meta-analyzed across OC and HT studies and subsequently meta-analyzed to obtain an overall estimate. The primary analysis was agnostic GWAS and interrogated all imputed genotypes using a P value threshold of <5.0 × 10-8; secondary analyses were candidate-based. RESULTS: The VTE case-only OC meta-analysis included 2895 OC users and 6607 nonusers; the case-only HT meta-analysis included 2434 HT users and 12 793 nonusers. In primary GWAS meta-analyses, no variant reached genome-wide significance, but the smallest P value approached statistical significance: rs9386463 (P = 5.03 × 10-8). We tested associations for 138 candidate variants and identified 2 that exceeded statistical significance (0.05/138 = 3.62 × 10-4): F5 rs6025 (P = 1.87 × 10-5; SI, 1.29; previously observed) and F11 rs2036914 (P = 2.0 × 10-4; SI, 0.91; new observation). CONCLUSION: The candidate variant approach to identify submultiplictive/supramultiplicative associations between genetic variation and OC and HT use identified a new association with common genetic variation in F11, while the agnostic interrogations did not yield new discoveries.

Genome-wide interaction study of dietary intake of fibre, fruits, and vegetables with risk of colorectal cancer.

BACKGROUND: Consumption of fibre, fruits and vegetables have been linked with lower colorectal cancer (CRC) risk. A genome-wide gene-environment (G × E) analysis was performed to test whether genetic variants modify these associations. METHODS: A pooled sample of 45 studies including up to 69,734 participants (cases: 29,896; controls: 39,838) of European ancestry were included. To identify G × E interactions, we used the traditional 1--degree-of-freedom (DF) G × E test and to improve power a 2-step procedure and a 3DF joint test that investigates the association between a genetic variant and dietary exposure, CRC risk and G × E interaction simultaneously. FINDINGS: The 3-DF joint test revealed two significant loci with p-value <5 × 10-8. Rs4730274 close to the SLC26A3 gene showed an association with fibre (p-value: 2.4 × 10-3) and G × fibre interaction with CRC (OR per quartile of fibre increase = 0.87, 0.80, and 0.75 for CC, TC, and TT genotype, respectively; G × E p-value: 1.8 × 10-7). Rs1620977 in the NEGR1 gene showed an association with fruit intake (p-value: 1.0 × 10-8) and G × fruit interaction with CRC (OR per quartile of fruit increase = 0.75, 0.65, and 0.56 for AA, AG, and GG genotype, respectively; G × E -p-value: 0.029). INTERPRETATION: We identified 2 loci associated with fibre and fruit intake that also modify the association of these dietary factors with CRC risk. Potential mechanisms include chronic inflammatory intestinal disorders, and gut function. However, further studies are needed for mechanistic validation and replication of findings. FUNDING: National Institutes of Health, National Cancer Institute. Full funding details for the individual consortia are provided in acknowledgments.

Association of MTHFR (C677T, A1298C) and MTRR A66G polymorphisms with fatty acids profile and risk of neural tube defects.

OBJECTIVE: This study aims to determine if 5,10-methylenetetrahydrofolate reductase (MTHFR C677T and A1298C) and methionine synthase reductase (MTRR A66G) gene polymorphisms were associated with fatty acid (FA) levels in mothers of fetuses with neural tube defects (NTDs) and whether these associations were modified by environmental factors. METHODS: Plasma FA composition was assessed using capillary gas chromatography. Concentrations of studied FA were compared between 42 mothers of NTDs fetuses and 30 controls as a function of each polymorphism by the Kruskal-Wallis nonparametric test. RESULTS: In MTHFR gene C677T polymorphism, cases with (CT + TT) genotype had lower monounsaturated FAs (MUFA) and omega-3 polyunsaturated FA (n-3 PUFA) levels, but higher omega-6 polyunsaturated FAs (n-6 PUFA) and omega-6 polyunsaturated FAs: omega-3 polyunsaturated FAs (n-6:n-3) ratio levels. In MTRR gene A66G polymorphism, cases with (AG + GG) genotype had lower MUFA levels, but higher PUFA and n-6 PUFA levels. Controls with (AG + GG) genotype had lower n-6 PUFA levels. In MTHFR gene C677T polymorphism, cases with smoking spouses and (CT + TT) genotype had lower MUFA and n-3 PUFA levels, but higher PUFA, n-6 PUFA, and n-6:n-3 ratio levels. Cases with (CT + TT) genotype and who used sauna during pregnancy had lower n-3 PUFA levels. In MTRR gene A66G polymorphism, cases with (AG + GG) genotype and who used sauna during pregnancy had higher PUFA and n-6 PUFA levels. CONCLUSIONS: Further research is required to clarify the association of FA metabolism and (MTHFR, MTRR) polymorphisms with NTDs.

No gene by stressful life events interaction on individual differences in adults' self-control.

Background: Difficulty with self-control, or the ability to alter impulses and behavior in a goal-directed way, predicts interpersonal conflict, lower socioeconomic attainments, and more adverse health outcomes. Etiological understanding, and intervention for low self-control is, therefore, a public health goal. A prominent developmental theory proposes that individuals with high genetic propensity for low self-control that are also exposed to stressful environments may be most at-risk of low levels of self-control. Here we examine if polygenic measures associated with behaviors marked by low self-control interact with stressful life events in predicting self-control. Methods: Leveraging molecular data from a large population-based Dutch sample (N = 7,090, Mage = 41.2) to test for effects of genetics (i.e., polygenic scores for ADHD and aggression), stressful life events (e.g., traffic accident, violent assault, financial problems), and a gene-by-stress interaction on self-control (measured with the ASEBA Self-Control Scale). Results: Both genetics (ß =.03 -.04, p <.001) and stressful life events (ß = .11 -.14, p <.001) were associated with individual differences in self-control. We find no evidence of a gene-by-stressful life events interaction on individual differences in adults' self-control. Conclusion: Our findings are consistent with the notion that genetic influences and stressful life events exert largely independent effects on adult self-control. However, the small effect sizes of polygenic scores increases the likelihood of null results. Genetically-informed longitudinal research in large samples can further inform the etiology of individual differences in self-control from early childhood into later adulthood and its downstream implications for public health.