Epigenome-wide association study of lung cancer among never smokers in two prospective cohorts in Shanghai, China.

BACKGROUND: The aetiology of lung cancer among individuals who never smoked remains elusive, despite 15% of lung cancer cases in men and 53% in women worldwide being unrelated to smoking. Epigenetic alterations, particularly DNA methylation (DNAm) changes, have emerged as potential drivers. Yet, few prospective epigenome-wide association studies (EWAS), primarily focusing on peripheral blood DNAm with limited representation of never smokers, have been conducted. METHODS: We conducted a nested case-control study of 80 never-smoking incident lung cancer cases and 83 never-smoking controls within the Shanghai Women's Health Study and Shanghai Men's Health Study. DNAm was measured in prediagnostic oral rinse samples using Illumina MethylationEPIC array. Initially, we conducted an EWAS to identify differentially methylated positions (DMPs) associated with lung cancer in the discovery sample of 101 subjects. The top 50 DMPs were further evaluated in a replication sample of 62 subjects, and results were pooled using fixed-effect meta-analysis. RESULTS: Our study identified three DMPs significantly associated with lung cancer at the epigenome-wide significance level of p<8.22×10-8. These DMPs were identified as cg09198866 (MYH9; TXN2), cg01411366 (SLC9A10) and cg12787323. Furthermore, examination of the top 1000 DMPs indicated significant enrichment in epithelial regulatory regions and their involvement in small GTPase-mediated signal transduction pathways. Additionally, GrimAge acceleration was identified as a risk factor for lung cancer (OR=1.19 per year; 95% CI 1.06 to 1.34). CONCLUSIONS: While replication in a larger sample size is necessary, our findings suggest that DNAm patterns in prediagnostic oral rinse samples could provide novel insights into the underlying mechanisms of lung cancer in never smokers.

Prospective epigenome and transcriptome analyses of cord and peripheral blood from preterm infants at risk of bronchopulmonary dysplasia.

Bronchopulmonary dysplasia (BPD) is a prevalent chronic lung disease of prematurity with limited treatment options. To uncover biomarkers of BPD risk, this study investigated epigenetic and transcriptomic signatures of prematurity at birth and during the neonatal period at day 14 and 28. Peripheral blood DNAs from preterm infants were applied to methylation arrays and cell-type composition was estimated by deconvolution. Covariate-adjusted robust linear regression elucidated BPD- and prolonged oxygen (≥ 14 days) exposure-associated CpGs. RNAs from cord and peripheral blood were sequenced, and differentially expressed genes (DEGs) for BPD or oxygen exposure were determined. Estimated neutrophil-lymphocyte ratios in peripheral blood at day 14 in BPD infants were significantly higher than nonBPD infants, suggesting an heightened inflammatory response in developing BPD. BPD-DEGs in cord blood indicated lymphopoiesis inhibition, altered Th1/Th2 responses, DNA damage, and organ degeneration. On day 14, BPD-associated CpGs were highly enriched in neutrophil activation, infection, and CD4 + T cell quantity, and BPD-DEGs were involved in DNA damage, cellular senescence, T cell homeostasis, and hyper-cytokinesis. On day 28, BPD-associated CpGs along with BPD-DEGs were enriched for phagocytosis, neurological disorder, and nucleotide metabolism. Oxygen supplementation markedly downregulated mitochondrial biogenesis genes and altered CpGs annotated to developmental genes. Prematurity-altered DNA methylation could cause abnormal lymphopoiesis, cellular assembly and cell cycle progression to increase BPD risk. Similar pathways between epigenome and transcriptome networks suggest coordination of the two in dysregulating leukopoiesis, adaptive immunity, and innate immunity. The results provide molecular insights into biomarkers for early detection and prevention of BPD.

Whole-genome survey reveals extensive variation in genetic diversity and inbreeding levels among peregrine falcon subspecies.

In efforts to prevent extinction, resource managers are often tasked with increasing genetic diversity in a population of concern to prevent inbreeding depression or improve adaptive potential in a changing environment. The assumption that all small populations require measures to increase their genetic diversity may be unwarranted, and limited resources for conservation may be better utilized elsewhere. We test this assumption in a case study focused on the peregrine falcon (Falco peregrinus), a cosmopolitan circumpolar species with 19 named subspecies. We used whole-genome resequencing to generate over two million single nucleotide polymorphisms (SNPs) from multiple individuals of all peregrine falcon subspecies. Our analyses revealed extensive variation among subspecies, with many island-restricted and nonmigratory populations possessing lower overall genomic diversity, elevated inbreeding coefficients (F ROH)-among the highest reported, and extensive runs of homozygosity (ROH) compared to mainland and migratory populations. Similarly, the majority of subspecies that are either nonmigratory or restricted to islands show a much longer history of low effective population size (N e). While mutational load analyses indicated an increased proportion of homozygous-derived deleterious variants (i.e., drift load) among nonmigrant and island populations compared to those that are migrant or reside on the mainland, no significant differences in the proportion of heterozygous deleterious variants (i.e., inbreeding load) was observed. Our results provide evidence that high levels of inbreeding may not be an existential threat for some populations or taxa. Additional factors such as the timing and severity of population declines are important to consider in management decisions about extinction potential.

