Proteogenomic analysis of chemo-refractory high-grade serous ovarian cancer.

To improve the understanding of chemo-refractory high-grade serous ovarian cancers (HGSOCs), we characterized the proteogenomic landscape of 242 (refractory and sensitive) HGSOCs, representing one discovery and two validation cohorts across two biospecimen types (formalin-fixed paraffin-embedded and frozen). We identified a 64-protein signature that predicts with high specificity a subset of HGSOCs refractory to initial platinum-based therapy and is validated in two independent patient cohorts. We detected significant association between lack of Ch17 loss of heterozygosity (LOH) and chemo-refractoriness. Based on pathway protein expression, we identified 5 clusters of HGSOC, which validated across two independent patient cohorts and patient-derived xenograft (PDX) models. These clusters may represent different mechanisms of refractoriness and implicate putative therapeutic vulnerabilities.

A distinct core regulatory module enforces oncogene expression in KMT2A-rearranged leukemia.

Acute myeloid leukemia with KMT2A (MLL) rearrangements is characterized by specific patterns of gene expression and enhancer architecture, implying unique core transcriptional regulatory circuitry. Here, we identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. Rapid transcription factor degradation followed by measurements of genome-wide transcription rates and superresolution microscopy revealed that MEF2D and IRF8 form a distinct core regulatory module with a narrow direct transcriptional program that includes activation of the key oncogenes MYC, HOXA9, and BCL2. Our study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Extensive Recovery of Embryonic Enhancer and Gene Memory Stored in Hypomethylated Enhancer DNA.

Developing and adult tissues use different cis-regulatory elements. Although DNA at some decommissioned embryonic enhancers is hypomethylated in adult cells, it is unknown whether this putative epigenetic memory is complete and recoverable. We find that, in adult mouse cells, hypomethylated CpG dinucleotides preserve a nearly complete archive of tissue-specific developmental enhancers. Sites that carry the active histone mark H3K4me1, and are therefore considered "primed," are mainly cis elements that act late in organogenesis. In contrast, sites decommissioned early in development retain hypomethylated DNA as a singular property. In adult intestinal and blood cells, sustained absence of polycomb repressive complex 2 indirectly reactivates most-and only-hypomethylated developmental enhancers. Embryonic and fetal transcriptional programs re-emerge as a result, in reverse chronology to cis element inactivation during development. Thus, hypomethylated DNA in adult cells preserves a "fossil record" of tissue-specific developmental enhancers, stably marking decommissioned sites and enabling recovery of this epigenetic memory.

The KDM5A/RBP2 histone demethylase represses NOTCH signaling to sustain neuroendocrine differentiation and promote small cell lung cancer tumorigenesis.

More than 90% of small cell lung cancers (SCLCs) harbor loss-of-function mutations in the tumor suppressor gene RB1 The canonical function of the RB1 gene product, pRB, is to repress the E2F transcription factor family, but pRB also functions to regulate cellular differentiation in part through its binding to the histone demethylase KDM5A (also known as RBP2 or JARID1A). We show that KDM5A promotes SCLC proliferation and SCLC's neuroendocrine differentiation phenotype in part by sustaining expression of the neuroendocrine transcription factor ASCL1. Mechanistically, we found that KDM5A sustains ASCL1 levels and neuroendocrine differentiation by repressing NOTCH2 and NOTCH target genes. To test the role of KDM5A in SCLC tumorigenesis in vivo, we developed a CRISPR/Cas9-based mouse model of SCLC by delivering an adenovirus (or an adeno-associated virus [AAV]) that expresses Cre recombinase and sgRNAs targeting Rb1, Tp53, and Rbl2 into the lungs of Lox-Stop-Lox Cas9 mice. Coinclusion of a KDM5A sgRNA decreased SCLC tumorigenesis and metastasis, and the SCLCs that formed despite the absence of KDM5A had higher NOTCH activity compared to KDM5A+/+ SCLCs. This work establishes a role for KDM5A in SCLC tumorigenesis and suggests that KDM5 inhibitors should be explored as treatments for SCLC.

