Genome-wide analysis of mammalian promoter architecture and evolution

Core promoterome of barley embryo

Precise localization and dissection of gene promoters are key to understanding transcriptional gene regulation and to successful bioengineering applications. The core RNA polymerase II initiation machinery is highly conserved among eukaryotes, leading to a general expectation of equivalent underlying mechanisms. Still, less is known about promoters in the plant kingdom. In this study, we employed cap analysis of gene expression (CAGE) at three embryonic developmental stages in barley to accurately map, annotate, and quantify transcription initiation events. Unsupervised discovery of de novo sequence clusters grouped promoters based on characteristic initiator and position-specific core-promoter motifs. This grouping was complemented by the annotation of transcription factor binding site (TFBS) motifs. Integration with genome-wide epigenomic data sets and gene ontology (GO) enrichment analysis further delineated the chromatin environments and functional roles of genes associated with distinct promoter categories. The TATA-box presence governs all features explored, supporting the general model of two separate genomic regulatory environments. We describe the extent and implications of alternative transcription initiation events, including those that are specific to developmental stages, which can affect the protein sequence or the presence of regions that regulate translation. The generated promoterome dataset provides a valuable genomic resource for enhancing the functional annotation of the barley genome. It also offers insights into the transcriptional regulation of individual genes and presents opportunities for the informed manipulation of promoter architecture, with the aim of enhancing traits of agronomic importance.

Engineered mouse H1 promoter mutants with superior RNA polymerase III activity

Vectors incorporating the human H1 (hH1) promoter are being applied for RNA interference (RNAi) experiments and genome editing. Although extensive studies have been conducted on the hH1 promoter, our understanding of the mouse H1 promoter remains limited. In this study, we predicted the 163 bp mouse H1 (mH1) promoter and 84 bp mouse H1 core (mH1 core) promoter through global alignment and detected its RNA polymerase II (Pol II) and III activities through the expression of the EGFP and the abundance of artificial sequence, which were generally slightly weaker than those of the hH1 promoter. Furthermore, to boost its Pol III activity, we engineered various promoter mutants by introducing mutations or systematically swapping elements. Surprisingly, the Pol II activity of mH1 core mut5 with AT stretch was at least 2-fold greater than that of the wild type, making it a potential candidate for target protein expression purposes. Fortunately, the Pol III activities of mH1 mut1 and mH1 core mut5 were at least 1.5 times stronger than those of the parental promoters in human and mouse cell lines on account of AT stretch, as did the mH1 mut4 with AT stretch and proximal sequence element (PSE) and TATA box insertion mutations. We highly recommend these three promoters as valuable supplements to the type 3 Pol III promoter toolbox.

Cut from the same cloth: RNAs transcribed from regulatory elements

A certain degree of chromatin openness is necessary for the activity of transcription-regulating regions within the genome, facilitating accessibility to RNA polymerases and subsequent synthesis of regulatory element RNAs (regRNAs) from these regions. The rapidly increasing number of studies underscores the significance of regRNAs across diverse cellular processes and diseases, challenging the paradigm that these transcripts are non-functional transcriptional noise. This review explores the multifaceted roles of regRNAs in human cells, encompassing rather well-studied entities such as promoter RNAs and enhancer RNAs (eRNAs), while also providing insights into overshadowed silencer RNAs and insulator RNAs. Furthermore, we assess notable examples of shorter regRNAs, like miRNAs, snRNAs, and snoRNAs, playing important roles. Expanding our discourse, we deliberate on the potential usage of regRNAs as biomarkers and novel targets for cancer and other human diseases.

The GC-content at the 5' ends of human protein-coding genes is undergoing mutational decay

BACKGROUND

In vertebrates, most protein-coding genes have a peak of GC-content near their 5' transcriptional start site (TSS). This feature promotes both the efficient nuclear export and translation of mRNAs. Despite the importance of GC-content for RNA metabolism, its general features, origin, and maintenance remain mysterious. We investigate the evolutionary forces shaping GC-content at the transcriptional start site (TSS) of genes through both comparative genomic analysis of nucleotide substitution rates between different species and by examining human de novo mutations.

RESULTS

Our data suggests that GC-peaks at TSSs were present in the last common ancestor of amniotes, and likely that of vertebrates. We observe that in apes and rodents, where recombination is directed away from TSSs by PRDM9, GC-content at the 5' end of protein-coding gene is currently undergoing mutational decay. In canids, which lack PRDM9 and perform recombination at TSSs, GC-content at the 5' end of protein-coding is increasing. We show that these patterns extend into the 5' end of the open reading frame, thus impacting synonymous codon position choices.

CONCLUSIONS

Our results indicate that the dynamics of this GC-peak in amniotes is largely shaped by historic patterns of recombination. Since decay of GC-content towards the mutation rate equilibrium is the default state for non-functional DNA, the observed decrease in GC-content at TSSs in apes and rodents indicates that the GC-peak is not being maintained by selection on most protein-coding genes in those species.

ZC3H4/Restrictor Exerts a Stranglehold on Pervasive Transcription

The regulation of transcription by RNA polymerase II (RNAPII) underpins all cellular processes and is perturbed in thousands of diseases. In humans, RNAPII transcribes ∼20000 protein-coding genes and engages in apparently futile non-coding transcription at thousands of other sites. Despite being so ubiquitous, this transcription is usually attenuated soon after initiation and the resulting products are immediately degraded by the nuclear exosome. We and others have recently described a new complex, "Restrictor", which appears to control such unproductive transcription. Underpinned by the RNA binding protein, ZC3H4, Restrictor curtails unproductive/pervasive transcription genome-wide. Here, we discuss these recent discoveries and speculate on some of the many unknowns regarding Restrictor function and mechanism.

Regional-specific calibration enables application of computational evidence for clinical classification of 5' cis-regulatory variants in Mendelian disease

