Exploiting meta-analysis of genome-wide interaction with serum 25-hydroxyvitamin D to identify novel genetic loci associated with pulmonary function.

BACKGROUND: Higher 25-hydroxyvitamin D (25(OH)D) concentrations in serum has a positive association with pulmonary function. Investigating genome-wide interactions with 25(OH)D may reveal new biological insights into pulmonary function. OBJECTIVES: We aimed to identify novel genetic variants associated with pulmonary function by accounting for 25(OH)D interactions. METHODS: We included 211,264 participants from the observational United Kingdom Biobank study with pulmonary function tests (PFTs), genome-wide genotypes, and 25(OH)D concentrations from 4 ancestral backgrounds-European, African, East Asian, and South Asian. Among PFTs, we focused on forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) because both were previously associated with 25(OH)D. We performed genome-wide association study (GWAS) analyses that accounted for variant×25(OH)D interaction using the joint 2 degree-of-freedom (2df) method, stratified by participants' smoking history and ancestry, and meta-analyzed results. We evaluated interaction effects to determine how variant-PFT associations were modified by 25(OH)D concentrations and conducted pathway enrichment analysis to examine the biological relevance of our findings. RESULTS: Our GWAS meta-analyses, accounting for interaction with 25(OH)D, revealed 30 genetic variants significantly associated with FEV1 or FVC (P2df <5.00×10-8) that were not previously reported for PFT-related traits. These novel variant signals were enriched in lung function-relevant pathways, including the p38 MAPK pathway. Among variants with genome-wide-significant 2df results, smoking-stratified meta-analyses identified 5 variants with 25(OH)D interactions that influenced FEV1 in both smoking groups (never smokers P1df interaction<2.65×10-4; ever smokers P1df interaction<1.71×10-5); rs3130553, rs2894186, rs79277477, and rs3130929 associations were only evident in never smokers, and the rs4678408 association was only found in ever smokers. CONCLUSION: Genetic variant associations with lung function can be modified by 25(OH)D, and smoking history can further modify variant×25(OH)D interactions. These results expand the known genetic architecture of pulmonary function and add evidence that gene-environment interactions, including with 25(OH)D and smoking, influence lung function.

Genetic and transcriptomic landscape of colonic diverticulosis.

OBJECTIVE: Colonic diverticulosis is a prevalent condition among older adults, marked by the presence of thin-walled pockets in the colon wall that can become inflamed, infected, haemorrhage or rupture. We present a case-control genetic and transcriptomic study aimed at identifying the genetic and cellular determinants underlying this condition and the relationship with other gastrointestinal disorders. DESIGN: We conducted DNA and RNA sequencing on colonic tissue from 404 patients with (N=172) and without (N=232) diverticulosis. We investigated variation in the transcriptome associated with diverticulosis and further integrated this variation with single-cell RNA-seq data from the human intestine. We also integrated our expression quantitative trait loci with genome-wide association study using Mendelian randomisation (MR). Furthermore, a Polygenic Risk Score analysis gauged associations between diverticulosis severity and other gastrointestinal disorders. RESULTS: We discerned 38 genes with differential expression and 17 with varied transcript usage linked to diverticulosis, indicating tissue remodelling as a primary diverticula formation mechanism. Diverticula formation was primarily linked to stromal and epithelial cells in the colon including endothelial cells, myofibroblasts, fibroblasts, goblet, tuft, enterocytes, neurons and glia. MR highlighted five genes including CCN3, CRISPLD2, ENTPD7, PHGR1 and TNFSF13, with potential causal effects on diverticulosis. Notably, ENTPD7 upregulation was confirmed in diverticulosis cases. Additionally, diverticulosis severity was positively correlated with genetic predisposition to diverticulitis. CONCLUSION: Our results suggest that tissue remodelling is a primary mechanism for diverticula formation. Individuals with an increased genetic proclivity to diverticulitis exhibit a larger numbers of diverticula on colonoscopy.

Investigating Associations of Omega-3 Fatty Acids, Lung Function Decline, and Airway Obstruction.

Rationale: Inflammation contributes to lung function decline and the development of chronic obstructive pulmonary disease. Omega-3 fatty acids have antiinflammatory properties and may benefit lung health. Objectives: To investigate associations of omega-3 fatty acids with lung function decline and incident airway obstruction in a diverse sample of adults from general-population cohorts. Methods: Complementary study designs: 1) longitudinal study of plasma phospholipid omega-3 fatty acids and repeated FEV1 and FVC measures in the NHLBI Pooled Cohorts Study and 2) two-sample Mendelian randomization (MR) study of genetically predicted omega-3 fatty acids and lung function parameters. Measurements and Main Results: The longitudinal study found that higher omega-3 fatty acid levels were associated with attenuated lung function decline in 15,063 participants, with the largest effect sizes for the most metabolically downstream omega-3 fatty acid, docosahexaenoic acid (DHA). An increase in DHA of 1% of total fatty acids was associated with attenuations of 1.4 ml/yr for FEV1 (95% confidence interval [CI], 1.1-1.8) and 2.0 ml/yr for FVC (95% CI, 1.6-2.4) and a 7% lower incidence of spirometry-defined airway obstruction (95% CI, 0.89-0.97). DHA associations persisted across sexes and smoking histories and in Black, White, and Hispanic participants, with associations of the largest magnitude in former smokers and Hispanic participants. The MR study showed similar trends toward positive associations of genetically predicted downstream omega-3 fatty acids with FEV1 and FVC. Conclusions: The longitudinal and MR studies provide evidence supporting beneficial effects of higher levels of downstream omega-3 fatty acids, especially DHA, on lung health.