Comprehensive genetic and phenotype analysis of 95 individuals with mosaic tuberous sclerosis complex.

Tuberous sclerosis complex (TSC) is a neurogenetic disorder due to loss-of-function TSC1 or TSC2 variants, characterized by tumors affecting multiple organs, including skin, brain, heart, lung, and kidney. Mosaicism for TSC1 or TSC2 variants occurs in 10%-15% of individuals diagnosed with TSC. Here, we report comprehensive characterization of TSC mosaicism by using massively parallel sequencing (MPS) of 330 TSC samples from a variety of tissues and fluids from a cohort of 95 individuals with mosaic TSC. TSC1 variants in individuals with mosaic TSC are much less common (9%) than in germline TSC overall (26%) (p < 0.0001). The mosaic variant allele frequency (VAF) is significantly higher in TSC1 than in TSC2, in both blood and saliva (median VAF: TSC1, 4.91%; TSC2, 1.93%; p = 0.036) and facial angiofibromas (median VAF: TSC1, 7.7%; TSC2 3.7%; p = 0.004), while the number of TSC clinical features in individuals with TSC1 and TSC2 mosaicism was similar. The distribution of mosaic variants across TSC1 and TSC2 is similar to that for pathogenic germline variants in general TSC. The systemic mosaic variant was not present in blood in 14 of 76 (18%) individuals with TSC, highlighting the value of analysis of multiple samples from each individual. A detailed comparison revealed that nearly all TSC clinical features are less common in individuals with mosaic versus germline TSC. A large number of previously unreported TSC1 and TSC2 variants, including intronic and large rearrangements (n = 11), were also identified.

Impairment of gamma-glutamyl transferase 1 activity in the metabolic pathogenesis of chromophobe renal cell carcinoma.

Chromophobe renal cell carcinoma (ChRCC) accounts for 5% of all sporadic renal cancers and can also occur in genetic syndromes including Birt-Hogg-Dube (BHD) and tuberous sclerosis complex (TSC). ChRCC has a distinct accumulation of abnormal mitochondria, accompanied by characteristic chromosomal imbalances and relatively few "driver" mutations. Metabolomic profiling of ChRCC and oncocytomas (benign renal tumors that share pathological features with ChRCC) revealed both similarities and differences between these tumor types, with principal component analysis (PCA) showing a distinct separation. ChRCC have a striking decrease in intermediates of the glutathione salvage pathway (also known as the gamma-glutamyl cycle) compared with adjacent normal kidney, as well as significant changes in glycolytic and pentose phosphate pathway intermediates. We also found that gamma glutamyl transferase 1 (GGT1), the key enzyme of the gamma-glutamyl cycle, is expressed at ∼100-fold lower levels in ChRCC compared with normal kidney, while no change in GGT1 expression was found in clear cell RCC (ccRCC). Significant differences in specific metabolite abundance were found in ChRCC vs. ccRCC, including the oxidative stress marker ophthalmate. Down-regulation of GGT1 enhanced the sensitivity to oxidative stress and treatment with buthionine sulfoximine (BSO), which was associated with changes in glutathione-pathway metabolites. These data indicate that impairment of the glutathione salvage pathway, associated with enhanced oxidative stress, may have key therapeutic implications for this rare tumor type for which there are currently no specific targeted therapies.

VIPER: Visualization Pipeline for RNA-seq, a Snakemake workflow for efficient and complete RNA-seq analysis.

