Global Run-On Sequencing (GRO-Seq)

HIV-1 Vpr induces ciTRAN to prevent transcriptional repression of the provirus

The functional consequences of circular RNA (circRNA) expression on HIV-1 replication are largely unknown. Using a customized protocol involving direct RNA nanopore sequencing, here, we captured circRNAs from HIV-1-infected T cells and identified ciTRAN, a circRNA that modulates HIV-1 transcription. We found that HIV-1 infection induces ciTRAN expression in a Vpr-dependent manner and that ciTRAN interacts with SRSF1, a protein known to repress HIV-1 transcription. Our results suggest that HIV-1 hijacks ciTRAN to exclude serine/arginine-rich splicing factor 1 (SRSF1) from the viral transcriptional complex, thereby promoting efficient viral transcription. In addition, we demonstrate that an SRSF1-inspired mimic can inhibit viral transcription regardless of ciTRAN induction. The hijacking of a host circRNA thus represents a previously unknown facet of primate lentiviruses in overcoming transmission bottlenecks.

A transcription-independent mechanism determines rapid periodic fluctuations of BRCA1 expression

BRCA1 expression is highly regulated to prevent genomic instability and tumorigenesis. Dysregulation of BRCA1 expression correlates closely with sporadic basal-like breast cancer and ovarian cancer. The most significant characteristic of BRCA1 regulation is periodic expression fluctuation throughout the cell cycle, which is important for the orderly progression of different DNA repair pathways throughout the various cell cycle phases and for further genomic stability. However, the underlying mechanism driving this phenomenon is poorly understood. Here, we demonstrate that RBM10-mediated RNA alternative splicing coupled to nonsense-mediated mRNA decay (AS-NMD), rather than transcription, determines the periodic fluctuations in G1/S-phase BRCA1 expression. Furthermore, AS-NMD broadly regulates the expression of period genes, such as DNA replication-related genes, in an uneconomical but more rapid manner. In summary, we identified an unexpected posttranscriptional mechanism distinct from canonical processes that mediates the rapid regulation of BRCA1 as well as other period gene expression during the G1/S-phase transition and provided insights into potential targets for cancer therapy.

Predictive model of transcriptional elongation control identifies trans regulatory factors from chromatin signatures

Promoter-proximal Polymerase II (Pol II) pausing is a key rate-limiting step for gene expression. DNA and RNA-binding trans-acting factors regulating the extent of pausing have been identified. However, we lack a quantitative model of how interactions of these factors determine pausing, therefore the relative importance of implicated factors is unknown. Moreover, previously unknown regulators might exist. Here we address this gap with a machine learning model that accurately predicts the extent of promoter-proximal Pol II pausing from large-scale genome and transcriptome binding maps and gene annotation and sequence composition features. We demonstrate high accuracy and generalizability of the model by validation on an independent cell line which reveals the model's cell line agnostic character. Model interpretation in light of prior knowledge about molecular functions of regulatory factors confirms the interconnection of pausing with other RNA processing steps. Harnessing underlying feature contributions, we assess the relative importance of each factor, quantify their predictive effects and systematically identify previously unknown regulators of pausing. We additionally identify 16 previously unknown 7SK ncRNA interacting RNA-binding proteins predictive of pausing. Our work provides a framework to further our understanding of the regulation of the critical early steps in transcriptional elongation.

The metazoan landscape of mitochondrial DNA gene order and content is shaped by selection and affects mitochondrial transcription

Mitochondrial DNA (mtDNA) harbors essential genes in most metazoans, yet the regulatory impact of the multiple evolutionary mtDNA rearrangements has been overlooked. Here, by analyzing mtDNAs from ~8000 metazoans we found high gene content conservation (especially of protein and rRNA genes), and codon preferences for mtDNA-encoded tRNAs across most metazoans. In contrast, mtDNA gene order (MGO) was selectively constrained within but not between phyla, yet certain gene stretches (ATP8-ATP6, ND4-ND4L) were highly conserved across metazoans. Since certain metazoans with different MGOs diverge in mtDNA transcription, we hypothesized that evolutionary mtDNA rearrangements affected mtDNA transcriptional patterns. As a first step to test this hypothesis, we analyzed available RNA-seq data from 53 metazoans. Since polycistron mtDNA transcripts constitute a small fraction of the steady-state RNA, we enriched for polycistronic boundaries by calculating RNA-seq read densities across junctions between gene couples encoded either by the same strand (SSJ) or by different strands (DSJ). We found that organisms whose mtDNA is organized in alternating reverse-strand/forward-strand gene blocks (mostly arthropods), displayed significantly reduced DSJ read counts, in contrast to organisms whose mtDNA genes are preferentially encoded by one strand (all chordates). Our findings suggest that mtDNA rearrangements are selectively constrained and likely impact mtDNA regulation.

Non-coding RNAs in human health and disease: potential function as biomarkers and therapeutic targets

Human diseases have been a critical threat from the beginning of human history. Knowing the origin, course of action and treatment of any disease state is essential. A microscopic approach to the molecular field is a more coherent and accurate way to explore the mechanism, progression, and therapy with the introduction and evolution of technology than a macroscopic approach. Non-coding RNAs (ncRNAs) play increasingly important roles in detecting, developing, and treating all abnormalities related to physiology, pathology, genetics, epigenetics, cancer, and developmental diseases. Noncoding RNAs are becoming increasingly crucial as powerful, multipurpose regulators of all biological processes. Parallel to this, a rising amount of scientific information has revealed links between abnormal noncoding RNA expression and human disorders. Numerous non-coding transcripts with unknown functions have been found in addition to advancements in RNA-sequencing methods. Non-coding linear RNAs come in a variety of forms, including circular RNAs with a continuous closed loop (circRNA), long non-coding RNAs (lncRNA), and microRNAs (miRNA). This comprises specific information on their biogenesis, mode of action, physiological function, and significance concerning disease (such as cancer or cardiovascular diseases and others). This study review focuses on non-coding RNA as specific biomarkers and novel therapeutic targets.

