An abundance of bidirectional promoters in the human genome

The evolutionary conservation of the bidirectional activity of the NWC gene promoter in jawed vertebrates and the domestication of the RAG transposon

The RAG-1 and RAG-2 genes form a recombinase complex that is indispensable for V(D)J recombination, which generates the diversity of immunoglobulins and T-cell receptors. It is widely accepted that the presence of RAGs in the genomes of jawed vertebrates and other lineages is a result of the horizontal transfer of a mobile genetic element. While a substantial amount of evidence has been gathered that clarifies the nature of the RAG transposon, far less attention has been paid to the genomic site of its integration in various host organisms. In all genomes of the jawed vertebrates that have been studied to date, the RAG genes are located in close proximity to the NWC gene. We have previously shown that the promoter of the murine NWC genes exhibits a bidirectional activity, which may have facilitated the integration and survival of the RAG transposon in the host genome. In this study, we characterise the promoters of the NWC homologues that are present in the representatives of other jawed vertebrates (H. sapiens, X. tropicalis and D. rerio). We show that the features that are characteristic for promoters as the hosts of a successful transposon integration (in terms of the arrangement, bidirectional and constitutive activity and the involvement of the Zfp143 transcription factor in the promoter regulation) are evolutionarily conserved, which indicates that the presence of RAG genes in jawed vertebrates is a direct result of a successful transposon integration into the NWC locus.

Tissue-specific impact of FADS cluster variants on FADS1 and FADS2 gene expression

Omega-6 (n-6) and omega-3 (n-3) long (≥ 20 carbon) chain polyunsaturated fatty acids (LC-PUFAs) play a critical role in human health and disease. Biosynthesis of LC-PUFAs from dietary 18 carbon PUFAs in tissues such as the liver is highly associated with genetic variation within the fatty acid desaturase (FADS) gene cluster, containing FADS1 and FADS2 that encode the rate-limiting desaturation enzymes in the LC-PUFA biosynthesis pathway. However, the molecular mechanisms by which FADS genetic variants affect LC-PUFA biosynthesis, and in which tissues, are unclear. The current study examined associations between common single nucleotide polymorphisms (SNPs) within the FADS gene cluster and FADS1 and FADS2 gene expression in 44 different human tissues (sample sizes ranging 70-361) from the Genotype-Tissue Expression (GTEx) Project. FADS1 and FADS2 expression were detected in all 44 tissues. Significant cis-eQTLs (within 1 megabase of each gene, False Discovery Rate, FDR<0.05, as defined by GTEx) were identified in 12 tissues for FADS1 gene expression and 23 tissues for FADS2 gene expression. Six tissues had significant (FDR< 0.05) eQTLs associated with both FADS1 and FADS2 (including artery, esophagus, heart, muscle, nerve, and thyroid). Interestingly, the identified eQTLs were consistently found to be associated in opposite directions for FADS1 and FADS2 expression. Taken together, findings from this study suggest common SNPs within the FADS gene cluster impact the transcription of FADS1 and FADS2 in numerous tissues and raise important questions about how the inverse expression of these two genes impact intermediate molecular (such a LC-PUFA and LC-PUFA-containing glycerolipid levels) and ultimately clinical phenotypes associated with inflammatory diseases and brain health.

Identification and functional analysis of long non-coding RNAs in human and mouse early embryos based on single-cell transcriptome data

Epigenetics regulations have an important role in fertilization and proper embryonic development, and several human diseases are associated with epigenetic modification disorders, such as Rett syndrome, Beckwith-Wiedemann syndrome and Angelman syndrome. However, the dynamics and functions of long non-coding RNAs (lncRNAs), one type of epigenetic regulators, in human pre-implantation development have not yet been demonstrated. In this study, a comprehensive analysis of human and mouse early-stage embryonic lncRNAs was performed based on public single-cell RNA sequencing data. Expression profile analysis revealed that lncRNAs are expressed in a developmental stage-specific manner during human early-stage embryonic development, whereas a more temporal-specific expression pattern was identified in mouse embryos. Weighted gene co-expression network analysis suggested that lncRNAs involved in human early-stage embryonic development are associated with several important functions and processes, such as oocyte maturation, zygotic genome activation and mitochondrial functions. We also found that the network of lncRNAs involved in zygotic genome activation was highly preservative between human and mouse embryos, whereas in other stages no strong correlation between human and mouse embryo was observed. This study provides insight into the molecular mechanism underlying lncRNA involvement in human pre-implantation embryonic development.

Genome-wide Analysis of Transcriptional Variability in a Large Maize-Teosinte Population

Gene expression regulation plays an important role in controlling plant phenotypes and adaptation. Here, we report a comprehensive assessment of gene expression variation through the transcriptome analyses of a large maize-teosinte experimental population. Genome-wide mapping identified 25 660 expression quantitative trait loci (eQTL) for 17 311 genes, capturing an unprecedented range of expression variation. We found that local eQTL were more frequently mapped to adjacent genes, displaying a mode of expression piggybacking, which consequently created co-regulated gene clusters. Genes within the co-regulated gene clusters tend to have relevant functions and shared chromatin modifications. Distant eQTL formed 125 significant distant eQTL hotspots with their targets significantly enriched in specific functional categories. By integrating different sources of information, we identified putative trans- regulators for a variety of metabolic pathways. We demonstrated that the bHLH transcription factor R1 and hexokinase HEX9 might act as crucial regulators for flavonoid biosynthesis and glycolysis, respectively. Moreover, we showed that domestication or improvement has significantly affected global gene expression, with many genes targeted by selection. Of particular interest, the Bx genes for benzoxazinoid biosynthesis may have undergone coordinated cis-regulatory divergence between maize and teosinte, and a transposon insertion that inactivates Bx12 was under strong selection as maize spread into temperate environments with a distinct herbivore community.

Transcription start site profiling uncovers divergent transcription and enhancer-associated RNAs in Drosophila melanogaster

Background

High-resolution transcription start site (TSS) mapping in D. melanogaster embryos and cell lines has revealed a rich and detailed landscape of both cis- and trans-regulatory elements and factors. However, TSS profiling has not been investigated in an orthogonal in vivo setting. Here, we present a comprehensive dataset that links TSS dynamics with nucleosome occupancy and gene expression in the wandering third instar larva, a developmental stage characterized by large-scale shifts in transcriptional programs in preparation for metamorphosis.

Results

The data recapitulate major regulatory classes of TSSs, based on peak width, promoter-proximal polymerase pausing, and cis-regulatory element density. We confirm the paucity of divergent transcription units in D. melanogaster, but also identify notable exceptions. Furthermore, we identify thousands of novel initiation events occurring at unannotated TSSs that can be classified into functional categories by their local density of histone modifications. Interestingly, a sub-class of these unannotated TSSs overlaps with functionally validated enhancer elements, consistent with a regulatory role for “enhancer RNAs” (eRNAs) in defining developmental transcription programs.

Conclusions

High-depth TSS mapping is a powerful strategy for identifying and characterizing low-abundance and/or low-stability RNAs. Global analysis of transcription initiation patterns in a developing organism reveals a vast number of novel initiation events that identify potential eRNAs as well as other non-coding transcripts critical for animal development.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4510-7) contains supplementary material, which is available to authorized users.

Intron retention enhances gene regulatory complexity in vertebrates

Background

