GM-CSF-dependent CD301b+ lung dendritic cells confer tolerance to inhaled allergens.

The severity of allergic asthma is driven by the balance between allergen-specific T regulatory (Treg) and T helper (Th)2 cells. However, it is unclear whether specific subsets of conventional dendritic cells (cDCs) promote the differentiation of these two T cell lineaeges. We have identified a subset of lung resident type 2 cDCs (cDC2s) that display high levels of CD301b and have potent Treg-inducing activity ex vivo. Single cell RNA sequencing and adoptive transfer experiments show that during allergic sensitization, many CD301b+ cDC2s transition in a stepwise manner to CD200+ cDC2s that selectively promote Th2 differentiation. GM-CSF augments the development and maintenance of CD301b+ cDC2s in vivo, and also selectively expands Treg-inducing CD301b+ cDC2s derived from bone marrow. Upon their adoptive transfer to recipient mice, lung-derived CD301b+ cDC2s confer immunological tolerance to inhaled allergens. Thus, GM-CSF maintains lung homeostasis by increasing numbers of Treg-inducing CD301b+ cDC2s.

GLIS3 expression in the thyroid gland in relation to TSH signaling and regulation of gene expression.

Loss of GLI-Similar 3 (GLIS3) function in mice and humans causes congenital hypothyroidism (CH). In this study, we demonstrate that GLIS3 protein is first detectable at E15.5 of murine thyroid development, a time at which GLIS3 target genes, such as Slc5a5 (Nis), become expressed. This, together with observations showing that ubiquitous Glis3KO mice do not display major changes in prenatal thyroid gland morphology, indicated that CH in Glis3KO mice is due to dyshormonogenesis rather than thyroid dysgenesis. Analysis of GLIS3 in postnatal thyroid suggested a link between GLIS3 protein expression and blood TSH levels. This was supported by data showing that treatment with TSH, cAMP, or adenylyl cyclase activators or expression of constitutively active PKA enhanced GLIS3 protein stability and transcriptional activity, indicating that GLIS3 activity is regulated at least in part by TSH/TSHR-mediated activation of PKA. The TSH-dependent increase in GLIS3 transcriptional activity would be critical for the induction of GLIS3 target gene expression, including several thyroid hormone (TH) biosynthetic genes, in thyroid follicular cells of mice fed a low iodine diet (LID) when blood TSH levels are highly elevated. Like TH biosynthetic genes, the expression of cell cycle genes is suppressed in ubiquitous Glis3KO mice fed a LID; however, in thyroid-specific Glis3 knockout mice, the expression of cell cycle genes was not repressed, in contrast to TH biosynthetic genes. This indicated that the inhibition of cell cycle genes in ubiquitous Glis3KO mice is dependent on changes in gene expression in GLIS3 target tissues other than the thyroid.

Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4.

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here, we determine the roles of CHD4 in enhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility. Its depletion leads to redistribution of transcription factors to previously unoccupied sites. During cellular reprogramming induced by the pioneer factor GATA3, CHD4 activity is necessary to prevent inappropriate chromatin opening. Mechanistically, CHD4 promotes nucleosome positioning over GATA3 binding motifs to compete with transcription factor-DNA interaction. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents unnecessary gene expression by editing chromatin binding activities of transcription factors.

Chemokine CCL19 promotes type 2 T-cell differentiation and allergic airway inflammation.