Non-coding RNAs in immunoregulation and autoimmunity: Technological advances and critical limitations

Immune cell function is critically dependent on precise control over transcriptional output from the genome. In this respect, integration of environmental signals that regulate gene expression, specifically by transcription factors, enhancer DNA elements, genome topography and non-coding RNAs (ncRNAs), are key components. The first three have been extensively investigated. Even though non-coding RNAs represent the vast majority of cellular RNA species, this class of RNA remains historically understudied. This is partly because of a lag in technological and bioinformatic innovations specifically capable of identifying and accurately measuring their expression. Nevertheless, recent progress in this domain has enabled a profusion of publications identifying novel sub-types of ncRNAs and studies directly addressing the function of ncRNAs in human health and disease. Many ncRNAs, including circular and enhancer RNAs, have now been demonstrated to play key functions in the regulation of immune cells and to show associations with immune-mediated diseases. Some ncRNAs may function as biomarkers of disease, aiding in diagnostics and in estimating response to treatment, while others may play a direct role in the pathogenesis of disease. Importantly, some are relatively stable and are amenable to therapeutic targeting, for example through gene therapy. Here, we provide an overview of ncRNAs and review technological advances that enable their study and hold substantial promise for the future. We provide context-specific examples by examining the associations of ncRNAs with four prototypical human autoimmune diseases, specifically rheumatoid arthritis, psoriasis, inflammatory bowel disease and multiple sclerosis. We anticipate that the utility and mechanistic roles of these ncRNAs in autoimmunity will be further elucidated in the near future.

TRIM24 is an insulin-responsive regulator of P-bodies

Insulin is a potent inducer of mRNA transcription and translation, contributing to metabolic regulation. Insulin has also been suggested to regulate mRNA stability through the processing body (P-body) molecular machinery. However, whether and how insulin regulates mRNA stability via P-bodies is not clear. Here we show that the E3-ligase TRIM24 is a critical factor linking insulin signalling to P-bodies. Upon insulin stimulation, protein kinase B (PKB, also known as Akt) phosphorylates TRIM24 and stimulates its shuttling from the nucleus into the cytoplasm. TRIM24 interacts with several critical components of P-bodies in the cytoplasm, promoting their polyubiquitylation, which consequently stabilises Pparγ mRNA. Inactivation of TRIM24 E3-ligase activity or prevention of its phosphorylation via knockin mutations in mice promotes hepatic Pparγ degradation via P-bodies. Consequently, both knockin mutations alleviate hepatosteatosis in mice fed on a high-fat diet. Our results demonstrate the critical role of TRIM24 in linking insulin signalling to P-bodies and have therapeutic implications for the treatment of hepatosteatosis.

Hypoxia classifier for transcriptome datasets

Molecular gene signatures are useful tools to characterize the physiological state of cell populations, but most have developed under a narrow range of conditions and cell types and are often restricted to a set of gene identities. Focusing on the transcriptional response to hypoxia, we aimed to generate widely applicable classifiers sourced from the results of a meta-analysis of 69 differential expression datasets which included 425 individual RNA-seq experiments from 33 different human cell types exposed to different degrees of hypoxia (0.1-5%[Formula: see text]) for 2-48 h. The resulting decision trees include both gene identities and quantitative boundaries, allowing for easy classification of individual samples without control or normoxic reference. Each tree is composed of 3-5 genes mostly drawn from a small set of just 8 genes (EGLN1, MIR210HG, NDRG1, ANKRD37, TCAF2, PFKFB3, BHLHE40, and MAFF). In spite of their simplicity, these classifiers achieve over 95% accuracy in cross validation and over 80% accuracy when applied to additional challenging datasets. Our results indicate that the classifiers are able to identify hypoxic tumor samples from bulk RNAseq and hypoxic regions within tumor from spatially resolved transcriptomics datasets. Moreover, application of the classifiers to histological sections from normal tissues suggest the presence of a hypoxic gene expression pattern in the kidney cortex not observed in other normoxic organs. Finally, tree classifiers described herein outperform traditional hypoxic gene signatures when compared against a wide range of datasets. This work describes a set of hypoxic gene signatures, structured as simple decision tress, that identify hypoxic samples and regions with high accuracy and can be applied to a broad variety of gene expression datasets and formats.

Functional annotation of breast cancer risk loci: current progress and future directions

Genome-wide association studies coupled with large-scale replication and fine-scale mapping studies have identified more than 150 genomic regions that are associated with breast cancer risk. Here, we review efforts to translate these findings into a greater understanding of disease mechanism. Our review comes in the context of a recently published fine-scale mapping analysis of these regions, which reported 352 independent signals and a total of 13,367 credible causal variants. The vast majority of credible causal variants map to noncoding DNA, implicating regulation of gene expression as the mechanism by which functional variants influence risk. Accordingly, we review methods for defining candidate-regulatory sequences, methods for identifying putative target genes and methods for linking candidate-regulatory sequences to putative target genes. We provide a summary of available data resources and identify gaps in these resources. We conclude that while much work has been done, there is still much to do. There are, however, grounds for optimism; combining statistical data from fine-scale mapping with functional data that are more representative of the normal "at risk" breast, generated using new technologies, should lead to a greater understanding of the mechanisms that influence an individual woman's risk of breast cancer.

HOTTIP-dependent R-loop formation regulates CTCF boundary activity and TAD integrity in leukemia

