The nomination of candidate genes underlying complex traits is often focused on genetic variations that alter mRNA abundance or result in non-conservative changes in amino acids. Although inconspicuous in complex trait analysis, genetic variants that affect splicing or RNA editing can also generate proteomic diversity and impact genetic traits. Indeed, it is known that splicing and RNA editing modulate several traits in humans and model organisms. Using high-throughput RNA sequencing (RNA-seq) analysis, it is now possible to integrate the genetics of transcript abundance, alternative splicing (AS) and editing with the analysis of complex traits. We recently demonstrated that both AS and mRNA editing are modulated by genetic and environmental factors, and potentially engender phenotypic diversity in a genetically segregating mouse population. Therefore, the analysis of splicing and RNA editing can expand not only the regulatory landscape of transcriptome and proteome complexity, but also the repertoire of candidate genes for complex traits.

Alternative splicing and mRNA editing are known to contribute to transcriptome diversity. Although alternative splicing is pervasive and contributes to a variety of pathologies, including cancer, the genetic context for individual differences in isoform usage is still evolving. Similarly, although mRNA editing is ubiquitous and associated with important biological processes such as intracellular viral replication and cancer development, individual variations in mRNA editing and the genetic transmissibility of mRNA editing are equivocal. Here, we have used linkage analysis to show that both mRNA editing and alternative splicing are regulated by the macrophage genetic background and environmental cues. We show that distinct loci, potentially harboring variable splice factors, regulate the splicing of multiple transcripts. Additionally, we show that individual genetic variability at the Apobec1 locus results in differential rates of C-to-U(T) editing in murine macrophages; with mouse strains expressing mostly a truncated alternative transcript isoform of Apobec1 exhibiting lower rates of editing. As a proof of concept, we have used linkage analysis to identify 36 high-confidence novel edited sites. These results provide a novel and complementary method that can be used to identify C-to-U editing sites in individuals segregating at specific loci and show that, beyond DNA sequence and structural changes, differential isoform usage and mRNA editing can contribute to intra-species genomic and phenotypic diversity.

Background

Accurate gene model predictions and annotation of alternative splicing events are imperative for genomic studies in organisms that contain genes with multiple exons. Currently most gene models for the intracellular parasite, Toxoplasma gondii, are based on computer model predictions without cDNA sequence verification. Additionally, the nature and extent of alternative splicing in Toxoplasma gondii is unknown. In this study, we used de novo transcript assembly and the published type II (ME49) genomic sequence to quantify the extent of alternative splicing in Toxoplasma and to improve the current Toxoplasma gene annotations.

Results

We used high-throughput RNA-sequencing data to assemble full-length transcripts, independently of a reference genome, followed by gene annotation based on the ME49 genome. We assembled 13,533 transcripts overlapping with known ME49 genes in ToxoDB and then used this set to; a) improve the annotation in the untranslated regions of ToxoDB genes, b) identify novel exons within protein-coding ToxoDB genes, and c) report on 50 previously unidentified alternatively spliced transcripts. Additionally, we assembled a set of 2,930 transcripts not overlapping with any known ME49 genes in ToxoDB. From this set, we have identified 118 new ME49 genes, 18 novel Toxoplasma genes, and putative non-coding RNAs.

Conclusion

RNA-seq data and de novo transcript assembly provide a robust way to update incompletely annotated genomes, like the Toxoplasma genome. We have used RNA-seq to improve the annotation of several Toxoplasma genes, identify alternatively spliced genes, novel genes, novel exons, and putative non-coding RNAs.

The severe virulence of Toxoplasma gondii in classical laboratory inbred mouse strains contradicts the hypothesis that house mice (Mus musculus) are the most important intermediate hosts for its transmission and evolution because death of the mouse before parasite transmission equals death of the parasite. However, the classical laboratory inbred mouse strains (Mus musculus domesticus), commonly used to test Toxoplasma strain differences in virulence, do not capture the genetic diversity within Mus musculus. Thus, it is possible that Toxoplasma strains that are severely virulent in laboratory inbred mice are avirulent in some other mouse sub-species. Here, we present insight into the responses of individual mouse strains, representing strains of the genetically divergent Mus musculus musculus, Mus musculus castaneus and Mus musculus domesticus, to infection with individual clonal and atypical Toxoplasma strains. We observed that, unlike M. m. domesticus, M. m. musculus and M. m. castaneus are resistant to the clonal Toxoplasma strains. For M. m. musculus, we show that this is due to a locus on chromosome 11 that includes the genes that encode the interferon gamma (IFNG)-inducible immunity-related GTPases (Irgs) that can kill the parasite by localising and subsequently vesiculating the parasitophorous vacuole membrane. However, despite the localization of known effector Irgs to the Toxoplasma parasitophorous vacuole membrane, we observed that some atypical Toxoplasma strains are virulent in all the mouse strains tested. The virulence of these atypical strains in M. m. musculus could not be attributed to individual rhoptry protein 5 (ROP5) alleles, a secreted parasite pseudokinase that antagonises the canonical effector Irgs and is indispensable for parasite virulence in laboratory inbred mice (M. m. domesticus). We conclude that murine resistance to Toxoplasma is modulated by complex interactions between host and parasite genotypes and may be independent of known effector Irgs on murine chromosome 11.

