Epigenomic Profiling of Human CD4+ T Cells Supports a Linear Differentiation Model and Highlights Molecular Regulators of Memory Development.

The impact of epigenetics on the differentiation of memory T (Tmem) cells is poorly defined. We generated deep epigenomes comprising genome-wide profiles of DNA methylation, histone modifications, DNA accessibility, and coding and non-coding RNA expression in naive, central-, effector-, and terminally differentiated CD45RA+ CD4+ Tmem cells from blood and CD69+ Tmem cells from bone marrow (BM-Tmem). We observed a progressive and proliferation-associated global loss of DNA methylation in heterochromatic parts of the genome during Tmem cell differentiation. Furthermore, distinct gradually changing signatures in the epigenome and the transcriptome supported a linear model of memory development in circulating T cells, while tissue-resident BM-Tmem branched off with a unique epigenetic profile. Integrative analyses identified candidate master regulators of Tmem cell differentiation, including the transcription factor FOXP1. This study highlights the importance of epigenomic changes for Tmem cell biology and demonstrates the value of epigenetic data for the identification of lineage regulators.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states.

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

Resident memory CD4+ T lymphocytes mobilize from bone marrow to contribute to a systemic secondary immune reaction.

Resident memory T lymphocytes (TRM ) of epithelial tissues and the Bm protect their host tissue. To what extent these cells are mobilized and contribute to systemic immune reactions is less clear. Here, we show that in secondary immune reactions to the measles-mumps-rubella (MMR) vaccine, CD4+ TRM are mobilized into the blood within 16 to 48 h after immunization in humans. This mobilization of TRM is cognate: TRM recognizing other antigens are not mobilized, unless they cross-react with the vaccine. We also demonstrate through methylome analyses that TRM are mobilized from the Bm. These mobilized cells make significant contribution to the systemic immune reaction, as evidenced by their T-cell receptor Vß clonotypes represented among the newly generated circulating memory T-cells, 14 days after vaccination. Thus, TRM of the Bm confer not only local, but also systemic immune memory.

Epigenetic control of region-specific transcriptional programs in mouse cerebellar and cortical astrocytes.

Astrocytes from the cerebral cortex (CTX) and cerebellum (CB) share basic molecular programs, but also form distinct spatial and functional subtypes. The regulatory epigenetic layers controlling such regional diversity have not been comprehensively investigated so far. Here, we present an integrated epigenome analysis of methylomes, open chromatin, and transcriptomes of astroglia populations isolated from the cortex or cerebellum of young adult mice. Besides a basic overall similarity in their epigenomic programs, cortical astrocytes and cerebellar astrocytes exhibit substantial differences in their overall open chromatin structure and in gene-specific DNA methylation. Regional epigenetic differences are linked to differences in transcriptional programs encompassing genes of region-specific transcription factor networks centered around Lhx2/Foxg1 in CTX astrocytes and the Zic/Irx families in CB astrocytes. The distinct epigenetic signatures around these transcription factor networks point to a complex interconnected and combinatorial regulation of region-specific transcriptomes. These findings suggest that key transcription factors, previously linked to temporal, regional, and spatial control of neurogenesis, also form combinatorial networks important for astrocytes. Our study provides a valuable resource for the molecular basis of regional astrocyte identity and physiology.

Differential degradation of RNA species by autophagy-related pathways in Arabidopsis.

The plant vacuole recycles proteins and RNA delivered to it by autophagy. In this study, by isolating intact vacuoles from Arabidopsis plants, followed by subsequent RNA purification, and deep sequencing, we provide a comprehensive characterization of Arabidopsis vacuolar RNAome. In the vacuolar RNAome, we detected ribosomal RNAs, transfer RNAs, including those of chloroplast origin, and in addition small RNA types. As autophagy is a main mechanism for the transport of RNA to the vacuole, atg5-1 mutants deficient in autophagy were included in our analysis. We observed severely reduced amounts of most chloroplast-derived RNA species in these mutants. Comparisons with cellular RNA composition provided an indication of possible up-regulation of alternative RNA breakdown pathways. By contrast, vacuolar RNA processing and composition in plants lacking vacuolar ribonuclease 2, involved in cellular RNA homeostasis, only showed minor alterations, possibly because of the presence of further so far unknown vacuolar RNase species. Among the small RNA types, we detected mature miRNAs in all vacuolar preparations but at much lower frequency in atg5-1, raising the possibility of a biological role for vacuolar miRNAs.