HOTTIP lncRNA is highly expressed in acute myeloid leukemia (AML) driven by MLL rearrangements or NPM1 mutations to mediate HOXA topologically associated domain (TAD) formation and drive aberrant transcription. However, the mechanism through which HOTTIP accesses CCCTC-binding factor (CTCF) chromatin boundaries and regulates CTCF-mediated genome topology remains unknown. Here, we show that HOTTIP directly interacts with and regulates a fraction of CTCF-binding sites (CBSs) in the AML genome by recruiting CTCF/cohesin complex and R-loop-associated regulators to form R-loops. HOTTIP-mediated R-loops reinforce the CTCF boundary and facilitate formation of TADs to drive gene transcription. Either deleting CBS or targeting RNase H to eliminate R-loops in the boundary CBS of β-catenin TAD impaired CTCF boundary activity, inhibited promoter/enhancer interactions, reduced β-catenin target expression, and mitigated leukemogenesis in xenograft mouse models with aberrant HOTTIP expression. Thus, HOTTIP-mediated R-loop formation directly reinforces CTCF chromatin boundary activity and TAD integrity to drive oncogene transcription and leukemia development.

CDK12 inhibition reduces abnormalities in cells from patients with myotonic dystrophy and in a mouse model

Myotonic dystrophy type 1 (DM1) is an RNA-based disease with no current treatment. It is caused by a transcribed CTG repeat expansion within the 3' untranslated region of the dystrophia myotonica protein kinase (DMPK) gene. Mutant repeat expansion transcripts remain in the nuclei of patients' cells, forming distinct microscopically detectable foci that contribute substantially to the pathophysiology of the condition. Here, we report small-molecule inhibitors that remove nuclear foci and have beneficial effects in the HSA^LR mouse model, reducing transgene expression, leading to improvements in myotonia, splicing, and centralized nuclei. Using chemoproteomics in combination with cell-based assays, we identify cyclin-dependent kinase 12 (CDK12) as a druggable target for this condition. CDK12 is a protein elevated in DM1 cell lines and patient muscle biopsies, and our results showed that its inhibition led to reduced expression of repeat expansion RNA. Some of the inhibitors identified in this study are currently the subject of clinical trials for other indications and provide valuable starting points for a drug development program in DM1.

Targeting Polyadenylation for Retention of RNA at Chromatin

The various steps of RNA polymerase II transcription, including transcription initiation, splicing, and termination, are interlinked and tightly coordinated. Efficient 3'end processing is defined by sequence motifs emerging in the nascent-transcribed RNA strand and the cotranscriptional binding of regulatory proteins. The processing of a mature 3'end consists of cleavage and polyadenylation and is coupled with RNA polymerase II transcription termination and the dissociation of the nascent RNA transcript from the chromatin-associated transcriptional template. The subcellular and subnuclear topological specificity of the various RNA species is important for their functions. For instance, the formation of RNA-binding protein interactions, critical for the final outcome of gene expression, may require the nucleoplasmic fully spliced and polyadenylated form of an RNA transcript. Thus, interfering with the critical step of transcription termination and 3'end formation provides a means for assaying the functional potential of a given RNA of interest.In this protocol, we describe a method for blocking 3'end processing of the nascent RNA transcript, by using RNase H-inactive antisense oligonucleotides targeting cleavage and polyadenylation, delivered via transient transfection in cell culture.

Approaches to Identify and Characterise the Post-Transcriptional Roles of lncRNAs in Cancer

It is becoming increasingly evident that the non-coding genome and transcriptome exert great influence over their coding counterparts through complex molecular interactions. Among non-coding RNAs (ncRNA), long non-coding RNAs (lncRNAs) in particular present increased potential to participate in dysregulation of post-transcriptional processes through both RNA and protein interactions. Since such processes can play key roles in contributing to cancer progression, it is desirable to continue expanding the search for lncRNAs impacting cancer through post-transcriptional mechanisms. The sheer diversity of mechanisms requires diverse resources and methods that have been developed and refined over the past decade. We provide an overview of computational resources as well as proven low-to-high throughput techniques to enable identification and characterisation of lncRNAs in their complex interactive contexts. As more cancer research strategies evolve to explore the non-coding genome and transcriptome, we anticipate this will provide a valuable primer and perspective of how these technologies have matured and will continue to evolve to assist researchers in elucidating post-transcriptional roles of lncRNAs in cancer.

Fetal sex modulates placental microRNA expression, potential microRNA-mRNA interactions, and levels of amino acid transporter expression and substrates: INFAT study subpopulation analysis of n-3 LCPUFA intervention during pregnancy and associations with offspring body composition

Background

Previously, we revealed sexually dimorphic mRNA expression and responsiveness to maternal dietary supplementation with n-3 long-chain polyunsaturated fatty acids (LCPUFA) in placentas from a defined INFAT study subpopulation. Here, we extended these analyses and explored the respective placental microRNA expression, putative microRNA-mRNA interactions, and downstream target processes as well as their associations with INFAT offspring body composition.

Results

We performed explorative placental microRNA profiling, predicted microRNA-mRNA interactions by bioinformatics, validated placental target microRNAs and their putative targets by RT-qPCR and western blotting, and measured amino acid levels in maternal and offspring cord blood plasma and placenta. microRNA, mRNA, protein, and amino acid levels were associated with each other and with offspring body composition from birth to 5 years of age. Forty-six differentially regulated microRNAs were found. Validations identified differential expression for microRNA-99a (miR-99a) and its predicted target genes mTOR, SLC7A5, encoding L-type amino acid transporter 1 (LAT1), and SLC6A6, encoding taurine transporter (TauT), and their prevailing significant sexually dimorphic regulation. Target mRNA levels were mostly higher in placentas from control male than from female offspring, whereas respective n-3 LCPUFA responsive target upregulation was predominantly found in female placentas, explaining the rather balanced expression levels between the sexes present only in the intervention group. LAT1 and TauT substrates tryptophan and taurine, respectively, were significantly altered in both maternal plasma at 32 weeks’ gestation and cord plasma following intervention, but not in the placenta. Several significant associations were observed for miR-99a, mTOR mRNA, SLC7A5 mRNA, and taurine and tryptophan in maternal and cord plasma with offspring body composition at birth, 1 year, 3 and 5 years of age.

