The functional consequences of alternative promoter use in mammalian genomes

CLN3 transcript complexity revealed by long-read RNA sequencing analysis

BACKGROUND

Batten disease is a group of rare inherited neurodegenerative diseases. Juvenile CLN3 disease is the most prevalent type, and the most common pathogenic variant shared by most patients is the "1-kb" deletion which removes two internal coding exons (7 and 8) in CLN3. Previously, we identified two transcripts in patient fibroblasts homozygous for the 1-kb deletion: the 'major' and 'minor' transcripts. To understand the full variety of disease transcripts and their role in disease pathogenesis, it is necessary to first investigate CLN3 transcription in "healthy" samples without juvenile CLN3 disease.

METHODS

We leveraged PacBio long-read RNA sequencing datasets from ENCODE to investigate the full range of CLN3 transcripts across various tissues and cell types in human control samples. Then we sought to validate their existence using data from different sources.

RESULTS

We found that a readthrough gene affects the quantification and annotation of CLN3. After taking this into account, we detected over 100 novel CLN3 transcripts, with no dominantly expressed CLN3 transcript. The most abundant transcript has median usage of 42.9%. Surprisingly, the known disease-associated 'major' transcripts are detected. Together, they have median usage of 1.5% across 22 samples. Furthermore, we identified 48 CLN3 ORFs, of which 26 are novel. The predominant ORF that encodes the canonical CLN3 protein isoform has median usage of 66.7%, meaning around one-third of CLN3 transcripts encode protein isoforms with different stretches of amino acids. The same ORFs could be found with alternative UTRs. Moreover, we were able to validate the translational potential of certain transcripts using public mass spectrometry data.

CONCLUSION

Overall, these findings provide valuable insights into the complexity of CLN3 transcription, highlighting the importance of studying both canonical and non-canonical CLN3 protein isoforms as well as the regulatory role of UTRs to fully comprehend the regulation and function(s) of CLN3. This knowledge is essential for investigating the impact of the 1-kb deletion and rare pathogenic variants on CLN3 transcription and disease pathogenesis.

Protein isoform-centric therapeutics: expanding targets and increasing specificity

Most protein-coding genes produce multiple protein isoforms; however, these isoforms are commonly neglected in drug discovery. The expression of protein isoforms can be specific to a disease, tissue and/or developmental stage, and this specific expression can be harnessed to achieve greater drug specificity than pan-targeting of all gene products and to enable improved treatments for diseases caused by aberrant protein isoform production. In recent years, several protein isoform-centric therapeutics have been developed. Here, we collate these studies and clinical trials to highlight three distinct but overlapping modes of action for protein isoform-centric drugs: isoform switching, isoform introduction or depletion, and modulation of isoform activity. In addition, we discuss how protein isoforms can be used clinically as targets for cell type-specific drug delivery and immunotherapy, diagnostic biomarkers and sources of cancer neoantigens. Collectively, we emphasize the value of a focus on isoforms as a route to discovering drugs with greater specificity and fewer adverse effects. This approach could enable the targeting of proteins for which pan-inhibition of all isoforms is toxic and poorly tolerated.

TransTEx: novel tissue-specificity scoring method for grouping human transcriptome into different expression groups

MOTIVATION

Although human tissues carry out common molecular processes, gene expression patterns can distinguish different tissues. Traditional informatics methods, primarily at the gene level, overlook the complexity of alternative transcript variants and protein isoforms produced by most genes, changes in which are linked to disease prognosis and drug resistance.

RESULTS

We developed TransTEx (Transcript-level Tissue Expression), a novel tissue-specificity scoring method, for grouping transcripts into four expression groups. TransTEx applies sequential cut-offs to tissue-wise transcript probability estimates, subsampling-based P-values and fold-change estimates. Application of TransTEx on GTEx mRNA-seq data divided 199 166 human transcripts into different groups as 17 999 tissue-specific (TSp), 7436 tissue-enhanced, 36 783 widely expressed (Wide), 79 191 lowly expressed (Low), and 57 757 no expression (Null) transcripts. Testis has the most (13 466) TSp isoforms followed by liver (890), brain (701), pituitary (435), and muscle (420). We found that the tissue specificity of alternative transcripts of a gene is predominantly influenced by alternate promoter usage. By overlapping brain-specific transcripts with the cell-type gene-markers in scBrainMap database, we found that 63% of the brain-specific transcripts were enriched in nonneuronal cell types, predominantly astrocytes followed by endothelial cells and oligodendrocytes. In addition, we found 61 brain cell-type marker genes encoding a total of 176 alternative transcripts as brain-specific and 22 alternative transcripts as testis-specific, highlighting the complex TSp and cell-type specific gene regulation and expression at isoform-level. TransTEx can be adopted to the analysis of bulk RNA-seq or scRNA-seq datasets to find tissue- and/or cell-type specific isoform-level gene markers.

AVAILABILITY AND IMPLEMENTATION

TransTEx database: https://bmi.cewit.stonybrook.edu/transtexdb/ and the R package is available via GitHub: https://github.com/pallavisurana1/TransTEx.

Integrative analysis of ultra-deep RNA-seq reveals alternative promoter usage as a mechanism of activating oncogenic programmes during prostate cancer progression

Transcription factor (TF) proteins regulate gene activity by binding to regulatory regions, most importantly at gene promoters. Many genes have alternative promoters (APs) bound by distinct TFs. The role of differential TF activity at APs during tumour development is poorly understood. Here we show, using deep RNA sequencing in 274 biopsies of benign prostate tissue, localized prostate tumours and metastatic castration-resistant prostate cancer, that AP usage increases as tumours progress and APs are responsible for a disproportionate amount of tumour transcriptional activity. Expression of the androgen receptor (AR), the key driver of prostate tumour activity, is correlated with elevated AP usage. We identified AR, FOXA1 and MYC as potential drivers of AP activation. DNA methylation is a likely mechanism for AP activation during tumour progression and lineage plasticity. Our data suggest that prostate tumours activate APs to magnify the transcriptional impact of tumour drivers, including AR and MYC.

DNA Hypomethylation Underlies Epigenetic Swapping between AGO1 and AGO1-V2 Isoforms in Tumors

Human tumors progress in part by accumulating epigenetic alterations, which include gains and losses of DNA methylation in different parts of the cancer cell genome. Recent work has revealed a link between these two opposite alterations by showing that DNA hypomethylation in tumors can induce the expression of transcripts that overlap downstream gene promoters and thereby induce their hypermethylation. Preliminary in silico evidence prompted us to investigate if this mechanism applies to the locus harboring AGO1, a gene that plays a central role in miRNA biogenesis and RNA interference. Inspection of public RNA-Seq datasets and RT-qPCR experiments show that an alternative transcript starting 13.4 kb upstream of AGO1 (AGO1-V2) is expressed specifically in testicular germ cells, and becomes aberrantly activated in different types of tumors, particularly in tumors of the esophagus, stomach, and lung. This expression pattern classifies AGO1-V2 into the group of "Cancer-Germline" (CG) genes. Analysis of transcriptomic and methylomic datasets provided evidence that transcriptional activation of AGO1-V2 depends on DNA demethylation of its promoter region. Western blot experiments revealed that AGO1-V2 encodes a shortened isoform of AGO1, corresponding to a truncation of 75 aa in the N-terminal domain, and which we therefore referred to as "∆NAGO1". Interestingly, significant correlations between hypomethylation/activation of AGO1-V2 and hypermethylation/repression of AGO1 were observed upon examination of tumor cell lines and tissue datasets. Overall, our study reveals the existence of a process of interdependent epigenetic alterations in the AGO1 locus, which promotes swapping between two AGO1 protein-coding mRNA isoforms in tumors.