While intron retention (IR) is now widely accepted as an important mechanism of mammalian gene expression control, it remains the least studied form of alternative splicing. To delineate conserved features of IR, we performed an exhaustive phylogenetic analysis in a highly purified and functionally defined cell type comprising neutrophilic granulocytes from five vertebrate species spanning 430 million years of evolution.

Results

Our RNA-sequencing-based analysis suggests that IR increases gene regulatory complexity, which is indicated by a strong anti-correlation between the number of genes affected by IR and the number of protein-coding genes in the genome of individual species. Our results confirm that IR affects many orthologous or functionally related genes in granulocytes. Further analysis uncovers new and unanticipated conserved characteristics of intron-retaining transcripts. We find that intron-retaining genes are transcriptionally co-regulated from bidirectional promoters. Intron-retaining genes have significantly longer 3′ UTR sequences, with a corresponding increase in microRNA binding sites, some of which include highly conserved sequence motifs. This suggests that intron-retaining genes are highly regulated post-transcriptionally.

Conclusions

Our study provides unique insights concerning the role of IR as a robust and evolutionarily conserved mechanism of gene expression regulation. Our findings enhance our understanding of gene regulatory complexity by adding another contributor to evolutionary adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1339-3) contains supplementary material, which is available to authorized users.

Catsper1 promoter is bidirectional and regulates the expression of a novel lncRNA

The Catsper1 gene, whose expression is restricted to male germ cells, has great importance in reproductive biology because of its function in sperm motility and fertilization. We previously reported that the promoter of this gene has transcriptional activity in either direction in a heterologous system. In the present study, we found that the Catsper1 promoter has in vitro transcriptional activity in either orientation in GC-1 spg mouse spermatogonial cells. The results also showed that this promoter regulates the expression of a new divergent Catsper1 gene named Catsper1au (Catsper1 antisense upstream transcript). Catsper1au is expressed in adult male mouse testis and liver tissues but not in female mouse liver or ovary tissues. In the testis, Catsper1au is expressed in embryos at 11.5 days post-coitum and from newborns to adults. This gene is also expressed in 1- to 3-week postnatal hearts and in 1-week to adult stage livers. The analysis of the 1402 bp whole genome sequence revealed that Catsper1au is an intronless and polyadenylated lncRNA, located in the nuclei of Sertoli and spermatogenic cells from adult testis. These data indicate that Catsper1au is divergently expressed from the Catsper1 promoter and could regulate gene expression during spermatogenesis.

It's easier to get along with the quiet neighbours

Why are neighbouring genes co‐expressed at the RNA level? While it is tempting to think that this is to enable coordination of functionally related proteins, analysis of human proteomics data by Rappsilber and colleagues (Kustatscher et al, 2017) suggests this is the exception not the rule. Rather it might either be just something that happens or, in some instances, an epiphenomenon of coordination of expression to enable reduced gene expression noise.

Chicken CCDC152 shares an NFYB-regulated bidirectional promoter with a growth hormone receptor antisense transcript and inhibits cells proliferation and migration

The chicken coiled-coil domain-containing protein 152 (CCDC152) recently has been identified as a novel one implicated in cell cycle regulation, cellular proliferation and migration by us. Here we demonstrate that CCDC152 is oriented in a head-to-head configuration with the antisense transcript of growth hormone receptor (GHR) gene. Through serial luciferase reporter assays, we firstly identified a minimal 102 bp intergenic region as a core bidirectional promoter to drive basal transcription in divergent orientations. And site mutation and transient transfected assays showed that nuclear transcription factor Y subunit beta (NFYB) could bind to the CCAAT box and directly transactivate this bidirectional promoter. SiRNA-mediated NFYB depletion could significantly down-regulate the expression of both GHR-AS-I6 and CCDC152. Additionally, the expression of GHR-AS-I6 was significantly up-regulated after CCDC152 overexpression. Overexpression of CCDC152 remarkably reduced cell proliferation and migration through JAK2/STAT signaling pathway. Thus, the GHR-AS-I6-CCDC152 bidirectional transcription unit, as a novel direct target of NFYB, is possibly essential for the accelerated proliferation and motility of different cells.

Identification of a novel metastasis inducing lncRNA which suppresses the KAI1/CD82 metastasis suppressor gene and is upregulated in triple-negative breast cancer

Inactivation of tumor/metastasis suppressor genes via epigenetic silencing is a frequent event in human cancers. KAI1/CD82 is a metastasis suppressor gene whose normal protecting activity is deficient in twelve different solid malignancies. Here we have identified and characterized a primarily nuclear non-polyadenylated, antisense (as)-lncRNA, initiating upstream of the KAI1 human metastasis suppressor gene transcription start site; and elongating in the opposite direction to KAI1 mRNA. We show that the KAI1 promoter is bi-directional giving rise to KAI1 mRNA and its as-lncRNA. Moreover, expression of this lncRNA transcript emerges to be inversely related to the KAI1 mRNA expression, and in direct relationship to the invasiveness level of human breast cancer derived cell lines. Importantly, knockdown of the KAI1 as-lncRNA in the triple-negative breast cancer cell line MDA-MB-231 have led to increased KAI1 mRNA and protein expression, manifested in stronger adhesion to fibronectin, retardation of cell migration and reduced cell invasion in vitro. Accordingly we have named this lncRNA, SKAI1BC, standing for "Suppressor of KAI1 in Breast Cancer". These results uncover a potential way to harness tumor metastasis via targeting SKAI1BC in human breast cancer, and perhaps also in other KAI1-deficient human malignancies.

Analysis of the relationship between coexpression domains and chromatin 3D organization

Gene order is not random in eukaryotic chromosomes, and co-regulated genes tend to be clustered. The mechanisms that determine co-regulation of large regions of the genome and its connection with chromatin three-dimensional (3D) organization are still unclear however. Here we have adapted a recently described method for identifying chromatin topologically associating domains (TADs) to identify coexpression domains (which we term “CODs”). Using human normal breast and breast cancer RNA-seq data, we have identified approximately 500 CODs. CODs in the normal and breast cancer genomes share similar characteristics but differ in their gene composition. COD genes have a greater tendency to be coexpressed with genes that reside in other CODs than with non-COD genes. Such inter-COD coexpression is maintained over large chromosomal distances in the normal genome but is partially lost in the cancer genome. Analyzing the relationship between CODs and chromatin 3D organization using Hi-C contact data, we find that CODs do not correspond to TADs. In fact, intra-TAD gene coexpression is the same as random for most chromosomes. However, the contact profile is similar between gene pairs that reside either in the same COD or in coexpressed CODs. These data indicate that co-regulated genes in the genome present similar patterns of contacts irrespective of the frequency of physical chromatin contacts between them.

Prokaryotic operons normally comprise functionally related genes whose expression is coordinated. Even though operons do not exist in most eukaryotes, results from the last fifteen years indicate that gene order is nonetheless not random in eukaryotes, and that coexpressed genes tend to be grouped in the genome. We identify here about 500 coexpression domain (CODs) in normal breast tissue. Interestingly, we find that genes within CODs often are coexpressed with other genes that reside in other CODs placed very far away in the same chromosome, which is indicative of long-range inter-COD co-regulation. Furthermore, we find that coexpressed genes within CODs or within co-regulated CODs display similar three-dimensional chromatin contacts, suggesting a spatial coordination of CODs.

RNA biology of disease-associated microsatellite repeat expansions