Conclusions

Our data suggest that the analyzed targets may be part of a sexually dimorphic molecular regulatory network in the placenta, possibly modulating gene expression per se and/or counteracting n-3 LCPUFA responsive changes, and thereby stabilizing respective placental and fetal amino acid levels. Our data propose placental miR-99, SLC7A5 mRNA, and taurine and tryptophan levels in maternal and fetal plasma as potentially predictive biomarkers for offspring body composition.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00345-x.

Male gametophyte development in flowering plants: A story of quarantine and sacrifice

Over 160 years ago, scientists made the first microscopic observations of angiosperm pollen. Unlike in animals, male meiosis in angiosperms produces a haploid microspore that undergoes one asymmetric division to form a vegetative cell and a generative cell. These two cells have distinct fates: the vegetative cell exits the cell cycle and elongates to form a tip-growing pollen tube; the generative cell divides once more in the pollen grain or within the growing pollen tube to form a pair of sperm cells. The concept that male germ cells are less active than the vegetative cell came from biochemical analyses and pollen structure anatomy early in the last century and is supported by the pollen transcriptome data of the last decade. However, the mechanism of how and when the transcriptional repression in male germ cells occurs is still not fully understood. In this review, we provide a brief account of the cytological and metabolic differentiation between the vegetative cell and male germ cells, with emphasis on the role of temporary callose walls, dynamic nuclear pore density, transcription repression, and histone variants. We further discuss the intercellular movement of small interfering RNA (siRNA) derived from transposable elements (TEs) and reexamine the function of TE expression in male germ cells.

New tools for 'ZEBRA-FISHING'

Zebrafish has been established as a classical model for developmental studies, yet in the past years, with the explosion of novel technological methods, the use of zebrafish as a model has expanded. One of the prominent fields that took advantage of zebrafish as a model organism early on is hematopoiesis, the process of blood cell generation from hematopoietic stem and progenitor cells (HSPCs). In zebrafish, HSPCs are born early during development in the aorta-gonad-mesonephros region and then translocate to the caudal hematopoietic tissue, where they expand and finally take residence in the kidney marrow. This journey is tightly regulated at multiple levels from extracellular signals to chromatin. In order to delineate the mechanistic underpinnings of this process, next-generation sequencing techniques could be an important ally. Here, we describe genome-wide approaches that have been undertaken to delineate zebrafish hematopoiesis.

An optimised method for intact nuclei isolation from diatoms

Due to their abundance in the oceans, their extraordinary biodiversity and the increasing use for biotech applications, the study of diatom biology is receiving more and more attention in the recent years. One of the limitations in developing molecular tools for diatoms lies in the peculiar nature of their cell wall, that is made of silica and organic molecules and that hinders the application of standard methods for cell lysis required, for example, to extract organelles. In this study we present a protocol for intact nuclei isolation from diatoms that was successfully applied to three different species: two pennates, Pseudo-nitzschia multistriata and Phaeodactylum tricornutum, and one centric diatom species, Chaetoceros diadema. Intact nuclei were extracted by treatment with acidified NH4F solution combined to low intensity sonication pulses and separated from cell debris via FAC-sorting upon incubation with SYBR Green. Microscopy observations confirmed the integrity of isolated nuclei and high sensitivity DNA electrophoresis showed that genomic DNA extracted from isolated nuclei has low degree of fragmentation. This protocol has proved to be a flexible and versatile method to obtain intact nuclei preparations from different diatom species and it has the potential to speed up applications such as epigenetic explorations as well as single cell ("single nuclei") genomics, transcriptomics and proteomics in different diatom species.

Nucleosome landscape reflects phenotypic differences in Trypanosoma cruzi life forms

Trypanosoma cruzi alternates between replicative and nonreplicative life forms, accompanied by a shift in global transcription levels and by changes in the nuclear architecture, the chromatin proteome and histone posttranslational modifications. To gain further insights into the epigenetic regulation that accompanies life form changes, we performed genome-wide high-resolution nucleosome mapping using two T. cruzi life forms (epimastigotes and cellular trypomastigotes). By combining a powerful pipeline that allowed us to faithfully compare nucleosome positioning and occupancy, more than 125 thousand nucleosomes were mapped, and approximately 20% of them differed between replicative and nonreplicative forms. The nonreplicative forms have less dynamic nucleosomes, possibly reflecting their lower global transcription levels and DNA replication arrest. However, dynamic nucleosomes are enriched at nonreplicative regulatory transcription initiation regions and at multigenic family members, which are associated with infective-stage and virulence factors. Strikingly, dynamic nucleosome regions are associated with GO terms related to nuclear division, translation, gene regulation and metabolism and, notably, associated with transcripts with different expression levels among life forms. Finally, the nucleosome landscape reflects the steady-state transcription expression: more abundant genes have a more deeply nucleosome-depleted region at putative 5' splice sites, likely associated with trans-splicing efficiency. Taken together, our results indicate that chromatin architecture, defined primarily by nucleosome positioning and occupancy, reflects the phenotypic differences found among T. cruzi life forms despite the lack of a canonical transcriptional control context.

Rapid and Scalable Profiling of Nascent RNA with fastGRO