Epigenomic profiling of isolated blood cell types reveals highly specific B cell smoking signatures and links to disease risk.

BACKGROUND: Tobacco smoking alters the DNA methylation profiles of immune cells which may underpin some of the pathogenesis of smoking-associated diseases. To link smoking-driven epigenetic effects in specific immune cell types with disease risk, we isolated six leukocyte subtypes, CD14+ monocytes, CD15+ granulocytes, CD19+ B cells, CD4+ T cells, CD8+ T cells, and CD56+ natural killer cells, from whole blood of 67 healthy adult smokers and 74 nonsmokers for epigenome-wide association study (EWAS) using Illumina 450k and EPIC methylation arrays. RESULTS: Numbers of smoking-associated differentially methylated sites (smCpGs) at genome-wide significance (p < 1.2 × 10-7) varied widely across cell types, from 5 smCpGs in CD8+ T cells to 111 smCpGs in CD19+ B cells. We found unique smoking effects in each cell type, some of which were not apparent in whole blood. Methylation-based deconvolution to estimate B cell subtypes revealed that smokers had 7.2% (p = 0.033) less naïve B cells. Adjusting for naïve and memory B cell proportions in EWAS and RNA-seq allowed the identification of genes enriched for B cell activation-related cytokine signaling pathways, Th1/Th2 responses, and hematopoietic cancers. Integrating with large-scale public datasets, 62 smCpGs were among CpGs associated with health-relevant EWASs. Furthermore, 74 smCpGs had reproducible methylation quantitative trait loci single nucleotide polymorphisms (SNPs) that were in complete linkage disequilibrium with genome-wide association study SNPs, associating with lung function, disease risks, and other traits. CONCLUSIONS: We observed blood cell-type-specific smCpGs, a naïve-to-memory shift among B cells, and by integrating genome-wide datasets, we identified their potential links to disease risks and health traits.

Epigenome-wide association study of bronchopulmonary dysplasia in preterm infants: results from the discovery-BPD program.

BACKGROUND: Bronchopulmonary dysplasia (BPD) is a lung disease in premature infants caused by therapeutic oxygen supplemental and characterized by impaired pulmonary development which persists into later life. While advances in neonatal care have improved survival rates of premature infants, cases of BPD have been increasing with limited therapeutic options for prevention and treatment. This study was designed to explore the relationship between gestational age (GA), birth weight, and estimated blood cell-type composition in premature infants and to elucidate early epigenetic biomarkers associated with BPD. METHODS: Cord blood DNA from preterm neonates that went on to develop BPD (n = 14) or not (non-BPD, n = 93) was applied to Illumina 450 K methylation arrays. Blood cell-type compositions were estimated using DNA methylation profiles. Multivariable robust regression analysis elucidated CpGs associated with BPD risk. cDNA microarray analysis of cord blood RNA identified differentially expressed genes in neonates who later developed BPD. RESULTS: The development of BPD and the need for oxygen supplementation were strongly associated with GA (BPD, p < 1.0E-04; O2 supplementation, p < 1.0E-09) and birth weight (BPD, p < 1.0E-02; O2 supplementation, p < 1.0E-07). The estimated nucleated red blood cell (NRBC) percent was negatively associated with birth weight and GA, positively associated with hypomethylation of the tobacco smoke exposure biomarker cg05575921, and high-NRBC blood samples displayed a hypomethylation profile. Epigenome-wide association study (EWAS) identified 38 (Bonferroni) and 275 (false discovery rate 1%) differentially methylated CpGs associated with BPD. BPD-associated CpGs in cord blood were enriched for lung maturation and hematopoiesis pathways. Stochastic epigenetic mutation burden at birth was significantly elevated among those who developed BPD (adjusted p = 0.02). Transcriptome changes in cord blood cells reflected cell cycle, development, and pulmonary disorder events in BPD. CONCLUSIONS: While results must be interpreted with caution because of the small size of this study, NRBC content strongly impacted DNA methylation profiles in preterm cord blood and EWAS analysis revealed potential insights into biological pathways involved in BPD pathogenesis.