BACKGROUND: Allergic asthma is driven largely by allergen-specific TH2 cells, which develop in regional lymph nodes on the interaction of naive CD4+ T cells with allergen-bearing dendritic cells that migrate from the lung. This migration event is dependent on CCR7 and its chemokine ligand, CCL21. However, is has been unclear whether the other CCR7 ligand, CCL19, has a role in allergic airway disease. OBJECTIVE: This study sought to define the role of CCL19 in TH2 differentiation and allergic airway disease. METHODS: Ccl19-deficient mice were studied in an animal model of allergic asthma. Dendritic cells or fibroblastic reticular cells from wild-type and Ccl19-deficient mice were cultured with naive CD4+ T cells, and cytokine production was measured by ELISA. Recombinant CCL19 was added to CD4+ T-cell cultures, and gene expression was assessed by RNA-sequencing and quantitative PCR. Transcription factor activation was assessed by flow cytometry. RESULTS: Lungs of Ccl19-deficient mice had less allergic airway inflammation, reduced airway hyperresponsiveness, and less IL-4 and IL-13 production compared with lungs of Ccl19-sufficient animals. Naive CD4+ T cells cocultured with Ccl19-deficient dendritic cells or fibroblastic reticular cells produced lower amounts of type 2 cytokines than did T cells cocultured with their wild-type counterparts. Recombinant CCL19 increased phosphorylation of STAT5 and induced expression of genes associated with TH2 cell and IL-2 signaling pathways. CONCLUSIONS: These results reveal a novel, TH2 cell-inducing function of CCL19 in allergic airway disease and suggest that strategies to block this pathway might help to reduce the incidence or severity of allergic asthma.

Evaluating Stacked Methylation Markers for Blood-Based Multicancer Detection.

The ability to detect several types of cancer using a non-invasive, blood-based test holds the potential to revolutionize oncology screening. We mined tumor methylation array data from the Cancer Genome Atlas (TCGA) covering 14 cancer types and identified two novel, broadly-occurring methylation markers at TLX1 and GALR1. To evaluate their performance as a generalized blood-based screening approach, along with our previously reported methylation biomarker, ZNF154, we rigorously assessed each marker individually or combined. Utilizing TCGA methylation data and applying logistic regression models within each individual cancer type, we found that the three-marker combination significantly increased the average area under the ROC curve (AUC) across the 14 tumor types compared to single markers (p = 1.158 × 10-10; Friedman test). Furthermore, we simulated dilutions of tumor DNA into healthy blood cell DNA and demonstrated increased AUC of combined markers across all dilution levels. Finally, we evaluated assay performance in bisulfite sequenced DNA from patient tumors and plasma, including early-stage samples. When combining all three markers, the assay correctly identified nine out of nine lung cancer plasma samples. In patient plasma from hepatocellular carcinoma, ZNF154 alone yielded the highest combined sensitivity and specificity values averaging 68% and 72%, whereas multiple markers could achieve higher sensitivity or specificity, but not both. Altogether, this study presents a comprehensive pipeline for the identification, testing, and validation of multi-cancer methylation biomarkers with a considerable potential for detecting a broad range of cancer types in patient blood samples.

DNASE1L3 enhances antitumor immunity and suppresses tumor progression in colon cancer.

DNASE1L3, an enzyme highly expressed in DCs, is functionally important for regulating autoimmune responses to self-DNA and chromatin. Deficiency of DNASE1L3 leads to development of autoimmune diseases in both humans and mice. However, despite the well-established causal relationship between DNASE1L3 and immunity, little is known about the involvement of DNASE1L3 in regulation of antitumor immunity, the foundation of modern antitumor immunotherapy. In this study, we identify DNASE1L3 as a potentially new regulator of antitumor immunity and a tumor suppressor in colon cancer. In humans, DNASE1L3 is downregulated in tumor-infiltrating DCs, and this downregulation is associated with poor patient prognosis and reduced tumor immune cell infiltration in many cancer types. In mice, Dnase1l3 deficiency in the tumor microenvironment enhances tumor formation and growth in several colon cancer models. Notably, the increased tumor formation and growth in Dnase1l3-deficient mice are associated with impaired antitumor immunity, as evidenced by a substantial reduction of cytotoxic T cells and a unique subset of DCs. Consistently, Dnase1l3-deficient DCs directly modulate cytotoxic T cells in vitro. To our knowledge, our study unveils a previously unknown link between DNASE1L3 and antitumor immunity and further suggests that restoration of DNASE1L3 activity may represent a potential therapeutic approach for anticancer therapy.

Role of GLIS3 in thyroid development and in the regulation of gene expression in thyroid specific Glis3KO mice.

