Evolutionary divergence of Firre localization and expression.

Long noncoding RNAs (lncRNAs) are rapidly evolving and thus typically poorly conserved in their sequences. How these sequence differences affect the characteristics and potential functions of lncRNAs with shared synteny remains unclear. Here we show that the syntenically conserved lncRNA Firre displays distinct expression and localization patterns in human and mouse. Single molecule RNA FISH reveals that in a range of cell lines, mouse Firre (mFirre) is predominantly nuclear, while human FIRRE (hFIRRE) is distributed between the cytoplasm and nucleus. This localization pattern is maintained in human/mouse hybrid cells expressing both human and mouse Firre, implying that the localization of the lncRNA is species autonomous. We find that the majority of hFIRRE transcripts in the cytoplasm are comprised of isoforms that are enriched in RRD repeats. We furthermore determine that in various tissues, mFirre is more highly expressed than its human counterpart. Our data illustrate that the rapid evolution of syntenic lncRNAs can lead to variations in lncRNA localization and abundance, which in turn may result in disparate lncRNA functions even in closely related species.

Stimulated emission depletion (STED) super resolution imaging of RNA- and protein-containing domains in fixed cells.

Super resolution microscopy has changed our capability to visualize and understand spatial arrangements of RNA- and protein-containing domains in individual cells. In a previous study, we described a novel lncRNA, Tumor-associated NBL2 transcript (TNBL), which originates from a primate specific macrosatellite repeat. We aimed to describe several aspects of TNBL lncRNA, with one focus being pinpointing its precise location in the nucleus, as well as visualizing its interactions with proteins to deduce its functionality. Using a combination of STimulated Emission Depletion (STED) super resolution microscopy, single molecule RNA (smRNA) FISH against TNBL, and immunofluorescence against SAM68 perinucleolar body, we resolved the spatial complexity of the interaction between TNBL aggregates and SAM68 bodies at the perinucleolar region. Here, we describe protocols for a step-by-step optimized smRNA FISH/IF and STED imaging, detailing parameter settings, and three-dimensional data analysis of spatial positioning of subnuclear structures. These protocols can be employed for single-cell imaging of complex nuclear RNA-protein structures.

Single-cell imaging reveals unexpected heterogeneity of telomerase reverse transcriptase expression across human cancer cell lines.

Telomerase is pathologically reactivated in most human cancers, where it maintains chromosomal telomeres and allows immortalization. Because telomerase reverse transcriptase (TERT) is usually the limiting component for telomerase activation, numerous studies have measured TERT mRNA levels in populations of cells or in tissues. In comparison, little is known about TERT expression at the single-cell and single-molecule level. To address this, we analyzed TERT expression across 10 human cancer lines using single-molecule RNA fluorescent in situ hybridization (FISH) and made several unexpected findings. First, there was substantial cell-to-cell variation in number of transcription sites and ratio of transcription sites to gene copies. Second, previous classification of lines as having monoallelic or biallelic TERT expression was found to be inadequate for capturing the TERT gene expression patterns. Finally, spliced TERT mRNA had primarily nuclear localization in cancer cells and induced pluripotent stem cells (iPSCs), in stark contrast to the expectation that spliced mRNA should be predominantly cytoplasmic. These data reveal unappreciated heterogeneity, complexity, and unconventionality in TERT expression across human cancer cells.

A novel long non-coding RNA from NBL2 pericentromeric macrosatellite forms a perinucleolar aggregate structure in colon cancer.

Primate-specific NBL2 macrosatellite is hypomethylated in several types of tumors, yet the consequences of this DNA hypomethylation remain unknown. We show that NBL2 conserved repeats are close to the centromeres of most acrocentric chromosomes. NBL2 associates with the perinucleolar region and undergoes severe demethylation in a subset of colorectal cancer (CRC). Upon DNA hypomethylation and histone acetylation, NBL2 repeats are transcribed in tumor cell lines and primary CRCs. NBL2 monomers exhibit promoter activity, and are contained within novel, non-polyA antisense lncRNAs, which we designated TNBL (Tumor-associated NBL2 transcript). TNBL is stable throughout the mitotic cycle, and in interphase nuclei preferentially forms a perinucleolar aggregate in the proximity of a subset of NBL2 loci. TNBL aggregates interact with the SAM68 perinucleolar body in a mirror-image cancer specific perinucleolar structure. TNBL binds with high affinity to several proteins involved in nuclear functions and RNA metabolism, such as CELF1 and NPM1. Our data unveil novel DNA and RNA structural features of a non-coding macrosatellite frequently altered in cancer.