BACKGROUND: RNA sequencing has become a ubiquitous technology used throughout life sciences as an effective method of measuring RNA abundance quantitatively in tissues and cells. The increase in use of RNA-seq technology has led to the continuous development of new tools for every step of analysis from alignment to downstream pathway analysis. However, effectively using these analysis tools in a scalable and reproducible way can be challenging, especially for non-experts. RESULTS: Using the workflow management system Snakemake we have developed a user friendly, fast, efficient, and comprehensive pipeline for RNA-seq analysis. VIPER (Visualization Pipeline for RNA-seq analysis) is an analysis workflow that combines some of the most popular tools to take RNA-seq analysis from raw sequencing data, through alignment and quality control, into downstream differential expression and pathway analysis. VIPER has been created in a modular fashion to allow for the rapid incorporation of new tools to expand the capabilities. This capacity has already been exploited to include very recently developed tools that explore immune infiltrate and T-cell CDR (Complementarity-Determining Regions) reconstruction abilities. The pipeline has been conveniently packaged such that minimal computational skills are required to download and install the dozens of software packages that VIPER uses. CONCLUSIONS: VIPER is a comprehensive solution that performs most standard RNA-seq analyses quickly and effectively with a built-in capacity for customization and expansion.

Cross-site comparison of ribosomal depletion kits for Illumina RNAseq library construction.

BACKGROUND: Ribosomal RNA (rRNA) comprises at least 90% of total RNA extracted from mammalian tissue or cell line samples. Informative transcriptional profiling using massively parallel sequencing technologies requires either enrichment of mature poly-adenylated transcripts or targeted depletion of the rRNA fraction. The latter method is of particular interest because it is compatible with degraded samples such as those extracted from FFPE and also captures transcripts that are not poly-adenylated such as some non-coding RNAs. Here we provide a cross-site study that evaluates the performance of ribosomal RNA removal kits from Illumina, Takara/Clontech, Kapa Biosystems, Lexogen, New England Biolabs and Qiagen on intact and degraded RNA samples. RESULTS: We find that all of the kits are capable of performing significant ribosomal depletion, though there are differences in their ease of use. All kits were able to remove ribosomal RNA to below 20% with intact RNA and identify ~ 14,000 protein coding genes from the Universal Human Reference RNA sample at >1FPKM. Analysis of differentially detected genes between kits suggests that transcript length may be a key factor in library production efficiency. CONCLUSIONS: These results provide a roadmap for labs on the strengths of each of these methods and how best to utilize them.

The association between germline BRCA2 variants and sensitivity to platinum-based chemotherapy among men with metastatic prostate cancer.

BACKGROUND: Breast cancer 2 (BRCA2)-associated breast and ovarian cancers are sensitive to platinum-based chemotherapy. It is unknown whether BRCA2-associated prostate cancer responds favorably to such treatment. METHODS: A retrospective analysis of a single-institution cohort of men with castration-resistant, metastatic prostate cancer was performed to determine the association between carrier status of pathogenic BRCA2 germline variants and prostate-specific antigen response to carboplatin-based chemotherapy. From 2001 through 2015, 8081 adult men with prostate cancer who had a consultation and/or underwent treatment at Dana-Farber Cancer Institute provided blood samples and consented to analyses of biologic material and clinical records. A subgroup of 141 men received at least 2 doses of carboplatin and docetaxel for castration-resistant disease (94% were also taxane refractory). These patients were categorized according to the absence or presence of pathogenic germline mutations in BRCA2 based on DNA sequencing from whole blood. The primary outcome was the response rate to carboplatin/docetaxel chemotherapy, defined according to a decline in prostate-specific antigen that exceeded 50% within 12 weeks of initiating this regimen. Associations between BRCA2 mutation status and response to carboplatin-based chemotherapy were tested using the Fisher exact test, with a 2-sided P value < .05 as the threshold for significance. RESULTS: Pathogenic germline BRCA2 variants were observed in 8 of 141 men (5.7%; 95% confidence interval, 2.5%-10.9%). Six of 8 BRCA2 carriers (75%) experienced prostate-specific antigen declines >50% within 12 weeks, compared with 23 of 133 noncarriers (17%; absolute difference, 58%; 95% confidence interval, 27%-88%; P < .001). Prostate cancer cell lines functionally corroborated these clinical findings. CONCLUSIONS: BRCA2-associated, castration-resistant prostate cancer is associated with a higher likelihood of response to carboplatin-based chemotherapy than non-BRCA2-associated prostate cancer. Cancer 2017;123:3532-9. © 2017 American Cancer Society.

