DNA methylation and allergic sensitizations: A genome-scale longitudinal study during adolescence.

BACKGROUND: The presence of allergic sensitization has a major influence on the development and course of common childhood conditions such as asthma and rhinitis. The etiology of allergic sensitization is poorly understood, and its underlying biological mechanisms are not well established. Several studies showed that DNA methylation (DNAm) at some CpGs is associated with allergic sensitization. However, no studies have focused on the critical adolescence period. METHODS: We assessed the association of pre- and postadolescence genome-wide DNAm with allergic sensitization against indoor, outdoor and food allergens, using linear mixed models. We hypothesized that DNAm is associated with sensitization in general, and with poly-sensitization status, and these associations are age- and gender-specific. We tested these hypotheses in the IoW cohort (n = 376) and examined the findings in the BAMSE cohort (n = 267). RESULTS: Via linear mixed models, we identified 35 CpGs in IoW associated with allergic sensitization (at false discovery rate of 0.05), of which 33 were available in BAMSE and replicated with respect to the direction of associations with allergic sensitization. At the 35 CpGs except for cg19210306 on C13orf27, a reduction in methylation among atopic subjects was observed, most notably for cg21220721 and cg11699125 (ACOT7). DNAm at cg10159529 was strongly correlated with expression of IL5RA in peripheral blood (P-value = 6.76 × 10-20 ). Three CpGs (cg14121142, cg23842695, and cg26496795) were identified in IoW with age-specific association between DNAm and allergic sensitization. CONCLUSION: In adolescence, the status of allergic sensitization was associated with DNAm differentiation and at some CpGs the association is likely to be age-specific.

Changes in DNA Methylation from Age 18 to Pregnancy in Type 1, 2, and 17 T Helper and Regulatory T-Cells Pathway Genes.

To succeed, pregnancies need to initiate immune biases towards T helper 2 (Th2) responses, yet little is known about what establishes this bias. Using the Illumina 450 K platform, we explored changes in DNA methylation (DNAm) of Th1, Th2, Th17, and regulatory T cell pathway genes before and during pregnancy. Female participants were recruited at birth (1989), and followed through age 18 years and their pregnancy (2011-2015). Peripheral blood DNAm was measured in 245 girls at 18 years; from among these girls, the DNAm of 54 women was repeatedly measured in the first (weeks 8-21, n = 39) and second (weeks 22-38, n = 35) halves of pregnancy, respectively. M-values (logit-transformed ß-values of DNAm) were analyzed: First, with repeated measurement models, cytosine-phosphate-guanine sites (CpGs) of pathway genes in pregnancy and at age 18 (nonpregnant) were compared for changes (p ≤ 0.05). Second, we tested how many of the 348 pathway-related CpGs changed compared to 10 randomly selected subsets of all other CpGs and compared to 10 randomly selected subsets of other CD4+-related CpGs (348 in each subset). Contrasted to the nonpregnant state, 27.7% of Th1-related CpGs changed in the first and 36.1% in the second half of pregnancy. Among the Th2 pathway CpGs, proportions of changes were 35.1% (first) and 33.8% (second half). The methylation changes suggest involvement of both Th1 and Th2 pathway CpGs in the immune bias during pregnancy. Changes in regulatory T cell and Th17 pathways need further exploration.

Comparison of different cell type correction methods for genome-scale epigenetics studies.