To date, clinical genetic testing for Mendelian disease variants has focused heavily on exonic coding and intronic gene regions. This multi-step study was undertaken to provide an evidence base for selecting and applying computational approaches for use in clinical classification of 5' cis-regulatory region variants. Curated datasets of clinically reported disease-causing 5' cis-regulatory region variants and variants from matched genomic regions in population controls were used to calibrate six bioinformatic tools as predictors of variant pathogenicity. Likelihood ratio estimates were aligned to code weights following ClinGen recommendations for application of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) classification scheme. Considering code assignment across all reference dataset variants, performance was best for CADD (81.2%) and REMM (81.5%). Optimized thresholds provided moderate evidence toward pathogenicity (CADD, REMM) and moderate (CADD) or supporting (REMM) evidence against pathogenicity. Both sensitivity and specificity of prediction were improved when further categorizing variants based on location in an EPDnew-defined promoter region. Combining predictions (CADD, REMM, and location in a promoter region) increased specificity at the expense of sensitivity. Importantly, the optimal CADD thresholds for assigning ACMG/AMP codes PP3 (≥10) and BP4 (≤8) were vastly different from recommendations for protein-coding variants (PP3 ≥25.3; BP4 ≤22.7); CADD <22.7 would incorrectly assign BP4 for >90% of reported disease-causing cis-regulatory region variants. Our results demonstrate the need to consider a tiered approach and tailored score thresholds to optimize bioinformatic impact prediction for clinical classification of 5' cis-regulatory region variants.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

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

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

Alternative TSS use is widespread in Cryptococcus fungi in response to environmental cues and regulated genome-wide by the transcription factor Tur1

Alternative transcription start site (TSS) usage regulation has been identified as a major means of gene expression regulation in metazoans. However, in fungi, its impact remains elusive as its study has thus far been restricted to model yeasts. Here, we first re-analyzed TSS-seq data to define genuine TSS clusters in 2 species of pathogenic Cryptococcus. We identified 2 types of TSS clusters associated with specific DNA sequence motifs. Our analysis also revealed that alternative TSS usage regulation in response to environmental cues is widespread in Cryptococcus, altering gene expression and protein targeting. Importantly, we performed a forward genetic screen to identify a unique transcription factor (TF) named Tur1, which regulates alternative TSS (altTSS) usage genome-wide when cells switch from exponential phase to stationary phase. ChiP-Seq and DamID-Seq analyses suggest that at some loci, the role of Tur1 might be direct. Tur1 has been previously shown to be essential for virulence in C. neoformans. We demonstrated here that a tur1Δ mutant strain is more sensitive to superoxide stress and phagocytosed more efficiently by macrophages than the wild-type (WT) strain.

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.

CapTrap-seq: a platform-agnostic and quantitative approach for high-fidelity full-length RNA sequencing

Long-read RNA sequencing is essential to produce accurate and exhaustive annotation of eukaryotic genomes. Despite advancements in throughput and accuracy, achieving reliable end-to-end identification of RNA transcripts remains a challenge for long-read sequencing methods. To address this limitation, we develop CapTrap-seq, a cDNA library preparation method, which combines the Cap-trapping strategy with oligo(dT) priming to detect 5' capped, full-length transcripts. In our study, we evaluate the performance of CapTrap-seq alongside other widely used RNA-seq library preparation protocols in human and mouse tissues, employing both ONT and PacBio sequencing technologies. To explore the quantitative capabilities of CapTrap-seq and its accuracy in reconstructing full-length RNA molecules, we implement a capping strategy for synthetic RNA spike-in sequences that mimics the natural 5'cap formation. Our benchmarks, incorporating the Long-read RNA-seq Genome Annotation Assessment Project (LRGASP) data, demonstrate that CapTrap-seq is a competitive, platform-agnostic RNA library preparation method for generating full-length transcript sequences.

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.

Methods for Functional Characterization of Genetic Polymorphisms of Non-Coding Regulatory Regions of the Human Genome

Currently, numerous associations between genetic polymorphisms and various diseases have been characterized through the Genome-Wide Association Studies. Majority of the clinically significant polymorphisms are localized in non-coding regions of the genome. While modern bioinformatic resources make it possible to predict molecular mechanisms that explain influence of the non-coding polymorphisms on gene expression, such hypotheses require experimental verification. This review discusses the methods for elucidating molecular mechanisms underlying dependence of the disease pathogenesis on specific genetic variants within the non-coding sequences. A particular focus is on the methods for identification of transcription factors with binding efficiency dependent on polymorphic variations. Despite remarkable progress in bioinformatic resources enabling prediction of the impact of polymorphisms on the disease pathogenesis, there is still the need for experimental approaches to investigate this issue.

Lineage-determining transcription factor-driven promoters regulate cell type-specific macrophage gene expression

Mammalian promoters consist of multifarious elements, which make them unique and support the selection of the proper transcript variants required under diverse conditions in distinct cell types. However, their direct DNA-transcription factor (TF) interactions are mostly unidentified. Murine bone marrow-derived macrophages (BMDMs) are a widely used model for studying gene expression regulation. Thus, this model serves as a rich source of various next-generation sequencing data sets, including a large number of TF cistromes. By processing and integrating the available cistromic, epigenomic and transcriptomic data from BMDMs, we characterized the macrophage-specific direct DNA-TF interactions, with a particular emphasis on those specific for promoters. Whilst active promoters are enriched for certain types of typically methylatable elements, more than half of them contain non-methylatable and prototypically promoter-distal elements. In addition, circa 14% of promoters-including that of Csf1r-are composed exclusively of 'distal' elements that provide cell type-specific gene regulation by specialized TFs. Similar to CG-rich promoters, these also contain methylatable CG sites that are demethylated in a significant portion and show high polymerase activity. We conclude that this unusual class of promoters regulates cell type-specific gene expression in macrophages, and such a mechanism might exist in other cell types too.

How DNA encodes the start of transcription

A deep-learning model reveals the rules that define transcription initiation.

Sequence basis of transcription initiation in the human genome