VarID2 quantifies gene expression noise dynamics and unveils functional heterogeneity of ageing hematopoietic stem cells.

Variability of gene expression due to stochasticity of transcription or variation of extrinsic signals, termed biological noise, is a potential driving force of cellular differentiation. Utilizing single-cell RNA-sequencing, we develop VarID2 for the quantification of biological noise at single-cell resolution. VarID2 reveals enhanced nuclear versus cytoplasmic noise, and distinct regulatory modes stratified by correlation between noise, expression, and chromatin accessibility. Noise levels are minimal in murine hematopoietic stem cells (HSCs) and increase during differentiation and ageing. Differential noise identifies myeloid-biased Dlk1+ long-term HSCs in aged mice with enhanced quiescence and self-renewal capacity. VarID2 reveals noise dynamics invisible to conventional single-cell transcriptome analysis.

The lncRNA Sweetheart regulates compensatory cardiac hypertrophy after myocardial injury in murine males.

After myocardial infarction in the adult heart the remaining, non-infarcted tissue adapts to compensate the loss of functional tissue. This adaptation requires changes in gene expression networks, which are mostly controlled by transcription regulating proteins. Long non-coding transcripts (lncRNAs) are taking part in fine-tuning such gene programs. We describe and characterize the cardiomyocyte specific lncRNA Sweetheart RNA (Swhtr), an approximately 10 kb long transcript divergently expressed from the cardiac core transcription factor coding gene Nkx2-5. We show that Swhtr is dispensable for normal heart development and function but becomes essential for the tissue adaptation process after myocardial infarction in murine males. Re-expressing Swhtr from an exogenous locus rescues the Swhtr null phenotype. Genes that depend on Swhtr after cardiac stress are significantly occupied and therefore most likely regulated by NKX2-5. The Swhtr transcript interacts with NKX2-5 and disperses upon hypoxic stress in cardiomyocytes, indicating an auxiliary role of Swhtr for NKX2-5 function in tissue adaptation after myocardial injury.

Mosaic cis-regulatory evolution drives transcriptional partitioning of HERVH endogenous retrovirus in the human embryo.

The human endogenous retrovirus type-H (HERVH) family is expressed in the preimplantation embryo. A subset of these elements are specifically transcribed in pluripotent stem cells where they appear to exert regulatory activities promoting self-renewal and pluripotency. How HERVH elements achieve such transcriptional specificity remains poorly understood. To uncover the sequence features underlying HERVH transcriptional activity, we performed a phyloregulatory analysis of the long terminal repeats (LTR7) of the HERVH family, which harbor its promoter, using a wealth of regulatory genomics data. We found that the family includes at least eight previously unrecognized subfamilies that have been active at different timepoints in primate evolution and display distinct expression patterns during human embryonic development. Notably, nearly all HERVH elements transcribed in ESCs belong to one of the youngest subfamilies we dubbed LTR7up. LTR7 sequence evolution was driven by a mixture of mutational processes, including point mutations, duplications, and multiple recombination events between subfamilies, that led to transcription factor binding motif modules characteristic of each subfamily. Using a reporter assay, we show that one such motif, a predicted SOX2/3 binding site unique to LTR7up, is essential for robust promoter activity in induced pluripotent stem cells. Together these findings illuminate the mechanisms by which HERVH diversified its expression pattern during evolution to colonize distinct cellular niches within the human embryo.

Unravelling the impact of aging on the human endothelial lncRNA transcriptome.