Microsatellites, or simple tandem repeat sequences, occur naturally in the human genome and have important roles in genome evolution and function. However, the expansion of microsatellites is associated with over two dozen neurological diseases. A common denominator among the majority of these disorders is the expression of expanded tandem repeat-containing RNA, referred to as xtrRNA in this review, which can mediate molecular disease pathology in multiple ways. This review focuses on the potential impact that simple tandem repeat expansions can have on the biology and metabolism of RNA that contain them and underscores important gaps in understanding. Merging the molecular biology of repeat expansion disorders with the current understanding of RNA biology, including splicing, transcription, transport, turnover and translation, will help clarify mechanisms of disease and improve therapeutic development.

Pervasive coexpression of spatially proximal genes is buffered at the protein level

Genes are not randomly distributed in the genome. In humans, 10% of protein-coding genes are transcribed from bidirectional promoters and many more are organised in larger clusters. Intriguingly, neighbouring genes are frequently coexpressed but rarely functionally related. Here we show that coexpression of bidirectional gene pairs, and closeby genes in general, is buffered at the protein level. Taking into account the 3D architecture of the genome, we find that co-regulation of spatially close, functionally unrelated genes is pervasive at the transcriptome level, but does not extend to the proteome. We present evidence that non-functional mRNA coexpression in human cells arises from stochastic chromatin fluctuations and direct regulatory interference between spatially close genes. Protein-level buffering likely reflects a lack of coordination of post-transcriptional regulation of functionally unrelated genes. Grouping human genes together along the genome sequence, or through long-range chromosome folding, is associated with reduced expression noise. Our results support the hypothesis that the selection for noise reduction is a major driver of the evolution of genome organisation.

Genome-wide identification and expression profiling of long non-coding RNAs in auditory and vestibular systems

Mammalian genomes encode multiple layers of regulation, including a class of RNA molecules known as long non-coding RNAs (lncRNAs). These are >200 nucleotides in length and similar to mRNAs, they are capped, polyadenylated, and spliced. In contrast to mRNAs, lncRNAs are less abundant and have higher tissue specificity, and have been linked to development, epigenetic processes, and disease. However, little is known about lncRNA function in the auditory and vestibular systems, or how they play a role in deafness and vestibular dysfunction. To help address this need, we performed a whole-genome identification of lncRNAs using RNA-seq at two developmental stages of the mouse inner ear sensory epithelium of the cochlea and vestibule. We identified 3,239 lncRNA genes, most of which were intergenic (lincRNAs) and 721 are novel. We examined temporal and tissue specificity by analyzing the developmental profiles on embryonic day 16.5 and at birth. The spatial and temporal patterns of three lncRNAs, two of which are in proximity to genes associated with hearing and deafness, were explored further. Our findings indicate that lncRNAs are prevalent in the sensory epithelium of the mouse inner ear and are likely to play key roles in regulating critical pathways for hearing and balance.

Comprehensive bioinformatics analysis of the characterization and determination underlying mechanisms of over-expression and co-expression of genes residing on 20q in colorectal cancer

The Long arm of chromosome 20 (20q) is closely related to the development of colorectal cancer, so identifying the expression profile of genes on 20q through a comprehensive overview is indispensable. In this article, preliminar experimental data, several available databases and bioinformatics tools such as the Cancer Genome Atlas, the Encyclopedia of DNA Elements, the JASPAR database and starBase were combined to analyze the correlation between genes and chromosomal aberrations, microRNA and transcription factors, as well as to explore the expression feature and potential regulative mechanism. The results showed that the most frequently unregulated genes in colorectal cancer arelocated on chromosome 20q, present a significant CNA–mRNA correlation.Furthermore, the genes with mRNA overexpression showed co-expression features and tended to be clustered within the same genomic neighborhoods. Then, several genes were selected to carry out further analysis and demonstrated that shared transcription factors, a conserved bidirectional promoter, and competition for a limited pool of microRNAin the 3’UTR of mRNA may be the underlying mechanisms behind the co-expression of physically adjacent genes.Finally, the databases, Lentivirus shRNA, and qPCR were used to find that these adjacent genes with co-expression cooperatively participated in the same biological pathways associated with the pathogenesis and development of colorectal cancer.

Sound of silence: the beauvericin cluster in Fusarium fujikuroi is controlled by cluster-specific and global regulators mediated by H3K27 modification

In this study, we compared the secondary metabolite profile of Fusarium fujikuroi and the histone deacetylase mutant ΔHDA1. We identified a novel peak in ΔHDA1, which was identified as beauvericin (BEA). Going in line with a 1000-fold increased BEA production, the respective non-ribosomal peptide synthetase (NRPS)-encoding gene (BEA1), as well as two adjacent genes (BEA2-BEA3), were significantly up-regulated in ΔHDA1 compared to the wild type. A special role was revealed for the ABC transporter Bea3: deletion of the encoding gene resulted in significant up-regulation of BEA1 and BEA2 and drastically elevated product yields. Furthermore, mutation of a conserved sequence motif in the promoter of BEA1 released BEA repression and resulted in elevated product levels. Candidate transcription factors (TFs) that could bind to this motif are the cluster-specific TF Bea4 as well as a homolog of the global mammalian Kruppel-like TF Yin Yang 1 (Yy1), both acting as repressors of BEA biosynthesis. In addition to Hda1, BEA biosynthesis is repressed by the activity of the H3K27 methyltransferase Kmt6. Consistently, Western blot analyses revealed a genome-wide enrichment of H3K27 acetylation (H3K27ac) in the ΔHDA1 and KMT6 knock-down mutants. Subsequent chromatin immunoprecipitation (ChIP) experiments showed elevated H3K27ac modification levels at the BEA cluster.

RNA Polymerase Collision versus DNA Structural Distortion: Twists and Turns Can Cause Break Failure

The twisting of DNA due to the movement of RNA polymerases is the basis of numerous classic experiments in molecular biology. Recent mouse genetic models indicate that chromosomal breakage is common at sites of transcriptional turbulence. Two key studies on this point mapped breakpoints to sites of either convergent or divergent transcription but arrived at different conclusions as to which is more detrimental and why. The issue hinges on whether DNA strand separation is the basis for the chromosomal instability or collision of RNA polymerases.

High resolution temporal transcriptomics of mouse embryoid body development reveals complex expression dynamics of coding and noncoding loci

Cellular responses to stimuli are rapid and continuous and yet the vast majority of investigations of transcriptional responses during developmental transitions typically use long interval time courses; limiting the available interpretive power. Moreover, such experiments typically focus on protein-coding transcripts, ignoring the important impact of long noncoding RNAs. We therefore evaluated coding and noncoding expression dynamics at unprecedented temporal resolution (6-hourly) in differentiating mouse embryonic stem cells and report new insight into molecular processes and genome organization. We present a highly resolved differentiation cascade that exhibits coding and noncoding transcriptional alterations, transcription factor network interactions and alternative splicing events, little of which can be resolved by long-interval developmental time-courses. We describe novel short lived and cycling patterns of gene expression and dissect temporally ordered gene expression changes in response to transcription factors. We elucidate patterns in gene co-expression across the genome, describe asynchronous transcription at bidirectional promoters and functionally annotate known and novel regulatory lncRNAs. These findings highlight the complex and dynamic molecular events underlying mammalian differentiation that can only be observed though a temporally resolved time course.

