An abundance of bidirectional promoters in the human genome

Genetic analysis of cognitive preservation in the midwestern Amish reveals a novel locus on chromosome 2

INTRODUCTION

Alzheimer's disease (AD) remains a debilitating condition with limited treatments and additional therapeutic targets needed. Identifying AD protective genetic loci may identify new targets and accelerate identification of therapeutic treatments. We examined a founder population to identify loci associated with cognitive preservation into advanced age.

METHODS

Genome-wide association and linkage analyses were performed on 946 examined and sampled Amish individuals, aged 76-95, who were either cognitively unimpaired (CU) or impaired (CI).

RESULTS

A total of 12 single nucleotide polymorphisms (SNPs) demonstrated suggestive association (P ≤ 5 × 10-4) with cognitive preservation. Genetic linkage analyses identified > 100 significant (logarithm of the odds [LOD] ≥ 3.3) SNPs, some which overlapped with the association results. Only one locus on chromosome 2 retained significance across multiple analyses.

DISCUSSION

A novel significant result for cognitive preservation on chromosome 2 includes the genes LRRTM4 and CTNNA2. Additionally, the lead SNP, rs1402906, impacts the POU3F2 transcription factor binding affinity, which regulates LRRTM4 and CTNNA2.

HIGHLIGHTS

GWAS and linkage identified over 100 loci associated with cognitive preservation. One locus on Chromosome 2 retained significance over multiple analyses. Predicted TFBSs near rs1402906 regulate genes associated with neurocognition.

Identification and functional characterization of bidirectional gene pairs and their intergenic regions in cotton

BACKGROUND

In research to improve the quality of transgenic crops, it is often necessary to introduce multiple functionally related genes into recipient plants simultaneously to improve crop genetic traits effectively. Compared with unidirectional promoters, bidirectional promoters simultaneously regulate the expression of multiple genes and improve the efficiency of biotechnology. Therefore, in this study, bidirectional gene pairs were systematically analyzed in Gossypium hirsutum TM-1, and the structure, function and evolutionary relationships of the bidirectional genes were analyzed. The endogenous bidirectional promoters of cotton were mined, and their specific regulatory elements and biological functions were explored to provide useful promoter resources and a theoretical basis for cultivating new cotton germplasms with excellent fiber quality.

RESULTS

Using an improved search model, a total of 1,383 bidirectional transcript pairs were identified in the Gossypium hirsutum TM-1 genome, and their gene structure and functional annotations were systematically analyzed. Thirty bidirectional intergenic sequences were randomly screened for promoter activity analysis via a transient expression system, and 25 intergenic sequences were found to have bidirectional promoter activity. Comparative analysis of the bidirectional gene profiles of the four cotton subspecies revealed that these subspecies presented abundant bidirectional gene pairs with high homology and that the bidirectional genes in the cotton subspecies were more similar in terms of their molecular functions, cellular components and biological processes. In addition, parallel analysis of bidirectional genes in dicotyledons and monocotyledons revealed that abundant bidirectional gene pairs exist in different species. Although the total number of orthologous bidirectional genes was similar, there was a significant difference in the number of orthologous bidirectional gene pairs between dicotyledons and monocotyledons. This evolutionary analysis of the function and structure of homologous bidirectional gene pairs in different varieties and different subspecies of the same species revealed potential pathways by which these gene pairs originated, which may be necessary for the evolution of a new species.

CONCLUSION

In this study, many bidirectional gene pairs in Gossypium hirsutum TM-1 were identified using computer programming, and systematic analysis was conducted to explore their functions and evolutionary relationships. In addition, the promoter activity of the bidirectional intergenic sequences was verified. The combination of computer programming screening, experimental validation and other methods is expected to provide preferred bidirectional promoters for transgenic breeding work via multigene cotransformation methods, and this information is valuable for genetic engineering research and applications.

Transcription Attenuation in Synthetic Promoters in Nonoverlapping Tandem Formation

Closely spaced promoters are ubiquitous in prokaryotic and eukaryotic genomes. How their structure and dynamics relate remains unclear, particularly for tandem formations. To study their transcriptional interference, we engineered two pairs and one trio of synthetic promoters in nonoverlapping, tandem formation, in single-copy plasmids transformed into Escherichia coli cells. From in vivo measurements, we found that these promoters in tandem formation can have attenuated transcription rates. The attenuation strength can be widely fine-tuned by the promoters' positioning, natural regulatory mechanisms, and other factors, including the antibiotic rifampicin, which is known to hamper RNAP promoter escape. From this, and supported by in silico models, we concluded that the attenuation in these constructs emerges from premature terminations generated by collisions between RNAPs elongating from upstream promoters and RNAPs occupying downstream promoters. Moreover, we found that these collisions can cause one or both RNAPs to falloff. Finally, the broad spectrum of possible, externally regulated, attenuation strengths observed in our synthetic tandem promoters suggests that they could become useful as externally controllable regulators of future synthetic circuits.

Molecular models of bidirectional promoter regulation