Circulating MicroRNAs, Polychlorinated Biphenyls, and Environmental Liver Disease in the Anniston Community Health Survey.

BACKGROUND: Polychlorinated biphenyl (PCB) exposures have been associated with liver injury in human cohorts, and steatohepatitis with liver necrosis in model systems. MicroRNAs (miRs) maintain cellular homeostasis and may regulate the response to environmental stress. OBJECTIVES: We tested the hypothesis that specific miRs are associated with liver disease and PCB exposures in a residential cohort. METHODS: Sixty-eight targeted hepatotoxicity miRs were measured in archived serum from 734 PCB-exposed participants in the cross-sectional Anniston Community Health Survey. Necrotic and other liver disease categories were defined by serum keratin 18 (K18) biomarkers. Associations were determined between exposure biomarkers (35 ortho-substituted PCB congeners) and disease biomarkers (highly expressed miRs or previously measured cytokines), and Ingenuity Pathway Analysis was performed. RESULTS: The necrotic liver disease category was associated with four up-regulated miRs (miR-99a-5p, miR-122-5p, miR-192-5p, and miR-320a) and five down-regulated miRs (let-7d-5p, miR-17-5p, miR-24-3p, miR-197-3p, and miR-221-3p). Twenty-two miRs were associated with the other liver disease category or with K18 measurements. Eleven miRs were associated with 24 PCBs, most commonly congeners with anti-estrogenic activities. Most of the exposure-associated miRs were associated with at least one serum hepatocyte death, pro-inflammatory cytokine or insulin resistance bioarker, or with both. Within each biomarker category, associations were strongest for the liver-specific miR-122-5p. Pathways of liver toxicity that were identified included inflammation/hepatitis, hyperplasia/hyperproliferation, cirrhosis, and hepatocellular carcinoma. Tumor protein p53 and tumor necrosis factor α were well integrated within the top identified networks. DISCUSSION: These results support the human hepatotoxicity of environmental PCB exposures while elucidating potential modes of PCB action. The MiR-derived liquid liver biopsy represents a promising new technique for environmental hepatology cohort studies. https://doi.org/10.1289/EHP9467.

Age-specific survival rates, causes of death, and allowable take of golden eagles in the western United States.

In the United States, the Bald and Golden Eagle Protection Act prohibits take of golden eagles (Aquila chrysaetos) unless authorized by permit, and stipulates that all permitted take must be sustainable. Golden eagles are unintentionally killed in conjunction with many lawful activities (e.g., electrocution on power poles, collision with wind turbines). Managers who issue permits for incidental take of golden eagles must determine allowable take levels and manage permitted take accordingly. To aid managers in making these decisions in the western United States, we used an integrated population model to obtain estimates of golden eagle vital rates and population size, and then used those estimates in a prescribed take level (PTL) model to estimate the allowable take level. Estimated mean annual survival rates for golden eagles ranged from 0.70 (95% credible interval = 0.66-0.74) for first-year birds to 0.90 (0.88-0.91) for adults. Models suggested a high proportion of adult female golden eagles attempted to breed and breeding pairs fledged a mean of 0.53 (0.39-0.72) young annually. Population size in the coterminous western United States has averaged ~31,800 individuals for several decades, with λ = 1.0 (0.96-1.05). The PTL model estimated a median allowable take limit of ~2227 (708-4182) individuals annually given a management objective of maintaining a stable population. We estimate that take averaged 2572 out of 4373 (59%) deaths annually, based on a representative sample of transmitter-tagged golden eagles. For the subset of golden eagles that were recovered and a cause of death determined, anthropogenic mortality accounted for an average of 74% of deaths after their first year; leading forms of take over all age classes were shooting (~670 per year), collisions (~611), electrocutions (~506), and poisoning (~427). Although observed take overlapped the credible interval of our allowable take estimate and the population overall has been stable, our findings indicate that additional take, unless mitigated for, may not be sustainable. Our analysis demonstrates the utility of the joint application of integrated population and prescribed take level models to management of incidental take of a protected species.

Does taxonomic and numerical resolution affect the assessment of invertebrate community structure in New World freshwater wetlands?

