Bidirectional gene organization: a common architectural feature of the human genome

A Bidirectional Non-Coding RNA Promoter Mediates Long-Range Gene Expression Regulation

Recent evidence suggests that human gene promoters display gene expression regulatory mechanisms beyond the typical single gene local transcription modulation. In mammalian genomes, genes with an associated bidirectional promoter are abundant; bidirectional promoter architecture serves as a regulatory hub for a gene pair expression. However, it has been suggested that its contribution to transcriptional regulation might exceed local transcription initiation modulation. Despite their abundance, the functional consequences of bidirectional promoter architecture remain largely unexplored. This work studies the long-range gene expression regulatory role of a long non-coding RNA gene promoter using chromosome conformation capture methods. We found that this particular bidirectional promoter contributes to distal gene expression regulation in a target-specific manner by establishing promoter-promoter interactions. In particular, we validated that the promoter-promoter interactions of this regulatory element with the promoter of distal gene BBX contribute to modulating the transcription rate of this gene; removing the bidirectional promoter from its genomic context leads to a rearrangement of BBX promoter-enhancer interactions and to increased gene expression. Moreover, long-range regulatory functionality is not directly dependent on its associated non-coding gene pair expression levels.

Economy of Effort or Sophisticated Programming? The Prevalence of Bidirectional Promoter Complexes in the Human Genome

An abundance of antisense promoters in the vicinity of the transcriptional start site of coding genes suggests that they play an important role in gene regulation. The divergent transcription of housekeeping genes by a common central promoter region allows for coordinated regulation of genes in related pathways and is also linked to higher promoter activity. However, closely positioned transcription start sites can also result in competition between overlapping promoter elements and generate a binary switch element. Furthermore, the direct competition resulting from the presence of an antisense promoter immediately downstream of the transcription start site of the gene produces an element that can exist in only one of two stable transcriptional states: sense or antisense. In this review, we summarize analyses of the prevalence of antisense transcription in higher eukaryotes and viruses, with a focus on the antisense promoters competing with the promoters of coding genes. The structures of bidirectional promoters driving the simultaneous expression of housekeeping genes are compared with examples of human bidirectional elements that have been shown to act as switches. Since many bidirectional elements contain a noncoding RNA as the divergent transcript, we describe examples of functional noncoding antisense transcripts that affect the epigenetic landscape and alter the expression of their host gene. Finally, we discuss opportunities for additional research on competing sense/antisense promoters, uncovering their potential role in programming cell differentiation.

Cell cycle regulation of the psoriasis associated gene CCHCR1 by transcription factor E2F1

The coiled-coil alpha-helical rod protein 1 (CCHCR1) was first identified as a candidate gene in psoriasis and has lately been found to be associated with a wide range of clinical conditions including COVID-19. CCHCR1 is located within P-bodies and centrosomes, but its exact role in these two subcellular structures and its transcriptional control remain largely unknown. Here, we showed that CCHCR1 shares a bidirectional promoter with its neighboring gene, TCF19. This bidirectional promoter is activated by the G1/S-regulatory transcription factor E2F1, and both genes are co-induced during the G1/S transition of the cell cycle. A luciferase reporter assay suggests that the short intergenic sequence, only 287 bp in length, is sufficient for the G1/S induction of both genes, but the expression of CCHCR1 is further enhanced by the presence of exon 1 from both TCF19 and CCHCR1. This research uncovers the transcriptional regulation of the CCHCR1 gene, offering new perspectives on its function. These findings contribute to the broader understanding of diseases associated with CCHCR1 and may serve as a foundational benchmark for future research in these vital medical fields.

Abundant binary promoter switches in lineage-determining transcription factors indicate a digital component of cell fate determination

Previous studies of the murine Ly49 and human KIR gene clusters implicated competing sense and antisense promoters in the control of variegated gene expression. In the current study, an examination of transcription factor genes defines an abundance of convergent and divergent sense/antisense promoter pairs, suggesting that competing promoters may control cell fate determination. Differentiation of CD34+ hematopoietic progenitors in vitro shows that cells with GATA1 antisense transcription have enhanced GATA2 transcription and a mast cell phenotype, whereas cells with GATA2 antisense transcription have increased GATA1 transcripts and an erythroblast phenotype. Detailed analyses of the AHR and RORC genes demonstrate the ability of competing promoters to act as binary switches and the association of antisense transcription with an immature/progenitor cell phenotype. These data indicate that alternative cell fates generated by promoter competition in lineage-determining transcription factors contribute to the programming of cell differentiation.

Molecular cloning of the gene promoter encoding the human CaVγ2/Stargazin divergent transcript (CACNG2-DT): characterization and regulation by the cAMP-PKA/CREB signaling pathway

CaVγ2 (Stargazin or TARPγ2) is a protein expressed in various types of neurons whose function was initially associated with a decrease in the functional expression of voltage-gated presynaptic Ca2+ channels (CaV) and which is now known to promote the trafficking of the postsynaptic α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptors (AMPAR) towards the cell membrane. Alterations in CaVγ2 expression has been associated with several neurological disorders, such as absence epilepsy. However, its regulation at the transcriptional level has not been intensively addressed. It has been reported that the promoter of the Cacng2 gene, encoding the rat CaVγ2, is bidirectional and regulates the transcription of a long non-coding RNA (lncRNA) in the antisense direction. Here, we investigate the proximal promoter region of the human CACNG2 gene in the antisense direction and show that this region includes two functional cAMP response elements that regulate the expression of a lncRNA called CACNG2-DT. The activity of these sites is significantly enhanced by forskolin, an adenylate cyclase activator, and inhibited by H89, a protein kinase A (PKA) antagonist. Therefore, this regulatory mechanism implies the activation of G protein-coupled receptors and downstream phosphorylation. Interestingly, we also found that the expression of CACNG2-DT may increase the levels of the CaVγ2 subunit. Together, these data provide novel information on the organization of the human CACNG2-DT gene promoter, describe modulatory domains and mechanisms that can mediate various regulatory inputs, and provide initial information on the molecular mechanisms that regulate the functional expression of the CaVγ2 protein.

CCIVR2 facilitates comprehensive identification of both overlapping and non-overlapping antisense transcripts within specified regions

