Complex controls: the role of alternative promoters in mammalian genomes

Hidden secrets of the cancer genome: unlocking the impact of non-coding mutations in gene regulatory elements

Discoveries in the field of genomics have revealed that non-coding genomic regions are not merely "junk DNA", but rather comprise critical elements involved in gene expression. These gene regulatory elements (GREs) include enhancers, insulators, silencers, and gene promoters. Notably, new evidence shows how mutations within these regions substantially influence gene expression programs, especially in the context of cancer. Advances in high-throughput sequencing technologies have accelerated the identification of somatic and germline single nucleotide mutations in non-coding genomic regions. This review provides an overview of somatic and germline non-coding single nucleotide alterations affecting transcription factor binding sites in GREs, specifically involved in cancer biology. It also summarizes the technologies available for exploring GREs and the challenges associated with studying and characterizing non-coding single nucleotide mutations. Understanding the role of GRE alterations in cancer is essential for improving diagnostic and prognostic capabilities in the precision medicine era, leading to enhanced patient-centered clinical outcomes.

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.

Iterative oxidation by TET1 is required for reprogramming of imprinting control regions and patterning of mouse sperm hypomethylated regions

Ten-eleven translocation (TET) enzymes iteratively oxidize 5-methylcytosine (5mC) to generate 5-hydroxymethylcytosine (5hmC), 5-formylcytosine, and 5-carboxylcytosine to facilitate active genome demethylation. Whether these bases are required to promote replication-coupled dilution or activate base excision repair during mammalian germline reprogramming remains unresolved due to the inability to decouple TET activities. Here, we generated two mouse lines expressing catalytically inactive TET1 (Tet1-HxD) and TET1 that stalls oxidation at 5hmC (Tet1-V). Tet1 knockout and catalytic mutant primordial germ cells (PGCs) fail to erase methylation at select imprinting control regions and promoters of meiosis-associated genes, validating the requirement for the iterative oxidation of 5mC for complete germline reprogramming. TET1V and TET1HxD rescue most hypermethylation of Tet1-/- sperm, suggesting the role of TET1 beyond its oxidative capability. We additionally identify a broader class of hypermethylated regions in Tet1 mutant mouse sperm that depend on TET oxidation for reprogramming. Our study demonstrates the link between TET1-mediated germline reprogramming and sperm methylome patterning.

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.

Transcriptional determinants of lipid mobilization in human adipocytes

Defects in adipocyte lipolysis drive multiple aspects of cardiometabolic disease, but the transcriptional framework controlling this process has not been established. To address this, we performed a targeted perturbation screen in primary human adipocytes. Our analyses identified 37 transcriptional regulators of lipid mobilization, which we classified as (i) transcription factors, (ii) histone chaperones, and (iii) mRNA processing proteins. On the basis of its strong relationship with multiple readouts of lipolysis in patient samples, we performed mechanistic studies on one hit, ZNF189, which encodes the zinc finger protein 189. Using mass spectrometry and chromatin profiling techniques, we show that ZNF189 interacts with the tripartite motif family member TRIM28 and represses the transcription of an adipocyte-specific isoform of phosphodiesterase 1B (PDE1B2). The regulation of lipid mobilization by ZNF189 requires PDE1B2, and the overexpression of PDE1B2 is sufficient to attenuate hormone-stimulated lipolysis. Thus, our work identifies the ZNF189-PDE1B2 axis as a determinant of human adipocyte lipolysis and highlights a link between chromatin architecture and lipid mobilization.

Characterization of NFE2L1-616, an isoform of nuclear factor-erythroid-2 related transcription factor-1 that activates antioxidant response element-regulated genes

The NFE2L1 transcription factor (aka Nrf1) is a basic leucine zipper protein that performs a critical role in the cellular stress response pathway. Here, we characterized a novel variant of NFE2L1 referred to as NFE2L1-616. The transcript encoding NFE2L1-616 is derived from an intronic promoter, and it has a distinct first exon than other reported full-length NFE2L1 isoforms. The NFE2L1-616 protein constitutively localizes in the nucleus as it lacks the N-terminal amino acid residues that targets other full-length NFE2L1 isoforms to the endoplasmic reticulum. The expression level of NFE2L1-616 is lower than other NFE2L1 isoforms. It is widely expressed across different cell lines and tissues that were examined. NFE2L1-616 showed strong transcriptional activity driving luciferase reporter expression from a promoter containing antioxidant response element. Together, the results suggest that NFE2L1-616 variant can function as a positive regulator in the transcriptional regulation of NFE2L1 responsive genes.

Alternative splicing in Atlantic salmon head kidney and SHK-1 cell line during the Piscirickettsia salmonis infection: A comparative transcriptome survey

Piscirickettsia salmonis, an intracellular bacterium in salmon aquaculture, is a big challenge because it is responsible for 54.2% of Atlantic salmon mortalities. In recent years, the high relevance of Alternative Splicing (AS) as a molecular mechanism associated with infectious conditions and host-pathogen interaction processes, especially in host immune activation, has been observed. Several studies have highlighted the role of AS in the host's immune response during viral, bacterial, and endoparasite infection. In the present study, we evaluated AS transcriptome profiles during P. salmonis infection in the two most used study models, SHK-1 cell line and salmon head kidney tissue. First, the SHK-1 cell line was exposed to P. salmonis infection at 0-, 7-, and 14-days post-infection (dpi). Following, total RNA was extracted for Illumina sequencing. On the other hand, RNA-Seq datasets of Atlantic salmon head kidney infected with the same P. salmonis strayingwase used. For both study models, the highest number of differentially alternative splicing (DAS) events was observed at 7 dpi, 16,830 DAS events derived from 9213 DAS genes in SHK-1 cells, and 13,820 DAS events from 7684 DAS genes in salmon HK. Alternative first exon (AF) was the most abundant AS type in the three infection times analyzed, representing 31% in SHK-1 cells and 228.6 in salmon HK; meanwhile, mutually exclusive exon (MX) was the least abundant. Notably, functional annotation of DAS genes in SHK-1 cells infected with P. salmonis showed a high presence of genes related to nucleotide metabolism. In contrast, the salmon head kidney exhibited many GO terms associated with immune response. Our findings reported the role of AS during P. salmonis infection in Atlantic salmon. These studies would contribute to a better understanding of the molecular bases that support the pathogen-host interaction, evidencing the contribution of AS regulating the transcriptional host response.