The incidence and prevalence of cardiovascular disease is highest among the elderly. There is a need to further understand the mechanisms behind endothelial cell aging in order to achieve vascular rejuvenation and minimize the onset of age-related vascular diseases. Long non-coding RNAs (lncRNAs) have been proposed to regulate numerous processes in the human genome, yet their function in vascular aging and their therapeutic potential remain largely unknown. This is primarily because the majority of studies investigating the impact of aging on lncRNA expression heavily rely on in vitro studies based on replicative senescence. Here, using a unique collection of young and aged endothelial cells isolated from native human arteries, we sought to characterize the age-related alterations in lncRNA expression profiles. We were able to detect a total of 4463 lncRNAs expressed in the human endothelium from which ∼17% (798) were altered in advanced age. One of the most affected lncRNAs in aging was the primate-specific, Prostate Cancer Associated Transcript (PCAT) 14. In our follow up analysis, using single molecule RNA FISH, we showed that PCAT14 is relatively abundant, localized almost exclusively in the nucleus of young endothelial cells, and silenced in the aged endothelium. Functionally, our studies proposed that downregulation of PCAT14 alters endothelial cell transcription profile and cell functions including endothelial cell migration, sprouting and inflammatory responses in vitro. Taken together, our data highlight that endothelial cell aging correlates with altered expression of lncRNAs, which could impair the endothelial regenerative capacity and enhance inflammatory phenotypes.

Memory B-cell like chronic lymphocytic leukaemia is associated with specific methylation profile of WNT5A promoter and undetectable expression of WNT5A gene.

Genome methylation profiles define naïve-like (n-CLL), memory-like (m-CLL), and intermediate (i-CLL) subsets of chronic lymphocytic leukaemia (CLL). The profiles can be easily determined by the analysis of the five-CpG signature. m-CLL, i-CLL, and n-CLL with the good, intermediate, and poor prognoses, respectively, differ by the somatic hypermutation status of the immunoglobulin heavy chain variable gene (IGHV), a widely used prognostic predictor in CLL. We have previously shown that the expression of WNT5A, encoding a ROR1 ligand, distinguishes patients with the worse outcome within the prognostically favourable IGHV-mutated subgroup. To analyse the mechanisms controlling WNT5A expression, we investigated the methylation status of 54 CpG sites within the WNT5A promoter and its relation to the WNT5A gene expression. In a cohort of 59 CLL patients balanced for combinations of IGHV and WNT5A statuses, we identified three promoter CpG sites whose methylation level correlated with the WNT5A expression within the IGHV-mutated subgroup. Further, we complemented our data with the methylation status of the five-CpG signature. IGHV-mutated/WNT5A-negative and IGHV-mutated/WNT5A-positive cases overlapped with m­CLL and i­CLL methylation subgroups, respectively, while most IGHV­unmutated samples were assigned to n-CLL. Median methylation levels of all the three CpG sites in the WNT5A promoter were lowest in i-CLL. Finally, a detailed analysis of m-CLL and i-CLL showed that undetectable WNT5A expression predicts longer treatment-free survival with higher statistical significance than the classification according to the five-CpG signature. To conclude, a favourable m-CLL subgroup is associated with mutated IGHV and undetectable WNT5A expression due to its promoter methylation.

Nuclear compartmentalization of TERT mRNA and TUG1 lncRNA is driven by intron retention.

The spatial partitioning of the transcriptome in the cell is an important form of gene-expression regulation. Here, we address how intron retention influences the spatio-temporal dynamics of transcripts from two clinically relevant genes: TERT (Telomerase Reverse Transcriptase) pre-mRNA and TUG1 (Taurine-Upregulated Gene 1) lncRNA. Single molecule RNA FISH reveals that nuclear TERT transcripts uniformly and robustly retain specific introns. Our data suggest that the splicing of TERT retained introns occurs during mitosis. In contrast, TUG1 has a bimodal distribution of fully spliced cytoplasmic and intron-retained nuclear transcripts. We further test the functionality of intron-retention events using RNA-targeting thiomorpholino antisense oligonucleotides to block intron excision. We show that intron retention is the driving force for the nuclear compartmentalization of these RNAs. For both RNAs, altering this splicing-driven subcellular distribution has significant effects on cell viability. Together, these findings show that stable retention of specific introns can orchestrate spatial compartmentalization of these RNAs within the cell. This process reveals that modulating RNA localization via targeted intron retention can be utilized for RNA-based therapies.