Pairs of sense and antisense transcriptions that are adjacent at their 5' and 3' regions are called divergent and convergent transcription, respectively. However, the structural properties of divergent/convergent transcription in different species or RNA biotypes are poorly characterized. Here, we developed CCIVR2, a program that facilitates identification of both overlapping and non-overlapping antisense transcripts produced from divergent/convergent transcription whose transcription start sites (TSS) or transcript end sites (TES) are located within a specified region. We used CCIVR2 to analyze antisense transcripts starting around the sense TSS (from divergent transcription) or ending around the sense TES (from convergent transcription) in 11 different species and found species- and RNA biotype-specific features of divergent/convergent transcription. Furthermore, we confirmed that CCIVR2 enables the identification of multiple sense/antisense transcript pairs from divergent transcription, including those with known functions in processes such as embryonic stem cell differentiation and TGFβ stimulation. CCIVR2 is therefore a valuable bioinformatics tool that facilitates the characterization of divergent/convergent transcription in different species and aids the identification of functional sense/antisense transcript pairs from divergent transcription in specified biological processes.

Long noncoding RNAs: biogenesis, regulation, function, and their emerging significance in toxicology

The repertoire of regulatory noncoding RNAs (ncRNAs) has been enriched by the inclusion of long noncoding RNA (lncRNA) that are longer than 200 nt. Some of the currently known lncRNAs, were reported in the 1990s before the term lncRNA was introduced. These lncRNAs have diverse regulatory functions including regulation of transcription via interactions with proteins and RNAs, chromatin remodeling, translation, posttranslational protein modification, protein trafficking and cell signaling. Predictably, the dysregulation of lncRNA expression due to exposure to toxicants may precipitate adverse health consequences. Dysregulation of lncRNAs has also been implicated in various adverse human health outcomes. There is an increasing agreement that lncRNA expression profiling data needs to be closely examined to determine whether their altered expression can be used as biomarkers of toxicity as well as adverse human health outcomes. This review summarizes the biogenesis, regulation, function of lncRNA and their emerging significance in toxicology and disease conditions. Because our understanding of the lncRNA-toxicity relationship is still evolving, this review discusses this developing field using some examples.

Strand Asymmetries Across Genomic Processes

Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.

Transforming Diagnosis and Therapeutics Using Cancer Genomics

In past quarter of the century, much has been understood about the genetic variation and abnormal genes that activate cancer in humans. All the cancers somehow possess alterations in the DNA sequence of cancer cell's genome. In present, we are heading toward the era where it is possible to obtain complete genome of the cancer cells for their better diagnosis, categorization and to explore treatment options.

SETD1A regulates transcriptional pause release of heme biosynthesis genes in leukemia

Histone methyltransferase SETD1A is critical for acute myeloid leukemia (AML) cell survival, but the molecular mechanism driving SETD1A gene regulation remains elusive. To delineate the role of SETD1A, we utilize a protein degrader technology to induce rapid SETD1A degradation in AML cell lines. SETD1A degradation results in immediate downregulation of transcripts associated with DNA repair and heme biosynthesis pathways. CRISPR-based functional analyses and metabolomics reveal an essential role of SETD1A to maintain mitochondrial respiration in AML cells. These SETD1A targets are enriched in head-to-head (H2H) genes. SETD1A degradation disrupts a non-enzymatic SETD1A domain-dependent cyclin K function, increases the Ser5P RNA polymerase II (RNAPII) at the transcriptional start site (TSS), and induces the promoter-proximal pausing of RNAPII in a strand-specific manner. This study reveals a non-enzymatic role for SETD1A in transcriptional pause release and provides insight into the mechanism of RNAPII pausing and its function in cancer.

Core circadian clock and light signaling genes brought into genetic linkage across the green lineage

The circadian clock is conserved at both the level of transcriptional networks as well as core genes in plants, ensuring that biological processes are phased to the correct time of day. In the model plant Arabidopsis (Arabidopsis thaliana), the core circadian SHAQKYF-type-MYB (sMYB) genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and REVEILLE (RVE4) show genetic linkage with PSEUDO-RESPONSE REGULATOR 9 (PRR9) and PRR7, respectively. Leveraging chromosome-resolved plant genomes and syntenic ortholog analysis enabled tracing this genetic linkage back to Amborella trichopoda, a sister lineage to the angiosperm, and identifying an additional evolutionarily conserved genetic linkage in light signaling genes. The LHY/CCA1-PRR5/9, RVE4/8-PRR3/7, and PIF3-PHYA genetic linkages emerged in the bryophyte lineage and progressively moved within several genes of each other across an array of angiosperm families representing distinct whole-genome duplication and fractionation events. Soybean (Glycine max) maintained all but two genetic linkages, and expression analysis revealed the PIF3-PHYA linkage overlapping with the E4 maturity group locus was the only pair to robustly cycle with an evening phase, in contrast to the sMYB-PRR morning and midday phase. While most monocots maintain the genetic linkages, they have been lost in the economically important grasses (Poaceae), such as maize (Zea mays), where the genes have been fractionated to separate chromosomes and presence/absence variation results in the segregation of PRR7 paralogs across heterotic groups. The environmental robustness model is put forward, suggesting that evolutionarily conserved genetic linkages ensure superior microhabitat pollinator synchrony, while wide-hybrids or unlinking the genes, as seen in the grasses, result in heterosis, adaptation, and colonization of new ecological niches.

Pan-Cancer Analysis Reveals the Prognostic Potential of the THAP9/THAP9-AS1 Sense-Antisense Gene Pair in Human Cancers

Human THAP9, which encodes a domesticated transposase of unknown function, and lncRNA THAP9-AS1 (THAP9-antisense1) are arranged head-to-head on opposite DNA strands, forming a sense and antisense gene pair. We predict that there is a bidirectional promoter that potentially regulates the expression of THAP9 and THAP9-AS1. Although both THAP9 and THAP9-AS1 are reported to be involved in various cancers, their correlative roles on each other’s expression has not been explored. We analyzed the expression levels, prognosis, and predicted biological functions of the two genes across different cancer datasets (TCGA, GTEx). We observed that although the expression levels of the two genes, THAP9 and THAP9-AS1, varied in different tumors, the expression of the gene pair was strongly correlated with patient prognosis; higher expression of the gene pair was usually linked to poor overall and disease-free survival. Thus, THAP9 and THAP9-AS1 may serve as potential clinical biomarkers of tumor prognosis. Further, we performed a gene co-expression analysis (using WGCNA) followed by a differential gene correlation analysis (DGCA) across 22 cancers to identify genes that share the expression pattern of THAP9 and THAP9-AS1. Interestingly, in both normal and cancer samples, THAP9 and THAP9-AS1 often co-express; moreover, their expression is positively correlated in each cancer type, suggesting the coordinated regulation of this H2H gene pair.