Macrophages display flexible activation states that range between pro-inflammatory (classical activation) and anti-inflammatory (alternative activation). These macrophage polarization states contribute to a variety of organismal phenotypes such as tissue remodeling and susceptibility to infectious and inflammatory diseases. Several macrophage- or immune-related genes have been shown to modulate infectious and inflammatory disease pathogenesis. However, the potential role that differences in macrophage activation phenotypes play in modulating differences in susceptibility to infectious and inflammatory disease is just emerging. We integrated transcriptional profiling and linkage analyses to determine the genetic basis for the differential murine macrophage response to inflammatory stimuli and to infection with the obligate intracellular parasite Toxoplasma gondii. We show that specific transcriptional programs, defined by distinct genomic loci, modulate macrophage activation phenotypes. In addition, we show that the difference between AJ and C57BL/6J macrophages in controlling Toxoplasma growth after stimulation with interferon gamma and tumor necrosis factor alpha mapped to chromosome 3, proximal to the Guanylate binding protein (Gbp) locus that is known to modulate the murine macrophage response to Toxoplasma. Using an shRNA-knockdown strategy, we show that the transcript levels of an RNA helicase, Ddx1, regulates strain differences in the amount of nitric oxide produced by macrophage after stimulation with interferon gamma and tumor necrosis factor. Our results provide a template for discovering candidate genes that modulate macrophage-mediated complex traits.

Unlabelled

Differences among individuals in susceptibility to infectious diseases can be modulated by host genetics. Much of the research in this field has aimed to identify loci within the host genome that are associated with these differences. In mice, A/J (AJ) and C57BL/6J (B6) mice show differential susceptibilities to various pathogens, including the intracellular pathogen Francisella tularensis. Because macrophages are the main initial target during F. tularensis infection, we explored early interactions of macrophages from these two mouse strains with F. tularensis as well as the genetic factors underlying these interactions. Our results indicate that bacterial interactions with bone marrow-derived macrophages (BMDMs) during early stages of infection are different in the AJ and B6 strains. During these early stages, bacteria are more numerous in B6 than in AJ macrophages and display differences in trafficking and early transcriptional response within these macrophages. To determine the genetic basis for these differences, we infected BMDMs isolated from recombinant inbred (RI) mice derived from reciprocal crosses between AJ and B6, and we followed early bacterial counts within these macrophages. Quantitative trait locus (QTL) analysis revealed a locus on chromosome 19 that is associated with early differences in bacterial counts in AJ versus B6 macrophages. QTL analysis of published data that measured the differential susceptibilities of the same RI mice to an in vivo challenge with F. tularensis confirmed the F. tularensis susceptibility QTL on chromosome 19. Overall, our results show that early interactions of macrophages with F. tularensis are dependent on the macrophage genetic background.

Importance

Francisella tularensis is a highly pathogenic bacterium with a very low infectious dose in humans. Some mechanisms of bacterial virulence have been elucidated, but the host genetic factors that contribute to host resistance or susceptibility are largely unknown. In this work, we have undertaken a genetic approach to assess what these factors are in mice. Analyzing early interactions of macrophages with the bacteria as well as data on overall susceptibility to infection revealed a locus on chromosome 19 that is associated with both phenotypes. In addition, our work revealed differences in the early macrophage response between macrophages with different genetic backgrounds. Overall, this work suggests some intriguing links between in vitro and in vivo infection models and should aid in further elucidating the genetic circuits behind the host response to Francisella tularensis infection.

Toxoplasma gondii is an intracellular parasite that infects a wide range of warm-blooded species. Rats vary in their susceptibility to this parasite. The Toxo1 locus conferring Toxoplasma resistance in rats was previously mapped to a region of chromosome 10 containing Nlrp1. This gene encodes an inflammasome sensor controlling macrophage sensitivity to anthrax lethal toxin (LT) induced rapid cell death (pyroptosis). We show here that rat strain differences in Toxoplasma infected macrophage sensitivity to pyroptosis, IL-1β/IL-18 processing, and inhibition of parasite proliferation are perfectly correlated with NLRP1 sequence, while inversely correlated with sensitivity to anthrax LT-induced cell death. Using recombinant inbred rats, SNP analyses and whole transcriptome gene expression studies, we narrowed the candidate genes for control of Toxoplasma-mediated rat macrophage pyroptosis to four genes, one of which was Nlrp1. Knockdown of Nlrp1 in pyroptosis-sensitive macrophages resulted in higher parasite replication and protection from cell death. Reciprocally, overexpression of the NLRP1 variant from Toxoplasma-sensitive macrophages in pyroptosis-resistant cells led to sensitization of these resistant macrophages. Our findings reveal Toxoplasma as a novel activator of the NLRP1 inflammasome in rat macrophages.

