Mapping the splicing landscape of the human immune system.

Most human genes code for more than one transcript. Different ratios of transcripts of the same gene can be found in different cell types or states, indicating differential use of transcription start sites or differential splicing. Such differential transcript use (DTUs) events provide an additional layer of regulation and protein diversity. With the exceptions of PTPRC and CIITA, there are very few reported cases of DTU events in the immune system. To rigorously map DTUs between different human immune cell types, we leveraged four publicly available RNA sequencing datasets. We identified 282 DTU events between five human healthy immune cell types that appear in at least two datasets. The patterns of the DTU events were mostly cell-type-specific or lineage-specific, in the context of the five cell types tested. DTUs correlated with the expression pattern of potential regulators, namely, splicing factors and transcription factors. Of the several immune related conditions studied, only sepsis affected the splicing of more than a few genes and only in innate immune cells. Taken together, we map the DTUs landscape in human peripheral blood immune cell types, and present hundreds of genes whose transcript use changes between cell types or upon activation.

The interweaved signatures of common-gamma-chain cytokines across immunologic lineages.

"γc" cytokines are a family whose receptors share a "common-gamma-chain" signaling moiety, and play central roles in differentiation, homeostasis, and communications of all immunocyte lineages. As a resource to better understand their range and specificity of action, we profiled by RNAseq the immediate-early responses to the main γc cytokines across all immunocyte lineages. The results reveal an unprecedented landscape: broader, with extensive overlap between cytokines (one cytokine doing in one cell what another does elsewhere) and essentially no effects unique to any one cytokine. Responses include a major downregulation component and a broad Myc-controlled resetting of biosynthetic and metabolic pathways. Various mechanisms appear involved: fast transcriptional activation, chromatin remodeling, and mRNA destabilization. Other surprises were uncovered: IL2 effects in mast cells, shifts between follicular and marginal zone B cells, paradoxical and cell-specific cross-talk between interferon and γc signatures, or an NKT-like program induced by IL21 in CD8+ T cells.

The transcriptional landscape of the giant freshwater prawn: Embryonic development and early sexual differentiation mechanisms.

The giant freshwater prawn pjMacrobrachium rosenbergii is one of the best studied species in aquaculture. However, the transcriptional changes associated with embryonic development and the sexual differentiation mechanism of M. rosenbergii remain to be elucidated. To characterize the embryonic development of this prawn and to determine whether differential expression and differential splicing play roles in the early sexual differentiation of M. rosenbergii, we profiled five developmental days of male and female embryos by RNA sequencing. We identified modules of co-expressed genes representing waves of transcription that correspond to physiological processes in early embryonic development (such as the maternal-to-zygotic transition) up to preparation for life outside the egg (development of muscles, cuticle etc.). Additionally, we found that hundreds of genes are differentially expressed between sexes, most of them uncharacterized, suggesting that the sex differentiation mechanism of M. rosenbergii might contain clade-specific elements. The resulting first-of-a-kind transcriptional map of embryonic development of male and female M. rosenbergii will guide future studies to reveal the roles of specific genes and splicing isoforms in the embryonic development and sexual differentiation process of M. rosenbergii.

Different sea urchin RAG-like genes were domesticated to carry out different functions.

The closely linked recombination activating genes (RAG1 and RAG2) in vertebrates encode the core of the RAG recombinase that mediates the V(D)J recombination of the immunoglobulin and T-cell receptor genes. RAG1 and RAG2 homologues (RAG1L and RAG2L) are present in multiple invertebrate phyla, including mollusks, nemerteans, cnidarians, and sea urchins. However, the function of the invertebrates' RAGL proteins is yet unknown. The sea urchins contain multiple RAGL genes that presumably originated in a common ancestral transposon. In this study, we demonstrated that two different RAG1L genes in the sea urchin Paracentrutus lividus (PlRAG1La and PlRAG1Lb) lost their mobility and, along with PlRAG2L, were fully domesticated to carry out different functions. We found that the examined echinoid RAGL homologues have distinct expression profiles in early developmental stages and in adult tissues. Moreover, the predicted structure of the proteins suggests that while PlRAG1La could maintain its endonuclease activity and create a heterotetramer with PlRAG2L, the PlRAG1Lb adopted a different function that does not include an interaction with DNA nor a collaboration with PlRAG2L. By characterizing the different RAG homologues in the echinoid lineage, we hope to increase the knowledge about the evolution of these genes and shed light on their domestication processes.