Loss of GLI-Similar 3 (GLIS3) function in mice and humans causes congenital hypothyroidism (CH). In this study, we demonstrate that GLIS3 protein is first detectable at E15.5 of murine thyroid development, a time when GLIS3 target genes, such as Slc5a5 (Nis), become also expressed. We further show that Glis3KO mice do not display any major changes in prenatal thyroid gland morphology indicating that CH in Glis3KO mice is due to dyshormonogenesis rather than thyroid dysgenesis. Analysis of thyroid-specific Glis3 knockout (Glis3-Pax8Cre) mice fed either a normal or low-iodine diet (ND or LID) revealed that, in contrast to ubiquitous Glis3KO mice, thyroid follicular cell proliferation and the expression of cell cycle genes were not repressed suggesting that the inhibition of thyroid follicular cell proliferation in ubiquitous Glis3KO mice is related to loss of GLIS3 function in other cell types. However, the expression of several thyroid hormone biosynthesis-, extracellular matrix (ECM)-, and inflammation-related genes was still suppressed in Glis3-Pax8Cre mice particularly under conditions of high blood levels of thyroid stimulating hormone (TSH). We further demonstrate that treatment with TSH, protein kinase A (PKA) or adenylyl cyclase activators or expression of constitutively active PKA enhances GLIS3 protein and activity, suggesting that GLIS3 transcriptional activity is regulated in part by TSH/TSHR-mediated activation of the PKA pathway. This mechanism of regulation provides an explanation for the dramatic increase in GLIS3 protein expression and the subsequent induction of GLIS3 target genes, including several thyroid hormone biosynthetic genes, in thyroid follicular cells of mice fed a LID.

Beadchip technology to detect DNA methylation in mouse faithfully recapitulates whole-genome bisulfite sequencing.

Aim: To facilitate wide-scale implementation of Illumina Mouse Methylation BeadChip (MMB) technology, array-based measurement of cytosine methylation was compared with the gold-standard assessment of DNA methylation by whole-genome bisulfite sequencing (WGBS). Methods: DNA methylation across two mouse strains (C57B6 and C3H) and both sexes was assessed using the MMB and compared with previously existing deep-coverage WGBS of mice of the same strain and sex. Results & conclusion: The findings demonstrated that 93.3-99.2% of sites had similar measurements of methylation across technologies and that differentially methylated cytosines and regions identified by each technology overlap and enrich for similar biological functions, suggesting that the MMB faithfully recapitulates the findings of WGBS.

Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4.

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here we determine the roles of CHD4 to enhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility at transcription factor binding sites; its depletion leads to altered motif scanning and redistribution of transcription factors to sites not previously occupied. During GATA3-mediated cellular reprogramming, CHD4 activity is necessary to prevent inappropriate chromatin opening and enhancer licensing. Mechanistically, CHD4 competes with transcription factor-DNA interaction by promoting nucleosome positioning over binding motifs. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents inappropriate gene expression by editing binding site selection by transcription factors.

GLIS3 regulates transcription of thyroid hormone biosynthetic genes in coordination with other thyroid transcription factors.

BACKGROUND: Loss of the transcription factor GLI-Similar 3 (GLIS3) function causes congenital hypothyroidism (CH) in both humans and mice due to decreased expression of several thyroid hormone (TH) biosynthetic genes in thyroid follicular cells. Whether and to what extent, GLIS3 regulates thyroid gene transcription in coordination with other thyroid transcriptional factors (TFs), such as PAX8, NKX2.1 and FOXE1, is poorly understood. METHODS: PAX8, NKX2.1, and FOXE1 ChIP-Seq analysis with mouse thyroid glands and rat thyrocyte PCCl3 cells was performed and compared to that of GLIS3 to analyze the co-regulation of gene transcription in thyroid follicular cells by these TFs. RESULTS: Analysis of the PAX8, NKX2.1, and FOXE1 cistromes identified extensive overlaps between these TF binding loci and those of GLIS3 indicating that GLIS3 shares many of the same regulatory regions with PAX8, NKX2.1, and FOXE1, particularly in genes associated with TH biosynthesis, induced by thyroid stimulating hormone (TSH), and suppressed in Glis3KO thyroid glands, including Slc5a5 (Nis), Slc26a4, Cdh16, and Adm2. ChIP-QPCR analysis showed that loss of GLIS3 did not significantly affect PAX8 or NKX2.1 binding and did not cause major alterations in H3K4me3 and H3K27me3 epigenetic signals. CONCLUSIONS: Our study indicates that GLIS3 regulates transcription of TH biosynthetic and TSH-inducible genes in thyroid follicular cells in coordination with PAX8, NKX2.1, and FOXE1 by binding within the same regulatory hub. GLIS3 does not cause major changes in chromatin structure at these common regulatory regions. GLIS3 may induce transcriptional activation by enhancing the interaction of these regulatory regions with other enhancers and/or RNA Polymerase II (Pol II) complexes.