Transcription initiation is a process that is essential to ensuring the proper function of any gene, yet we still lack a unified understanding of sequence patterns and rules that explain most transcription start sites in the human genome. By predicting transcription initiation at base-pair resolution from sequences with a deep learning-inspired explainable model called Puffin, we show that a small set of simple rules can explain transcription initiation at most human promoters. We identify key sequence patterns that contribute to human promoter activity, each activating transcription with distinct position-specific effects. Furthermore, we explain the sequence basis of bidirectional transcription at promoters, identify the links between promoter sequence and gene expression variation across cell types, and explore the conservation of sequence determinants of transcription initiation across mammalian species.

Natural antisense transcripts as versatile regulators of gene expression

Long non-coding RNAs (lncRNAs) are emerging as a major class of gene products that have central roles in cell and developmental biology. Natural antisense transcripts (NATs) are an important subset of lncRNAs that are expressed from the opposite strand of protein-coding and non-coding genes and are a genome-wide phenomenon in both eukaryotes and prokaryotes. In eukaryotes, a myriad of NATs participate in regulatory pathways that affect expression of their cognate sense genes. Recent developments in the study of NATs and lncRNAs and large-scale sequencing and bioinformatics projects suggest that whether NATs regulate expression, splicing, stability or translation of the sense transcript is influenced by the pattern and degrees of overlap between the sense-antisense pair. Moreover, epigenetic gene regulatory mechanisms prevail in somatic cells whereas mechanisms dependent on the formation of double-stranded RNA intermediates are prevalent in germ cells. The modulating effects of NATs on sense transcript expression make NATs rational targets for therapeutic interventions.

Dissection of core promoter syntax through single nucleotide resolution modeling of transcription initiation

Our understanding of how the DNA sequences of cis-regulatory elements encode transcription initiation patterns remains limited. Here we introduce CLIPNET, a deep learning model trained on population-scale PRO-cap data that accurately predicts the position and quantity of transcription initiation with single nucleotide resolution from DNA sequence. Interpretation of CLIPNET revealed a complex regulatory syntax consisting of DNA-protein interactions in five major positions between -200 and +50 bp relative to the transcription start site, as well as more subtle positional preferences among different transcriptional activators. Transcriptional activator and core promoter motifs occupy different positions and play distinct roles in regulating initiation, with the former driving initiation quantity and the latter initiation position. We identified core promoter motifs that explain initiation patterns in the majority of promoters and enhancers, including DPR motifs and AT-rich TBP binding sequences in TATA-less promoters. Our results provide insights into the sequence architecture governing transcription initiation.

OVCH1 Antisense RNA 1 is differentially expressed between non-frail and frail old adults

While some old adults stay healthy and non-frail up to late in life, others experience multimorbidity and frailty often accompanied by a pro-inflammatory state. The underlying molecular mechanisms for those differences are still obscure. Here, we used gene expression analysis to understand the molecular underpinning between non-frail and frail individuals in old age. Twenty-four adults (50% non-frail and 50% frail) from InCHIANTI study were included. Total RNA extracted from whole blood was analyzed by Cap Analysis of Gene Expression (CAGE). CAGE identified transcription start site (TSS) and active enhancer regions. We identified a set of differentially expressed (DE) TSS and enhancer between non-frail and frail and male and female participants. Several DE TSSs were annotated as lncRNA (XIST and TTTY14) and antisense RNAs (ZFX-AS1 and OVCH1 Antisense RNA 1). The promoter region chr6:366,786,54-366,787,97;+ was DE and overlapping the longevity CDKN1A gene. GWAS-LD enrichment analysis identifies overlapping LD-blocks with the DE regions with reported traits in GWAS catalog (isovolumetric relaxation time and urinary tract infection frequency). Furthermore, we used weighted gene co-expression network analysis (WGCNA) to identify changes of gene expression associated with clinical traits and identify key gene modules. We performed functional enrichment analysis of the gene modules with significant trait/module correlation. One gene module is showing a very distinct pattern in hub genes. Glycogen Phosphorylase L (PYGL) was the top ranked hub gene between non-frail and frail. We predicted transcription factor binding sites (TFBS) and motif activity. TF involved in age-related pathways (e.g., FOXO3 and MYC) shows different expression patterns between non-frail and frail participants. Expanding the study of OVCH1 Antisense RNA 1 and PYGL may help understand the mechanisms leading to loss of homeostasis that ultimately causes frailty.

A single-cell genomic strategy for alternative transcript start sites identification

Alternative transcription start sites (TSSs) usage plays a critical role in gene transcription regulation in mammals. However, precisely identifying alternative TSSs remains challenging at the genome-wide level. We report a single-cell genomic technology for alternative TSSs annotation and cell heterogeneity detection. In the method, we utilize Fluidigm C1 system to capture individual cells of interest, SMARTer cDNA synthesis kit to recover full-length cDNAs, then dual priming oligonucleotide system to specifically enrich TSSs for genomic analysis. We apply this method to a genome-wide study of alternative TSSs identification in two different IFN-β stimulated mouse embryonic fibroblasts (MEFs). The data clearly discriminate two IFN-β stimulated MEFs. Moreover, our results indicate 81% expressed genes in these two cell types containing multiple TSSs, which is much higher than previous predictions based on Cap-Analysis Gene Expression (CAGE) (58%) or empirical determination (54%) in various cell types. This indicates that alternative TSSs are more pervasive than expected and implies our strategy could position them at an unprecedented sensitivity. It would be helpful for elucidating their biological insights in future.

Pre-piRNA trimming safeguards piRNAs against erroneous targeting by RNA-dependent RNA polymerase

