The alternative role of DNA methylation in splicing regulation

Epigenetic Dysregulation in Alcohol-Associated Behaviors: Preclinical and Clinical Evidence

Alcohol use disorder (AUD) is characterized by loss of control over intake and drinking despite harmful consequences. At a molecular level, AUD is associated with long-term neuroadaptations in key brain regions that are involved in reward processing and decision-making. Over the last decades, a great effort has been made to understand the neurobiological basis underlying AUD. Epigenetic mechanisms have emerged as an important mechanism in the regulation of long-term alcohol-induced gene expression changes. Here, we review the literature supporting a role for epigenetic processes in AUD. We particularly focused on the three most studied epigenetic mechanisms: DNA methylation, Histone modification and non-coding RNAs. Clinical studies indicate an association between AUD and DNA methylation both at the gene and global levels. Using behavioral paradigms that mimic some of the characteristics of AUD, preclinical studies demonstrate that changes in epigenetic mechanisms can functionally impact alcohol-associated behaviors. While many studies support a therapeutic potential for targeting epigenetic enzymes, more research is needed to fully understand their role in AUD. Identification of brain circuits underlying alcohol-associated behaviors has made major advances in recent years. However, there are very few studies that investigate how epigenetic mechanisms can affect these circuits or impact the neuronal ensembles that promote alcohol-associated behaviors. Studies that focus on the role of circuit-specific and cell-specific epigenetic changes for clinically relevant alcohol behaviors may provide new insights on the functional role of epigenetic processes in AUD.

Sex-specific changes in autosomal methylation rate in ageing common terns

Senescence, an age-related decline in survival and/or reproductive performance, occurs in species across the tree of life. Molecular mechanisms underlying this within-individual phenomenon are still largely unknown, but DNA methylation changes with age are among the candidates. Using a longitudinal approach, we investigated age-specific changes in autosomal methylation of common terns, relatively long-lived migratory seabirds known to show senescence. We collected blood at 1-, 3- and/or 4-year intervals, extracted DNA from the erythrocytes and estimated autosomal DNA methylation by mapping Reduced Representative Bisulfite Sequencing reads to a de novo assembled reference genome. We found autosomal methylation levels to decrease with age within females, but not males, and no evidence for selective (dis)appearance of birds of either sex in relation to their methylation level. Moreover, although we found positions in the genome to consistently vary in their methylation levels, individuals did not show such strong consistent variance. These results pave the way for studies at the level of genome features or specific positions, which should elucidate the functional consequences of the patterns observed, and how they translate to the ageing phenotype.

Patient-specific identification of genome-wide DNA-methylation differences between intracranial and extracranial melanoma metastases

Melanomas frequently metastasize to distant organs and especially intracranial metastases still represent a major clinical challenge. Epigenetic reprogramming of intracranial metastases is thought to be involved in therapy failure, but so far only little is known about patient-specific DNA-methylation differences between intra- and extracranial melanoma metastases. Hierarchical clustering of the methylomes of 24 patient-matched intra- and extracranial melanoma metastases pairs revealed that intra- and extracranial metastases of individual patients were more similar to each other than to metastases in the same tissue from other patients. Therefore, a personalized analysis of each metastases pair was done by a Hidden Markov Model to classify methylation levels of individual CpGs as decreased, unchanged or increased in the intra- compared to the extracranial metastasis. The predicted DNA-methylation alterations were highly patient-specific differing in the number and methylation states of altered CpGs. Nevertheless, four important general observations were made: (i) intracranial metastases of most patients mainly showed a reduction of DNA-methylation, (ii) cytokine signaling was most frequently affected by differential methylation in individual metastases pairs, but also MAPK, PI3K/Akt and ECM signaling were often altered, (iii) frequently affected genes were mainly involved in signaling, growth, adhesion or apoptosis, and (iv) an enrichment of functional terms related to channel and transporter activities supports previous findings for a brain-like phenotype. In addition, the derived set of 17 signaling pathway genes that distinguished intra- from extracranial metastases in more than 50% of patients included well-known oncogenes (e.g. PRKCA, DUSP6, BMP4) and several other genes known from neuronal disorders (e.g. EIF4B, SGK1, CACNG8). Moreover, associations of gene body methylation alterations with corresponding gene expression changes revealed that especially the three signaling pathway genes JAK3, MECOM, and TNXB differ strongly in their expression between patient-matched intra- and extracranial metastases. Our analysis contributes to an in-depth characterization of DNA-methylation differences between patient-matched intra- and extracranial melanoma metastases and may provide a basis for future experimental studies to identify targets for new therapeutic approaches.

Synonymous alterations of cancer-associated Trp53 CpG mutational hotspots cause fatal developmental jaw malocclusions but no tumors in knock-in mice

Intragenic CpG dinucleotides are tightly conserved in evolution yet are also vulnerable to methylation-dependent mutation, raising the question as to why these functionally critical sites have not been deselected by more stable coding sequences. We previously showed in cell lines that altered exonic CpG methylation can modify promoter start sites, and hence protein isoform expression, for the human TP53 tumor suppressor gene. Here we extend this work to the in vivo setting by testing whether synonymous germline modifications of exonic CpG sites affect murine development, fertility, longevity, or cancer incidence. We substituted the DNA-binding exons 5-8 of Trp53, the mouse ortholog of human TP53, with variant-CpG (either CpG-depleted or -enriched) sequences predicted to encode the normal p53 amino acid sequence; a control construct was also created in which all non-CpG sites were synonymously substituted. Homozygous Trp53-null mice were the only genotype to develop tumors. Mice with variant-CpG Trp53 sequences remained tumor-free, but were uniquely prone to dental anomalies causing jaw malocclusion (p < .0001). Since the latter phenotype also characterises murine Rett syndrome due to dysfunction of the trans-repressive MeCP2 methyl-CpG-binding protein, we hypothesise that CpG sites may exert non-coding phenotypic effects via pre-translational cis-interactions of 5-methylcytosine with methyl-binding proteins which regulate mRNA transcript initiation, expression or splicing, although direct effects on mRNA structure or translation are also possible.

Interactive DNA Methylation Array Analysis with ShinyÉPICo

Arrays provide a cost-effective platform for the analysis of human DNA methylation. ShinyÉPICo is an interactive, web-based, and graphical tool that allows the user to analyze Illumina DNA methylation arrays (450 k and EPIC), from the user's own computer or from a server. This tool covers the analysis entirely, from the raw data input to the final list of differentially methylated positions or regions. Here, we describe the steps of the analysis, the different parameters available, and useful information to understand and select the best options in each step.

Interactive regulation of DNA demethylase gene TET1 and m6A methyltransferase gene METTL3 in myoblast differentiation

DNA methylation (5mC) and mRNA N⁶-methyladenosine (m⁶A) play an essential role in gene transcriptional regulation. DNA methylation has been well established to be involved in skeletal muscle development. Interacting regulatory mechanisms between DNA methylation and mRNA m⁶A modification have been identified in a variety of biological processes. However, the effect of m⁶A on skeletal muscle differentiation and the underlying mechanisms are still unclear. It is also unknown whether there is an interaction between DNA methylation and mRNA m⁶A modification in skeletal myogenesis. In the present study, we used m⁶A-IP-qPCR, LC-MS/MS and dot blot assays to determine that the DNA demethylase gene, TET1, exhibited increased m⁶A levels and decreased mRNA expression during bovine skeletal myoblast differentiation. Dual-luciferase reporter assays and RIP experiments demonstrated that METTL3 suppressed TET1 expression by regulating TET1 mRNA stability in a m⁶A-YTHDF2-dependent manner. Furthermore, TET1 mediated DNA demethylation of itself, MYOD1 and MYOG, thereby stimulating their expression to promote myogenic differentiation. Ectopic expression of TET1 rescued the effect of METTL3 knockdown on reduced myotubes. In contrast, TET1 knockdown impaired the myogenic differentiation promoted by METTL3 overexpression. Moreover, ChIP experiments found that TET1 could bind and demethylate METTL3 DNA, which enhanced METTL3 expression. In addition, TET1 knockdown decreased m⁶A levels. ChIP assays also showed that TET1 knockdown contributed to the binding of H3K4me3 and H3K27me3 to METTL3 DNA. Our results revealed a negative feedback regulatory loop between TET1 and METTL3 in myoblast differentiation, which unveiled the interplay among DNA methylation, RNA methylation and histone methylation in skeletal myogenesis.