Neighbours matter: Effects of genomic organization on gene expression plasticity in response to environmental stresses during biological invasions

Gene expression regulation has been widely recognized as an important molecular mechanism underlying phenotypic plasticity in environmental adaptation. However, it remains largely unexplored on the effects of genomic organization on gene expression plasticity under environmental stresses during biological invasions. Here, we use an invasive model ascidian, Ciona robusta, to investigate how genomic organization affects gene expression in response to salinity stresses during range expansions. Our study showed that neighboring genes were co-expressed and approximately 30% of stress responsive genes were physically clustered on chromosomes. Such coordinated expression was substantially affected by the physical distance and orientation of genes. Interestingly, the overall expression correlation of neighboring genes was significantly decreased under high salinity stresses, illustrating that the co-expression regulation could be disrupted by salinity challenges. Furthermore, the clustering of genes was associated with their function constraints and expression patterns - operon genes enriched in gene expression machinery had the highest transcriptional activity and expression stability. Notably, our analyses showed that the tail-to-tail genes, mainly involved in biological functions related to phosphorylation, homeostatic process, and ion transport, exhibited higher intrinsic expression variability and greater response to salinity challenges. Altogether, the results obtained here provide new insights into the effects of gene organization on gene expression plasticity under environmental challenges, hence improving our knowledge on mechanisms of rapid environmental adaptation during biological invasions.

A versatile transposon-based technology to generate loss- and gain-of-function phenotypes in the mouse liver

BACKGROUND

Understanding the contribution of gene function in distinct organ systems to the pathogenesis of human diseases in biomedical research requires modifying gene expression through the generation of gain- and loss-of-function phenotypes in model organisms, for instance, the mouse. However, methods to modify both germline and somatic genomes have important limitations that prevent easy, strong, and stable expression of transgenes. For instance, while the liver is remarkably easy to target, nucleic acids introduced to modify the genome of hepatocytes are rapidly lost, or the transgene expression they mediate becomes inhibited due to the action of effector pathways for the elimination of exogenous DNA. Novel methods are required to overcome these challenges, and here we develop a somatic gene delivery technology enabling long-lasting high-level transgene expression in the entire hepatocyte population of mice.

RESULTS

We exploit the fumarylacetoacetate hydrolase (Fah) gene correction-induced regeneration in Fah-deficient livers, to demonstrate that such approach stabilizes luciferase expression more than 5000-fold above the level detected in WT animals, following plasmid DNA introduction complemented by transposon-mediated chromosomal gene transfer. Building on this advancement, we created a versatile technology platform for performing gene function analysis in vivo in the mouse liver. Our technology allows the tag-free expression of proteins of interest and silencing of any arbitrary gene in the mouse genome. This was achieved by applying the HADHA/B endogenous bidirectional promoter capable of driving well-balanced bidirectional expression and by optimizing in vivo intronic artificial microRNA-based gene silencing. We demonstrated the particular usefulness of the technology in cancer research by creating a p53-silenced and hRas G12V-overexpressing tumor model.

CONCLUSIONS

We developed a versatile technology platform for in vivo somatic genome editing in the mouse liver, which meets multiple requirements for long-lasting high-level transgene expression. We believe that this technology will contribute to the development of a more accurate new generation of tools for gene function analysis in mice.

Analytical kinetic model of native tandem promoters in E. coli

Closely spaced promoters in tandem formation are abundant in bacteria. We investigated the evolutionary conservation, biological functions, and the RNA and single-cell protein expression of genes regulated by tandem promoters in E. coli. We also studied the sequence (distance between transcription start sites ‘d_TSS’, pause sequences, and distances from oriC) and potential influence of the input transcription factors of these promoters. From this, we propose an analytical model of gene expression based on measured expression dynamics, where RNAP-promoter occupancy times and d_TSS are the key regulators of transcription interference due to TSS occlusion by RNAP at one of the promoters (when d_TSS ≤ 35 bp) and RNAP occupancy of the downstream promoter (when d_TSS > 35 bp). Occlusion and downstream promoter occupancy are modeled as linear functions of occupancy time, while the influence of d_TSS is implemented by a continuous step function, fit to in vivo data on mean single-cell protein numbers of 30 natural genes controlled by tandem promoters. The best-fitting step is at 35 bp, matching the length of DNA occupied by RNAP in the open complex formation. This model accurately predicts the squared coefficient of variation and skewness of the natural single-cell protein numbers as a function of d_TSS. Additional predictions suggest that promoters in tandem formation can cover a wide range of transcription dynamics within realistic intervals of parameter values. By accurately capturing the dynamics of these promoters, this model can be helpful to predict the dynamics of new promoters and contribute to the expansion of the repertoire of expression dynamics available to synthetic genetic constructs.

Tandem promoters are common in nature, but investigations on their dynamics have so far largely relied on synthetic constructs. Thus, their regulation and potentially unique dynamics remain unexplored. We first performed a comprehensive exploration of the conservation of genes regulated by these promoters in E. coli and the properties of their input transcription factors. We then measured protein and RNA levels expressed by 30 Escherichia coli tandem promoters, to establish an analytical model of the expression dynamics of genes controlled by such promoters. We show that start site occlusion and downstream RNAP occupancy can be realistically captured by a model with RNAP binding affinity, the time length of open complex formation, and the nucleotide distance between transcription start sites. This study contributes to a better understanding of the unique dynamics tandem promoters can bring to the dynamics of gene networks and will assist in their use in synthetic genetic circuits.

Regulatory mechanisms ensuring coordinated expression of functionally related genes

Coordinated spatiotemporal expression of large sets of genes is required for the development and homeostasis of organisms. To achieve this goal, organisms use myriad strategies where they form operons, utilize bidirectional promoters, cluster genes, share enhancers among genes by DNA looping, and form topologically associated domains and transcriptional condensates. Coexpression achieved by these different strategies is hypothesized to have functional importance in minimizing gene expression variability, establishing dosage balance to ensure stoichiometry of protein complexes, and minimizing accumulation of toxic intermediate metabolites. By combining gene-editing tools with computational modeling, recent studies tested the advantages of adjacent genes located in pairs and clusters. We propose that with the advancement of gene editing, single-cell sequencing, and imaging tools, one could readily test the functional importance of different coexpression strategies in a variety of biological processes.