Dysregulation of COVID-19 related gene expression in the COPD lung.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) patients are at increased risk of poor outcome from Coronavirus disease (COVID-19). Early data suggest elevated Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) receptor angiotensin converting enzyme 2 (ACE2) expression, but relationships to disease phenotype and downstream regulators of inflammation in the Renin-Angiotensin system (RAS) are unknown. We aimed to determine the relationship between RAS gene expression relevant to SARS-CoV-2 infection in the lung with disease characteristics in COPD, and the regulation of newly identified SARS-CoV-2 receptors and spike-cleaving proteases, important for SARS-CoV-2 infection. METHODS: We quantified gene expression using RNA sequencing of epithelial brushings and bronchial biopsies from 31 COPD and 37 control subjects. RESULTS: ACE2 gene expression (log2-fold change (FC)) was increased in COPD compared to ex-smoking (HV-ES) controls in epithelial brushings (0.25, p = 0.042) and bronchial biopsies (0.23, p = 0.050), and correlated with worse lung function (r = - 0.28, p = 0.0090). ACE2 was further increased in frequent exacerbators compared to infrequent exacerbators (0.51, p = 0.00045) and associated with use of ACE inhibitors (ACEi) (0.50, p = 0.0034), having cardiovascular disease (0.23, p = 0.048) or hypertension (0.34, p = 0.0089), and inhaled corticosteroid use in COPD subjects in bronchial biopsies (0.33, p = 0.049). Angiotensin II receptor type (AGTR)1 and 2 expression was decreased in COPD bronchial biopsies compared to HV-ES controls with log2FC of -0.26 (p = 0.033) and - 0.40, (p = 0.0010), respectively. However, the AGTR1:2 ratio was increased in COPD subjects compared with HV-ES controls, log2FC of 0.57 (p = 0.0051). Basigin, a newly identified potential SARS-CoV-2 receptor was also upregulated in both brushes, log2FC of 0.17 (p = 0.0040), and bronchial biopsies, (log2FC of 0.18 (p = 0.017), in COPD vs HV-ES. Transmembrane protease, serine (TMPRSS)2 was not differentially regulated between control and COPD. However, various other spike-cleaving proteases were, including TMPRSS4 and Cathepsin B, in both epithelial brushes (log2FC of 0.25 (p = 0.0012) and log2FC of 0.56 (p = 5.49E-06), respectively) and bronchial biopsies (log2FC of 0.49 (p = 0.00021) and log2FC of 0.246 (p = 0.028), respectively). CONCLUSION: This study identifies key differences in expression of genes related to susceptibility and aetiology of COVID-19 within the COPD lung. Further studies to understand the impact on clinical course of disease are now required.

Neuronal deficiency of p38α-MAPK ameliorates symptoms and pathology of APP or Tau-transgenic Alzheimer's mouse models.

Alzheimer's disease (AD) is the leading cause of dementia with very limited therapeutic options. Amyloid ß (Aß) and phosphorylated Tau (p-Tau) are key pathogenic molecules in AD. P38α-MAPK is specifically activated in AD lesion sites. However, its effects on AD pathogenesis, especially on p-Tau-associated brain pathology, and the underlying molecular mechanisms remain unclear. We mated human APP-transgenic mice and human P301S Tau-transgenic mice with mapk14-floxed and neuron-specific Cre-knock-in mice. We observed that deletion of p38α-MAPK specifically in neurons improves the cognitive function of both 9-month-old APP and Tau-transgenic AD mice, which is associated with decreased Aß and p-Tau load in the brain. We further used next-generation sequencing to analyze the gene transcription in brains of p38α-MAPK deficient and wild-type APP-transgenic mice, which indicated that deletion of p38α-MAPK regulates the transcription of calcium homeostasis-related genes, especially downregulates the expression of grin2a, a gene encoding NMDAR subunit NR2A. Cell culture experiments further verified that deletion of p38α-MAPK inhibits NMDA-triggered calcium influx and neuronal apoptosis. Our systemic studies of AD pathogenic mechanisms using both APP- and Tau-transgenic mice suggested that deletion of neuronal p38α-MAPK attenuates AD-associated brain pathology and protects neurons in AD pathogenesis. This study supports p38α-MAPK as a novel target for AD therapy.

Exogenous RNAi mechanisms contribute to transcriptome adaptation by phased siRNA clusters in Paramecium.

Extensive research has characterized distinct exogenous RNAi pathways interfering in gene expression during vegetative growth of the unicellular model ciliate Paramecium. However, role of RNAi in endogenous transcriptome regulation, and environmental adaptation is unknown. Here, we describe the first genome-wide profiling of endogenous sRNAs in context of different transcriptomic states (serotypes). We developed a pipeline to identify, and characterize 2602 siRNA producing clusters (SRCs). Our data show no evidence that SRCs produce miRNAs, and in contrast to other species, no preference for strand specificity of siRNAs. Interestingly, most SRCs overlap coding genes and a separate group show siRNA phasing along the entire open reading frame, suggesting that the mRNA transcript serves as a source for siRNAs. Integrative analysis of siRNA abundance and gene expression levels revealed surprisingly that mRNA and siRNA show negative as well as positive associations. Two RNA-dependent RNA Polymerase mutants, RDR1 and RDR2, show a drastic loss of siRNAs especially in phased SRCs accompanied with increased mRNA levels. Importantly, most SRCs depend on both RDRs, reminiscent to primary siRNAs in the RNAi against exogenous RNA, indicating mechanistic overlaps between exogenous and endogenous RNAi contributing to flexible transcriptome adaptation.

Unique and assay specific features of NOMe-, ATAC- and DNase I-seq data.

Chromatin accessibility maps are important for the functional interpretation of the genome. Here, we systematically analysed assay specific differences between DNase I-seq, ATAC-seq and NOMe-seq in a side by side experimental and bioinformatic setup. We observe that most prominent nucleosome depleted regions (NDRs, e.g. in promoters) are roboustly called by all three or at least two assays. However, we also find a high proportion of assay specific NDRs that are often 'called' by only one of the assays. We show evidence that these assay specific NDRs are indeed genuine open chromatin sites and contribute important information for accurate gene expression prediction. While technically ATAC-seq and DNase I-seq provide a superb high NDR calling rate for relatively low sequencing costs in comparison to NOMe-seq, NOMe-seq singles out for its genome-wide coverage allowing to not only detect NDRs but also endogenous DNA methylation and as we show here genome wide segmentation into heterochromatic B domains and local phasing of nucleosomes outside of NDRs. In summary, our comparisons strongly suggest to consider assay specific differences for the experimental design and for generalized and comparative functional interpretations.

Genome-Wide Analysis of the Nucleosome Landscape in Individuals with Coffin-Siris Syndrome.

The switch/sucrose non-fermenting (SWI/SNF) complex is an ATP-dependent chromatin remodeller that regulates the spacing of nucleosomes and thereby controls gene expression. Heterozygous mutations in genes encoding subunits of the SWI/SNF complex have been reported in individuals with Coffin-Siris syndrome (CSS), with the majority of the mutations in ARID1B. CSS is a rare congenital disorder characterized by facial dysmorphisms, digital anomalies, and variable intellectual disability. We hypothesized that mutations in genes encoding subunits of the ubiquitously expressed SWI/SNF complex may lead to alterations of the nucleosome profiles in different cell types. We performed the first study on CSS-patient samples and investigated the nucleosome landscapes of cell-free DNA (cfDNA) isolated from blood plasma by whole-genome sequencing. In addition, we studied the nucleosome landscapes of CD14+ monocytes from CSS-affected individuals by nucleosome occupancy and methylome-sequencing (NOMe-seq) as well as their expression profiles. In cfDNA of CSS-affected individuals with heterozygous ARID1B mutations, we did not observe major changes in the nucleosome profile around transcription start sites. In CD14+ monocytes, we found few genomic regions with different nucleosome occupancy when compared to controls. RNA-seq analysis of CD14+ monocytes of these individuals detected only few differentially expressed genes, which were not in proximity to any of the identified differential nucleosome-depleted regions. In conclusion, we show that heterozygous mutations in the human SWI/SNF subunit ARID1B do not have a major impact on the nucleosome landscape or gene expression in blood cells. This might be due to functional redundancy, cell-type specificity, or alternative functions of ARID1B.

Expression of retinoic acid signaling components ADH7 and ALDH1A1 is reduced in aniridia limbal epithelial cells and a siRNA primary cell based aniridia model.

PAX6-related Aniridia is a sight-threatening disease involving progression of secondary glaucoma and aniridia related keratopathy (ARK). Change or loss of limbal epithelial progenitors causes epithelial surface defects. We analyzed the effect of PAX6 on mRNA expression changes with a two-step approach, as follows. First, we sequenced mRNA from limbal epithelial cells isolated from controls and aniridia patients. Second, we confirmed the bioinformatics and literature-based result list for a small interfering RNA (siRNA)-based primary aniridia cell model (PAX6 knockdown). With this approach, we expected that the genes directly influenced by PAX6 would be distinguishable from those affected secondarily by the ARK disease state. Therefore, epithelial cells were isolated from the limbus region of two patients with aniridia and cultured in keratinocyte serum-free medium. Normal control cells were obtained from the limbus region of corneal donors. For the siRNA-based aniridia cell model, cells were transfected with Lipofectamine and 5 nM siRNA against PAX6 or control treatment. All cells were lysed to yield DNA, RNA, and protein. Reduction of PAX6 protein was assessed by western blot. Aniridia and control Poly-A-enriched RNA libraries were subjected to next-generation sequencing. The differential analysis was a combination of quantification with RSEM and differential tests with edgeR. Gene lists were filtered by comparison to NCBI GEO datasets, annotated with DAVID, and manually annotated using a literature search. Based on the resulting filtered gene list, qPCR primers were purchased, and candidate genes (TP63, ABCG2, ADH7, ALDH1A1, PITX1, DKK1, DSG1, KRT12, KRT3, KRT13, SPINK6, SPINK7, CTSV, SERPINB1) were verified by qPCR on the siRNA-based aniridia cell model. We identified genes that might be regulated by PAX6 and showed that SPINK7 mRNA, which codes for a protease inhibitor, is downregulated in patients as well as in our primary aniridia cell model. ALDH1A1 and AHD7 mRNA levels were reduced in limbal epithelial cells of aniridia patients, and both transcripts were downregulated by PAX6 knockdown in our cell model. This siRNA-based aniridia cell model is a valuable tool for confirming identified PAX6-affected genes that might promote ARK pathogenesis. The model recapitulated expression changes for SPINK7, ADH7, and ALDH1A1 that were also observed in patient samples. These results provide evidence that PAX6 might drive corneal epithelial differentiation by direct or indirect control of retinoic acid signaling processes through ADH7 and ALDH1A1.

Combining transcription factor binding affinities with open-chromatin data for accurate gene expression prediction.

The binding and contribution of transcription factors (TF) to cell specific gene expression is often deduced from open-chromatin measurements to avoid costly TF ChIP-seq assays. Thus, it is important to develop computational methods for accurate TF binding prediction in open-chromatin regions (OCRs). Here, we report a novel segmentation-based method, TEPIC, to predict TF binding by combining sets of OCRs with position weight matrices. TEPIC can be applied to various open-chromatin data, e.g. DNaseI-seq and NOMe-seq. Additionally, Histone-Marks (HMs) can be used to identify candidate TF binding sites. TEPIC computes TF affinities and uses open-chromatin/HM signal intensity as quantitative measures of TF binding strength. Using machine learning, we find low affinity binding sites to improve our ability to explain gene expression variability compared to the standard presence/absence classification of binding sites. Further, we show that both footprints and peaks capture essential TF binding events and lead to a good prediction performance. In our application, gene-based scores computed by TEPIC with one open-chromatin assay nearly reach the quality of several TF ChIP-seq data sets. Finally, these scores correctly predict known transcriptional regulators as illustrated by the application to novel DNaseI-seq and NOMe-seq data for primary human hepatocytes and CD4+ T-cells, respectively.

Two sets of RNAi components are required for heterochromatin formation in trans triggered by truncated transgenes.

Across kingdoms, RNA interference (RNAi) has been shown to control gene expression at the transcriptional- or the post-transcriptional level. Here, we describe a mechanism which involves both aspects: truncated transgenes, which fail to produce intact mRNA, induce siRNA accumulation and silencing of homologous loci in trans in the ciliate Paramecium We show that silencing is achieved by co-transcriptional silencing, associated with repressive histone marks at the endogenous gene. This is accompanied by secondary siRNA accumulation, strictly limited to the open reading frame of the remote locus. Our data shows that in this mechanism, heterochromatic marks depend on a variety of RNAi components. These include RDR3 and PTIWI14 as well as a second set of components, which are also involved in post-transcriptional silencing: RDR2, PTIWI13, DCR1 and CID2. Our data indicates differential processing of nascent un-spliced and long, spliced transcripts thus suggesting a hitherto-unrecognized functional interaction between post-transcriptional and co-transcriptional RNAi. Both sets of RNAi components are required for efficient trans-acting RNAi at the chromatin level and our data indicates similar mechanisms contributing to genome wide regulation of gene expression by epigenetic mechanisms.

Evening expression of arabidopsis GIGANTEA is controlled by combinatorial interactions among evolutionarily conserved regulatory motifs.

Diurnal patterns of gene transcription are often conferred by complex interactions between circadian clock control and acute responses to environmental cues. Arabidopsis thaliana GIGANTEA (GI) contributes to photoperiodic flowering, circadian clock control, and photoreceptor signaling, and its transcription is regulated by the circadian clock and light. We used phylogenetic shadowing to identify three evolutionarily constrained regions (conserved regulatory modules [CRMs]) within the GI promoter and show that CRM2 is sufficient to confer a similar transcriptional pattern as the full-length promoter. Dissection of CRM2 showed that one subfragment (CRM2-A) contributes light inducibility, while another (CRM2-B) exhibits a diurnal response. Mutational analysis showed that three ABA RESPONSE ELEMENT LIKE (ABREL) motifs in CRM2-A and three EVENING ELEMENTs (EEs) in CRM2-B are essential in combination to confer a high amplitude diurnal pattern of expression. Genome-wide analysis identified characteristic spacing patterns of EEs and 71 A. thaliana promoters containing three EEs. Among these promoters, that of FLAVIN BINDING KELCH REPEAT F-BOX1 was analyzed in detail and shown to harbor a CRM functionally related to GI CRM2. Thus, combinatorial interactions among EEs and ABRELs confer diurnal patterns of transcription via an evolutionarily conserved module present in GI and other evening-expressed genes.

Human cytomegalovirus drives epigenetic imprinting of the IFNG locus in NKG2Chi natural killer cells.

Memory type 1 T helper (T(H)1) cells are characterized by the stable expression of interferon (IFN)-γ as well as by the epigenetic imprinting of the IFNG locus. Among innate cells, NK cells play a crucial role in the defense against cytomegalovirus (CMV) and represent the main source of IFN-γ. Recently, it was shown that memory-like features can be observed in NK cell subsets after CMV infection. However, the molecular mechanisms underlying NK cell adaptive properties have not been completely defined. In the present study, we demonstrated that only NKG2Chi NK cells expanded in human CMV (HCMV) seropositive individuals underwent epigenetic remodeling of the IFNG conserved non-coding sequence (CNS) 1, similar to memory CD8(+) T cells or T(H)1 cells. The accessibility of the CNS1 was required to enhance IFN-γ transcriptional activity in response to NKG2C and 2B4 engagement, which led to consistent IFN-γ production in NKG2C(hi) NK cells. Thus, our data identify epigenetic imprinting of the IFNG locus as selective hallmark and crucial mechanism driving strong and stable IFN-γ expression in HCMV-specific NK cell expansions, providing a molecular basis for the regulation of adaptive features in innate cells.

Elevated salicylic acid levels conferred by increased expression of ISOCHORISMATE SYNTHASE 1 contribute to hyperaccumulation of SUMO1 conjugates in the Arabidopsis mutant early in short days 4.

Post-translational modification of proteins by attachment of small ubiquitin-like modifier (SUMO) is essential for plant growth and development. Mutations in the SUMO protease early in short days 4 (ESD4) cause hyperaccumulation of conjugates formed between SUMO and its substrates, and phenotypically are associated with extreme early flowering and impaired growth. We performed a suppressor mutagenesis screen of esd4 and identified a series of mutants called suppressor of esd4 (sed), which delay flowering, enhance growth and reduce hyperaccumulation of SUMO conjugates. Genetic mapping and genome sequencing indicated that one of these mutations (sed111) is in the gene salicylic acid induction-deficient 2 (SID2), which encodes ISOCHORISMATE SYNTHASE I, an enzyme required for biosynthesis of salicylic acid (SA). Analyses showed that compared with wild-type plants, esd4 contains higher levels of SID2 mRNA and about threefold more SA, whereas sed111 contains lower SA levels. Other sed mutants also contain lower SA levels but are not mutant for SID2, although most reduce SID2 mRNA levels. Therefore, higher SA levels contribute to the small size, early flowering and elevated SUMO conjugate levels of esd4. Our results support previous data indicating that SUMO homeostasis influences SA biosynthesis in wild-type plants, and also demonstrate that elevated levels of SA strongly increase the abundance of SUMO conjugates.

Analysis of the Arabidopsis shoot meristem transcriptome during floral transition identifies distinct regulatory patterns and a leucine-rich repeat protein that promotes flowering.

Flowering of Arabidopsis thaliana is induced by exposure to long days (LDs). During this process, the shoot apical meristem is converted to an inflorescence meristem that forms flowers, and this transition is maintained even if plants are returned to short days (SDs). We show that exposure to five LDs is sufficient to commit the meristem of SD-grown plants to flower as if they were exposed to continuous LDs. The MADS box proteins SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) and FRUITFULL (FUL) play essential roles in this commitment process and in the induction of flowering downstream of the transmissible FLOWERING LOCUS T (FT) signal. We exploited laser microdissection and Solexa sequencing to identify 202 genes whose transcripts increase in the meristem during floral commitment. Expression of six of these transcripts was tested in different mutants, allowing them to be assigned to FT-dependent or FT-independent pathways. Most, but not all, of those dependent on FT and its paralog TWIN SISTER OF FT (TSF) also relied on SOC1 and FUL. However, this dependency on FT and TSF or SOC1 and FUL was often bypassed in the presence of the short vegetative phase mutation. FLOR1, which encodes a leucine-rich repeat protein, was induced in the early inflorescence meristem, and flor1 mutations delayed flowering. Our data contribute to the definition of LD-dependent pathways downstream and in parallel to FT.