BACKGROUND: Whole blood is frequently utilized in genome-wide association studies of DNA methylation patterns in relation to environmental exposures or clinical outcomes. These associations can be confounded by cellular heterogeneity. Algorithms have been developed to measure or adjust for this heterogeneity, and some have been compared in the literature. However, with new methods available, it is unknown whether the findings will be consistent, if not which method(s) perform better. RESULTS: Methods: We compared eight cell-type correction methods including the method in the minfi R package, the method by Houseman et al., the Removing unwanted variation (RUV) approach, the methods in FaST-LMM-EWASher, ReFACTor, RefFreeEWAS, and RefFreeCellMix R programs, along with one approach utilizing surrogate variables (SVAs). We first evaluated the association of DNA methylation at each CpG across the whole genome with prenatal arsenic exposure levels and with cancer status, adjusted for estimated cell-type information obtained from different methods. We then compared CpGs showing statistical significance from different approaches. For the methods implemented in minfi and proposed by Houseman et al., we utilized homogeneous data with composition of some blood cells available and compared them with the estimated cell compositions. Finally, for methods not explicitly estimating cell compositions, we evaluated their performance using simulated DNA methylation data with a set of latent variables representing "cell types". RESULTS: Results from the SVA-based method overall showed the highest agreement with all other methods except for FaST-LMM-EWASher. Using homogeneous data, minfi provided better estimations on cell types compared to the originally proposed method by Houseman et al. Further simulation studies on methods free of reference data revealed that SVA provided good sensitivities and specificities, RefFreeCellMix in general produced high sensitivities but specificities tended to be low when confounding is present, and FaST-LMM-EWASher gave the lowest sensitivity but highest specificity. CONCLUSIONS: Results from real data and simulations indicated that SVA is recommended when the focus is on the identification of informative CpGs. When appropriate reference data are available, the method implemented in the minfi package is recommended. However, if no such reference data are available or if the focus is not on estimating cell proportions, the SVA method is suggested.

Genome-wide DNA methylation at birth in relation to in utero arsenic exposure and the associated health in later life.

BACKGROUND: In utero arsenic exposure may alter fetal developmental programming by altering DNA methylation, which may result in a higher risk of disease in later life. We evaluated the association between in utero arsenic exposure and DNA methylation (DNAm) in cord blood and its influence in later life. METHODS: Genome-wide DNA methylation in cord blood from 64 subjects in the Taiwanese maternal infant and birth cohort was analyzed. Robust regressions were applied to assess the association of DNA methylation with in utero arsenic exposure. Multiple testing was adjusted by controlling false discovery rate (FDR) of 0.05. The DAVID bioinformatics tool was implemented for functional annotation analyses on the detected CpGs. The identified CpGs were further tested in an independent cohort. For the CpGs replicated in the independent cohort, linear mixed models were applied to assess the association of DNA methylation with low-density lipoprotein (LDL) at different ages (2, 5, 8, 11 and 14 years). RESULTS: In total, 579 out of 385,183 CpGs were identified after adjusting for multiple testing (FDR = 0.05), of which ~60% were positively associated with arsenic exposure. Functional annotation analysis on these CpGs detected 17 KEGG pathways (FDR = 0.05) including pathways for cardiovascular diseases (CVD) and diabetes mellitus. In the independent cohort, about 46% (252 out of 553 CpGs) of the identified CpGs showed associations consistent with those in the study cohort. In total, 11 CpGs replicated in the independent cohort were in the pathways related to CVD and diabetes mellitus. Via longitudinal analyses, we found at 5 out of the 11 CpGs methylation was associated with LDL over time and interactions between DNA methylation and time were observed at 4 of the 5 CpGs, cg25189764 (coeff = 0.157, p-value = 0.047), cg04986899 (coeff. For interaction [coeff.int] = 0.030, p-value = 0.024), cg04903360 (coeff.int = 0.026, p-value = 0.032), cg08198265 (coeff.int = -0.063, p-value = 0.0021), cg10473311 (coeff.int = -0.021, p-value = 0.027). CONCLUSION: In utero arsenic exposure was associated with cord blood DNA methylation at various CpGs. The identified CpGs may help determine pathological epigenetic mechanisms linked to in utero arsenic exposure. Five CpGs (cg25189764, cg04986899, cg04903360, cg08198265 and cg10473311) may serve as epigenetic markers for changes in LDL later in life.

Changes in Bone Marrow Tumor and Immune Cells Correlate with Durability of Remissions Following BCMA CAR T Therapy in Myeloma.

