Event Analysis: Using Transcript Events To Improve Estimates of Abundance in RNA-seq Data

Shifts in isoform usage underlie transcriptional differences in regulatory T cells in type 1 diabetes

Genome-wide association studies have identified numerous loci with allelic associations to Type 1 Diabetes (T1D) risk. Most disease-associated variants are enriched in regulatory sequences active in lymphoid cell types, suggesting that lymphocyte gene expression is altered in T1D. Here we assay gene expression between T1D cases and healthy controls in two autoimmunity-relevant lymphocyte cell types, memory CD4+/CD25+ regulatory T cells (Treg) and memory CD4+/CD25- T cells, using a splicing event-based approach to characterize tissue-specific transcriptomes. Limited differences in isoform usage between T1D cases and controls are observed in memory CD4+/CD25- T-cells. In Tregs, 402 genes demonstrate differences in isoform usage between cases and controls, particularly RNA recognition and splicing factor genes. Many of these genes are regulated by the variable inclusion of exons that can trigger nonsense mediated decay. Our results suggest that dysregulation of gene expression, through shifts in alternative splicing in Tregs, contributes to T1D pathophysiology.

UBASH3A Interacts with PTPN22 to Regulate IL2 Expression and Risk for Type 1 Diabetes

UBASH3A is a negative regulator of T cell activation and IL-2 production and plays key roles in autoimmunity. Although previous studies revealed the individual effects of UBASH3A on risk for type 1 diabetes (T1D; a common autoimmune disease), the relationship of UBASH3A with other T1D risk factors remains largely unknown. Given that another well-known T1D risk factor, PTPN22, also inhibits T cell activation and IL-2 production, we investigated the relationship between UBASH3A and PTPN22. We found that UBASH3A, via its Src homology 3 (SH3) domain, physically interacts with PTPN22 in T cells, and that this interaction is not altered by the T1D risk coding variant rs2476601 in PTPN22. Furthermore, our analysis of RNA-seq data from T1D cases showed that the amounts of UBASH3A and PTPN22 transcripts exert a cooperative effect on IL2 expression in human primary CD8⁺ T cells. Finally, our genetic association analyses revealed that two independent T1D risk variants, rs11203203 in UBASH3A and rs2476601 in PTPN22, interact statistically, jointly affecting risk for T1D. In summary, our study reveals novel interactions, both biochemical and statistical, between two independent T1D risk loci, and suggests how these interactions may affect T cell function and increase risk for T1D.

Sex-Biased Expression Is Associated With Chromatin State in Drosophila melanogaster and Drosophila simulans

In Drosophila melanogaster and D. simulans head tissue, 60% of orthologous genes show evidence of sex-biased expression in at least one species. Of these, ∼39% (2,192) are conserved in direction. We hypothesize enrichment of open chromatin in the sex where we see expression bias and closed chromatin in the opposite sex. Male-biased orthologs are significantly enriched for H3K4me3 marks in males of both species (∼89% of male-biased orthologs vs. ∼76% of unbiased orthologs). Similarly, female-biased orthologs are significantly enriched for H3K4me3 marks in females of both species (∼90% of female-biased orthologs vs. ∼73% of unbiased orthologs). The sex-bias ratio in female-biased orthologs was similar in magnitude between the two species, regardless of the closed chromatin (H3K27me2me3) marks in males. However, in male-biased orthologs, the presence of H3K27me2me3 in both species significantly reduced the correlation between D. melanogaster sex-bias ratio and the D. simulans sex-bias ratio. Male-biased orthologs are enriched for evidence of positive selection in the D. melanogaster group. There are more male-biased genes than female-biased genes in both species. For orthologs with gains/losses of sex-bias between the two species, there is an excess of male-bias compared to female-bias, but there is no consistent pattern in the relationship between H3K4me3 or H3K27me2me3 chromatin marks and expression. These data suggest chromatin state is a component of the maintenance of sex-biased expression and divergence of sex-bias between species is reflected in the complexity of the chromatin status.

Variation in leaf transcriptome responses to elevated ozone corresponds with physiological sensitivity to ozone across maize inbred lines

We examine the impact of sustained elevated ozone concentration on the leaf transcriptome of 5 diverse maize inbred genotypes, which vary in physiological sensitivity to ozone (B73, Mo17, Hp301, C123, NC338), using long reads to assemble transcripts and short reads to quantify expression of these transcripts. More than 99% of the long reads, 99% of the assembled transcripts, and 97% of the short reads map to both B73 and Mo17 reference genomes. Approximately 95% of the genes with assembled transcripts belong to known B73-Mo17 syntenic loci and 94% of genes with assembled transcripts are present in all temperate lines in the NAM pan-genome. While there is limited evidence for alternative splicing in response to ozone stress, there is a difference in the magnitude of differential expression among the 5 genotypes. The transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes) was modest, while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. There is the potential for tandem duplication in 30% of genes with assembled transcripts, but there is no obvious association between potential tandem duplication and differential expression. Genes with a common response across the 5 genotypes (83 genes) were associated with photosynthesis, in particular photosystem I. The functional annotation of genes not differentially expressed in B73 but responsive in the other 4 genotypes (789) identifies reactive oxygen species. This suggests that B73 has a different response to long term ozone exposure than the other 4 genotypes. The relative magnitude of the genotypic response to ozone, and the enrichment analyses are consistent regardless of whether aligning short reads to: long read assembled transcripts; the B73 reference; the Mo17 reference. We find that prolonged ozone exposure directly impacts the photosynthetic machinery of the leaf.

Genetic Control of Splicing at SIRPG Modulates Risk of Type 1 Diabetes

Signal regulatory protein SIRPγ (CD172G) is expressed on the surface of lymphocytes, where it acts by engaging its ligand, CD47. SIRPG, which encodes SIRPγ, contains a nonsynonymous coding variant, rs6043409, which is significantly associated with risk for type 1 diabetes. SIRPG produces multiple transcript isoforms via alternative splicing, all encoding potentially functional proteins. We show that rs6043409 alters a predicted exonic splicing enhancer, resulting in significant shifts in the distribution of SIRPG transcript isoforms. All of these transcript isoforms produced protein upon transient expression in vitro. However, CRISPR/Cas9 targeting of one of the alternatively spliced exons in SIRPG eliminated all SIRPγ expression in Jurkat T cells. These targeted cells formed fewer cell-cell conjugates with each other than with wild-type Jurkat cells, expressed reduced levels of genes associated with CD47 signaling, and had significantly increased levels of cell-surface CD47. In primary CD4+ and CD8+ T cells, cell-surface SIRPγ levels in response to anti-CD3 stimulation varied quantitatively by rs6043409 genotype. Our results suggest that SIRPG is the most likely causative gene for type 1 diabetes risk in the 20p13 region and highlight the role of alternative splicing in lymphocytes in mediating the genetic risk for autoimmunity.

tappAS: a comprehensive computational framework for the analysis of the functional impact of differential splicing

Recent advances in long-read sequencing solve inaccuracies in alternative transcript identification of full-length transcripts in short-read RNA-Seq data, which encourages the development of methods for isoform-centered functional analysis. Here, we present tappAS, the first framework to enable a comprehensive Functional Iso-Transcriptomics (FIT) analysis, which is effective at revealing the functional impact of context-specific post-transcriptional regulation. tappAS uses isoform-resolved annotation of coding and non-coding functional domains, motifs, and sites, in combination with novel analysis methods to interrogate different aspects of the functional readout of transcript variants and isoform regulation. tappAS software and documentation are available at https://app.tappas.org.