Epigenetics of Skeletal Muscle-Associated Genes in the ASB, LRRC, TMEM, and OSBPL Gene Families

Histopathologic image-based deep learning classifier for predicting platinum-based treatment responses in high-grade serous ovarian cancer

Platinum-based chemotherapy is the cornerstone treatment for female high-grade serous ovarian carcinoma (HGSOC), but choosing an appropriate treatment for patients hinges on their responsiveness to it. Currently, no available biomarkers can promptly predict responses to platinum-based treatment. Therefore, we developed the Pathologic Risk Classifier for HGSOC (PathoRiCH), a histopathologic image-based classifier. PathoRiCH was trained on an in-house cohort (n = 394) and validated on two independent external cohorts (n = 284 and n = 136). The PathoRiCH-predicted favorable and poor response groups show significantly different platinum-free intervals in all three cohorts. Combining PathoRiCH with molecular biomarkers provides an even more powerful tool for the risk stratification of patients. The decisions of PathoRiCH are explained through visualization and a transcriptomic analysis, which bolster the reliability of our model's decisions. PathoRiCH exhibits better predictive performance than current molecular biomarkers. PathoRiCH will provide a solid foundation for developing an innovative tool to transform the current diagnostic pipeline for HGSOC.

Late-Stage Skeletal Muscle Transcriptome in Duchenne muscular dystrophy shows a BMP4-Induced Molecular Signature

Duchenne muscular dystrophy (DMD) is a fatal X-linked recessive disease due to loss-of-function mutations in the DYSTROPHIN gene. DMD-related skeletal muscle wasting is typified by an aberrant immune response involving upregulation of TGFβ family of cytokines. We previously demonstrated that bone morphogenetic protein 4 (BMP4) is increased in DMD and BMP4 stimulation induces a 20-fold upregulation of Smad8 transcription. However, the role of BMP4 in severely affected DMD skeletal muscle is unknown. We hypothesized that transcriptomic signatures in severely affected human DMD skeletal muscle are driven by BMP4 signaling. Transcriptomes from skeletal muscle biopsies of late-stage DMD vs. non-DMD controls and C2C12 muscle cells with or without BMP4 stimulation were generated by RNA-Seq and analyzed for single transcript differential expression as well as by Ingenuity Pathway Analysis and weighted gene co-expression network analyses. A total of 2,328 and 5,291 transcripts in the human muscle and C2C12 muscle cells, respectively, were differentially expressed. We identified an overlapping molecular signature of 1,027 genes dysregulated in DMD muscle that were induced in BMP4-stimulated C2C12 muscle cells. Highly upregulated DMD transcripts that overlapped with BMP4-stimulated C2C12 muscle cells included ADAMTS3, HCAR2, SERPING1, SMAD8 , and UNC13C. The DMD transcriptome was characterized by dysregulation of pathways involving immune function, extracellular matrix remodeling, and metabolic/mitochondrial function. In summary, we define a late-stage DMD skeletal muscle transcriptome that substantially overlaps with the BMP4-induced molecular signature in C2C12 muscle cells. This supports BMP4 as a disease-driving regulator of transcriptomic changes in late-stage DMD skeletal muscle and expands our understanding of the evolution of dystrophic signaling pathways and their associated gene networks that could be explored for therapeutic development.

Restoration of epigenetic impairment in the skeletal muscle and chronic inflammation resolution as a therapeutic approach in sarcopenia

Sarcopenia is an age-associated loss of skeletal muscle mass, strength, and function, accompanied by severe adverse health outcomes, such as falls and fractures, functional decline, high health costs, and mortality. Hence, its prevention and treatment have become increasingly urgent. However, despite the wide prevalence and extensive research on sarcopenia, no FDA-approved disease-modifying drugs exist. This is probably due to a poor understanding of the mechanisms underlying its pathophysiology. Recent evidence demonstrate that sarcopenia development is characterized by two key elements: (i) epigenetic dysregulation of multiple molecular pathways associated with sarcopenia pathogenesis, such as protein remodeling, insulin resistance, mitochondria impairments, and (ii) the creation of a systemic, chronic, low-grade inflammation (SCLGI). In this review, we focus on the epigenetic regulators that have been implicated in skeletal muscle deterioration, their individual roles, and possible crosstalk. We also discuss epidrugs, which are the pharmaceuticals with the potential to restore the epigenetic mechanisms deregulated in sarcopenia. In addition, we discuss the mechanisms underlying failed SCLGI resolution in sarcopenia and the potential application of pro-resolving molecules, comprising specialized pro-resolving mediators (SPMs) and their stable mimetics and receptor agonists. These compounds, as well as epidrugs, reveal beneficial effects in preclinical studies related to sarcopenia. Based on these encouraging observations, we propose the combination of epidrugs with SCLI-resolving agents as a new therapeutic approach for sarcopenia that can effectively attenuate of its manifestations.