Identification of Splice Variants and Isoforms in Transcriptomics and Proteomics

Alternative splicing is pivotal to the regulation of gene expression and protein diversity in eukaryotic cells. The detection of alternative splicing events requires specific omics technologies. Although short-read RNA sequencing has successfully supported a plethora of investigations on alternative splicing, the emerging technologies of long-read RNA sequencing and top-down mass spectrometry open new opportunities to identify alternative splicing and protein isoforms with less ambiguity. Here, we summarize improvements in short-read RNA sequencing for alternative splicing analysis, including percent splicing index estimation and differential analysis. We also review the computational methods used in top-down proteomics analysis regarding proteoform identification, including the construction of databases of protein isoforms and statistical analyses of search results. While many improvements in sequencing and computational methods will result from emerging technologies, there should be future endeavors to increase the effectiveness, integration, and proteome coverage of alternative splicing events.

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.

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.

TET1 Catalytic Activity is Required for Reprogramming of Imprinting Control Regions and Patterning of Sperm-Specific Hypomethylated Regions

DNA methylation erasure is required for mammalian primordial germ cell reprogramming. TET enzymes iteratively oxidize 5-methylcytosine to generate 5-hyroxymethylcytosine (5hmC), 5-formylcytosine, and 5-carboxycytosine to facilitate active genome demethylation. Whether these bases are required to promote replication-coupled dilution or activate base excision repair during germline reprogramming remains unresolved due to the lack of genetic models that decouple TET activities. Here, we generated two mouse lines expressing catalytically inactive TET1 ( Tet1-HxD ) and TET1 that stalls oxidation at 5hmC ( Tet1-V ). Tet1 -/- , Tet1 V/V , and Tet1 HxD/HxD sperm methylomes show that TET1 V and TET1 HxD rescue most Tet1 -/- hypermethylated regions, demonstrating the importance of TET1’s extra-catalytic functions. Imprinted regions, in contrast, require iterative oxidation. We further reveal a broader class of hypermethylated regions in sperm of Tet1 mutant mice that are excluded from de novo methylation during male germline development and depend on TET oxidation for reprogramming. Our study underscores the link between TET1-mediated demethylation during reprogramming and sperm methylome patterning.

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.

Integrative analysis of Iso-Seq and RNA-seq reveals dynamic changes of alternative promoter, alternative splicing and alternative polyadenylation during Angiotensin II-induced senescence in rat primary aortic endothelial cells

In eukaryotes, alternative promoter (AP), alternative splicing (AS), and alternative polyadenylation (APA) are three crucial regulatory mechanisms that modulate message RNA (mRNA) diversity. Although AP, AS and APA are involved in diverse biological processess, whether they have dynamic changes in Angiotensin II (Ang II) induced senescence in rat primary aortic endothelial cells (RAECs), an important cellular model for studying cardiovascular disease, remains unclear. Here we integrated both PacBio single-molecule long-read isoform sequencing (Iso-Seq) and Illumina short-read RNA sequencing (RNA-seq) to analyze the changes of AP, AS and APA in Ang II-induced senescent RAECs. Iso-Seq generated 36,278 isoforms from 10,145 gene loci and 65.81% of these isoforms are novel, which were further cross-validated by public data obtained by other techonologies such as CAGE, PolyA-Seq and 3'READS. APA contributed most to novel isoforms, followed by AS and AP. Further investigation showed that AP, AS and APA could all contribute to the regulation of isoform, but AS has more dynamic changes compared to AP and APA upon Ang II stimulation. Genes undergoing AP, AS and APA in Ang II-treated cells are enriched in various pathways related to aging or senescence, suggesting that these molecular changes are involved in functional alterations during Ang II-induced senescence. Together, the present study largely improved the annotation of rat genome and revealed gene expression changes at isoform level, extending the understanding of the complexity of gene regulation in Ang II-treated RAECs, and also provided novel clues for discovering the regulatory mechanism undelying Ang II caused vascular senescence and diseases.

LAFITE Reveals the Complexity of Transcript Isoforms in Subcellular Fractions

Characterization of the subcellular distribution of RNA is essential for understanding the molecular basis of biological processes. Here, the subcellular nanopore direct RNA-sequencing (DRS) of four lung cancer cell lines (A549, H1975, H358, and HCC4006) is performed, coupled with a computational pipeline, Low-abundance Aware Full-length Isoform clusTEr (LAFITE), to comprehensively analyze the full-length cytoplasmic and nuclear transcriptome. Using additional DRS and orthogonal data sets, it is shown that LAFITE outperforms current methods for detecting full-length transcripts, particularly for low-abundance isoforms that are usually overlooked due to poor read coverage. Experimental validation of six novel isoforms exclusively identified by LAFITE further confirms the reliability of this pipeline. By applying LAFITE to subcellular DRS data, the complexity of the nuclear transcriptome is revealed in terms of isoform diversity, 3'-UTR usage, m6A modification patterns, and intron retention. Overall, LAFITE provides enhanced full-length isoform identification and enables a high-resolution view of the RNA landscape at the isoform level.

Hippocampal ceRNA networks from chronic intermittent ethanol vapor-exposed male mice and functional analysis of top-ranked lncRNA genes for ethanol drinking phenotypes