Genome-wide profiling of nascent RNA has become a fundamental tool to study transcription regulation. Unlike steady-state RNA-sequencing (RNA-seq), nascent RNA profiling mirrors real-time activity of RNA polymerases and provides an accurate readout of transcriptome-wide variations. Some species of nuclear RNAs (i.e., large intergenic noncoding RNAs [lincRNAs] and eRNAs) have a short half-life and can only be accurately gauged by nascent RNA techniques. Furthermore, nascent RNA-seq detects post-cleavage RNA at termination sites and promoter-associated antisense RNAs, providing insights into RNA polymerase II (RNAPII) dynamics and processivity. Here, we present a run-on assay with 4-thio ribonucleotide (4-S-UTP) labeling, followed by reversible biotinylation and affinity purification via streptavidin. Our protocol allows streamlined sample preparation within less than 3 days. We named the technique fastGRO (fast Global Run-On). We show that fastGRO is highly reproducible and yields a more complete and extensive coverage of nascent RNA than comparable techniques can. Importantly, we demonstrate that fastGRO is scalable and can be performed with as few as 0.5 × 10⁶ cells.

Barbieri et al. developed fastGRO, a nascent RNA-sequencing technique based on nuclear run-on. Using a streamlined, under-3-days protocol, fastGRO tracks the activity of RNA polymerase for differential gene expression analysis, polymerase kinetic studies, and profiling of lowly expressed and unstable RNA species. A low-input fastGRO protocol profiles nascent RNA in as little as 0.5 × 10⁶ cells.

SIX6 is a TAL1-regulated transcription factor in T-ALL and associated with inferior outcome

T-cell acute lymphoblastic leukemia (T-ALL) is a hematological malignancy driven by abnormal activity of transcription factors. Here we report an aberrant expression of the developmental transcription factor SIX6 in the TAL1-subtype of T-ALL. Our results demonstrate that the binding of TAL1 and GATA3 transcription factors into an upstream enhancer element directly regulates SIX6 expression. High expression of SIX6 was associated with inferior event-free survival within three independent patient cohorts. At a functional level, CRISPR-Cas9-mediated knockout of the SIX6 gene in TAL1 positive Jurkat cells induced changes in genes associated with the mTOR-, K-RAS-, and TNFα-related molecular signatures but did not impair cell proliferation or viability. There was also no acceleration of T-ALL development within a Myc driven zebrafish tumor model in vivo. Taken together, our results show that SIX6 belongs to the TAL1 regulatory gene network in T-ALL but is alone insufficient to influence the development or maintenance of T-ALL.

Noncoding RNAs in Cardiovascular Disease: Current Knowledge, Tools and Technologies for Investigation, and Future Directions: A Scientific Statement From the American Heart Association

BACKGROUND

The discovery that much of the non-protein-coding genome is transcribed and plays a diverse functional role in fundamental cellular processes has led to an explosion in the development of tools and technologies to investigate the role of these noncoding RNAs in cardiovascular health. Furthermore, identifying noncoding RNAs for targeted therapeutics to treat cardiovascular disease is an emerging area of research. The purpose of this statement is to review existing literature, offer guidance on tools and technologies currently available to study noncoding RNAs, and identify areas of unmet need.

METHODS

The writing group used systematic literature reviews (including MEDLINE, Web of Science through 2018), expert opinion/statements, analyses of databases and computational tools/algorithms, and review of current clinical trials to provide a broad consensus on the current state of the art in noncoding RNA in cardiovascular disease.

RESULTS

Significant progress has been made since the initial studies focusing on the role of miRNAs (microRNAs) in cardiovascular development and disease. Notably, recent progress on understanding the role of novel types of noncoding small RNAs such as snoRNAs (small nucleolar RNAs), tRNA (transfer RNA) fragments, and Y-RNAs in cellular processes has revealed a noncanonical function for many of these molecules. Similarly, the identification of long noncoding RNAs that appear to play an important role in cardiovascular disease processes, coupled with the development of tools to characterize their interacting partners, has led to significant mechanistic insight. Finally, recent work has characterized the unique role of extracellular RNAs in mediating intercellular communication and their potential role as biomarkers.

CONCLUSIONS

The rapid expansion of tools and pipelines for isolating, measuring, and annotating these entities suggests that caution in interpreting results is warranted until these methodologies are rigorously validated. Most investigators have focused on investigating the functional role of single RNA entities, but studies suggest complex interaction between different RNA molecules. The use of network approaches and advanced computational tools to understand the interaction of different noncoding RNA species to mediate a particular phenotype may be required to fully comprehend the function of noncoding RNAs in mediating disease phenotypes.

An Overview of Non-coding RNAs and Cardiovascular System

Cardiovascular disease management and timely diagnosis remain a major dilemma. Delineating molecular mechanisms of cardiovascular diseases is opening horizon in the field of molecular medicines and in the development of early diagnostic markers. Non-coding RNAs are the highly functional and vibrant nucleic acids and are known to be involved in the regulation of endothelial cells, vascular and smooth muscles cells, cardiac metabolism, ischemia, inflammation and many processes in cardiovascular system. This chapter is comprehensively focusing on the overview of the non-coding RNAs including their discovery, generation, classification and functional regulation. In addition, overview regarding different non-coding RNAs as long non-coding, siRNAs and miRNAs involvement in the cardiovascular diseases is also addressed. Detailed functional analysis of this vast group of highly regulatory molecules will be promising for shaping future drug discoveries.

EPIGENE: genome-wide transcription unit annotation using a multivariate probabilistic model of histone modifications

BACKGROUND

Understanding the transcriptome is critical for explaining the functional as well as regulatory roles of genomic regions. Current methods for the identification of transcription units (TUs) use RNA-seq that, however, require large quantities of mRNA rendering the identification of inherently unstable TUs, e.g. miRNA precursors, difficult. This problem can be alleviated by chromatin-based approaches due to a correlation between histone modifications and transcription.