Lineage-specific canonical and non-canonical activity of EZH2 in advanced prostate cancer subtypes.

Enhancer of zeste homolog 2 (EZH2) is a histone methyltransferase and emerging therapeutic target that is overexpressed in most castration-resistant prostate cancers and implicated as a driver of disease progression and resistance to hormonal therapies. Here we define the lineage-specific action and differential activity of EZH2 in both prostate adenocarcinoma and neuroendocrine prostate cancer (NEPC) subtypes of advanced prostate cancer to better understand the role of EZH2 in modulating differentiation, lineage plasticity, and to identify mediators of response and resistance to EZH2 inhibitor therapy. Mechanistically, EZH2 modulates bivalent genes that results in upregulation of NEPC-associated transcriptional drivers (e.g., ASCL1) and neuronal gene programs in NEPC, and leads to forward differentiation after targeting EZH2 in NEPC. Subtype-specific downstream effects of EZH2 inhibition on cell cycle genes support the potential rationale for co-targeting cyclin/CDK to overcome resistance to EZH2 inhibition.

Lineage-specific canonical and non-canonical activity of EZH2 in advanced prostate cancer subtypes.

Enhancer of zeste homolog 2 (EZH2) is a histone methyltransferase and emerging therapeutic target that is overexpressed in most castration-resistant prostate cancers and implicated as a driver of disease progression and resistance to hormonal therapies. Here we define the lineage-specific action and differential activity of EZH2 in both prostate adenocarcinoma (PRAD) and neuroendocrine prostate cancer (NEPC) subtypes of advanced prostate cancer to better understand the role of EZH2 in modulating differentiation, lineage plasticity, and to identify mediators of response and resistance to EZH2 inhibitor therapy. Mechanistically, EZH2 modulates bivalent genes that results in upregulation of NEPC-associated transcriptional drivers (e.g., ASCL1) and neuronal gene programs, and leads to forward differentiation after targeting EZH2 in NEPC. Subtype-specific downstream effects of EZH2 inhibition on cell cycle genes support the potential rationale for co-targeting cyclin/CDK to overcome resistance to EZH2 inhibition.

Telomerase reactivation induces progression of mouse Braf V600E -driven thyroid cancers without telomere lengthening.

Mutations in the promoter of the telomerase reverse transcriptase ( TERT ) gene are the paradigm of a cross-cancer alteration in a non-coding region. TERT promoter mutations (TPMs) are biomarkers of poor prognosis in several tumors, including thyroid cancers. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert -123C>T ) and crossed it with thyroid-specific Braf V600E -mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all Braf V600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of Braf V600E +Tert -123C>T and Braf V600E +K5-Tert mice progressed to poorly differentiated thyroid cancers at week 20, respectively. Braf+Tert tumors showed increased mitosis and necrosis in areas of solid growth, and older animals from these cohorts displayed anaplastic-like features, i.e., spindle cells and macrophage infiltration. Murine Tert promoter mutation increased Tert transcription in vitro and in vivo , but temporal and intra-tumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine and chemokine signaling, were overactivated. Braf+Tert animals remained responsive to MAPK pathway inhibitors. These models constitute useful pre-clinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs.

Telomerase Upregulation Induces Progression of Mouse BrafV600E-Driven Thyroid Cancers and Triggers Nontelomeric Effects.