Identification of differentially methylated regions in rare diseases from a single-patient perspective

BACKGROUND

DNA methylation (5-mC) is being widely recognized as an alternative in the detection of sequence variants in the diagnosis of some rare neurodevelopmental and imprinting disorders. Identification of alterations in DNA methylation plays an important role in the diagnosis and understanding of the etiology of those disorders. Canonical pipelines for the detection of differentially methylated regions (DMRs) usually rely on inter-group (e.g., case versus control) comparisons. However, these tools might perform suboptimally in the context of rare diseases and multilocus imprinting disturbances due to small cohort sizes and inter-patient heterogeneity. Therefore, there is a need to provide a simple but statistically robust pipeline for scientists and clinicians to perform differential methylation analyses at the single patient level as well as to evaluate how parameter fine-tuning may affect differentially methylated region detection.

RESULT

We implemented an improved statistical method to detect differentially methylated regions in correlated datasets based on the Z-score and empirical Brown aggregation methods from a single-patient perspective. To accurately assess the predictive power of our method, we generated semi-simulated data using a public control population of 521 samples and investigated how the size of the control population, methylation difference, and region size affect DMR detection. In addition, we validated the detection of methylation events in patients suffering from rare multi-locus imprinting disturbance and evaluated how this method could complement existing tools in the context of clinical diagnosis.

CONCLUSION

In this study, we present a robust statistical method to perform differential methylation analysis at the single patient level and describe its optimal parameters to increase DMRs identification performance. Finally, we show its diagnostic utility when applied to rare disorders.

A seven-transmembrane protein-TM7SF3, resides in nuclear speckles and regulates alternative splicing

The seven-transmembrane superfamily member 3 protein (TM7SF3) is a p53-regulated homeostatic factor that attenuates cellular stress and the unfolded protein response. Here we show that TM7SF3 localizes to nuclear speckles; eukaryotic nuclear bodies enriched in splicing factors. This unexpected location for a trans -membranal protein enables formation of stable complexes between TM7SF3 and pre-mRNA splicing factors including DHX15, LARP7, HNRNPU, RBM14, and HNRNPK. Indeed, TM7SF3 regulates alternative splicing of >330 genes, mainly at the 3'end of introns by directly modulating the activity of splicing factors such as HNRNPK. These effects are observed both in cell lines and primary human pancreatic islets. Accordingly, silencing of TM7SF3 results in differential expression of 1465 genes (about 7% of the human genome); with 844 and 621 genes being up- or down-regulated, respectively. Our findings implicate TM7SF3, as a resident protein of nuclear speckles and suggest a role for seven-transmembrane proteins as regulators of alternative splicing.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Genome-wide alternation and effect of DNA methylation in the impairments of steroidogenesis and spermatogenesis after PM2.5 exposure

The effects of ambient fine particles on male reproductive health have raised widespread concern. The particular underlying mechanisms of the damage remain largely unclear and demand more research in new directions. Previous research has revealed that DNA methylation plays an important role in male reproductive development and is also vulnerable to environmental influences. However, there hasn't been enough investigation into the involvement of DNA methylation in PM_2.5-induced male reproductive toxicity. Here, we establish a real-time PM_2.5 exposure model and revealed that PM_2.5 exposure could lead to testicular dysfunction including spermatogenesis impairment and steroid hormone dysfunction. In particular, the decrease in the testicular global level of 5-methylcytosine (5mC) indicated a possible association of DNA methylation with testicular injury induced by PM_2.5 exposure. Further genome-wide methylation analysis revealed genomic hypomethylation of testicular DNA and identified more than 1000 differentially methylated regions in both CAP and UA versus FA, indicating that PM_2.5 exposure, even low-dose, could modulate the testicular methylome. Furthermore, integrated analysis of methylome and transcriptome identified some key methylated genes and networks, which may be involved in spermatogenesis and synthesis of steroid hormone. The testicular methylation levels of key genes especially Cyp11a1 and Pax8 raised, and their consequent reduced expression may impair the testosterone and sperm production process. Our research provides fundamental knowledge as well as novel insights into the possible involvement of DNA methylation in PM_2.5-induced male reproductive harm.

Adapting to climate with limited genetic diversity: Nucleotide, DNA methylation and microbiome variation among populations of the social spider Stegodyphus dumicola

Understanding the role of genetic and nongenetic variants in modulating phenotypes is central to our knowledge of adaptive responses to local conditions and environmental change, particularly in species with such low population genetic diversity that it is likely to limit their evolutionary potential. A first step towards uncovering the molecular mechanisms underlying population-specific responses to the environment is to carry out environmental association studies. We associated climatic variation with genetic, epigenetic and microbiome variation in populations of a social spider with extremely low standing genetic diversity. We identified genetic variants that are associated strongly with environmental variation, particularly with average temperature, a pattern consistent with local adaptation. Variation in DNA methylation in many genes was strongly correlated with a wide set of climate parameters, thereby revealing a different pattern of associations than that of genetic variants, which show strong correlations to a more restricted range of climate parameters. DNA methylation levels were largely independent of cis-genetic variation and of overall genetic population structure, suggesting that DNA methylation can work as an independent mechanism. Microbiome composition also correlated with environmental variation, but most strong associations were with precipitation-related climatic factors. Our results suggest a role for both genetic and nongenetic mechanisms in shaping phenotypic responses to local environments.

PARP1′s Involvement in RNA Polymerase II Elongation: Pausing and Releasing Regulation through the Integrator and Super Elongation Complex

RNA polymerase elongation along the gene body is tightly regulated to ensure proper transcription and alternative splicing events. Understanding the mechanism and factors critical in regulating the rate of RNA polymerase II elongation and processivity is clearly important. Recently we showed that PARP1, a well-known DNA repair protein, when bound to chromatin, regulates RNA polymerase II elongation. However, the mechanism by which it does so is not known. In the current study, we aimed to tease out how PARP1 regulates RNAPII elongation. We show, both in vivo and in vitro, that PARP1 binds directly to the Integrator subunit 3 (IntS3), a member of the elongation Integrator complex. The association between the two proteins is mediated via the C-terminal domain of PARP1 to the C-terminal domain of IntS3. Interestingly, the occupancy of IntS3 along two PARP1 target genes mimicked that of PARP1, suggesting a role in its recruitment/assembly of elongation factors. Indeed, the knockdown of PARP1 resulted in differential chromatin association and gene occupancy of IntS3 and other key elongation factors. Most of these PARP1-mediated effects were due to the physical presence of PARP1 rather than its PARylation activity. These studies argue that PARP1 controls the progressive RNAPII elongation complexes. In summary, we present a platform to begin to decipher PARP1′s role in recruiting/scaffolding elongation factors along the gene body regions during RNA polymerase II elongation and gene regulation.

Premature transcription termination at the expanded GAA repeats and aberrant alternative polyadenylation contributes to the Frataxin transcriptional deficit in Friedreich's ataxia

Frataxin deficiency in Friedreich's ataxia results from transcriptional downregulation of the FXN gene caused by expansion of the intronic trinucleotide guanine-adenine-adenine (GAA) repeats. We used multiple transcriptomic approaches to determine the molecular mechanism of transcription inhibition caused by long GAAs. We uncovered that transcription of FXN in patient cells is prematurely terminated upstream of the expanded repeats leading to the formation of a novel, truncated and stable RNA. This FXN early terminated transcript (FXN-ett) undergoes alternative, non-productive splicing and does not contribute to the synthesis of functional frataxin. The level the FXN-ett RNA directly correlates with the length of the longer of the two expanded GAA tracts. Targeting GAAs with antisense oligonucleotides or excision of the repeats eliminates the transcription impediment, diminishes expression of the aberrant FXN-ett, while increasing levels of FXN mRNA and frataxin. Non-productive transcription may represent a common phenomenon and attractive therapeutic target in diseases caused by repeat-mediated transcription aberrations.