RESULTS

Here, we introduce EPIGENE, a novel chromatin segmentation method for the identification of active TUs using transcription-associated histone modifications. Unlike the existing chromatin segmentation approaches, EPIGENE uses a constrained, semi-supervised multivariate hidden Markov model (HMM) that models the observed combination of histone modifications using a product of independent Bernoulli random variables, to identify active TUs. Our results show that EPIGENE can identify genome-wide TUs in an unbiased manner. EPIGENE-predicted TUs show an enrichment of RNA Polymerase II at the transcription start site and in gene body indicating that they are indeed transcribed. Comprehensive validation using existing annotations revealed that 93% of EPIGENE TUs can be explained by existing gene annotations and 5% of EPIGENE TUs in HepG2 can be explained by microRNA annotations. EPIGENE outperformed the existing RNA-seq-based approaches in TU prediction precision across human cell lines. Finally, we identified 232 novel TUs in K562 and 43 novel cell-specific TUs all of which were supported by RNA Polymerase II ChIP-seq and Nascent RNA-seq data.

CONCLUSION

We demonstrate the applicability of EPIGENE to identify genome-wide active TUs and to provide valuable information about unannotated TUs. EPIGENE is an open-source method and is freely available at: https://github.com/imbbLab/EPIGENE.

Retinoic Acid-Regulated Target Genes During Development: Integrative Genomics Analysis

Retinoic acid (RA), a major natural active metabolite of vitamin A (VA) is well known to play critical roles in embryonic development. The effects of RA are mediated by nuclear receptors (RARs), which regulate the expression of gene batteries involved in cell growth and differentiation. Since the early 1990s several laboratories have focused on understanding how RA-regulated genes and RAR binding sites operate by studying the differentiation of embryonal carcinoma cells and embryonic stem cells. The development of hybridization-based microarray technology and high performance software analysis programs has allowed the characterization of thousands of RA-regulated genes. During the two last decades, publication of the genome sequence of various organisms has allowed advances in massive parallel sequencing and bioinformatics analysis of genome-wide data sets. These new generation sequencing (NGS) technologies have revolutionized the field by providing a global integrated picture of RA-regulated gene networks and the regulatory programs involved in cell fate decisions during embryonal carcinoma and embryonic stem cells differentiation. Now the challenge is to reconstruct the RA-regulated gene networks at the single cell level during the development of specialized embryonic tissues.

The Tumor Suppressor ARID1A Controls Global Transcription via Pausing of RNA Polymerase II

AT-rich interactive domain-containing proteins 1A and 1B (ARID1A and ARID1B) are mutually exclusive subunits of the chromatin remodeler SWI/SNF. ARID1A is the most frequently mutated chromatin regulator across all cancers, and ovarian clear cell carcinoma (OCCC) carries the highest prevalence of ARID1A mutations (∼57%). Despite evidence implicating ARID1A in tumorigenesis, the mechanism remains elusive. Here, we demonstrate that ARID1A binds active regulatory elements in OCCC. Depletion of ARID1A represses RNA polymerase II (RNAPII) transcription but results in modest changes to accessibility. Specifically, pausing of RNAPII is severely impaired after loss of ARID1A. Compromised pausing results in transcriptional dysregulation of active genes, which is compensated by upregulation of ARID1B. However, a subset of ARID1A-dependent genes is not rescued by ARID1B, including many p53 and estrogen receptor (ESR1) targets. Our results provide insight into ARID1A-mediated tumorigenesis and unveil functions of SWI/SNF in modulating RNAPII dynamics.

Long Noncoding RNA Discovery in Cardiovascular Disease: Decoding Form to Function

Despite significant improvements during the past 3 decades, cardiovascular disease remains a leading worldwide health epidemic. The recent identification of a fascinating group of mediators known as long noncoding RNAs (lncRNAs) has provided a wealth of new biology to explore for cardiovascular risk mitigation. lncRNAs are expressed in a highly context-specific fashion, and multiple lines of evidence implicated them in diverse biological processes. Indeed, abnormalities of lncRNAs have been directly linked with human ailments, including cardiovascular biology and disease. Of particular interest to the cardiovascular research community, dysregulation in lncRNA regulatory circuits have been associated with cardiac pathological hypertrophy, vascular disease, cell fate programming and development, atherosclerosis, dyslipidemia, and metabolic syndrome. Although techniques in interrogating noncoding RNAs are rapidly evolving, a major challenge in studying lncRNAs remains navigating through multiple technical constraints. In this review, we provide a road map for lncRNA discovery and interrogation in biological systems relevant to cardiovascular disease and highlight approaches to decipher their modes of action.

Recruitment of BRCA1 limits MYCN-driven accumulation of stalled RNA polymerase

MYC is an oncogenic transcription factor that binds globally to active promoters and promotes transcriptional elongation by RNA polymerase II (RNAPII)^1,2. Deregulated expression of the paralogous protein MYCN drives the development of neuronal and neuroendocrine tumours and is often associated with a particularly poor prognosis³. Here we show that, similar to MYC, activation of MYCN in human neuroblastoma cells induces escape of RNAPII from promoters. If the release of RNAPII from transcriptional pause sites (pause release) fails, MYCN recruits BRCA1 to promoter-proximal regions. Recruitment of BRCA1 prevents MYCN-dependent accumulation of stalled RNAPII and enhances transcriptional activation by MYCN. Mechanistically, BRCA1 stabilizes mRNA decapping complexes and enables MYCN to suppress R-loop formation in promoter-proximal regions. Recruitment of BRCA1 requires the ubiquitin-specific protease USP11, which binds specifically to MYCN when MYCN is dephosphorylated at Thr58. USP11, BRCA1 and MYCN stabilize each other on chromatin, preventing proteasomal turnover of MYCN. Because BRCA1 is highly expressed in neuronal progenitor cells during early development⁴ and MYC is less efficient than MYCN in recruiting BRCA1, our findings indicate that a cell-lineage-specific stress response enables MYCN-driven tumours to cope with deregulated RNAPII function.