Using antisense oligonucleotides for the physiological modulation of the alternative splicing of NF1 exon 23a during PC12 neuronal differentiation.

Neurofibromatosis Type 1 (NF1) is a genetic condition affecting approximately 1:3500 persons worldwide. The NF1 gene codes for neurofibromin protein, a GTPase activating protein (GAP) and a negative regulator of RAS. The NF1 gene undergoes alternative splicing of exon 23a (E23a) that codes for 21 amino acids placed at the center of the GAP related domain (GRD). E23a-containing type II neurofibromin exhibits a weaker Ras-GAP activity compared to E23a-less type I isoform. Exon E23a has been related with the cognitive impairment present in NF1 individuals. We designed antisense Phosphorodiamidate Morpholino Oligomers (PMOs) to modulate E23a alternative splicing at physiological conditions of gene expression and tested their impact during PC12 cell line neuronal differentiation. Results show that any dynamic modification of the natural ratio between type I and type II isoforms disturbed neuronal differentiation, altering the proper formation of neurites and deregulating both the MAPK/ERK and cAMP/PKA signaling pathways. Our results suggest an opposite regulation of these pathways by neurofibromin and the possible existence of a feedback loop sensing neurofibromin-related signaling. The present work illustrates the utility of PMOs to study alternative splicing that could be applied to other alternatively spliced genes in vitro and in vivo.

Identification of human long noncoding RNAs associated with nonalcoholic fatty liver disease and metabolic homeostasis.

A growing number of long noncoding RNAs (lncRNAs) have emerged as vital metabolic regulators. However, most human lncRNAs are nonconserved and highly tissue specific, vastly limiting our ability to identify human lncRNA metabolic regulators (hLMRs). In this study, we established a pipeline to identify putative hLMRs that are metabolically sensitive, disease relevant, and population applicable. We first progressively processed multilevel human transcriptome data to select liver lncRNAs that exhibit highly dynamic expression in the general population, show differential expression in a nonalcoholic fatty liver disease (NAFLD) population, and respond to dietary intervention in a small NAFLD cohort. We then experimentally demonstrated the responsiveness of selected hepatic lncRNAs to defined metabolic milieus in a liver-specific humanized mouse model. Furthermore, by extracting a concise list of protein-coding genes that are persistently correlated with lncRNAs in general and NAFLD populations, we predicted the specific function for each hLMR. Using gain- and loss-of-function approaches in humanized mice as well as ectopic expression in conventional mice, we validated the regulatory role of one nonconserved hLMR in cholesterol metabolism by coordinating with an RNA-binding protein, PTBP1, to modulate the transcription of cholesterol synthesis genes. Our work overcame the heterogeneity intrinsic to human data to enable the efficient identification and functional definition of disease-relevant human lncRNAs in metabolic homeostasis.

The Tug1 lncRNA locus is essential for male fertility.

BACKGROUND: Several long noncoding RNAs (lncRNAs) have been shown to function as components of molecular machines that play fundamental roles in biology. While the number of annotated lncRNAs in mammalian genomes has greatly expanded, studying lncRNA function has been a challenge due to their diverse biological roles and because lncRNA loci can contain multiple molecular modes that may exert function. RESULTS: We previously generated and characterized a cohort of 20 lncRNA loci knockout mice. Here, we extend this initial study and provide a more detailed analysis of the highly conserved lncRNA locus, taurine-upregulated gene 1 (Tug1). We report that Tug1-knockout male mice are sterile with underlying defects including a low number of sperm and abnormal sperm morphology. Because lncRNA loci can contain multiple modes of action, we wanted to determine which, if any, potential elements contained in the Tug1 genomic region have any activity. Using engineered mouse models and cell-based assays, we provide evidence that the Tug1 locus harbors two distinct noncoding regulatory activities, as a cis-DNA repressor that regulates neighboring genes and as a lncRNA that can regulate genes by a trans-based function. We also show that Tug1 contains an evolutionary conserved open reading frame that when overexpressed produces a stable protein which impacts mitochondrial membrane potential, suggesting a potential third coding function. CONCLUSIONS: Our results reveal an essential role for the Tug1 locus in male fertility and uncover evidence for distinct molecular modes in the Tug1 locus, thus highlighting the complexity present at lncRNA loci.