Mutations in the promoter of the telomerase reverse transcriptase (TERT) gene are the paradigm of a cross-cancer alteration in a noncoding region. TERT promoter mutations (TPM) are biomarkers of poor prognosis in cancer, including thyroid tumors. TPMs enhance TERT transcription, which is otherwise silenced in adult tissues, thus reactivating a bona fide oncoprotein. To study TERT deregulation and its downstream consequences, we generated a Tert mutant promoter mouse model via CRISPR/Cas9 engineering of the murine equivalent locus (Tert-123C>T) and crossed it with thyroid-specific BrafV600E-mutant mice. We also employed an alternative model of Tert overexpression (K5-Tert). Whereas all BrafV600E animals developed well-differentiated papillary thyroid tumors, 29% and 36% of BrafV600E+Tert-123C>T and BrafV600E+K5-Tert mice progressed to poorly differentiated cancers at week 20, respectively. Tert-upregulated tumors showed increased mitosis and necrosis in areas of solid growth, and older animals displayed anaplastic-like features, that is, spindle cells and macrophage infiltration. Murine TPM increased Tert transcription in vitro and in vivo, but temporal and intratumoral heterogeneity was observed. RNA-sequencing of thyroid tumor cells showed that processes other than the canonical Tert-mediated telomere maintenance role operate in these specimens. Pathway analysis showed that MAPK and PI3K/AKT signaling, as well as processes not previously associated with this tumor etiology, involving cytokine, and chemokine signaling, were overactivated. These models constitute useful preclinical tools to understand the cell-autonomous and microenvironment-related consequences of Tert-mediated progression in advanced thyroid cancers and other aggressive tumors carrying TPMs. IMPLICATIONS: Telomerase-driven cancer progression activates pathways that can be dissected and perhaps therapeutically exploited.

Uncovering transcriptomic landscape alterations of CAN-2409 in in vitro and in vivo glioma models.

Rationale: CAN-2409 is a locally delivered oncolytic therapy, which results in vaccination against the injected tumor. CAN-2409 consists of a non-replicating adenovirus armed with the Herpes virus thymidine kinase, which metabolizes ganciclovir into a phosphorylated nucleotide that is incorporated into the tumor cell's genome, thereby inflicting immunogenic cancer cell death. While CAN-2409's immunological impact has been well characterized, its effects on the tumor cells transcriptome remains unknown. We compared the transcriptomic landscape after treatment of glioblastoma models with CAN-2409 in vitro and in vivo to assess how the interplay with the tumor microenvironment influences CAN-2409-mediated transcriptome alterations. Methods: We performed RNA-Seq with CAN-2409 treated patient-derived glioma stem-like cells and tumors of C57/BL6 mice and compared KEGG pathway usage and differential gene expression focusing on immune cell and cytokine profiles. T-cell -killing assays were performed to assess candidate effectors. Results: PCA analysis showed distinct clustering of control and CAN-2409 samples under both conditions. KEGG pathway analysis revealed significant enrichment for p53 signaling and cell cycle pathway, with similar dynamics for key regulators of both pathways in vitro and in vivo, including MYC, CCNB1, PLK1 and CDC20. Selected alterations (PLK1 and CCNB1) were validated at the protein level. Cytokine expression analysis revealed upregulation of pro-inflammatory IL12a under both conditions; immune cell gene profiling showed reduction of myeloid associated genes. T-cell-killing assays showed increased killing in the presence of IL-12. Conclusion: CAN-2409 significantly alters the transcriptome both in vitro and in vivo. Comparison of pathway enrichment revealed mutual and differential utilization of pathways under both conditions, suggesting a modulating influence on the cell cycle in tumor cells, and of the tumor microenvironment on the transcriptome in vivo. IL-12 synthesis likely depends on interactions with the tumor microenvironment, and it facilitates CAN-2409 cell killing. This dataset provides potential to understand resistance mechanisms and identify potential biomarkers for future studies.

Androgen receptor inhibition suppresses anti-tumor neutrophil response against bone metastatic prostate cancer via regulation of TßRI expression.