The molecular mechanisms regulating the development and progression of alcohol use disorder (AUD) are largely unknown. While noncoding RNAs have previously been implicated as playing key roles in AUD, long-noncoding RNA (lncRNA) remains understudied in relation to AUD. In this study, we first identified ethanol-responsive lncRNAs in the mouse hippocampus that are transcriptional network hub genes. Microarray analysis of lncRNA, miRNA, circular RNA, and protein coding gene expression in the hippocampus from chronic intermittent ethanol vapor- or air- (control) exposed mice was used to identify ethanol-responsive competing endogenous RNA (ceRNA) networks. Highly interconnected lncRNAs (genes that had the strongest overall correlation to all other dysregulated genes identified) were ranked. The top four lncRNAs were novel, previously uncharacterized genes named Gm42575, 4930413E15Rik, Gm15767, and Gm33447, hereafter referred to as Pitt1, Pitt2, Pitt3, and Pitt4, respectively. We subsequently tested the hypothesis that CRISPR/Cas9 mutagenesis of the putative promoter and first exon of these lncRNAs in C57BL/6J mice would alter ethanol drinking behavior. The Drinking in the Dark (DID) assay was used to examine binge-like drinking behavior, and the Every-Other-Day Two-Bottle Choice (EOD-2BC) assay was used to examine intermittent ethanol consumption and preference. No significant differences between control and mutant mice were observed in the DID assay. Female-specific reductions in ethanol consumption were observed in the EOD-2BC assay for Pitt1, Pitt3, and Pitt4 mutant mice compared to controls. Male-specific alterations in ethanol preference were observed for Pitt1 and Pitt2. Female-specific increases in ethanol preference were observed for Pitt3 and Pitt4. Total fluid consumption was reduced in Pitt1 and Pitt2 mutants at 15% v/v ethanol and in Pitt3 and Pitt4 at 20% v/v ethanol in females only. We conclude that all lncRNAs targeted altered ethanol drinking behavior, and that lncRNAs Pitt1, Pitt3, and Pitt4 influenced ethanol consumption in a sex-specific manner. Further research is necessary to elucidate the biological mechanisms for these effects. These findings add to the literature implicating noncoding RNAs in AUD and suggest lncRNAs also play an important regulatory role in the disease.

Z-DNA-Containing Long Terminal Repeats of Human Endogenous Retrovirus Families Provide Alternative Promoters for Human Functional Genes

Transposable elements (TEs) account for approximately 45% of the human genome. TEs have proliferated randomly and integrated into functional genes during hominoid radiation. They appear as right-handed B-DNA double helices and slightly elongated left-handed Z-DNAs. Human endogenous retrovirus (HERV) families are widely distributed in human chromosomes at a ratio of 8%. They contain a 5'-long terminal repeat (LTR)-gag-pol-env-3'-LTR structure. LTRs contain the U3 enhancer and promoter region, transcribed R region, and U5 region. LTRs can influence host gene expression by acting as regulatory elements. In this review, we describe the alternative promoters derived from LTR elements that overlap Z-DNA by comparing Z-hunt and DeepZ data for human functional genes. We also present evidence showing the regulatory activity of LTR elements containing Z-DNA in GSDML. Taken together, the regulatory activity of LTR elements with Z-DNA allows us to understand gene function in relation to various human diseases.

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.

SMAD4-201 transcript as a putative biomarker in colorectal cancer

Background

Transcripts with alternative 5′-untranslated regions (UTRs) result from the activity of alternative promoters and they can determine gene expression by influencing its stability and translational efficiency, thus executing complex regulation of developmental, physiological and pathological processes. Transcriptional regulation of human SMAD4, a key tumor suppressor deregulated in most gastrointestinal cancers, entails four alternative promoters. These promoters and alternative transcripts they generate remain unexplored as contributors to the SMAD4 deregulation in cancer. The aim of this study was to investigate the relative abundance of the transcript SMAD4–201 in colorectal cell lines and tissues in order to establish if its fluctuations may be associated with colorectal cancer (CRC).

Methods

Relative abundance of SMAD4–201 in total SMAD4 mRNA was analyzed using quantitative PCR in a set of permanent human colon cell lines and tumor and corresponding healthy tissue samples from patients with CRC.

Results

The relative abundance of SMAD4–201 in analyzed cell lines varied between 16 and 47%. A similar relative abundance of SMAD4–201 transcript was found in the majority of analyzed human tumor tissue samples, and it was averagely 20% lower in non-malignant in comparison to malignant tissue samples (p = 0.001). Transcript SMAD4–202 was not detectable in any of the analyzed samples, so the observed fluctuations in the composition of SMAD4 transcripts can be attributed to transcripts other than SMAD4–201 and SMAD4–202.

Conclusion

The expression profile of SMAD4–201 in human tumor and non-tumor tissue samples may indicate the translational potential of this molecule in CRC, but further research is needed to clarify its usability as a potential biomarker for early diagnosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09186-z.

Functional Annotation of Custom Transcriptomes

Many eukaryotic genes can give rise to different alternative transcripts depending on stage of development, cell type, and physiological cues. Current transcriptome-wide sequencing technologies highlight the remarkable extent of this regulation in metazoans and allow for RNA isoforms to be profiled in increasingly small biological samples and with a growing confidence. Understanding biological functions of sample-specific transcripts is a major challenge in genomics and RNA processing fields. Here we describe simple bioinformatics workflows that facilitate this task by streamlining reference-guided annotation of novel transcripts. A key part of our protocol is the R package factR that rapidly matches custom-assembled transcripts to their likely host genes, deduces the sequence and domain structure of novel protein products, and predicts sensitivity of newly identified RNA isoforms to nonsense-mediated decay.

Two promoters integrate multiple enhancer inputs to drive wild-type knirps expression in the Drosophila melanogaster embryo