The Piwi/Piwi-interacting RNA (piRNA) pathway protects genome integrity in animal germ lines. Maturation of piRNAs involves nucleolytic processing at both 5' and 3' ends. The ribonuclease PARN-1 and its orthologs mediate piRNA 3' trimming in worms, insects, and mammals. However, the significance of this evolutionarily conserved processing step is not fully understood. Employing C. elegans as a model, we recently discovered that 3' trimming protects piRNAs against non-templated nucleotide additions and degradation. Here, we find that worms lacking PARN-1 accumulate an uncharacterized RNA species termed anti-piRNAs, which are antisense to piRNAs. Anti-piRNAs associate with Piwi proteins, are 17-19 nucleotides long, and begin with 5' guanine or adenine. Untrimmed pre-piRNAs are misdirected by the terminal nucleotidyltransferase RDE-3 and RNA-dependent RNA polymerase EGO-1, leading to the formation of anti-piRNAs. This work identifies a class of small RNAs in parn-1 mutants and provides insight into the activities of RDE-3, EGO-1, and Piwi proteins.

Distinct roles of LARP1 and 4EBP1/2 in regulating translation and stability of 5'TOP mRNAs

A central mechanism of mTOR complex 1 (mTORC1) signaling is the coordinated translation of ribosomal protein and translation factor mRNAs mediated by the 5'-terminal oligopyrimidine motif (5'TOP). Recently, La-related protein 1 (LARP1) was proposed to be the specific regulator of 5'TOP mRNA translation downstream of mTORC1, while eIF4E-binding proteins (4EBP1/2) were suggested to have a general role in translational repression of all transcripts. Here, we use single-molecule translation site imaging of 5'TOP and canonical mRNAs to study the translation of single mRNAs in living cells. Our data reveal that 4EBP1/2 has a dominant role in repression of translation of both 5'TOP and canonical mRNAs during pharmacological inhibition of mTOR. In contrast, we find that LARP1 selectively protects 5'TOP mRNAs from degradation in a transcriptome-wide analysis of mRNA half-lives. Our results clarify the roles of 4EBP1/2 and LARP1 in regulating 5'TOP mRNAs and provide a framework to further study how these factors control cell growth during development and disease.

Interplay between the transcription preinitiation complex and the +1 nucleosome

Eukaryotic transcription starts with the assembly of a preinitiation complex (PIC) on core promoters. Flanking this region is the +1 nucleosome, the first nucleosome downstream of the core promoter. While this nucleosome is rich in epigenetic marks and plays a key role in transcription regulation, how the +1 nucleosome interacts with the transcription machinery has been a long-standing question. Here, we summarize recent structural and functional studies of the +1 nucleosome in complex with the PIC. We specifically focus on how differently organized promoter-nucleosome templates affect the assembly of the PIC and PIC-Mediator on chromatin and result in distinct transcription initiation.

Molecular analysis of the human cytoglobin mRNA isoforms

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

Quantitative analysis of transcription start site selection reveals control by DNA sequence, RNA polymerase II activity and NTP levels

Transcription start site (TSS) selection is a key step in gene expression and occurs at many promoter positions over a wide range of efficiencies. Here we develop a massively parallel reporter assay to quantitatively dissect contributions of promoter sequence, nucleoside triphosphate substrate levels and RNA polymerase II (Pol II) activity to TSS selection by 'promoter scanning' in Saccharomyces cerevisiae (Pol II MAssively Systematic Transcript End Readout, 'Pol II MASTER'). Using Pol II MASTER, we measure the efficiency of Pol II initiation at 1,000,000 individual TSS sequences in a defined promoter context. Pol II MASTER confirms proposed critical qualities of S. cerevisiae TSS -8, -1 and +1 positions, quantitatively, in a controlled promoter context. Pol II MASTER extends quantitative analysis to surrounding sequences and determines that they tune initiation over a wide range of efficiencies. These results enabled the development of a predictive model for initiation efficiency based on sequence. We show that genetic perturbation of Pol II catalytic activity alters initiation efficiency mostly independently of TSS sequence, but selectively modulates preference for the initiating nucleotide. Intriguingly, we find that Pol II initiation efficiency is directly sensitive to guanosine-5'-triphosphate levels at the first five transcript positions and to cytosine-5'-triphosphate and uridine-5'-triphosphate levels at the second position genome wide. These results suggest individual nucleoside triphosphate levels can have transcript-specific effects on initiation, representing a cryptic layer of potential regulation at the level of Pol II biochemical properties. The results establish Pol II MASTER as a method for quantitative dissection of transcription initiation in eukaryotes.

RNA-based Therapeutics: Past, Present and Future Prospects, Challenges in Cancer Treatment

It is becoming more and harder in today's climate to disregard the impact of cancer on social health. Even though a significant amount of money is spent annually on cancer research, it still ranks as the second leading cause of death worldwide. Additionally, only about half of the patients suffering from complex forms of cancer survive a year after receiving traditional cancer therapies. A method for silencing genes is called RNA interference (RNAi). Such a method is very effective in focusing on genes linked to cancer. Most gene products implicated in cancer have recently been used as RNA interference (RNAi) therapeutic targets. According to the findings from this research, RNAi application is necessary for today's cancer treatment to target functioning carcinogenic molecules and tumor resistance to chemotherapy and radiation. Proapoptotic and antiproliferative activity has been reported from previous research studies on cell culture systems, animal models, and clinical trials through the knockdown of gene products from RNAi technology. Numerous novel RNAi-based medications are now in the clinical trial stages thanks to the discovery of the RNAi mechanism and advancements in the area. In the future, genomic-based personalized medicines can be developed through this RNAi therapy. Hopefully, cancer sufferers will find this sort of therapy to be one of the most effective ones. Various kinds of RNA-based treatments, such as aptamers, small interfering RNAs, microRNAs, antisense oligonucleotides, and messenger RNA, are covered in broad terms in this study. We also present an overview of the RNA-based therapies that have received regulatory approval in the past or are now undergoing clinical studies.

Cis-Regulatory Elements in Mammals

In cis-regulatory elements, enhancers and promoters with complex molecular interactions are used to coordinate gene transcription through physical proximity and chemical modifications. These processes subsequently influence the phenotypic characteristics of an organism. An in-depth exploration of enhancers and promoters can substantially enhance our understanding of gene regulatory networks, shedding new light on mammalian development, evolution and disease pathways. In this review, we provide a comprehensive overview of the intrinsic structural attributes, detection methodologies as well as the operational mechanisms of enhancers and promoters, coupled with the relevant novel and innovative investigative techniques used to explore their actions. We further elucidated the state-of-the-art research on the roles of enhancers and promoters in the realms of mammalian development, evolution and disease, and we conclude with forward-looking insights into prospective research avenues.