TAD hierarchy restricts poised LTR activation and loss of TAD hierarchy promotes LTR co-option in cancer

Transposable elements (TEs) are abundant in the human genome, and they provide the sources for genetic and functional diversity. The regulation of TEs expression and their functional consequences in physiological conditions and cancer development remain to be fully elucidated. Previous studies suggested TEs are repressed by DNA methylation and chromatin modifications. The effect of 3D chromatin topology on TE regulation remains elusive. Here, by integrating transcriptome and 3D genome architecture studies, we showed that haploinsufficient loss of NIPBL selectively activates alternative promoters at the long terminal repeats (LTRs) of the TE subclasses. This activation occurs through the reorganization of topologically associating domain (TAD) hierarchical structures and recruitment of proximal enhancers. These observations indicate that TAD hierarchy restricts transcriptional activation of LTRs that already possess open chromatin features. In cancer, perturbation of the hierarchical chromatin topology can lead to co-option of LTRs as functional alternative promoters in a context-dependent manner and drive aberrant transcriptional activation of novel oncogenes and other divergent transcripts. These data uncovered a new layer of regulatory mechanism of TE expression beyond DNA and chromatin modification in human genome. They also posit the TAD hierarchy dysregulation as a novel mechanism for alternative promoter-mediated oncogene activation and transcriptional diversity in cancer, which may be exploited therapeutically.

CYP19A1 Expression Is Controlled by mRNA Stability of the Upstream Transcription Factor AP-2γ in Placental JEG3 Cells

CYP19A1 encodes aromatase, which converts testosterone to estrogen, and is induced during placental maturation. To elucidate the molecular mechanism underlying this function, histone methylation was analyzed using the placental cytotrophoblast cell line, JEG3. Treatment of JEG3 cells with 3-deazaneplanocin A, an inhibitor of several methyltransferases, resulted in increased CYP19A1 expression, accompanied by removal of the repressive mark H3K27me3 from the CYP19A1 promoter. However, this increase was not observed in cells treated with GSK126, another specific inhibitor for H3K27me3 methylation. Expression of TFAP2C, which encodes AP-2γ, a transcription factor that regulates CYP19A1, was also elevated on 3-deazaneplanocin A treatment. Interestingly, TFAP2C messenger RNA (mRNA) was readily degraded in JEG3 cells but protected from degradation in the presence of 3-deazaneplanocin A. TFAP2C mRNA contained N6-methyladenosines, which were reduced on drug treatment. These observations indicate that the TFAP2C mRNA undergoes adenosine methylation and rapid degradation, whereas 3-deazaneplanocin A suppresses methylation, resulting in an increase in AP-2γ levels. We conclude that the increase in AP-2γ expression via stabilization of the TFAP2C mRNA is likely to underlie the increased CYP19A1 expression.

Evolution of the enhancer-rich regulatory region of the gene for the cell-type specific transcription factor POU1F1

Precise spatio-temporal expression of genes in organogenesis is regulated by the coordinated interplay of DNA elements such as promoter and enhancers present in the regulatory region of a given locus. POU1F1 transcription factor plays a crucial role in the development of somatotrophs, lactotrophs and thyrotrophs in the anterior pituitary gland, and in maintaining high expression of growth hormone, prolactin and TSH. In mouse, expression of POU1F1 is controlled by a region fenced by two CTCF sites, containing 5 upstream enhancer elements, designated E-A (5' to 3'). Elements C, B and A correspond to elements shown previously to play a role in pituitary development and hormonal expression; functional roles for elements E and D have not been reported. We performed comparative sequence analysis of this regulatory region and discovered that three elements, B, C and E, are present in all vertebrate groups except Agnatha. One very long (>2 kb) element (A) is unique to mammals suggesting a specific change in regulation of the gene in this group. Using DNA accessibility assay (ATAC-seq) we showed that conserved elements in anterior pituitary of four non-mammals are open, suggesting functionality as regulatory elements. We showed that, in many non-mammalian vertebrates, an additional upstream exon closely follows element E, leading to alternatively spliced transcripts. Here, element E functions as an alternative promoter, but in mammals this feature is lost, suggesting conversion of alternative promoter to enhancer. Our work shows that regulation of POU1F1 changed markedly during the course of vertebrate evolution, use of a low number of enhancer elements combined with alternative promoters in non-mammalian vertebrates being replaced by use of a unique combination of regulatory units in mammals. Most importantly, our work suggests that evolutionary conversion of alternate promoter to enhancer could be one of the evolutionary mechanisms of enhancer birth.

Detecting differential transcript usage in complex diseases with SPIT

Differential transcript usage (DTU) plays a crucial role in determining how gene expression differs among cells, tissues, and developmental stages, contributing to the complexity and diversity of biological systems. In abnormal cells, it can also lead to deficiencies in protein function and underpin disease pathogenesis. Analyzing DTU via RNA sequencing (RNA-seq) data is vital, but the genetic heterogeneity in populations with complex diseases presents an intricate challenge due to diverse causal events and undetermined subtypes. Although the majority of common diseases in humans are categorized as complex, state-of-the-art DTU analysis methods often overlook this heterogeneity in their models. We therefore developed SPIT, a statistical tool that identifies predominant subgroups in transcript usage within a population along with their distinctive sets of DTU events. This study provides comprehensive assessments of SPIT's methodology and applies it to analyze brain samples from individuals with schizophrenia, revealing previously unreported DTU events in six candidate genes.

Genome-wide analyses reveal the regulatory roles of DNA methylation-regulated alternative promoter transcripts in breast cancer

A certain proportion of genes are regulated by multiple, distinct promoters, revealing a dynamic landscape of the cancer transcriptome. However, the contribution of alternative promoters (APs) in breast cancer (BRCA) remains largely unexplored. Here, we identified 3654 genes with multiple promoters in BRCA patients, and 53 of them could generate distinct AP transcripts that are dysregulated and prognosis-related in BRCA, namely prognosis-related dysregulated AP (prdeAP) transcripts. Interestingly, when we searched for the genomic signatures of these prdeAP genes, we found that the promoter regions of 92% of the prdeAP genes were enriched with abundant DNA methylation signals. Through further bioinformatic analysis and experimental validation, we showed that AP selections of TANK, UNKL, CCL28, and MAP1LC3A were regulated by DNA methylation upon their corresponding promoter regions. Functionally, by overexpressing AP variants of TANK, we found that TANK|55731 could dramatically suppress MDA-MB-231 cell proliferation and migration. Meanwhile, pan-cancer survival analyses suggested that AP variants of TANK provided more accurate prognostic predictive ability than TANK gene in a variety of tumor types, including BRCA. Together, by uncovering the DNA methylation-regulated AP transcripts with tumor prognostic features, our work revealed a novel layer of regulators in BRCA progression and provided potential targets that served as effective biomarkers for anti-BRCA treatment.

