RESUMO
Alternative pre-mRNA splicing has long been proposed to contribute greatly to proteome complexity. However, the extent to which mature mRNA isoforms are successfully translated into protein remains controversial. Here, we used high-throughput RNA sequencing and mass spectrometry (MS)-based proteomics to better evaluate the translation of alternatively spliced mRNAs. To increase proteome coverage and improve protein quantitation, we optimized cell fractionation and sample processing steps at both the protein and peptide level. Furthermore, we generated a custom peptide database trained on analysis of RNA-seq data with MAJIQ, an algorithm optimized to detect and quantify differential and unannotated splice junction usage. We matched tandem mass spectra acquired by data-dependent acquisition (DDA) against our custom RNA-seq based database, as well as SWISS-PROT and RefSeq databases to improve identification of splicing-derived proteoforms by 28% compared with use of the SWISS-PROT database alone. Altogether, we identified peptide evidence for 554 alternate proteoforms corresponding to 274 genes. Our increased depth and detection of proteins also allowed us to track changes in the transcriptome and proteome induced by T-cell stimulation, as well as fluctuations in protein subcellular localization. In sum, our data here confirm that use of generic databases in proteomic studies underestimates the number of spliced mRNA isoforms that are translated into protein and provides a workflow that improves isoform detection in large-scale proteomic experiments.
Assuntos
Algoritmos , Processamento Alternativo , Bases de Dados de Ácidos Nucleicos , Perfilação da Expressão Gênica , Sequenciamento de Nucleotídeos em Larga Escala , Peptídeos , Isoformas de RNA , Humanos , Peptídeos/genética , Peptídeos/metabolismo , Isoformas de Proteínas/biossíntese , Isoformas de Proteínas/genética , Isoformas de RNA/biossíntese , Isoformas de RNA/genética , Espectrometria de Massas em TandemRESUMO
Human CD4+ T cells are often subdivided into distinct subtypes, including Th1, Th2, Th17, and Treg cells, that are thought to carry out distinct functions in the body. Typically, these T-cell subpopulations are defined by the expression of distinct gene repertoires; however, there is variability between studies regarding the methods used for isolation and the markers used to define each T-cell subtype. Therefore, how reliably studies can be compared to one another remains an open question. Moreover, previous analysis of gene expression in CD4+ T-cell subsets has largely focused on gene expression rather than alternative splicing. Here we take a meta-analysis approach, comparing eleven independent RNA-seq studies of human Th1, Th2, Th17, and/or Treg cells to determine the consistency in gene expression and splicing within each subtype across studies. We find that known master-regulators are consistently enriched in the appropriate subtype; however, cytokines and other genes often used as markers are more variable. Importantly, we also identify previously unknown transcriptomic markers that appear to consistently differentiate between subsets, including a few Treg-specific splicing patterns. Together this work highlights the heterogeneity in gene expression between samples designated as the same subtype, but also suggests additional markers that can be used to define functional groupings.
Assuntos
Linfócitos T CD4-Positivos/fisiologia , Expressão Gênica/genética , Splicing de RNA/genética , Subpopulações de Linfócitos T/fisiologia , Transcriptoma/genética , Adulto , Células Cultivadas , Citocinas/genética , Feminino , Humanos , Masculino , Pessoa de Meia-IdadeRESUMO
Cells adjust to hypoxic stress within the tumor microenvironment by downregulating energy-consuming processes including translation. To delineate mechanisms of cellular adaptation to hypoxia, we performed RNA-Seq of normoxic and hypoxic head and neck cancer cells. These data revealed a significant down regulation of genes known to regulate RNA processing and splicing. Exon-level analyses classified > 1,000 mRNAs as alternatively spliced under hypoxia and uncovered a unique retained intron (RI) in the master regulator of translation initiation, EIF2B5. Notably, this intron was expressed in solid tumors in a stage-dependent manner. We investigated the biological consequence of this RI and demonstrate that its inclusion creates a premature termination codon (PTC), that leads to a 65kDa truncated protein isoform that opposes full-length eIF2Bε to inhibit global translation. Furthermore, expression of 65kDa eIF2Bε led to increased survival of head and neck cancer cells under hypoxia, providing evidence that this isoform enables cells to adapt to conditions of low oxygen. Additional work to uncover -cis and -trans regulators of EIF2B5 splicing identified several factors that influence intron retention in EIF2B5: a weak splicing potential at the RI, hypoxia-induced expression and binding of the splicing factor SRSF3, and increased binding of total and phospho-Ser2 RNA polymerase II specifically at the intron retained under hypoxia. Altogether, these data reveal differential splicing as a previously uncharacterized mode of translational control under hypoxia and are supported by a model in which hypoxia-induced changes to cotranscriptional processing lead to selective retention of a PTC-containing intron in EIF2B5.
Assuntos
Fator de Iniciação 2B em Eucariotos/genética , Perfilação da Expressão Gênica/métodos , Íntrons/genética , Biossíntese de Proteínas/genética , Hipóxia Tumoral/genética , Processamento Alternativo/genética , Sequência de Bases , Linhagem Celular Tumoral , Regulação Neoplásica da Expressão Gênica , Loci Gênicos , Neoplasias de Cabeça e Pescoço/genética , Neoplasias de Cabeça e Pescoço/patologia , Humanos , Modelos Biológicos , Motivos de Nucleotídeos/genética , Fosforilação , Reação em Cadeia da Polimerase , Ligação Proteica , Isoformas de Proteínas/genética , Isoformas de Proteínas/metabolismo , RNA Polimerase II/metabolismo , RNA Mensageiro/genética , RNA Mensageiro/metabolismo , Sequências Reguladoras de Ácido Nucleico/genética , Reprodutibilidade dos TestesRESUMO
Glycogen synthase kinase-3 (GSK-3) is a constitutively active, ubiquitously expressed protein kinase that regulates multiple signaling pathways. In vitro kinase assays and genetic and pharmacological manipulations of GSK-3 have identified more than 100 putative GSK-3 substrates in diverse cell types. Many more have been predicted on the basis of a recurrent GSK-3 consensus motif ((pS/pT)XXX(S/T)), but this prediction has not been tested by analyzing the GSK-3 phosphoproteome. Using stable isotope labeling of amino acids in culture (SILAC) and MS techniques to analyze the repertoire of GSK-3-dependent phosphorylation in mouse embryonic stem cells (ESCs), we found that â¼2.4% of (pS/pT)XXX(S/T) sites are phosphorylated in a GSK-3-dependent manner. A comparison of WT and Gsk3a;Gsk3b knock-out (Gsk3 DKO) ESCs revealed prominent GSK-3-dependent phosphorylation of multiple splicing factors and regulators of RNA biosynthesis as well as proteins that regulate transcription, translation, and cell division. Gsk3 DKO reduced phosphorylation of the splicing factors RBM8A, SRSF9, and PSF as well as the nucleolar proteins NPM1 and PHF6, and recombinant GSK-3ß phosphorylated these proteins in vitro RNA-Seq of WT and Gsk3 DKO ESCs identified â¼190 genes that are alternatively spliced in a GSK-3-dependent manner, supporting a broad role for GSK-3 in regulating alternative splicing. The MS data also identified posttranscriptional regulation of protein abundance by GSK-3, with â¼47 proteins (1.4%) whose levels increased and â¼78 (2.4%) whose levels decreased in the absence of GSK-3. This study provides the first unbiased analysis of the GSK-3 phosphoproteome and strong evidence that GSK-3 broadly regulates alternative splicing.
Assuntos
Glicogênio Sintase Quinase 3 beta/metabolismo , Quinase 3 da Glicogênio Sintase/metabolismo , Proteínas de Homeodomínio/metabolismo , Proteínas Nucleares/metabolismo , Processamento de Proteína Pós-Traducional , Proteínas de Ligação a RNA/metabolismo , Fatores de Processamento de Serina-Arginina/metabolismo , Processamento Alternativo , Animais , Isótopos de Carbono , Linhagem Celular , Proteínas de Ligação a DNA/química , Proteínas de Ligação a DNA/metabolismo , Células-Tronco Embrionárias/enzimologia , Células-Tronco Embrionárias/metabolismo , Técnicas de Inativação de Genes , Quinase 3 da Glicogênio Sintase/genética , Glicogênio Sintase Quinase 3 beta/genética , Proteínas de Homeodomínio/química , Camundongos , Isótopos de Nitrogênio , Proteínas Nucleares/química , Nucleofosmina , Mapeamento de Peptídeos , Fosforilação , Estabilidade Proteica , Proteômica/métodos , Proteínas de Ligação a RNA/química , Proteínas Recombinantes/metabolismo , Proteínas Repressoras , Fatores de Processamento de Serina-Arginina/química , Especificidade por SubstratoRESUMO
Microbial colonization of the infant gastrointestinal tract (GIT) begins at birth, is shaped by the maternal microbiota, and is profoundly altered by antibiotic treatment. Antibiotic treatment of mothers during pregnancy influences colonization of the GIT microbiota of their infants. The role of the GIT microbiota in regulating adaptive immune function against systemic viral infections during infancy remains undefined. We used a mouse model of perinatal antibiotic exposure to examine the effect of GIT microbial dysbiosis on infant CD8(+) T cell-mediated antiviral immunity. Maternal antibiotic treatment/treated (MAT) during pregnancy and lactation resulted in profound alterations in the composition of the GIT microbiota in mothers and infants. Streptococcus spp. dominated the GIT microbiota of MAT mothers, whereas Enterococcus faecalis predominated within the MAT infant GIT. MAT infant mice subsequently exhibited increased and accelerated mortality following vaccinia virus infection. Ag-specific IFN-γ-producing CD8(+) T cells were reduced in sublethally infected MAT infant mice. MAT CD8(+) T cells from uninfected infant mice also demonstrated a reduced capacity to sustain IFN-γ production following in vitro activation. We additionally determined that control infant mice became more susceptible to infection if they were born in an animal facility using stricter standards of hygiene. These data indicate that undisturbed colonization and progression of the GIT microbiota during infancy are necessary to promote robust adaptive antiviral immune responses.
Assuntos
Linfócitos T CD8-Positivos/imunologia , Enterococcus faecalis/fisiologia , Microbioma Gastrointestinal , Streptococcus/fisiologia , Vaccinia virus/imunologia , Vacínia/microbiologia , Imunidade Adaptativa , Animais , Animais Recém-Nascidos , Antibacterianos/administração & dosagem , Antibacterianos/efeitos adversos , Células Cultivadas , Feminino , Interferon gama/metabolismo , Exposição Materna/efeitos adversos , Camundongos , Camundongos Endogâmicos C57BL , Vacínia/imunologiaRESUMO
GB virus C (GBV-C) is a non-pathogenic flavivirus that may play a role in modulating HIV disease. Multiple genotypes of GBV-C that have been identified to date that may differentially regulate HIV; however, the number of complete GBV-C sequences published to date is very limited. We sequenced full-length GBV-C genomes from four individuals with HIV/HCV co-infection in the United States. Intergenotypic recombination was evident in two of these individuals. Evaluation of additional full-length GBV-C genomes would facilitate the creation of full-length, replication-competent molecular clones of GBV-C to evaluate the phenotypic diversity of GBV-C genotypes and provide important molecular data on this understudied virus.
Assuntos
Vírus GB C/genética , Vírus GB C/isolamento & purificação , Genoma Viral/genética , Recombinação Genética/genética , Sequência de Aminoácidos , Sequência de Bases , Coinfecção , Infecções por Flaviviridae/virologia , Humanos , Filogenia , Estudos Prospectivos , RNA Viral/genética , Análise de Sequência de RNA , Estados UnidosRESUMO
Relapsed or refractory B-cell acute lymphoblastic leukemia (B-ALL) is a major cause of pediatric cancer-related deaths. Relapse-specific mutations do not account for all chemotherapy failures in B- ALL patients, suggesting additional mechanisms of resistance. By mining RNA-seq datasets of paired diagnostic/relapse pediatric B-ALL samples, we discovered pervasive alternative splicing (AS) patterns linked to relapse and affecting drivers of resistance to glucocorticoids, anti-folates, and thiopurines. Most splicing variations represented cassette exon skipping, "poison" exon inclusion, and intron retention, phenocopying well-documented loss-of-function mutations. In contrast, relapse-associated AS of NT5C2 mRNA yielded an isoform with the functionally uncharacterized in-frame exon 6a. Incorporation of the 8-amino acid sequence SQVAVQKR into this enzyme created a putative phosphorylation site and resulted in elevated nucleosidase activity, which is a known consequence of gain-of-function mutations in NT5C2 and a common determinant of 6-mercaptopurine (6-MP) resistance. Consistent with this finding, NT5C2ex6a and the R238W hotspot variant conferred comparable levels of resistance to 6-MP in B-ALL cells both in vitro and in vivo. Furthermore, both the NT5C2ex6a and R238W variants induced collateral sensitivity to the inosine monophosphate dehydrogenase (IMPDH) inhibitor mizoribine. These results ascribe an important role for splicing perturbations in chemotherapy resistance in relapsed B-ALL and suggest that IMPDH inhibitors, including the commonly used immunosuppressive agent mycophenolate mofetil, could be a valuable therapeutic option for treating thiopurine-resistant leukemias.
RESUMO
Relapsed or refractory B-cell acute lymphoblastic leukemia (B-ALL) is a major cause of pediatric cancer-related deaths. Relapse-specific mutations do not account for all chemotherapy failures in B-ALL patients, suggesting additional mechanisms of resistance. By mining RNA-seq datasets of paired diagnostic/relapse pediatric B-ALL samples, we discovered pervasive alternative splicing (AS) patterns linked to relapse and affecting drivers of resistance to glucocorticoids, anti-folates, and thiopurines. Most splicing variations represented cassette exon skipping, "poison" exon inclusion, and intron retention, phenocopying well-documented loss-of-function mutations. In contrast, relapse-associated AS of NT5C2 mRNA yielded an isoform with the functionally uncharacterized in-frame exon 6a. Incorporation of the 8-amino acid sequence SQVAVQKR into this enzyme created a putative phosphorylation site and resulted in elevated nucleosidase activity, which is a known consequence of gain-of-function mutations in NT5C2 and a common determinant of 6-mercaptopurine (6-MP) resistance. Consistent with this finding, NT5C2ex6a and the R238W hotspot variant conferred comparable levels of resistance to 6-MP in B-ALL cells both in vitro and in vivo. Furthermore, both the NT5C2ex6a and R238W variants induced collateral sensitivity to the inosine monophosphate dehydrogenase (IMPDH) inhibitor mizoribine. These results ascribe an important role for splicing perturbations in chemotherapy resistance in relapsed B-ALL and suggest that IMPDH inhibitors, including the commonly used immunosuppressive agent mycophenolate mofetil, could be a valuable therapeutic option for treating thiopurine-resistant leukemias.
RESUMO
The ubiquity of RNA-seq has led to many methods that use RNA-seq data to analyze variations in RNA splicing. However, available methods are not well suited for handling heterogeneous and large datasets. Such datasets scale to thousands of samples across dozens of experimental conditions, exhibit increased variability compared to biological replicates, and involve thousands of unannotated splice variants resulting in increased transcriptome complexity. We describe here a suite of algorithms and tools implemented in the MAJIQ v2 package to address challenges in detection, quantification, and visualization of splicing variations from such datasets. Using both large scale synthetic data and GTEx v8 as benchmark datasets, we assess the advantages of MAJIQ v2 compared to existing methods. We then apply MAJIQ v2 package to analyze differential splicing across 2,335 samples from 13 brain subregions, demonstrating its ability to offer insights into brain subregion-specific splicing regulation.
Assuntos
Algoritmos , Splicing de RNA , RNA-Seq , Benchmarking , EncéfaloRESUMO
Alternative splicing occurs in the vast majority of human genes, giving rise to distinct mRNA and protein isoforms. We, and others, have previously identified hundreds of genes that change their isoform expression upon T cell activation via alternative splicing; however, how these changes link activation input with functional output remains largely unknown. Here, we investigate how costimulation of T cells through the CD28 receptor impacts alternative splicing in T cells activated through the T cell receptor (TCR, CD3) and find that while CD28 signaling alone has minimal impact on splicing, it enhances the extent of change for up to 20% of TCR-induced alternative splicing events. Interestingly, a set of CD28-enhanced splicing events occur within genes encoding key components of the apoptotic signaling pathway; namely caspase-9, Bax, and Bim. Using both CRISPR-edited cells and antisense oligos to force expression of specific isoforms, we show for all three of these genes that the isoform induced by CD3/CD28 costimulation promotes resistance to apoptosis, and that changes in all three genes together function combinatorially to further promote cell viability. Finally, we show that the JNK signaling pathway, induced downstream of CD3/CD28 costimulation, is required for each of these splicing events, further highlighting their co-regulation. Together, these findings demonstrate that alternative splicing is a key mechanism by which costimulation of CD28 promotes viability of activated T cells.
Assuntos
Antígenos CD28 , Linfócitos T , Humanos , Linfócitos T/metabolismo , Antígenos CD28/metabolismo , Processamento Alternativo , Sobrevivência Celular , Receptores de Antígenos de Linfócitos T/metabolismo , ApoptoseRESUMO
The effects of confounding factors on gene expression analysis have been extensively studied following the introduction of high-throughput microarrays and subsequently RNA sequencing. In contrast, there is a lack of equivalent analysis and tools for RNA splicing. Here we first assess the effect of confounders on both expression and splicing quantifications in two large public RNA-Seq datasets (TARGET, ENCODE). We show quantification of splicing variations are affected at least as much as those of gene expression, revealing unwanted sources of variations in both datasets. Next, we develop MOCCASIN, a method to correct the effect of both known and unknown confounders on RNA splicing quantification and demonstrate MOCCASIN's effectiveness on both synthetic and real data. Code, synthetic and corrected datasets are all made available as resources.
Assuntos
Algoritmos , Biologia Computacional/métodos , Perfilação da Expressão Gênica/métodos , Sequenciamento de Nucleotídeos em Larga Escala/métodos , Splicing de RNA , Bases de Dados Genéticas , Humanos , Células K562 , RNA-Seq/métodos , Reprodutibilidade dos Testes , SoftwareRESUMO
Alternative splicing contributes to gene expression dynamics in many tissues, yet its role in auditory development remains unclear. We performed whole-exome sequencing in individuals with sensorineural hearing loss (SNHL) and identified pathogenic mutations in Epithelial Splicing-Regulatory Protein 1 (ESRP1). Patient-derived induced pluripotent stem cells showed alternative splicing defects that were restored upon repair of an ESRP1 mutant allele. To determine how ESRP1 mutations cause hearing loss, we evaluated Esrp1-/- mouse embryos and uncovered alterations in cochlear morphogenesis, auditory hair cell differentiation, and cell fate specification. Transcriptome analysis revealed impaired expression and splicing of genes with essential roles in cochlea development and auditory function. Aberrant splicing of Fgfr2 blocked stria vascularis formation due to erroneous ligand usage, which was corrected by reducing Fgf9 gene dosage. These findings implicate mutations in ESRP1 as a cause of SNHL and demonstrate the complex interplay between alternative splicing, inner ear development, and auditory function.
Assuntos
Processamento Alternativo/genética , Cóclea/embriologia , Perda Auditiva/genética , Mutação/genética , Proteínas de Ligação a RNA/genética , Animais , Diferenciação Celular/genética , Cóclea/metabolismo , Camundongos KnockoutRESUMO
Genetically engineered microbial biosensors have yet to realize commercial success in environmental applications due, in part, to difficulties associated with transducing and transmitting traditional bioluminescent information. Bioelectrochemical systems (BESs) output a direct electric signal that can be incorporated into devices for remote environmental monitoring. Here, we describe a BES-based biosensor with genetically encoded specificity for a toxic metal. By placing an essential component of the metal reduction (Mtr) pathway of Shewanella oneidensis under the control of an arsenic-sensitive promoter, we have genetically engineered a strain that produces increased current in response to arsenic when inoculated into a BES. Our BES-based biosensor has a detection limit of ~40 µM arsenite with a linear range up to 100 µM arsenite. Because our transcriptional circuit relies on the activation of a single promoter, similar sensing systems may be developed to detect other analytes by the swap of a single genetic part.