The efficiency of biodiversity assessments and biomonitoring studies is commonly challenged by limitations in taxonomic identification and quantification approaches. In this study, we assessed the effects of different taxonomic and numerical resolutions on a range of community structure metrics in invertebrate compositional data sets from six regions distributed across North and South America. We specifically assessed the degree of similarity in the metrics (richness, equitability, beta diversity, heterogeneity in community composition and congruence) for data sets identified to a coarse resolution (usually family level) and the finest taxonomic resolution practical (usually genus level, sometimes species or morphospecies) and by presence-absence and relative abundance numerical resolutions. Spearman correlations showed highly significant and positive associations between univariate metrics (richness and equitability) calculated for coarse- and finest-resolution datasets. Procrustes analysis detected significant congruence between composition datasets. Higher correlation coefficients were found for datasets with the same numerical resolutions regardless of the taxonomic level (about 90%), while the correlations for comparisons across numerical resolutions were consistently lower. Our findings indicate that family-level resolution can be used as a surrogate of finer taxonomic resolutions to calculate a range of biodiversity metrics commonly used to describe invertebrate community structure patterns in New World freshwater wetlands without significant loss of information. However, conclusions on biodiversity patterns derived from datasets with different numerical resolutions should be critically considered in studies on wetland invertebrates.

Germline and Somatic Genetic Variants in the p53 Pathway Interact to Affect Cancer Risk, Progression, and Drug Response.

Insights into oncogenesis derived from cancer susceptibility loci (SNP) hold the potential to facilitate better cancer management and treatment through precision oncology. However, therapeutic insights have thus far been limited by our current lack of understanding regarding both interactions of these loci with somatic cancer driver mutations and their influence on tumorigenesis. For example, although both germline and somatic genetic variation to the p53 tumor suppressor pathway are known to promote tumorigenesis, little is known about the extent to which such variants cooperate to alter pathway activity. Here we hypothesize that cancer risk-associated germline variants interact with somatic TP53 mutational status to modify cancer risk, progression, and response to therapy. Focusing on a cancer risk SNP (rs78378222) with a well-documented ability to directly influence p53 activity as well as integration of germline datasets relating to cancer susceptibility with tumor data capturing somatically-acquired genetic variation provided supportive evidence for this hypothesis. Integration of germline and somatic genetic data enabled identification of a novel entry point for therapeutic manipulation of p53 activities. A cluster of cancer risk SNPs resulted in increased expression of prosurvival p53 target gene KITLG and attenuation of p53-mediated responses to genotoxic therapies, which were reversed by pharmacologic inhibition of the prosurvival c-KIT signal. Together, our results offer evidence of how cancer susceptibility SNPs can interact with cancer driver genes to affect cancer progression and identify novel combinatorial therapies. SIGNIFICANCE: These results offer evidence of how cancer susceptibility SNPs can interact with cancer driver genes to affect cancer progression and present novel therapeutic targets.

Single-cell analyses identify dysfunctional CD16+ CD8 T cells in smokers.

Tobacco smoke exposure contributes to the global burden of communicable and chronic diseases. To identify immune cells affected by smoking, we use single-cell RNA sequencing on peripheral blood from smokers and nonsmokers. Transcriptomes reveal a subpopulation of FCGR3A (CD16)-expressing Natural Killer (NK)-like CD8 T lymphocytes that increase in smokers. Mass cytometry confirms elevated CD16+ CD8 T cells in smokers. Inferred as highly differentiated by pseudotime analysis, NK-like CD8 T cells express markers characteristic of effector memory re-expressing CD45RA T (TEMRA) cells. Indicative of immune aging, smokers' CD8 T cells are biased toward differentiated cells and smokers have fewer naïve cells than nonsmokers. DNA methylation-based models show that smoking dose is associated with accelerated aging and decreased telomere length, a biomarker of T cell senescence. Immune aging accompanies T cell senescence, which can ultimately lead to impaired immune function. This suggests a role for smoking-induced, senescence-associated immune dysregulation in smoking-mediated pathologies.

Mining a human transcriptome database for chemical modulators of NRF2.

Nuclear factor erythroid-2 related factor 2 (NRF2) encoded by the NFE2L2 gene is a transcription factor critical for protecting cells from chemically-induced oxidative stress. We developed computational procedures to identify chemical modulators of NRF2 in a large database of human microarray data. A gene expression biomarker was built from statistically-filtered gene lists derived from microarray experiments in primary human hepatocytes and cancer cell lines exposed to NRF2-activating chemicals (oltipraz, sulforaphane, CDDO-Im) or in which the NRF2 suppressor Keap1 was knocked down by siRNA. Directionally consistent biomarker genes were further filtered for those dependent on NRF2 using a microarray dataset from cells after NFE2L2 siRNA knockdown. The resulting 143-gene biomarker was evaluated as a predictive tool using the correlation-based Running Fisher algorithm. Using 59 gene expression comparisons from chemically-treated cells with known NRF2 activating potential, the biomarker gave a balanced accuracy of 93%. The biomarker was comprised of many well-known NRF2 target genes (AKR1B10, AKR1C1, NQO1, TXNRD1, SRXN1, GCLC, GCLM), 69% of which were found to be bound directly by NRF2 using ChIP-Seq. NRF2 activity was assessed across ~9840 microarray comparisons from ~1460 studies examining the effects of ~2260 chemicals in human cell lines. A total of 260 and 43 chemicals were found to activate or suppress NRF2, respectively, most of which have not been previously reported to modulate NRF2 activity. Using a NRF2-responsive reporter gene in HepG2 cells, we confirmed the activity of a set of chemicals predicted using the biomarker. The biomarker will be useful for future gene expression screening studies of environmentally-relevant chemicals.