Pan-cancer analysis of promoter activity quantitative trait loci

Altered promoter activity has been generally observed in diverse biological processes, including tumorigenesis. Accumulating evidence suggests that employing a quantitative trait locus mapping approach is effective in comprehending the genetic basis of promoter activity. By utilizing genotype data from The Cancer Genome Atlas and calculating corresponding promoter activity values using proActiv, we systematically evaluated the impact of genetic variants on promoter activity and identified >1.0 million promoter activity quantitative trait loci (paQTLs) as both cis- and trans-acting. Additionally, leveraging data from the genome-wide association study (GWAS) catalog, we discovered >1.3 million paQTLs that overlap with known GWAS linkage disequilibrium regions. Remarkably, ∼9324 paQTLs exhibited significant associations with patient prognosis. Moreover, investigating the impact of promoter activity on >1000 imputed antitumor therapy responses among pan-cancer patients revealed >43 000 million significant associations. Furthermore, ∼25 000 significant associations were identified between promoter activity and immune cell abundance. Finally, a user-friendly data portal, Pancan-paQTL (https://www.hbpding.com/PancanPaQTL/), was constructed for users to browse, search and download data of interest. Pancan-paQTL serves as a comprehensive multidimensional database, enabling functional and clinical investigations into genetic variants associated with promoter activity, drug responses and immune infiltration across multiple cancer types.

SARS-CoV-2 hijacks a cell damage response, which induces transcription of a more efficient Spike S-acyltransferase

SARS-CoV-2 infection requires Spike protein-mediated fusion between the viral and cellular membranes. The fusogenic activity of Spike depends on its post-translational lipid modification by host S-acyltransferases, predominantly ZDHHC20. Previous observations indicate that SARS-CoV-2 infection augments the S-acylation of Spike when compared to mere Spike transfection. Here, we find that SARS-CoV-2 infection triggers a change in the transcriptional start site of the zdhhc20 gene, both in cells and in an in vivo infection model, resulting in a 67-amino-acid-long N-terminally extended protein with approx. 40 times higher Spike acylating activity, resulting in enhanced fusion of viruses with host cells. Furthermore, we observed the same induced transcriptional change in response to other challenges, such as chemically induced colitis and pore-forming toxins, indicating that SARS-CoV-2 hijacks an existing cell damage response pathway to optimize it fusion glycoprotein.

Identifying promoter sequence architectures via a chunking-based algorithm using non-negative matrix factorisation

Core promoters are stretches of DNA at the beginning of genes that contain information that facilitates the binding of transcription initiation complexes. Different functional subsets of genes have core promoters with distinct architectures and characteristic motifs. Some of these motifs inform the selection of transcription start sites (TSS). By discovering motifs with fixed distances from known TSS positions, we could in principle classify promoters into different functional groups. Due to the variability and overlap of architectures, promoter classification is a difficult task that requires new approaches. In this study, we present a new method based on non-negative matrix factorisation (NMF) and the associated software called seqArchR that clusters promoter sequences based on their motifs at near-fixed distances from a reference point, such as TSS. When combined with experimental data from CAGE, seqArchR can efficiently identify TSS-directing motifs, including known ones like TATA, DPE, and nucleosome positioning signal, as well as novel lineage-specific motifs and the function of genes associated with them. By using seqArchR on developmental time courses, we reveal how relative use of promoter architectures changes over time with stage-specific expression. seqArchR is a powerful tool for initial genome-wide classification and functional characterisation of promoters. Its use cases are more general: it can also be used to discover any motifs at near-fixed distances from a reference point, even if they are present in only a small subset of sequences.

Causal role of a promoter polymorphism in natural variation of the Arabidopsis floral repressor gene FLC

Noncoding polymorphism frequently associates with phenotypic variation, but causation and mechanism are rarely established. Noncoding single-nucleotide polymorphisms (SNPs) characterize the major haplotypes of the Arabidopsis thaliana floral repressor gene FLOWERING LOCUS C (FLC). This noncoding polymorphism generates a range of FLC expression levels, determining the requirement for and the response to winter cold. The major adaptive determinant of these FLC haplotypes was shown to be the autumnal levels of FLC expression. Here, we investigate how noncoding SNPs influence FLC transcriptional output. We identify an upstream transcription start site (uTSS) cluster at FLC, whose usage is increased by an A variant at the promoter SNP-230. This variant is present in relatively few Arabidopsis accessions, with the majority containing G at this site. We demonstrate a causal role for the A variant at -230 in reduced FLC transcriptional output. The G variant upregulates FLC expression redundantly with the major transcriptional activator FRIGIDA (FRI). We demonstrate an additive interaction of SNP-230 with an intronic SNP+259, which also differentially influences uTSS usage. Combinatorial interactions between noncoding SNPs and transcriptional activators thus generate quantitative variation in FLC transcription that has facilitated the adaptation of Arabidopsis accessions to distinct climates.

Chromatin analysis of adult pluripotent stem cells reveals a unique stemness maintenance strategy

Many highly regenerative organisms maintain adult pluripotent stem cells throughout their life, but how the long-term maintenance of pluripotency is accomplished is unclear. To decipher the regulatory logic of adult pluripotent stem cells, we analyzed the chromatin organization of stem cell genes in the planarian Schmidtea mediterranea. We identify a special chromatin state of stem cell genes, which is distinct from that of tissue-specific genes and resembles constitutive genes. Where tissue-specific promoters have detectable transcription factor binding sites, the promoters of stem cell-specific genes instead have sequence features that broadly decrease nucleosome binding affinity. This genic organization makes pluripotency-related gene expression the default state in these cells, which is maintained by the activity of chromatin remodelers ISWI and SNF2 in the stem cells.

Transcription factor IID parks and drives preinitiation complexes at sharp or broad promoters

Core promoters are sites where transcriptional regulatory inputs of a gene are integrated to direct the assembly of the preinitiation complex (PIC) and RNA polymerase II (Pol II) transcription output. Until now, core promoter functions have been investigated by distinct methods, including Pol II transcription initiation site mappings and structural characterization of PICs on distinct promoters. Here, we bring together these previously unconnected observations and hypothesize how, on metazoan TATA promoters, the precisely structured building up of transcription factor (TF) IID-based PICs results in sharp transcription start site (TSS) selection; or, in contrast, how the less strictly controlled positioning of the TATA-less promoter DNA relative to TFIID-core PIC components results in alternative broad TSS selections by Pol II.

Transcription start site heterogeneity and its role in RNA fate determination distinguish HIV-1 from other retroviruses and are mediated by core promoter elements

HIV-1 uses heterogeneous transcription start sites (TSSs) to generate two RNA 5´ isoforms that adopt radically different structures and perform distinct replication functions. Although these RNAs differ in length by only two bases, exclusively, the shorter RNA is encapsidated while the longer RNA is excluded from virions and provides intracellular functions. The current study examined TSS usage and packaging selectivity for a broad range of retroviruses and found that heterogeneous TSS usage was a conserved feature of all tested HIV-1 strains, but all other retroviruses examined displayed unique TSSs. Phylogenetic comparisons and chimeric viruses' properties provided evidence that this mechanism of RNA fate determination was an innovation of the HIV-1 lineage, with determinants mapping to core promoter elements. Fine-tuning differences between HIV-1 and HIV-2, which uses a unique TSS, implicated purine residue positioning plus a specific TSS-adjacent dinucleotide in specifying multiplicity of TSS usage. Based on these findings, HIV-1 expression constructs were generated that differed from the parental strain by only two point mutations yet each expressed only one of HIV-1's two RNAs. Replication defects of the variant with only the presumptive founder TSS were less severe than those for the virus with only the secondary start site. IMPORTANCE Retroviruses use RNA both to encode their proteins and to serve in place of DNA as their genomes. A recent surprising discovery was that the genomic RNAs and messenger RNAs of HIV-1 are not identical but instead differ subtly on one of their ends. These differences enable the functional separation of HIV-1 RNAs into genome and messenger roles. In this report, we examined a broad collection of HIV-1-related viruses and discovered that each produced only one end class of RNA, and thus must differ from HIV-1 in how they specify RNA fates. By comparing regulatory signals, we generated virus variants that pinpointed the determinants of HIV-1 RNA fates, as well as HIV-1 variants that produced only one or the other functional class of RNA. Competition and replication assays confirmed that HIV-1 has evolved to rely on the coordinated actions of both its RNA forms.

Pre-piRNA trimming safeguards piRNAs against erroneous targeting by RNA-dependent RNA Polymerase

In animal germ lines, The Piwi/piRNA pathway plays a crucial role in safeguarding genome integrity and promoting fertility. Following transcription from discrete genomic loci, piRNA precursors undergo nucleolytic processing at both 5' and 3' ends. The ribonuclease PARN-1 and its orthologs mediate piRNA 3' trimming in worms, insects and mammals. Yet, the significance of this evolutionarily conserved processing step is not well understood. Employing C. elegans as a model organism, our recent work has demonstrated that 3' trimming protects piRNAs against non-templated nucleotide additions and degradation. In this study, we present an unexpected finding that C. elegans deficient for PARN-1 accumulate a heretofore uncharacterized RNA species termed anti-piRNAs, which are antisense to piRNAs. These anti-piRNAs associate with Piwi proteins and display the propensity for a length of 17-19 nucleotides and 5' guanine and adenine residues. We show that untrimmed pre-piRNAs in parn-1 mutants are modified by the terminal nucleotidyltransferase RDE-3 and erroneously targeted by the RNA-dependent RNA polymerase EGO-1, thereby giving rise to anti-piRNAs. Taken together, our work identifies a previously unknown class of small RNAs upon loss of parn-1 and provides mechanistic insight to activities of RDE-3, EGO-1 and Piwi proteins.

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

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

Flipons and small RNAs accentuate the asymmetries of pervasive transcription by the reset and sequence-specific microcoding of promoter conformation

The role of alternate DNA conformations such as Z-DNA in the regulation of transcription is currently underappreciated. These structures are encoded by sequences called flipons, many of which are enriched in promoter and enhancer regions. Through a change in their conformation, flipons provide a tunable mechanism to mechanically reset promoters for the next round of transcription. They act as actuators that capture and release energy to ensure that the turnover of the proteins at promoters is optimized to cell state. Likewise, the single-stranded DNA formed as flipons cycle facilitates the docking of RNAs that are able to microcode promoter conformations and canalize the pervasive transcription commonly observed in metazoan genomes. The strand-specific nature of the interaction between RNA and DNA likely accounts for the known asymmetry of epigenetic marks present on the histone tetramers that pair to form nucleosomes. The role of these supercoil-dependent processes in promoter choice and transcriptional interference is reviewed. The evolutionary implications are examined: the resilience and canalization of flipon-dependent gene regulation is contrasted with the rapid adaptation enabled by the spread of flipon repeats throughout the genome. Overall, the current findings underscore the important role of flipons in modulating the readout of genetic information and how little we know about their biology.

Locus-resolution analysis of L1 regulation and retrotransposition potential in mouse embryonic development

Mice harbor ∼2800 intact copies of the retrotransposon Long Interspersed Element 1 (L1). The in vivo retrotransposition capacity of an L1 copy is defined by both its sequence integrity and epigenetic status, including DNA methylation of the monomeric units constituting young mouse L1 promoters. Locus-specific L1 methylation dynamics during development may therefore elucidate and explain spatiotemporal niches of endogenous retrotransposition but remain unresolved. Here, we interrogate the retrotransposition efficiency and epigenetic fate of source (donor) L1s, identified as mobile in vivo. We show that promoter monomer loss consistently attenuates the relative retrotransposition potential of their offspring (daughter) L1 insertions. We also observe that most donor/daughter L1 pairs are efficiently methylated upon differentiation in vivo and in vitro. We use Oxford Nanopore Technologies (ONT) long-read sequencing to resolve L1 methylation genome-wide and at individual L1 loci, revealing a distinctive "smile" pattern in methylation levels across the L1 promoter region. Using Pacific Biosciences (PacBio) SMRT sequencing of L1 5' RACE products, we then examine DNA methylation dynamics at the mouse L1 promoter in parallel with transcription start site (TSS) distribution at locus-specific resolution. Together, our results offer a novel perspective on the interplay between epigenetic repression, L1 evolution, and genome stability.

Mapping the splicing landscape of the human immune system

Most human genes code for more than one transcript. Different ratios of transcripts of the same gene can be found in different cell types or states, indicating differential use of transcription start sites or differential splicing. Such differential transcript use (DTUs) events provide an additional layer of regulation and protein diversity. With the exceptions of PTPRC and CIITA, there are very few reported cases of DTU events in the immune system. To rigorously map DTUs between different human immune cell types, we leveraged four publicly available RNA sequencing datasets. We identified 282 DTU events between five human healthy immune cell types that appear in at least two datasets. The patterns of the DTU events were mostly cell-type-specific or lineage-specific, in the context of the five cell types tested. DTUs correlated with the expression pattern of potential regulators, namely, splicing factors and transcription factors. Of the several immune related conditions studied, only sepsis affected the splicing of more than a few genes and only in innate immune cells. Taken together, we map the DTUs landscape in human peripheral blood immune cell types, and present hundreds of genes whose transcript use changes between cell types or upon activation.

Combining TSS-MPRA and sensitive TSS profile dissimilarity scoring to study the sequence determinants of transcription initiation

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

Sequence basis of transcription initiation in human genome

Transcription initiation is an essential process for ensuring proper function of any gene, however, a unified understanding of sequence patterns and rules that determine transcription initiation sites in human genome remains elusive. By explaining transcription initiation at basepair resolution from sequence with a deep learning-inspired explainable modeling approach, here we show that simple rules can explain the vast majority of human promoters. We identified key sequence patterns that contribute to human promoter function, each activating transcription with a distinct position-specific effect curve that likely reflects its mechanism of promoting transcription initiation. Most of these position-specific effects have not been previously characterized, and we verified them using experimental perturbations of transcription factors and sequences. We revealed the sequence basis of bidirectional transcription at promoters and links between promoter selectivity and gene expression variation across cell types. Additionally, by analyzing 241 mammalian genomes and mouse transcription initiation site data, we showed that the sequence determinants are conserved across mammalian species. Taken together, we provide a unified model of the sequence basis of transcription initiation at the basepair level that is broadly applicable across mammalian species, and shed new light on basic questions related to promoter sequence and function.

CapTrap-Seq: A platform-agnostic and quantitative approach for high-fidelity full-length RNA transcript sequencing

Long-read RNA sequencing is essential to produce accurate and exhaustive annotation of eukaryotic genomes. Despite advancements in throughput and accuracy, achieving reliable end-to-end identification of RNA transcripts remains a challenge for long-read sequencing methods. To address this limitation, we developed CapTrap-seq, a cDNA library preparation method, which combines the Cap-trapping strategy with oligo(dT) priming to detect 5'capped, full-length transcripts, together with the data processing pipeline LyRic. We benchmarked CapTrap-seq and other popular RNA-seq library preparation protocols in a number of human tissues using both ONT and PacBio sequencing. To assess the accuracy of the transcript models produced, we introduced a capping strategy for synthetic RNA spike-in sequences that mimics the natural 5'cap formation in RNA spike-in molecules. We found that the vast majority (up to 90%) of transcript models that LyRic derives from CapTrap-seq reads are full-length. This makes it possible to produce highly accurate annotations with minimal human intervention.

Investigating the NRAS 5' UTR as a target for small molecules

Neuroblastoma RAS (NRAS) is an oncogene that is deregulated and highly mutated in cancers including melanomas and acute myeloid leukemias. The 5' untranslated region (UTR) (5' UTR) of the NRAS mRNA contains a G-quadruplex (G4) that regulates translation. Here we report a novel class of small molecule that binds to the G4 structure located in the 5' UTR of the NRAS mRNA. We used a small molecule microarray screen to identify molecules that selectively bind to the NRAS-G4 with submicromolar affinity. One compound inhibits the translation of NRAS in vitro but showed only moderate effects on the NRAS levels in cellulo. Rapid Amplification of cDNA Ends and RT-PCR analysis revealed that the predominant NRAS transcript does not possess the G4 structure. Thus, although NRAS transcripts lack a G4 in many cell lines the concept of targeting folded regions within 5' UTRs to control translation remains a highly attractive strategy.