Contribution of the Wolffian duct mesenchyme to the formation of the female reproductive tract.

The female reproductive tract develops from its embryonic precursor, the Müllerian duct. In close proximity to the Müllerian duct lies the precursor for the male reproductive tract, the Wolffian duct, which is eliminated in the female embryo during sexual differentiation. We discovered that a component of the Wolffian duct, its mesenchyme, is not eliminated after sexual differentiation. Instead, the Wolffian duct mesenchyme underwent changes in transcriptome and chromatin accessibility from male tract to female tract identity, and became a unique mesenchymal population in the female reproductive tract with localization and transcriptome distinct from the mesenchyme derived from the Müllerian duct. Partial ablation of the Wolffian duct mesenchyme stunted the growth of the fetal female reproductive tract in ex vivo organ culture. These findings reveal a new fetal origin of mesenchymal tissues for female reproductive tract formation and reshape our understanding of sexual differentiation of reproductive tracts.

Low-input ATAC&mRNA-seq protocol for simultaneous profiling of chromatin accessibility and gene expression.

We present a simple, fast, and robust protocol (low-input ATAC&mRNA-seq) to simultaneously generate ATAC-seq and mRNA-seq libraries from the same cells in limited cell numbers by coupling a simplified ATAC procedure using whole cells with a novel mRNA-seq approach that features a seamless on-bead process including direct mRNA isolation from the cell lysate, solid-phase cDNA synthesis, and direct tagmentation of mRNA/cDNA hybrids for library preparation. It enables dual-omics profiling from limited material when joint epigenome and transcriptome analyses are needed. For complete details on the use and execution of this protocol, please refer to Li et al. (2021).

CD11b+ lung dendritic cells at different stages of maturation induce Th17 or Th2 differentiation.

Dendritic cells (DC) in the lung that induce Th17 differentiation remain incompletely understood, in part because conventional CD11b+ DCs (cDC2) are heterogeneous. Here, we report a population of cDCs that rapidly accumulates in lungs of mice following house dust extract inhalation. These cells are Ly-6C+, are developmentally and phenotypically similar to cDC2, and strongly promote Th17 differentiation ex vivo. Single cell RNA-sequencing (scRNA-Seq) of lung cDC2 indicates 5 distinct clusters. Pseudotime analysis of scRNA-Seq data and adoptive transfer experiments with purified cDC2 subpopulations suggest stepwise developmental progression of immature Ly-6C+Ly-6A/E+ cDC2 to mature Ly-6C-CD301b+ lung resident cDC2 lacking Ccr7 expression, which then further mature into CD200+ migratory cDC2 expressing Ccr7. Partially mature Ly-6C+Ly-6A/E-CD301b- cDC2, which express Il1b, promote Th17 differentiation. By contrast, CD200+ mature cDC2 strongly induce Th2, but not Th17, differentiation. Thus, Th17 and Th2 differentiation are promoted by lung cDC2 at distinct stages of maturation.

A simple and robust method for simultaneous dual-omics profiling with limited numbers of cells.

Deciphering epigenetic regulation of gene expression requires measuring the epigenome and transcriptome jointly. Single-cell multi-omics technologies have been developed for concurrent profiling of chromatin accessibility and gene expression. However, multi-omics profiling of low-input bulk samples remains challenging. Therefore, we developed low-input ATAC&mRNA-seq, a simple and robust method for studying the role of chromatin structure in gene regulation in a single experiment with thousands of cells, to maximize insights from limited input material by obtaining ATAC-seq and mRNA-seq data simultaneously from the same cells with data quality comparable to that of conventional mono-omics assays. Integrative data analysis revealed similar strong association between promoter accessibility and gene expression when using the data of low-input ATAC&mRNA-seq as when using single-assay data, underscoring the accuracy and reliability of our dual-omics assay to generate both datum types simultaneously with just thousands of cells. We envision our method to be widely applied in many biological disciplines with limited materials.