Proteogenomic, Epigenetic, and Clinical Implications of Recurrent Aberrant Splice Variants in Clear Cell Renal Cell Carcinoma

BACKGROUND

Alternative mRNA splicing can be dysregulated in cancer, resulting in the generation of aberrant splice variants (SVs). Given the paucity of actionable genomic mutations in clear cell renal cell carcinoma (ccRCC), aberrant SVs may be an avenue to novel mechanisms of pathogenesis.

OBJECTIVE

To identify and characterize aberrant SVs enriched in ccRCC.

DESIGN, SETTING, AND PARTICIPANTS

Using RNA-seq data from the Cancer Cell Line Encyclopedia, we identified neojunctions uniquely expressed in ccRCC. Candidate SVs were then checked for expression across normal tissue in the Genotype-Tissue Expression Project and primary tumor tissue from The Cancer Genome Atlas (TCGA), Clinical Proteomic Tumor Analysis Consortium (CPTAC), and our institutional Total Cancer Care database.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Clinicopathologic, genomic, and survival data were available for all cohorts. Epigenetic data were available for the TCGA and CPTAC cohorts. Proteomic data were available for the CPTAC cohort. The association of aberrant SV expression with these variables was examined using the Kruskal-Wallis test, pairwise t test, Spearman correlation test, and Cox regression analysis.

RESULTS AND LIMITATIONS

Our pipeline identified 16 ccRCC-enriched SVs. EGFR, HPCAL1-SV and RNASET2-SV expression was negatively correlated with gene-specific CpG methylation. We derived a survival risk score based primarily on the expression of five SVs (RNASET2, FGD1, PDZD2, COBLL1, and PTPN14), which was consistent and applicable across multiple cohorts on multivariate analysis. The splicing factor RBM4, which modulates splicing of HIF-1α, exhibited significantly lower expression at the protein level in the high-risk group, as defined by our SV-based score.

CONCLUSIONS

We describe 16 aberrant SVs enriched in ccRCC, many of which are associated with disease biology and/or clinical outcomes. This study provides a novel strategy for identifying and characterizing disease-specific aberrant SVs.

PATIENT SUMMARY

We describe a method to identify disease targets and biomarkers using transcriptomic analysis beyond somatic mutations or gene expression. Kidney tumors express unique splice variants that may provide additional prognostic information following surgery.

Early Development of Cardiac Fibrosis in Young Old-Father Offspring

Cardiac aging is characterized by progressive fibrosis. Epidemiological studies have found that advanced paternal age is associated with an increased risk of heart failure in the next generation. This study is aimed at evaluating the effect of paternal age, in the young male rat progeny, on cardiac phenotype under circulatory stress conditions. Offspring rats were obtained by mating old males (24 months old) with young females (two months old) and by mating young males (two months old) with the same young females. Hypertension was induced in old father offspring (OFO) rats and young old father (YFO) offspring rats using L-NAME (N(ω)-nitro-L-arginine methyl ester). The OFO L-NAME rats showed a high blood pressure phenotype associated with substantial cardiac hypertrophy and an exacerbation of cardiac fibrosis compared to the YFO L-NAME rats. Histological analysis of heart tissue showed an expansion of the extracellular matrix, with fibroblasts displaying markers of epicardial origin (Tcf21, Tbx18, and Wt1) in the OFO group. Moreover, western blot and protein phosphorylation antibody array identified the TGF-β2 receptor pathway as preferentially activated in aged hearts as well as in OFO cardiac tissue treated with L-NAME. In addition, old father offspring rats (OFO+OFO L-NAME) had increased cardiac DNA methylation. In young hypertensive progeny, advanced paternal age at conception may be a risk factor for early progression towards cardiac fibrosis. An intergenerational transmission may be behind the paternal age-related cardiac remodeling in the young offspring.

Epigenetic Dysregulations in Arsenic-Induced Carcinogenesis

Arsenic is a crucial environmental metalloid whose high toxicity levels negatively impact human health. It poses significant health concerns to millions of people in developed and developing countries such as the USA, Canada, Bangladesh, India, China, and Mexico by enhancing sensitivity to various types of diseases, including cancers. However, how arsenic causes changes in gene expression that results in heinous conditions remains elusive. One of the proposed essential mechanisms that still has seen limited research with regard to causing disease upon arsenic exposure is the dysregulation of epigenetic components. In this review, we have extensively summarized current discoveries in arsenic-induced epigenetic modifications in carcinogenesis and angiogenesis. Importantly, we highlight the possible mechanisms underlying epigenetic reprogramming through arsenic exposure that cause changes in cell signaling and dysfunctions of different epigenetic elements.

Multi-omics analysis reveals RNA splicing alterations and their biological and clinical implications in lung adenocarcinoma

Alternative RNA splicing is one of the most important mechanisms of posttranscriptional gene regulation, which contributes to protein diversity in eukaryotes. It is well known that RNA splicing dysregulation is a critical mechanism in tumor pathogenesis and the rationale for the promising splice-switching therapeutics for cancer treatment. Although we have a comprehensive understanding of DNA mutations, abnormal gene expression profiles, epigenomics, and proteomics in lung adenocarcinoma (LUAD), little is known about its aberrant alternative splicing profiles. Here, based on the multi-omics data generated from over 1000 samples, we systematically studied the RNA splicing alterations in LUAD and revealed their biological and clinical implications. We identified 3688 aberrant alternative splicing events (AASEs) in LUAD, most of which were alternative promoter and exon skip. The specific regulatory roles of RNA binding proteins, somatic mutations, and DNA methylations on AASEs were comprehensively interrogated. We dissected the functional implications of AASEs and concluded that AASEs mainly affected biological processes related to tumor proliferation and metastasis. We also found that one subtype of LUAD with a particular AASEs pattern was immunogenic and had a better prognosis and response rate to immunotherapy. These findings revealed novel events related to tumorigenesis and tumor immune microenvironment and laid the foundation for the development of splice-switching therapies for LUAD.

Association of DNA methylation with steroidogenic enzymes in Cushing's adenoma

DNA methylation and demethylation regulate the transcription of genes. DNA methylation-associated gene expression of adrenal steroidogenic enzymes may regulate cortisol production in cortisol-producing adenoma (CPA). We aimed to determine the DNA methylation levels of all genes encoding steroidogenic enzymes involved in CPA. Additionally, the aims were to clarify the DNA methylation-associated gene expression and evaluate the difference of CPA genotype from others using DNA methylation data. Twenty-five adrenal CPA and six nonfunctioning adrenocortical adenoma (NFA) samples were analyzed. RNA sequencing and DNA methylation array were performed. The methylation levels at 118 methylation sites of the genes were investigated, and their methylation and mRNA levels were subsequently integrated. Among all the steroidogenic enzyme genes studied, CYP17A1 gene was mainly found to be hypomethylated in CPA compared to that in NFA, and the Benjamini-Hochberg procedure demonstrated that methylation levels at two sites in the CYP17A1 gene body were statistically significant. PRKACA mutant CPAs predominantly exhibited hypomethylation of CYP17A1 gene compared with the GNAS mutant CPAs. Inverse associations between CYP17A1 methylation in three regions of the gene body and its mRNA levels were observed in the NFAs and CPAs. In applying clustering analysis using CYP17A1 methylation and mRNA levels, CPAs with PRKACA mutation were differentiated from NFAs and CPAs with a GNAS mutation. We demonstrated that CPAs exhibited hypomethylation of the CYP17A1 gene body in CPA, especially in the PRKACA mutant CPAs. Methylation of CYP17A1 gene may influence its transcription levels.

