Alternative splicing (AS) allows increased diversity and orthogonal regulation of the transcriptional products of mammalian genomes. To assess the distribution and variation of alternative splicing across cell lineages of the immune system, we comprehensively analyzed RNA sequencing and microarray data generated by the Immunological Genome Project Consortium. AS is pervasive: 60% of genes showed frequent AS isoforms in T or B lymphocytes, with 7,599 previously unreported isoforms. Distinct cell specificity was observed, with differential exon skipping in 5% of genes otherwise coexpressed in both B and T cells. The distribution of isoforms was mostly all or none, suggesting on/off switching as a frequent mode of AS regulation in lymphocytes. From the identification of differential exon use in the microarray data, clustering of exon inclusion/exclusion patterns across all Immunological Genome Project cell types showed that ∼70% of AS exons are distributed along a common pattern linked to lineage differentiation and cell cycling. Other AS events distinguished myeloid from lymphoid cells or affected only a small set of exons without clear lineage specificity (e.g., Ptprc). Computational analysis predicted specific associations between AS exons and splicing regulators, which were verified by detection of the hnRPLL/Ptprc connection.

To determine the breadth and underpinning of changes in immunocyte gene expression due to genetic variation in mice, we performed, as part of the Immunological Genome Project, gene expression profiling for CD4(+) T cells and neutrophils purified from 39 inbred strains of the Mouse Phenome Database. Considering both cell types, a large number of transcripts showed significant variation across the inbred strains, with 22% of the transcriptome varying by 2-fold or more. These included 119 loci with apparent complete loss of function, where the corresponding transcript was not expressed in some of the strains, representing a useful resource of "natural knockouts." We identified 1222 cis-expression quantitative trait loci (cis-eQTL) that control some of this variation. Most (60%) cis-eQTLs were shared between T cells and neutrophils, but a significant portion uniquely impacted one of the cell types, suggesting cell type-specific regulatory mechanisms. Using a conditional regression algorithm, we predicted regulatory interactions between transcription factors and potential targets, and we demonstrated that these predictions overlap with regulatory interactions inferred from transcriptional changes during immunocyte differentiation. Finally, comparison of these and parallel data from CD4(+) T cells of healthy humans demonstrated intriguing similarities in variability of a gene's expression: the most variable genes tended to be the same in both species, and there was an overlap in genes subject to strong cis-acting genetic variants. We speculate that this "conservation of variation" reflects a differential constraint on intraspecies variation in expression levels of different genes, either through lower pressure for some genes, or by favoring variability for others.