Endothelial gene regulatory elements associated with cardiopharyngeal lineage differentiation

Endothelial cells (EC) differentiate from multiple sources, including the cardiopharyngeal mesoderm, which gives rise also to cardiac and branchiomeric muscles. The enhancers activated during endothelial differentiation within the cardiopharyngeal mesoderm are not completely known. Here, we use a cardiogenic mesoderm differentiation model that activates an endothelial transcription program to identify endothelial regulatory elements activated in early cardiogenic mesoderm. Integrating chromatin remodeling and gene expression data with available single-cell RNA-seq data from mouse embryos, we identify 101 putative regulatory elements of EC genes. We then apply a machine-learning strategy, trained on validated enhancers, to predict enhancers. Using this computational assay, we determine that 50% of these sequences are likely enhancers, some of which are already reported. We also identify a smaller set of regulatory elements of well-known EC genes and validate them using genetic and epigenetic perturbation. Finally, we integrate multiple data sources and computational tools to search for transcriptional factor binding motifs. In conclusion, we show EC regulatory sequences with a high likelihood to be enhancers, and we validate a subset of them using computational and cell culture models. Motif analyses show that the core EC transcription factors GATA/ETS/FOS is a likely driver of EC regulation in cardiopharyngeal mesoderm.

Epigenetics of Genes Preferentially Expressed in Dissimilar Cell Populations: Myoblasts and Cerebellum

While studying myoblast methylomes and transcriptomes, we found that CDH15 had a remarkable preference for expression in both myoblasts and cerebellum. To understand how widespread such a relationship was and its epigenetic and biological correlates, we systematically looked for genes with similar transcription profiles and analyzed their DNA methylation and chromatin state and accessibility profiles in many different cell populations. Twenty genes were expressed preferentially in myoblasts and cerebellum (Myob/Cbl genes). Some shared DNA hypo- or hypermethylated regions in myoblasts and cerebellum. Particularly striking was ZNF556, whose promoter is hypomethylated in expressing cells but highly methylated in the many cell populations that do not express the gene. In reporter gene assays, we demonstrated that its promoter's activity is methylation sensitive. The atypical epigenetics of ZNF556 may have originated from its promoter's hypomethylation and selective activation in sperm progenitors and oocytes. Five of the Myob/Cbl genes (KCNJ12, ST8SIA5, ZIC1, VAX2, and EN2) have much higher RNA levels in cerebellum than in myoblasts and displayed myoblast-specific hypermethylation upstream and/or downstream of their promoters that may downmodulate expression. Differential DNA methylation was associated with alternative promoter usage for Myob/Cbl genes MCF2L, DOK7, CNPY1, and ANK1. Myob/Cbl genes PAX3, LBX1, ZNF556, ZIC1, EN2, and VAX2 encode sequence-specific transcription factors, which likely help drive the myoblast and cerebellum specificity of other Myob/Cbl genes. This study extends our understanding of epigenetic/transcription associations related to differentiation and may help elucidate relationships between epigenetic signatures and muscular dystrophies or cerebellar-linked neuropathologies.

Implications of tRNA abundance on translation elongation across bovine tissues