miR778 mediates gene expression, histone modification, and DNA methylation during cyst nematode parasitism

Despite the known critical regulatory functions of microRNAs, histone modifications, and DNA methylation in reprograming plant epigenomes in response to pathogen infection, the molecular mechanisms underlying the tight coordination of these components remain poorly understood. Here, we show how Arabidopsis (Arabidopsis thaliana) miR778 coordinately modulates the root transcriptome, histone methylation, and DNA methylation via post-transcriptional regulation of the H3K9 methyltransferases SU(var)3-9 homolog 5 (SUVH5) and SUVH6 upon infection by the beet cyst nematode Heterodera schachtii. miR778 post-transcriptionally silences SUVH5 and SUVH6 upon nematode infection. Manipulation of the expression of miR778 and its two target genes significantly altered plant susceptibility to H. schachtii. RNA-seq analysis revealed a key role of SUVH5 and SUVH6 in reprograming the transcriptome of Arabidopsis roots upon H. schachtii infection. In addition, chromatin immunoprecipitation (ChIP)-seq analysis established SUVH5 and SUVH6 as the main enzymes mediating H3K9me2 deposition in Arabidopsis roots in response to nematode infection. ChIP-seq analysis also showed that these methyltransferases possess distinct DNA binding preferences in that they are targeting transposable elements under noninfected conditions and protein-coding genes in infected plants. Further analyses indicated that H3K9me2 deposition directed by SUVH5 and SUVH6 contributes to gene expression changes both in roots and in nematode feeding sites and preferentially associates with CG DNA methylation. Together, our results uncovered multi-layered epigenetic regulatory mechanisms coordinated by miR778 during Arabidopsis-H. schachtii interactions.

Maternal attachment insecurity, maltreatment history, and depressive symptoms are associated with broad DNA methylation signatures in infants

The early environment, including maternal characteristics, provides many cues to young organisms that shape their long-term physical and mental health. Identifying the earliest molecular events that precede observable developmental outcomes could help identify children in need of support prior to the onset of physical and mental health difficulties. In this study, we examined whether mothers' attachment insecurity, maltreatment history, and depressive symptoms were associated with alterations in DNA methylation patterns in their infants, and whether these correlates in the infant epigenome were associated with socioemotional and behavioral functioning in toddlerhood. We recruited 156 women oversampled for histories of depression, who completed psychiatric interviews and depression screening during pregnancy, then provided follow-up behavioral data on their children at 18 months. Buccal cell DNA was obtained from 32 of their infants for a large-scale analysis of methylation patterns across 5 × 106 individual CpG dinucleotides, using clustering-based significance criteria to control for multiple comparisons. We found that tens of thousands of individual infant CpGs were alternatively methylated in association with maternal attachment insecurity, maltreatment in childhood, and antenatal and postpartum depressive symptoms, including genes implicated in developmental patterning, cell-cell communication, hormonal regulation, immune function/inflammatory response, and neurotransmission. Density of DNA methylation at selected genes from the result set was also significantly associated with toddler socioemotional and behavioral problems. This is the first report to identify novel regions of the human infant genome at which DNA methylation patterns are associated longitudinally both with maternal characteristics and with offspring socioemotional and behavioral problems in toddlerhood.

FBXL7 Body Hypomethylation Is Frequent in Tumors from the Digestive and Respiratory Tracts and Is Associated with Risk-Factor Exposure

Squamous cell carcinoma is the main histological tumor type in the upper aerodigestive tract (UADT), including the esophagus (ESCC) and the head and neck sites, as well as the oral cavity (OCSCC), larynx (LSCC) and oropharynx (OPSCC). These tumors are induced by alcohol and tobacco exposure, with the exception of a subgroup of OPSCC linked to human papillomavirus (HPV) infection. Few genes are frequently mutated in UADT tumors, pointing to other molecular mechanisms being involved during carcinogenesis. The F-box and leucine-rich repeat protein 7 (FBXL7) is a potential tumor-suppressing gene, one that is frequently hypermethylated in pancreatic cancer and where the encoded protein promotes the degradation of AURKA, BIRC5 and c-SRC. Thus, the aim of this study was to evaluate the methylation and expression profile of FBXL7 in the UADT and the gene’s association with the clinical, etiological and pathological characteristics of patients, as well as the expression of its degradation targets. Here we show that the FBXL7 gene’s body is hypomethylated in the UADT, independently of histology, but not in virus-associated tumors. FBXL7 body methylation and gene expression levels were correlated in the ESCC, LSCC, OCSCC and OPSCC. Immunohistochemistry analysis showed that FBXL7 protein levels are not correlated with the levels of its degradation targets, AURKA and BIRC5, in the UADT. The high discriminatory potential of FBXL7 body hypomethylation between non-tumor and tumor tissues makes it a promising biomarker.

Reversal of Epigenetic Peroxisome Proliferator-Activated Receptor-γ Suppression by Diacerein Alleviates Oxidative Stress and Osteoarthritis in Mice

Aims: The pathogenesis of osteoarthritis (OA) is characterized by oxidative stress (OS) and sustained inflammation that are substantially associated with epigenetic DNA methylation alterations of osteogenic gene expression. Diacerein as an anthraquinone anti-OA drug exhibits multiple chondroprotective properties, but less clarified pharmacological actions. Since anthraquinone contain an epigenetic modulating property, in this study we investigate whether the anti-OA functions of diacerein involve DNA methylation modulation and antioxidant signaling. Results: The OA mice incurred by destabilization of medial meniscus exhibited marked suppression of peroxisome proliferator-activated receptor-gamma (PPARγ), a chondroprotective transcription factor with anti-inflammation and OS-balancing properties, aberrant upregulations of DNA methyltransferase (DNMT)1/3a, and PPARγ promoter hypermethylation in knee joint cartilage. Diacerein treatment mitigated the cartilage damage and significantly inhibited the DNMT1/3a upregulation, the PPARγ promoter hypermethylation, and the PPARγ loss, and it effectively corrected the adverse expression of antioxidant enzymes and inflammatory cytokines. In cultured chondrocytes, diacerein reduced the interleukin-1β-induced PPARγ suppression and the abnormal expression of its downstream antioxidant enzymes in a gain of DNMT and PPARγ inhibition-sensitive manner, and in PPARγ knockout mice, the anti-OA effects of diacerein were significantly reduced. Innovation: Our work reveals a novel anti-OA pharmacological property of diacerein and identifies the aberrant DNMT elevation and the resultant PPARγ suppression as an important epigenetic pathway that mediates diacerein's anti-OA activities. Conclusion: DNA methylation aberration and the resultant PPARγ suppression contribute significantly to epigenetic OA pathogenesis, and targeting PPARγ suppression via DNA demethylation is an important component of diacerein's anti-OA functions. Antioxid. Redox Signal. 37, 40-53.

Alternative Splicing and Its Roles in Plant Metabolism

Plant metabolism, including primary metabolism such as tricarboxylic acid cycle, glycolysis, shikimate and amino acid pathways as well as specialized metabolism such as biosynthesis of phenolics, alkaloids and saponins, contributes to plant survival, growth, development and interactions with the environment. To this end, these metabolic processes are tightly and finely regulated transcriptionally, post-transcriptionally, translationally and post-translationally in response to different growth and developmental stages as well as the constantly changing environment. In this review, we summarize and describe the current knowledge of the regulation of plant metabolism by alternative splicing, a post-transcriptional regulatory mechanism that generates multiple protein isoforms from a single gene by using alternative splice sites during splicing. Numerous genes in plant metabolism have been shown to be alternatively spliced under different developmental stages and stress conditions. In particular, alternative splicing serves as a regulatory mechanism to fine-tune plant metabolism by altering biochemical activities, interaction and subcellular localization of proteins encoded by splice isoforms of various genes.

The potential consequences of bidirectional promoter methylation on GLA and HNRNPH2 expression in Fabry disease phenotypes in a family of patients carrying a GLA deletion variant