Integrative Epigenetic and Molecular Analysis Reveals a Novel Promoter for a New Isoform of the Transcription Factor TEAD4

TEAD4 is a transcription factor that plays a crucial role in the Hippo pathway by regulating the expression of genes related to proliferation and apoptosis. It is also involved in the maintenance and differentiation of the trophectoderm during pre- and post-implantation embryonic development. An alternative promoter for the TEAD4 gene was identified through epigenetic profile analysis, and a new transcript from the intronic region of TEAD4 was discovered using the 5'RACE method. The transcript of the novel promoter encodes a TEAD4 isoform (TEAD4-ΔN) that lacks the DNA-binding domain but retains the C-terminal protein-protein interaction domain. Gene expression studies, including end-point PCR and Western blotting, showed that full-length TEAD4 was present in all investigated tissues. However, TEAD4-ΔN was only detectable in certain cell types. The TEAD4-ΔN promoter is conserved throughout evolution and demonstrates transcriptional activity in transient-expression experiments. Our study reveals that TEAD4 interacts with the alternative promoter and increases the expression of the truncated isoform. DNA methylation plays a crucial function in the restricted expression of the TEAD4-ΔN isoform in specific tissues, including the umbilical cord and the placenta. The data presented indicate that the DNA-methylation status of the TEAD4-ΔN promoter plays a critical role in regulating organ size, cancer development, and placenta differentiation.

Transcript CD81-215 may be a long noncoding RNA of stromal origin with tumor-promoting role in colon cancer

The role of tetraspanin CD81 in malignant transformation is best studied in colorectal cancer, and it appears that other transcripts beside the fully coding mRNA may also be dysregulated in malignant cells. Recent data from a comprehensive pan-cancer transcriptome analysis demonstrated differential activity of two alternative CD81 gene promoters in malignant versus nonmalignant gut mucosa. The promoter active in gut mucosa gives rise to transcripts CD81-203 and CD81-213, while the promoter active in colon and rectal cancer gives rise to transcripts CD81-205 and CD81-215. Our study aimed to explore the biomarker potential of the transcripts from the alternative CD81 gene promoters in colon cancer, as well as to investigate their structure and potential function using in silico tools. The analysis of the transcripts' expression in several colon cell lines cultivated in 2D and 3D and a set of colon cancer and healthy gut mucosa samples by qPCR and RNA sequencing suggested their low expression and stromal origin. Expression patterns in tumor and nontumor tissue along with in silico data suppose that the transcript CD81-215 may be a noncoding RNA of stromal origin with possible involvement in signaling related to malignant transformation.

Epigenetic Landscape and Therapeutic Implication of Gene Isoforms of Doublecortin-Like Kinase 1 for Cancer Stem Cells

While significant strides have been made in understanding cancer biology, the enhancement in patient survival is limited, underscoring the urgency for innovative strategies. Epigenetic modifications characterized by hereditary shifts in gene expression without changes to the DNA sequence play a critical role in producing alternative gene isoforms. When these processes go awry, they influence cancer onset, growth, spread, and cancer stemness. In this review, we delve into the epigenetic and isoform nuances of the protein kinase, doublecortin-like kinase 1 (DCLK1). Recognized as a hallmark of tumor stemness, DCLK1 plays a pivotal role in tumorigenesis, and DCLK1 isoforms, shaped by alternative promoter usage and splicing, can reveal potential therapeutic touchpoints. Our discussion centers on recent findings pertaining to the specific functions of DCLK1 isoforms and the prevailing understanding of its epigenetic regulation via its two distinct promoters. It is noteworthy that all DCLK1 isoforms retain their kinase domain, suggesting that their unique functionalities arise from non-kinase mechanisms. Consequently, our research has pivoted to drugs that specifically influence the epigenetic generation of these DCLK1 isoforms. We posit that a combined therapeutic approach, harnessing both the epigenetic regulators of specific DCLK1 isoforms and DCLK1-targeted drugs, may prove more effective than therapies that solely target DCLK1.

Alternative Transcripts Diversify Genome Function for Phenome Relevance to Health and Diseases

Manipulation using alternative exon splicing (AES), alternative transcription start (ATS), and alternative polyadenylation (APA) sites are key to transcript diversity underlying health and disease. All three are pervasive in organisms, present in at least 50% of human protein-coding genes. In fact, ATS and APA site use has the highest impact on protein identity, with their ability to alter which first and last exons are utilized as well as impacting stability and translation efficiency. These RNA variants have been shown to be highly specific, both in tissue type and stage, with demonstrated importance to cell proliferation, differentiation and the transition from fetal to adult cells. While alternative exon splicing has a limited effect on protein identity, its ubiquity highlights the importance of these minor alterations, which can alter other features such as localization. The three processes are also highly interwoven, with overlapping, complementary, and competing factors, RNA polymerase II and its CTD (C-terminal domain) chief among them. Their role in development means dysregulation leads to a wide variety of disorders and cancers, with some forms of disease disproportionately affected by specific mechanisms (AES, ATS, or APA). Challenges associated with the genome-wide profiling of RNA variants and their potential solutions are also discussed in this review.

The Role of NCS1 in Immunotherapy and Prognosis of Human Cancer

The Neural Calcium Sensor1 (NCS1) is a crucial protein that binds to Ca2+ and is believed to play a role in regulating tumor invasion and cell proliferation. However, the role of NCS1 in immune infiltration and cancer prognosis is still unknown. Our study aimed to explore the expression profile, immune infiltration pattern, prognostic value, biological function, and potential compounds targeting NCS1 using public databases. High expression of NCS1 was detected by immune histochemical staining in LIHC (Liver hepatocellular carcinoma), BRCA (Breast invasive carcinoma), KIRC (Kidney renal clear cell carcinoma), and SKCM (Skin Cutaneous Melanoma). The expression of NCS1 in cancer was determined by TCGA (The Cancer Genome Atlas Program), GTEx (The Genotype-Tissue Expression), the Kaplan-Meier plotter, GEO (Gene Expression Omnibus), GEPIA2.0 (Gene Expression Profiling Interactive Analysis 2.0), HPA (The Human Protein Atlas), UALCAN, TIMER2.0, TISIDB, Metascape, Drugbank, chEMBL, and ICSDB databases. NCS1 has genomic mutations as well as aberrant DNA methylation in multiple cancers compared to normal tissues. Also, NCS1 was significantly different in the immune microenvironment, tumor mutational burden (TMB), microsatellite instability (MSI), and immune infiltrate-associated cells in different cancers, which could be used for the typing of immune and molecular subtypes of cancer and the presence of immune checkpoint resistance in several cancers. Univariate regression analysis, multivariate regression analysis, and gene enrichment analysis to construct prognostic models revealed that NCS1 is involved in immune regulation and can be used as a prognostic biomarker for SKCM, LIHC, BRCA, COAD, and KIRC. These results provide clues from a bioinformatic perspective and highlight the importance of NCS1 in a variety of cancers.

The single-cell landscape of alternative transcription start sites of diabetic retina