From intronization to intron loss: How the interplay between mRNA-associated processes can shape the architecture and the expression of eukaryotic genes

Transcription-coupled processes such as capping, splicing, and cleavage/polyadenylation participate in the journey from genes to proteins. Although they are traditionally thought to serve only as steps in the generation of mature mRNAs, a synthesis of available data indicates that these processes could also act as a driving force for the evolution of eukaryotic genes. A theoretical framework for how mRNA-associated processes may shape gene structure and expression has recently been proposed. Factors that promote splicing and cleavage/polyadenylation in this framework compete for access to overlapping or neighboring signals throughout the transcription cycle. These antagonistic interactions allow mechanisms for intron gain and splice site recognition as well as common trends in eukaryotic gene structure and expression to be coherently integrated. Here, I extend this framework further. Observations that largely (but not exclusively) revolve around the formation of DNA-RNA hybrid structures, called R loops, and promoter directionality are integrated. Additionally, the interplay between splicing factors and cleavage/polyadenylation factors is theorized to also affect the formation of intragenic DNA double-stranded breaks thereby contributing to intron loss. The most notable prediction in this proposition is that RNA molecules can mediate intron loss by serving as a template to repair DNA double-stranded breaks. The framework presented here leverages a vast body of empirical observations, logically extending previous suggestions, and generating verifiable predictions to further substantiate the view that the intracellular environment plays an active role in shaping the structure and the expression of eukaryotic genes.

Analysis of Promoters of Arabidopsis thaliana Divergent Gene Pair SERAT3;2 and IDH-III Shows SERAT3;2 Promoter is Nested Within the IDH-III Promoter

Intergenic regions of divergent gene pairs show bidirectional promoter activity but whether regulatory sequences for gene expression in opposite directions are shared is not established. In this study, promoters of divergently arranged gene pair At4g35640-At4g35650 (SERAT3;2-IDH-III) of Arabidopsis thaliana were analyzed to identify overlapping regulatory regions. Both genes showed the highest expression in flower buds and flowers. 5' RACE experiments extended the intergenic region from 161 bp shown in TAIR annotation to 512 bp. GUS analysis of transgenic A. thaliana plants carrying the 691 bp fragment (512 bp intergenic region plus 5' UTR of both the genes) linked to uidA gene revealed that SERAT3;2 promoter drives gene expression in the tapetum, whereas IDH-III promoter functions specifically in microspores/pollen. Serial 5' deletion of the 691 bp fragment showed SERAT3;2 promoter extends up to -355 position, whereas IDH-III promoter encompasses the 512 bp intergenic region. In transgenics, uidA transcript levels were lower than native SERAT3;2 and IDH-III transcripts indicating presence of additional cis regulatory elements beyond the 691 bp fragment. The present study demonstrated for the first time occurrence of a nested promoter in plants and identified a novel bidirectional promoter capable of driving gene expression in tapetum and microspores/pollen.

Long non-coding RNA expression profile in cervical cancer tissues

Cervical cancer (CC), one of the most common types of cancer of the female population, presents an enormous challenge in diagnosis and treatment. Long non-coding (lnc)RNAs, non-coding (nc)RNAs with length >200 nucleotides, have been identified to be associated with multiple types of cancer, including CC. This class of nc transcripts serves an important role in tumor suppression and oncogenic signaling pathways. In the present study, the microarray method was used to obtain the expression profile of lncRNAs and protein-coding mRNAs and to compare the expression of lncRNAs between CC tissues and corresponding adjacent non-cancerous tissues in order to screen potential lncRNAs for associations with CC. Overall, 3356 lncRNAs with significantly different expression pattern in CC tissues compared with adjacent non-cancerous tissues were identified, while 1,857 of them were upregulated. These differentially expressed lncRNAs were additionally classified into 5 subgroups. Reverse transcription quantitative polymerase chain reactions were performed to validate the expression pattern of 5 random selected lncRNAs, and 2lncRNAs were identified to have significantly different expression in CC samples compared with adjacent non-cancerous tissues. This finding suggests that those lncRNAs with different expression may serve important roles in the development of CC, and the expression data may provide information for additional study on the involvement of lncRNAs in CC.

Genetic variability of the activity of bidirectional promoters: a pilot study in bovine muscle

Bidirectional promoters are regulatory regions co-regulating the expression of two neighbouring genes organized in a head-to-head orientation. In recent years, these regulatory regions have been studied in many organisms; however, no investigation to date has been done to analyse the genetic variation of the activity of this type of promoter regions. In our study, we conducted an investigation to first identify bidirectional promoters sharing genes expressed in bovine Longissimus thoracis and then to find genetic variants affecting the activity of some of these bidirectional promoters. Combining bovine gene information and expression data obtained using RNA-Seq, we identified 120 putative bidirectional promoters active in bovine muscle. We experimentally validated in vitro 16 of these bidirectional promoters. Finally, using gene expression and whole-genome genotyping data, we explored the variability of the activity in muscle of the identified bidirectional promoters and discovered genetic variants affecting their activity. We found that the expression level of 77 genes is correlated with the activity of 12 bidirectional promoters. We also identified 57 single nucleotide polymorphisms associated with the activity of 5 bidirectional promoters. To our knowledge, our study is the first analysis in any species of the genetic variability of the activity of bidirectional promoters.

Developmental transitions in Arabidopsis are regulated by antisense RNAs resulting from bidirectionally transcribed genes

Transcription terminators are DNA elements located at the 3' end of genes that ensure efficient cleavage of nascent RNA generating the 3' end of mRNA, as well as facilitating disengagement of elongating DNA-dependent RNA polymerase II. Surprisingly, terminators are also a potent source of antisense transcription. We have recently described an Arabidopsis antisense transcript originating from the 3' end of a master regulator of Arabidopsis thaliana seed dormancy DOG1. In this review, we discuss the broader implications of our discovery in light of recent developments in yeast and Arabidopsis. We show that, surprisingly, the key features of terminators that give rise to antisense transcription are preserved between Arabidopsis and yeast, suggesting a conserved mechanism. We also compare our discovery to known antisense-based regulatory mechanisms, highlighting the link between antisense-based gene expression regulation and major developmental transitions in plants.

Evolutionary acquisition of promoter-associated non-coding RNA (pancRNA) repertoires diversifies species-dependent gene activation mechanisms in mammals

Background

Recent transcriptome analyses have shown that long non-coding RNAs (ncRNAs) play extensive roles in transcriptional regulation. In particular, we have reported that promoter-associated ncRNAs (pancRNAs) activate the partner gene expression via local epigenetic changes.

Results

Here, we identify thousands of genes under pancRNA-mediated transcriptional activation in five mammalian species in common. In the mouse, 1) pancRNA-partnered genes confined their expression pattern to certain tissues compared to pancRNA-lacking genes, 2) expression of pancRNAs was significantly correlated with the enrichment of active chromatin marks, H3K4 trimethylation and H3K27 acetylation, at the promoter regions of the partner genes, 3) H3K4me1 marked the pancRNA-partnered genes regardless of their expression level, and 4) C- or G-skewed motifs were exclusively overrepresented between−200 and−1 bp relative to the transcription start sites of the pancRNA-partnered genes. More importantly, the comparative transcriptome analysis among five different mammalian species using a total of 25 counterpart tissues showed that the overall pancRNA expression profile exhibited extremely high species-specificity compared to that of total mRNA, suggesting that interspecies difference in pancRNA repertoires might lead to the diversification of mRNA expression profiles.