RN7SK small nuclear RNA controls bidirectional transcription of highly expressed gene pairs in skin

Pausing of RNA polymerase II (Pol II) close to promoters is a common regulatory step in RNA synthesis, and is coordinated by a ribonucleoprotein complex scaffolded by the noncoding RNA RN7SK. The function of RN7SK-regulated gene transcription in adult tissue homoeostasis is currently unknown. Here, we deplete RN7SK during mouse and human epidermal stem cell differentiation. Unexpectedly, loss of this small nuclear RNA specifically reduces transcription of numerous cell cycle regulators leading to cell cycle exit and differentiation. Mechanistically, we show that RN7SK is required for efficient transcription of highly expressed gene pairs with bidirectional promoters, which in the epidermis co-regulated cell cycle and chromosome organization. The reduction in transcription involves impaired splicing and RNA decay, but occurs in the absence of chromatin remodelling at promoters and putative enhancers. Thus, RN7SK is directly required for efficient Pol II transcription of highly transcribed bidirectional gene pairs, and thereby exerts tissue-specific functions, such as maintaining a cycling cell population in the epidermis.

The noncoding RNA RN7SK regulates RNA polymerase II pausing and splicing. Here the authors deplete RN7SK in mouse and human during epidermal stem cell differentiation and reveal a novel role in orchestrating bidirectional transcription of highly expressed gene pairs.

Bidirectional promoters: an enigmatic genome architecture and their roles in cancers

Bidirectional promoters are the transcription regulatory regions of genes positioned head-to-head on opposite strands. Specific sequence signals, chromatin modifications and three-dimensional structures of the transcription site facilitate the unconventional yet tightly regulated transcription proceeding in both directions from these promoters. Mutations or aberrant epigenetic changes can lead to abnormal enhanced or reduced expression from either of the bidirectionally transcribed genes resulting in tumorigenesis. Moreover, bidirectionally transcribed genes might also contribute towards the immune regulation in tumor microenvironment. In this review, we aimed to expound the characteristic features of bidirectional promoters alongside their transcriptional regulations, and ultimately, the association of these enigmatic genomic elements in different cancers.

High-Resolution Mapping of Transcription Initiation in the Asexual Stages of Toxoplasma gondii

Toxoplasma gondii is a common parasite of humans and animals, causing life-threatening disease in the immunocompromized, fetal abnormalities when contracted during gestation, and recurrent ocular lesions in some patients. Central to the prevalence and pathogenicity of this protozoan is its ability to adapt to a broad range of environments, and to differentiate between acute and chronic stages. These processes are underpinned by a major rewiring of gene expression, yet the mechanisms that regulate transcription in this parasite are only partially characterized. Deciphering these mechanisms requires a precise and comprehensive map of transcription start sites (TSSs); however, Toxoplasma TSSs have remained incompletely defined. To address this challenge, we used 5'-end RNA sequencing to genomically assess transcription initiation in both acute and chronic stages of Toxoplasma. Here, we report an in-depth analysis of transcription initiation at promoters, and provide empirically-defined TSSs for 7603 (91%) protein-coding genes, of which only 1840 concur with existing gene models. Comparing data from acute and chronic stages, we identified instances of stage-specific alternative TSSs that putatively generate mRNA isoforms with distinct 5' termini. Analysis of the nucleotide content and nucleosome occupancy around TSSs allowed us to examine the determinants of TSS choice, and outline features of Toxoplasma promoter architecture. We also found pervasive divergent transcription at Toxoplasma promoters, clustered within the nucleosomes of highly-symmetrical phased arrays, underscoring chromatin contributions to transcription initiation. Corroborating previous observations, we asserted that Toxoplasma 5' leaders are among the longest of any eukaryote studied thus far, displaying a median length of approximately 800 nucleotides. Further highlighting the utility of a precise TSS map, we pinpointed motifs associated with transcription initiation, including the binding sites of the master regulator of chronic-stage differentiation, BFD1, and a novel motif with a similar positional arrangement present at 44% of Toxoplasma promoters. This work provides a critical resource for functional genomics in Toxoplasma, and lays down a foundation to study the interactions between genomic sequences and the regulatory factors that control transcription in this parasite.

Pan-Cancer Analysis of Head-to-Head Gene Pairs in Terms of Transcriptional Activity, Co-expression and Regulation

Background

Head-to-Head (H2H) gene pairs are regulated by bidirectional promoters and divergently transcribed from opposite DNA strands with transcription start sites (TSSs) separated within 1 kb. H2H organization is ancient and conserved, and H2H pairs tend to exhibit similar expression patterns. Although some H2H genes have been reported to be associated with disease and cancer, there is a lack of systematic studies on H2H organization in the scenario of cancer development.

Methods

Human H2H gene pairs were identified based on GENCODE hg19 and the functional relevance of H2H pairs was explored through function enrichment and semantic similarity analysis. To investigate the association between H2H organization and carcinogenesis, pan-cancer differential analysis of H2H genes about transcriptional activity, co-expression and transcriptional regulation by transcription factors and enhancers were performed based on data from The Cancer Genome Atlas. Cox proportional hazards regression model and log-rank test were used to determine the prognostic powers of H2H pairs.

Results

In the present study, we first updated H2H genes from 1,447 to 3,150 pairs, from which the peak group with TSS distance of 1–100 was observed as expected in our previous work. It was found that housekeeping genes, mitochondrial-functional associated genes and cancer genes tend to be organized in H2H arrangement. Pan-cancer analysis indicates that H2H genes are transcriptionally active than random genes in both normal and cancer tissues, but H2H pairs display higher correlation in cancer than in normal. Particularly, housekeeping H2H pairs are differentially correlated much more significantly than non-housekeeping H2H pairs are. Some of differentially correlated H2H pairs were found to be associated with prognosis. The alteration of TF similarity seems to contribute to differential co-expression of H2H pairs during carcinogenesis; meanwhile remote enhancers also at least partly explain the differential co-expression and co-regulation of H2H pairs.

Conclusion