Chimeric antigen-receptor (CAR) T cells lead to high response rates in myeloma, but most patients experience recurrent disease. We combined several high-dimensional approaches to study tumor/immune cells in the tumor microenvironment (TME) of myeloma patients pre- and post-B-cell maturation antigen (BCMA)-specific CAR T therapy. Lower diversity of pretherapy T-cell receptor (TCR) repertoire, presence of hyperexpanded clones with exhaustion phenotype, and BAFF+PD-L1+ myeloid cells in the marrow correlated with shorter progression-free survival (PFS) following CAR T therapy. In contrast, longer PFS was associated with an increased proportion of CLEC9A+ dendritic cells (DC), CD27+TCF1+ T cells with diverse T-cell receptors, and emergence of T cells expressing marrow-residence genes. Residual tumor cells at initial response express stemlike genes, and tumor recurrence was associated with the emergence of new dominant clones. These data illustrate a dynamic interplay between endogenous T, CAR T, myeloid/DC, and tumor compartments that affects the durability of response following CAR T therapy in myeloma. SIGNIFICANCE: There is an unmet need to identify determinants of durable responses following BCMA CAR T therapy of myeloma. High-dimensional analysis of the TME was performed to identify features of immune and tumor cells that correlate with survival and suggest several strategies to improve outcomes following CAR T therapy. See related commentary by Graham and Maus, p. 478. This article is highlighted in the In This Issue feature, p. 476.

DNA Methylation at Birth is Associated with Childhood Serum Immunoglobulin E Levels.

Immunoglobulin E (IgE) is known to play an important role in allergic diseases. Epigenetic traits acquired due to modification of deoxyribonucleic acid (DNA) methylation (DNAm) in early life may have phenotypic consequences through their role in transcriptional regulation with relevance to the developmental origins of diseases including allergy. However, epigenome-scale studies on the longitudinal association of cord blood DNAm with IgE over time are lacking. Our study aimed to examine the association of DNAm at birth with childhood serum IgE levels during early life. Genome-scale DNAm and total serum IgE measured at birth, 5, 8, and 11 years of children in the Taiwan Maternal and Infant Cohort Study were included in the study in the discovery stage. Linear mixed models were implemented to assess the association between cord blood DNAm at ~310K 5'-cytosine-phosphate-guanine-3' (CpG) sites with repeated IgE measurements, adjusting for cord blood IgE. Identified statistically significant CpGs (at a false discovery rate, FDR, of 0.05) were further tested in an independent replication cohort, the Isle of Wight (IoW) birth cohort. We mapped replicated CpGs to genes and conducted gene ontology analysis using ToppFun to identify significantly enriched pathways and biological processes of the genes. Cord blood DNAm of 273 CpG sites were significantly (FDR = 0.05) associated with IgE levels longitudinally. Among the identified CpGs available in both cohorts (184 CpGs), 92 CpGs (50%) were replicated in the IoW in terms of consistency in direction of associations between DNA methylation and IgE levels later in life, and 16 of the 92 CpGs showed statistically significant associations (P < .05). Gene ontology analysis identified 4 pathways (FDR = 0.05). The identified 16 CpG sites had the potential to serve as epigenetic markers associated with later IgE production, beneficial to allergic disease prevention and intervention.

Aberrant Extrafollicular B Cells, Immune Dysfunction, Myeloid Inflammation, and MyD88-Mutant Progenitors Precede Waldenstrom Macroglobulinemia.

Waldenstrom macroglobulinemia (WM) and its precursor IgM gammopathy are distinct disorders characterized by clonal mature IgM-expressing B-cell outgrowth in the bone marrow. Here, we show by high-dimensional single-cell immunogenomic profiling of patient samples that these disorders originate in the setting of global B-cell compartment alterations, characterized by expansion of genomically aberrant extrafollicular B cells of the nonmalignant clonotype. Alterations in the immune microenvironment preceding malignant clonal expansion include myeloid inflammation and naïve B- and T-cell depletion. Host response to these early lesions involves clone-specific T-cell immunity that may include MYD88 mutation-specific responses. Hematopoietic progenitors carry the oncogenic MYD88 mutations characteristic of the malignant WM clone. These data support a model for WM pathogenesis wherein oncogenic alterations and signaling in progenitors, myeloid inflammation, and global alterations in extrafollicular B cells create the milieu promoting extranodal pattern of growth in differentiated malignant cells. SIGNIFICANCE: These data provide evidence that growth of the malignant clone in WM is preceded by expansion of extrafollicular B cells, myeloid inflammation, and immune dysfunction in the preneoplastic phase. These changes may be related in part to MYD88 oncogenic signaling in pre-B progenitor cells and suggest a novel model for WM pathogenesis. This article is highlighted in the In This Issue feature, p. 549.

