Expression and regulatory roles of lncRNAs in G-CIMP-low vs G-CIMP-high Glioma: an in-silico analysis

BACKGROUND: Clinically relevant glioma subtypes, such as the glioma-CpG island methylator phenotype (G-CIMP), have been defined by epigenetics. In this study, the role of long non-coding RNAs in association with the poor-prognosis G-CMIP-low phenotype and the good-prognosis G-CMIP-high phenotype was investigated. Functional associations of lncRNAs with mRNAs and miRNAs were examined to hypothesize influencing factors of the aggressive phenotype. METHODS: RNA-seq data on 250 samples from TCGA\u27s Pan-Glioma study, quantified for lncRNA and mRNAs (GENCODE v28), were analyzed for differential expression between G-CIMP-low and G-CIMP-high phenotypes. Functional interpretation of the differential lncRNAs was performed by Ingenuity Pathway Analysis. Spearman rank order correlation estimates between lncRNA, miRNA, and mRNA nominated differential lncRNA with a likely miRNA sponge function. RESULTS: We identified 4371 differentially expressed features (mRNA = 3705; lncRNA = 666; FDR ≤ 5%). From these, the protein-coding gene TP53 was identified as an upstream regulator of differential lncRNAs PANDAR and PVT1 (p = 0.0237) and enrichment was detected in the development of carcinoma (p = 0.0176). Two lncRNAs (HCG11, PART1) were positively correlated with 342 mRNAs, and their correlation estimates diminish after adjusting for either of the target miRNAs: hsa-miR-490-3p, hsa-miR-129-5p. This suggests a likely sponge function for HCG11 and PART1. CONCLUSIONS: These findings identify differential lncRNAs with oncogenic features that are associated with G-CIMP phenotypes. Further investigation with controlled experiments is needed to confirm the molecular relationships

Noninvasive approaches to detect methylation-based markers to monitor gliomas

In this review, we summarize the current approaches used to detect glioma tissue-derived DNA methylation markers in liquid biopsy specimens with the aim to diagnose, prognosticate and potentially track treatment response and evolution of patients with gliomas

An early requirement for maternal FoxH1 during zebrafish gastrulation

AbstractThe Forkhead Box H1 (FoxH1) protein is a co-transcription factor recruited by phosphorylated Smad2 downstream of several TGFβs, including Nodal-related proteins. We have reassessed the function of zebrafish FoxH1 using antisense morpholino oligonucleotides (MOs). MOs targeting translation of foxH1 disrupt embryonic epiboly movements during gastrulation and cause death on the first day of development. The FoxH1 morphant phenotype is much more severe than that of zebrafish carrying foxh1/schmalspur (sur) DNA-binding domain mutations, FoxH1 splice-blocking morphants or other Nodal pathway mutants, and it cannot be altered by concomitant perturbations in Nodal signaling. Apart from disrupting epiboly, FoxH1 MO treatment disrupts convergence and internalization movements. Late gastrula-stage FoxH1 morphants exhibit delayed mesoderm and endoderm marker gene expression and failed patterning of the central nervous system. Probing FoxH1 morphant RNA by microarray, we identified a cohort of five keratin genes – cyt1, cyt2, krt4, krt8 and krt18 – that are normally transcribed in the embryo's enveloping layer (EVL) and which have significantly reduced expression in FoxH1-depleted embryos. Simultaneously disrupting these keratins with a mixture of MOs reproduces the FoxH1 morphant phenotype. Our studies thus point to an essential role for maternal FoxH1 and downstream keratins during gastrulation that is epistatic to Nodal signaling

Optimizing exosomal RNA isolation for RNA-Seq analyses of archival sera specimens