Approximately 11% of human genes are transcribed by a bidirectional promoter (BDP), defined as two genes with <1 kb between their transcription start sites. Despite their evolutionary conservation and enrichment for housekeeping genes and oncogenes, the regulatory role of BDPs remains unclear. BDPs have been suggested to facilitate gene coregulation and/or decrease expression noise. This review discusses these potential regulatory functions through the context of six prospective underlying mechanistic models: a single nucleosome free region, shared transcription factor/regulator binding, cooperative negative supercoiling, bimodal histone marks, joint activation by enhancer(s), and RNA-mediated recruitment of regulators. These molecular mechanisms may act independently and/or cooperatively to facilitate the coregulation and/or decreased expression noise predicted of BDPs.

Novel insights on the positive correlation between sense and antisense pairs on gene expression

A considerable proportion of the eukaryotic genome undergoes transcription, leading to the generation of noncoding RNA molecules that lack protein-coding information and are not subjected to translation. These noncoding RNAs (ncRNAs) are well recognized to have essential roles in several biological processes. Long noncoding RNAs (lncRNAs) represent the most extensive category of ncRNAs found in the human genome. Much research has focused on investigating the roles of cis-acting lncRNAs in the regulation of specific target gene expression. In the majority of instances, the regulation of sense gene expression by its corresponding antisense pair occurs in a negative (discordant) manner, resulting in the suppression of the target genes. The notion that a negative correlation exists between sense and antisense pairings is, however, not universally valid. In fact, several recent studies have reported a positive relationship between corresponding cis antisense pairs within plants, budding yeast, and mammalian cancer cells. The positive (concordant) correlation between anti-sense and sense transcripts leads to an increase in the level of the sense transcript within the same genomic loci. In addition, mechanisms such as altering chromatin structure, the formation of R loops, and the recruitment of transcription factors can either enhance transcription or stabilize sense transcripts through their antisense pairs. The primary objective of this work is to provide a comprehensive understanding of both aspects of antisense regulation, specifically focusing on the positive correlation between sense and antisense transcripts in the context of eukaryotic gene expression, including its implications towards cancer progression. This article is categorized under: RNA Processing > 3' End Processing Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

Cellular control of protein levels: A systems biology perspective

How cells regulate protein levels is a central question of biology. Over the past decades, molecular biology research has provided profound insights into the mechanisms and the molecular machinery governing each step of the gene expression process, from transcription to protein degradation. Recent advances in transcriptomics and proteomics have complemented our understanding of these fundamental cellular processes with a quantitative, systems-level perspective. Multi-omic studies revealed significant quantitative, kinetic and functional differences between the genome, transcriptome and proteome. While protein levels often correlate with mRNA levels, quantitative investigations have demonstrated a substantial impact of translation and protein degradation on protein expression control. In addition, protein-level regulation appears to play a crucial role in buffering protein abundances against undesirable mRNA expression variation. These findings have practical implications for many fields, including gene function prediction and precision medicine.

Development of a compact bidirectional promoter-driven dual chimeric antigen receptor (CAR) construct targeting CD19 and CD20 in the Sleeping Beauty (SB) transposon system

BACKGROUND

A bidirectional promoter-driven chimeric antigen receptor (CAR) cassette provides the simultaneous expression of two CARs, which significantly enhances dual antigen-targeted CAR T-cell therapy.

METHODS

We developed a second-generation CAR directing CD19 and CD20 antigens, incorporating them in a head-to-head orientation from a bidirectional promoter using a single Sleeping Beauty transposon system. The efficacy of bidirectional promoter-driven dual CD19 and CD20 CAR T cells was determined in vitro against cell lines expressing either, or both, CD19 and CD20 antigens. In vivo antitumor activity was tested in Raji lymphoma-bearing immunodeficient NOD-scid IL2Rgammanull (NSG) mice.

RESULTS

Of all tested promoters, the bidirectional EF-1α promoter optimally expressed transcripts from both sense (CD19-CAR) and antisense (GFP.CD20-CAR) directions. Superior cytotoxicity, cytokine production and antigen-specific activation were observed in vitro in the bidirectional EF-1α promoter-driven CD19/CD20 CAR T cells. In contrast, a unidirectional construct driven by the EF-1α promoter, but using self-cleaving peptide-linked CD19 and CD20 CARs, showed inferior expression and in vitro function. Treatment of mice bearing advanced Raji lymphomas with bidirectional EF-1α promoter-driven CD19/CD20 CAR T cells effectively controlled tumor growth and extended the survival of mice compared with group treated with single antigen targeted CAR T cells.

CONCLUSION

The use of bidirectional promoters in a single vector offers advantages of size and robust CAR expression with the potential to expand use in other forms of gene therapies like CAR T cells.

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.

Proteasome Inhibition Reprograms Chromatin Landscape in Breast Cancer