Bone metastatic disease of prostate cancer (PCa) is incurable and progression in bone is largely dictated by tumor-stromal interactions in the bone microenvironment. We showed previously that bone neutrophils initially inhibit bone metastatic PCa growth yet metastatic PCa becomes resistant to neutrophil response. Further, neutrophils isolated from tumor-bone lost their ability to suppress tumor growth through unknown mechanisms. With this study, our goal was to define the impact of metastatic PCa on neutrophil function throughout tumor progression and to determine the potential of neutrophils as predictive biomarkers of metastatic disease. Using patient peripheral blood polymorphonuclear neutrophils (PMNs), we identified that PCa progression dictates PMN cell surface markers and gene expression, but not cytotoxicity against PCa. Importantly, we also identified a novel phenomenon in which second generation androgen deprivation therapy (ADT) suppresses PMN cytotoxicity via increased transforming growth factor beta receptor I (TßRI). High dose testosterone and genetic or pharmacologic TßRI inhibition rescued androgen receptor-mediated neutrophil suppression and restored neutrophil anti-tumor immune response. These studies highlight the ability to leverage standard-care ADT to generate neutrophil anti-tumor responses against bone metastatic PCa.

Heterogeneity and transcriptional drivers of triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is a heterogeneous disease with limited treatment options. To characterize TNBC heterogeneity, we defined transcriptional, epigenetic, and metabolic subtypes and subtype-driving super-enhancers and transcription factors by combining functional and molecular profiling with computational analyses. Single-cell RNA sequencing revealed relative homogeneity of the major transcriptional subtypes (luminal, basal, and mesenchymal) within samples. We found that mesenchymal TNBCs share features with mesenchymal neuroblastoma and rhabdoid tumors and that the PRRX1 transcription factor is a key driver of these tumors. PRRX1 is sufficient for inducing mesenchymal features in basal but not in luminal TNBC cells via reprogramming super-enhancer landscapes, but it is not required for mesenchymal state maintenance or for cellular viability. Our comprehensive, large-scale, multiplatform, multiomics study of both experimental and clinical TNBC is an important resource for the scientific and clinical research communities and opens venues for future investigation.

Splicing factor SRSF1 controls distinct molecular programs in regulatory and effector T cells implicated in systemic autoimmune disease.

Systemic autoimmune diseases are characterized by hyperactive effector T cells (Teffs), aberrant cytokines and chemokines, and dysfunctional regulatory T cells (Tregs). We previously uncovered new roles for serine/arginine-rich splicing factor 1 (SRSF1) in the control of genes involved in T cell signaling and cytokine production in human T cells. SRSF1 levels are decreased in T cells from patients with systemic lupus erythematosus (SLE), and low levels correlate with severe disease. Moreover, T cell-conditional Srsf1-deficient mice recapitulate the autoimmune phenotype, exhibiting CD4 T cell hyperactivity, dysfunctional Tregs, systemic autoimmunity, and tissue inflammation. However, the role of SRSF1 in controlling molecular programs in Teffs and Tregs and how these pathways are implicated in autoimmunity is not known. Here, by comparative bioinformatics analysis, we demonstrate that SRSF1 controls largely distinct gene programs in Tregs and Teffs in vivo. SRSF1 regulates 189 differentially expressed genes (DEGs) unique to Tregs, 582 DEGs unique to Teffs, and 29 DEGs shared between both. Shared genes included IL-17A, IL-17F, CSF1, CXCL10, and CXCR4, and were highly enriched for inflammatory response and cytokine-cytokine receptor interaction pathways. SRSF1 controls distinct pathways in Tregs, which include chemokine signaling and immune cell differentiation, compared with pathways in Teffs, which include cytokine production, T cell homeostasis, and activation. We identified putative mRNA binding targets of SRSF1 which include CSF1, CXCL10, and IL-17F. Finally, comparisons with transcriptomics profiles from lupus-prone MRL/lpr mice reveal that SRSF1 controls genes and pathways implicated in autoimmune disease. The target genes of SRSF1 and putative binding targets we discovered, have known roles in systemic autoimmunity. Our findings suggest that SRSF1 controls distinct molecular pathways in Tregs and Teffs and aberrant SRSF1 levels may contribute to their dysfunction and immunopathogenesis of systemic autoimmune disease.