Proper development depends on precise spatiotemporal gene expression patterns. Most developmental genes are regulated by multiple enhancers and often by multiple core promoters that generate similar transcripts. We hypothesize that multiple promoters may be required either because enhancers prefer a specific promoter or because multiple promoters serve as a redundancy mechanism. To test these hypotheses, we studied the expression of the knirps locus in the early Drosophila melanogaster embryo, which is mediated by multiple enhancers and core promoters. We found that one of these promoters resembles a typical "sharp" developmental promoter, while the other resembles a "broad" promoter usually associated with housekeeping genes. Using synthetic reporter constructs, we found that some, but not all, enhancers in the locus show a preference for one promoter, indicating that promoters provide both redundancy and specificity. By analyzing the reporter dynamics, we identified specific burst properties during the transcription process, namely burst size and frequency, that are most strongly tuned by the combination of promoter and enhancer. Using locus-sized reporters, we discovered that enhancers with no promoter preference in a synthetic setting have a preference in the locus context. Our results suggest that the presence of multiple promoters in a locus is due both to enhancer preference and a need for redundancy and that "broad" promoters with dispersed transcription start sites are common among developmental genes. They also imply that it can be difficult to extrapolate expression measurements from synthetic reporters to the locus context, where other variables shape a gene's overall expression pattern.

Canonical and Divergent N-Terminal HBx Isoform Proteins Unveiled: Characteristics and Roles during HBV Replication

Hepatitis B virus (HBV) X protein (HBx) is a viral regulatory and multifunctional protein. It is well-known that the canonical HBx reading frame bears two phylogenetically conserved internal in-frame translational initiation codons at Met2 and Met3, thus possibly generating divergent N-terminal smaller isoforms during translation. Here, we demonstrate that the three distinct HBx isoforms are generated from the ectopically expressed HBV HBx gene, named XF (full-length), XM (medium-length), and XS (short-length); they display different subcellular localizations when expressed individually in cultured hepatoma cells. Particularly, the smallest HBx isoform, XS, displayed a predominantly cytoplasmic localization. To study HBx proteins during viral replication, we performed site-directed mutagenesis to target the individual or combinatorial expression of the HBx isoforms within the HBV viral backbone (full viral genome). Our results indicate that of all HBx isoforms, only the smallest HBx isoform, XS, can restore WT levels of HBV replication, and bind to the viral mini chromosome, thereby establishing an active chromatin state, highlighting its crucial activities during HBV replication. Intriguingly, we found that sequences of HBV HBx genotype H are devoid of the conserved Met3 position, and therefore HBV genotype H infection is naturally silent for the expression of the HBx XS isoform. Finally, we found that the HBx XM (medium-length) isoform shares significant sequence similarity with the N-terminus domain of the COMMD8 protein, a member of the copper metabolism MURR1 domain-containing (COMMD) protein family. This novel finding might facilitate studies on the phylogenetic origin of the HBV X protein. The identification and functional characterization of its isoforms will shift the paradigm by changing the concept of HBx from being a unique, canonical, and multifunctional protein toward the occurrence of different HBx isoforms, carrying out different overlapping functions at different subcellular localizations during HBV genome replication. Significantly, our current work unveils new crucial HBV targets to study for potential antiviral research, and human virus pathogenesis.

Heterozygous Tropomodulin 3 mice have improved lung vascularization after chronic hypoxia

The molecular mechanisms leading to high altitude pulmonary hypertension (HAPH) remains poorly understood. We previously analyzed the whole genome sequence of Kyrgyz highland population and identified eight genomic intervals having a potential role in HAPH. Tropomodulin 3 gene (TMOD3) which encodes a protein that binds and caps the pointed ends of actin filaments and inhibits cell migration, was one of the top candidates. Here we systematically sought additional evidence to validate the functional role of TMOD3. In-silico analysis reveals that some of the SNPs in HAPH associated genomic intervals were positioned in a regulatory region that could result in alternative splicing of TMOD3. In order to functionally validate the role of TMOD3 in HAPH, we exposed Tmod3-/+ mice to 4 weeks of constant hypoxia, i.e. 10% O2 and analyzed both functional (hemodynamic measurements) and structural (angiography) parameters related to HAPH. The hemodynamic measurements, such as right ventricular systolic pressure, a surrogate measure for pulmonary arterial systolic pressure, and right ventricular contractility (RV- ± dP/dt), increases with hypoxia did not separate between Tmod3-/+ and control mice. Remarkably, there was a significant increase in the number of lung vascular branches and total length of pulmonary vascular branches (p < 0.001) in Tmod3-/+ after 4 weeks of constant hypoxia as compared to controls. Notably, the Tmod3-/+ endothelial cells migration was also significantly higher than that from the wild-type littermates. Our results indicate that, under chronic hypoxia, lower levels of Tmod3 play an important role in the maintenance or neo-vascularization of pulmonary arteries.

The ASIP gene in the llama (Lama glama): Alternative transcripts, expression and relation with color phenotypes

The Agouti gene (ASIP) is one of the most important genes for coat color determination in mammals. It has a complex structure with several promoters and alternative non-coding first exons that are transcribed into mRNAs with different 5'UTR. These mRNA isoforms regulate the temporal and spatial expression of the gene, producing diverse pigmentation patterns. Here, we studied ASIP transcriptional variants and their expression in the skin of llamas with different coat color phenotypes. We also described the ASIP locus, including promoter usage and the splicing events that originate each transcript variant. Using 5'RACE-PCR we isolated seven ASIP transcripts with alternative 5'UTR, where exons 1A, 1A', 1C, 1D, and a novel non-coding exon 1A" were identified. Additionally, new alternative spliced forms were found. The diversity of ASIP 5'UTRs is originated by a complex pattern of alternative promoter usage, multiple transcription start sites and splicing events that include exon skipping and alternative 3' splicing site selection. We found that ASIP was highly expressed in llamas with white and brown phenotypes while black animals presented very low expression. The main responsible for this difference was a fusion transcript between ASIP and NCOA6 genes, which was present in the skin of white and brown llamas but not in the black ones. The rest of ASIP transcripts presented very low expression in the skin, indicating that the main regulation point for ASIP gene expression is at the transcriptional level. Nevertheless, the characteristics of the 5'UTRs sequences suggest that alternative transcripts could be regulated differently at the protein synthesis level.