Introduction: Translation is a crucial stage of gene expression. It may also act as an additional layer of regulation that plays an important role in gene expression and function. Highly expressed genes are believed to be codon-biased to support increased protein production, in which quickly translated codons correspond to highly abundant tRNAs. Synonymous SNPs, considered to be silent due to the degeneracy of the genetic code, may shift protein abundance and function through alterations in translational efficiency and suboptimal pairing to lowly abundant tRNAs. Methods: Here, we applied Quantitative Mature tRNA sequencing (QuantM-tRNAseq) and ribosome profiling across bovine tissues in order to investigate the relationship between tRNA expression and slowed translation. Results: Moreover, we have identified genes modulated at transcriptional and/or translational levels underlying tissue-specific biological processes. We have also successfully defined pausing sites that depict the regulatory information encoded within the open reading frame of transcripts, which could be related to translation rate and facilitate proper protein folding. This work offers an atlas of distinctive pausing sites across three bovine tissues, which provides an opportunity to predict codon optimality and understand tissue-specific mechanisms of regulating protein synthesis.

Establishment of a 6-signature risk model associated with cellular senescence for predicting the prognosis of breast cancer

This study focused on screening novel markers associated with cellular senescence for predicting the prognosis of breast cancer. The RNA-seq expression profile of BRCA and clinical data were obtained from TCGA. The pam algorithm was used to cluster patients based on senescence-related genes. The weighted gene co-expression network analysis was used to identify co-expressed genes, and LASSO-Cox analysis was performed to build a risk prognosis model. The performance of the model was also evaluated. We additionally explored the role of senescence in cancer development and possible regulatory mechanism. The patients were clustered into 2 subtypes. A total of 5259 genes significantly related to senescence were identified by weighted gene co-expression network analysis. LASSO-Cox finally established a 6-signature risk model (ADAMTS8, DCAF12L2, PCDHA10, PGK1, SLC16A2, and TMEM233) that exhibited favorable and stable performance in our training, validation, and whole BRCA datasets. Furthermore, the superiority of our model was also observed after comparing it to other published models. The 6-signature was proved to be an independent risk factor for prognosis. In addition, mechanism prediction implied the activation of glycometabolism processes such as glycolysis and TCA cycle under the condition of senescence. Glycometabolism pathways were further found to negatively correlate with the infiltration level of CD8 T-cells and natural killer cells but positively correlate with M2 macrophage infiltration and expressions of tissue degeneration biomarkers, which suggested the deficit immune surveillance and risk of tumor migration. The constructed 6-gene model based on cellular senescence could be an effective indicator for predicting the prognosis of BRCA.

asb5a/asb5b Double Knockout Affects Zebrafish Cardiac Contractile Function

Ankyrin repeat and suppression-of-cytokine-signaling box (Asb) proteins, a subset of ubiquitin ligase E3, include Asb5 with six ankyrin-repeat domains. Zebrafish harbor two asb5 gene isoforms, asb5a and asb5b. Currently, the effects of asb5 gene inactivation on zebrafish embryonic development and heart function are unknown. Using CRISPR/Cas9, we generated asb5a-knockout zebrafish, revealing no abnormal phenotypes at 48 h post-fertilization (hpf). In situ hybridization showed similar asb5a and asb5b expression patterns, indicating the functional redundancy of these isoforms. Morpholino interference was used to target asb5b in wild-type and asb5a-knockout zebrafish. Knocking down asb5b in the wild-type had no phenotypic impact, but simultaneous asb5b knockdown in asb5a-knockout homozygotes led to severe pericardial cavity enlargement and atrial dilation. RNA-seq and cluster analyses identified significantly enriched cardiac muscle contraction genes in the double-knockout at 48 hpf. Moreover, semi-automatic heartbeat analysis demonstrated significant changes in various heart function indicators. STRING database/Cytoscape analyses confirmed that 11 cardiac-contraction-related hub genes exhibited disrupted expression, with three modules containing these genes potentially regulating cardiac contractile function through calcium ion channels. This study reveals functional redundancy in asb5a and asb5b, with simultaneous knockout significantly impacting zebrafish early heart development and contraction, providing key insights into asb5's mechanism.

Exploring tumor-normal cross-talk with TranNet: Role of the environment in tumor progression