Subtype and grade-dependent spatial heterogeneity of T-cell infiltration in pediatric glioma.

Brain tumors are the leading cause of cancer-related mortality in children and have distinct genomic and molecular features compared with adult glioma. However, the properties of immune cells in these tumors has been vastly understudied compared with their adult counterparts. We combined multiplex immunofluorescence immunohistochemistry coupled with machine learning and single-cell mass cytometry to evaluate T-cells infiltrating pediatric glial tumors. We show that low-grade tumors are characterized by greater T-cell density compared with high-grade glioma (HGG). However, even among low-grade tumors, T-cell infiltration can be highly variable and subtype-dependent, with greater T-cell density in pleomorphic xanthoastrocytoma and ganglioglioma. CD3+ T-cell infiltration correlates inversely with the expression of SOX2, an embryonal stem cell marker commonly expressed by glial tumors. T-cells within both HGG and low-grade glioma (LGG) exhibit phenotypic heterogeneity and tissue-resident memory T-cells consist of distinct subsets of CD103+ and TCF1+ cells that exhibit distinct spatial localization patterns. TCF1+ T-cells are located closer to the vessels while CD103+ resident T-cells reside within the tumor further away from the vasculature. Recurrent tumors are characterized by a decline in CD103+ tumor-infiltrating T-cells. BRAFV600E mutation is immunogenic in children with LGG and may serve as a target for immune therapy. These data provide several novel insights into the subtype-dependent and grade-dependent changes in immune architecture in pediatric gliomas and suggest that harnessing tumor-resident T-cells may be essential to improve immune control in glioma.

Risk-associated alterations in marrow T cells in pediatric leukemia.

Current management of childhood leukemia is tailored based on disease risk determined by clinical features at presentation. Whether properties of the host immune response impact disease risk and outcome is not known. Here, we combine mass cytometry, single cell genomics, and functional studies to characterize the BM immune environment in children with B cell acute lymphoblastic leukemia and acute myelogenous leukemia at presentation. T cells in leukemia marrow demonstrate evidence of chronic immune activation and exhaustion/dysfunction, with attrition of naive T cells and TCF1+ stem-like memory T cells and accumulation of terminally differentiated effector T cells. Marrow-infiltrating NK cells also exhibit evidence of dysfunction, particularly in myeloid leukemia. Properties of immune cells identified distinct immune phenotype-based clusters correlating with disease risk in acute lymphoblastic leukemia. High-risk immune signatures were associated with expression of stem-like genes on tumor cells. These data provide a comprehensive assessment of the immune landscape of childhood leukemias and identify targets potentially amenable to therapeutic intervention. These studies also suggest that properties of the host response with depletion of naive T cells and accumulation of terminal-effector T cells may contribute to the biologic basis of disease risk. Properties of immune microenvironment identified here may also impact optimal application of immune therapies, including T cell-redirection approaches in childhood leukemia.

Hepatic Tumor Formation in Adult Mice Developmentally Exposed to Organotin.

BACKGROUND: Tributyltin (TBT) is a persistent and bioaccumulative environmental toxicant. Developmental exposure to TBT has been shown to cause fatty liver disease (steatosis), as well as increased adiposity in many species, leading to its characterization as an obesogen. OBJECTIVE: We aimed to determine the long-term effects of developmental TBT exposure on the liver. METHODS: C57BL/6J mice were exposed to a dose of TBT (0.5mg/kg body weight per day; 3.07µM) below the current developmental no observed adverse effect level (NOAEL) via drinking water, or drinking water alone, provided to the dam from preconception through lactation. Sires were exposed during breeding and lactation. Pups from two parity cycles were included in this study. Animals were followed longitudinally, and livers of offspring were analyzed by pathological evaluation, immunohistochemistry, immunoblotting, and RNA sequencing. RESULTS: Developmental exposure to TBT led to increased adiposity and hepatic steatosis at 14 and 20 weeks of age and increased liver adenomas at 45 weeks of age in male offspring. Female offspring displayed increased adiposity as compared with males, but TBT did not lead to an increase in fatty liver or tumor development in female offspring. Liver tumors in male mice were enriched in pathways and gene signatures associated with human and rodent nonalcoholic fatty liver disease (NAFLD) and hepatocellular carcinoma (HCC). This includes down-regulation of growth hormone receptor (GHR) and of STAT5 signaling, which occurred in response to TBT exposure and preceded liver tumor development. CONCLUSIONS: These data reveal a previously unappreciated ability of TBT to increase risk for liver tumorigenesis in mice in a sex-specific manner. Taken together, these findings provide new insights into how early life environmental exposures contribute to liver disease in adulthood. https://doi.org/10.1289/EHP5414.