Viral proteogenomic and expression profiling during productive replication of a skin-tropic herpesvirus in the natural host

Efficient transmission of herpesviruses is essential for dissemination in host populations; however, little is known about the viral genes that mediate transmission, mostly due to a lack of natural virus-host model systems. Marek's disease is a devastating herpesviral disease of chickens caused by Marek's disease virus (MDV) and an excellent natural model to study skin-tropic herpesviruses and transmission. Like varicella zoster virus that causes chicken pox in humans, the only site where infectious cell-free MD virions are efficiently produced is in epithelial skin cells, a requirement for host-to-host transmission. Here, we enriched for heavily infected feather follicle epithelial skin cells of live chickens to measure both viral transcription and protein expression using combined short- and long-read RNA sequencing and LC/MS-MS bottom-up proteomics. Enrichment produced a previously unseen breadth and depth of viral peptide sequencing. We confirmed protein translation for 84 viral genes at high confidence (1% FDR) and correlated relative protein abundance with RNA expression levels. Using a proteogenomic approach, we confirmed translation of most well-characterized spliced viral transcripts and identified a novel, abundant isoform of the 14 kDa transcript family via IsoSeq transcripts, short-read intron-spanning sequencing reads, and a high-quality junction-spanning peptide identification. We identified peptides representing alternative start codon usage in several genes and putative novel microORFs at the 5' ends of two core herpesviral genes, pUL47 and ICP4, along with strong evidence of independent transcription and translation of the capsid scaffold protein pUL26.5. Using a natural animal host model system to examine viral gene expression provides a robust, efficient, and meaningful way of validating results gathered from cell culture systems.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single-cell resolution