There is a growing awareness that tumor-adjacent normal tissues used as control samples in cancer studies do not represent fully healthy tissues. Instead, they are intermediates between healthy tissues and tumors. The factors that contribute to the deviation of such control samples from healthy state include exposure to the tumor-promoting factors, tumor-related immune response, and other aspects of tumor microenvironment. Characterizing the relation between gene expression of tumor-adjacent control samples and tumors is fundamental for understanding roles of microenvironment in tumor initiation and progression, as well as for identification of diagnostic and prognostic biomarkers for cancers. To address the demand, we developed and validated TranNet, a computational approach that utilizes gene expression in matched control and tumor samples to study the relation between their gene expression profiles. TranNet infers a sparse weighted bipartite graph from gene expression profiles of matched control samples to tumors. The results allow us to identify predictors (potential regulators) of this transition. To our knowledge, TranNet is the first computational method to infer such dependencies. We applied TranNet to the data of several cancer types and their matched control samples from The Cancer Genome Atlas (TCGA). Many predictors identified by TranNet are genes associated with regulation by the tumor microenvironment as they are enriched in G-protein coupled receptor signaling, cell-to-cell communication, immune processes, and cell adhesion. Correspondingly, targets of inferred predictors are enriched in pathways related to tissue remodelling (including the epithelial-mesenchymal Transition (EMT)), immune response, and cell proliferation. This implies that the predictors are markers and potential stromal facilitators of tumor progression. Our results provide new insights into the relationships between tumor adjacent control sample, tumor and the tumor environment. Moreover, the set of predictors identified by TranNet will provide a valuable resource for future investigations.

TMEM200A is a potential prognostic biomarker and correlated with immune infiltrates in gastric cancer

Background

Gastric cancer (GC) is one of the most common malignant tumors in the digestive system. Several transmembrane (TMEM) proteins are defined as tumor suppressors or oncogenes. However, the role and underlying mechanism of TMEM200A in GC remain unclear.

Methods

We analyzed the expression of TMEM200A in GC. Furthermore, the influence of TMEM200A on survival of GC patients was evaluated. The correlations between the clinical information and TMEM200A expression were analyzed using chi-square test and logistic regression. Relevant prognostic factors were identified performing univariate and multivariate analysis. Gene set enrichment analysis (GSEA) was performed based on the TCGA dataset. Finally, we explore the relationship between TMEM200A expression and cancer immune infiltrates using CIBERSORT.

Results

TMEM200A was up-regulated in GC tissues than that in adjacent non-tumor tissues based on TCGA database. Meta-analysis and RT-qPCR validated the difference in TMEM200A expression. Kaplan-Meier curves suggested the increased TMEM200A had a poor prognosis in GC patients. The chi-square test and logistic regression analyses showed that the TMEM200A expression correlates significantly with T stage. Multivariate analysis showed that TMEM200A expression might be an important independent predictor of poor overall survival in GC patients. GSEA identified five immune-related signaling pathways and five tumor-related signaling pathways significantly enriched in the high TMEM200A expression phenotype pathway. Finally, we found CD8+ T cells is apparently decreased in high TMEM200A expression group. Conversely, eosinophils is increased in high expression group compared with low expression group.

Conclusion

TMEM200A is a potential prognostic biomarker and correlated with immune infiltrates in GC.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

Mining of chicken muscle growth genes and the function of important candidate gene RPL3L in muscle development

The birth weight of chickens does not significantly affect the weight at slaughter, while the different growth rate after birth was one of the important reasons for the difference in slaughter weight. Also, the increase in chickens' postnatal skeletal muscle weight is the main cause of the slaughter weight gain, but which genes are involved in this biological process is still unclear. In this study, by integrating four transcriptome datasets containing chicken muscles at different developmental times or different chicken tissues in public databases, a total of nine candidate genes that may be related to postnatal muscle development in chickens were obtained, including RPL3L, FBP2, ASB4, ASB15, CKMT2, PGAM1, YIPF7, PFKM, and LDHA. One of these candidate genes is RPL3L, whose 42 bp insertion/deletion (indel) mutation significantly correlated with multiple carcass traits in the F2 resource population from Xinghua chickens crossing with White Recessive Rock (WRR) chickens, including live weight, carcass weight, half eviscerated weight, eviscerated weight, breast meat weight, wing weight, leg muscle shear force, and breast muscle shear force. Also, there was a very significant difference between different genotypes of the RPL3L 42 bp indel mutation in these trains. Further experiments showed that RPL3L was highly expressed in chicken skeletal muscle, and its overexpression could promote the proliferation and inhibit the differentiation of chicken myoblasts by regulating ASB4 and ASB15 expression. Our findings demonstrated that the RPL3L 42 bp indel may be one of the molecular markers of chicken weight-related traits.