Dioxin-like compound exposures and DNA methylation in the Anniston Community Health Survey Phase II.

The Anniston Community Health Survey (ACHS-I) was initially conducted from 2005 to 2007 to assess polychlorinated biphenyl (PCB) exposures in Anniston, Alabama residents. In 2014, a follow-up study (ACHS-II) was conducted to measure the same PCBs as in ACHS-I and additional compounds e.g., polychlorinated dibenzo-p-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs), and dioxin-like non-ortho (cPCBs) substituted PCBs. In this epigenome-wide association study (EWAS), we examined the associations between PCDD, PCDF, and PCB exposures and DNA methylation. Whole blood DNA methylation was measured using Illumina EPIC arrays (n=292). We modeled lipid-adjusted toxic equivalencies (TEQs) for: ΣDioxins (sum of 28 PCDDs, PCDFs, cPCBs, and mPCBs), PCDDs, PCDFs, cPCBs, and mPCBs using robust multivariable linear regression adjusting for age, race, sex, smoking, bisulfite conversion batch, and estimated percentages of six blood cell types. Among all exposures we identified 10 genome-wide (Bonferroni p≤6.74E-08) and 116 FDR (p≤5.00E-02) significant associations representing 10 and 113 unique CpGs, respectively. Of the 10 genome-wide associations, seven (70%) occurred in the PCDDs and four (40%) of these associations had an absolute differential methylation ≥1.00%, based on the methylation difference between the highest and lowest exposure quartiles. Most of the associations (six, 60%) represented hypomethylation changes. Of the 10 unique CpGs, eight (80%) were in genes shown to be associated with dioxins and/or PCBs based on data from the 2019 Comparative Toxicogenomics Database. In this study, we have identified a set of CpGs in blood DNA that may be particularly susceptible to dioxin, furan, and dioxin-like PCB exposures.

Polychlorinated biphenyl exposure and DNA methylation in the Anniston Community Health Survey.

Anniston, Alabama was home to a major polychlorinated biphenyl (PCB) production facility from 1929 until 1971. The Anniston Community Health Survey I and II (ACHS-I 2005-2007, ACHS-II 2013-2014) were conducted to explore the effects of PCB exposures. In this report we examined associations between PCB exposure and DNA methylation in whole blood using EPIC arrays (ACHS-I, n = 518; ACHS-II, n = 299). For both cohorts, 35 PCBs were measured in serum. We modelled methylation versus PCB wet-weight concentrations for: the sum of 35 PCBs, mono-ortho substituted PCBs, di-ortho substituted PCBs, tri/tetra-ortho substituted PCBs, oestrogenic PCBs, and antiestrogenic PCBs. Using robust multivariable linear regression, we adjusted for age, race, sex, smoking, total lipids, and six blood cell-type percentages. We carried out a two-stage analysis; discovery in ACHS-I followed by replication in ACHS-II. In ACHS-I, we identified 28 associations (17 unique CpGs) at p ≤ 6.70E-08 and 369 associations (286 unique CpGs) at FDR p ≤ 5.00E-02. A large proportion of the genes have been observed to interact with PCBs or dioxins in model studies. Among the 28 genome-wide significant CpG/PCB associations, 14 displayed replicated directional effects in ACHS-II; however, only one in ACHS-II was statistically significant at p ≤ 1.70E-04. While we identified many novel CpGs significantly associated with PCB exposures in ACHS-I, the differential methylation was modest and the effect was attenuated seven years later in ACHS-II, suggesting a lack of persistence of the associations between PCB exposures and altered DNA methylation in blood cells.

Elevated Heterozygosity in Adults Relative to Juveniles Provides Evidence of Viability Selection on Eagles and Falcons.