A CpG island promoter drives the CXXC5 gene expression

CXXC5 is a member of the zinc-finger CXXC family that binds to unmethylated CpG dinucleotides. CXXC5 modulates gene expressions resulting in diverse cellular events mediated by distinct signaling pathways. However, the mechanism responsible for CXXC5 expression remains largely unknown. We found here that of the 14 annotated CXXC5 transcripts with distinct 5' untranslated regions encoding the same protein, transcript variant 2 with the highest expression level among variants represents the main transcript in cell models. The DNA segment in and at the immediate 5'-sequences of the first exon of variant 2 contains a core promoter within which multiple transcription start sites are present. Residing in a region with high G-C nucleotide content and CpG repeats, the core promoter is unmethylated, deficient in nucleosomes, and associated with active RNA polymerase-II. These findings suggest that a CpG island promoter drives CXXC5 expression. Promoter pull-down revealed the association of various transcription factors (TFs) and transcription co-regulatory proteins, as well as proteins involved in histone/chromatin, DNA, and RNA processing with the core promoter. Of the TFs, we verified that ELF1 and MAZ contribute to CXXC5 expression. Moreover, the first exon of variant 2 may contain a G-quadruplex forming region that could modulate CXXC5 expression.

Inflammation drives alternative first exon usage to regulate immune genes including a novel iron-regulated isoform of Aim2

Determining the layers of gene regulation within the innate immune response is critical to our understanding of the cellular responses to infection and dysregulation in disease. We identified a conserved mechanism of gene regulation in human and mouse via changes in alternative first exon (AFE) usage following inflammation, resulting in changes to the isoforms produced. Of these AFE events, we identified 95 unannotated transcription start sites in mice using a de novo transcriptome generated by long-read native RNA-sequencing, one of which is in the cytosolic receptor for dsDNA and known inflammatory inducible gene, Aim2. We show that this unannotated AFE isoform of Aim2 is the predominant isoform expressed during inflammation and contains an iron-responsive element in its 5′UTR enabling mRNA translation to be regulated by iron levels. This work highlights the importance of examining alternative isoform changes and translational regulation in the innate immune response and uncovers novel regulatory mechanisms of Aim2.

Structure and function of DEAH-box helicase 32 and its role in cancer

DEAH-box helicase 32 (DHX32) is an RNA helicase with unique structural characteristics that is involved in numerous biological processes associated with RNA, including ribosome biosynthesis, transcription, mRNA splicing and translation. Increasing evidence suggests that abnormal DHX32 expression contributes to cancer initiation and development, due to dysregulated cell proliferation, differentiation, apoptosis and other processes. In the current review, the discovery, structure and function of DHX32, as well as the association between abnormal DHX32 expression and tumors are discussed. DHX32 expression is downregulated in acute lymphoblastic leukemia, but upregulated in solid tumors, including colorectal and breast cancer. Furthermore, DHX32 expression levels are associated with the pathological and clinical features of the cancer. Therefore, DHX32 may serve as a novel liquid biopsy marker for auxiliary diagnosis and prognosis screening, as well as a possible target for cancer therapy. The molecular mechanism underlying the contribution of DHX32 towards the initiation and development of cancer requires further investigation for the development of anticancer treatments based on manipulating DHX32 expression and function.

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.

Deletion of a non-canonical regulatory sequence causes loss of Scn1a expression and epileptic phenotypes in mice

BACKGROUND

Genes with multiple co-active promoters appear common in brain, yet little is known about functional requirements for these potentially redundant genomic regulatory elements. SCN1A, which encodes the NaV1.1 sodium channel alpha subunit, is one such gene with two co-active promoters. Mutations in SCN1A are associated with epilepsy, including Dravet syndrome (DS). The majority of DS patients harbor coding mutations causing SCN1A haploinsufficiency; however, putative causal non-coding promoter mutations have been identified.

METHODS

To determine the functional role of one of these potentially redundant Scn1a promoters, we focused on the non-coding Scn1a 1b regulatory region, previously described as a non-canonical alternative transcriptional start site. We generated a transgenic mouse line with deletion of the extended evolutionarily conserved 1b non-coding interval and characterized changes in gene and protein expression, and assessed seizure activity and alterations in behavior.

RESULTS

Mice harboring a deletion of the 1b non-coding interval exhibited surprisingly severe reductions of Scn1a and NaV1.1 expression throughout the brain. This was accompanied by electroencephalographic and thermal-evoked seizures, and behavioral deficits.

CONCLUSIONS

This work contributes to functional dissection of the regulatory wiring of a major epilepsy risk gene, SCN1A. We identified the 1b region as a critical disease-relevant regulatory element and provide evidence that non-canonical and seemingly redundant promoters can have essential function.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

Overview of PAX gene family: analysis of human tissue-specific variant expression and involvement in human disease

Paired-box (PAX) genes encode a family of highly conserved transcription factors found in vertebrates and invertebrates. PAX proteins are defined by the presence of a paired domain that is evolutionarily conserved across phylogenies. Inclusion of a homeodomain and/or an octapeptide linker subdivides PAX proteins into four groups. Often termed "master regulators", PAX proteins orchestrate tissue and organ development throughout cell differentiation and lineage determination, and are essential for tissue structure and function through maintenance of cell identity. Mutations in PAX genes are associated with myriad human diseases (e.g., microphthalmia, anophthalmia, coloboma, hypothyroidism, acute lymphoblastic leukemia). Transcriptional regulation by PAX proteins is, in part, modulated by expression of alternatively spliced transcripts. Herein, we provide a genomics update on the nine human PAX family members and PAX homologs in 16 additional species. We also present a comprehensive summary of human tissue-specific PAX transcript variant expression and describe potential functional significance of PAX isoforms. While the functional roles of PAX proteins in developmental diseases and cancer are well characterized, much remains to be understood regarding the functional roles of PAX isoforms in human health. We anticipate the analysis of tissue-specific PAX transcript variant expression presented herein can serve as a starting point for such research endeavors.

