Alternative Transcription Start Sites Research Articles

A Bayesian Change Point Model for Dynamic Alternative Transcription Start Site Usage During Cellular Differentiation.

ABSTRACT An alternative transcription start site (ATSS) is a major driving force for increasing the complexity of transcripts in human tissues. As a transcriptional regulatory mechanism, ATSS has biological significance. Many studies have confirmed that ATSS plays an important role in diseases and cell development and differentiation. However, exploration of its dynamic mechanisms remains insufficient. Identifying ATSS change points during cell differentiation is critical for elucidating potential dynamic mechanisms. For relative ATSS usage as percentage data, the existing methods lack sensitivity to detect the change point for ATSS longitudinal data. In addition, some methods have strict requirements for data distribution and cannot be applied to deal with this problem. In this study, the Bayesian change point detection model was first constructed using reparameterization techniques for two parameters of a beta distribution for the percentage data type, and the posterior distributions of parameters and change points were obtained using Markov Chain Monte Carlo (MCMC) sampling. With comprehensive simulation studies, the performance of the Bayesian change point detection model is found to be consistently powerful and robust across most scenarios with different sample sizes and beta distributions. Second, differential ATSS events in the real data, whose change points were identified using our method, were clustered according to their change points. Last, for each change point, pathway and transcription factor motif analyses were performed on its differential ATSS events. The results of our analyses demonstrated the effectiveness of the Bayesian change point detection model and provided biological insights into cell differentiation.

Characterization and expression of the microtubule associated protein tau gene isoforms and their impact on clinical outcomes in glioma patients

Gliomas are the most common tumors found in the central nervous system. They are classified based on their histopathological malignancy. Low-grade gliomas (LGG) can turn into glioblastomas (GBM), the most aggressive form of brain cancer with the worst prognosis. Previous studies have shown that microtubule-associated protein tau (MAPT) is increased in GBM and that reducing it can improve survival in animal models. On the other hand, elevated MAPT expression has been linked to increased survival in LGG, indicating its potential as a survival biomarker. However, previous studies have primarily focused on the canonical protein-coding transcripts and gene-level expression of MAPT, disregarding the diversity of its many isoforms. Many diseases in almost every tissue have been associated with a dysregulation of alternative splicing, either through an altered ratio of the expressed isoforms or by expressing aberrant transcripts with different alternative transcription start and termination sites. This study investigates the association of specific MAPT transcripts with survival in LGG and GBM using publicly available genomic datasets and performing statistical analyses. The results of the study indicated coding and several non-coding isoforms associated with overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI). After including the relevant predictors as covariates, it was observed that among the protein-coding isoforms, MAPT-205 and MAPT-201 were significantly associated with clinical survival endpoints, but with opposite trends. These findings highlight the importance of studying the role of MAPT isoforms in glioma cellular phenotypes, aggressiveness, and patient survival.

Differences in 5'untranslated regions highlight the importance of translational regulation of dosage sensitive genes

BackgroundUntranslated regions (UTRs) are important mediators of post-transcriptional regulation. The length of UTRs and the composition of regulatory elements within them are known to vary substantially across genes, but little is known about the reasons for this variation in humans. Here, we set out to determine whether this variation, specifically in 5’UTRs, correlates with gene dosage sensitivity.ResultsWe investigate 5’UTR length, the number of alternative transcription start sites, the potential for alternative splicing, the number and type of upstream open reading frames (uORFs) and the propensity of 5’UTRs to form secondary structures. We explore how these elements vary by gene tolerance to loss-of-function (LoF; using the LOEUF metric), and in genes where changes in dosage are known to cause disease. We show that LOEUF correlates with 5’UTR length and complexity. Genes that are most intolerant to LoF have longer 5’UTRs, greater TSS diversity, and more upstream regulatory elements than their LoF tolerant counterparts. We show that these differences are evident in disease gene-sets, but not in recessive developmental disorder genes where LoF of a single allele is tolerated.ConclusionsOur results confirm the importance of post-transcriptional regulation through 5'UTRs in tight regulation of mRNA and protein levels, particularly for genes where changes in dosage are deleterious and lead to disease. Finally, to support gene-based investigation we release a web-based browser tool, VuTR, that supports exploration of the composition of individual 5'UTRs and the impact of genetic variation within them.

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.