Conclusions

The present study raises the interesting possibility that the gain and/or loss of gene-activation-associated pancRNA repertoires, caused by formation or disruption of the genomic GC-skewed structure in the course of evolution, finely shape the tissue-specific pattern of gene expression according to a given species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3662-1) contains supplementary material, which is available to authorized users.

The Evolution and Expression Pattern of Human Overlapping lncRNA and Protein-coding Gene Pairs

Long non-coding RNA overlapping with protein-coding gene (lncRNA-coding pair) is a special type of overlapping genes. Protein-coding overlapping genes have been well studied and increasing attention has been paid to lncRNAs. By studying lncRNA-coding pairs in human genome, we showed that lncRNA-coding pairs were more likely to be generated by overprinting and retaining genes in lncRNA-coding pairs were given higher priority than non-overlapping genes. Besides, the preference of overlapping configurations preserved during evolution was based on the origin of lncRNA-coding pairs. Further investigations showed that lncRNAs promoting the splicing of their embedded protein-coding partners was a unilateral interaction, but the existence of overlapping partners improving the gene expression was bidirectional and the effect was decreased with the increased evolutionary age of genes. Additionally, the expression of lncRNA-coding pairs showed an overall positive correlation and the expression correlation was associated with their overlapping configurations, local genomic environment and evolutionary age of genes. Comparison of the expression correlation of lncRNA-coding pairs between normal and cancer samples found that the lineage-specific pairs including old protein-coding genes may play an important role in tumorigenesis. This work presents a systematically comprehensive understanding of the evolution and the expression pattern of human lncRNA-coding pairs.

The PRR11-SKA2 Bidirectional Transcription Unit Is Negatively Regulated by p53 through NF-Y in Lung Cancer Cells

We previously identified proline-rich protein 11 (PRR11) as a novel cancer-related gene that is implicated in the regulation of cell cycle and tumorigenesis. Our recent study demonstrated that PRR11 and its adjacent gene, kinetochore associated 2 (SKA2), constitute a classic head-to-head gene pair that is coordinately regulated by nuclear factor Y (NF-Y). In the present study, we further show that the PRR11-SKA2 bidirectional transcription unit is an indirect target of the tumor suppressor p53. A luciferase reporter assay revealed that overexpression of wild type p53, but not mutant p53, significantly represses the basal activity and NF-Y mediated transactivation of the PRR11-SKA2 bidirectional promoter. Deletion and mutation analysis of the PRR11-SKA2 promoter revealed that p53-mediated PRR11-SKA2 repression is dependent on the presence of functional NF-Y binding sites. Furthermore, a co-immunoprecipitation assay revealed that p53 associates with NF-Y in lung cancer cells, and a chromatin immunoprecipitation assay showed that p53 represses PRR11-SKA2 transcription by reducing the binding amount of NF-Y in the PRR11-SKA2 promoter region. Consistently, the ability of p53 to downregulate PRR11-SKA2 transcription was significantly attenuated upon siRNA-mediated depletion of nuclear factor Y subunit beta (NF-YB). Notably, lung cancer patients with lower expression of either PRR11 or SKA2 along with wild type p53 exhibited the best overall survival compared with others with p53 mutation and/or higher expression of either PRR11 or SKA2. Taken together, our results demonstrate that p53 negatively regulates the expression of the PRR11-SKA2 bidirectional transcription unit through NF-Y, suggesting that the inability to repress the PRR11-SKA2 bidirectional transcription unit after loss of p53 might contribute to tumorigenesis.

Long noncoding RNA repertoire in chicken liver and adipose tissue

BACKGROUND

Improving functional annotation of the chicken genome is a key challenge in bridging the gap between genotype and phenotype. Among all transcribed regions, long noncoding RNAs (lncRNAs) are a major component of the transcriptome and its regulation, and whole-transcriptome sequencing (RNA-Seq) has greatly improved their identification and characterization. We performed an extensive profiling of the lncRNA transcriptome in the chicken liver and adipose tissue by RNA-Seq. We focused on these two tissues because of their importance in various economical traits for which energy storage and mobilization play key roles and also because of their high cell homogeneity. To predict lncRNAs, we used a recently developed tool called FEELnc, which also classifies them with respect to their distance and strand orientation to the closest protein-coding genes. Moreover, to confidently identify the genes/transcripts expressed in each tissue (a complex task for weakly expressed molecules such as lncRNAs), we probed a particularly large number of biological replicates (16 per tissue) compared to common multi-tissue studies with a larger set of tissues but less sampling.

RESULTS

We predicted 2193 lncRNA genes, among which 1670 were robustly expressed across replicates in the liver and/or adipose tissue and which were classified into 1493 intergenic and 177 intragenic lncRNAs located between and within protein-coding genes, respectively. We observed similar structural features between chickens and mammals, with strong synteny conservation but without sequence conservation. As previously reported, we confirm that lncRNAs have a lower and more tissue-specific expression than mRNAs. Finally, we showed that adjacent lncRNA-mRNA genes in divergent orientation have a higher co-expression level when separated by less than 1 kb compared to more distant divergent pairs. Among these, we highlighted for the first time a novel lncRNA candidate involved in lipid metabolism, lnc_DHCR24, which is highly correlated with the DHCR24 gene that encodes a key enzyme of cholesterol biosynthesis.

CONCLUSIONS

We provide a comprehensive lncRNA repertoire in the chicken liver and adipose tissue, which shows interesting patterns of co-expression between mRNAs and lncRNAs. It contributes to improving the structural and functional annotation of the chicken genome and provides a basis for further studies on energy storage and mobilization traits in the chicken.

MYH9 binds to lncRNA gene PTCSC2 and regulates FOXE1 in the 9q22 thyroid cancer risk locus

A locus on chromosome 9q22 harbors a SNP (rs965513) firmly associated with risk of papillary thyroid carcinoma (PTC). The locus also comprises the forkhead box E1 (FOXE1) gene, which is implicated in thyroid development, and a long noncoding RNA (lncRNA) gene, papillary thyroid cancer susceptibility candidate 2 (PTCSC2). How these might interact is not known. Here we report that PTCSC2 binds myosin-9 (MYH9). In a bidirectional promoter shared by FOXE1 and PTCSC2, MYH9 inhibits the promoter activity in both directions. This inhibition can be reversed by PTCSC2, which acts as a suppressor. RNA knockdown of FOXE1 in primary thyroid cells profoundly interferes with the p53 pathway. We propose that the interaction between the lncRNA, its binding protein MYH9, and the coding gene FOXE1 underlies the predisposition to PTC triggered by rs965513.

A light-switchable bidirectional expression system in filamentous fungus Trichoderma reesei

The filamentous fungi Trichoderma reesei is widely used in the production of cellulolytic enzymes and recombinant proteins. However, only moderate success has been achieved in expressing heterologous proteins in T. reesei. Light-dependent control of DNA transcription, and protein expression have been demonstrated in bacteria, fungi, and mammalian cells. In this study, light inducible transactivators, a "light-on" bidirectional promoter and a "light-off" promoter were constructed successfully in T. reesei for the first time. Our light inducible transactivators can homodimerize and bind to the upstream region of artificial promoters to activate or repress genes transcription. Additionally, we upgraded the light-inducible system to on-off system that can simultaneously control the expression of multiple heterologous proteins in T. reesei. Moreover, a native cellulase-free background for the expression of heterologous proteins was achieved by knocking out the genes involved in transcriptional regulation and encoding of cellulases: xyr1, cbh1, and cbh2. Our light-switchable system showed a very little background protein expression and robust activation in the blue light with significantly improved heterologous protein expression. We demonstrate that our light-switchable system has a potential application as an on/off "switch" that can simultaneously regulate the expression of multiple genes in T. reesei under native cellulase-free background.