GLIS1 regulates trabecular meshwork function and intraocular pressure and is associated with glaucoma in humans.

Chronically elevated intraocular pressure (IOP) is the major risk factor of primary open-angle glaucoma, a leading cause of blindness. Dysfunction of the trabecular meshwork (TM), which controls the outflow of aqueous humor (AqH) from the anterior chamber, is the major cause of elevated IOP. Here, we demonstrate that mice deficient in the Krüppel-like zinc finger transcriptional factor GLI-similar-1 (GLIS1) develop chronically elevated IOP. Magnetic resonance imaging and histopathological analysis reveal that deficiency in GLIS1 expression induces progressive degeneration of the TM, leading to inefficient AqH drainage from the anterior chamber and elevated IOP. Transcriptome and cistrome analyses identified several glaucoma- and extracellular matrix-associated genes as direct transcriptional targets of GLIS1. We also identified a significant association between GLIS1 variant rs941125 and glaucoma in humans (P = 4.73 × 10-6), further supporting a role for GLIS1 into glaucoma etiology. Our study identifies GLIS1 as a critical regulator of TM function and maintenance, AqH dynamics, and IOP.

DNMT1 reads heterochromatic H4K20me3 to reinforce LINE-1 DNA methylation.

DNA methylation and trimethylated histone H4 Lysine 20 (H4K20me3) constitute two important heterochromatin-enriched marks that frequently cooperate in silencing repetitive elements of the mammalian genome. However, it remains elusive how these two chromatin modifications crosstalk. Here, we report that DNA methyltransferase 1 (DNMT1) specifically 'recognizes' H4K20me3 via its first bromo-adjacent-homology domain (DNMT1BAH1). Engagement of DNMT1BAH1-H4K20me3 ensures heterochromatin targeting of DNMT1 and DNA methylation at LINE-1 retrotransposons, and cooperates with the previously reported readout of histone H3 tail modifications (i.e., H3K9me3 and H3 ubiquitylation) by the RFTS domain to allosterically regulate DNMT1's activity. Interplay between RFTS and BAH1 domains of DNMT1 profoundly impacts DNA methylation at both global and focal levels and genomic resistance to radiation-induced damage. Together, our study establishes a direct link between H4K20me3 and DNA methylation, providing a mechanism in which multivalent recognition of repressive histone modifications by DNMT1 ensures appropriate DNA methylation patterning and genomic stability.

UDP-glucose and P2Y14 receptor amplify allergen-induced airway eosinophilia.

Airway eosinophilia is a hallmark of allergic asthma and is associated with mucus production, airway hyperresponsiveness, and shortness of breath. Although glucocorticoids are widely used to treat asthma, their prolonged use is associated with several side effects. Furthermore, many individuals with eosinophilic asthma are resistant to glucocorticoid treatment, and they have an unmet need for novel therapies. Here, we show that UDP-glucose (UDP-G), a nucleotide sugar, is selectively released into the airways of allergen-sensitized mice upon their subsequent challenge with that same allergen. Mice lacking P2Y14R, the receptor for UDP-G, had decreased airway eosinophilia and airway hyperresponsiveness compared with wild-type mice in a protease-mediated model of asthma. P2Y14R was dispensable for allergic sensitization and for the production of type 2 cytokines in the lung after challenge. However, UDP-G increased chemokinesis in eosinophils and enhanced their response to the eosinophil chemoattractant, CCL24. In turn, eosinophils triggered the release of UDP-G into the airway, thereby amplifying eosinophilic recruitment. This positive feedback loop was sensitive to therapeutic intervention, as a small molecule antagonist of P2Y14R inhibited airway eosinophilia. These findings thus reveal a pathway that can be therapeutically targeted to treat asthma exacerbations and glucocorticoid-resistant forms of this disease.