Post-sepsis immunosuppression depends on NKT cell regulation of mTOR/IFN-γ in NK cells.

As treatment of the early, inflammatory phase of sepsis improves, post-sepsis immunosuppression and secondary infection have increased in importance. How early inflammation drives immunosuppression remains unclear. Although IFN-γ typically helps microbial clearance, we found that increased plasma IFN-γ in early clinical sepsis was associated with the later development of secondary Candida infection. Consistent with this observation, we found that exogenous IFN-γ suppressed macrophage phagocytosis of zymosan in vivo, and antibody blockade of IFN-γ after endotoxemia improved survival of secondary candidemia. Transcriptomic analysis of innate lymphocytes during endotoxemia suggested that NKT cells drove IFN-γ production by NK cells via mTORC1. Activation of invariant NKT (iNKT) cells with glycolipid antigen drove immunosuppression. Deletion of iNKT cells in Cd1d-/- mice or inhibition of mTOR by rapamycin reduced immunosuppression and susceptibility to secondary Candida infection. Thus, although rapamycin is typically an immunosuppressive medication, in the context of sepsis, rapamycin has the opposite effect. These results implicated an NKT cell/mTOR/IFN-γ axis in immunosuppression following endotoxemia or sepsis. In summary, in vivo iNKT cells activated mTORC1 in NK cells to produce IFN-γ, which worsened macrophage phagocytosis, clearance of secondary Candida infection, and mortality.

ImmGen report: sexual dimorphism in the immune system transcriptome.

Sexual dimorphism in the mammalian immune system is manifested as more frequent and severe infectious diseases in males and, on the other hand, higher rates of autoimmune disease in females, yet insights underlying those differences are still lacking. Here we characterize sex differences in the immune system by RNA and ATAC sequence profiling of untreated and interferon-induced immune cell types in male and female mice. We detect very few differentially expressed genes between male and female immune cells except in macrophages from three different tissues. Accordingly, very few genomic regions display differences in accessibility between sexes. Transcriptional sexual dimorphism in macrophages is mediated by genes of innate immune pathways, and increases after interferon stimulation. Thus, the stronger immune response of females may be due to more activated innate immune pathways prior to pathogen invasion.

JingleBells: A Repository of Immune-Related Single-Cell RNA-Sequencing Datasets.