RNAirport: a deep neural network–based database characterizing representative gene models in plants

A 5′-leader, known initially as the 5′-untranslated region, contains multiple isoforms due to alternative splicing (aS) and alternative transcription start site (aTSS). Therefore, a representative 5′-leader is demanded to examine the embedded RNA regulatory elements in controlling translation efficiency. Here, we develop a ranking algorithm and a deep-learning model to annotate representative 5′-leaders for five plant species. We rank the intra-sample and inter-sample frequency of aS-mediated transcript isoforms using the Kruskal–Wallis test–based algorithm and identify the representative aS-5′-leader. To further assign a representative 5′-end, we train the deep-learning model 5′leaderP to learn aTSS-mediated 5′-end distribution patterns from cap-analysis gene expression data. The model accurately predicts the 5′-end, confirmed experimentally in Arabidopsis and rice. The representative 5′-leader-contained gene models and 5′leaderP can be accessed at RNAirport (http://www.rnairport.com/leader5P/). This stage 1 5′-leader annotation records 5′-leader diversity and will pave the way to Ribo-Seq open-reading frame annotation, identical to the project recently initiated by human GENCODE.

BORIS/CTCFL epigenetically reprograms clustered CTCF binding sites into alternative transcriptional start sites

BackgroundPervasive usage of alternative promoters leads to the deregulation of gene expression in carcinogenesis and may drive the emergence of new genes in spermatogenesis. However, little is known regarding the mechanisms underpinning the activation of alternative promoters.ResultsHere we describe how alternative cancer-testis-specific transcription is activated. We show that intergenic and intronic CTCF binding sites, which are transcriptionally inert in normal somatic cells, could be epigenetically reprogrammed into active de novo promoters in germ and cancer cells. BORIS/CTCFL, the testis-specific paralog of the ubiquitously expressed CTCF, triggers the epigenetic reprogramming of CTCF sites into units of active transcription. BORIS binding initiates the recruitment of the chromatin remodeling factor, SRCAP, followed by the replacement of H2A histone with H2A.Z, resulting in a more relaxed chromatin state in the nucleosomes flanking the CTCF binding sites. The relaxation of chromatin around CTCF binding sites facilitates the recruitment of multiple additional transcription factors, thereby activating transcription from a given binding site. We demonstrate that the epigenetically reprogrammed CTCF binding sites can drive the expression of cancer-testis genes, long noncoding RNAs, retro-pseudogenes, and dormant transposable elements.ConclusionsThus, BORIS functions as a transcription factor that epigenetically reprograms clustered CTCF binding sites into transcriptional start sites, promoting transcription from alternative promoters in both germ cells and cancer cells.

Synaptopodin-2 Isoforms Have Specific Binding Partners and Display Distinct, Muscle Cell Type-Specific Expression Patterns.

Synaptopodin-2 (SYNPO2) is a protein associated with the Z-disc in striated muscle cells. It interacts with α-actinin and filamin C, playing a role in Z-disc maintenance under stress by chaperone-assisted selective autophagy (CASA). In smooth muscle cells, SYNPO2 is a component of dense bodies. Furthermore, it has been proposed to play a role in tumor cell proliferation and metastasis in many different kinds of cancers. Alternative transcription start sites and alternative splicing predict the expression of six putative SYNPO2 isoforms differing by extended amino- and/or carboxy-termini. Our analyses at mRNA and protein levels revealed differential expression of SYNPO2 isoforms in cardiac, skeletal and smooth muscle cells. We identified synemin, an intermediate filament protein, as a novel binding partner of the PDZ-domain in the amino-terminal extension of the isoforms mainly expressed in cardiac and smooth muscle cells, and demonstrated colocalization of SYNPO2 and synemin in both cell types. A carboxy-terminal extension, mainly expressed in smooth muscle cells, is sufficient for association with dense bodies and interacts with α-actinin. SYNPO2 therefore represents an additional and novel link between intermediate filaments and the Z-discs in cardiomyocytes and dense bodies in smooth muscle cells, respectively. In pathological skeletal muscle samples, we identified SYNPO2 in the central and intermediate zones of target fibers of patients with neurogenic muscular atrophy, and in nemaline bodies. Our findings help to understand distinct functions of individual SYNPO2 isoforms in different muscle tissues, but also in tumor pathology.

CamoTSS: analysis of alternative transcription start sites for cellular phenotypes and regulatory patterns from 5' scRNA-seq data

Five-prime single-cell RNA-seq (scRNA-seq) has been widely employed to profile cellular transcriptomes, however, its power of analysing transcription start sites (TSS) has not been fully utilised. Here, we present a computational method suite, CamoTSS, to precisely identify TSS and quantify its expression by leveraging the cDNA on read 1, which enables effective detection of alternative TSS usage. With various experimental data sets, we have demonstrated that CamoTSS can accurately identify TSS and the detected alternative TSS usages showed strong specificity in different biological processes, including cell types across human organs, the development of human thymus, and cancer conditions. As evidenced in nasopharyngeal cancer, alternative TSS usage can also reveal regulatory patterns including systematic TSS dysregulations.

Transcriptional diversity of the oxytocin receptor in prairie voles: mechanistic implications for behavioral neuroscience and maternal physiology.

The neurohormone oxytocin regulates many aspects of physiology primarily by binding to its receptor, the oxytocin receptor. The oxytocin receptor gene (Oxtr) has been shown to have alternative transcripts in the mouse brain which may each have different biological functions or be used in specific contexts. A popular animal model for studying oxytocin-dependent social behaviors is the prairie vole, a biparental and monogamous rodent. Alternative transcriptional capacity of Oxtr in prairie voles is unknown. We used 5' rapid amplification of cDNA ends to identify alternative Oxtr transcription start sites in prairie vole brain tissue and uterine tissue. We then validated expression of specific transcripts in fetal brains and assessed the impact of exogenous oxytocin administration in utero on offspring brain development. We identified seven distinct Oxtr transcripts, all of which are present in both brain and uterine tissue. We then demonstrated that maternal oxytocin administration alters expression of a specific subset of Oxtr transcripts and that these different transcripts are under unique epigenetic regulation, such that in the perinatal period only one of the alternative transcripts is associated with DNA methylation in the Oxtr promoter. These data establish the existence of multiple Oxtr transcripts in prairie vole brain and uterine tissue and implicate oxytocin in the regulation of alternative transcript expression. These data have significant implications for our understanding of null mutant models in both mice and voles and translation in human birth and behavior.

Multi-landmark alignment of genomic signals reveals conserved expression patterns across transcription start sites

The prevalent one-dimensional alignment of genomic signals to a reference landmark is a cornerstone of current methods to study transcription and its DNA-dependent processes but it is prone to mask potential relations among multiple DNA elements. We developed a systematic approach to align genomic signals to multiple locations simultaneously by expanding the dimensionality of the genomic-coordinate space. We analyzed transcription in human and uncovered a complex dependence on the relative position of neighboring transcription start sites (TSSs) that is consistently conserved among cell types. The dependence ranges from enhancement to suppression of transcription depending on the relative distances to the TSSs, their intragenic position, and the transcriptional activity of the gene. Our results reveal a conserved hierarchy of alternative TSS usage within a previously unrecognized level of genomic organization and provide a general methodology to analyze complex functional relationships among multiple types of DNA elements.

The promoter as a trip navigator: Guiding alternative polyadenylation site destinations

Alfonso-Gonzalez etal.1 present an innovative combination of long-read-sequencing approaches that reveal coupling of alternative transcription start sites and alternative polyadenylation site usage on a global level.

Long-read direct RNA sequencing reveals epigenetic regulation of chimeric gene-transposon transcripts in Arabidopsis thaliana

Transposable elements (TEs) are accumulated in both intergenic and intragenic regions in plant genomes. Intragenic TEs often act as regulatory elements of associated genes and are also co-transcribed with genes, generating chimeric TE-gene transcripts. Despite the potential impact on mRNA regulation and gene function, the prevalence and transcriptional regulation of TE-gene transcripts are poorly understood. By long-read direct RNA sequencing and a dedicated bioinformatics pipeline, ParasiTE, we investigated the transcription and RNA processing of TE-gene transcripts in Arabidopsis thaliana. We identified a global production of TE-gene transcripts in thousands of A. thaliana gene loci, with TE sequences often being associated with alternative transcription start sites or transcription termination sites. The epigenetic state of intragenic TEs affects RNAPII elongation and usage of alternative poly(A) signals within TE sequences, regulating alternative TE-gene isoform production. Co-transcription and inclusion of TE-derived sequences into gene transcripts impact regulation of RNA stability and environmental responses of some loci. Our study provides insights into TE-gene interactions that contributes to mRNA regulation, transcriptome diversity, and environmental responses in plants.

Unraveling the influences of sequence and position on yeast uORF activity using massively parallel reporter systems and machine learning.

Upstream open-reading frames (uORFs) are potent cis-acting regulators of mRNA translation and nonsense-mediated decay (NMD). While both AUG- and non-AUG initiated uORFs are ubiquitous in ribosome profiling studies, few uORFs have been experimentally tested. Consequently, the relative influences of sequence, structural, and positional features on uORF activity have not been determined. We quantified thousands of yeast uORFs using massively parallel reporter assays in wildtype and ∆upf1 yeast. While nearly all AUG uORFs were robust repressors, most non-AUG uORFs had relatively weak impacts on expression. Machine learning regression modeling revealed that both uORF sequences and locations within transcript leaders predict their effect on gene expression. Indeed, alternative transcription start sites highly influenced uORF activity. These results define the scope of natural uORF activity, identify features associated with translational repression and NMD, and suggest that the locations of uORFs in transcript leaders are nearly as predictive as uORF sequences.

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

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

PO-01-145 DYNAMIC UTR COMPOSITION IMPACTS RNA BINDING PROTEIN COMPLEXING TO CONTROL TRANSLATION AND LIMIT CARDIAC GAP JUNCTION FORMATION DURING STRESS

Translational control has emerged as a critical yet dynamic process, representing one of the most rapid means by which the cardiomyocyte can alter its proteome. Several key cardiac mRNAs, including connexin43 (Cx43, gene name GJA1) are understood to undergo alternative translation initiation with resulting products (ie, GJA1-20k) impacting cardiac function. We previously found alterations to GJA1 mRNA untranslated region (UTR) length occur in aged and stressed heart tissue, suppressing internal translation of GJA1-20k, and limiting Cx43 gap junction formation. Having therefore identified dynamic mRNA cis element-based regulation of Cx43/GJA1-20k translation, we next asked what trans-acting factors are at play to alter GJA1 mRNA translation during stress. The IMP family of RNA binding proteins (RBPs) regulate alternative and canonical translation initiation of several genes, which led us to investigate how they may impact GJA1 mRNA. Using GJA1 as a model RNA, we test the hypothesis that alternative transcription start site usage controls inclusion of cis-acting elements in the mRNA 5’UTR affecting RBP binding and regulation of translation. Pathological hypertrophy was induced by isoproterenol deliverey to adult mice via implanted minipumps. Mouse models are complemented with matured human IPSC-derived cardiomyocytes, primary mouse ventricular cardiomyocytes, exposed to chronic hypoxia or TGF-β in vitro. Single molecule fluorescence in situ hybridization (smFISH) coupled with super-resolution microscopy localize GJA1 mRNA and RBPs. Hypertrophic hearts show reduced GJA1-20k and IMP1 expression. Using smFISH to label GJA1 mRNA, we find IMP1 complexes with GJA1 mRNA predominantly in comparison to IMP2 and 3, and that this complexing is reduced during stress, concomitant with loss of GJA1-20k internal translation. Biochemical assays confirm GJA1 mRNA/IMP1 binding by RNA-IP and synthetic reduction of IMP1 expression by RNAi or CRISPR/Cas9 alters GJA1 mRNA localization, reduces internal translation of GJA1-20k, and limits gap junction formation. We find comparable results during hypoxia in primary mouse cardiomyocytes and human induced pluripotent stem cell derived cardiomyocytes. Dynamic GJA1-5’UTR composition therefore represents a powerful mechanism connecting alterations in transcription to gene-specific translational regulation by RBPs such as IMP1, disruption of intercellular communication during stress, and arrhythmogenesis.

Abstract LB270: LINGO2-microRNA cluster controls oral cancer stemness ability and cancer progression

Abstract Cancer recurrence and metastasis are the primary reasons for treatment failure in late-stage oral cancer. Cancer stem cells are the root of cancer recurrence and metastasis. By using the microRNAome analysis of Taiwan OSCC cohort, we found miR-876-3p was highly correlated to OSCC recurrence. The precursor miR-876 promoted in vitro OSCC cell proliferation, migration, and cancer stemness and regulated several cancer stemness-relative genes such as NANOG, OCT4, and LIN28A expression. miR-876 and miR-873 are located in the second intron of the 5’-untransalation region (5’-UTR) at the alternative variant leucine-rich repeat and Ig domain containing 2 (LINGO2). Interestingly, only OSCC cells presented LINGO2 protein expression but did not show in primary oral keratinocytes (HOK). The unique CpG islands at the alternative transcription start site of LINGO2 hinted at crucial epigenetic regulation that controls LINGO2 activation. Silencing of LINGO2 expression reduced OSCC cell growth and cancer stemness abilities. LINGO2 expression was higher in OSCC tumor tissue than in adjacent normal tissue in both Taiwan and The Cancer Genome Atlas (TCGA) head and neck cancer cohorts. Moreover, LINGO2 also served as a poor prognostic marker in cancer patients' overall and recurrent free survival time. In summary, we found a novel LINGO2-microRNA cluster that regulates OSCC stemness and metastasis abilities. Those cancer specific expressions may provide a new therapy niche in OSCC cancer therapy. Citation Format: Wei-Min Chang, Li-Jie Li, Che-Hsuan Lin, Michael Hsiao, Peter Mu-Hsin Chang. LINGO2-microRNA cluster controls oral cancer stemness ability and cancer progression [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 2 (Clinical Trials and Late-Breaking Research); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(8_Suppl):Abstract nr LB270.

An epigenetic switch controls an alternative NR2F2 isoform that unleashes a metastatic program in melanoma

Metastatic melanoma develops once transformed melanocytic cells begin to de-differentiate into migratory and invasive melanoma cells with neural crest cell (NCC)-like and epithelial-to-mesenchymal transition (EMT)-like features. However, it is still unclear how transformed melanocytes assume a metastatic melanoma cell state. Here, we define DNA methylation changes that accompany metastatic progression in melanoma patients and discover Nuclear Receptor Subfamily 2 Group F, Member 2 – isoform 2 (NR2F2-Iso2) as an epigenetically regulated metastasis driver. NR2F2-Iso2 is transcribed from an alternative transcriptional start site (TSS) and it is truncated at the N-terminal end which encodes the NR2F2 DNA-binding domain. We find that NR2F2-Iso2 expression is turned off by DNA methylation when NCCs differentiate into melanocytes. Conversely, this process is reversed during metastatic melanoma progression, when NR2F2-Iso2 becomes increasingly hypomethylated and re-expressed. Our functional and molecular studies suggest that NR2F2-Iso2 drives metastatic melanoma progression by modulating the activity of full-length NR2F2 (Isoform 1) over EMT- and NCC-associated target genes. Our findings indicate that DNA methylation changes play a crucial role during metastatic melanoma progression, and their control of NR2F2 activity allows transformed melanocytes to acquire NCC-like and EMT-like features. This epigenetically regulated transcriptional plasticity facilitates cell state transitions and metastatic spread.

Smar2C2: A Simple and Efficient Protocol for the Identification of Transcription Start Sites.

Promoters and the noncoding sequences that drive their function are fundamental aspects of genes that are critical to their regulation. The transcription preinitiation complex binds and assembles on promoters where it facilitates transcription. The transcription start site (TSS) is located downstream of the promoter sequence and is defined as the location in the genome where polymerase begins transcribing DNA into RNA. Knowing the location of TSSs is useful for annotation of genes, identification of non-coding sequences important to gene regulation, detection of alternative TSSs, and understanding of 5' UTR content. Several existing techniques make it possible to accurately identify TSSs, but are often difficult to perform experimentally, require large amounts of input RNA, or are unable to identify a large number of TSSs from a single sample. Many of these protocols take advantage of template switching reverse transcriptases (TSRTs), which reliably place an adaptor at the 5' end of a first strand synthesis of cDNA. Here, we introduce a protocol that exploits TSRT activity combined with rolling circle amplification to identify TSSs with several unique advantages over existing methods. Sequence adaptors are placed on the 5' and 3' end of the full-length cDNA copy of a transcript. A splint compatible with those adaptors is then used to circularize the full-length cDNA. Linear DNA containing concatemers of the cDNA are generated using rolling circle amplification, and a sequencing library is formed by fragmenting the concatemers. This protocol is straightforward to execute, requiring limited bench time with relatively stable reagents. Using extremely low amounts of RNA input, this protocol produces large numbers of accurate, deduplicated TSSs genome wide. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Splint generation Basic Protocol 2: RNA extraction Basic Protocol 3: cDNA synthesis Basic Protocol 4: cDNA circularization and amplification Basic Protocol 5: Library generation.

Nanopore-Based Direct RNA Sequencing of the Trypanosoma brucei Transcriptome Identifies Novel lncRNAs.

Trypanosomatids are single-cell eukaryotic parasites. Unlike higher eukaryotes, they control gene expression post-transcriptionally and not at the level of transcription initiation. This involves all known cellular RNA circuits, from mRNA processing to mRNA decay, to translation, in addition to a large panel of RNA-interacting proteins that modulate mRNA abundance. However, other forms of gene regulation, for example by lncRNAs, cannot be excluded. LncRNAs are poorly studied in trypanosomatids, with only a single lncRNA characterized to date. Furthermore, it is not clear whether the complete inventory of trypanosomatid lncRNAs is known, because of the inherent cDNA-recoding and DNA-amplification limitations of short-read RNA sequencing. Here, we overcome these limitations by using long-read direct RNA sequencing (DRS) on nanopore arrays. We analyze the native RNA pool of the two main lifecycle stages of the African trypanosome Trypanosoma brucei, with a special emphasis on the inventory of lncRNAs. We identify 207 previously unknown lncRNAs, 32 of which are stage-specifically expressed. We also present insights into the complexity of the T. brucei transcriptome, including alternative transcriptional start and stop sites and potential transcript isoforms, to provide a bias-free understanding of the intricate RNA landscape in T. brucei.

Identifying transcript-level differential expression in primary human immune cells

BackgroundMultipotential hematopoietic stem cells differentiate into a wide variety of immune cells with a diversity of functions, including the ability to respond to a variety of stimuli. Importantly, numerous studies have demonstrated the importance of gene transcription in defining cell identity and functions. While these studies have primarily been performed at the level of the gene, it is known that key immune genes such as CD44 and CD45 generate multiple different transcripts that are differentially expressed across different immune cells, and that encode proteins with different sequences and functions. Prior genomic surveys have shown that the mechanisms for generating diversity in expressed transcripts (alternate splicing, alternate transcription start sites, etc.) are very active in immune cells, but have been lacking in terms of identifying genes with multiple transcripts, that are differentially expressed, and likely to affect cell functions. MethodsWe first identified the set of genes that had at least two transcripts expressed in our RNA sequencing dataset generated from purified populations of neutrophils, monocytes and five lymphocyte populations (B, NK, γδ T, CD4 + T and CD8 + T) from twelve healthy donors. Next, we developed a heuristic approach to identify genes where two or more transcripts have distinct expression patterns across lymphoid and/or myeloid populations. We then focused our annotation and interpretation on differentially expressed transcripts that affect the coding sequence. This process was repeated to identify transcripts that were differentially expressed between monocytes and populations of macrophages and LPS-stimulated macrophages derived from these monocytes in vitro. ResultsWe found that over 55 % of genes had two or more expressed transcripts, with an average ∼3 transcripts per gene, and that 70 % of these had at least two of the transcripts that encoded proteins with different sequences. As expected, we identified a complex pattern of differential expression for multiple transcripts encoding the CD45 transmembrane protein, but we also found similar evidence for ten other genes (CD300A, FYB1, GPI, LITAF, PSMA1, PTMA, RPL32, SEPTIN9, SH3BP2, SH3KBP1) when comparing the expression patterns of transcripts within myeloid and lymphoid cells. We also identified five genes with differentially expressed transcripts associated with the transition from monocytes to macrophages (FNBP1, KLF6, and SEPTIN9) or between macrophages and LPS-stimulated macrophages (CD44, OAZ2, and SEPTIN9). For the most part, we found that the different transcripts of these genes are expected to impact specific biological functions, for example the different transcripts of SEPTIN9 likely regulate the cytoskeleton in immune cells via their interactions with actins filaments and microtubules. ConclusionsThis analytic approach successfully identified multi-transcript genes that are differentially expressed across immune cells and could be applied to other transcriptomic data. Data availability statementResearchers can request access to the individual-level data from the current study by contacting the Montreal Heart Institute ethics committee at the following institutional email address: cer.icm@icm-mhi.org.