Fabry disease (FD) is a rare inherited disease characterized by a wide range of symptoms attributed to GLA mutations resulting in defective α-galactosidase A (α-Gal A) and accumulation of glycosphingolipids. The GLA locus is paired in a divergent manner with the heterogeneous nuclear ribonucleoprotein HNRNPH2 locus mapped in the RPL36A-HNRNPH2 readthrough locus. As a follow-up to our recent finding of the co-regulation of GLA and HNRNPH2 via a bidirectional promoter (BDP) in normal kidney and skin cells, the potential accumulative influence of BDP methylation and GLA mutation on the severity of FD in patients from the same family, two males and two females carrying a GLA deletion mutation, c.1033_1034delTC (p.Ser345Argfs) was addressed in the present study. The molecular analyses of the FD patients compared with the control revealed that the expression of GLA was significantly low (P<0.05), and HNRNPH2 showed a tendency of low expression (P=0.1) when BDP methylation was elevated in FD patients, compared with low BDP methylation and high GLA expression (P<0.05), and a high trend of HNRNPH2 expression in normal individuals. The accumulative effects of the mutation and BDP methylation with the severity of the disease were observed in three patients. One male FD patient, a member of the FD family diagnosed with progressive loss of kidney function, hypertension, and eventually a stroke, and the lowest level of α-Gal A enzyme activity showed the highest BDP DNA methylation level. It is concluded that the DNA methylation of GLA-HNRNPH2 BDP may serve a role in diagnosing and treating FD.

Relevance and Regulation of Alternative Splicing in Plant Heat Stress Response: Current Understanding and Future Directions

Alternative splicing (AS) is a major mechanism for gene expression in eukaryotes, increasing proteome diversity but also regulating transcriptome abundance. High temperatures have a strong impact on the splicing profile of many genes and therefore AS is considered as an integral part of heat stress response. While many studies have established a detailed description of the diversity of the RNAome under heat stress in different plant species and stress regimes, little is known on the underlying mechanisms that control this temperature-sensitive process. AS is mainly regulated by the activity of splicing regulators. Changes in the abundance of these proteins through transcription and AS, post-translational modifications and interactions with exonic and intronic cis-elements and core elements of the spliceosomes modulate the outcome of pre-mRNA splicing. As a major part of pre-mRNAs are spliced co-transcriptionally, the chromatin environment along with the RNA polymerase II elongation play a major role in the regulation of pre-mRNA splicing under heat stress conditions. Despite its importance, our understanding on the regulation of heat stress sensitive AS in plants is scarce. In this review, we summarize the current status of knowledge on the regulation of AS in plants under heat stress conditions. We discuss possible implications of different pathways based on results from non-plant systems to provide a perspective for researchers who aim to elucidate the molecular basis of AS under high temperatures.

Paradigm shifts in animal epigenetics: Research on non-model species leads to new insights into dependencies, functions and inheritance of DNA methylation

Recent investigations with non-model species and whole-genome approaches have challenged several paradigms in animal epigenetics. They revealed that epigenetic variation in populations is not the mere consequence of genetic variation, but is a semi-independent or independent source of phenotypic variation, depending on mode of reproduction. DNA methylation is not positively correlated with genome size and phylogenetic position as earlier believed, but has evolved differently between and within higher taxa. Epigenetic marks are usually not completely erased in the zygote and germ cells as generalized from mouse, but often persist and can be transgenerationally inherited, making them evolutionarily relevant. Gene body methylation and promoter methylation are similar in vertebrates and invertebrates with well methylated genomes but transposon silencing through methylation is variable. The new data also suggest that animals use epigenetic mechanisms to cope with rapid environmental changes and to adapt to new environments. The main benefiters are asexual populations, invaders, sessile taxa and long-lived species.

LARGE-SCALE INTEGRATION OF DNA METHYLATION AND GENE EXPRESSION ARRAY PLATFORMS IDENTIFIES BOTH cis AND trans RELATIONSHIPS

Although epigenome-wide association studies (EWAS) have been successful in identifying DNA methylation (DNAm) patterns associated with disease states, any further characterization of etiologic mechanisms underlying disease remains elusive. This knowledge gap does not originate from a lack of DNAm–trait associations, but rather stems from study design issues that affect the interpretability of EWAS results. Despite known limitations in predicting the function of a particular CpG site, most EWAS maintain the broad assumption that altered DNAm results in a concomitant change of transcription at the most proximal gene. This study integrated DNAm and gene expression (GE) measurements in two cohorts, the Adolescent and Young Adult Twin Study (AYATS) and the Pregnancy, Race, Environment, Genes (PREG) study, to improve the understanding of epigenomic regulatory mechanisms. CpG sites associated with GE in cis were enriched in areas of transcription factor binding and areas of intermediate-to-low CpG density. CpG sites associated with trans GE were also enriched in areas of known regulatory significance, including enhancer regions. These results highlight issues with restricting DNAm-transcript annotations to small genomic intervals and question the validity of assuming a cis DNAm–GE pathway. Based on these findings, the interpretation of EWAS results is limited in studies without multi-omic support and further research should identify genomic regions in which GE-associated DNAm is overrepresented. An in-depth characterization of GE-associated CpG sites could improve predictions of the downstream functional impact of altered DNAm and inform best practices for interpreting DNAm–trait associations generated by EWAS.

LDICDL: LncRNA-Disease Association Identification Based on Collaborative Deep Learning

It has been proved that long noncoding RNA (lncRNA) plays critical roles in many human diseases. Therefore, inferring associations between lncRNAs and diseases can contribute to disease diagnosis, prognosis and treatment. To overcome the limitation of traditional experimental methods such as expensive and time-consuming, several computational methods have been proposed to predict lncRNA-disease associations by fusing different biological data. However, the prediction performance of lncRNA-disease associations identification needs to be improved. In this study, we propose a computational model (named LDICDL) to identify lncRNA-disease associations based on collaborative deep learning. It uses an automatic encoder to denoise multiple lncRNA feature information and multiple disease feature information, respectively. Then, the matrix decomposition algorithm is employed to predict the potential lncRNA-disease associations. In addition, to overcome the limitation of matrix decomposition, the hybrid model is developed to predict associations between new lncRNA (or disease) and diseases (or lncRNA). The ten-fold cross validation and de novo test are applied to evaluate the performance of method. The experimental results show LDICDL outperforms than other state-of-the-art methods in prediction performance.

Gene-body methylation in plants: mechanisms, functions and important implications for understanding evolutionary processes

Gene body methylation (gbM) is an epigenetic mark where gene exons are methylated in the CG context only, as opposed to CHG and CHH contexts (where H stands for A, C, or T). CG methylation is transmitted transgenerationally in plants, opening the possibility that gbM may be shaped by adaptation. This presupposes, however, that gbM has a function that affects phenotype, which has been a topic of debate in the literature. Here, we review our current knowledge of gbM in plants. We start by presenting the well-elucidated mechanisms of plant gbM establishment and maintenance. We then review more controversial topics: the evolution of gbM and the potential selective pressures that act on it. Finally, we discuss the potential functions of gbM that may affect organismal phenotypes: gene expression stabilization and upregulation, inhibition of aberrant transcription (reverse and internal), prevention of aberrant intron retention, and protection against TE insertions. To bolster the review of these topics, we include novel analyses to assess the effect of gbM on transcripts. Overall, a growing body of literature finds that gbM correlates with levels and patterns of gene expression. It is not clear, however, if this is a causal relationship. Altogether, functional work suggests that the effects of gbM, if any, must be relatively small, but there is nonetheless evidence that it is shaped by natural selection. We conclude by discussing the potential adaptive character of gbM and its implications for an updated view of the mechanisms of adaptation in plants.

Genome-Wide DNA Methylation and Its Effect on Gene Expression During Subclinical Mastitis in Water Buffalo