Emerging roles of macrosatellite repeats in genome organization and disease development.

Abundant repetitive DNA sequences are an enigmatic part of the human genome. Despite increasing evidence on the functionality of DNA repeats, their biologic role is still elusive and under frequent debate. Macrosatellites are the largest of the tandem DNA repeats, located on one or multiple chromosomes. The contribution of macrosatellites to genome regulation and human health was demonstrated for the D4Z4 macrosatellite repeat array on chromosome 4q35. Reduced copy number of D4Z4 repeats is associated with local euchromatinization and the onset of facioscapulohumeral muscular dystrophy. Although the role other macrosatellite families may play remains rather obscure, their diverse functionalities within the genome are being gradually revealed. In this review, we will outline structural and functional features of coding and noncoding macrosatellite repeats, and highlight recent findings that bring these sequences into the spotlight of genome organization and disease development.

Barcelona conference on epigenetics and cancer 2016 - beyond cancer genomes.

The Barcelona Conference on Epigenetics and Cancer (BCEC) entitled "Beyond Cancer Genomes" took place October 13th and 14th 2016 in Barcelona. The 2016 BCEC was the fourth edition of a series of annual conferences coordinated by Marcus Buschbeck and subsequently organized by leading research centers in Barcelona together with B•DEBATE, a joint initiative of BIOCAT and "La Caixa" Foundation. Salvador Aznar-Benitah, Eduard Batlle, and Raúl Méndez from the Institute for Research in Biomedicine in Barcelona selected the 2016 BCEC panel of speakers. As the title indicates, this year's conference expanded the epigenetic focus to include additional cancer-relevant topics, such as tumor heterogeneity and RNA regulation. Methods to develop therapeutic approaches on the basis of novel insights have been discussed in great detail. The conference has attracted 217 participants from 11 countries.

Helicase Lymphoid-Specific Enzyme Contributes to the Maintenance of Methylation of SST1 Pericentromeric Repeats That Are Frequently Demethylated in Colon Cancer and Associate with Genomic Damage.

DNA hypomethylation at repetitive elements accounts for the genome-wide DNA hypomethylation common in cancer, including colorectal cancer (CRC). We identified a pericentromeric repeat element called SST1 frequently hypomethylated (>5% demethylation compared with matched normal tissue) in several cancers, including 28 of 128 (22%) CRCs. SST1 somatic demethylation associated with genome damage, especially in tumors with wild-type TP53. Seven percent of the 128 CRCs exhibited a higher ("severe") level of demethylation (≥10%) that co-occurred with TP53 mutations. SST1 demethylation correlated with distinct histone marks in CRC cell lines and primary tumors: demethylated SST1 associated with high levels of the repressive histone 3 lysine 27 trimethylation (H3K27me3) mark and lower levels of histone 3 lysine 9 trimethylation (H3K9me3). Furthermore, induced demethylation of SST1 by 5-aza-dC led to increased H3K27me3 and reduced H3K9me3. Thus, in some CRCs, SST1 demethylation reflects an epigenetic reprogramming associated with changes in chromatin structure that may affect chromosomal integrity. The chromatin remodeler factor, the helicase lymphoid-specific (HELLS) enzyme, called the "epigenetic guardian of repetitive elements", interacted with SST1 as shown by chromatin immunoprecipitation, and down-regulation of HELLS by shRNA resulted in demethylation of SST1 in vitro. Altogether these results suggest that HELLS contributes to SST1 methylation maintenance. Alterations in HELLS recruitment and function could contribute to the somatic demethylation of SST1 repeat elements undergone before and/or during CRC pathogenesis.