H2H pairs tend to show much stronger positive expression correlation in cancer than in normal due to differential regulation of bidirectional promoters. The study provides insights into the significance of H2H organization in carcinogenesis and the underlying dysfunctional regulation mechanisms. Those differentially correlated H2H pairs associated with survival have the potential to be prognostic biomarkers and therapeutic targets for cancer.

Pairing of segmentation clock genes drives robust pattern formation

Gene expression is an inherently stochastic process ^1,2; however, organismal development and homeostasis require that cells spatiotemporally coordinate the expression of large sets of genes. Coexpressed gene pairs in metazoans often reside in the same chromosomal neighborhood, with gene pairs representing 10% - 50% of all genes depending on species ^3–6. As shared upstream regulators can ensure correlated gene expression, the selective advantage of maintaining adjacent gene pairs remains unknown ⁶. Here, using two linked zebrafish segmentation clock genes, her1 and her7, and combining single-cell transcript counting, genetic engineering, real-time imaging and computational modeling, we reveal that gene pairing boosts correlated transcription and provides phenotypic robustness for developmental pattern formation. Our results demonstrate that disrupting gene pairing disrupts oscillations and segmentation, identifying the selective pressure retaining correlated transcription to sustain a robust and rapid developmental clock. We anticipate that these findings will inspire investigating advantages of gene pairing in other systems and engineering precise synthetic clocks in embryos and organoids.

Therapeutic Angiogenesis by a "Dynamic Duo": Simultaneous Expression of HGF and VEGF165 by Novel Bicistronic Plasmid Restores Blood Flow in Ischemic Skeletal Muscle

Therapeutic angiogenesis is a promising strategy for relief of ischemic conditions, and gene delivery was used to stimulate blood vessels’ formation and growth. We have previously shown that intramuscular injection of a mixture containing plasmids encoding vascular endothelial growth factor (VEGF)165 and hepatocyte growth factor (HGF) leads to restoration of blood flow in mouse ischemic limb, and efficacy of combined delivery was superior to each plasmid administered alone. In this work, we evaluated different approaches for co-expression of HGF and VEGF165 genes in a panel of candidate plasmid DNAs (pDNAs) with internal ribosome entry sites (IRESs), a bidirectional promoter or two independent promoters for each gene of interest. Studies in HEK293T culture showed that all plasmids provided synthesis of HGF and VEGF165 proteins and stimulated capillary formation by human umbilical vein endothelial cells (HUVEC), indicating the biological potency of expressed factors. Tests in skeletal muscle explants showed a dramatic difference and most plasmids failed to express HGF and VEGF165 in a significant quantity. However, a bicistronic plasmid with two independent promoters (cytomegalovirus (CMV) for HGF and chicken b-actin (CAG) for VEGF165) provided expression of both grow factors in skeletal muscle at an equimolar ratio. Efficacy tests of bicistronic plasmid were performed in a mouse model of hind limb ischemia. Intramuscular administration of plasmid induced significant restoration of perfusion compared to an empty vector and saline. These findings were supported by increased CD31+ capillary density in animals that received pHGF/VEGF. Overall, our study reports a first-in-class candidate gene therapy drug to deliver two pivotal angiogenic growth factors (HGF and VEGF165) with properties that provide basis for future development of treatment for an unmet medical need—peripheral artery disease and associated limb ischemia.

LncRNA-MAP3K4 regulates vascular inflammation through the p38 MAPK signaling pathway and cis-modulation of MAP3K4

Chronic vascular inflammation plays a key role in the pathogenesis of atherosclerosis. Long non-coding RNAs (lncRNAs) have emerged as essential inflammation regulators. We identify a novel lncRNA termed lncRNA-MAP3K4 that is enriched in the vessel wall and regulates vascular inflammation. In the aortic intima, lncRNA-MAP3K4 expression was reduced by 50% during the progression of atherosclerosis (chronic inflammation) and 70% during endotoxemia (acute inflammation). lncRNA-MAP3K4 knockdown reduced the expression of key inflammatory factors (eg, ICAM-1, E-selectin, MCP-1) in endothelial cells or vascular smooth muscle cells and decreased monocytes adhesion to endothelium, as well as reducing TNF-α, IL-1β, COX2 expression in macrophages. Mechanistically, lncRNA-MAP3K4 regulates inflammation through the p38 MAPK signaling pathway. lncRNA-MAP3K4 shares a bidirectional promoter with MAP3K4, an upstream regulator of the MAPK signaling pathway, and regulates its transcription in cis. lncRNA-MAP3K4 and MAP3K4 show coordinated expression in response to inflammation in vivo and in vitro. Similar to lncRNA-MAP3K4, MAP3K4 knockdown reduced the expression of inflammatory factors in several different vascular cells. Furthermore, lncRNA-MAP3K4 and MAP3K4 knockdown showed cooperativity in reducing inflammation in endothelial cells. Collectively, these findings unveil the role of a novel lncRNA in vascular inflammation by cis-regulating MAP3K4 via a p38 MAPK pathway.

Transcription closed and open complex formation coordinate expression of genes with a shared promoter region

Many genes are spaced closely, allowing coordination without explicit control through shared regulatory elements and molecular interactions. We study the dynamics of a stochastic model of a gene-pair in a head-to-head configuration, sharing promoter elements, which accounts for the rate-limiting steps in transcription initiation. We find that only in specific regions of the parameter space of the rate-limiting steps is orderly coexpression exhibited, suggesting that successful cooperation between closely spaced genes requires the coevolution of compatible rate-limiting step configuration. The model predictions are validated using in vivo single-cell, single-RNA measurements of the dynamics of pairs of genes sharing promoter elements. Our results suggest that, in E. coli, the kinetics of the rate-limiting steps in active transcription can play a central role in shaping the dynamics of gene-pairs sharing promoter elements.

Characterization of bidirectional gene pairs in The Cancer Genome Atlas (TCGA) dataset