Investigating associations of omega-3 fatty acids, lung function decline, and airway obstruction.

Rationale: Inflammation contributes to lung function decline and the development of chronic obstructive pulmonary disease. Omega-3 fatty acids have anti-inflammatory properties and may benefit lung health. Objectives: Investigate associations of omega-3 fatty acids with lung function decline and incident airway obstruction in adults of diverse races/ethnicities from general population cohorts. Methods: Complementary study designs: (1) longitudinal study of plasma phospholipid omega-3 fatty acids and repeated FEV 1 and FVC measures in the National Heart, Lung, and Blood Institute Pooled Cohorts Study, and (2) two-sample Mendelian Randomization (MR) study of genetically predicted omega-3 fatty acids and lung function parameters. Measurements and Main Results: The longitudinal study found that higher omega-3 fatty acid concentrations were associated with attenuated lung function decline in 15,063 participants, with the largest effect sizes for docosahexaenoic acid (DHA). One standard deviation higher DHA was associated with an attenuation of 1.8 mL/year for FEV 1 (95% confidence interval [CI] 1.3-2.2) and 2.4 mL/year for FVC (95% CI 1.9-3.0). One standard deviation higher DHA was also associated with a 9% lower incidence of spirometry-defined airway obstruction (95% CI 0.86-0.97). DHA associations persisted across sexes, smoking histories, and Black, white and Hispanic participants, with the largest magnitude associations in former smokers and Hispanics. The MR study showed positive associations of genetically predicted omega-3 fatty acids with FEV 1 and FVC, with statistically significant findings across multiple MR methods. Conclusions: The longitudinal and MR studies provide evidence supporting beneficial effects of higher circulating omega-3 fatty acids, especially DHA, on lung health.

Correcting signal biases and detecting regulatory elements in STARR-seq data.

High-throughput reporter assays such as self-transcribing active regulatory region sequencing (STARR-seq) have made it possible to measure regulatory element activity across the entire human genome at once. The resulting data, however, present substantial analytical challenges. Here, we identify technical biases that explain most of the variance in STARR-seq data. We then develop a statistical model to correct those biases and to improve detection of regulatory elements. This approach substantially improves precision and recall over current methods, improves detection of both activating and repressive regulatory elements, and controls for false discoveries despite strong local correlations in signal.

AP-1 subunits converge promiscuously at enhancers to potentiate transcription.

The AP-1 transcription factor (TF) dimer contributes to many biological processes and environmental responses. AP-1 can be composed of many interchangeable subunits. Unambiguously determining the binding locations of these subunits in the human genome is challenging because of variable antibody specificity and affinity. Here, we definitively establish the genome-wide binding patterns of five AP-1 subunits by using CRISPR to introduce a common antibody tag on each subunit. We find limited evidence for strong dimerization preferences between subunits at steady state and find that, under a stimulus, dimerization patterns reflect changes in the transcriptome. Further, our analysis suggests that canonical AP-1 motifs indiscriminately recruit all AP-1 subunits to genomic sites, which we term AP-1 hotspots. We find that AP-1 hotspots are predictive of cell type-specific gene expression and of genomic responses to glucocorticoid signaling (more so than super-enhancers) and are significantly enriched in disease-associated genetic variants. Together, these results support a model where promiscuous binding of many AP-1 subunits to the same genomic location play a key role in regulating cell type-specific gene expression and environmental responses.

Genome wide association analysis in a mouse advanced intercross line.

The LG/J x SM/J advanced intercross line of mice (LG x SM AIL) is a multigenerational outbred population. High minor allele frequencies, a simple genetic background, and the fully sequenced LG and SM genomes make it a powerful population for genome-wide association studies. Here we use 1,063 AIL mice to identify 126 significant associations for 50 traits relevant to human health and disease. We also identify thousands of cis- and trans-eQTLs in the hippocampus, striatum, and prefrontal cortex of ~200 mice. We replicate an association between locomotor activity and Csmd1, which we identified in an earlier generation of this AIL, and show that Csmd1 mutant mice recapitulate the locomotor phenotype. Our results demonstrate the utility of the LG x SM AIL as a mapping population, identify numerous novel associations, and shed light on the genetic architecture of mammalian behavior.

Quantifying the Impact of Non-coding Variants on Transcription Factor-DNA Binding.

Many recent studies have emphasized the importance of genetic variants and mutations in cancer and other complex human diseases. The overwhelming majority of these variants occur in non-coding portions of the genome, where they can have a functional impact by disrupting regulatory interactions between transcription factors (TFs) and DNA. Here, we present a method for assessing the impact of non-coding mutations on TF-DNA interactions, based on regression models of DNA-binding specificity trained on high-throughput in vitro data. We use ordinary least squares (OLS) to estimate the parameters of the binding model for each TF, and we show that our predictions of TF-binding changes due to DNA mutations correlate well with measured changes in gene expression. In addition, by leveraging distributional results associated with OLS estimation, for each predicted change in TF binding we also compute a normalized score (z-score) and a significance value (p-value) reflecting our confidence that the mutation affects TF binding. We use this approach to analyze a large set of pathogenic non-coding variants, and we show that these variants lead to significant differences in TF binding between alleles, compared to a control set of common variants. Thus, our results indicate that there is a strong regulatory component to the pathogenic non-coding variants identified thus far.