Splicing factor SRSF1 controls autoimmune-related molecular pathways in regulatory T cells distinct from FoxP3.

Regulatory T cells (Tregs) are vital for maintaining immune self-tolerance, and their impaired function leads to autoimmune disease. Mutations in FoxP3, the master transcriptional regulator of Tregs, leads to immune dysregulation, polyendocrinopathy, enteropathy, X-linked (IPEX) syndrome in humans and the early lethal "scurfy" phenotype with multi-organ autoimmune disease in mice. We recently identified serine/arginine-rich splicing factor 1 (SRSF1) as an indispensable regulator of Treg homeostasis and function. Intriguingly, Treg-conditional SRSF1-deficient mice exhibit early lethal systemic autoimmunity with multi-organ inflammation reminiscent of the scurfy mice. Importantly, SRSF1 is decreased in T cells from patients with the autoimmune disease systemic lupus erythematosus (SLE), and low SRSF1 levels inversely correlate with disease severity. Given that the Treg-specific deficiency of SRSF1 causes similarly profound autoimmune disease outcomes in mice as the deficiency/mutation in FoxP3, we aimed to evaluate the genes and molecular pathways controlled by these two indispensable regulatory proteins. We performed comparative bioinformatic analyses of transcriptomic profiles of Tregs from Srsf1-knockout mice and two Foxp3 mutant mice--the FoxP3-deficient ΔFoxp3 and the Foxp3 M370I mutant mice. We identified 132 differentially expressed genes (DEGs) unique to Srsf1-ko Tregs, 503 DEGs unique to Foxp3 M370I Tregs, and 1367 DEGs unique to ΔFoxp3 Tregs. Gene set enrichment and pathway analysis of DEGs unique to Srsf1-ko Tregs indicate that SRSF1 controls cytokine and immune response pathways. Conversely, FoxP3 controls pathways involved in DNA replication and cell cycle. Besides the distinct gene signatures, we identified only 30 shared genes between all three Treg mutants, mostly contributing to cytokine and immune defense pathways. Prominent genes included the chemokines CXCR6 and CCL1 and the checkpoint inhibitors FASLG and PDCD1. Thus, we demonstrate that SRSF1 and FoxP3 control common and distinct molecular pathways implicated in autoimmunity. Our analyses suggest that SRSF1 controls crucial immune functions in Tregs contributing to immune tolerance, and perturbations in its levels lead to systemic autoimmunity via mechanisms that are largely distinct from FoxP3.

Splicing factor SRSF1 is essential for CD8 T cell function and host antigen-specific viral immunity.

Cytotoxic CD8 T cells are crucial for the host antigen-specific immune response to viral pathogens. Here we report the identification of an essential role for the serine/arginine-rich splicing factor (SRSF) 1 in CD8 T cell homeostasis and function. Specifically, SRSF1 is necessary for the maintenance of normal CD8 T lymphocyte numbers in the lymphoid compartment, and for the proliferative capacity and cytotoxic function of CD8 T cells. Furthermore, SRSF1 is required for antigen-specific IFN-γ cytokine responses in a viral infection challenge in mice. Transcriptomics analyses of Srsf1-deficient T cells reveal that SRSF1 controls proliferation, MAP kinase signaling and IFN signaling pathways. Mechanistically, SRSF1 controls the expression and activity of the Mnk2/p38-MAPK axis at the molecular level. Our findings reveal previously unrecognized roles for SRSF1 in the physiology and function of cytotoxic CD8 T lymphocytes and a potential molecular mechanism in viral immunopathogenesis.