Long-read transcriptome sequencing reveals abundant promoter diversity in distinct molecular subtypes of gastric cancer

BACKGROUND

Deregulated gene expression is a hallmark of cancer; however, most studies to date have analyzed short-read RNA sequencing data with inherent limitations. Here, we combine PacBio long-read isoform sequencing (Iso-Seq) and Illumina paired-end short-read RNA sequencing to comprehensively survey the transcriptome of gastric cancer (GC), a leading cause of global cancer mortality.

RESULTS

We performed full-length transcriptome analysis across 10 GC cell lines covering four major GC molecular subtypes (chromosomal unstable, Epstein-Barr positive, genome stable and microsatellite unstable). We identify 60,239 non-redundant full-length transcripts, of which > 66% are novel compared to current transcriptome databases. Novel isoforms are more likely to be cell line and subtype specific, expressed at lower levels with larger number of exons, with longer isoform/coding sequence lengths. Most novel isoforms utilize an alternate first exon, and compared to other alternative splicing categories, are expressed at higher levels and exhibit higher variability. Collectively, we observe alternate promoter usage in 25% of detected genes, with the majority (84.2%) of known/novel promoter pairs exhibiting potential changes in their coding sequences. Mapping these alternate promoters to TCGA GC samples, we identify several cancer-associated isoforms, including novel variants of oncogenes. Tumor-specific transcript isoforms tend to alter protein coding sequences to a larger extent than other isoforms. Analysis of outcome data suggests that novel isoforms may impart additional prognostic information.

CONCLUSIONS

Our results provide a rich resource of full-length transcriptome data for deeper studies of GC and other gastrointestinal malignancies.

Getting started: altering promoter choice as a mechanism for cell type differentiation

In this issue of Genes & Development, Lu and colleagues (pp. 663-677) have discovered a key new mechanism of alternative promoter choice that is involved in differentiation of spermatocytes. Promoter choice has strong potential as mechanism for differentiation of many different cell types.

RNA splicing alteration in the response to platinum chemotherapy in ovarian cancer: A possible biomarker and therapeutic target

Since its discovery, alternative splicing has been recognized as a powerful way for a cell to amplify the genetic information and for a living organism to adapt, evolve, and survive. We now know that a very high number of genes are regulated by alternative splicing and that alterations of splicing have been observed in different types of human diseases, including cancer. Here, we review the accumulating knowledge that links the regulation of alternative splicing to the response to chemotherapy, focusing our attention on ovarian cancer and platinum-based treatments. Moreover, we discuss how expanding information could be exploited to identify new possible biomarkers of platinum response, to better select patients, and/or to design new therapies able to overcome platinum resistance.

Modulating gene regulation to treat genetic disorders

Over a thousand diseases are caused by mutations that alter gene expression levels. The potential of nuclease-deficient zinc fingers, TALEs or CRISPR fusion systems to treat these diseases by modulating gene expression has recently emerged. These systems can be applied to modify the activity of gene-regulatory elements - promoters, enhancers, silencers and insulators, subsequently changing their target gene expression levels to achieve therapeutic benefits - an approach termed cis-regulation therapy (CRT). Here, we review emerging CRT technologies and assess their therapeutic potential for treating a wide range of diseases caused by abnormal gene dosage. The challenges facing the translation of CRT into the clinic are discussed.

Prognostic and predictive value of PD-L2 DNA methylation and mRNA expression in melanoma

BACKGROUND

PD-L1 (programmed cell death 1 ligand 1) expression in melanoma has been associated with a better response to anti-PD-1 (programmed cell death 1) therapy. However, patients with PD-L1-negative melanomas can respond to anti-PD-1 blockade, suggesting that the other PD-1 ligand, PD-L2 (programmed cell death 1 ligand 2), might also be relevant for efficacy of PD-1 inhibition. We investigated PD-L2 expression and methylation as a prognostic and predictive biomarker in melanoma.

METHODS

DNA methylation at five CpG loci and gene expression of PD-L2 were evaluated with regard to survival in 470 melanomas from The Cancer Genome Atlas. PD-L2 promoter methylation in correlation with PD-L2 mRNA and protein expression was analyzed in human melanoma cell lines. Prognostic and predictive value of PD-L2 methylation was validated using quantitative methylation-specific PCR in a multicenter cohort of 129 melanoma patients receiving anti-PD-1 therapy. mRNA sequencing data of 121 melanoma patients receiving anti-PD-1 therapy provided by Liu et al. were analyzed for PD-L2 mRNA expression.

RESULTS

We found significant correlations between PD-L2 methylation and mRNA expression levels in melanoma tissues and cell lines. Interferon-γ inducible PD-L2 protein expression correlated with PD-L2 promoter methylation in melanoma cells. PD-L2 DNA promoter hypomethylation and high mRNA expression were found to be strong predictors of prolonged overall survival. In pre-treatment melanoma samples from patients receiving anti-PD-1 therapy, low PD-L2 DNA methylation and high PD-L2 mRNA expression predicted longer progression-free survival.

CONCLUSION

PD-L2 expression seems to be regulated via DNA promoter methylation. PD-L2 DNA methylation and mRNA expression may predict progression-free survival in melanoma patients receiving anti-PD-1 immunotherapy. Assessment of PD-L2 should be included in further clinical trials with anti-PD-1 antibodies.

Methylation-Based Therapies for Colorectal Cancer

Colorectal carcinogenesis (CRC) is caused by the gradual long-term accumulation of both genetic and epigenetic changes. Recently, epigenetic alterations have been included in the classification of the CRC molecular subtype, and this points out their prognostic impact. As epigenetic modifications are reversible, they may represent relevant therapeutic targets. DNA methylation, catalyzed by DNA methyltransferases (DNMTs), regulates gene expression. For many years, the deregulation of DNA methylation has been considered to play a substantial part in CRC etiology and evolution. Despite considerable advances in CRC treatment, patient therapy response persists as limited, and their profit from systemic therapies are often hampered by the introduction of chemoresistance. In addition, inter-individual changes in therapy response in CRC patients can arise from their specific (epi)genetic compositions. In this review article, we summarize the options of CRC treatment based on DNA methylation status for their predictive value. This review also includes the therapy outcomes based on the patient's methylation status in CRC patients. In addition, the current challenge of research is to develop therapeutic inhibitors of DNMT. Based on the essential role of DNA methylation in CRC development, the application of DNMT inhibitors was recently proposed for the treatment of CRC patients, especially in patients with DNA hypermethylation.