The 26S proteasome is the major protein degradation machinery in cells. Cancer cells use the proteasome to modulate gene expression networks that promote tumor growth. Proteasome inhibitors have emerged as effective cancer therapeutics, but how they work mechanistically remains unclear. Here, using integrative genomic analysis, we discovered unexpected reprogramming of the chromatin landscape and RNA polymerase II (RNAPII) transcription initiation in breast cancer cells treated with the proteasome inhibitor MG132. The cells acquired dynamic changes in chromatin accessibility at specific genomic loci termed differentially open chromatin regions (DOCR). DOCRs with decreased accessibility were promoter proximal and exhibited unique chromatin architecture associated with divergent RNAPII transcription. Conversely, DOCRs with increased accessibility were primarily distal to transcription start sites and enriched in oncogenic superenhancers predominantly accessible in non-basal breast tumor subtypes. These findings describe the mechanisms by which the proteasome modulates the expression of gene networks intrinsic to breast cancer biology.

SIGNIFICANCE

Our study provides a strong basis for understanding the mechanisms by which proteasome inhibitors exert anticancer effects. We find open chromatin regions that change during proteasome inhibition, are typically accessible in non-basal breast cancers.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II often initiates transcription in both directions at the promoter transcription start site, but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements such as enhancers. These bidirectional RNA Polymerase II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the transcription start site or cis-regulatory elements. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both transcription start sites and cis-regulatory elements in A. thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional RNA Polymerase II initiation at both promoters and gene-proximal cis-regulatory elements of A. thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both D. melanogaster and H. sapiens, as well as cis-regulatory elements in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around transcription start sites in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

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.

Molecular analysis of the human cytoglobin mRNA isoforms

Multiple functions have been proposed for the ubiquitously expressed vertebrate globin cytoglobin (Cygb), including nitric oxide (NO) metabolism, lipid peroxidation/signalling, superoxide dismutase activity, reactive oxygen/nitrogen species (RONS) scavenging, regulation of blood pressure, antifibrosis, and both tumour suppressor and oncogenic effects. Since alternative splicing can expand the biological roles of a gene, we investigated whether this mechanism contributes to the functional diversity of Cygb. By mining of cDNA data and molecular analysis, we identified five alternative mRNA isoforms for the human CYGB gene (V-1 to V-5). Comprehensive RNA-seq analyses of public datasets from human tissues and cells confirmed that the canonical CYGB V-1 isoform is the primary CYGB transcript in the majority of analysed datasets. Interestingly, we revealed that isoform V-3 represented the predominant CYGB variant in hepatoblastoma (HB) cell lines and in the majority of analysed normal and HB liver tissues. CYGB V-3 mRNA is transcribed from an alternate upstream promoter and hypothetically encodes a N-terminally truncated CYGB protein, which is not recognized by some antibodies used in published studies. Little to no transcriptional evidence was found for the other CYGB isoforms. Comparative transcriptomics and flow cytometry on CYGB+/+ and gene-edited CYGB-/- HepG2 HB cells did not unveil a knockout phenotype and, thus, a potential function for CYGB V-3. Our study reveals that the CYGB gene is transcriptionally more complex than previously described as it expresses alternative mRNA isoforms of unknown function. Additional experimental data are needed to clarify the biological meaning of those alternative CYGB transcripts.

A plant cell death-inducing protein from litchi interacts with Peronophythora litchii pectate lyase and enhances plant resistance

Cell wall degrading enzymes, including pectate lyases (PeLs), released by plant pathogens, break down protective barriers and/or activate host immunity. The direct interactions between PeLs and plant immune-related proteins remain unclear. We identify two PeLs, PlPeL1 and PlPeL1-like, critical for full virulence of Peronophythora litchii on litchi (Litchi chinensis). These proteins enhance plant susceptibility to oomycete pathogens in a PeL enzymatic activity-dependent manner. However, LcPIP1, a plant immune regulator secreted by litchi, binds to PlPeL1/PlPeL1-like, and attenuates PlPeL1/PlPeL1-like induced plant susceptibility to Phytophthora capsici. LcPIP1 also induces cell death and various immune responses in Nicotiana benthamiana. Conserved in plants, LcPIP1 homologs bear a conserved "VDMASG" motif and exhibit immunity-inducing activity. Furthermore, SERK3 interacts with LcPIP1 and is required for LcPIP1-induced cell death. NbPIP1 participates in immune responses triggered by the PAMP protein INF1. In summary, our study reveals the dual roles of PlPeL1/PlPeL1-like in plant-pathogen interactions: enhancing pathogen virulence through PeL enzymatic activity while also being targeted by LcPIP1, thus enhancing plant immunity.

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.

Genetic analysis of cognitive preservation in the midwestern Amish reveals a novel locus on chromosome 2

INTRODUCTION

Alzheimer disease (AD) remains a debilitating condition with limited treatments and additional therapeutic targets needed. Identifying AD protective genetic loci may identify new targets and accelerate identification of therapeutic treatments. We examined a founder population to identify loci associated with cognitive preservation into advanced age.

METHODS

Genome-wide association and linkage analyses were performed on 946 examined and sampled Amish individuals, aged 76-95, who were either cognitively unimpaired (CU) or impaired (CI).

RESULTS