Exosomes are endosome-derived membrane vesicles that contain proteins, lipids, and nucleic acids. The exosomal transcriptome mediates intercellular communication, and represents an understudied reservoir of novel biomarkers for human diseases. Next-generation sequencing enables complex quantitative characterization of exosomal RNAs from diverse sources. However, detailed protocols describing exosome purification for preparation of exosomal RNA-sequence (RNA-Seq) libraries are lacking. Here we compared methods for isolation of exosomes and extraction of exosomal RNA from human cell-free serum, as well as strategies for attaining equal representation of samples within pooled RNA-Seq libraries. We compared commercial precipitation with ultracentrifugation for exosome purification and confirmed the presence of exosomes via both transmission electron microscopy and immunoblotting. Exosomal RNA extraction was compared using four different RNA purification methods. We determined the minimal starting volume of serum required for exosome preparation and showed that high quality exosomal RNA can be isolated from sera stored for over a decade. Finally, RNA-Seq libraries were successfully prepared with exosomal RNAs extracted from human cell-free serum, cataloguing both coding and non-coding exosomal transcripts. This method provides researchers with strategic options to prepare RNA-Seq libraries and compare RNA-Seq data quantitatively from minimal volumes of fresh and archival human cell-free serum for disease biomarker discovery

Nucleosome positioning and histone modifications define relationships between regulatory elements and nearby gene expression in breast epithelial cells