Integration of RNA-seq and ATAC-seq identifies muscle-regulated hub genes in cattle

As the main product of livestock, muscle itself plays an irreplaceable role in maintaining animal body movement and regulating metabolism. Therefore, it is of great significance to explore its growth, development and regeneration to improve the meat yield and quality of livestock. In this study, we attempted to use RNA-seq and ATAC-seq techniques to identify differentially expressed genes (DEGs) specifically expressed in bovine skeletal muscle as potential candidates for studying the regulatory mechanisms of muscle development. Microarray data from 8 tissue samples were selected from the GEO database for analysis. First, we obtained gene modules related to each tissue through WGCNA analysis. Through Gene Ontology (GO) functional annotation, the module of lightyellow (ME_lightyellow) was closely related to muscle development, and 213 hub genes were screened as follow-up research targets. Further, the difference analysis showed that, except for PREB, all other candidate hub genes were up-regulated (muscle group vs. other-group). ATAC-seq analysis showed that muscle-specific accessible chromatin regions were mainly located in promoter of genes related to muscle structure development (GO:0061061), muscle cell development (GO:0055001) and muscle system process (GO:0003012), which were involved in cAMP, CGMP-PKG, MAPK, and other signaling pathways. Next, we integrated the results of RNA-seq and ATAC-seq analysis, and 54 of the 212 candidate hub genes were identified as key regulatory genes in skeletal muscle development. Finally, through motif analysis, 22 of the 54 key genes were found to be potential target genes of transcription factor MEF2C. Including CAPN3, ACTN2, MB, MYOM3, SRL, CKM, ALPK3, MAP3K20, UBE2G1, NEURL2, CAND2, DOT1L, HRC, MAMSTR, FSD2, LRRC2, LSMEM1, SLC29A2, FHL3, KLHL41, ATXN7L2, and PDRG1. This provides a potential reference for studying the molecular mechanism of skeletal muscle development in mammals.

Squalene through Its Post-Squalene Metabolites Is a Modulator of Hepatic Transcriptome in Rabbits

Squalene is a natural bioactive triterpene and an important intermediate in the biosynthesis of sterols. To assess the effect of this compound on the hepatic transcriptome, RNA-sequencing was carried out in two groups of male New Zealand rabbits fed either a diet enriched with 1% sunflower oil or the same diet with 0.5% squalene for 4 weeks. Hepatic lipids, lipid droplet area, squalene, and sterols were also monitored. The Squalene administration downregulated 9 transcripts and upregulated 13 transcripts. The gene ontology of transcripts fitted into the following main categories: transporter of proteins and sterols, lipid metabolism, lipogenesis, anti-inflammatory and anti-cancer properties. When the results were confirmed by RT-qPCR, rabbits receiving squalene displayed significant hepatic expression changes of LOC100344884 (PNPLA3), GCK, TFCP2L1, ASCL1, ACSS2, OST4, FAM91A1, MYH6, LRRC39, LOC108176846, GLT1D1 and TREH. A squalene-enriched diet increased hepatic levels of squalene, lanosterol, dihydrolanosterol, lathosterol, zymostenol and desmosterol. Strong correlations were found among specific sterols and some squalene-changed transcripts. Incubation of the murine AML12 hepatic cell line in the presence of lanosterol, dihydrolanosterol, zymostenol and desmosterol reproduced the observed changes in the expressions of Acss2, Fam91a1 and Pnpla3. In conclusion, these findings indicate that the squalene and post-squalene metabolites play important roles in hepatic transcriptional changes required to protect the liver against malfunction.

Risks and rewards of big-data in epigenomics research: an interview with Melanie Ehrlich