Macrophages play an essential role in the early immune response against Toxoplasma and are the cell type preferentially infected by the parasite in vivo. Interferon gamma (IFNγ) elicits a variety of anti-Toxoplasma activities in macrophages. Using a genome-wide CRISPR screen we identify 353 Toxoplasma genes that determine parasite fitness in naїve or IFNγ-activated murine macrophages, seven of which are further confirmed. We show that one of these genes encodes dense granule protein GRA45, which has a chaperone-like domain, is critical for correct localization of GRAs into the PVM and secretion of GRA effectors into the host cytoplasm. Parasites lacking GRA45 are more susceptible to IFNγ-mediated growth inhibition and have reduced virulence in mice. Together, we identify and characterize an important chaperone-like GRA in Toxoplasma and provide a resource for the community to further explore the function of Toxoplasma genes that determine fitness in IFNγ-activated macrophages.

Upon invasion of Lewis rat macrophages, Toxoplasma rapidly induces programmed cell death (pyroptosis), which prevents Toxoplasma replication, possibly explaining the resistance of the Lewis rat to Toxoplasma Using a chemical mutagenesis screen, we identified Toxoplasma mutants that no longer induced pyroptosis. Whole-genome sequencing led to the identification of three Toxoplasma parasitophorous vacuole-localized dense granule proteins, GRA35, GRA42, and GRA43, that are individually required for induction of Lewis rat macrophage pyroptosis. Macrophage infection with Δgra35, Δgra42, and Δgra43 parasites led to greatly reduced cell death rates and enhanced parasite replication. Lewis rat macrophages infected with parasites containing a single, double, or triple deletion of these GRAs showed similar levels of cell viability, suggesting that the three GRAs function in the same pathway. Deletion of GRA42 or GRA43 resulted in GRA35 (and other GRAs) being retained inside the parasitophorous vacuole instead of being localized to the parasitophorous vacuole membrane. Despite having greatly enhanced replication in Lewis rat macrophages in vitro, Δgra35, Δgra42, and Δgra43 parasites did not establish a chronic infection in Lewis rats. Toxoplasma did not induce F344 rat macrophage pyroptosis, but F344 rats infected with Δgra35, Δgra42, and Δgra43 parasites had reduced cyst numbers. Thus, these GRAs determined parasite in vivo fitness in F344 rats. Overall, our data suggest that these three Toxoplasma dense granule proteins play a critical role in establishing a chronic infection in vivo, independently of their role in mediating macrophage pyroptosis, likely due to their importance in regulating protein localization to the parasitophorous vacuole membrane.IMPORTANCE Inflammasomes are major components of the innate immune system and are responsible for detecting various microbial and environmental danger signals. Upon invasion of Lewis rat macrophages, the parasite rapidly activates the NLRP1 inflammasome, resulting in pyroptosis and elimination of the parasite's replication niche. The work reported here revealed that Toxoplasma GRA35, GRA42, and GRA43 are required for induction of Lewis rat macrophage pyroptosis. GRA42 and GRA43 mediate the correct localization of other GRAs, including GRA35, to the parasitophorous vacuole membrane. These three GRAs were also found to be important for parasite in vivo fitness in a Toxoplasma-susceptible rat strain, independently of their role in NLRP1 inflammasome activation, suggesting that they perform other important functions. Thus, this study identified three GRAs that mediate the induction of Lewis rat macrophage pyroptosis and are required for pathogenesis of the parasite.

Most isolates of Toxoplasma from Europe and North America fall into one of three genetically distinct clonal lineages, the type I, II and III lineages. However, in South America these strains are rarely isolated and instead a great variety of other strains are found. T. gondii strains differ widely in a number of phenotypes in mice, such as virulence, persistence, oral infectivity, migratory capacity, induction of cytokine expression and modulation of host gene expression. The outcome of toxoplasmosis in patients is also variable and we hypothesize that, besides host and environmental factors, the genotype of the parasite strain plays a major role. The molecular basis for these differences in pathogenesis, especially in strains other than the clonal lineages, remains largely unexplored. Macrophages play an essential role in the early immune response against T. gondii and are also the cell type preferentially infected in vivo. To determine if non-canonical Toxoplasma strains have unique interactions with the host cell, we infected murine macrophages with 29 different Toxoplasma strains, representing global diversity, and used RNA-sequencing to determine host and parasite transcriptomes. We identified large differences between strains in the expression level of known parasite effectors and large chromosomal structural variation in some strains. We also identified novel strain-specifically regulated host pathways, including the regulation of the type I interferon response by some atypical strains. IFNβ production by infected cells was associated with parasite killing, independent of interferon gamma activation, and dependent on endosomal Toll-like receptors in macrophages and the cytoplasmic receptor retinoic acid-inducible gene 1 (RIG-I) in fibroblasts.