Transcriptional Interference Promotes Rapid Expression Divergence of Drosophila Nested Genes

Nested genes are the most common form of protein-coding overlap in eukaryotic genomes. Previous studies have shown that nested genes accumulate rapidly over evolutionary time, typically via the insertion of short young duplicate genes into long introns. However, the evolutionary relationship between nested genes remains unclear. Here, I compare RNA-seq expression profiles of nested, proximal intra-chromosomal, intermediate intra-chromosomal, distant intra-chromosomal, and inter-chromosomal gene pairs in two Drosophila species. I find that expression profiles of nested genes are more divergent than those of any other class of genes, supporting the hypothesis that concurrent expression of nested genes is deleterious due to transcriptional interference. Further analysis reveals that expression profiles of derived nested genes are more divergent than those of their ancestral un-nested orthologs, which are more divergent than those of un-nested genes with similar genomic features. Thus, gene expression divergence between nested genes is likely caused by selection against nesting of genes with insufficiently divergent expression profiles, as well as by continued expression divergence after nesting. Moreover, expression divergence and sequence evolutionary rates are elevated in young nested genes and reduced in old nested genes, indicating that a burst of rapid evolution occurs after nesting. Together, these findings suggest that similarity between expression profiles of nested genes is deleterious due to transcriptional interference, and that natural selection addresses this problem both by eradicating highly deleterious nestings and by enabling rapid expression divergence of surviving nested genes, thereby quickly limiting or abolishing transcriptional interference.

Histone modifications facilitate the coexpression of bidirectional promoters in rice

Background

Bidirectional gene pairs are highly abundant and mostly co-regulated in eukaryotic genomes. The structural features of bidirectional promoters (BDPs) have been well studied in yeast, humans and plants. However, the underlying mechanisms responsible for the coexpression of BDPs remain understudied, especially in plants.

Results

Here, we characterized chromatin features associated with rice BDPs. Several unique chromatin features were present in rice BDPs but were missing from unidirectional promoters (UDPs), including overrepresented active histone marks, canonical nucleosomes and underrepresented H3K27me3. In particular, overrepresented active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me2 and H3K36me3) were truly overrepresented in type I BDPs but not in the other two BDPs, based on a Kolmogorov-Smirnov test.

Conclusions

Our analyses indicate that active marks (H3K4ac, H4K12ac, H4K16ac, H3K4me3, H3K9ac and H3K27ac) may coordinate with repressive marks (H3K27me3 and H3K9me1/3) to build a unique chromatin structure that favors the coregulation of bidirectional gene pairs. Thus, our findings help to enhance the understanding of unique epigenetic mechanisms that regulate bidirectional gene pairs and may improve the manipulation of gene pairs for crop bioengineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3125-0) contains supplementary material, which is available to authorized users.

Avian Leukosis Virus Activation of an Antisense RNA Upstream of TERT in B-Cell Lymphomas

Avian leukosis virus (ALV) induces tumors by integrating its proviral DNA into the chicken genome and altering the expression of nearby genes via strong promoter and enhancer elements. Viral integration sites that contribute to oncogenesis are selected in tumor cells. Deep-sequencing analysis of B-cell lymphoma DNA confirmed that the telomerase reverse transcriptase (TERT) gene promoter is a common ALV integration target. Twenty-six unique proviral integration sites were mapped between 46 and 3,552 nucleotides (nt) upstream of the TERT transcription start site, predominantly in the opposite transcriptional orientation to TERT Transcriptome-sequencing (RNA-seq) analysis of normal bursa revealed a transcribed region upstream of TERT in the opposite orientation, suggesting the TERT promoter is bidirectional. This transcript appears to be an uncharacterized antisense RNA. We have previously shown that TERT expression is upregulated in tumors with integrations in the TERT promoter region. We now report that the viral promoter drives the expression of a chimeric transcript containing viral sequences spliced to exons 4 through 7 of this antisense RNA. Clonal expansion of cells with ALV integrations driving overexpression of the TERT antisense RNA suggest it may have a role in tumorigenesis.

IMPORTANCE

The data suggest that ALV integrations in the TERT promoter region drive the overexpression of a novel antisense RNA and contribute to the development of lymphomas.

Involvement of FOXL2 and RSPO1 in Ovarian Determination, Development, and Maintenance in Mammals

In mammals, sex determination is a process through which the gonad is committed to differentiate into a testis or an ovary. This process relies on a delicate balance between genetic pathways that promote one fate and inhibit the other. Once the gonad is committed to the female pathway, ovarian differentiation begins and, depending on the species, is completed during gestation or shortly after birth. During this step, granulosa cell precursors, steroidogenic cells, and primordial germ cells start to express female-specific markers in a sex-dimorphic manner. The germ cells then arrest at prophase I of meiosis and, together with somatic cells, assemble into functional structures. This organization gives the ovary its definitive morphology and functionality during folliculogenesis. Until now, 2 main genetic cascades have been shown to be involved in female sex differentiation. The first is driven by FOXL2, a transcription factor that also plays a crucial role in folliculogenesis and ovarian fate maintenance in adults. The other operates through the WNT/CTNNB1 canonical pathway and is regulated primarily by R-spondin1. Here, we discuss the roles of FOXL2 and RSPO1/WNT/ CTNNB1 during ovarian development and homeostasis in different models, such as humans, goats, and rodents.

A potent bidirectional promoter from the monocot cereal Eleusine coracana

Ragi bifunctional α-amylase-trypsin inhibitor (RBI) of Eleusine coracana (L.) Gaertn. (finger millet) simultaneously inhibits α-amylase and trypsin. In continuation of previous work on the cloning, expression and characterization of RBI, a bidirectional promoter from finger millet was explored on the basis of experimental observations. Two trypsin inhibitors were identified while purifying RBI from a trypsin-Sepharose column eluent. Using an FPLC gel filtration column, these three inhibitors were purified to homogeneity and subjected to MALDI-TOF-TOF-MS/MS analysis and N-terminal sequencing. Both ragi trypsin inhibitors (RTIs) showed the same N-terminal sequence and considerable sequence similarity to RBI, indicating the presence of a multigene protease inhibitor family in finger millet. To gain insight into the evolution of these genes, the upstream region of RBI was explored by Genome Walking. Interestingly, on sequencing, a genome walking product of ∼1 Kb showed presence of an N-terminal RBI specific primer sequence twice but in opposite directions and leaving an intervening region of ∼0.9 Kb. The intervening region was presumed to represent an E. coracana bidirectional promoter (EcBDP), intuitively having a divergent RBI-RTI gene pair at two sides. For assaying the bidirectionality of promoter activity, a dual reporter GUS-GFP vector construct was made for plant expression containing the reporter genes at two ends of EcBDP, which was used to transform Agrobacterium tumefaciens LBA 4404. Transient plant transformation by recombinant Agrobacterium cells was carried out in onion scale epidermal cells and finger millet seedling leaves. Simultaneous expression of GUS and GFP under EcBDP established it as a potent natural bidirectional promoter from monocot origin, thereby potentially having vast application in cereal gene manipulations. In addition, inducibility of the EcBDP by either abscisic acid or cold treatment, as determined by transient transformation in onion, would substantiate more precise control of gene expression to mitigate the effects of adverse environmental conditions.