Recruitment of BRCA1 limits MYCN-driven accumulation of stalled RNA polymerase

MYC is an oncogenic transcription factor that binds globally to active promoters and promotes transcriptional elongation by RNA polymerase II (RNAPII)^1,2. Deregulated expression of the paralogous protein MYCN drives the development of neuronal and neuroendocrine tumours and is often associated with a particularly poor prognosis³. Here we show that, similar to MYC, activation of MYCN in human neuroblastoma cells induces escape of RNAPII from promoters. If the release of RNAPII from transcriptional pause sites (pause release) fails, MYCN recruits BRCA1 to promoter-proximal regions. Recruitment of BRCA1 prevents MYCN-dependent accumulation of stalled RNAPII and enhances transcriptional activation by MYCN. Mechanistically, BRCA1 stabilizes mRNA decapping complexes and enables MYCN to suppress R-loop formation in promoter-proximal regions. Recruitment of BRCA1 requires the ubiquitin-specific protease USP11, which binds specifically to MYCN when MYCN is dephosphorylated at Thr58. USP11, BRCA1 and MYCN stabilize each other on chromatin, preventing proteasomal turnover of MYCN. Because BRCA1 is highly expressed in neuronal progenitor cells during early development⁴ and MYC is less efficient than MYCN in recruiting BRCA1, our findings indicate that a cell-lineage-specific stress response enables MYCN-driven tumours to cope with deregulated RNAPII function.

Interactions between the transcription and replication machineries regulate the RNA and DNA synthesis in the herpesviruses

The temporal coordination of viral gene expression is imperative for the regulation of the herpesvirus replication cycle. While the main factors of this transcriptional coordination are known, the subtler control mechanisms of gene expression remain elusive. Recent long read sequencing-based approached have revealed an intricate meshwork of overlaps between the herpesvirus transcripts and the overlap of the replication origins with noncoding RNAs. It has been shown that the transcriptional apparatuses can physically interfere with one another while transcribing overlapping regions. We hypothesize that transcriptional interference regulates the global gene expression across the herpesvirus genome. Additionally, an overall decrease in transcriptional activity in individual viral genes has been observed following the onset of DNA replication. An overlap of the replication origins with specific transcripts has also been described in several herpesviruses. The genome-wide interactions between the transcriptional apparatuses and between the replication and transcriptional machineries suggest the existence of novel layers of genetic regulation.

An Overview of Methodologies in Studying lncRNAs in the High-Throughput Era: When Acronyms ATTACK!

The discovery of pervasive transcription in eukaryotic genomes provided one of many surprising (and perhaps most surprising) findings of the genomic era and led to the uncovering of a large number of previously unstudied transcriptional events. This pervasive transcription leads to the production of large numbers of noncoding RNAs (ncRNAs) and thus opened the window to study these diverse, abundant transcripts of unclear relevance and unknown function. Since that discovery, recent advances in high-throughput sequencing technologies have identified a large collection of ncRNAs, from microRNAs to long noncoding RNAs (lncRNAs). Subsequent discoveries have shown that many lncRNAs play important roles in various eukaryotic processes; these discoveries have profoundly altered our understanding of the regulation of eukaryotic gene expression. Although the identification of ncRNAs has become a standard experimental approach, the functional characterization of these diverse ncRNAs remains a major challenge. In this chapter, we highlight recent progress in the methods to identify lncRNAs and the techniques to study the molecular function of these lncRNAs and the application of these techniques to the study of plant lncRNAs.

Saccharomyces cerevisiae Metabolic Labeling with 4-thiouracil and the Quantification of Newly Synthesized mRNA As a Proxy for RNA Polymerase II Activity

Global defects in RNA polymerase II transcription might be overlooked by transcriptomic studies analyzing steady-state RNA. Indeed, the global decrease in mRNA synthesis has been shown to be compensated by a simultaneous decrease in mRNA degradation to restore normal steady-state levels. Hence, the genome-wide quantification of mRNA synthesis, independently from mRNA decay, is the best direct reflection of RNA polymerase II transcriptional activity. Here, we discuss a method using non-perturbing metabolic labeling of nascent RNAs in Saccharomyces cerevisiae (S. cerevisiae). Specifically, the cells are cultured for 6 min with a uracil analog, 4-thiouracil, and the labeled newly transcribed RNAs are purified and quantified to determine the synthesis rates of all individual mRNA. Moreover, using labeled Schizosaccharomyces pombe cells as internal standard allows comparing mRNA synthesis in different S. cerevisiae strains. Using this protocol and fitting the data with a dynamic kinetic model, the corresponding mRNA decay rates can be determined.

Long Non-Coding RNAs Guide the Fine-Tuning of Gene Regulation in B-Cell Development and Malignancy