Dynamic regulation of chromatin organizer SATB1 via TCR-induced alternative promoter switch during T-cell development

The chromatin organizer SATB1 is highly enriched in thymocytes and is essential for T-cell development. Although SATB1 regulates a large number of genes important for T-cell development, the mechanism(s) regulating expression of SATB1 during this process remain elusive. Using chromatin immune precipitation-seq-based occupancy profiles of H3K4me3 and H3Kme1 at Satb1 gene locus, we predicted four different alternative promoters of Satb1 in mouse thymocytes and characterized them. The expression of Satb1 transcript variants with distinct 5′ UTRs occurs in a stage-specific manner during T-cell development and is dependent on TCR signaling. The observed discrepancy between the expression levels of SATB1 mRNA and protein in developing thymocytes can be explained by the differential translatability of Satb1 transcript variants as confirmed by polysome profiling and in vitro translation assay. We show that Satb1 alternative promoters exhibit lineage-specific chromatin accessibility during T-cell development from progenitors. Furthermore, TCF1 regulates the Satb1 P2 promoter switch during CD4SP development, via direct binding to the Satb1 P2 promoter. CD4SP T cells from TCF1 KO mice exhibit downregulation of P2 transcript variant expression as well as low levels of SATB1 protein. Collectively, these results provide unequivocal evidence toward alternative promoter switch-mediated developmental stage-specific regulation of SATB1 in thymocytes.

Promoter DNA Hypermethylation and Paradoxical Gene Activation

DNA methylation is a stable epigenetic modification that contributes to the spatiotemporal regulation of gene expression. The manner in which DNA methylation contributes to transcriptional control is dependent on the biological context, including physiological state and the properties of the DNA itself. Classically, dense promoter DNA methylation is associated with transcriptional repression. However, growing evidence suggests that this association may not always hold true, and promoter hypermethylation now also appears to be associated with high transcriptional activity. Furthermore, in a selection of contexts, increasing levels of promoter methylation correlate directly with increased gene expression. These findings postulate a context-dependent model whereby epigenetic contributions to transcriptional regulation occur in a more complex and dynamic manner. We present current evidence documenting promoter hypermethylation and high levels of gene expression, offer insights into the possible mechanisms by which this occurs, and discuss the potential implications for both research and clinical applications.

Exploring 3D chromatin contacts in gene regulation: The evolution of approaches for the identification of functional enhancer-promoter interaction

Mechanisms underlying gene regulation are key to understand how multicellular organisms with various cell types develop from the same genetic blueprint. Dynamic interactions between enhancers and genes are revealed to play central roles in controlling gene transcription, but the determinants to link functional enhancer-promoter pairs remain elusive. A major challenge is the lack of reliable approach to detect and verify functional enhancer-promoter interactions (EPIs). In this review, we summarized the current methods for detecting EPIs and described how developing techniques facilitate the identification of EPI through assessing the merits and drawbacks of these methods. We also reviewed recent state-of-art EPI prediction methods in terms of their rationale, data usage and characterization. Furthermore, we briefly discussed the evolved strategies for validating functional EPIs.

Delineating the role of MITF isoforms in pigmentation and tissue homeostasis

MITF, a gene that is mutated in familial melanoma and Waardenburg syndrome, encodes multiple isoforms expressed from alternative promoters that share common coding exons but have unique amino termini. It is not completely understood how these isoforms influence pigmentation in different tissues and how the expression of these independent isoforms of MITF is regulated. Here, we show that melanocytes express two isoforms of MITF, MITF-A and MITF-M. The expression of MITF-A is partially regulated by a newly identified retinoid enhancer element located upstream of the MITF-A promoter. Mitf-A knockout mice have only subtle changes in melanin accumulation in the hair and reduced Tyr expression in the eye. In contrast, Mitf-M-null mice have enlarged kidneys, lack neural crest-derived melanocytes in the skin, choroid, and iris stroma, yet maintain pigmentation within the retinal pigment epithelium and iris pigment epithelium of the eye. Taken together, these studies identify a critical role for MITF-M in melanocytes, a minor role for MITF-A in regulating pigmentation in the hair and Tyr expression in the eye, and a novel role for MITF-M in size control of the kidney.

Evidence that alternative transcriptional initiation is largely nonadaptive

Alternative transcriptional initiation (ATI) refers to the frequent observation that one gene has multiple transcription start sites (TSSs). Although this phenomenon is thought to be adaptive, the specific advantage is rarely known. Here, we propose that each gene has one optimal TSS and that ATI arises primarily from imprecise transcriptional initiation that could be deleterious. This error hypothesis predicts that (i) the TSS diversity of a gene reduces with its expression level; (ii) the fractional use of the major TSS increases, but that of each minor TSS decreases, with the gene expression level; and (iii) cis-elements for major TSSs are selectively constrained, while those for minor TSSs are not. By contrast, the adaptive hypothesis does not make these predictions a priori. Our analysis of human and mouse transcriptomes confirms each of the three predictions. These and other findings strongly suggest that ATI predominantly results from molecular errors, requiring a major revision of our understanding of the precision and regulation of transcription.

The transcription of a gene may start from one of several transcription start sites, a phenomenon known as alternative transcriptional initiation. Contrary to common belief, this study shows that variation of the transcription start site of a given gene is nonadaptive and is largely attributable to transcriptional initiation error that is typically deleterious.

