Sites of transcription initiation drive mRNA isoform selection.

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

Pervasive, Coordinated Protein-Level Changes Driven by Transcript Isoform Switching during Meiosis.

To better understand the gene regulatory mechanisms that program developmental processes, we carried out simultaneous genome-wide measurements of mRNA, translation, and protein through meiotic differentiation in budding yeast. Surprisingly, we observed that the levels of several hundred mRNAs are anti-correlated with their corresponding protein products. We show that rather than arising from canonical forms of gene regulatory control, the regulation of at least 380 such cases, or over 8% of all measured genes, involves temporally regulated switching between production of a canonical, translatable transcript and a 5' extended isoform that is not efficiently translated into protein. By this pervasive mechanism for the modulation of protein levels through a natural developmental program, a single transcription factor can coordinately activate and repress protein synthesis for distinct sets of genes. The distinction is not based on whether or not an mRNA is induced but rather on the type of transcript produced.

Efficient, specific, and combinatorial control of endogenous exon splicing with dCasRx-RBM25.

Efficient targeted control of splicing is a major goal of functional genomics and therapeutic applications. Guide (g)RNA-directed, deactivated (d)Cas CRISPR enzymes fused to splicing effectors represent a promising strategy due to the flexibility of these systems. However, efficient, specific, and generalizable activation of endogenous exons using this approach has not been previously reported. By screening over 300 dCasRx-splicing factor fusion proteins tethered to splicing reporters, we identify dCasRx-RBM25 as a potent activator of exons. Moreover, dCasRx-RBM25 efficiently activates the splicing of ∼90% of targeted endogenous alternative exons and displays high on-target specificity. Using gRNA arrays for combinatorial targeting, we demonstrate that dCasRx-RBM25 enables multiplexed activation and repression of exons. Using this feature, the targeting of neural-regulated exons in Ptpb1 and Puf60 in embryonic stem cells reveals combinatorial effects on downstream alternative splicing events controlled by these factors. Collectively, our results enable versatile, combinatorial exon-resolution functional assays and splicing-directed therapeutic applications.

Integrated genomic analysis reveals key features of long undecoded transcript isoform-based gene repression.

Long undecoded transcript isoforms (LUTIs) represent a class of non-canonical mRNAs that downregulate gene expression through the combined act of transcriptional and translational repression. While single gene studies revealed important aspects of LUTI-based repression, how these features affect gene regulation on a global scale is unknown. Using transcript leader and direct RNA sequencing, here, we identify 74 LUTI candidates that are specifically induced in meiotic prophase. Translational repression of these candidates appears to be ubiquitous and is dependent on upstream open reading frames. However, LUTI-based transcriptional repression is variable. In only 50% of the cases, LUTI transcription causes downregulation of the protein-coding transcript isoform. Higher LUTI expression, enrichment of histone 3 lysine 36 trimethylation, and changes in nucleosome position are the strongest predictors of LUTI-based transcriptional repression. We conclude that LUTIs downregulate gene expression in a manner that integrates translational repression, chromatin state changes, and the magnitude of LUTI expression.

Distinct Structures and Dynamics of Chromatosomes with Different Human Linker Histone Isoforms.

The repeating structural unit of metazoan chromatin is the chromatosome, a nucleosome bound to a linker histone, H1. There are 11 human H1 isoforms with diverse cellular functions, but how they interact with the nucleosome remains elusive. Here, we determined the cryoelectron microscopy (cryo-EM) structures of chromatosomes containing 197 bp DNA and three different human H1 isoforms, respectively. The globular domains of all three H1 isoforms bound to the nucleosome dyad. However, the flanking/linker DNAs displayed substantial distinct dynamic conformations. Nuclear magnetic resonance (NMR) and H1 tail-swapping cryo-EM experiments revealed that the C-terminal tails of the H1 isoforms mainly controlled the flanking DNA orientations. We also observed partial ordering of the core histone H2A C-terminal and H3 N-terminal tails in the chromatosomes. Our results provide insights into the structures and dynamics of the chromatosomes and have implications for the structure and function of chromatin.

The clinical utility and diagnostic implementation of human subject cell transdifferentiation followed by RNA sequencing.