More than half of mammalian protein-coding genes have multiple transcription start sites. Alternative transcription start site (TSS) modulate mRNA stability, localization, and translation efficiency on post-transcription level, and even generate novel protein isoforms. However, differential TSS usage among cell types in healthy and diabetic retina remains poorly characterized. In this study, by using 5'-tag-based single-cell RNA sequencing, we identified cell type-specific alternative TSS events and key transcription factors for each of retinal cell types. We observed that lengthening of 5'- UTRs in retinal cell types are enriched for multiple RNA binding protein binding sites, including splicing regulators Rbfox1/2/3 and Nova1. Furthermore, by comparing TSS expression between healthy and diabetic retina, we identified elevated apoptosis signal in Müller glia and microglia, which can be served as a putative early sign of diabetic retinopathy. By measuring 5'UTR isoforms in retinal single-cell dataset, our work provides a comprehensive panorama of alternative TSS and its potential consequence related to post-transcriptional regulation. We anticipate our assay can not only provide insights into cellular heterogeneity driven by transcriptional initiation, but also open up the perspectives for identification of novel diagnostic indexes for diabetic retinopathy.

Detecting differential transcript usage in complex diseases with SPIT

Differential transcript usage (DTU) plays a crucial role in determining how gene expression differs among cells, tissues, and different developmental stages, thereby contributing to the complexity and diversity of biological systems. In abnormal cells, it can also lead to deficiencies in protein function, potentially leading to pathogenesis of diseases. Detecting such events for single-gene genetic traits is relatively uncomplicated; however, the heterogeneity of populations with complex diseases presents an intricate challenge due to the presence of diverse causal events and undetermined subtypes. SPIT is the first statistical tool that quantifies the heterogeneity in transcript usage within a population and identifies predominant subgroups along with their distinctive sets of DTU events. We provide comprehensive assessments of SPIT's methodology in both single-gene and complex traits and report the results of applying SPIT to analyze brain samples from individuals with schizophrenia. Our analysis reveals previously unreported DTU events in six candidate genes.

Integrative genomic analyses of promoter G-quadruplexes reveal their selective constraint and association with gene activation

G-quadruplexes (G4s) regulate DNA replication and gene transcription, and are enriched in promoters without fully appreciated functional relevance. Here we show high selection pressure on putative G4 (pG4) forming sequences in promoters through investigating genetic and genomic data. Analyses of 76,156 whole-genome sequences reveal that G-tracts and connecting loops in promoter pG4s display lower or higher allele frequencies, respectively, than pG4-flanking regions, and central guanines (Gs) in G-tracts show higher selection pressure than other Gs. Additionally, pG4-promoters produce over 72.4% of transcripts, and promoter G4-containing genes are expressed at relatively high levels. Most genes repressed by TMPyP4, a G4-ligand, regulate epigenetic processes, and promoter G4s are enriched with gene activation histone marks, chromatin remodeler and transcription factor binding sites. Consistently, cis-expression quantitative trait loci (cis-eQTLs) are enriched in promoter pG4s and their G-tracts. Overall, our study demonstrates selective constraint of promoter G4s and reinforces their stimulative role in gene expression.

Pre- and Post-Transcriptional Control of HBV Gene Expression: The Road Traveled towards the New Paradigm of HBx, Its Isoforms, and Their Diverse Functions

Hepatitis B virus (HBV) is an enveloped DNA human virus belonging to the Hepadnaviridae family. Perhaps its main distinguishable characteristic is the replication of its genome through a reverse transcription process. The HBV circular genome encodes only four overlapping reading frames, encoding for the main canonical proteins named core, P, surface, and X (or HBx protein). However, pre- and post-transcriptional gene regulation diversifies the full HBV proteome into diverse isoform proteins. In line with this, hepatitis B virus X protein (HBx) is a viral multifunctional and regulatory protein of 16.5 kDa, whose canonical reading frame presents two phylogenetically conserved internal in-frame translational initiation codons, and which results as well in the expression of two divergent N-terminal smaller isoforms of 8.6 and 5.8 kDa, during translation. The canonical HBx, as well as the smaller isoform proteins, displays different roles during viral replication and subcellular localizations. In this article, we reviewed the different mechanisms of pre- and post-transcriptional regulation of protein expression that take place during viral replication. We also investigated all the past and recent evidence about HBV HBx gene regulation and its divergent N-terminal isoform proteins. Evidence has been collected for over 30 years. The accumulated evidence simply strengthens the concept of a new paradigm of the canonical HBx, and its smaller divergent N-terminal isoform proteins, not only during viral replication, but also throughout cell pathogenesis.

Reciprocal regulation of flower induction by ELF3α and ELF3β generated via alternative promoter usage

Flowering is critical for sexual reproduction and fruit production. Several pear (Pyrus sp.) varieties produce few flower buds, but the underlying mechanisms are unknown. The circadian clock regulator EARLY FLOWERING3 (ELF3) serves as a scaffold protein in the evening complex that controls flowering. Here, we report that the absence of a 58-bp sequence in the 2nd intron of PbELF3 is genetically associated with the production of fewer flower buds in pear. From rapid amplification of cDNA ends sequencing results, we identified a short, previously unknown transcript from the PbELF3 locus, which we termed PbELF3β, whose transcript level was significantly lower in pear cultivars that lacked the 58-bp region. The heterologous expression of PbELF3β in Arabidopsis (Arabidopsis thaliana) accelerated flowering, whereas the heterologous expression of the full-length transcript PbELF3α caused late flowering. Notably, ELF3β was functionally conserved in other plants. Deletion of the 2nd intron reduced AtELF3β expression and caused delayed flowering time in Arabidopsis. AtELF3β physically interacted with AtELF3α, disrupting the formation of the evening complex and consequently releasing its repression of flower induction genes such as GIGANTEA (GI). AtELF3β had no effect in the absence of AtELF3α, supporting the idea that AtELF3β promotes flower induction by blocking AtELF3α function. Our findings show that alternative promoter usage at the ELF3 locus allows plants to fine-tune flower induction.

CAGE sequencing reveals CFTR-dependent dysregulation of type I IFN signaling in activated cystic fibrosis macrophages

An intense, nonresolving airway inflammatory response leads to destructive lung disease in cystic fibrosis (CF). Dysregulation of macrophage immune function may be a key facet governing the progression of CF lung disease, but the underlying mechanisms are not fully understood. We used 5' end centered transcriptome sequencing to profile P. aeruginosa LPS-activated human CF macrophages, showing that CF and non-CF macrophages deploy substantially distinct transcriptional programs at baseline and following activation. This includes a significantly blunted type I IFN signaling response in activated patient cells relative to healthy controls that was reversible upon in vitro treatment with CFTR modulators in patient cells and by CRISPR-Cas9 gene editing to correct the F508del mutation in patient-derived iPSC macrophages. These findings illustrate a previously unidentified immune defect in human CF macrophages that is CFTR dependent and reversible with CFTR modulators, thus providing new avenues in the search for effective anti-inflammatory interventions in CF.

TP53BP1, a dual-coding gene, uses promoter switching and translational reinitiation to express a smORF protein