Synteny-based mapping-by-sequencing enabled by targeted enrichment.

Mapping-by-sequencing, as implemented in SHOREmap ('SHOREmapping'), is greatly accelerating the identification of causal mutations. The original SHOREmap approach based on resequencing of bulked segregants required a highly accurate and complete reference sequence. However, current whole-genome or transcriptome assemblies from next-generation sequencing data of non-model organisms do not produce chromosome-length scaffolds. We have therefore developed a method that exploits synteny with a related genome for genetic mapping. We first demonstrate how mapping-by-sequencing can be performed using a reduced number of markers, and how the associated decrease in the number of markers can be compensated for by enrichment of marker sequences. As proof of concept, we apply this method to Arabidopsis thaliana gene models ordered by synteny with the genome sequence of the distant relative Brassica rapa, whose genome has several large-scale rearrangements relative to A. thaliana. Our approach provides an alternative method for high-resolution genetic mapping in species that lack finished genome reference sequences or for which only RNA-seq assemblies are available. Finally, for improved identification of causal mutations by fine-mapping, we introduce a new likelihood ratio test statistic, transforming local allele frequency estimations into a confidence interval similar to conventional mapping intervals.

Molecular control of endurance training adaptation in male mouse skeletal muscle.

Skeletal muscle has an enormous plastic potential to adapt to various external and internal perturbations. Although morphological changes in endurance-trained muscles are well described, the molecular underpinnings of training adaptation are poorly understood. We therefore aimed to elucidate the molecular signature of muscles of trained male mice and unravel the training status-dependent responses to an acute bout of exercise. Our results reveal that, even though at baseline an unexpectedly low number of genes define the trained muscle, training status substantially affects the transcriptional response to an acute challenge, both quantitatively and qualitatively, in part associated with epigenetic modifications. Finally, transiently activated factors such as the peroxisome proliferator-activated receptor-γ coactivator 1α are indispensable for normal training adaptation. Together, these results provide a molecular framework of the temporal and training status-dependent exercise response that underpins muscle plasticity in training.

Alterations in the hepatocyte epigenetic landscape in steatosis.

Fatty liver disease or the accumulation of fat in the liver, has been reported to affect the global population. This comes with an increased risk for the development of fibrosis, cirrhosis, and hepatocellular carcinoma. Yet, little is known about the effects of a diet containing high fat and alcohol towards epigenetic aging, with respect to changes in transcriptional and epigenomic profiles. In this study, we took up a multi-omics approach and integrated gene expression, methylation signals, and chromatin signals to study the epigenomic effects of a high-fat and alcohol-containing diet on mouse hepatocytes. We identified four relevant gene network clusters that were associated with relevant pathways that promote steatosis. Using a machine learning approach, we predict specific transcription factors that might be responsible to modulate the functionally relevant clusters. Finally, we discover four additional CpG loci and validate aging-related differential CpG methylation. Differential CpG methylation linked to aging showed minimal overlap with altered methylation in steatosis.