The “bidirectional gene pair” indicates a particular head-to-head gene organization in which transcription start sites of two genes are located on opposite strands of genomic DNA within a region of one kb. Despite bidirectional gene pairs are well characterized, little is known about their expression profiles and regulation features in tumorigenesis. We used RNA-seq data from The Cancer Genome Atlas (TCGA) dataset for a systematic analysis of the expression profiles of bidirectional gene pairs in 13 cancer datasets. Gene pairs on the opposite strand with transcription end site distance within one kb or on the same strand with the distance of two genes between 1–10 kb and gene pairs comprising two randomly chosen genes were used as control gene pairs (CG1, CG2, and random). We identified and characterized up-/down-regulated genes by comparing the expression level between tumors and adjacent normal tissues in 13 TCGA datasets. There were no consistently significant difference in the percentage of up-/down-regulated genes between bidirectional and control/random genes in most of TCGA datasets. However, the percentage of bidirectional gene pairs comprising two up- or two down-regulated genes was significantly higher than gene pairs from CG1/2 in 12/11 analyzed TCGA datasets and the random gene pairs in all 13 TCGA datasets. Then we identified the methylation correlated bidirectional genes to explore the regulatory mechanism of bidirectional genes. Like the differentially expressed gene pairs, the bidirectional genes in a pair were significantly prone to be both hypo- or hyper-methylation correlated genes in 12/13 TCGA datasets when comparing to the CG2/random gene pairs despite no significant difference between the percentages of hypo-/hyper-methylation correlated genes in bidirectional and CG2/random genes in most of TCGA datasets. Finally, we explored the correlation between bidirectional genes and patient’s survival, identifying prognostic bidirectional genes and prognostic bidirectional gene pairs in each TCGA dataset. Remarkably, we found a group of prognostic bidirectional gene pairs in which the combination of two protein coding genes with different expression level correlated with different survival prognosis in survival analysis for OS. The percentage of these gene pairs in bidirectional gene pair were significantly higher than the gene pairs in controls in COAD datasets and lower in none of 13 TCGA datasets.

Regulatory Architecture of the Neuronal Cacng2/Tarpγ2 Gene Promoter: Multiple Repressive Domains, a Polymorphic Regulatory Short Tandem Repeat, and Bidirectional Organization with Co-regulated lncRNAs

CACNG2 (TARPγ2, Stargazin) is a multi-functional regulator of excitatory neurotransmission and has been implicated in the pathological processes of several brain diseases. Cacng2 function is dependent upon expression level, but currently, little is known about the molecular mechanisms that control expression of this gene. To address this deficit and investigate disease-related gene variants, we have cloned and characterized the rat Cacng2 promoter and have defined three major features: (i) multiple repressive domains that include an array of RE-1 silencing transcription factor (REST) elements, and a calcium regulatory element-binding factor (CaRF) element, (ii) a (poly-GA) short tandem repeat (STR), and (iii) bidirectional organization with expressed lncRNAs. Functional activity of the promoter was demonstrated in transfected neuronal cell lines (HT22 and PC12), but although selective removal of REST and CaRF domains was shown to enhance promoter-driven transcription, the enhanced Cacng2 promoter constructs were still about fivefold weaker than a comparable rat Synapsin-1 promoter sequence. Direct evidence of REST activity at the Cacng2 promoter was obtained through co-transfection with an established dominant-negative REST (DNR) construct. Investigation of the GA-repeat STR revealed polymorphism across both animal strains and species, and size variation was also observed in absence epilepsy disease model cohorts (Genetic Absence Epilepsy Rats, Strasbourg [GAERS] and non-epileptic control [NEC] rats). These data provide evidence of a genotype (STR)-phenotype correlation that may be unique with respect to proximal gene regulatory sequence in the demonstrated absence of other promoter, or 3' UTR variants in GAERS rats. However, although transcriptional regulatory activity of the STR was demonstrated in further transfection studies, we did not find a GAERS vs. NEC difference, indicating that this specific STR length variation may only be relevant in the context of other (Cacna1h and Kcnk9) gene variants in this disease model. Additional studies revealed further (bidirectional) complexity at the Cacng2 promoter, and we identified novel, co-regulated, antisense rat lncRNAs that are paired with Cacng2 mRNA. These studies have provided novel insights into the organization of a synaptic protein gene promoter, describing multiple repressive and modulatory domains that can mediate diverse regulatory inputs.

Synthesis of Seed-Specific Bidirectional Promoters for Metabolic Engineering of Anthocyanin-Rich Maize

Tissue-specific promoters play an important role in plant molecular farming. Here, we describe a strategy to modify the tissue specificity of a maize embryo-specific bidirectional promoter PZmBD1. Six types of cis-elements, i.e. RY repeats (R), GCN4 (G), the prolamin box (P), Skn-1 (S), and the ACGT and AACA (A) motifs, were collected and fused to PZmBD1 to generate eight chimeric putative bidirectional promoters. Qualitative and quantitative analysis of reporter genes driven by the promoters showed that two promoters exhibited high seed-specific bidirectional activity in maize transient and stable transformed systems. The stronger one was chosen and fused to the intergenic region of two gene clusters consisting of four anthocyanin biosynthesis-related genes (ZmBz1, ZmBz2, ZmC1 and ZmR2) and seven reporter genes, resulting in the first embryo and endosperm anthocyanin-rich purple maize. Anthocyanin analysis showed that the total anthocyanin content reaches 2,910 mg kg-1 DW in transgenic maize and cyanidin is the major anthocyanin in transgenic maize, as in natural varieties. The expression profile analysis of endogenous genes showed that the anthocyanin biosynthesis pathway was activated by two transgenic transcription factor genes ZmC1 and ZmR2. Our results indicate that both the modification strategy and these functionally characterized tissue-specific bidirectional promoters generated could be used for genetic research and development of plant biotechnology products. The anthocyanin-rich purple maize could provide economic natural colorants for the food and beverage industry, and valuable germplasm for developing anthocyanin-rich fresh corn.

Engineered bidirectional promoters enable rapid multi-gene co-expression optimization

Numerous synthetic biology endeavors require well-tuned co-expression of functional components for success. Classically, monodirectional promoters (MDPs) have been used for such applications, but MDPs are limited in terms of multi-gene co-expression capabilities. Consequently, there is a pressing need for new tools with improved flexibility in terms of genetic circuit design, metabolic pathway assembly, and optimization. Here, motivated by nature's use of bidirectional promoters (BDPs) as a solution for efficient gene co-expression, we generate a library of 168 synthetic BDPs in the yeast Komagataella phaffii (syn. Pichia pastoris), leveraging naturally occurring BDPs as a parts repository. This library of synthetic BDPs allows for rapid screening of diverse expression profiles and ratios to optimize gene co-expression, including for metabolic pathways (taxadiene, β-carotene). The modular design strategies applied for creating the BDP library could be relevant in other eukaryotic hosts, enabling a myriad of metabolic engineering and synthetic biology applications.