Viability selection yields adult populations that are more genetically variable than those of juveniles, producing a positive correlation between heterozygosity and survival. Viability selection could be the result of decreased heterozygosity across many loci in inbred individuals and a subsequent decrease in survivorship resulting from the expression of the deleterious alleles. Alternatively, locus-specific differences in genetic variability between adults and juveniles may be driven by forms of balancing selection, including heterozygote advantage, frequency-dependent selection, or selection across temporal and spatial scales. We use a pooled-sequencing approach to compare genome-wide and locus-specific genetic variability between 74 golden eagle (Aquila chrysaetos), 62 imperial eagle (Aquila heliaca), and 69 prairie falcon (Falco mexicanus) juveniles and adults. Although genome-wide genetic variability is comparable between juvenile and adult golden eagles and prairie falcons, imperial eagle adults are significantly more heterozygous than juveniles. This evidence of viability selection may stem from a relatively smaller imperial eagle effective population size and potentially greater genetic load. We additionally identify ~2000 single-nucleotide polymorphisms across the 3 species with extreme differences in heterozygosity between juveniles and adults. Many of these markers are associated with genes implicated in immune function or olfaction. These loci represent potential targets for studies of how heterozygote advantage, frequency-dependent selection, and selection over spatial and temporal scales influence survivorship in avian species. Overall, our genome-wide data extend previous studies that used allozyme or microsatellite markers and indicate that viability selection may be a more common evolutionary phenomenon than often appreciated.

Smoking-associated AHRR demethylation in cord blood DNA: impact of CD235a+ nucleated red blood cells.

BACKGROUND: Numerous studies have demonstrated that DNA methylation levels in the aryl hydrocarbon receptor repressor (AHRR) gene measured in cord blood are significantly associated with prenatal tobacco smoke exposure and can be used as a fetal exposure biomarker. The mechanism driving this demethylation has not been determined and it is unclear if all cord blood cell types are impacted. Nucleated red blood cells (nRBCs/CD235a+ cells) are developmentally immature RBCs that display genome-wide hypomethylation and are observed at increased frequency in the cord blood of smoking mothers. We tested if AHRR methylation levels in CD235a+ nRBCs or nRBC counts influenced AHRR methylation in whole cord blood. METHODS: Cord blood was collected from smoking (n = 34) and nonsmoking (n = 19) mothers and DNA was prepared from whole cord blood, isolated CD235a+ nRBCs, and CD14+ monocytes. AHRR methylation in cord blood DNA was measured using Illumina 850K arrays (cg05575921, chr5:373378). Pyrosequencing was used to compare methylation levels among cord blood, CD235a+, and CD14+ cells. We measured nRBC percentages using conventional complete blood counts and estimated percent nRBCs by a deconvolution model. RESULTS: Methylation levels in AHRR were significantly lower in nRBCs relative to whole cord blood and CD14+ monocytes. While AHRR methylation levels in the cell types were significantly correlated across all subjects, methylation values at the chr5:373378 CpG averaged 14.6% lower in nRBCs (range 0.4 to 24.8%; p = 3.8E-13) relative to CD14+, with nonsmokers showing a significantly greater hypomethylation (- 4.1%, p = 1.8E-02). Methylation level at the AHRR chr5:373378 CpG was strongly associated with self-reported smoking in both CD14+ monocytes (t test p = 5.7E-09) and nRBCs (p = 4.8E-08), as well as cotinine levels (regression p = 1.1E-07 and p = 3.6E-04, respectively). For subjects with whole blood 850K data, robust linear regression models adjusting for estimated cell type composition, either including nRBCs counts or estimates, modestly increased the association between smoking and cg05575921 methylation. CONCLUSIONS: Prenatal smoke exposure was highly significantly associated with AHRR methylation in cord blood, CD14+ monocytes, and CD235a+ nRBCs. AHRR methylation levels in nRBCs and nRBC counts had minimal effect on cord blood methylation measurements. However, regression models using estimated nRBCs or actual nRBC counts outperformed those lacking these covariates.

Associations between Maternal Tobacco Smoke Exposure and the Cord Blood [Formula: see text] DNA Methylome.