Subclinical mastitis (SCM) in buffalo is one of the most challenging paradoxes for the dairy sector with very significant milk production losses and poses an imminent danger to milch animal’s milk-producing ability. We present here the genome-wide methylation specific to SCM in water buffalo and its consequential effect on the gene expression landscape for the first time. Whole-genome DNA methylation profiles from peripheral blood lymphocytes and gene expression profiles from milk somatic cells of healthy and SCM cases were catalogued from the MeDIP-Seq and RNA-Seq data. The average methylation in healthy buffaloes was found to be higher than that in the SCM-infected buffaloes. DNA methylation was abundant in the intergenic region followed by the intronic region in both healthy control and SCM groups. A total of 3,950 differentially methylated regions (DMRs) were identified and annotated to 370 differentially methylated genes (DMGs), most of which were enriched in the promoter region. Several important pathways were activated due to hypomethylation and belonged to the Staphylococcus aureus infection, Th17 cell differentiation, and antigen processing and presentation pathways along with others of defense responses. DNA methylome was compared with transcriptome to understand the regulatory role of DNA methylation on gene expression specific to SCM in buffaloes. A total of 4,778 significant differentially expressed genes (DEGs) were extracted in response to SCM, out of which 67 DMGs were also found to be differentially expressed, suggesting that during SCM, DNA methylation could be one of the epigenetic regulatory mechanisms of gene expression. Genes like CSF2RB, LOC102408349, C3 and PZP like, and CPAMD8 were found to be downregulated in our study, which are known to be involved in the immune response to SCM. Association of DNA methylation with transposable elements, miRNAs, and lncRNAs was also studied. The present study reports a buffalo SCM web resource (BSCM2TDb) available at http://webtom.cabgrid.res.in/BSCM2TDb that catalogues all the mastitis-related information of the analyses results of this study in a single place. This will be of immense use to buffalo researchers to understand the host–pathogen interaction involving SCM, which is required in endeavors of mastitis control and management.

The Epigenetics of Psychosis: A Structured Review with Representative Loci

The evidence for an environmental component in chronic psychotic disorders is strong and research on the epigenetic manifestations of these environmental impacts has commenced in earnest. In reviewing this research, the focus is on three genes as models for differential methylation, MCHR1, AKT1 and TDO2, each of which have been investigated for genetic association with psychotic disorders. Environmental factors associated with psychotic disorders, and which interact with these model genes, are explored in depth. The location of transcription factor motifs relative to key methylation sites is evaluated for predicted gene expression results, and for other sites, evidence is presented for methylation directing alternative splicing. Experimental results from key studies show differential methylation: for MCHR1, in psychosis cases versus controls; for AKT1, as a pre-existing methylation pattern influencing brain activation following acute administration of a psychosis-eliciting environmental stimulus; and for TDO2, in a pattern associated with a developmental factor of risk for psychosis, in all cases the predicted expression impact being highly dependent on location. Methylation induced by smoking, a confounding variable, exhibits an intriguing pattern for all three genes. Finally, how differential methylation meshes with Darwinian principles is examined, in particular as it relates to the “flexible stem” theory of evolution.

Overactivation of IL6 cis‑signaling in leukocytes is an inflammatory hallmark of deep vein thrombosis

Inflammation is a protective response of the body to various injuries, which is strictly regulated by a variety of factors, including immune cells and soluble mediators. However, dysfunction of this defensive mechanism often results in inflammation‑driven diseases, such as deep vein thrombosis (DVT). The complex relationship between inflammatory cell activity and DVT has not been fully elucidated. The present study aimed to investigate the role of interleukin‑6 (IL6) signaling transduction in DVT. To this aim, the expression levels of transmembrane isoforms of the IL6 receptor (IL6R) and the glycoprotein 130 responsible for the IL6 cis‑signaling were evaluated in the peripheral blood mononuclear cells of patients with DVT and of healthy controls. The results indicated that leukocytes from patients with DVT exhibited overexpression of both IL6R and gp130 membrane isoforms and that these were strongly associated with the occurrence of DVT. Overall, the present findings indicated that IL6 cis‑signaling may have a direct involvement in the leukocyte activation in DVT and may serve as a predictive biomarker of DVT development.

Transcriptional Alterations in X-Linked Dystonia–Parkinsonism Caused by the SVA Retrotransposon

X-linked dystonia–parkinsonism (XDP) is a severe neurodegenerative disorder that manifests as adult-onset dystonia combined with parkinsonism. A SINE-VNTR-Alu (SVA) retrotransposon inserted in an intron of the TAF1 gene reduces its expression and alters splicing in XDP patient-derived cells. As a consequence, increased levels of the TAF1 intron retention transcript TAF1-32i can be found in XDP cells as compared to healthy controls. Here, we investigate the sequence of the deep intronic region included in this transcript and show that it is also present in cells from healthy individuals, albeit in lower amounts than in XDP cells, and that it undergoes degradation by nonsense-mediated mRNA decay. Furthermore, we investigate epigenetic marks (e.g., DNA methylation and histone modifications) present in this intronic region and the spanning sequence. Finally, we show that the SVA evinces regulatory potential, as demonstrated by its ability to repress the TAF1 promoter in vitro. Our results enable a better understanding of the disease mechanisms underlying XDP and transcriptional alterations caused by SVA retrotransposons.

Methionine- and Choline-Deficient Diet Identifies an Essential Role for DNA Methylation in Plasmacytoid Dendritic Cell Biology

Diet plays an important role in lifestyle disorders associated with the disturbed immune system. During the study of methionine- and choline-deficient diet-induced nonalcoholic fatty liver disease, we observed a specific decrease in the plasmacytoid dendritic cell (pDC) fraction from murine spleens. While delineating the role for individual components, we identified that l-methionine supplementation correlates with representation of the pDC fraction. S-adenosylmethionine (SAM) is a key methyl donor, and we demonstrate that supplementation of methionine-deficient medium with SAM but not homocysteine reverses the defect in pDC development. l-Methionine has been implicated in maintenance of methylation status in the cell. Based on our observed effect of SAM and zebularine on DC subset development, we sought to clarify the role of DNA methylation in pDC biology. Whole-genome bisulfite sequencing analysis from the splenic DC subsets identified that pDCs display differentially hypermethylated regions in comparison with classical DC (cDC) subsets, whereas cDC1 and cDC2 exhibited comparable methylated regions, serving as a control in our study. We validated differentially methylated regions in the sorted pDC, CD8α⁺ cDC1, and CD4⁺ cDC2 subsets from spleens as well as FL-BMDC cultures. Upon analysis of genes linked with differentially methylated regions, we identified that differential DNA methylation is associated with the MAPK pathway such that its inhibition guides DC development toward the pDC subtype. Overall, our study identifies an important role for methionine in pDC biology.

Gene architecture directs splicing outcome in separate nuclear spatial regions

How the splicing machinery defines exons or introns as the spliced unit has remained a puzzle for 30 years. Here, we demonstrate that peripheral and central regions of the nucleus harbor genes with two distinct exon-intron GC content architectures that differ in the splicing outcome. Genes with low GC content exons, flanked by long introns with lower GC content, are localized in the periphery, and the exons are defined as the spliced unit. Alternative splicing of these genes results in exon skipping. In contrast, the nuclear center contains genes with a high GC content in the exons and short flanking introns. Most splicing of these genes occurs via intron definition, and aberrant splicing leads to intron retention. We demonstrate that the nuclear periphery and center generate different environments for the regulation of alternative splicing and that two sets of splicing factors form discrete regulatory subnetworks for the two gene architectures. Our study connects 3D genome organization and splicing, thus demonstrating that exon and intron definition modes of splicing occur in different nuclear regions.

A Genome-Wide Investigation of Effects of Aberrant DNA Methylation on the Usage of Alternative Promoters in Hepatocellular Carcinoma

BACKGROUND