Engineered bidirectional promoters enable rapid multi-gene co-expression optimization

Numerous synthetic biology endeavors require well-tuned co-expression of functional components for success. Classically, monodirectional promoters (MDPs) have been used for such applications, but MDPs are limited in terms of multi-gene co-expression capabilities. Consequently, there is a pressing need for new tools with improved flexibility in terms of genetic circuit design, metabolic pathway assembly, and optimization. Here, motivated by nature’s use of bidirectional promoters (BDPs) as a solution for efficient gene co-expression, we generate a library of 168 synthetic BDPs in the yeast Komagataella phaffii (syn. Pichia pastoris), leveraging naturally occurring BDPs as a parts repository. This library of synthetic BDPs allows for rapid screening of diverse expression profiles and ratios to optimize gene co-expression, including for metabolic pathways (taxadiene, β-carotene). The modular design strategies applied for creating the BDP library could be relevant in other eukaryotic hosts, enabling a myriad of metabolic engineering and synthetic biology applications.

Classic monodirectional promoters are of limited use for multiple gene co-expression. Here the authors generate a library of 168 bidirectional promoters for the yeast K. phaffii (syn. P. pastoris) with diverse expression profiles to optimize metabolic pathway design.

Microbes involving in carcinogenesis; growing state of the art

Lateral gene transfer (LGT) has been demonstrated as a transfer process of novel genes between different species. LGT proceedings are occurring between microbes and plants, as well as between microbes and animals. New evidence demonstrates that bacterial insertional mutagenesis may occur in cancer cells. Due to the important role of genetic changes in the increase of cell proliferation and cancer development, we reviewed the effects of microbial-animal LGT in human oncogenesis. In addition, viral DNA can induce cancer development by random insertion into cancer-related genes or by inducing translocations. In conclusion, growing evidence shows the contribution of the microbial genome in cancer and autoimmune disease.

Control of human gene expression: High abundance of divergent transcription in genes containing both INR and BRE elements in the core promoter

Background

DNA sequence elements in the core promoter can play a central role in regulation of gene expression. Core elements (e.g. INR and TATA box) are located within ~50bp of the transcription start site and both upstream and downstream elements are known. Although all can affect the level of gene expression, their mechanism of action has yet to be fully defined. The studies described here are focused on two core promoter elements, INR and BRE, in the human genome. The locations of the two elements were determined in a large number of human promoters and the results were interpreted in terms of overall promoter function.

Results

A total of 13,406 promoters were collected from the reference version of the human genome and found to contain 62,891 INR sequences and 32,290 BRE. An INR sequence was found in the core region of 1231 (9.2%) promoters and a BRE in 2592 (19.3%); 121 promoters (0.9%) have both INR and BRE elements. Counts support the view that most human promoters lack an INR or BRE element in the core promoter. Further analysis was carried out with the aligned aggregate of promoters from each chromosome. The results showed distinct INR distributions in separate chromosome groups indicating a degree of chromosome specificity to the way core promoter elements are deployed in the genome. The rare promoters with both INR and BRE elements were found to be enriched among the genes with divergent transcription. Enrichment raises the possibility that core promoter elements can have a function in chromosome organization as well as in initiation of transcription.

Long Non-Coding RNA-Mediated Regulation of the Interferon Response: A New Perspective on a Familiar Theme

The interferon (IFN) response is a critical and ubiquitous component of the innate immune response to pathogens. Detailed studies in the last decades have elucidated the function of a large number of proteins that mediate the complex signaling pathways and gene expression programs involved in the interferon response. The recent discovery of the long non-coding RNAs (lncRNAs) as a new category of cellular effectors has led to studies aiming to understand the role of these transcripts in the IFN response. Several high throughput studies have shown that a large number of lncRNAs are differentially expressed following IFN stimulation and/or viral infections. In-depth study of a very small fraction of the identified lncRNAs has revealed critical roles for this class of transcripts in the regulation of multiple steps of the IFN response, and pointed to the presence of an extensive RNA-mediated regulatory network during the antiviral response. As the vast majority of the identified potential regulatory lncRNAs remain unstudied, it is highly likely that future studies will reveal a completely new perspective on the regulation of the IFN response, with lncRNA- and protein-mediated regulatory networks coordinating the duration, magnitude, and character of this aspect of the innate immune response. In addition to providing a more complete picture of the IFN response, these studies will likely identify new therapeutic targets that in the long term may impact the therapeutic options available against microbial infections and diseases of the immune system.

Pan-Cancer Analysis Reveals Differential Susceptibility of Bidirectional Gene Promoters to DNA Methylation, Somatic Mutations, and Copy Number Alterations

Bidirectional gene promoters affect the transcription of two genes, leading to the hypothesis that they should exhibit protection against genetic or epigenetic changes in cancer. Therefore, they provide an excellent opportunity to learn about promoter susceptibility to somatic alteration in tumors. We tested this hypothesis using data from genome-scale DNA methylation (14 cancer types), simple somatic mutation (10 cancer types), and copy number variation profiling (14 cancer types). For DNA methylation, the difference in rank differential methylation between tumor and tumor-adjacent normal matched samples based on promoter type was tested by the Wilcoxon rank sum test. Logistic regression was used to compare differences in simple somatic mutations. For copy number alteration, a mixed effects logistic regression model was used. The change in methylation between non-diseased tissues and their tumor counterparts was significantly greater in single compared to bidirectional promoters across all 14 cancer types examined. Similarly, the extent of copy number alteration was greater in single gene compared to bidirectional promoters for all 14 cancer types. Furthermore, among 10 cancer types with available simple somatic mutation data, bidirectional promoters were slightly more susceptible. These results suggest that selective pressures related with specific functional impacts during carcinogenesis drive the susceptibility of promoter regions to somatic alteration.

Bidirectional promoters exhibit characteristic chromatin modification signature associated with transcription elongation in both sense and antisense directions

Background

In contrast to unidirectional promoters wherein antisense transcription results in short transcripts which are rapidly degraded, bidirectional promoters produce mature transcripts in both sense and antisense orientation. To understand the molecular mechanism of how productive bidirectional transcription is regulated, we focused on delineating the chromatin signature of bidirectional promoters.

Results