RNA sequencing (RNA-seq) has recently been used in translational research settings to facilitate diagnoses of Mendelian disorders. A significant obstacle for clinical laboratories in adopting RNA-seq is the low or absent expression of a significant number of disease-associated genes/transcripts in clinically accessible samples. As this is especially problematic in neurological diseases, we developed a clinical diagnostic approach that enhanced the detection and evaluation of tissue-specific genes/transcripts through fibroblast-to-neuron cell transdifferentiation. The approach is designed specifically to suit clinical implementation, emphasizing simplicity, cost effectiveness, turnaround time, and reproducibility. For clinical validation, we generated induced neurons (iNeurons) from 71 individuals with primary neurological phenotypes recruited to the Undiagnosed Diseases Network. The overall diagnostic yield was 25.4%. Over a quarter of the diagnostic findings benefited from transdifferentiation and could not be achieved by fibroblast RNA-seq alone. This iNeuron transcriptomic approach can be effectively integrated into diagnostic whole-transcriptome evaluation of individuals with genetic disorders.

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease.

A major fraction of loci identified by genome-wide association studies (GWASs) mediate alternative splicing, but mechanistic interpretation is hindered by the technical limitations of short-read RNA sequencing (RNA-seq), which cannot directly link splicing events to full-length protein isoforms. Long-read RNA-seq represents a powerful tool to characterize transcript isoforms, and recently, infer protein isoform existence. Here, we present an approach that integrates information from GWASs, splicing quantitative trait loci (sQTLs), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes that colocalized with BMD associations (H4PP ≥ 0.75). We generated PacBio Iso-Seq data (N = ∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were unannotated. By casting the sQTLs onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense-mediated decay and 190 that potentially resulted in the expression of unannotated protein isoforms. Finally, we functionally validated colocalizing sQTLs in TPM2, in which siRNA-mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization but exhibited no effect upon knockdown of the entire gene. Our approach should be to generalize across diverse clinical traits and to provide insights into protein isoform activities modulated by GWAS loci.

Oncofetal IGF2BP3-mediated control of microRNA structural diversity in the malignancy of early-stage lung adenocarcinoma.

The nature of microRNA (miRNA) dysfunction in carcinogenesis remains controversial because of the complex connection between miRNA structural diversity and biological processes. Here, we found that oncofetal IGF2BP3 regulates the selective production of a subset of 3'-isoforms (3'-isomiRs), including miR-21-5p and Let-7 family, which induces significant changes in their cellular seed occupancy and structural components, establishing a cancer-specific gene expression profile. The D-score, reflecting dominant production of a representative miR-21-5p+C (a 3'-isomiR), discriminated between clinical early-stage lung adenocarcinoma (LUAD) cases with low and high recurrence risks, and was associated with molecular features of cell cycle progression, epithelial-mesenchymal transition pressure, and immune evasion. We found that IGF2BP3 controls the production of miR-21-5p+C by directing the nuclear Drosha complex to select the cleavage site. IGF2BP3 was also involved in the production of 3'-isomiRs of miR-425-5p and miR-454-3p. IGF2BP3-regulated these three miRNAs are suggested to be associated with the regulation of p53, TGF-ß, and TNF pathways in LUAD. Knockdown of IGF2BP3 also induced a selective upregulation of Let-7 3'-isomiRs, leading to increased cellular Let-7 seed occupancy and broad repression of its target genes encoding cell cycle regulators. The D-score is an index that reflects this cellular situation. Our results suggest that the aberrant regulation of miRNA structural diversity is a critical component for controlling cellular networks, thus supporting the establishment of a malignant gene expression profile in early stage LUAD.

Quantitative Control for Stoichiometric Protein Synthesis.

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production.

Extensive Structural Differences of Closely Related 3' mRNA Isoforms: Links to Pab1 Binding and mRNA Stability.

Alternative polyadenylation generates numerous 3' mRNA isoforms that can vary in biological properties, such as stability and localization. We developed methods to obtain transcriptome-scale structural information and protein binding on individual 3' mRNA isoforms in vivo. Strikingly, near-identical mRNA isoforms can possess dramatically different structures throughout the 3' UTR. Analyses of identical mRNAs in different species or refolded in vitro indicate that structural differences in vivo are often due to trans-acting factors. The level of Pab1 binding to poly(A)-containing isoforms is surprisingly variable, and differences in Pab1 binding correlate with the extent of structural variation for closely spaced isoforms. A pattern encompassing single-strandedness near the 3' terminus, double-strandedness of the poly(A) tail, and low Pab1 binding is associated with mRNA stability. Thus, individual 3' mRNA isoforms can be remarkably different physical entities in vivo. Sequences responsible for isoform-specific structures, differential Pab1 binding, and mRNA stability are evolutionarily conserved, indicating biological function.

A vertebrate Vangl2 translational variant required for planar cell polarity.

First described in the milkweed bug Oncopeltus fasciatus, planar cell polarity (PCP) is a developmental process essential for embryogenesis and development of polarized structures in Metazoans. This signaling pathway involves a set of evolutionarily conserved genes encoding transmembrane (Vangl, Frizzled, Celsr) and cytoplasmic (Prickle, Dishevelled) molecules. Vangl2 is of major importance in embryonic development as illustrated by its pivotal role during neural tube closure in human, mouse, Xenopus, and zebrafish embryos. Here, we report on the molecular and functional characterization of a Vangl2 isoform, Vangl2-Long, containing an N-terminal extension of about 50 aa, which arises from an alternative near-cognate AUA translation initiation site, lying upstream of the conventional start codon. While missing in Vangl1 paralogs and in all invertebrates, including Drosophila, this N-terminal extension is conserved in all vertebrate Vangl2 sequences. We show that Vangl2-Long belongs to a multimeric complex with Vangl1 and Vangl2. Using morpholino oligonucleotides to specifically knockdown Vangl2-Long in Xenopus, we found that this isoform is functional and required for embryo extension and neural tube closure. Furthermore, both Vangl2 and Vangl2-Long must be correctly expressed for the polarized distribution of the PCP molecules Pk2 and Dvl1 and for centriole rotational polarity in ciliated epidermal cells. Altogether, our study suggests that Vangl2-Long significantly contributes to the pool of Vangl2 molecules present at the plasma membrane to maintain PCP in vertebrate tissues.

Myosin-1C differentially displaces tropomyosin isoforms altering their inhibition of motility.

Force generation and motility by actomyosin in nonmuscle cells are spatially regulated by ∼40 tropomyosin (Tpm) isoforms. The means by which Tpms are targeted to specific cellular regions and the mechanisms that result in differential activity of myosin paralogs are unknown. We show that Tpm3.1 and Tpm1.7 inhibit Myosin-IC (Myo1C), with Tpm1.7 more effectively reducing the number of gliding filaments than Tpm3.1. Strikingly, cosedimentation and fluorescence microscopy assays revealed that Tpm3.1 is displaced from actin by Myo1C and not by myosin-II. In contrast, Tpm1.7 is only weakly displaced by Myo1C. Unlike other characterized myosins, Myo1C motility is inhibited by Tpm when the Tpm-actin filament is activated by myosin-II. These results point to a mechanism for the exclusion of myosin-I paralogs from cellular Tpm-decorated actin filaments that are activated by other myosins. Additionally, our results suggest a potential mechanism for myosin-induced Tpm sorting in cells.

The structures of salt-inducible kinase 3 in complex with inhibitors reveal determinants for binding and selectivity.

The salt-inducible kinases (SIKs) 1 to 3, belonging to the AMPK-related kinase family, serve as master regulators orchestrating a diverse set of physiological processes such as metabolism, bone formation, immune response, oncogenesis, and cardiac rhythm. Owing to its key regulatory role, the SIK kinases have emerged as compelling targets for pharmacological intervention across a diverse set of indications. Therefore, there is interest in developing SIK inhibitors with defined selectivity profiles both to further dissect the downstream biology and for treating disease. However, despite a large pharmaceutical interest in the SIKs, experimental structures of SIK kinases are scarce. This is likely due to the challenges associated with the generation of proteins suitable for structural studies. By adopting a rational approach to construct design and protein purification, we successfully crystallized and subsequently solved the structure of SIK3 in complex with HG-9-91-01, a potent SIK inhibitor. To enable further SIK3-inhibitor complex structures we identified an antibody fragment that facilitated crystallization and enabled a robust protocol suitable for structure-based drug design. The structures reveal SIK3 in an active conformation, where the ubiquitin-associated domain is shown to provide further stabilization to this active conformation. We present four pharmacologically relevant and distinct SIK3-inhibitor complexes. These detail the key interaction for each ligand and reveal how different regions of the ATP site are engaged by the different inhibitors to achieve high affinity. Notably, the structure of SIK3 in complex with a SIK3 specific inhibitor offers insights into isoform selectivity.

Nanoblot: an R-package for visualization of RNA isoforms from long-read RNA-sequencing data.

RT-PCR and northern blots have long been used to study RNA isoforms usage for single genes. Recent advancements in long-read sequencing have yielded unprecedented information about the usage and abundance of these RNA isoforms. However, visualization of long-read sequencing data remains challenging due to the high information density. To alleviate these issues, we have developed NanoBlot, an open-source R-package that generates northern blot and RT-PCR-like images from long-read sequencing data. NanoBlot requires aligned, positionally sorted and indexed BAM files. Plotting is based around ggplot2 and is easily customizable. Advantages of NanoBlot include a robust system for designing probes to visualize isoforms including excluding reads based on the presence or absence of a specified region, an elegant solution to representing isoforms with continuous variations in length, and the ability to overlay multiple genes in the same plot using different colors. We present examples of nanoblots compared to actual northern blot data. In addition to traditional gel-like images, the NanoBlot package can also output other visualizations such as violin plots and 3'-RACE-like plots focused on 3'-end isoforms visualization. The use of the NanoBlot package should provide a simple answer to some of the challenges of visualizing long-read RNA-sequencing data.

Network integration and protein structural binding analysis of neurodegeneration-related interactome.

Neurodegenerative diseases (NDs) usually connect with aggregation and molecular interactions of pathological proteins. The integration of accumulative data from clinical and biomedical research will allow for the excavation of pathological proteins and related interactors. It is also important to systematically study their interacting proteins in order to find more related proteins and potential therapeutic targets. Understanding binding regions in protein interactions will help functional proteomics and provide an alternative method for predicting novel interactions. This study integrated data from biomedical research to achieve systematic mining and analysis of pathogenic proteins and their interaction network. A workflow has been built as a solution for the collective information of proteins involved in NDs, related protein-protein interactions (PPIs) and interactive visualizations. It also included protein isoforms and mapped them in a disease-related PPI network to illuminate the impact of alternative splicing on protein binding. The interacting proteins enriched by diseases and biological processes (BPs) revealed possible regulatory modules. A high-resolution network with structural affinity information was generated. Finally, Neurodegenerative Disease Atlas (NDAtlas) was constructed with an interactive and intuitive view of protein docking with 3D molecular graphics beyond the traditional 2D network. NDAtlas is available at http://bis.zju.edu.cn/ndatlas.

MiR-30a-5p isoform -1|1 promotes the progression of gastric cancer by inhibiting TMEM66 and reducing intratumoral cytotoxic T cells.

Gastric cancer is histologically classified into the intestinal subtype, which forms tubular structures, and the aggressive diffuse subtype, characterized by rapid invasion and poor prognosis. The variety and quantity of miRNA isoforms between different histological subtypes of gastric cancer were unknown. Through systematic filtering, we found that more diverse miR-30a-5p isoforms was present in the diffuse subtype of gastric cancer, and was associated with patients' worse survival independent of tumor stage based on the TCGA miRNA-seq data. Among all nine isoforms of miR-30a-5p, miR-30a-5p -1|1 was more abundant than the archetype of miR-30a-5p. Higher expression of miR-30a-5p -1|1 was observed in patients with advanced tumor stage and poor survival. Furthermore, miR-30a-5p -1|1 could promote the metastasis of gastric cancer cells both in vitro and in vivo by down-regulating TMEM66. In clinical samples, decreased expression of TMEM66 was characteristic of gastric cancer, and the low level of TMEM66 correlated with deceased CD8 positive cells in the tumor microenvironment probably due to decreased cytokines production. In conclusion, the variety of miR-30a-5p isoforms correlates with worse survival in gastric cancer patients. Moreover, miR-30a-5p -1|1 could promote gastric cancer metastasis by inhibiting TMEM66 and the infiltration of intratumoral CD8 positive cells.

AhR promotes phosphorylation of ARNT isoform 1 in human T cell malignancies as a switch for optimal AhR activity.

The aryl hydrocarbon receptor nuclear translocator (ARNT) is a transcription factor present in immune cells as a long and short isoform, referred to as isoforms 1 and 3, respectively. However, investigation into potential ARNT isoform­specific immune functions is lacking despite the well-established heterodimerization requirement of ARNT with, and for the activity of, the aryl hydrocarbon receptor (AhR), a critical mediator of immune homeostasis. Here, using global and targeted transcriptomics analyses, we show that the relative ARNT isoform 1:3 ratio in human T cell lymphoma cells dictates the amplitude and direction of AhR target gene regulation. Specifically, shifting the ARNT isoform 1:3 ratio lower by suppressing isoform 1 enhances, or higher by suppressing isoform 3 abrogates, AhR responsiveness to ligand activation through preprograming a cellular genetic background that directs explicit gene expression patterns. Moreover, the fluctuations in gene expression patterns that accompany a decrease or increase in the ARNT isoform 1:3 ratio are associated with inflammation or immunosuppression, respectively. Molecular studies identified the unique casein kinase 2 (CK2) phosphorylation site within isoform 1 as an essential parameter to the mechanism of ARNT isoform­specific regulation of AhR signaling. Notably, CK2-mediated phosphorylation of ARNT isoform 1 is dependent on ligand-induced AhR nuclear translocation and is required for optimal AhR target gene regulation. These observations reveal ARNT as a central modulator of AhR activity predicated on the status of the ARNT isoform ratio and suggest that ARNT-based therapies are a viable option for tuning the immune system to target immune disorders.

Structure and Roles of V-type ATPases.

V-ATPases are membrane-embedded protein complexes that function as ATP hydrolysis-driven proton pumps. V-ATPases are the primary source of organellar acidification in all eukaryotes, making them essential for many fundamental cellular processes. Enzymatic activity can be modulated by regulated and reversible disassembly of the complex, and several subunits of mammalian V-ATPase have multiple isoforms that are differentially localized. Although the biochemical properties of the different isoforms are currently unknown, mutations in specific subunit isoforms have been associated with various diseases, making V-ATPases potential drug targets. V-ATPase structure and activity have been best characterized in Saccharomyces cerevisiae, where recent structures have revealed details about the dynamics of the enzyme, the proton translocation pathway, and conformational changes associated with regulated disassembly and autoinhibition.

ClusTrast: a short read de novo transcript isoform assembler guided by clustered contigs.

BACKGROUND: Transcriptome assembly from RNA-sequencing data in species without a reliable reference genome has to be performed de novo, but studies have shown that de novo methods often have inadequate ability to reconstruct transcript isoforms. We address this issue by constructing an assembly pipeline whose main purpose is to produce a comprehensive set of transcript isoforms. RESULTS: We present the de novo transcript isoform assembler ClusTrast, which takes short read RNA-seq data as input, assembles a primary assembly, clusters a set of guiding contigs, aligns the short reads to the guiding contigs, assembles each clustered set of short reads individually, and merges the primary and clusterwise assemblies into the final assembly. We tested ClusTrast on real datasets from six eukaryotic species, and showed that ClusTrast reconstructed more expressed known isoforms than any of the other tested de novo assemblers, at a moderate reduction in precision. For recall, ClusTrast was on top in the lower end of expression levels (<15% percentile) for all tested datasets, and over the entire range for almost all datasets. Reference transcripts were often (35-69% for the six datasets) reconstructed to at least 95% of their length by ClusTrast, and more than half of reference transcripts (58-81%) were reconstructed with contigs that exhibited polymorphism, measuring on a subset of reliably predicted contigs. ClusTrast recall increased when using a union of assembled transcripts from more than one assembly tool as primary assembly. CONCLUSION: We suggest that ClusTrast can be a useful tool for studying isoforms in species without a reliable reference genome, in particular when the goal is to produce a comprehensive transcriptome set with polymorphic variants.

Generation of B7-H3 isoform regulated by ANXA2/NSUN2/YBX1 axis in human glioma.

In recent years, in the development of emerging immunotherapy, B7-H3 is also termed as CD276 and has become a novel chimeric antigen receptor (CAR)-T target against glioma and other tumours, and aroused extensive attention. However, B7-H3 has three isoforms (2, 3 and 4Ig) with the controversial expression and elusive function in tumour especially glioma. The current study mainly focuses on the regulatory factors and related mechanisms of generation of different B7-H3 isoforms. First, we have determined that 2Ig is dominant in glioma with high malignancy, and 4Ig is widely expressed, whereas 3Ig shows negative expression in all glioma. Next, we have further found that RNA binding protein annexin A2 (ANXA2) is essential for B7-H3 isoform maintenance, but fail to determine the choice of 4Ig or 2Ig. RNA methyltransferase NOP2/Sun RNA methyltransferase 2 (NSUN2) and 5-methylcytosine reader Y-box binding protein 1 (YBX1) facilitate the production of 2Ig. Our findings have uncovered a series of factors (ANXA2/NSUN2/YBX1) that can determine the alternative generation of different isoforms of B7-H3 in glioma. Our result aims to help peers gain a clearer understanding of the expression and regulatory mechanisms of B7H3 in tumour patients, and to provide better strategies for designing B7H3 as a target in immunotherapy.