BACKGROUND: Maternal tobacco smoke exposure has been associated with altered DNA methylation. However, previous studies largely used methylation arrays, which cover a small fraction of CpGs, and focused on whole cord blood. OBJECTIVES: The current study examined the impact of in utero exposure to maternal tobacco smoke on the cord blood [Formula: see text] DNA methylome. METHODS: The methylomes of 20 Hispanic white newborns ([Formula: see text] exposed to any maternal tobacco smoke in pregnancy; [Formula: see text] unexposed) from the Maternal and Child Health Study (MACHS) were profiled by whole-genome bisulfite sequencing (median coverage: [Formula: see text]). Statistical analyses were conducted using the Regression Analysis of Differential Methylation (RADMeth) program because it performs well on low-coverage data (minimizes false positives and negatives). RESULTS: We found that 10,381 CpGs were differentially methylated by tobacco smoke exposure [neighbor-adjusted p-values that are additionally corrected for multiple testing based on the Benjamini-Hochberg method for controlling the false discovery rate (FDR) [Formula: see text]]. From these CpGs, RADMeth identified 557 differentially methylated regions (DMRs) that were overrepresented ([Formula: see text]) in important regulatory regions, including enhancers. Of nine DMRs that could be queried in a reduced representation bisulfite sequencing (RRBS) study of adult [Formula: see text] cells ([Formula: see text] smokers; [Formula: see text] nonsmokers), four replicated ([Formula: see text]). Additionally, a CpG in the promoter of SLC7A8 (percent methylation difference: [Formula: see text] comparing exposed to unexposed) replicated ([Formula: see text]) in an EPIC (Illumina) array study of cord blood [Formula: see text] cells ([Formula: see text] exposed to sustained maternal tobacco smoke; [Formula: see text] unexposed) and in a study of adult [Formula: see text] cells across two platforms (EPIC: [Formula: see text] smokers; [Formula: see text] nonsmokers; 450K: [Formula: see text] smokers; [Formula: see text] nonsmokers). CONCLUSIONS: Maternal tobacco smoke exposure in pregnancy is associated with cord blood [Formula: see text] DNA methylation in key regulatory regions, including enhancers. While we used a method that performs well on low-coverage data, we cannot exclude the possibility that some results may be false positives. However, we identified a differentially methylated CpG in amino acid transporter SLC7A8 that is highly reproducible, which may be sensitive to cigarette smoke in both cord blood and adult [Formula: see text] cells. https://doi.org/10.1289/EHP3398.

Crohn's disease IRGM risk alleles are associated with altered gene expression in human tissues.

Crohn's disease (CD) is a chronic inflammatory gastrointestinal disorder. Genetic association studies have implicated dysregulated autophagy in CD. Among risk loci identified are a promoter single nucleotide polymorphism (SNP)( rs13361189 ) and two intragenic SNPs ( rs9637876 , rs10065172 ) in immunity-related GTPase family M ( IRGM) a gene that encodes a protein of the autophagy initiation complex. All three SNPs have been proposed to modify IRGM expression, but reports have been divergent and largely derived from cell lines. Here, analyzing RNA-Sequencing data of human tissues from the Genotype-Tissue Expression Project, we found that rs13361189 minor allele carriers had reduced IRGM expression in whole blood and terminal ileum, and upregulation in ileum of ZNF300P1, a locus adjacent to IRGM on chromosome 5q33.1 that encodes a long noncoding RNA. Whole blood and ileum from minor allele carriers had altered expression of multiple additional genes that have previously been linked to colitis and/or autophagy. Notable among these was an increase in ileum of LTF (lactoferrin), an established fecal inflammatory biomarker of CD, and in whole blood of TNF, a key cytokine in CD pathogenesis. Last, we confirmed that risk alleles at all three loci associated with increased risk for CD but not ulcerative colitis in a case-control study. Taken together, our findings suggest that genetically encoded IRGM deficiency may predispose to CD through dysregulation of inflammatory gene networks. Gene expression profiling of disease target tissues in genetically susceptible populations is a promising strategy for revealing new leads for the study of molecular pathogenesis and, potentially, for precision medicine. NEW & NOTEWORTHY Single nucleotide polymorphisms in immunity-related GTPase family M ( IRGM), a gene that encodes an autophagy initiation protein, have been linked epidemiologically to increased risk for Crohn's disease (CD). Here, we show for the first time that subjects with risk alleles at two such loci, rs13361189 and rs10065172 , have reduced IRGM expression in whole blood and terminal ileum, as well as dysregulated expression of a wide array of additional genes that regulate inflammation and autophagy.

Sulforaphane enriched transcriptome of lung mitochondrial energy metabolism and provided pulmonary injury protection via Nrf2 in mice.