The alternative usage of promoters provides a way to regulate gene expression, has a significant influence on the transcriptome, and contributes to the cellular transformation of cancer. However, the function of alternative promoters (APs) in hepatocellular carcinoma (HCC) has not been systematically studied yet. In addition, the potential mechanism of regulation to the usage of APs remains unclear. DNA methylation, one of the most aberrant epigenetic modifications in cancers, is known to regulate transcriptional activity. Whether DNA methylation regulates the usage of APs needs to be explored. Here, we aim to investigate the effects of DNA methylation on usage of APs in HCC.

METHODS

Promoter activities were calculated based on RNA-seq data. Functional enrichment analysis was implemented to conduct GO terms. Correlation tests were used to detect the correlation between promoter activity and methylation status. The LASSO regression model was used to generate a diagnostic model. Kaplan-Meier analysis was used to compare the overall survival between high and low methylation groups. RNA-seq and whole-genome bisulfite sequencing (WGBS) in HCC samples were performed to validate the correlation of promoter activity and methylation.

RESULTS

We identified 855 APs in total, which could be well used to distinguish cancer from normal samples. The correlation of promoter activity and DNA methylation in APs was observed, and the APs with negative correlation were defined as methylation-regulated APs (mrAPs). Six mrAPs were identified to generate a diagnostic model with good performance (AUC = 0.97). Notably, the majority of mrAPs had CpG sites that could be used to predict clinical outcomes by methylation status. Finally, we verified 85.6% of promoter activity variation and 92.3% of methylation changes in our paired RNA-seq and WGBS samples, respectively. The negative correlation between promoter activity and methylation status was further confirmed in our HCC samples.

CONCLUSION

The aberrant methylation status plays a critical role in the precision usage of APs in HCC, which sheds light on the mechanism of cancer development and provides a new insight into cancer screening and treatment.

Whole-Genome DNA Methylation Dynamics of Sheep Preimplantation Embryo Investigated by Single-Cell DNA Methylome Sequencing

The early stages of mammalian embryonic development involve the participation and cooperation of numerous complex processes, including nutritional, genetic, and epigenetic mechanisms. However, in embryos cultured in vitro, a developmental block occurs that affects embryo development and the efficiency of culture. Although the block period is reported to involve the transcriptional repression of maternal genes and transcriptional activation of zygotic genes, how epigenetic factors regulate developmental block is still unclear. In this study, we systematically analyzed whole-genome methylation levels during five stages of sheep oocyte and preimplantation embryo development using single-cell level whole genome bisulphite sequencing (SC-WGBS) technology. Then, we examined several million CpG sites in individual cells at each evaluated developmental stage to identify the methylation changes that take place during the development of sheep preimplantation embryos. Our results showed that two strong waves of methylation changes occurred, namely, demethylation at the 8-cell to 16-cell stage and methylation at the 16-cell to 32-cell stage. Analysis of DNA methylation patterns in different functional regions revealed a stable hypermethylation status in 3'UTRs and gene bodies; however, significant differences were observed in intergenic and promoter regions at different developmental stages. Changes in methylation at different stages of preimplantation embryo development were also compared to investigate the molecular mechanisms involved in sheep embryo development at the methylation level. In conclusion, we report a detailed analysis of the DNA methylation dynamics during the development of sheep preimplantation embryos. Our results provide an explanation for the complex regulatory mechanisms underlying the embryo developmental block based on changes in DNA methylation levels.

Using Machine Learning to Predict Obesity Based on Genome-Wide and Epigenome-Wide Gene-Gene and Gene-Diet Interactions

Obesity is associated with many chronic diseases that impair healthy aging and is governed by genetic, epigenetic, and environmental factors and their complex interactions. This study aimed to develop a model that predicts an individual's risk of obesity by better characterizing these complex relations and interactions focusing on dietary factors. For this purpose, we conducted a combined genome-wide and epigenome-wide scan for body mass index (BMI) and up to three-way interactions among 402,793 single nucleotide polymorphisms (SNPs), 415,202 DNA methylation sites (DMSs), and 397 dietary and lifestyle factors using the generalized multifactor dimensionality reduction (GMDR) method. The training set consisted of 1,573 participants in exam 8 of the Framingham Offspring Study (FOS) cohort. After identifying genetic, epigenetic, and dietary factors that passed statistical significance, we applied machine learning (ML) algorithms to predict participants' obesity status in the test set, taken as a subset of independent samples (n = 394) from the same cohort. The quality and accuracy of prediction models were evaluated using the area under the receiver operating characteristic curve (ROC-AUC). GMDR identified 213 SNPs, 530 DMSs, and 49 dietary and lifestyle factors as significant predictors of obesity. Comparing several ML algorithms, we found that the stochastic gradient boosting model provided the best prediction accuracy for obesity with an overall accuracy of 70%, with ROC-AUC of 0.72 in test set samples. Top predictors of the best-fit model were 21 SNPs, 230 DMSs in genes such as CPT1A, ABCG1, SLC7A11, RNF145, and SREBF1, and 26 dietary factors, including processed meat, diet soda, French fries, high-fat dairy, artificial sweeteners, alcohol intake, and specific nutrients and food components, such as calcium and flavonols. In conclusion, we developed an integrated approach with ML to predict obesity using omics and dietary data. This extends our knowledge of the drivers of obesity, which can inform precision nutrition strategies for the prevention and treatment of obesity. Clinical Trial Registration: [www.ClinicalTrials.gov], the Framingham Heart Study (FHS), [NCT00005121].

The Neuroepigenetic Landscape of Vertebrate and Invertebrate Models of Neurodegenerative Diseases

Vertebrate and invertebrate models of neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis, have been paramount to our understanding of the pathophysiology of these conditions; however, the brain epigenetic landscape is less well established in these disease models. DNA methylation, histone modifications, and microRNAs are among commonly studied mechanisms of epigenetic regulation. Genome-wide studies and candidate studies of specific methylation marks, histone marks, and microRNAs have demonstrated the dysregulation of these mechanisms in models of neurodegenerative diseases; however, the studies to date are scarce and inconclusive and the implications of many of these changes are still not fully understood. In this review, we summarize epigenetic changes reported to date in the brain of vertebrate and invertebrate models used to study neurodegenerative diseases, specifically diseases affecting the aging population. We also discuss caveats of epigenetic research so far and the use of disease models to understand neurodegenerative diseases, with the aim of improving the use of model organisms in this context in future studies.

The Identification of Prognostic and Metastatic Alternative Splicing in Skin Cutaneous Melanoma

Skin cutaneous melanoma (SKCM) is a type of highly invasive cancer originated from melanocytes. It is reported that aberrant alternative splicing (AS) plays an important role in the neoplasia and metastasis of many types of cancer. Therefore, we investigated whether ASEs of pre-RNA have such an influence on the prognosis of SKCM and the related mechanism of ASEs in SKCM. The RNA-seq data and ASEs data for SKCM patients were obtained from the TCGA and TCGASpliceSeq database. The univariate Cox regression revealed 1265 overall survival-related splicing events (OS-SEs). Screened by Lasso regression, 4 OS-SEs were identified and used to construct an effective prediction model (AUC: .904), whose risk score was proved to be an independent prognostic factor. Furthermore, Kruskal-Wallis test and Mann-Whitney-Wilcoxon test showed that an aberrant splicing type of aminoacyl tRNA synthetase complex-interacting multifunctional protein 2 (AIMP2) regulated by CDC-like kinase 1 (CLK1) was associated with the metastasis and stage of SKCM. Besides, the overlapped signal pathway for AIMP2 was galactose metabolism identified by the co-expression analysis. External database validation also confirmed that AIMP2, CLK1, and the galactose metabolism were associated with the metastasis and stage of SKCM patients. ChIP-seq and ATAC-seq methods further confirmed the transcription regulation of CLK1, AIMP2, and other key genes, whose cellular expression was detected by Single Cell Sequencing. In conclusion, we proposed that CLK1-regulated AIMP2-78704-ES might play a critical role in the tumorigenesis and metastasis of SKCM via galactose metabolism. Besides, we established an effective model with MTMR14-63114-ES, URI1-48867-ES, BATF2-16724-AP, and MED22-88025-AP to predict the metastasis and prognosis of SKCM patients.