The complexity of the metazoan proteome is significantly increased by the expression of small proteins (<100 aa) derived from smORFs within lncRNAs, uORFs, 3' UTRs and, reading frames overlapping the CDS. These smORF encoded proteins (SEPs) have diverse roles, ranging from the regulation of cellular physiological to essential developmental functions. We report the characterization of a new member of this protein family, SEP53BP1, derived from a small internal ORF that overlaps the CDS encoding 53BP1. Its expression is coupled to the utilization of an alternative, cell-type specific promoter coupled to translational reinitiation events mediated by a uORF in the alternative 5' TL of the mRNA. This uORF-mediated reinitiation at an internal ORF is also observed in zebrafish. Interactome studies indicate that the human SEP53BP1 associates with components of the protein turnover pathway including the proteasome, and the TRiC/CCT chaperonin complex, suggesting that it may play a role in cellular proteostasis.

Condition-specific 3' mRNA isoform half-lives and stability elements in yeast

Alternative polyadenylation generates numerous 3' mRNA isoforms that can differ in their stability, structure, and function. These isoforms can be used to map mRNA stabilizing and destabilizing elements within 3' untranslated regions (3'UTRs). Here, we examine how environmental conditions affect 3' mRNA isoform turnover and structure in yeast cells on a transcriptome scale. Isoform stability broadly increases when cells grow more slowly, with relative half-lives of most isoforms being well correlated across multiple conditions. Surprisingly, dimethyl sulfate probing reveals that individual 3' isoforms have similar structures across different conditions, in contrast to the extensive structural differences that can exist between closely related isoforms in an individual condition. Unexpectedly, most mRNA stabilizing and destabilizing elements function only in a single growth condition. The genes associated with some classes of condition-specific stability elements are enriched for different functional categories, suggesting that regulated mRNA stability might contribute to adaptation to different growth environments. Condition-specific stability elements do not result in corresponding condition-specific changes in steady-state mRNA isoform levels. This observation is consistent with a compensatory mechanism between polyadenylation and stability, and it suggests that condition-specific mRNA stability elements might largely reflect condition-specific regulation of mRNA 3' end formation.

Genome-Wide Mapping of 5' Isoforms with 5'-Seq

The transcriptome is far more complex than previously assumed. Transcripts from the same gene can differ in terms of transcription start site, transcription end site, or pattern of splicing, and growing evidence supports the functional importance of these distinct transcript isoforms. Easily identifying these isoforms experimentally via library construction and high-throughput sequencing is crucial. Current library construction methods for identifying transcription start sites (5' transcript isoforms) involve large number of steps and (expensive) reagents, utilization of cDNA intermediates for adapter ligation, and are less suitable for studying low-abundance isoforms. Here, I describe a quick protocol for the generation of sequencing libraries to define capped 5' isoforms (5'-Seq) of various abundances in yeast and suggest a 5' isoform data analysis pipeline. The protocol relies on the utilization of a dephosphorylation-decapping method (oligo-capping) to generate a sequencing library from mRNA fragments and is a simplification of previously published 5' isoform protocols in terms of the handling steps, time, and cost. This method is exemplified using Saccharomyces cerevisiae mRNA, but it can be applied to various cellular conditions to study the effects of 5' transcript isoforms on transcriptional and/or translational regulation. © 2023 Wiley Periodicals LLC. Basic Protocol: Construction of a DNA sequencing library from capped 5' isoforms Support Protocol: Sequencing data analysis.

Sex chromosome mechanisms in cardiac development and disease

Many human diseases, including cardiovascular disease, show differences between men and women in pathology and treatment outcomes. In the case of cardiac disease, sex differences are exemplified by differences in the frequency of specific types of congenital and adult-onset heart disease. Clinical studies have suggested that gonadal hormones are a factor in sex bias. However, recent research has shown that gene and protein networks under non-hormonal control also account for cardiac sex differences. In this review, we describe the sex chromosome pathways that lead to sex differences in the development and function of the heart and highlight how these findings affect future care and treatment of cardiac disease.

Identification and characterization of the Doublesex gene and its mRNA isoforms in the brine shrimp Artemia franciscana

Doublesex (DSX) proteins are members of the Doublesex/mab-3-related (DMRT) protein family and play crucial roles in sex determination and differentiation among the animal kingdom. In the present study, we identified two Doublesex (Dsx)-like mRNA isoforms in the brine shrimp Artemia franciscana (Kellogg 1906), which are generated by the combination of alternative promoters, alternative splicing and alternative polyadenylation. The two transcripts exhibited sex-biased enrichment, which we termed AfrDsxM and AfrDsxF. They share a common region which encodes an identical N-terminal DNA-binding (DM) domain. RT-qPCR analyses showed that AfrDsxM is dominantly expressed in male Artemia while AfrDsxF is specifically expressed in females. Expression levels of both isoforms increased along with the developmental stages of their respective sexes. RNA interference with dsRNA showed that the knockdown of AfrDsxM in male larvae led to the appearance of female traits including an ovary-like structure in the original male reproductive system and an elevated expression of vitellogenin. However, silencing of AfrDsxF induced no clear phenotypic change in female Artemia. These results indicated that the male AfrDSXM may act as inhibiting regulator upon the default female developmental mode in Artemia. Furthermore, electrophoretic mobility shift assay analyses revealed that the unique DM domain of AfrDSXs can specifically bind to promoter segments of potential downstream target genes like AfrVtg. These data show that AfrDSXs play crucial roles in regulating sexual development in Artemia, and further provide insight into the evolution of sex determination/differentiation in sexual organisms.

Cloning and characterization of nitrate reductase gene in kelp Saccharina japonica (Laminariales, Phaeophyta)

BACKGROUND

Brown macroalgae dominate temperate coastal ecosystems, and their productivity is typically limited by nitrate availability. As an economically important kelp, Saccharina japonica is the most productive farmed seaweed and needs to be supplemented with sufficient nitrate throughout the cultivation process. However, molecular characterization of genes involved in nitrogen assimilation has not been conducted in brown macroalgae.

RESULTS

Here, we described the identification of the nitrate reductase (NR) gene from S. japonica (SjNR). Using two different cloning methods for SjNR, i.e. rapid amplification of cDNA ends (RACE) and cDNA cloning alone, a single fragment was obtained respectively. According to results of sequence analysis between these two fragments, the tentative coding sequence in two clones, SjNR-L and SjNR-S, were suggested to represent two transcripts of the single copy SjNR, and the ATG of SjNR-S was located inside the third exon of SjNR-L. In the 5' upstream sequence of each transcript, promoter core elements, response elements, especially multiple N response elements which occurred in microalgal NR, were all predicted. Further sequence analysis revealed that both transcripts encoded all five domains conserved in eukaryotic plant NRs. RT-qPCR results showed that the transcription level of SjNR in juvenile sporophytes could be significantly induced by nitrate and inhibited by ammonium, which was in line with plant NRs. The recombinant SjNR-L and SjNR-S were all proved to have NR activity, suggesting that the single-copy gene SjNR might be regulated on transcription level based on alternative promoters and multiple transcriptional start sites. Moreover, both NADH and NADPH were found to be able to act as electron donors for SjNR alone, which is the first confirmation that brown algal NR has a NAD(P)H-bispecific form.

CONCLUSION

These results will provide a scientific basis for understanding the N demand of kelp in various stages of cultivation and evaluating the environmental remediation potential of kelp in eutrophic sea areas.

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.

NF-κB regulates the expression of STING via alternative promoter usage

AIMS

Stimulator of interferon genes (STING) is a transmembrane protein in endoplasmic reticulum and plays crucial roles in autophagy, antiviral and anti-tumor responses. However, there are few studies on the transcriptional regulation mechanism of STING.

MAIN METHODS

The 5' RACE experiment was used to determine the location of STING promoters. Luciferase reporting assay confirmed the activity and core region of STING internal promoter. Site-directed mutagenesis confirmed that NF-κB regulates the activity of STING promoters. The regulation of NF-κB on STING was investigated by real-time quantitative PCR, western blot, chromatin immunoprecipitation assay and lipopolysaccharide (LPS) inflammatory cell model.

KEY FINDINGS

There was also a transcription start site at the 17 bp sequence upstream of STING second exon. STING-285 was the core region of the internal promoter. After NF-κB binding site mutation, the activity of STING internal promoter decreased significantly. In addition, we found that NF-κB can bind to the promoter region of wild-type STING. Overexpression of NF-κB significantly increased the activity of STING internal promoter and wild-type promoter, while knockdown of endogenous NF-κB significantly inhibited the activity of STING promoters. The binding of NF-κB to STING promoters in vivo were confirmed by chromatin immunoprecipitation assay. Meanwhile, we stimulated HeLa cells with LPS to activate the NF-κB pathway and found that STING expression was up-regulated.

SIGNIFICANCE

These results suggest that transcription factor NF-κB positively regulates the expression of STING via alternative promoter usage. This provides a new basis and potential drug targets for the clinical treatment of STING related diseases.

Applications of deep learning in understanding gene regulation

Gene regulation is a central topic in cell biology. Advances in omics technologies and the accumulation of omics data have provided better opportunities for gene regulation studies than ever before. For this reason deep learning, as a data-driven predictive modeling approach, has been successfully applied to this field during the past decade. In this article, we aim to give a brief yet comprehensive overview of representative deep-learning methods for gene regulation. Specifically, we discuss and compare the design principles and datasets used by each method, creating a reference for researchers who wish to replicate or improve existing methods. We also discuss the common problems of existing approaches and prospectively introduce the emerging deep-learning paradigms that will potentially alleviate them. We hope that this article will provide a rich and up-to-date resource and shed light on future research directions in this area.

A Guide to MethylationToActivity: A Deep Learning Framework That Reveals Promoter Activity Landscapes from DNA Methylomes in Individual Tumors

Genome-wide DNA methylomes have contributed greatly to tumor detection and subclassification. However, interpreting the biological impact of the DNA methylome at the individual gene level remains a challenge. MethylationToActivity (M2A) is a pipeline that uses convolutional neural networks to infer H3K4me3 and H3K27ac enrichment from DNA methylomes and thus infer promoter activity. It was shown to be highly accurate and robust in revealing promoter activity landscapes in various pediatric and adult cancers. The following will present a user-friendly guide through the model pipeline.

The potential of cold-shock promoters for the expression of recombinant proteins in microbes and mammalian cells

BACKGROUND

Low-temperature expression of recombinant proteins may be advantageous to support their proper folding and preserve bioactivity. The generation of expression vectors regulated under cold conditions can improve the expression of some target proteins that are difficult to express in different expression systems. The cspA encodes the major cold-shock protein from Escherichia coli (CspA). The promoter of cspA has been widely used to develop cold shock-inducible expression platforms in E. coli. Moreover, it is often necessary to employ expression systems other than bacteria, particularly when recombinant proteins require complex post-translational modifications. Currently, there are no commercial platforms available for expressing target genes by cold shock in eukaryotic cells. Consequently, genetic elements that respond to cold shock offer the possibility of developing novel cold-inducible expression platforms, particularly suitable for yeasts, and mammalian cells.

CONCLUSIONS

This review covers the importance of the cellular response to low temperatures and the prospective use of cold-sensitive promoters to direct the expression of recombinant proteins. This concept may contribute to renewing interest in applying white technologies to produce recombinant proteins that are difficult to express.

Mediators of Capillary-to-Venule Conversion in the Chronic Inflammatory Skin Disease Psoriasis

Psoriasis is a chronic inflammatory skin disease characterized by epidermal hyperplasia and hyperkeratosis, immune cell infiltration and vascular remodeling. Despite the emerging recognition of vascular normalization as a potential strategy for managing psoriasis, an in-depth delineation of the remodeled dermal vasculature has been missing. In this study, we exploited 5' single-cell RNA sequencing to investigate the transcriptomic alterations in different subpopulations of blood vascular and lymphatic endothelial cells directly isolated from psoriatic and healthy human skin. Individual subtypes of endothelial cells underwent specific molecular repatterning associated with cell adhesion and extracellular matrix organization. Blood capillaries, in particular, showed upregulation of the melanoma cell adhesion molecule as well as its binding partners and adopted postcapillary venule‒like characteristics during chronic inflammation that are more permissive to leukocyte transmigration. We also identified psoriasis-specific interactions between cis-regulatory enhancers and promoters for each endothelial cell subtype, revealing the dysregulated gene regulatory networks in psoriasis. Together, our results provide more insights into the specific transcriptional responses and epigenetic signatures of endothelial cells lining different vessel compartments in chronic skin inflammation.

Gene product diversity: adaptive or not?

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

Alternative Transcription Start Site Usage and Functional Implications in Pathogenic Fungi

Pathogenic fungi require delicate gene regulation mechanisms to adapt to diverse living environments and escape host immune systems. Recent advances in sequencing technology have exposed the complexity of the fungal genome, thus allowing the gradual disentanglement of multiple layers of gene expression control. Alternative transcription start site (aTSS) usage, previously reported to be prominent in mammals and to play important roles in physiopathology, is also present in fungi to fine-tune gene expression. Depending on the alteration in their sequences, RNA isoforms arising from aTSSs acquire different characteristics that significantly alter their stability and translational capacity as well as the properties and biologic functions of the resulting proteins. Disrupted control of aTSS usage has been reported to severely impair growth, virulence, and the infectious capacity of pathogenic fungi. Here, we discuss principle concepts, mechanisms, and the functional implication of aTSS usage in fungi.

An Epigenetic Role of Mitochondria in Cancer

Mitochondria are not only the main energy supplier but are also the cell metabolic center regulating multiple key metaborates that play pivotal roles in epigenetics regulation. These metabolites include acetyl-CoA, α-ketoglutarate (α-KG), S-adenosyl methionine (SAM), NAD⁺, and O-linked beta-N-acetylglucosamine (O-GlcNAc), which are the main substrates for DNA methylation and histone post-translation modifications, essential for gene transcriptional regulation and cell fate determination. Tumorigenesis is attributed to many factors, including gene mutations and tumor microenvironment. Mitochondria and epigenetics play essential roles in tumor initiation, evolution, metastasis, and recurrence. Targeting mitochondrial metabolism and epigenetics are promising therapeutic strategies for tumor treatment. In this review, we summarize the roles of mitochondria in key metabolites required for epigenetics modification and in cell fate regulation and discuss the current strategy in cancer therapies via targeting epigenetic modifiers and related enzymes in metabolic regulation. This review is an important contribution to the understanding of the current metabolic-epigenetic-tumorigenesis concept.

Uncovering the impacts of alternative splicing on the proteome with current omics techniques

The high-throughput sequencing of cellular RNAs has underscored a broad effect of isoform diversification through alternative splicing on the transcriptome. Moreover, the differential production of transcript isoforms from gene loci has been recognized as a critical mechanism in cell differentiation, organismal development, and disease. Yet, the extent of the impact of alternative splicing on protein production and cellular function remains a matter of debate. Multiple experimental and computational approaches have been developed in recent years to address this question. These studies have unveiled how molecular changes at different steps in the RNA processing pathway can lead to differences in protein production and have functional effects. New and emerging experimental technologies open exciting new opportunities to develop new methods to fully establish the connection between messenger RNA expression and protein production and to further investigate how RNA variation impacts the proteome and cell function. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing Translation > Regulation RNA Evolution and Genomics > Computational Analyses of RNA.

Epigenetic promoter alterations in GI tumour immune-editing and resistance to immune checkpoint inhibition

OBJECTIVES

Epigenomic alterations in cancer interact with the immune microenvironment to dictate tumour evolution and therapeutic response. We aimed to study the regulation of the tumour immune microenvironment through epigenetic alternate promoter use in gastric cancer and to expand our findings to other gastrointestinal tumours.

DESIGN

Alternate promoter burden (APB) was quantified using a novel bioinformatic algorithm (proActiv) to infer promoter activity from short-read RNA sequencing and samples categorised into APB_high, APB_int and APB_low. Single-cell RNA sequencing was performed to analyse the intratumour immune microenvironment. A humanised mouse cancer in vivo model was used to explore dynamic temporal interactions between tumour kinetics, alternate promoter usage and the human immune system. Multiple cohorts of gastrointestinal tumours treated with immunotherapy were assessed for correlation between APB and treatment outcomes.

RESULTS

APB_high gastric cancer tumours expressed decreased levels of T-cell cytolytic activity and exhibited signatures of immune depletion. Single-cell RNAsequencing analysis confirmed distinct immunological populations and lower T-cell proportions in APB_high tumours. Functional in vivo studies using 'humanised mice' harbouring an active human immune system revealed distinct temporal relationships between APB and tumour growth, with APB_high tumours having almost no human T-cell infiltration. Analysis of immunotherapy-treated patients with GI cancer confirmed resistance of APB_high tumours to immune checkpoint inhibition. APB_high gastric cancer exhibited significantly poorer progression-free survival compared with APB_low (median 55 days vs 121 days, HR 0.40, 95% CI 0.18 to 0.93, p=0.032).

CONCLUSION

These findings demonstrate an association between alternate promoter use and the tumour microenvironment, leading to immune evasion and immunotherapy resistance.

Differential regulation of alternative promoters emerges from unified kinetics of enhancer-promoter interaction

Many eukaryotic genes contain alternative promoters with distinct expression patterns. How these promoters are differentially regulated remains elusive. Here, we apply single-molecule imaging to quantify the transcriptional regulation of two alternative promoters (P1 and P2) of the Bicoid (Bcd) target gene hunchback in syncytial blastoderm Drosophila embryos. Contrary to the previous notion that Bcd only activates P2, we find that Bcd activates both promoters via the same two enhancers. P1 activation is less frequent and requires binding of more Bcd molecules than P2 activation. Using a theoretical model to relate promoter activity to enhancer states, we show that the two promoters follow common transcription kinetics driven by sequential Bcd binding at the two enhancers. Bcd binding at either enhancer primarily activates P2, while P1 activation relies more on Bcd binding at both enhancers. These results provide a quantitative framework for understanding the kinetic mechanisms of complex eukaryotic gene regulation.

Alternative promoters differ in their expression patterns, whose mechanisms are not well understood. Here the authors show that alternative promoters of a Drosophila embryonic gene hunchback are regulated by different action modes of two enhancers.

Defining the evolutionary and gene regulatory logic of vertebrate neuronal diversity

How is the exquisite cellular complexity of the vertebrate brain achieved? In this issue of Neuron, Closser et al. (2022) reveal that an expanded neuronal gene regulatory landscape may drive evolutionary cellular diversification by providing complex context and cell-specific control of effector genes.

Concordant Androgen-Regulated Expression of Divergent Rhox5 Promoters in Sertoli Cells

Concordant transcriptional regulation can generate multiple gene products that collaborate to achieve a common goal. Here we report a case of concordant transcriptional regulation that instead drives a single protein to be produced in the same cell type from divergent promoters. This gene product-the RHOX5 homeobox transcription factor-is translated from 2 different mRNAs with different 5' untranslated regions (UTRs) transcribed from alternative promoters. Despite the fact that these 2 promoters-the proximal promoter (Pp) and the distal promoter (Pd)-exhibit different patterns of tissue-specific activity, share no obvious sequence identity, and depend on distinct transcription factors for expression, they exhibit a remarkably similar expression pattern in the testes. In particular, both depend on androgen signaling for expression in the testes, where they are specifically expressed in Sertoli cells and have a similar stage-specific expression pattern during the seminiferous epithelial cycle. We report evidence for 3 mechanisms that collaborate to drive concordant Pp/Pd expression. First, both promoters have an intrinsic ability to respond to androgen receptor and androgen. Second, the Pp acts as an enhancer to promote androgen-dependent transcription from the Pd. Third, Pd transcription is positively autoregulated by the RHOX5 protein, which is first produced developmentally from the Pp. Together, our data support a model in which the Rhox5 homeobox gene evolved multiple mechanisms to activate both of its promoters in Sertoli cells to produce Rhox5 in an androgen-dependent manner during different phases of spermatogenesis.

Jpx RNA regulates CTCF anchor site selection and formation of chromosome loops

Chromosome loops shift dynamically during development, homeostasis, and disease. CCCTC-binding factor (CTCF) is known to anchor loops and construct 3D genomes, but how anchor sites are selected is not yet understood. Here, we unveil Jpx RNA as a determinant of anchor selectivity. Jpx RNA targets thousands of genomic sites, preferentially binding promoters of active genes. Depleting Jpx RNA causes ectopic CTCF binding, massive shifts in chromosome looping, and downregulation of >700 Jpx target genes. Without Jpx, thousands of lost loops are replaced by de novo loops anchored by ectopic CTCF sites. Although Jpx controls CTCF binding on a genome-wide basis, it acts selectively at the subset of developmentally sensitive CTCF sites. Specifically, Jpx targets low-affinity CTCF motifs and displaces CTCF protein through competitive inhibition. We conclude that Jpx acts as a CTCF release factor and shapes the 3D genome by regulating anchor site usage.

The Ends Dictate the Means: Promoter Switching in Herpesvirus Gene Expression

Herpesvirus gene expression is dynamic and complex, with distinct complements of viral genes expressed at specific times in different infection contexts. These complex patterns of viral gene expression arise in part from the integration of multiple cellular and viral signals that affect the transcription of viral genes. The use of alternative promoters provides an increased level of control, allowing different promoters to direct the transcription of the same gene in response to distinct temporal and contextual cues. While once considered rare, herpesvirus alternative promoter usage was recently found to be far more pervasive and impactful than previously thought. Here we review several examples of promoter switching in herpesviruses and discuss the functional consequences on the transcriptional and post-transcriptional regulation of viral gene expression.

Transcriptional Regulation of the Human IL5RA Gene through Alternative Promoter Usage during Eosinophil Development

Regulation of the IL-5 receptor alpha (IL5RA) gene is complicated, with two known promoters (P1 and P2) driving transcription, and two known isoforms (transmembrane and soluble) dichotomously affecting the signaling potential of the protein products. Here, we sought to determine the patterns of P1 and P2 promoter usage and transcription factor occupancy during primary human eosinophil development from CD34⁺ hematopoietic stem cell progenitors. We found that during eosinophilopoiesis, both promoters were active but subject to distinct temporal regulation, coincident with combinatorial interactions of transcription factors, including GATA-1, PU.1, and C/EBP family members. P1 displayed a relatively constant level of activity throughout eosinophil development, while P2 activity peaked early and waned thereafter. The soluble IL-5Rα mRNA peaked early and showed the greatest magnitude fold-induction, while the signaling-competent transmembrane isoform peaked moderately. Two human eosinophilic cell lines whose relative use of P1 and P2 were similar to eosinophils differentiated in culture were used to functionally test putative transcription factor binding sites. Transcription factor occupancy was then validated in primary cultures by ChIP. We conclude that IL-5-dependent generation of eosinophils from CD34⁺ precursors involves complex and dynamic activity including both promoters, several interacting transcription factors, and both signaling and antagonistic protein products.

A pan-cancer analysis of alternative splicing of splicing factors in 6904 patients

Great progress has been made in the investigation on mutation and expression of splicing factor. However, little is known on the role of alternative splicing of splicing factors across cancers. Here, we reported a pan-cancer analysis of alternative splicing of splicing factors spanning 6904 patients across 16 cancer types, and identified 167 splicing factors with implications regulating cancer-specific splicing patterns through alternative splicing. Furthermore, we found that abnormal splicing events of splicing factors could serve as potential common regulators for alternative splicing in different cancers. In addition, we developed a splicing-derived neoepitopes database (ASPNs), which provided the corresponding putative alternative splicing-derived neoepitopes of 16 cancer types. Our results suggested that alternative splicing of splicing factors involved in the pre-RNA splicing process was common across cancer types and may represent an underestimated hallmark of tumorigenesis.

DNABERT: pre-trained Bidirectional Encoder Representations from Transformers model for DNA-language in genome

MOTIVATION

Deciphering the language of non-coding DNA is one of the fundamental problems in genome research. Gene regulatory code is highly complex due to the existence of polysemy and distant semantic relationship, which previous informatics methods often fail to capture especially in data-scarce scenarios.

RESULTS

To address this challenge, we developed a novel pre-trained bidirectional encoder representation, named DNABERT, to capture global and transferrable understanding of genomic DNA sequences based on up and downstream nucleotide contexts. We compared DNABERT to the most widely used programs for genome-wide regulatory elements prediction and demonstrate its ease of use, accuracy and efficiency. We show that the single pre-trained transformers model can simultaneously achieve state-of-the-art performance on prediction of promoters, splice sites and transcription factor binding sites, after easy fine-tuning using small task-specific labeled data. Further, DNABERT enables direct visualization of nucleotide-level importance and semantic relationship within input sequences for better interpretability and accurate identification of conserved sequence motifs and functional genetic variant candidates. Finally, we demonstrate that pre-trained DNABERT with human genome can even be readily applied to other organisms with exceptional performance. We anticipate that the pre-trained DNABERT model can be fined tuned to many other sequence analyses tasks.

AVAILABILITY AND IMPLEMENTATION

The source code, pretrained and finetuned model for DNABERT are available at GitHub (https://github.com/jerryji1993/DNABERT).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evolution of tissue and developmental specificity of transcription start sites in Bos taurus indicus

To further the understanding of the evolution of transcriptional regulation, we profiled genome-wide transcriptional start sites (TSSs) in two sub-species, Bos taurus taurus and Bos taurus indicus, that diverged approximately 500,000 years ago. Evolutionary and developmental-stage differences in TSSs were detected across the sub-species, including translocation of dominant TSS and changes in TSS distribution. The 16% of all SNPs located in significant differentially used TSS clusters across sub-species had significant shifts in allele frequency (472 SNPs), indicating they may have been subject to selection. In spleen and muscle, a higher relative TSS expression was observed in Bos indicus than Bos taurus for all heat shock protein genes, which may be responsible for the tropical adaptation of Bos indicus.

Forutan et al. measure RNA expression and map transcription start sites (TSSs) between two sub-species of cattle in adult and fetal tissues and demonstrate translocation and changes in TSS distribution between the sub-species. This study provides insight into cattle adaptation and provides support of rapidly evolving TSS between two recently diverged subspecies.

Alternative Promoter Use Governs the Expression of IgLON Cell Adhesion Molecules in Histogenetic Fields of the Embryonic Mouse Brain

The members of the IgLON superfamily of cell adhesion molecules facilitate fundamental cellular communication during brain development, maintain functional brain circuitry, and are associated with several neuropsychiatric disorders such as depression, autism, schizophrenia, and intellectual disabilities. Usage of alternative promoter-specific 1a and 1b mRNA isoforms in Lsamp, Opcml, Ntm, and the single promoter of Negr1 in the mouse and human brain has been previously described. To determine the precise spatiotemporal expression dynamics of Lsamp, Opcml, Ntm isoforms, and Negr1, in the developing brain, we generated isoform-specific RNA probes and carried out in situ hybridization in the developing (embryonic, E10.5, E11.5, 13.5, 17; postnatal, P0) and adult mouse brains. We show that promoter-specific expression of IgLONs is established early during pallial development (at E10.5), where it remains throughout its differentiation through adulthood. In the diencephalon, midbrain, and hindbrain, strong expression patterns are initiated a few days later and begin fading after birth, being only faintly expressed during adulthood. Thus, the expression of specific IgLONs in the developing brain may provide the means for regionally specific functionality as well as for specific regional vulnerabilities. The current study will therefore improve the understanding of how IgLON genes are implicated in the development of neuropsychiatric disorders.

ExTraMapper: Exon- and Transcript-level mappings for orthologous gene pairs

MOTIVATION

Access to large-scale genomics and transcriptomics data from various tissues and cell lines allowed the discovery of wide-spread alternative splicing events and alternative promoter usage in mammalians. Between human and mouse, gene-level orthology is currently present for nearly 16k protein-coding genes spanning a diverse repertoire of over 200k total transcript isoforms.

RESULTS

Here, we describe a novel method, ExTraMapper, which leverages sequence conservation between exons of a pair of organisms and identifies a fine-scale orthology mapping at the exon and then transcript level. ExTraMapper identifies more than 350k exon mappings, as well as 30k transcript mappings between human and mouse using only sequence and gene annotation information. We demonstrate that ExTraMapper identifies a larger number of exon and transcript mappings compared to previous methods. Further, it identifies exon fusions, splits, and losses due to splice site mutations, and finds mappings between microexons that are previously missed. By reanalysis of RNA-seq data from 13 matched human and mouse tissues, we show that ExTraMapper improves the correlation of transcript-specific expression levels suggesting a more accurate mapping of human and mouse transcripts. We also applied the method to detect conserved exon and transcript pairs between human and rhesus macaque genomes to highlight the point that ExTraMapper is applicable to any pair of organisms that have orthologous gene pairs.

AVAILABILITY

The source code and the results are available at https://github.com/ay-lab/ExTraMapper and http://ay-lab-tools.lji.org/extramapper.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.