Abstract\ud \ud Background\ud The precise nature of how cell type specific chromatin structures at enhancer sites affect gene expression is largely unknown. Here we identified cell type specific enhancers coupled with gene expression in two different types of breast epithelial cells, HMEC (normal breast epithelial cells) and MDAMB231 (triple negative breast cancer cell line).\ud \ud \ud Results\ud Enhancers were defined by modified neighboring histones [using chromatin immunoprecipitation followed by sequencing (ChIP-seq)] and nucleosome depletion [using formaldehyde-assisted isolation of regulatory elements followed by sequencing (FAIRE-seq)]. Histone modifications at enhancers were related to the expression levels of nearby genes up to 750 kb away. These expression levels were correlated with enhancer status (poised or active), defined by surrounding histone marks. Furthermore, about fifty percent of poised and active enhancers contained nucleosome-depleted regions. We also identified response element motifs enriched at these enhancer sites that revealed key transcription factors (e.g. TP63) likely involved in regulating breast epithelial enhancer-mediated gene expression. By utilizing expression data, potential target genes of more than 600 active enhancers were identified. These genes were involved in proteolysis, epidermis development, cell adhesion, mitosis, cell cycle, and DNA replication.\ud \ud \ud Conclusions\ud These findings facilitate the understanding of epigenetic regulation specifically, such as the relationships between regulatory elements and gene expression and generally, how breast epithelial cellular phenotypes are determined by cell type specific enhancers.National Institutes of Health [R01 CA136924 to GAC, T32CA009320 to HN]Genetic Associations and Mechanisms in Oncology (GAME-ON

SIRT1 regulates Mxd1 during malignant melanoma progression

In a murine melanoma model, malignant transformation promoted by a sustained stress condition was causally related to increased levels of reactive oxygen species resulting in DNA damage and massive epigenetic alterations. Since the chromatin modifier Sirtuin-1 (SIRT1) is a protein attracted to double-stranded DNA break (DSB) sites and can recruit other components of the epigenetic machinery, we aimed to define the role of SIRT1 in melanomagenesis through our melanoma model. The DNA damage marker, gamma H2AX was found increased in melanocytes after 24 hours of deadhesion, accompanied by increased SIRT1 expression and decreased levels of its target, H4K16ac. Moreover, SIRT1 started to be associated to DNMT3B during the stress condition, and this complex was maintained along malignant progression. Mxd1 was identified by ChIP-seq among the DNA sequences differentially associated with SIRT1 during deadhesion and was shown to be a common target of both, SIRT1 and DNMT3B. In addition, Mxd1 was found downregulated from pre-malignant melanocytes to metastatic melanoma cells. Treatment with DNMT inhibitor 5AzaCdR reversed the Mxd1 expression. Sirt1 stable silencing increased Mxd1 mRNA expression and led to down-regulation of MYC targets, such as Cdkn1a, Bcl2 and Psen2, whose upregulation is associated with human melanoma aggressiveness and poor prognosis. We demonstrated a novel role of the stress responsive protein SIRT1 in malignant transformation of melanocytes associated with deadhesion. Mxd1 was identified as a new SIRT1 target gene. SIRT1 promoted Mxd1 silencing, which led to increased activity of MYC oncogene contributing to melanoma progression.FAPESP [2011/0166-38, 2011/12306-1, 2014/13663-0, 2015/07925-5, 2016/06488-3]DAAD [PKZ A/12/79134]FAPESP/BAYLAT [2012/51300-7]Univ Fed Sao Paulo UNIFESP, Dept Pharmacol, Ontogeny & Epigenet Lab, Sao Paulo, SP, BrazilUniv Sao Paulo, Ribeirao Preto Med Sch, Dept Genet, Ribeirao Preto, SP, BrazilFriedrich Alexander Univ Erlangen Nurnberg FAU, Inst Pathol, Expt Tumorpathol, Erlangen, GermanyFriedrich Alexander Univ Erlangen Nurnberg FAU, Dept Pediat & Adolescent Med, Erlangen, GermanyUniv Fed Sao Paulo UNIFESP, Dept Pharmacol, Ontogeny & Epigenet Lab, Sao Paulo, SP, BrazilFAPESP [2011/0166-38, 2011/12306-1, 2014/13663-0, 2015/07925-5, 2016/06488-3]DAAD [PKZ A/12/79134]FAPESP/BAYLAT [2012/51300-7]Web of Scienc

Super-Enhancer-Associated LncRNA UCA1 Interacts Directly with AMOT to Activate YAP Target Genes in Epithelial Ovarian Cancer

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of tumorigenesis, and yet their mechanistic roles remain challenging to characterize. Here, we integrate functional proteomics with lncRNA-interactome profiling to characterize Urothelial Cancer Associated 1 (UCA1), a candidate driver of ovarian cancer development. Reverse phase protein array (RPPA) analysis indicates that UCA1 activates transcription coactivator YAP and its target genes. In vivo RNA antisense purification (iRAP) of UCA1 interacting proteins identified angiomotin (AMOT), a known YAP regulator, as a direct binding partner. Loss-of-function experiments show that AMOT mediates YAP activation by UCA1, as UCA1 enhances the AMOT-YAP interaction to promote YAP dephosphorylation and nuclear translocation. Together, we characterize UCA1 as a lncRNA regulator of Hippo-YAP signaling and highlight the UCA1-AMOT-YAP signaling axis in ovarian cancer development

Comprehensive functional annotation of 77 prostate cancer risk loci.

Genome-wide association studies (GWAS) have revolutionized the field of cancer genetics, but the causal links between increased genetic risk and onset/progression of disease processes remain to be identified. Here we report the first step in such an endeavor for prostate cancer. We provide a comprehensive annotation of the 77 known risk loci, based upon highly correlated variants in biologically relevant chromatin annotations--we identified 727 such potentially functional SNPs. We also provide a detailed account of possible protein disruption, microRNA target sequence disruption and regulatory response element disruption of all correlated SNPs at r(2) ≥ 0.88%. 88% of the 727 SNPs fall within putative enhancers, and many alter critical residues in the response elements of transcription factors known to be involved in prostate biology. We define as risk enhancers those regions with enhancer chromatin biofeatures in prostate-derived cell lines with prostate-cancer correlated SNPs. To aid the identification of these enhancers, we performed genomewide ChIP-seq for H3K27-acetylation, a mark of actively engaged enhancers, as well as the transcription factor TCF7L2. We analyzed in depth three variants in risk enhancers, two of which show significantly altered androgen sensitivity in LNCaP cells. This includes rs4907792, that is in linkage disequilibrium (r(2) = 0.91) with an eQTL for NUDT11 (on the X chromosome) in prostate tissue, and rs10486567, the index SNP in intron 3 of the JAZF1 gene on chromosome 7. Rs4907792 is within a critical residue of a strong consensus androgen response element that is interrupted in the protective allele, resulting in a 56% decrease in its androgen sensitivity, whereas rs10486567 affects both NKX3-1 and FOXA-AR motifs where the risk allele results in a 39% increase in basal activity and a 28% fold-increase in androgen stimulated enhancer activity. Identification of such enhancer variants and their potential target genes represents a preliminary step in connecting risk to disease process