CUT&Tag-BS for simultaneous profiling of histone modification and DNA methylation with high efficiency and low cost.

It remains a challenge to decipher the complex relationship between DNA methylation, histone modification, and the underlying DNA sequence with limited input material. Here, we developed an efficient, low-input, and low-cost method for the simultaneous profiling of genomic localization of histone modification and methylation status of the underlying DNA at single-base resolution from the same cells in a single experiment by integrating cleavage under targets and tagmentation (CUT&Tag) with tagmentation-based bisulfite sequencing (CUT&Tag-BS). We demonstrated the validity of our method using representative histone modifications of euchromatin and constitutive and facultative heterochromatin (H3K4me1, H3K9me3, and H3K27me3, respectively). Similar histone modification enrichment patterns were observed in CUT&Tag-BS compared with non-bisulfite-treated control, and H3K4me1-marked regions were found to mostly be CpG poor, lack methylation concordance, and exhibit prevalent DNA methylation heterogeneity among mouse embryonic stem cells (mESCs). We anticipate that CUT&Tag-BS will be widely applied to directly address the genomic relationship between DNA methylation and histone modification, especially in low-input scenarios with precious biological samples.

The genomic regulatory elements for estrogen receptor alpha transactivation-function-1 regulated genes.

Estrogen receptor alpha (ERα) is a ligand-dependent transcription regulator, containing two transactivation functional domains, AF-1 and AF-2. The selective estrogen receptor modulators (SERMs), including 4-hydroxytamoxifen (4OHT), activate AF-1 preferentially rather than AF-2. However, it is unclear whether this specific function is related to the tissue-selective functionality of SERMs. Moreover, there is no information determining AF-1-dependent estrogenic-genes existing in tissues. We sought to identify AF-1-dependent estrogenic-genes using the AF-2 mutated knock-in (KI) mouse model, AF2ERKI. AF2ER is an AF-2 disrupted estradiol (E2)-insensitive mutant ERα, but AF-1-dependent transcription can be activated by the estrogen-antagonists, fulvestrant (ICI) and 4OHT. Gene profiling and ChIP-Seq analysis identified Klk1b21 as an ICI-inducible gene in AF2ERKI uterus. The regulatory activity was analyzed further using a cell-based reporter assay. The 5'-flanking 0.4k bp region of Klk1b21 gene responded as an ERα AF-1-dependent estrogen-responsive promoter. The 150 bp minimum ERα binding element (EBE) consists of three direct repeats. These three half-site sequences were essential for the ERα-dependent transactivation and were differentially recognized by E2 and 4OHT for the gene activation. This response was impaired when the minimum EBE was fused with a thymidine-kinase promoter but could be restored by fusion with the 100 bp minimum transcription initiation element (TIE) of Klk1b21, suggesting that the cooperative function of EBE and TIE is essential for mediating AF-1-dependent transactivation. These findings provide the first in vivo evidence that endogenous ERα AF-1 dominant estrogenic-genes exist in estrogen-responsive organs. Such findings will aid in understanding the mechanism of ERα-dependent tissue-selective activity of SERMs.

Single Nucleotide Resolution Analysis Reveals Pervasive, Long-Lasting DNA Methylation Changes by Developmental Exposure to a Mitochondrial Toxicant.

Mitochondrial-driven alterations of the epigenome have been reported, but whether they are relevant at the organismal level remains unknown. The viable yellow agouti mouse (Avy) is a powerful epigenetic biosensor model that reports on the DNA methylation status of the Avy locus, which is established prior to the three-germ-layer separation, through the coat color of the animals. Here we show that maternal exposure to rotenone, a potent mitochondrial complex I inhibitor, not only changes the DNA methylation status of the Avy locus in the skin but broadly affects the liver DNA methylome of the offspring. These effects are accompanied by altered gene expression programs that persist throughout life, and which associate with impairment of antioxidant activity and mitochondrial function in aged animals. These pervasive and lasting genomic effects suggest a putative role for mitochondria in regulating life-long gene expression programs through developmental nuclear epigenetic remodeling.