Analysis of Ly49 gene transcripts in mature NK cells supports a role for the Pro1 element in gene activation, not gene expression

The variegated expression of murine Ly49 loci has been associated with the probabilistic behavior of an upstream promoter active in immature cells, the Pro1 element. However, recent data suggest that Pro1 may be active in mature natural killer (NK) cells and function as an enhancer element. To assess directly if Pro1 transcripts are present in mature Ly49-expressing NK cells, RNA-sequencing of the total transcript pool was performed on freshly isolated splenic NK cells sorted for expression of either Ly49G or Ly49I. No Pro1 transcripts were detected from the Ly49a, Ly49c or Ly49i genes in mature Ly49(+) NK cells that contained high levels of Pro2 transcripts. Low levels of Ly49g Pro1 transcripts were found in both Ly49G(+) and Ly49G(-) populations, consistent with the presence of a small population of mature NK cells undergoing Ly49g gene activation, as previously demonstrated by culture of splenic NK cells in interleukin-2. Ly49 gene reporter constructs containing Pro1 failed to show any enhancer activity of Pro1 on Pro2 in a mature Ly49-expressing cell line. Taken together, the results are consistent with Pro1 transcription having a role in gene activation in developing NK, and argue against a role for Pro1 in Ly49 gene transcription by mature NK cells.

Functional characterization of open chromatin in bidirectional promoters of rice

Bidirectional gene pairs tend to be highly coregulated and function in similar biological processes in eukaryotic genomes. Structural features and functional consequences of bidirectional promoters (BDPs) have received considerable attention among diverse species. However, the underlying mechanisms responsible for the bidirectional transcription and coexpression of BDPs remain poorly understood in plants. In this study, we integrated DNase-seq, RNA-seq, ChIP-seq and MNase-seq data and investigated the effect of physical DNase I hypersensitive site (DHS) positions on the transcription of rice BDPs. We found that the physical position of a DHS relative to the TSS of bidirectional gene pairs can affect the expression of the corresponding genes: the closer a DHS is to the TSS, the higher is the expression level of the genes. Most importantly, we observed that the distribution of DHSs plays a significant role in the regulation of transcription and the coexpression of gene pairs, which are possibly mediated by orchestrating the positioning of histone marks and canonical nucleosomes around BDPs. Our results demonstrate that the combined actions of chromatin structures with DHSs, which contain functional cis-elements for interaction with transcriptional machinery, may play an important role in the regulation of the bidirectional transcription or coexpression in rice BDPs. Our findings may help to enhance the understanding of DHSs in the regulation of bidirectional gene pairs.

A low-temperature-responsive element involved in the regulation of the Arabidopsis thaliana At1g71850/At1g71860 divergent gene pair

The bidirectional promoter of the Arabidopsis thaliana gene pair At1g71850/At1g71860 harbors low-temperature-responsive elements, which participate in anti-correlated transcription regulation of the driving genes in response to environmental low temperature. A divergent gene pair is defined as two adjacent genes organized head to head in opposite orientation, sharing a common promoter region. Divergent gene pairs are mainly coexpressed, but some display opposite regulation. The mechanistic basis of such anti-correlated regulation is not well understood. Here, the regulation of the Arabidopsis thaliana gene pair At1g71850/At1g71860 was investigated. Semi-quantitative RT-PCR and Genevestigator analyses showed that while one of the pair was upregulated by exposure to low temperature, the same treatment downregulated the other. Promoter::GUS fusion transgenes were used to show that this behavior was driven by a bidirectional promoter, which harbored an as-1 motif, associated with the low-temperature response; mutation of this sequence produced a significant decrease in cold-responsive expression. With regard to the as-1 motif in the native orientation repressing the promoter's low-temperature responsiveness, the same as-1 motif introduced in the reverse direction showed a slight enhancement in the promoter's responsiveness to low-temperature exposure, indicating that the orientation of the motif was important for the promoter's activity. These findings provide new insights into the complex transcriptional regulation of bidirectional gene pairs as well as plant stress response.

Principles for RNA metabolism and alternative transcription initiation within closely spaced promoters

Mammalian transcriptomes are complex and formed by extensive promoter activity. In addition, gene promoters are largely divergent and initiate transcription of reverse-oriented promoter upstream transcripts (PROMPTs). Although PROMPTs are commonly terminated early, influenced by polyadenylation sites, promoters often cluster so that the divergent activity of one might impact another. Here we found that the distance between promoters strongly correlates with the expression, stability and length of their associated PROMPTs. Adjacent promoters driving divergent mRNA transcription support PROMPT formation, but owing to polyadenylation site constraints, these transcripts tend to spread into the neighboring mRNA on the same strand. This mechanism to derive new alternative mRNA transcription start sites (TSSs) is also evident at closely spaced promoters supporting convergent mRNA transcription. We suggest that basic building blocks of divergently transcribed core promoter pairs, in combination with the wealth of TSSs in mammalian genomes, provide a framework with which evolution shapes transcriptomes.

The intergenic region of the maize defensin-like protein genes Def1 and Def2 functions as an embryo-specific asymmetric bidirectional promoter

Bidirectional promoters are identified in diverse organisms with widely varied genome sizes, including bacteria, yeast, mammals, and plants. However, little research has been done on any individual endogenous bidirectional promoter from plants. Here, we describe a promoter positioned in the intergenic region of two defensin-like protein genes, Def1 and Def2 in maize (Zea mays). We examined the expression profiles of Def1 and Def2 in 14 maize tissues by qRT-PCR, and the results showed that this gene pair was expressed abundantly and specifically in seeds. When fused to either green fluorescent protein (GFP) or β-glucuronidase (GUS) reporter genes, P ZmBD1 , P ZmDef1 , and P ZmDef2 were active and reproduced the expression patterns of both Def1 and Def2 genes in transformed immature maize embryos, as well as in developing seeds of transgenic maize. Comparative analysis revealed that PZmBD1 shared most of the expression characteristics of the two polar promoters, but displayed more stringent embryo specificity, delayed expression initiation, and asymmetric promoter activity. Moreover, a truncated promoter study revealed that the core promoters only exhibit basic bidirectional activity, while interacting with necessary cis-elements, which leads to polarity and different strengths. The sophisticated interaction or counteraction between the core promoter and cis-elements may potentially regulate bidirectional promoters.

Long noncoding RNA variations in cardiometabolic diseases

Cardiometabolic diseases are characterized as a combination of multiple risk factors for cardiovascular disease (CVD) and metabolic diseases including diabetes mellitus, dyslipidemia, hypertension and abdominal obesity. This cluster of abnormalities individually and interdependently leads to atherosclerosis and CVD morbidity and mortality. In the past decade, genome-wide association studies (GWASs) have identified a series of cardiometabolic disease-associated variants that can collectively explain a small proportion of the variability. Intriguingly, the susceptibility variants imputed from GWASs usually do not reside in the coding regions, suggesting a crucial role of the noncoding elements of the genome. In recent years, emerging evidence suggests that noncoding RNA (ncRNA) is functional for physiology and pathophysiology of human diseases. These include microRNAs and long noncoding RNAs (lncRNAs) that are now implicated in human diseases. The ncRNAs can interact with each other and with proteins, to interfere gene expressions, leading to the development of many human disorders. Although evidence suggests the functional role of lncRNAs in cardiometabolic traits, the molecular mechanisms of gene regulation underlying cardiometabolic diseases remain to be better defined. Here, we summarize the recent discoveries of lncRNA variations in the context of cardiometabolic diseases.