Deciphering the synergistic and redundant roles of CG and non-CG DNA methylation in plant development and transposable element silencing

DNA methylation plays key roles in transposable element (TE) silencing and gene expression regulation. DNA methylation occurs at CG, CHG and CHH sequence contexts in plants. However, the synergistic and redundant roles of CG and non-CG methylation are poorly understood. By introducing CRISPR/Cas9-induced met1 mutation into the ddcc (drm1 drm2 cmt2 cmt3) mutant, we attempted to knock out all five DNA methyltransferases in Arabidopsis and then investigate the synergistic and redundant roles of CG and non-CG DNA methylation. We found that the homozygous ddcc met1 quintuple mutants are embryonically lethal, although met1 and ddcc mutants only display some developmental abnormalities. Unexpectedly, the ddcc met1 quintuple mutations only reduce transmission through the male gametophytes. The ddcc met1^+/- mutants show apparent size divergence, which is not associated with difference in DNA methylation patterns, but associated with the difference in the levels of DNA damage. Finally, we show that a group of TEs are specifically activated in the ddcc met1^+/- mutants. This work reveals that CG and non-CG DNA methylation synergistically and redundantly regulate plant reproductive development, vegetative development and TE silencing in Arabidopsis. Our findings provide insights into the roles of DNA methylation in plant development.

OUP accepted manuscript

Alternative splicing (AS) and transcription elongation are vital biological processes, and their dysregulation causes multiple diseases, including tumors. However, the coregulatory mechanism of AS and transcription elongation in tumors remains unclear. This study demonstrates a novel AS pattern of tight junction protein 1 (ZO1) regulated by the RNA polymerase II elongation rate in colorectal cancer (CRC). Glioma tumor suppressor candidate region gene 1 (GLTSCR1) decreases the transcription elongation rate of ZO1 to provide a time window for binding of the splicing factor HuR to the specific motif in intron 22 of ZO1 and spliceosome recognition of the weak 3′ and 5′ splice sites in exon 23 to promote exon 23 inclusion. Since exon 23 inclusion in ZO1 suppresses migration and invasion of CRC cells, our findings suggest a novel potential therapeutic target for CRC.

Epigenetic modifications of tumor necrosis factor-alpha in joint cartilage tissue from osteoarthritis patients - CONSORT

BACKGROUND

Osteoarthritis (OA) remains one of the most common osteopathy for centuries, which can be attributed to multiple risk factors including mechanical and biochemical ones. More and more studies verified that inflammatory cytokines play important roles in the progression of OA, such as tumor necrosis factor-alpha (TNF-α). In this study, we aimed to investigate the relationship between epigenetic manifestations of TNF-? and the pathogenesis of OA.

METHODS

Totally, 37 OA patients' cartilage was collected through the knee joint and 13 samples of articular cartilage as healthy control was collected through traumatic amputation. Real-time PCR, Western blot and ELISA analysis were performed to observe the expression of target genes and proteins in collected samples.

RESULTS

Compared with the healthy control group, TNF-? was over-expressing in cartilage which was collected from OA patients. DNA hypomethylation, histone hyperacetylation and histone methylation were observed in the TNF-? promoter in OA compared with normal patients, and we also studied series of enzymes associated with epigenetics. The results showed that by increasing DNA methylation and decreasing histone acetylation in the TNF-? promoter, and TNF-? over-expression in OA cartilage was suppressed, histone methylation has no significant correlation with OA.

CONCLUSION

In conclusion, the changes of epigenetic status regulate TNF-α expression in the cells, which are pivotal to the OA disease process. These results may give us a better understanding of OA and may provide new therapeutic options.

The Chinese pine genome and methylome unveil key features of conifer evolution

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

Epigenetic drug development for autoimmune and inflammatory diseases

ES-62 is the major secreted protein of the parasitic filarial nematode, Acanthocheilonema viteae. The molecule exists as a large tetramer (MW, ~240kD), which possesses immunomodulatory properties by virtue of multiple phosphorylcholine (PC) moieties attached to N-type glycans. By suppressing inflammatory immune responses, ES-62 can prevent disease development in certain mouse models of allergic and autoimmune conditions, including joint pathology in collagen-induced arthritis (CIA), a model of rheumatoid arthritis (RA). Such protection is associated with functional suppression of "pathogenic" hyper-responsive synovial fibroblasts (SFs), which exhibit an aggressive inflammatory and bone-damaging phenotype induced by their epigenetic rewiring in response to the inflammatory microenvironment of the arthritic joint. Critically, exposure to ES-62 in vivo induces a stably-imprinted CIA-SF phenotype that exhibits functional responses more typical of healthy, Naïve-SFs. Consistent with this, ES-62 "rewiring" of SFs away from the hyper-responsive phenotype is associated with suppression of ERK activation, STAT3 activation and miR-155 upregulation, signals widely associated with SF pathogenesis. Surprisingly however, DNA methylome analysis of Naïve-, CIA- and ES-62-CIA-SF cohorts reveals that rather than simply preventing pathogenic rewiring of SFs, ES-62 induces further changes in DNA methylation under the inflammatory conditions pertaining in the inflamed joint, including targeting genes associated with ciliogenesis, to programme a novel "resolving" CIA-SF phenotype. In addition to introducing a previously unsuspected aspect of ES-62's mechanism of action, such unique behaviour signposts the potential for developing DNA methylation signatures predictive of pathogenesis and its resolution and hence, candidate mechanisms by which novel therapeutic interventions could prevent SFs from perpetuating joint inflammation and destruction in RA. Pertinent to these translational aspects of ES-62-behavior, small molecule analogues (SMAs) based on ES-62's active PC-moieties mimic the rewiring of SFs as well as the protection against joint disease in CIA afforded by the parasitic worm product.

The pattern of DNA methylation alteration, and its association with the changes of gene expression and alternative splicing during phosphate starvation in tomato

DNA methylation is changed and associates with gene expression alterations in plant response to phosphate starvation (Pi-), a common stress that impacts plant growth and productivity. However, in the horticultural model species Solanum lycopersicum (tomato), the dynamics of DNA methylation and its relationship with changes in gene transcription and alternative splicing (AS) under Pi- are unknown. Here, we performed integrative methylome and transcriptome analyses of tomato seedlings under Pi-deficient and -sufficient conditions. We found Pi- caused a slight increase in the overall methylation level, with millions of differentially methylated cytosines (DmCs) and a few hundred differentially methylated regions (DMRs). We also identified thousands of differentially expressed (DE) and differential AS (DAS) genes induced by Pi-, and found that DmCs were more abundant in non-expressed genes than in DE or DAS genes. Moreover, DNA methylation alterations weakly correlated with transcription changes but not with DAS events, and hyper-CHH-DMRs overlapping with transposable elements (TEs) were enriched in a subset of Pi starvation response (PSR) genes. We propose that changes in DNA methylation may be associated with the differential expression of some PSR genes, but that most of these changes probably control the expression of nearby TEs, rather than directly affecting the transcription or AS of PSR genes. Besides, the pattern of methylation changes upon Pi- may largely be shaped by TE distributions. Together, our study provides comprehensive insights into the association of DNA methylation with gene transcription and AS under Pi- in tomato and may contribute to unveiling novel roles of epigenetic mechanisms in plant stress response.

Arsenic-induced epigenetic changes in cancer development

Arsenic is a ubiquitous metalloid whose high levels of toxicity pose major health concerns to millions of people worldwide by increasing susceptibility to various cancers and non-cancer illnesses. Since arsenic is not a mutagen, the mechanism by which it causes changes in gene expression and disease pathogenesis is not clear. One possible mechanism is through generation of reactive oxygen species. Another equally important mechanism still very much in its infancy is epigenetic dysregulation. In this review, we discuss recent discoveries underlying arsenic-induced epigenetic changes in cancer development. Importantly, we highlight the proposed mechanisms targeted by arsenic to drive oncogenic gene expression.