We report generation and utility of a reporter system that enables simultaneous scoring of transcriptional activity in opposite directions. Testing of putative bidirectional promoters in this system demonstrates no measurable bias towards any one direction of transcription. We analyzed the NUP26L-PIH1D3 bidirectional gene pair during Retinoic acid mediated differentiation of embryonic carcinoma cells. In their native context, we observed that the chromatin landscape at and around the transcription regulatory region between the pair of bidirectional genes is modulated in concordance with transcriptional activity of each gene in the pair. We then extended this analysis to 974 bidirectional gene pairs in two different cell lines, H1 human embryonic stem cells and CD4 positive T cells using publicly available ChIP-Seq and RNA-Seq data. Bidirectional gene pairs were classified based on the intergenic distance separating the two TSS of the transcripts analyzed as well as the relative expression of each transcript in a bidirectional gene pair. We report that for the entire range of intergenic distance separating bidirectional genes, the expression profile of such genes (symmetric or asymmetric) matches the histone modification profile of marks associated with active transcription initiation and elongation.

Conclusions

We demonstrate unique distribution of histone modification marks that correlate robustly with the transcription status of genes regulated by bidirectional promoters. These findings strongly imply that occurrence of these marks might signal the transcription machinery to drive maturation of antisense transcription from the bidirectional promoters.

Electronic supplementary material

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

CpG Islands: A Historical Perspective

The discovery of CpG islands (CGIs) and the study of their structure and properties run parallel to the development of molecular biology in the last two decades of the twentieth century and to the development of high-throughput genomic technologies at the turn of the millennium. First identified as discrete G + C-rich regions of unmethylated DNA in several vertebrates, CGIs were soon found to display additional distinctive chromatin features from the rest of the genome in terms of accessibility and of the epigenetic modifications of their histones. These features, together with their colocalization with promoters and with origins of DNA replication in mammals, highlighted their relevance in the regulation of genomic processes. Recent approaches have shown with unprecedented detail the dynamics and diversity of the epigenetic landscape of CGIs during normal development and under pathological conditions. Also, comparative analyses across species have started revealing how CGIs evolve and contribute to the evolution of the vertebrate genome.

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.

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.

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.

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.

A mutated dph3 gene causes sensitivity of Schizosaccharomyces pombe cells to cytotoxic agents

Dph3 is involved in diphthamide modification of the eukaryotic translation elongation factor eEF2 and in Elongator-mediated modifications of tRNAs, where a 5-methoxycarbonyl-methyl moiety is added to wobble uridines. Lack of such modifications affects protein synthesis due to inaccurate translation of mRNAs at ribosomes. We have discovered that integration of markers at the msh3 locus of Schizosaccharomyces pombe impaired the function of the nearby located dph3 gene. Such integrations rendered cells sensitive to the cytotoxic drugs hydroxyurea and methyl methanesulfonate. We constructed dph3 and msh3 strains with mutated ATG start codons (ATGmut), which allowed investigating drug sensitivity without potential interference by marker insertions. The dph3-ATGmut and a dph3::loxP-ura4-loxM gene disruption strain, but not msh3-ATGmut, turned out to be sensitive to hydroxyurea and methyl methanesulfonate, likewise the strains with cassettes integrated at the msh3 locus. The fungicide sordarin, which inhibits diphthamide modified eEF2 of Saccharomyces cerevisiae, barely affected survival of wild type and msh3Δ S. pombe cells, while the dph3Δ mutant was sensitive. The msh3-ATG mutation, but not dph3Δ or the dph3-ATG mutation caused a defect in mating-type switching, indicating that the ura4 marker at the dph3 locus did not interfere with Msh3 function. We conclude that Dph3 is required for cellular resistance to the fungicide sordarin and to the cytotoxic drugs hydroxyurea and methyl methanesulfonate. This is likely mediated by efficient translation of proteins in response to DNA damage and replication stress.

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.

Bidirectional transcription of Linc00441 and RB1 via H3K27 modification-dependent way promotes hepatocellular carcinoma

The retinoblastoma gene (RB1), a known tumor-suppressor gene (TSG), was decreased in multiple cancers including hepatocellular carcinoma (HCC). Here we focused on the bidirectional transcripted long noncoding RNA (Linc00441) with neighbor gene RB1 to investigate whether Linc00441 is involved in the suppression of RB1 in HCC. We found that aberrant upregulated intranuclear Linc00441 was reversely correlated with RB1 expression in human HCC samples. The gain- and loss-of-function investigation revealed that Linc00441 could promote the proliferation of HCC cells in vitro and in vivo with an apoptosis suppression and cell cycle rearrangement. Furthermore, RNA pull-down assay indicated the decreased level of RB1 induced by Linc00441 was associated with the incidental methylation by DNMT3A recruited by Linc00441. On the contrary, the transcription factor (TCF-4) enhanced H3K27 acetylation and direct transcription factor for Linc00441 was responsible for the upregulation of Linc00441 in HCC. In conclusion, the epigenetic interaction between Linc00441 and bidirectional transcripted neighbor RB1 may be a de novo theory cutting-point for the inactivation of RB1 in HCC and may serve as targeting site for tumor therapy in the future.

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.

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.

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.

Genes Linked to Endometriosis by GWAS Are Integral to Cytoskeleton Regulation and Suggests That Mesothelial Barrier Homeostasis Is a Factor in the Pathogenesis of Endometriosis

Endometriosis, defined by the presence of ectopic endometrial lesions, is a common disease in reproductive-age women that profoundly affects patients' quality of life. Various pathogenic models have been proposed, but the origin of endometriosis remains elusive. In this article, we propose that the mesothelial barrier, which protects the underlying stroma from endometrial transplants present in retrograde menstrual fluid, can be compromised by activation of the epithelial to mesenchymal transition (EMT) repair mechanism that lead to temporary loss of barrier integrity. Absent of the mesothelial barrier, endometrial cells can more readily adhere to the underlying peritoneal stroma and establish endometrial lesions. The hypothesis is based on the clinical and experimental observations that correlate the location of endometrial lesions with areas of mesothelial damage, together with genetic evidence that 4 genes associated with endometriosis are direct regulators of the actin-cytoskeleton, which coordinates mesothelial barrier integrity. It supports past observations that implicate the peritoneum in the pathogenesis of endometriosis and unifies previously disparate theories that endometriosis may be triggered by infection, mechanical damage, and inflammation since each of these mechanisms can induce EMT in the mesothelium. If the hypothesis is correct, inhibition of EMT in the mesothelial barrier provides a novel paradigm for the prevention and treatment of endometriosis.