12 SNPs demonstrated suggestive association (P≤5×10-4) with cognitive preservation. Genetic linkage analyses identified >100 significant (LOD≥3.3) SNPs, some which overlapped with the association results. Only one locus on chromosome 2 retained significance across multiple analyses.

DISCUSSION

A novel significant result for cognitive preservation on chromosome 2 includes the genes LRRTM4 and CTNNA2. Additionally, the lead SNP, rs1402906, impacts the POU3F2 transcription factor binding affinity, which regulates LRRTM4 and CTNNA2.

Generation of de novo miRNAs from template switching during DNA replication

The mechanisms generating novel genes and genetic information are poorly known, even for microRNA (miRNA) genes with an extremely constrained design. All miRNA primary transcripts need to fold into a stem-loop structure to yield short gene products ([Formula: see text]22 nt) that bind and repress their mRNA targets. While a substantial number of miRNA genes are ancient and highly conserved, short secondary structures coding for entirely novel miRNA genes have been shown to emerge in a lineage-specific manner. Template switching is a DNA-replication-related mutation mechanism that can introduce complex changes and generate perfect base pairing for entire hairpin structures in a single event. Here, we show that the template-switching mutations (TSMs) have participated in the emergence of over 6,000 suitable hairpin structures in the primate lineage to yield at least 18 new human miRNA genes, that is 26% of the miRNAs inferred to have arisen since the origin of primates. While the mechanism appears random, the TSM-generated miRNAs are enriched in introns where they can be expressed with their host genes. The high frequency of TSM events provides raw material for evolution. Being orders of magnitude faster than other mechanisms proposed for de novo creation of genes, TSM-generated miRNAs enable near-instant rewiring of genetic information and rapid adaptation to changing environments.

Exploring the role of transcriptional and post-transcriptional processes in mRNA co-expression

Co-expression of two or more genes at the single-cell level is usually associated with functional co-regulation. While mRNA co-expression-measured as the correlation in mRNA levels-can be influenced by both transcriptional and post-transcriptional events, transcriptional regulation is typically considered dominant. We review and connect the literature describing transcriptional and post-transcriptional regulation of co-expression. To enhance our understanding, we integrate four datasets spanning single-cell gene expression data, single-cell promoter activity data and individual transcript half-lives. Confirming expectations, we find that positive co-expression necessitates promoter coordination and similar mRNA half-lives. Surprisingly, negative co-expression is favored by differences in mRNA half-lives, contrary to initial predictions from stochastic simulations. Notably, this association manifests specifically within clusters of genes. We further observe a striking compensation between promoter coordination and mRNA half-lives, which additional stochastic simulations suggest might give rise to the observed co-expression patterns. These findings raise intriguing questions about the functional advantages conferred by this compensation between distal kinetic steps.

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.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II (Pol II) often initiates transcription in both directions at the promoter transcription start site (TSS), but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements (CREs) such as enhancers. These bidirectional Pol II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the TSS or CRE. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both TSSs and CREs in Arabidopsis thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional Pol II initiation at both promoters and gene-proximal CREs of Arabidopsis thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both Drosophila melanogaster and Homo sapiens, as well as CREs in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around TSSs in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

Complex origin, evolution, and diversification of non-canonically organized OVATE-OFP and OVATE-Like OFP gene pair across Embryophyta

Ovate Family Proteins (OFP) is a plant-specific gene family of negative transcriptional regulators. Till-date, a handful of in-silico studies have provided glimpses into family size, expansion patterns, and genic features across all major plant lineages. A major lacuna exists in understanding origin of organisation complexity of members such as those arranged in a head-to-head manner which may lead to transcriptional co-regulation via a common bi-directional promoter. To address this gap, we investigated the origin, organization and evolution of two head-to-head arranged gene pairs of homologs of AtOFP2-AtOFP17, and, AtOFP4-AtOFP20 across Archaeplastida. The ancestral forms of AtOFP2, AtOFP4, AtOFP17, and AtOFP20 are likely to have evolved in last common ancestors of Embryophyta (land plants) given their complete absence in Rhodophyta and Chlorophyta. The OFP gene family originated and expanded in Bryophyta, including protein variants with complete (OVATE-OFP) or partial (OVATE-Like OFP) OVATE domain; with head-to-head organization present only in Spermatophyta (gymnosperms and angiosperms). Ancestral State Reconstruction revealed the origin of head-to-head organized gene pair in gymnosperms, with both genes being OVATE-OFP (homologs of AtOFP2/4). Phylogenetic reconstruction and copy number analysis suggests the presence of a single copy of the head-to-head arranged pair of OFP2/4 (OVATE)-OFP17/20 (OVATE-Like) in all angiosperms except Brassicaceae, and a duplication event in last common ancestor of core Brassicaceae approximately 32-54 MYA leading to origin of AtOFP2-AtOFP17 and AtOFP4-AtOFP20 as paralogs. Synteny analysis of genomic regions harbouring homologs of AtOFP2-AtOFP17, AtOFP4-AtOFP20 and AtOFP2/4-AtOFP17/20 across angiosperms suggested ancestral nature of AtOFP2-AtOFP17 gene pair. The present study thus establishes the orthology and evolutionary history of two non-canonically organised gene pairs with variation in their OVATE domain. The non-canonical organisation, atleast in Brassicaceae, has the potential of generating complex transcriptional regulation mediated via a common bi-directional promoter. The study thus lays down a framework to understand evolution of gene and protein structure, transcriptional regulation and function across a phylogenetic lineage through comparative analyses.