Melanie Ehrlich, PhD, is a professor in the Tulane Cancer Center, the Tulane Center for Medical Bioinformatics and Genomics and the Hayward Human Genetics Program at Tulane Medical School, New Orleans, LA. She obtained her PhD in molecular biology in 1971 from the State University of New York at Stony Brook and completed postdoctoral research at Albert Einstein College of Medicine in 1972. She has been working on various aspects of epigenetics, starting with DNA methylation, since 1973. Her group made many first findings about DNA methylation (see below). For example, in 1982 and 1983, in collaboration with Charles Gehrke at the University of Missouri, she was the first to report tissue-specific and cancer-specific differences in overall DNA methylation in humans. In 1985, Xian-Yang Zhang and Richard Wang in her lab discovered a class of human DNA sequences specifically hypomethylated in sperm. In 1998, her group was the first to describe extensive losses of DNA methylation in pericentromeric and centromeric DNA repeats in human cancer. Her lab's many publications on the prevalence of both DNA hypermethylation and hypomethylation in the same cancers brought needed balance to our understanding of the epigenetics of cancer and to its clinical implications [1]. Besides working on cancer epigenetics, her research group has helped elucidate cytogenetic and gene expression abnormalities in the immunodeficiency, centromeric and facial anomalies (ICF) syndrome, a rare recessive disease often caused by mutations in DNMT3B. Her group also studied the epigenetics and transcriptomics of facioscapulohumeral muscular dystrophy (FSHD), whose disease locus is a tandem 3.3-kb repeat at subtelomeric 4q (that happens to be hypomethylated in ICF DNA [2]). Her study of FSHD has taken her in the direction of muscle (skeletal muscle, heart and aorta) epigenetics [3-6]. Recently, she has led research that applies epigenetics much more rigorously than usual to the evaluation of genetic variants from genome-wide association studies (GWAS) of osteoporosis and obesity. In continued collaboration with Sriharsa Pradhan at New England Biolabs and Michelle Lacey at Tulane University, she has compared 5-hydroxymethylcytosine and 5-methylcytosine clustering in various human tissues [7] and is studying myoblast methylomes that they generated by a new high-resolution enzymatic technique (enzymatic methyl-seq).

Epigenetics of Mitochondria-Associated Genes in Striated Muscle

Striated muscle has especially large energy demands. We identified 97 genes preferentially expressed in skeletal muscle and heart, but not in aorta, and found significant enrichment for mitochondrial associations among them. We compared the epigenomic and transcriptomic profiles of the 27 genes associated with striated muscle and mitochondria. Many showed strong correlations between their tissue-specific transcription levels, and their tissue-specific promoter, enhancer, or open chromatin as well as their DNA hypomethylation. Their striated muscle-specific enhancer chromatin was inside, upstream, or downstream of the gene, throughout much of the gene as a super-enhancer (CKMT2, SLC25A4, and ACO2), or even overlapping a neighboring gene (COX6A2, COX7A1, and COQ10A). Surprisingly, the 3' end of the 1.38 Mb PRKN (PARK2) gene (involved in mitophagy and linked to juvenile Parkinson's disease) displayed skeletal muscle/myoblast-specific enhancer chromatin, a myoblast-specific antisense RNA, as well as brain-specific enhancer chromatin. We also found novel tissue-specific RNAs in brain and embryonic stem cells within PPARGC1A (PGC-1α), which encodes a master transcriptional coregulator for mitochondrial formation and metabolism. The tissue specificity of this gene's four alternative promoters, including a muscle-associated promoter, correlated with nearby enhancer chromatin and open chromatin. Our in-depth epigenetic examination of these genes revealed previously undescribed tissue-specific enhancer chromatin, intragenic promoters, regions of DNA hypomethylation, and intragenic noncoding RNAs that give new insights into transcription control for this medically important set of genes.

Methylome of skeletal muscle tissue in patients with hypertension and diabetes undergoing cardiopulmonary bypass

Background: Epigenomic changes occurring during surgery have been neglected in research; diabetes and hypertension can affect the epigenome but little is known about the epigenetics of skeletal muscle (SKM). Methods: DNA methylation was profiled via Illumina MethylationEPIC arrays in SKM samples obtained at the beginning and end of heart surgery with cardiopulmonary bypass. Results: Methylation in patients with hypertension and diabetes was significantly different, more so for uncontrolled diabetes; hypertension alone produced minimal effect. The affected pathways involved IL-1, IL-12, IL-18, TNF-α, IFN-γ, VEGF, NF-κB and Wnt signaling, apoptosis and DNA damage response. Significant changes occurred during surgery and included loci in the Hippo-YAP/TAZ pathway. Conclusion: Cardiopulmonary bypass surgery affects the SKM methylome, and the combination of hypertension and diabetes induces changes in the SKM epigenome in contrast to hypertension alone.