Recent advances in single-cell RNA-sequencing (scRNA-seq) technology increase the understanding of immune differentiation and activation processes, as well as the heterogeneity of immune cell types. Although the number of available immune-related scRNA-seq datasets increases rapidly, their large size and various formats render them hard for the wider immunology community to use, and read-level data are practically inaccessible to the non-computational immunologist. To facilitate datasets reuse, we created the JingleBells repository for immune-related scRNA-seq datasets ready for analysis and visualization of reads at the single-cell level (http://jinglebells.bgu.ac.il/). To this end, we collected the raw data of publicly available immune-related scRNA-seq datasets, aligned the reads to the relevant genome, and saved aligned reads in a uniform format, annotated for cell of origin. We also added scripts and a step-by-step tutorial for visualizing each dataset at the single-cell level, through the commonly used Integrated Genome Viewer (www.broadinstitute.org/igv/). The uniform scRNA-seq format used in JingleBells can facilitate reuse of scRNA-seq data by computational biologists. It also enables immunologists who are interested in a specific gene to visualize the reads aligned to this gene to estimate cell-specific preferences for splicing, mutation load, or alleles. Thus JingleBells is a resource that will extend the usefulness of scRNA-seq datasets outside the programming aficionado realm.

Mapping Whole-Transcriptome Splicing in Mouse Hematopoietic Stem Cells.

Hematopoietic stem cells (HSCs) are rare cells that generate all the various types of blood and immune cells. High-quality transcriptome data have enabled the identification of significant genes for HSCs. However, most genes are expressed in various forms by alternative splicing (AS), extending transcriptome complexity. Here, we delineate AS to determine which isoforms are expressed in mouse HSCs. Our analysis of microarray and RNA-sequencing data includes differential expression of splicing factors that may regulate AS, and a complete map of splicing isoforms. Multiple types of isoforms for known HSC genes and unannotated splicing that may alter gene function are presented. Transcriptome-wide identification of genes and their respective isoforms in mouse HSCs will open another dimension for adult stem cells.

Tight coevolution of proliferating cell nuclear antigen (PCNA)-partner interaction networks in fungi leads to interspecies network incompatibility.

The structure and connectivity of protein-protein interaction (PPI) networks are maintained throughout evolution by coordinated changes (coevolution) of network proteins. Despite extensive research, relatively little is known regarding the molecular basis and functional implications of the coevolution of PPI networks. Here, we used proliferating cell nuclear antigen, a hub protein that mediates DNA replication and repair in eukaryotes, as a model system to study the coevolution of PPI networks in fungi. Using a combined bioinformatics and experimental approach, we discovered that PCNA-partner interactions tightly coevolved in fungal species, leading to specific modes of recognition. We found that fungal proliferating cell nuclear antigen-partner interaction networks diverged into two distinct groups as a result of such coevolution and that hybrid networks of these groups are functionally noncompatible in Saccharomyces cerevisiae. Our results indicate that the coevolution of PPI networks can form functional barriers between fungal species, and thus can promote and fix speciation.

Comparative cross-species alternative splicing in plants.

Alternative splicing (AS) can add significantly to genome complexity. Plants are thought to exhibit less AS than animals. An algorithm, based on expressed sequence tag (EST) pairs gapped alignment, was developed that takes advantage of the relatively small intron and exon size in plants and directly compares pairs of ESTs to search for AS. EST pairs gapped alignment was first evaluated in Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and tomato (Solanum lycopersicum) for which annotated genome sequence is available and was shown to accurately predict splicing events. The method was then applied to 11 plant species that include 17 cultivars for which enough ESTs are available. The results show a large, 3.7-fold difference in AS rates between plant species with Arabidopsis and rice in the lower range and lettuce (Lactuca sativa) and sorghum (Sorghum bicolor) in the upper range. Hence, compared to higher animals, plants show a much greater degree of variety in their AS rates and in some plant species the rates of animal and plant AS are comparable although the distribution of AS types may differ. In eudicots but not monocots, a correlation between genome size and AS rates was detected, implying that in eudicots the mechanisms that lead to larger genomes are a driving force for the evolution of AS.

Whole-genome microarray in Arabidopsis facilitates global analysis of retained introns.

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism that can increase protein diversity and affect mRNA stability. Different types of AS have been observed; these include exon skipping, alternative donor or acceptor site and intron retention. In humans, exon skipping is the most common type while intron retention is rare. In contrast, in Arabidopsis, intron retention is the most prevalent AS type (approximately 40%). Here we show that direct transcript expression analysis using high-density oligonucleotide-based whole-genome microarrays (WGAs) is particularly amenable for assessing global intron retention in Arabidopsis. By applying a novel algorithm retained introns are detected in 8% of the transcripts examined. A sampling of 14 transcripts showed that 86% can be confirmed by RT-PCR. This rate of detection predicts an overall total AS rate of 20% for Arabidopsis compared with 10-22% based on EST/cDNA-based analysis. These findings will facilitate monitoring constitutive and dynamic whole-genome splicing on the next generation WGA slides.

Intron retention is a major phenomenon in alternative splicing in Arabidopsis.

Alternative splicing (AS) combines different transcript splice junctions that result in transcripts with shuffled exons, alternative 5' or 3' splicing sites, retained introns and different transcript termini. In this way, multiple mRNA species and proteins can be created from a single gene expanding the potential informational content of eukaryotic genomes. Search algorithms of AS forms in a variety of Arabidopsis databases showed they contained an unusually high fraction of retained introns (above 30%), compared with 10% that was reported for humans. The preponderance of retained introns (65%) were either part of open reading frames, present in the UTR region or present as the last intron in the transcript, indicating that their occurrence would not participate in non-sense-mediated decay. Interestingly, the functional distribution of the transcripts with retained introns is skewed towards stress and external/internal stimuli-related functions. A sampling of the alternative transcripts with retained introns were confirmed by RT-PCR and were shown to co-purify with polyribosomes, indicating their nuclear export. Thus, retained introns are a prominent feature of AS in Arabidopsis and as such may play a regulatory function.