Proteasome inhibition reprograms chromatin landscape in breast cancer

The 26S proteasome is the major protein degradation machinery in cells. Cancer cells use the proteasome to modulate gene expression networks that promote tumor growth. Proteasome inhibitors have emerged as effective cancer therapeutics, but how they work mechanistically remains unclear. Here, using integrative genomic analysis, we discovered unexpected reprogramming of the chromatin landscape and RNAPII transcription initiation in breast cancer cells treated with the proteasome inhibitor MG132. The cells acquired dynamic changes in chromatin accessibility at specific genomic loci termed Differentially Open Chromatin Regions (DOCRs). DOCRs with decreased accessibility were promoter proximal and exhibited unique chromatin architecture associated with divergent RNAPII transcription. Conversely, DOCRs with increased accessibility were primarily distal to transcription start sites and enriched in oncogenic super enhancers predominantly accessible in non-basal breast tumor subtypes. These findings describe the mechanisms by which the proteasome modulates the expression of gene networks intrinsic to breast cancer biology.

Highlights

Proteasome inhibition uncovers de novo Differential Open Chromatin Regions (DOCRs) in breast cancer cells. Proteasome inhibitor sensitive promoters exhibit a distinctive chromatin architecture with discrete transcription initiation patterns.Proteasome inhibition reprograms accessibility of super enhancers.Proteasome inhibitor sensitive super enhancers distinguish basal from non-basal breast cancer subtypes.

A Small De Novo CNV Deletion of the Paternal Copy of FOXF1, Leaving lncRNA FENDRR Intact, Provides Insight into Their Bidirectional Promoter Region

Pathogenic single-nucleotide variants (SNVs) and copy-number variant (CNV) deletions involving the FOXF1 transcription factor gene or CNV deletions of its distant lung-specific enhancer are responsible for alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV), a rarely diagnosed lethal lung developmental disorder in neonates. In contrast to SNVs within FOXF1 and CNV deletions involving only the FOXF1 enhancer, larger-sized deletions involving FOXF1 and the adjacent, oppositely oriented lncRNA gene FENDRR have additionally been associated with hypoplastic left heart syndrome and single umbilical artery (SUA). Here, in an ACDMPV infant without any congenital heart defect or SUA, we identified a small 5 kb CNV deletion that removed the paternal allele of FOXF1 and its promoter, leaving FENDRR and its promoter intact. Reporter assay in the IMR-90 fetal cell line implied that the deletion may indeed not have significantly affected FENDRR expression. It also showed a polarization of the FOXF1-FENDRR inter-promoter region consisting of its ability to increase the transcription of FENDRR but not FOXF1. Interestingly, this transcription-stimulating activity was suppressed in the presence of the FOXF1 promoter. Our data shed more light on the interactions between neighboring promoters of FOXF1-FENDRR and possibly other divergently transcribed mRNA-lncRNA gene pairs.

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.

Higher-order modular regulation of the human proteome

Operons are transcriptional modules that allow bacteria to adapt to environmental changes by coordinately expressing the relevant set of genes. In humans, biological pathways and their regulation are more complex. If and how human cells coordinate the expression of entire biological processes is unclear. Here, we capture 31 higher-order co-regulation modules, which we term progulons, by help of supervised machine-learning on proteomics data. Progulons consist of dozens to hundreds of proteins that together mediate core cellular functions. They are not restricted to physical interactions or co-localisation. Progulon abundance changes are primarily controlled at the level of protein synthesis and degradation. Implemented as a web app at www.proteomehd.net/progulonFinder, our approach enables the targeted search for progulons of specific cellular processes. We use it to identify a DNA replication progulon and reveal multiple new replication factors, validated by extensive phenotyping of siRNA-induced knockdowns. Progulons provide a new entry point into the molecular understanding of biological processes.

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.

The Functional Meaning of 5'UTR in Protein-Coding Genes

As it is well known, messenger RNA has many regulatory regions along its sequence length. One of them is the 5' untranslated region (5'UTR), which itself contains many regulatory elements such as upstream ORFs (uORFs), internal ribosome entry sites (IRESs), microRNA binding sites, and structural components involved in the regulation of mRNA stability, pre-mRNA splicing, and translation initiation. Activation of the alternative, more upstream transcription start site leads to an extension of 5'UTR. One of the consequences of 5'UTRs extension may be head-to-head gene overlap. This review describes elements in 5'UTR of protein-coding transcripts and the functional significance of protein-coding genes 5' overlap with implications for transcription, translation, and disease.

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.