Nrf2 is essential to antioxidant response element (ARE)-mediated host defense. Sulforaphane (SFN) is a phytochemical antioxidant known to affect multiple cellular targets including Nrf2-ARE pathway in chemoprevention. However, the role of SFN in non-malignant airway disorders remain unclear. To test if pre-activation of Nrf2-ARE signaling protects lungs from oxidant-induced acute injury, wild-type (Nrf2+/+) and Nrf2-deficient (Nrf2-/-) mice were given SFN orally or as standardized broccoli sprout extract diet (SBE) before hyperoxia or air exposure. Hyperoxia-induced pulmonary injury and oxidation indices were significantly reduced by SFN or SBE in Nrf2+/+ mice but not in Nrf2-/- mice. SFN upregulated a large cluster of basal lung genes that are involved in mitochondrial oxidative phosphorylation, energy metabolism, and cardiovascular protection only in Nrf2+/+ mice. Bioinformatic analysis elucidated ARE-like motifs on these genes. Transcript abundance of the mitochondrial machinery genes remained significantly higher after hyperoxia exposure in SFN-treated Nrf2+/+ mice than in SFN-treated Nrf2-/- mice. Nuclear factor-κB was suggested to be a central molecule in transcriptome networks affected by SFN. Minor improvement of hyperoxia-caused lung histopathology and neutrophilia by SFN in Nrf2-/- mice implies Nrf2-independent or alternate effector mechanisms. In conclusion, SFN is suggested to be as a preventive intervention in a preclinical model of acute lung injury by linking mitochondria and Nrf2. Administration of SFN alleviated acute lung injury-like pathogenesis in a Nrf2-dependent manner. Potential AREs in the SFN-inducible transcriptome for mitochondria bioenergetics provided a new insight into the downstream mechanisms of Nrf2-mediated pulmonary protection.

Epigenome-wide meta-analysis of DNA methylation and childhood asthma.

BACKGROUND: Epigenetic mechanisms, including methylation, can contribute to childhood asthma. Identifying DNA methylation profiles in asthmatic patients can inform disease pathogenesis. OBJECTIVE: We sought to identify differential DNA methylation in newborns and children related to childhood asthma. METHODS: Within the Pregnancy And Childhood Epigenetics consortium, we performed epigenome-wide meta-analyses of school-age asthma in relation to CpG methylation (Illumina450K) in blood measured either in newborns, in prospective analyses, or cross-sectionally in school-aged children. We also identified differentially methylated regions. RESULTS: In newborns (8 cohorts, 668 cases), 9 CpGs (and 35 regions) were differentially methylated (epigenome-wide significance, false discovery rate < 0.05) in relation to asthma development. In a cross-sectional meta-analysis of asthma and methylation in children (9 cohorts, 631 cases), we identified 179 CpGs (false discovery rate < 0.05) and 36 differentially methylated regions. In replication studies of methylation in other tissues, most of the 179 CpGs discovered in blood replicated, despite smaller sample sizes, in studies of nasal respiratory epithelium or eosinophils. Pathway analyses highlighted enrichment for asthma-relevant immune processes and overlap in pathways enriched both in newborns and children. Gene expression correlated with methylation at most loci. Functional annotation supports a regulatory effect on gene expression at many asthma-associated CpGs. Several implicated genes are targets for approved or experimental drugs, including IL5RA and KCNH2. CONCLUSION: Novel loci differentially methylated in newborns represent potential biomarkers of risk of asthma by school age. Cross-sectional associations in children can reflect both risk for and effects of disease. Asthma-related differential methylation in blood in children was substantially replicated in eosinophils and respiratory epithelium.

A distinct class of antioxidant response elements is consistently activated in tumors with NRF2 mutations.

NRF2 is a redox-responsive transcription factor that regulates expression of cytoprotective genes via its interaction with DNA sequences known as antioxidant response elements (AREs). NRF2 activity is induced by oxidative stress, but oxidative stress is not the only context in which NRF2 can be activated. Mutations that disrupt the interaction between NRF2 and KEAP1, an inhibitor of NRF2, lead to NRF2 hyperactivation and promote oncogenesis. The mechanisms underlying NRF2's oncogenic properties remain unclear, but likely involve aberrant expression of select NRF2 target genes. We tested this possibility using an integrative genomics approach to get a precise view of the direct NRF2 target genes dysregulated in tumors with NRF2 hyperactivating mutations. This approach revealed a core set of 32 direct NRF2 targets that are consistently upregulated in NRF2 hyperactivated tumors. This set of NRF2 "cancer target genes" includes canonical redox-related NRF2 targets, as well as target genes that have not been previously linked to NRF2 activation. Importantly, NRF2-driven upregulation of this gene set is largely independent of the organ system where the tumor developed. One key distinguishing feature of these NRF2 cancer target genes is that they are regulated by high affinity AREs that fall within genomic regions possessing a ubiquitously permissive chromatin signature. This implies that these NRF2 cancer target genes are responsive to oncogenic NRF2 in most tissues because they lack the regulatory constraints that restrict expression of most other NRF2 target genes. This NRF2 cancer target gene set also serves as a reliable proxy for NRF2 activity, and high NRF2 activity is associated with significant decreases in survival in multiple cancer types. Overall, the pervasive upregulation of these NRF2 cancer targets across multiple cancers, and their association with negative outcomes, suggests that these will be central to dissecting the functional implications of NRF2 hyperactivation in several cancer contexts.