Altered expression of lncRNAs in autism spectrum disorder

Long non-coding RNAs (lncRNAs) have been recognized as an important epigenetic factor in the evolution of neuropsychiatric conditions. We have selected five lncRNAs (DISC2, PRKAR2A-AS1, LOC105375675, LRRC2-AS1, and LOC101928237) to measure their expression in blood samples of children with autism spectrum disorder (ASD) versus children with normal development. Expressions of DISC2, PRKAR2A-AS1 and LOC101928237 have been enhanced in ASD cases compared with healthy children (Posterior Beta = 2.508, P value<0.0001; Posterior Beta = 2.793, P value = 0.014 and Posterior Beta = 1.646, P value <0.0001, respectively). On the other hand, expression of LRRC2-AS1 has been lower in ASD patients compared with controls (Posterior Beta = -3.781, P value<0.0001). Remarkably, expression of DISC2 and PRKAR2A-AS1 have been lower in girls compared with boys (Posterior Beta = -0.982, P value<0.0001 and Posterior Beta = -0.135, P value<0.0001, respectively). In addition, expression of DISC2 has been lower in ASD cases aged more than 6 compared with those aged less than 6 years (Posterior Beta = -0.876, P value = 0.003). DISC2, LOC101928237, LRRC2-AS1, and PRKAR2A-AS1 had the area under curve (AUC) values of 0.76, 0.90, 0.92, and 0.79 in distinguishing between ASD and healthy children. Expression levels of none of DISC2, LOC101928237, LOC105375675, LRRC2-AS1, and PRKAR2A-AS1 were correlated with age of ASD cases or healthy controls. A significant correlation was detected between expressions of DISC2 and PRKAR2A-AS1. There were inverse correlations between the following pairs of lncRNAs: DISC2/LRRC2-AS1, DISC2/LOC101928237, LRRC2-AS1/PRKAR2A-AS1, LOC101928237/LRRC2-AS1, and LOC101928237 /LOC105375675. We conclude that DISC2, LOC101928237, LRRC2-AS1, and PRKAR2A-AS1 might be used as potential markers for this condition.

Epigenetics of Muscle- and Brain-Specific Expression of KLHL Family Genes

KLHL and the related KBTBD genes encode components of the Cullin-E3 ubiquitin ligase complex and typically target tissue-specific proteins for degradation, thereby affecting differentiation, homeostasis, metabolism, cell signaling, and the oxidative stress response. Despite their importance in cell function and disease (especially, KLHL40, KLHL41, KBTBD13, KEAP1, and ENC1), previous studies of epigenetic factors that affect transcription were predominantly limited to promoter DNA methylation. Using diverse tissue and cell culture whole-genome profiles, we examined 17 KLHL or KBTBD genes preferentially expressed in skeletal muscle or brain to identify tissue-specific enhancer and promoter chromatin, open chromatin (DNaseI hypersensitivity), and DNA hypomethylation. Sixteen of the 17 genes displayed muscle- or brain-specific enhancer chromatin in their gene bodies, and most exhibited specific intergenic enhancer chromatin as well. Seven genes were embedded in super-enhancers (particularly strong, tissue-specific clusters of enhancers). The enhancer chromatin regions typically displayed foci of DNA hypomethylation at peaks of open chromatin. In addition, we found evidence for an intragenic enhancer in one gene upregulating expression of its neighboring gene, specifically for KLHL40/HHATL and KLHL38/FBXO32 gene pairs. Many KLHL/KBTBD genes had tissue-specific promoter chromatin at their 5' ends, but surprisingly, two (KBTBD11 and KLHL31) had constitutively unmethylated promoter chromatin in their 3' exons that overlaps a retrotransposed KLHL gene. Our findings demonstrate the importance of expanding epigenetic analyses beyond the 5' ends of genes in studies of normal and abnormal gene regulation.