CTCF blocks antisense transcription initiation at divergent promoters

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense transcripts (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is not well understood. Here, by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters. Primary transcript RNA-FISH shows that CTCF lowers burst fraction but not burst intensity of uasTrx and that co-bursting of sense and antisense transcripts is disfavored. Genome editing, chromatin conformation studies and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its architectural function. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

CTCF blocks antisense transcription initiation at divergent promoters

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense transcripts (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is not well understood. Here, by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters. Primary transcript RNA-FISH shows that CTCF lowers burst fraction but not burst intensity of uasTrx and that co-bursting of sense and antisense transcripts is disfavored. Genome editing, chromatin conformation studies and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its architectural function. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

Navigating the Multiverse of Antisense RNAs: The Transcription- and RNA-Dependent Dimension

Evidence accumulated over the past decades shows that the number of identified antisense transcripts is continuously increasing, promoting them from transcriptional noise to real genes with specific functions. Indeed, recent studies have begun to unravel the complexity of the antisense RNA (asRNA) world, starting from the multidimensional mechanisms that they can exert in physiological and pathological conditions. In this review, we discuss the multiverse of the molecular functions of asRNAs, describing their action through transcription-dependent and RNA-dependent mechanisms. Then, we report the workflow and methodologies to study and functionally characterize single asRNA candidates.

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.

Database of Potential Promoter Sequences in the Capsicum annuum Genome

In this study, we used a mathematical method for the multiple alignment of highly divergent sequences (MAHDS) to create a database of potential promoter sequences (PPSs) in the Capsicum annuum genome. To search for PPSs, 20 statistically significant classes of sequences located in the range from −499 to +100 nucleotides near the annotated genes were calculated. For each class, a position–weight matrix (PWM) was computed and then used to identify PPSs in the C. annuum genome. In total, 825,136 PPSs were detected, with a false positive rate of 0.13%. The PPSs obtained with the MAHDS method were tested using TSSFinder, which detects transcription start sites. The databank of the found PPSs provides their coordinates in chromosomes, the alignment of each PPS with the PWM, and the level of statistical significance as a normal distribution argument, and can be used in genetic engineering and biotechnology.

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.

Histone Variants in the Specialization of Plant Chromatin

The basic unit of chromatin, the nucleosome, is an octamer of four core histone proteins (H2A, H2B, H3, and H4) and serves as a fundamental regulatory unit in all DNA-templated processes. The majority of nucleosome assembly occurs during DNA replication when these core histones are produced en masse to accommodate the nascent genome. In addition, there are a number of nonallelic sequence variants of H2A and H3 in particular, known as histone variants, that can be incorporated into nucleosomes in a targeted and replication-independent manner. By virtue of their sequence divergence from the replication-coupled histones, these histone variants can impart unique properties onto the nucleosomes they occupy and thereby influence transcription and epigenetic states, DNA repair, chromosome segregation, and other nuclear processes in ways that profoundly affect plant biology. In this review, we discuss the evolutionary origins of these variants in plants, their known roles in chromatin, and their impacts on plant development and stress responses. We focus on the individual and combined roles of histone variants in transcriptional regulation within euchromatic and heterochromatic genome regions. Finally, we highlight gaps in our understanding of plant variants at the molecular, cellular, and organismal levels, and we propose new directions for study in the field of plant histone variants.

Promoter antisense RNAs: beyond transcription by-products of active promoters

Omics-based technologies, which have developed rapidly over the last few decades, have generated increasing evidence demonstrating pervasive divergent transcription from RNA polymerase II (Pol II) promoters of eukaryotic genome, and indeed have raised considerable discussion as to their potential physiopathological function. Unlike many other long non-coding RNAs (lncRNAs), promoter antisense RNAs (PAS RNAs) were initially considered to be merely passive transcription by-products of active promoters. However, recent studies have begun to reveal their critical importance in a broad spectrum of biological processes. In this Review, I summarize recent technological advances that enable accurate detection of PAS RNA and discuss the mechanisms of PAS RNA biogenesis emphasizing the functional importance of its structure enabling the diverse functions of PAS RNA in transcription and chromatin regulation.

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.

Recent Advances in Gene Therapy for Cancer Theranostics

There is great interest in developing gene therapies for many disease indications, including cancer. However, successful delivery of nucleic acids to tumor cells is a major challenge, and in vivo efficacy is difficult to predict. Cancer theranostics is an approach combining anti-tumor therapy with imaging or diagnostic capabilities, with the goal of monitoring successful delivery and efficacy of a therapeutic agent in a tumor. Successful theranostics must maintain a high degree of anticancer targeting and efficacy while incorporating high-contrast imaging agents that are nontoxic and compatible with clinical imaging modalities. This review highlights recent advancements in theranostic strategies, including imaging technologies and genetic engineering approaches. Graphical Abstract.

Refining the genomic location of single nucleotide polymorphism variation affecting Atlantic salmon maturation timing at a key large-effect locus

Efforts to understand the genetic underpinnings of phenotypic variation are becoming more and more frequent in molecular ecology. Such efforts often lead to the identification of candidate regions showing signals of association and/or selection. These regions may contain multiple genes and therefore validation of which genes are actually responsible for the signal is required. In Atlantic salmon (Salmo salar), a large-effect locus for maturation timing, an ecologically important trait, occurs in a genomic region including two genes, vgll3 and akap11, but data for clearly determining which of the genes (or both) contribute to the association have been lacking. Here, we take advantage of natural recombination events detected between the two candidate genes in a salmon broodstock to reduce linkage disequilibrium at the locus, thus enabling delineation of the influence of variation at these two genes on early maturation. By rearing 5,895 males to maturation age, of which 81% had recombinant vgll3/akap11 allelic combinations, we found that vgll3 single nucleotide polymorphism (SNP) variation was strongly associated with early maturation, whereas there was little or no association between akap11 SNP variation and early maturation. These findings provide strong evidence supporting vgll3 as the primary candidate gene in the chromosome 25 locus for influencing early maturation. This will help guide future research for understanding the genetic processes controlling early maturation. This also exemplifies the utility of natural recombinants to more precisely map causal variation underlying ecologically important phenotypic diversity.

X-Chromosome Inactivation and Autosomal Random Monoallelic Expression as "Faux Amis"

X-chromosome inactivation (XCI) and random monoallelic expression of autosomal genes (RMAE) are two paradigms of gene expression regulation where, at the single cell level, genes can be expressed from either the maternal or paternal alleles. X-chromosome inactivation takes place in female marsupial and placental mammals, while RMAE has been described in mammals and also other species. Although the outcome of both processes results in random monoallelic expression and mosaicism at the cellular level, there are many important differences. We provide here a brief sketch of the history behind the discovery of XCI and RMAE. Moreover, we review some of the distinctive features of these two phenomena, with respect to when in development they are established, their roles in dosage compensation and cellular phenotypic diversity, and the molecular mechanisms underlying their initiation and stability.

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.

Transcriptomic Analysis in Multiple Myeloma and Primary Plasma Cell Leukemia with t(11;14) Reveals Different Expression Patterns with Biological Implications in Venetoclax Sensitivity

Simple Summary

The growing interest in BCL2 inhibitors for the treatment of multiple myeloma (MM) has led to the need for biomarkers that are able to predict patient’s sensitivity to the drug. The presence of the chromosomal translocation t(11;14) in MM is mainly associated with sensitivity to venetoclax and good prognosis. The incidence of t(11;14) largely increases in primary Plasma Cell Leukemia (pPCL) in association with an unfavorable outcome. Currently, data concerning pPCL sensitivity to venetoclax are virtually absent. In this context, we investigated the transcriptome of MM and pPCL with t(11;14), evidencing that the two clinical entities are likely responsive to venetoclax based on different molecular programs, thus prompting further studies to elucidate better novel potential predictive biomarkers.

Abstract

Mechanisms underlying the pathophysiology of primary Plasma Cell Leukemia (pPCL) and intramedullary multiple myeloma (MM) need to be further elucidated, being potentially relevant for improving therapeutic approaches. In such a context, the MM and pPCL subgroups characterized by t(11;14) deserve a focused investigation, as the presence of the translocation is mainly associated with sensitivity to venetoclax. Herein, we investigated a proprietary cohort of MM and pPCL patients, focusing on the transcriptional signature of samples carrying t(11;14), whose incidence increases in pPCL in association with an unfavorable outcome. In addition, we evaluated the expression levels of the BCL2-gene family members and of a panel of B-cell genes recently reported to be associated with sensitivity to venetoclax in MM. Moreover, transcriptional analysis of lncRNAs in the two clinical settings led to the identification of several differentially expressed transcripts, among which the SNGH6 deregulated lncRNA might be relevant in the pathogenesis and prognosis of pPCL with t(11;14). Overall, our data suggest that MMs and pPCLs with t(11;14) might be responsive to venetoclax based on different molecular programs, prompting further studies to elucidate better novel potential predictive biomarkers.

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.

Genetic and epigenetic features of promoters with ubiquitous chromatin accessibility support ubiquitous transcription of cell-essential genes

Gene expression is controlled by regulatory elements within accessible chromatin. Although most regulatory elements are cell type-specific, a subset is accessible in nearly all the 517 human and 94 mouse cell and tissue types assayed by the ENCODE consortium. We systematically analyzed 9000 human and 8000 mouse ubiquitously-accessible candidate cis-regulatory elements (cCREs) with promoter-like signatures (PLSs) from ENCODE, which we denote ubi-PLSs. These are more CpG-rich than non-ubi-PLSs and correspond to genes with ubiquitously high transcription, including a majority of cell-essential genes. ubi-PLSs are enriched with motifs of ubiquitously-expressed transcription factors and preferentially bound by transcriptional cofactors regulating ubiquitously-expressed genes. They are highly conserved between human and mouse at the synteny level but exhibit frequent turnover of motif sites; accordingly, ubi-PLSs show increased variation at their centers compared with flanking regions among the ∼186 thousand human genomes sequenced by the TOPMed project. Finally, ubi-PLSs are enriched in genes implicated in Mendelian diseases, especially diseases broadly impacting most cell types, such as deficiencies in mitochondrial functions. Thus, a set of roughly 9000 mammalian promoters are actively maintained in an accessible state across cell types by a distinct set of transcription factors and cofactors to ensure the transcriptional programs of cell-essential genes.

Promoter switching in response to changing environment and elevated expression of protein-coding genes overlapping at their 5' ends

Despite the number of studies focused on sense-antisense transcription, the key question of whether such organization evolved as a regulator of gene expression or if this is only a byproduct of other regulatory processes has not been elucidated to date. In this study, protein-coding sense-antisense gene pairs were analyzed with a particular focus on pairs overlapping at their 5' ends. Analyses were performed in 73 human transcription start site libraries. The results of our studies showed that the overlap between genes is not a stable feature and depends on which TSSs are utilized in a given cell type. An analysis of gene expression did not confirm that overlap between genes causes downregulation of their expression. This observation contradicts earlier findings. In addition, we showed that the switch from one promoter to another, leading to genes overlap, may occur in response to changing environment of a cell or tissue. We also demonstrated that in transfected and cancerous cells genes overlap is observed more often in comparison with normal tissues. Moreover, utilization of overlapping promoters depends on particular state of a cell and, at least in some groups of genes, is not merely coincidental.

Co-regulation and function of FOXM1/RHNO1 bidirectional genes in cancer

The FOXM1 transcription factor is an oncoprotein and a top biomarker of poor prognosis in human cancer. Overexpression and activation of FOXM1 is frequent in high-grade serous carcinoma (HGSC), the most common and lethal form of human ovarian cancer, and is linked to copy number gains at chromosome 12p13.33. We show that FOXM1 is co-amplified and co-expressed with RHNO1, a gene involved in the ATR-Chk1 signaling pathway that functions in the DNA replication stress response. We demonstrate that FOXM1 and RHNO1 are head-to-head (i.e., bidirectional) genes (BDG) regulated by a bidirectional promoter (BDP) (named F/R-BDP). FOXM1 and RHNO1 each promote oncogenic phenotypes in HGSC cells, including clonogenic growth, DNA homologous recombination repair, and poly-ADP ribosylase inhibitor resistance. FOXM1 and RHNO1 are one of the first examples of oncogenic BDG, and therapeutic targeting of FOXM1/RHNO1 BDG is a potential therapeutic approach for ovarian and other cancers.

Crystallin gene expression: Insights from studies of transcriptional bursting

Cellular differentiation is marked by temporally and spatially regulated gene expression. The ocular lens is one of the most powerful mammalian model system since it is composed from only two cell subtypes, called lens epithelial and fiber cells. Lens epithelial cells differentiate into fiber cells through a series of spatially and temporally orchestrated processes, including massive production of crystallins, cellular elongation and the coordinated degradation of nuclei and other organelles. Studies of transcriptional and posttranscriptional gene regulatory mechanisms in lens provide a wide range of opportunities to understand global molecular mechanisms of gene expression as steady-state levels of crystallin mRNAs reach very high levels comparable to globin genes in erythrocytes. Importantly, dysregulation of crystallin gene expression results in lens structural abnormalities and cataracts. The mRNA life cycle is comprised of multiple stages, including transcription, splicing, nuclear export into cytoplasm, stabilization, localization, translation and ultimate decay. In recent years, development of modern mRNA detection methods with single molecule and single cell resolution enabled transformative studies to visualize the mRNA life cycle to generate novel insights into the sequential regulatory mechanisms of gene expression during embryogenesis. This review is focused on recent major advancements in studies of transcriptional bursting in differentiating lens fiber cells, analysis of nascent mRNA expression from bi-directional promoters, transient nuclear accumulation of specific mRNAs, condensation of chromatin prior lens fiber cell denucleation, and outlines future studies to probe the interactions of individual mRNAs with specific RNA-binding proteins (RBPs) in the cytoplasm and regulation of translation and mRNA decay.

The Dm-Myb Oncoprotein Contributes to Insulator Function and Stabilizes Repressive H3K27me3 PcG Domains

Drosophila Myb (Dm-Myb) encodes a protein that plays a key role in regulation of mitotic phase genes. Here, we further refine its role in the context of a developing tissue as a potentiator of gene expression required for proper RNA polymerase II (RNA Pol II) function and efficient H3K4 methylation at promoters. In contrast to its role in gene activation, Myb is also required for repression of many genes, although no specific mechanism for this role has been proposed. We now reveal a critical role for Myb in contributing to insulator function, in part by promoting binding of insulator proteins BEAF-32 and CP190 and stabilizing H3K27me3 Polycomb-group (PcG) domains. In the absence of Myb, H3K27me3 is markedly reduced throughout the genome, leading to H3K4me3 spreading and gene derepression. Finally, Myb is enriched at boundaries that demarcate chromatin environments, including chromatin loop anchors. These results reveal functions of Myb that extend beyond transcriptional regulation.

Myb has been considered a transcriptional activator of primarily M phase genes. Here, Santana et al. show that Myb also contributes to insulator function, in part by promoting binding of insulator factors, and is required to stabilize repressive domains in the genome.