Changes in DNA methylation from pre- to post-adolescence are associated with pubertal exposures.

BACKGROUND: Adolescence is a period characterized by major biological development, which may be associated with changes in DNA methylation (DNA-M). However, it is unknown to what extent DNA-M varies from pre- to post-adolescence, whether the pattern of changes is different between females and males, and how adolescence-related factors are associated with changes in DNA-M. METHODS: Genome-scale DNA-M at ages 10 and 18 years in whole blood of 325 subjects (n = 140 females) in the Isle of Wight (IOW) birth cohort was analyzed using Illumina Infinium arrays (450K and EPIC). Linear mixed models were used to examine DNA-M changes between pre- and post-adolescence and whether the changes were gender-specific. Adolescence-related factors and environmental exposure factors were assessed on their association with DNA-M changes. Replication of findings was attempted in the comparable Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. RESULTS: In the IOW cohort, after controlling for technical variation and cell compositions at both pre- and post-adolescence, 15,532 cytosine-phosphate-guanine (CpG) sites (of 400,825 CpGs, 3.88%) showed statistically significant DNA-M changes from pre-adolescence to post-adolescence invariant to gender (false discovery rate (FDR) = 0.05). Of these 15,532 CpGs, 10,212 CpGs (66%) were replicated in the ALSPAC cohort. Pathway analysis using Ingenuity Pathway Analysis (IPA) identified significant biological pathways related to growth and development of the reproductive system, emphasizing the importance of this period of transition on epigenetic state of genes. In addition, in IOW, we identified 1179 CpGs with gender-specific DNA-M changes. In the IOW cohort, body mass index (BMI) at age 10 years, age of growth spurt, nonsteroidal drugs use, and current smoking status showed statistically significant associations with DNA-M changes at 15 CpGs on 14 genes such as the AHRR gene. For BMI at age 10 years, the association was gender-specific. Findings on current smoking status were replicated in the ALSPAC cohort. CONCLUSION: Adolescent transition is associated with changes in DNA-M at more than 15K CpGs. Identified pathways emphasize the importance of this period of transition on epigenetic state of genes relevant to cell growth and immune system development.

Epigenome-wide association study of asthma and wheeze characterizes loci within HK1.

BACKGROUND: To identify novel epigenetic markers of adolescent asthma and replicate findings in an independent cohort, then explore whether such markers are detectable at birth, predictive of early-life wheeze, and associated with gene expression in cord blood. METHODS: We performed epigenome-wide screening with recursive random forest feature selection and internal validation in the IOW birth cohort. We then tested whether we could replicate these findings in the independent cohort ALSPAC and followed-up our top finding with children of the IOW cohort. RESULTS: We identified 10 CpG sites associated with adolescent asthma at a 5% false discovery rate (IOW, n = 370), five of which exhibited evidence of associations in the replication study (ALSPAC, n = 720). One site, cg16658191, within HK1 displayed particularly strong associations after cellular heterogeneity adjustments in both cohorts (ORIOW = 0.17, 95% CI 0.04-0.57) (ORALSPAC = 0.57, 95% CI 0.38-0.87). Additionally, higher expression of HK1 (OR = 3.81, 95% CI 1.41-11.77) in cord blood was predictive of wheezing in infancy (n = 82). CONCLUSION: We identified novel associations between asthma and wheeze with methylation at cg16658191 and the expression of HK1, which may serve as markers of, predictors of, and potentially etiologic factors involved in asthma and early life wheeze.