With the introduction of next generation sequencing methods, such as RNA sequencing, it has become apparent that alterations in the non-coding regions of our genome are important in the development of cancer. Particularly interesting is the class of long non-coding RNAs (lncRNAs), including the recently described subclass of circular RNAs (circRNAs), which display tissue- and cell-type specific expression patterns and exert diverse regulatory functions in the cells. B-cells undergo complex and tightly regulated processes in order to develop from antigen naïve cells residing in the bone marrow to the highly diverse and competent effector cells circulating in peripheral blood. These processes include V(D)J recombination, rapid proliferation, somatic hypermutation and clonal selection, posing a risk of malignant transformation at each step. The aim of this review is to provide insight into how lncRNAs including circRNAs, participate in normal B-cell differentiation, and how deregulation of these molecules is involved in the development of B-cell malignancies. We describe the prognostic value and functional significance of specific deregulated lncRNAs in diseases such as acute lymphoblastic leukemia, chronic lymphocytic leukemia, mantle cell lymphoma, diffuse large B-cell lymphoma, follicular lymphoma, Burkitt lymphoma and multiple myeloma, and we provide an overview of the current knowledge on the role of circRNAs in these diseases.

Integrated omics approaches to characterize a nuclear receptor corepressor-associated histone deacetylase in mouse skeletal muscle

Nuclear receptors regulate gene expression by differentially binding to coactivators or corepressors in a ligand-dependent manner, which further recruits a set of epigenome-modifying enzymes that remodel chromatin conformation. Histone acetylation is a major epigenomic change controlled by histone acetyltransferases (HATs) and histone deacetylases (HDACs). HDAC3 is the only HDAC that confers the enzymatic activity to the complexes nucleated by nuclear receptor corepressors NCoR and SMRT. To address the metabolic function of HDAC3, we have deleted it specifically in mouse skeletal muscles. We have performed the following omics profiling in skeletal muscles of these mice: (1) RNA-seq profiling of total RNA; (2) Global nuclear run-on (GRO-seq) analysis of nascent RNAs; (3) Chromatin immuno-precipitation (ChIP-seq) of HDAC3 at both early evening and early morning; (4) proteomics profiling with mass spectrometry; (5) snap-shot metabolomics profiling of water-soluble metabolites at the basal condition; (6) snap-shot metabolomics profiling of lipid species at the basal condition; (7) kinetic fluxomics analysis of glucose utilization using ¹³C₆-glucose In vivo during treadmill running exercise. These approaches have provided several novel insights into how nuclear receptors regulate circadian rhythm of skeletal muscle fuel metabolism, which has been published elsewhere. Here we present the original datasets and technical details during the execution, analysis, and interpretation of these omics studies.

A critical role for nucleoporin 358 (Nup358) in transposon silencing and piRNA biogenesis in Drosophila

Piwi-interacting RNAs (piRNAs) are a class of small noncoding RNAs that bind Piwi proteins to silence transposons and to regulate gene expression. In Drosophila germ cells, the Aubergine (Aub)-Argonaute 3 (Ago3)-dependent ping-pong cycle generates most germline piRNAs. Loading of antisense piRNAs amplified by this cycle enables Piwi to enter the nucleus and silence transposons. Nuclear localization is crucial for Piwi function in transposon silencing, but how this process is regulated remains unknown. It is also not known whether any of the components of the nuclear pore complex (NPC) directly function in the piRNA pathway. Here, we show that nucleoporin 358 (Nup358) and Piwi interact with each other and that a germline knockdown (GLKD) of Nup358 with short hairpin RNA prevents Piwi entry into the nucleus. The Nup358 GLKD also activated transposons, increased genomic instability, and derailed piRNA biogenesis because of a combination of decreased piRNA precursor transcription and a collapse of the ping-pong cycle. Our results point to a critical role for Nup358 in the piRNA pathway, laying the foundation for future studies to fully elucidate the mechanisms by which Nup358 contributes to piRNA biogenesis and transposon silencing.

DNA methylation at enhancers identifies distinct breast cancer lineages

Breast cancers exhibit genome-wide aberrant DNA methylation patterns. To investigate how these affect the transcriptome and which changes are linked to transformation or progression, we apply genome-wide expression-methylation quantitative trait loci (emQTL) analysis between DNA methylation and gene expression. On a whole genome scale, in cis and in trans, DNA methylation and gene expression have remarkably and reproducibly conserved patterns of association in three breast cancer cohorts (n = 104, n = 253 and n = 277). The expression-methylation quantitative trait loci associations form two main clusters; one relates to tumor infiltrating immune cell signatures and the other to estrogen receptor signaling. In the estrogen related cluster, using ChromHMM segmentation and transcription factor chromatin immunoprecipitation sequencing data, we identify transcriptional networks regulated in a cell lineage-specific manner by DNA methylation at enhancers. These networks are strongly dominated by ERα, FOXA1 or GATA3 and their targets were functionally validated using knockdown by small interfering RNA or GRO-seq analysis after transcriptional stimulation with estrogen.

Emerging roles for long noncoding RNAs in B-cell development and malignancy

Long noncoding (lnc)RNAs have emerged as essential mediators of cellular biology, differentiation and malignant transformation. LncRNAs have a broad range of possible functions at the transcriptional, posttranscriptional and protein level and their aberrant expression significantly contributes to the hallmarks of cancer cell biology. In addition, their high tissue- and cell-type specificity makes lncRNAs especially interesting as biomarkers, prognostic factors or specific therapeutic targets. Here, we review current knowledge on lncRNA expression changes during normal B-cell development, indicating essential functions in the differentiation process. In addition we address lncRNA deregulation in B-cell malignancies, the putative prognostic value of this as well as the molecular functions of multiple deregulated lncRNAs. Altogether, the discussed work indicates major roles for lncRNAs in normal and malignant B cells affecting oncogenic pathways as well as the response to common therapeutics.

Thyroid Hormone Signaling Pathways: Time for a More Precise Nomenclature

Current literature makes a distinction between two pathways for thyroid hormone signaling: genomic and nongenomic. However, this classification is a source of confusion. We propose a clarification in the nomenclature that may help to avoid unproductive controversies and favor progress in this field of research. Four types of thyroid hormone signaling are defined, and the experimental criteria for classification are discussed.