Synthetic glucocorticoids, such as dexamethasone, have been used as a treatment for many immune conditions, such as asthma and, more recently, severe COVID-19. Single-cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single-cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated peripheral blood mononuclear cells from 96 African American children. We used novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we showed that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and showed widespread genetic regulation of the transcriptional dynamics of the gene expression response.

Transcription start site heterogeneity and its role in RNA fate determination distinguish HIV-1 from other retroviruses and are mediated by core promoter elements

HIV-1 uses heterogeneous transcription start sites (TSSs) to generate two RNA 5' isoforms that adopt radically different structures and perform distinct replication functions. Although these RNAs differ in length by only two bases, exclusively the shorter RNA is encapsidated while the longer RNA is excluded from virions and provides intracellular functions. The current study examined TSS usage and packaging selectivity for a broad range of retroviruses and found that heterogenous TSS usage was a conserved feature of all tested HIV-1 strains, but all other retroviruses examined displayed unique TSSs. Phylogenetic csomparisons and chimeric viruses' properties provided evidence that this mechanism of RNA fate determination was an innovation of the HIV-1 lineage, with determinants mapping to core promoter elements. Fine-tuning differences between HIV-1 and HIV-2, which uses a unique TSS, implicated purine residue positioning plus a specific TSS-adjacent dinucleotide in specifying multiplicity of TSS usage. Based on these findings, HIV-1 expression constructs were generated that differed from the parental strain by only two point mutations yet each expressed only one of HIV-1's two RNAs. Replication defects of the variant with only the presumptive founder TSS were less severe than those for the virus with only the secondary start site.