Multiple surveys of transcriptional initiation showed that mammalian genes typically have multiple transcription start sites such that transcription is initiated from any one of these sites. Many researchers believe that this phenomenon is adaptive because it allows production of multiple transcripts, from the same gene, that potentially vary in function or post-transcriptional regulation. Nevertheless, it is also possible that each gene has only one optimal transcription start site and that alternative transcriptional initiation arises primarily from molecular errors that are slightly deleterious. This error hypothesis makes a series of predictions about the amount of transcription start site diversity per gene, relative uses of the various start sites of a gene, among-tissue and across-species differences in start site usage, and the evolutionary conservation of cis-regulatory elements of various start sites, all of which are verified in our analyses of genome-wide transcription start site data from the human and mouse. These findings strongly suggest that alternative transcriptional initiation largely reflects molecular errors instead of molecular adaptations and require a rethink of the precision and regulation of transcription.

Revealing the alternative promoter usage of SAF/MAZ gene by bichromatic fluorescent reporter construct

The large-scale identification of putative alternative promoters study shows more than 52% of human genes are regulated by alternative promoters. The human myc-associated zinc finger protein (SAF/MAZ) gene have SAF-1 and SAF-3 variants transcripted from two transcription start sites (TSSs). By using SAF/MAZ promoter as a model, we set up an approach to probe how the alternative promoters are regulated in real time. We have constructed the bichromatic fluorescent reporter driven by SAF/MAZ 5'-proximal promoter plasmids from which transactivation status of SAF-1 and SAF-3 alternative promoter could be monitored by EGFP and DsRed expression respectively. The results showed that the SAF-3 expression is regulated by alternative promoters. When the bichromatic fluorescent reporter was driven by -1692/+277 or -1401/+277 SAF/MAZ promoter the dominant expression of SAF-3 would be observed in comparison with SAF-1 expression. We also identified that Elk-1 is an inhibitory transcription factor for SAF-3 expression. The temporal diversity of SAF-1 and SAF-3 expressions can be observed via bichromatic fluorescent reporters. These imply that the bichromatic fluorescent reporter driven by alternative promoter construct might be a useful tool for decoding the temporal regulatory repertoire of alternative promoter in human genes.

Molecular characterization of the giant freshwater prawn (Macrobrachium rosenbergii) beta-actin gene promoter

Constitutive promoters are important tools for gene function studies and transgenesis. The Beta-actin (actb1) gene promoter has been isolated from many species but remains to be cloned from the giant freshwater prawn (Macrobrachium rosenbergii). In this study, we cloned and characterized the Mractb1 gene promoter. Two alternative promoters were identified for the Mractb1 gene, which direct the generation of two transcripts with different 5′ untranslated regions. Three CpG islands were predicted in the upstream sequence, which are intimately related to transcription initiation and promoter activity. In addition to the CCAAT-box and the CArG-box, molecular dissection of the flanking sequence revealed the existence of one negative and two positive elements in the upstream region and the first intron. Finally, the Mractb1 promoter demonstrated comparative activity to the carp (Cyprinus carpio) actb1 promoter. Our investigations provide a valuable genetic tool for gene function studies and shed light on the regulation of the Mractb1 gene.

Human kallikrein-related peptidase 12 (KLK12) splice variants discriminate benign from cancerous breast tumors

OBJECTIVES

As kallikrein-related peptidase 12 (KLK12) has been implicated in the cancer progression and alternative splicing plays significant role in this disease, the aim of this study was to examine the expression profile and the clinical impact of the KLK12 splice variants in breast cancer.

DESIGN AND METHODS

Total RNA was isolated and reverse transcripted from 141 tissues. Afterwards, quantitative real-time PCR were conducted, followed by the performance of the comparative CT (2-ΔΔCT) method for relative quantification, whilst their correlation with the clinicopathological features of breast malignancies were assessed by statistical analysis.

RESULTS

Both KLK12sv1/2 and KLK12sv3 showed higher expression in non-cancerous than in cancerous samples. KLKsv1/2 (P = 0.001) upregulated and KLK12sv3 (P < 0.001) downregulated in the malignant compared to the benign tumors and their discriminative ability was verified by ROC curve analysis. Moreover, KLK12sv3 was associated with grade (P = 0.012) and hormonal receptor status (P = 0.001). Furthermore, Kaplan-Meier and Cox regression analyses showed that patients with positive KLK12sv1/2 and KLK12sv3 levels presented a significantly longer disease-free survival (P = 0.014 and P = 0.013, respectively) and overall survival (P = 0.062 and P = 0.004, respectively).

CONCLUSIONS

Our results demonstrate the discriminative value of KLK12sv1/2 and KLK12sv3 between benign and malignant breast tumors as well as their potential favorable prognostic significance in breast adenocarcinoma.

SEASTAR: systematic evaluation of alternative transcription start sites in RNA

Alternative first exons diversify the transcriptomes of eukaryotes by producing variants of the 5' Untranslated Regions (5'UTRs) and N-terminal coding sequences. Accurate transcriptome-wide detection of alternative first exons typically requires specialized experimental approaches that are designed to identify the 5' ends of transcripts. We developed a computational pipeline SEASTAR that identifies first exons from RNA-seq data alone then quantifies and compares alternative first exon usage across multiple biological conditions. The exons inferred by SEASTAR coincide with transcription start sites identified directly by CAGE experiments and bear epigenetic hallmarks of active promoters. To determine if differential usage of alternative first exons can yield insights into the mechanism controlling gene expression, we applied SEASTAR to an RNA-seq dataset that tracked the reprogramming of mouse fibroblasts into induced pluripotent stem cells. We observed dynamic temporal changes in the usage of alternative first exons, along with correlated changes in transcription factor expression. Using a combined sequence motif and gene set enrichment analysis we identified N-Myc as a regulator of alternative first exon usage in the pluripotent state. Our results demonstrate that SEASTAR can leverage the available RNA-seq data to gain insights into the control of gene expression and alternative transcript variation in eukaryotic transcriptomes.