The Determinants of Directionality in Transcriptional Initiation

A new paradigm has emerged in recent years characterizing transcription initiation as a bidirectional process encompassing a larger proportion of the genome than previously thought. Past concepts of coding genes thinly scattered among a vast background of transcriptionally inert noncoding DNA have been abandoned. A richer picture has taken shape, integrating transcription of coding genes, enhancer RNAs (eRNAs), and various other noncoding transcriptional events. In this review we give an overview of recent studies detailing the mechanisms of RNA polymerase II (RNA Pol II)-based transcriptional initiation and discuss the ways in which transcriptional direction is established as well as its functional implications.

Regulation of cell-to-cell variability in divergent gene expression

Cell-to-cell variability (noise) is an important feature of gene expression that impacts cell fitness and development. The regulatory mechanism of this variability is not fully understood. Here we investigate the effect on gene expression noise in divergent gene pairs (DGPs). We generated reporters driven by divergent promoters, rearranged their gene order, and probed their expressions using time-lapse fluorescence microscopy and single-molecule fluorescence in situ hybridization (smFISH). We show that two genes in a co-regulated DGP have higher expression covariance compared with the separate, tandem and convergent configurations, and this higher covariance is caused by more synchronized firing of the divergent transcriptions. For differentially regulated DGPs, the regulatory signal of one gene can stochastically ‘leak' to the other, causing increased gene expression noise. We propose that the DGPs' function in limiting or promoting gene expression noise may enhance or compromise cell fitness, providing an explanation for the conservation pattern of DGPs.

Gene expression noise affects cell fitness and development. Here, Yan et al. show that co-regulated divergent gene pairs (DGPs) suppress uncorrelated gene expression noise due to more synchronized transcription firing, and differentially regulated DGPs enhance gene expression noise due to transcription leakage.

Landscape and Dynamics of Transcription Initiation in the Malaria Parasite Plasmodium falciparum

A comprehensive map of transcription start sites (TSSs) across the highly AT-rich genome of P. falciparum would aid progress toward deciphering the molecular mechanisms that underlie the timely regulation of gene expression in this malaria parasite. Using high-throughput sequencing technologies, we generated a comprehensive atlas of transcription initiation events at single-nucleotide resolution during the parasite intra-erythrocytic developmental cycle. This detailed analysis of TSS usage enabled us to define architectural features of plasmodial promoters. We demonstrate that TSS selection and strength are constrained by local nucleotide composition. Furthermore, we provide evidence for coordinate and stage-specific TSS usage from distinct sites within the same transcription unit, thereby producing transcript isoforms, a subset of which are developmentally regulated. This work offers a framework for further investigations into the interactions between genomic sequences and regulatory factors governing the complex transcriptional program of this major human pathogen.

Mutational Strand Asymmetries in Cancer Genomes Reveal Mechanisms of DNA Damage and Repair

Mutational processes constantly shape the somatic genome, leading to immunity, aging, cancer, and other diseases. When cancer is the outcome, we are afforded a glimpse into these processes by the clonal expansion of the malignant cell. Here, we characterize a less explored layer of the mutational landscape of cancer: mutational asymmetries between the two DNA strands. Analyzing whole-genome sequences of 590 tumors from 14 different cancer types, we reveal widespread asymmetries across mutagenic processes, with transcriptional ("T-class") asymmetry dominating UV-, smoking-, and liver-cancer-associated mutations and replicative ("R-class") asymmetry dominating POLE-, APOBEC-, and MSI-associated mutations. We report a striking phenomenon of transcription-coupled damage (TCD) on the non-transcribed DNA strand and provide evidence that APOBEC mutagenesis occurs on the lagging-strand template during DNA replication. As more genomes are sequenced, studying and classifying their asymmetries will illuminate the underlying biological mechanisms of DNA damage and repair.

Dissecting the Roles of Divergent and Convergent Transcription in Chromosome Instability

The interplay of transcription, topological tension, and chromosome breakage is a subject of intense interest, but, with so many facets to the problem, it is difficult to test. Here, we vary the orientation of promoters relative to one another in a yeast system that permits sensitive detection of chromosome breaks. Interestingly, convergent transcription that would direct RNA polymerases into one another does not increase chromosome breakage. In contrast, divergent transcription that would create underwound and potentially single-stranded DNA does cause a marked increase in chromosome breakage. Furthermore, we examine the role that topoisomerases are playing in preventing genome instability at these promoters and find that Top2 is required to prevent instability at converging promoters.

Regulatory dissection of the CBX5 and hnRNPA1 bi-directional promoter in human breast cancer cells reveals novel transcript variants differentially associated with HP1α down-regulation in metastatic cells

BACKGROUND

The three members of the human heterochromatin protein 1 (HP1) family of proteins, HP1α, HP1β, and HPγ, are involved in chromatin packing and epigenetic gene regulation. HP1α is encoded from the CBX5 gene and is a suppressor of metastasis. CBX5 is down-regulated at the transcriptional and protein level in metastatic compared to non-metastatic breast cancer. CBX5 shares a bi-directional promoter structure with the hnRNPA1 gene. But whereas CBX5 expression is down-regulated in metastatic cells, hnRNAP1 expression is constant. Here, we address the regulation of CBX5 in human breast cancer.

METHODS

Transient transfection and transposon mediated integration of dual-reporter mini-genes containing the bi-directional hnRNPA1 and CBX5 promoter was performed to investigate transcriptional regulation in breast cancer cell lines. Bioinformatics and functional analysis were performed to characterize transcriptional events specifically regulating CBX5 expression. TSA treatment and Chromatin Immunoprecipitation (ChIP) were performed to investigate the chromatin structure along CBX5 in breast cancer cells. Finally, expression of hnRNPA1 and CBX5 mRNA isoforms were measured by quantitative reverse transcriptase PCR (qRT-PCR) in breast cancer tissue samples.

RESULTS

We demonstrate that an hnRNPA1 and CBX5 bi-directional core promoter fragment does not comprise intrinsic capacity for specific CBX5 down-regulation in metastatic cells. Characterization of transcriptional events in the 20 kb CBX5 intron 1 revealed existence of several novel CBX5 transcripts. Two of these encode consensus HP1α protein but used autonomous promoters in intron 1 by which HP1α expression could be de-coupled from the bi-directional promoter. In addition, another CBX5 transcriptional isoform, STET, was discovered. This transcript includes CBX5 exon 1 and part of intron 1 sequences but lacks inclusion of HP1α encoding exons. Inverse correlation between STET and HP1α coding CBX5 mRNA expression was observed in breast cancer cell lines and tissue samples from breast cancer patients.

CONCLUSION

We find that HP1α is down-regulated in a mechanism involving CBX5 promoter downstream sequences and that regulation through alternative polyadenylation and splicing generates a transcript, STET, with potential importance in carcinogenesis.