Mutational scans reveal differential evolvability of Drosophila promoters and enhancers

Rapid enhancer and slow promoter evolution have been demonstrated through comparative genomics. However, it is not clear how this information is encoded genetically and if this can be used to place evolution in a predictive context. Part of the challenge is that our understanding of the potential for regulatory evolution is biased primarily toward natural variation or limited experimental perturbations. Here, to explore the evolutionary capacity of promoter variation, we surveyed an unbiased mutation library for three promoters in Drosophila melanogaster. We found that mutations in promoters had limited to no effect on spatial patterns of gene expression. Compared to developmental enhancers, promoters are more robust to mutations and have more access to mutations that can increase gene expression, suggesting that their low activity might be a result of selection. Consistent with these observations, increasing the promoter activity at the endogenous locus of shavenbaby led to increased transcription yet limited phenotypic changes. Taken together, developmental promoters may encode robust transcriptional outputs allowing evolvability through the integration of diverse developmental enhancers.

This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

Decoding murine cytomegalovirus

The genomes of both human cytomegalovirus (HCMV) and murine cytomegalovirus (MCMV) were first sequenced over 20 years ago. Similar to HCMV, the MCMV genome had initially been proposed to harbor ≈170 open reading frames (ORFs). More recently, omics approaches revealed HCMV gene expression to be substantially more complex comprising several hundred viral ORFs. Here, we provide a state-of-the art reannotation of lytic MCMV gene expression based on integrative analysis of a large set of omics data. Our data reveal 365 viral transcription start sites (TiSS) that give rise to 380 and 454 viral transcripts and ORFs, respectively. The latter include >200 small ORFs, some of which represented the most highly expressed viral gene products. By combining TiSS profiling with metabolic RNA labelling and chemical nucleotide conversion sequencing (dSLAM-seq), we provide a detailed picture of the expression kinetics of viral transcription. This not only resulted in the identification of a novel MCMV immediate early transcript encoding the m166.5 ORF, which we termed ie4, but also revealed a group of well-expressed viral transcripts that are induced later than canonical true late genes and contain an initiator element (Inr) but no TATA- or TATT-box in their core promoters. We show that viral upstream ORFs (uORFs) tune gene expression of longer viral ORFs expressed in cis at translational level. Finally, we identify a truncated isoform of the viral NK-cell immune evasin m145 arising from a viral TiSS downstream of the canonical m145 mRNA. Despite being ≈5-fold more abundantly expressed than the canonical m145 protein it was not required for downregulating the NK cell ligand, MULT-I. In summary, our work will pave the way for future mechanistic studies on previously unknown cytomegalovirus gene products in an important virus animal model.

Alternative promoters in CpG depleted regions are prevalently associated with epigenetic misregulation of liver cancer transcriptomes

Transcriptional regulation is commonly governed by alternative promoters. However, the regulatory architecture in alternative and reference promoters, and how they differ, remains elusive. In 100 CAGE-seq libraries from hepatocellular carcinoma patients, here we annotate 4083 alternative promoters in 2926 multi-promoter genes, which are largely undetected in normal livers. These genes are enriched in oncogenic processes and predominantly show association with overall survival. Alternative promoters are narrow nucleosome depleted regions, CpG island depleted, and enriched for tissue-specific transcription factors. Globally tumors lose DNA methylation. We show hierarchical retention of intragenic DNA methylation with CG-poor regions rapidly losing methylation, while CG-rich regions retain it, a process mediated by differential SETD2, H3K36me3, DNMT3B, and TET1 binding. This mechanism is validated in SETD2 knockdown cells and SETD2-mutated patients. Selective DNA methylation loss in CG-poor regions makes the chromatin accessible for alternative transcription. We show alternative promoters can control tumor transcriptomes and their regulatory architecture.