Understanding host response to infectious salmon anaemia virus in an Atlantic salmon cell line using single-cell RNA sequencing.

BACKGROUND: Infectious Salmon Anaemia Virus (ISAV) is an Orthomixovirus that represents a large problem for salmonid aquaculture worldwide. Current prevention and treatment methods are only partially effective. Genetic selection and genome engineering have the potential to develop ISAV resistant salmon stocks. Both strategies can benefit from an improved understanding of the genomic regulation of ISAV pathogenesis. Here, we used single-cell RNA sequencing of an Atlantic salmon cell line to provide the first high dimensional insight into the transcriptional landscape that underpins host-virus interaction during early ISAV infection. RESULTS: Salmon head kidney (SHK-1) cells were single-cell RNA sequenced at 24, 48 and 96 h post-ISAV challenge. At 24 h post infection, cells showed expression signatures consistent with viral entry, with genes such as PI3K, FAK or JNK being upregulated relative to uninfected cells. At 48 and 96 h, infected cells showed a clear anti-viral response, characterised by the expression of IFNA2 or IRF2. Uninfected bystander cells at 48 and 96 h also showed clear transcriptional differences, potentially suggesting paracrine signalling from infected cells. These bystander cells expressed pathways such as mRNA sensing, RNA degradation, ubiquitination or proteasome; and up-regulation of mitochondrial ribosome genes also seemed to play a role in the host response to the infection. Correlation between viral and host genes revealed novel genes potentially key for this fish-virus interaction. CONCLUSIONS: This study has increased our understanding of the cellular response of Atlantic salmon during ISAV infection and revealed host-virus interactions at the cellular level. Our results highlight various potential key genes in this host-virus interaction, which can be manipulated in future functional studies to increase the resistance of Atlantic salmon to ISAV.

Stem cell-derived porcine macrophages as a new platform for studying host-pathogen interactions.

BACKGROUND: Infectious diseases of farmed and wild animals pose a recurrent threat to food security and human health. The macrophage, a key component of the innate immune system, is the first line of defence against many infectious agents and plays a major role in shaping the adaptive immune response. However, this phagocyte is a target and host for many pathogens. Understanding the molecular basis of interactions between macrophages and pathogens is therefore crucial for the development of effective strategies to combat important infectious diseases. RESULTS: We explored how porcine pluripotent stem cells (PSCs) can provide a limitless in vitro supply of genetically and experimentally tractable macrophages. Porcine PSC-derived macrophages (PSCdMs) exhibited molecular and functional characteristics of ex vivo primary macrophages and were productively infected by pig pathogens, including porcine reproductive and respiratory syndrome virus (PRRSV) and African swine fever virus (ASFV), two of the most economically important and devastating viruses in pig farming. Moreover, porcine PSCdMs were readily amenable to genetic modification by CRISPR/Cas9 gene editing applied either in parental stem cells or directly in the macrophages by lentiviral vector transduction. CONCLUSIONS: We show that porcine PSCdMs exhibit key macrophage characteristics, including infection by a range of commercially relevant pig pathogens. In addition, genetic engineering of PSCs and PSCdMs affords new opportunities for functional analysis of macrophage biology in an important livestock species. PSCs and differentiated derivatives should therefore represent a useful and ethical experimental platform to investigate the genetic and molecular basis of host-pathogen interactions in pigs, and also have wider applications in livestock.

Using machine learning to detect the differential usage of novel gene isoforms.

BACKGROUND: Differential isoform usage is an important driver of inter-individual phenotypic diversity and is linked to various diseases and traits. However, accurately detecting the differential usage of different gene transcripts between groups can be difficult, in particular in less well annotated genomes where the spectrum of transcript isoforms is largely unknown. RESULTS: We investigated whether machine learning approaches can detect differential isoform usage based purely on the distribution of reads across a gene region. We illustrate that gradient boosting and elastic net approaches can successfully identify large numbers of genes showing potential differential isoform usage between Europeans and Africans, that are enriched among relevant biological pathways and significantly overlap those identified by previous approaches. We demonstrate that diversity at the 3' and 5' ends of genes are primary drivers of these differences between populations. CONCLUSION: Machine learning methods can effectively detect differential isoform usage from read fraction data, and can provide novel insights into the biological differences between groups.

Transcriptional and Linkage Analyses Identify Loci that Mediate the Differential Macrophage Response to Inflammatory Stimuli and Infection.

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.

Surgical site infections after craniotomy: a matched health-care cost and length of stay study.

OBJECTIVE: To assess the extra health-care costs and length of stay resulting from surgical site infection (SSI), as well as to identify the most frequent aetiological microorganisms of SSIs among Jordanian craniotomy surgery patients. METHOD: A retrospective, descriptive, correlational and nested 1:1 matched case-control design was used. A computerised list of patients, who underwent surgery between May 2009 and March 2015, was generated in the targeted hospital. A final bill for every selected patient was also determined. Patients were divided equally into two groups: patients with an SSI and patients without an SSI. RESULTS: A total of 64 patients were recruited. The SSI-group had a significant higher mean health-care cost of $7,899.08 (p=0.001) and a longer stay in hospital (mean additional days: 23.17) than the non-SSI group. Furthermore, Acinetobacter baumannii and Staphylococcus aureus were determined as the most predominant causative agents of SSI, at 39.1% and 26.1% of SSI patients, respectively. CONCLUSION: The results of this study can be considered as a baseline for national benchmarking to evaluate the quality of care provided to targeted patients. This study should encourage nurse administrators to adopt protocols and strategies that promote infection control measures, as well as to develop new methods of surveillance on universal precautions adherence. This may limit pathogen contamination in the surgical wound, shorten length of stay and decrease health-care costs.

Comparative ribosome profiling uncovers a dominant role for translational control in Toxoplasma gondii.

BACKGROUND: The lytic cycle of the protozoan parasite Toxoplasma gondii, which involves a brief sojourn in the extracellular space, is characterized by defined transcriptional profiles. For an obligate intracellular parasite that is shielded from the cytosolic host immune factors by a parasitophorous vacuole, the brief entry into the extracellular space is likely to exert enormous stress. Due to its role in cellular stress response, we hypothesize that translational control plays an important role in regulating gene expression in Toxoplasma during the lytic cycle. Unlike transcriptional profiles, insights into genome-wide translational profiles of Toxoplasma gondii are lacking. METHODS: We have performed genome-wide ribosome profiling, coupled with high throughput RNA sequencing, in intracellular and extracellular Toxoplasma gondii parasites to investigate translational control during the lytic cycle. RESULTS: Although differences in transcript abundance were mostly mirrored at the translational level, we observed significant differences in the abundance of ribosome footprints between the two parasite stages. Furthermore, our data suggest that mRNA translation in the parasite is potentially regulated by mRNA secondary structure and upstream open reading frames. CONCLUSION: We show that most of the Toxoplasma genes that are dysregulated during the lytic cycle are translationally regulated.

The genetic basis for individual differences in mRNA splicing and APOBEC1 editing activity in murine macrophages.

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.

Inflammasome sensor NLRP1 controls rat macrophage susceptibility to Toxoplasma gondii.

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.

Incorporating alternative splicing and mRNA editing into the genetic analysis of complex traits.

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.

Transcriptional analysis of murine macrophages infected with different Toxoplasma strains identifies novel regulation of host signaling pathways.

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.

Comparative transcriptional analysis identifies genes associated with the attenuation of Theileria parva infected cells after long-term in vitro culture.

Autologous administration of attenuated Theileria parva-infected cells induces immunity to T. parva in cattle. The mechanism of attenuation, however, is largely unknown. Here, we used RNA sequencing of pathogenic and attenuated T. parva-infected T-cells to elucidate the transcriptional changes underpinning attenuation. We observed differential expression of several host genes, including TRAIL, PD-1, TGF-ß and granzymes that are known to regulate inflammation and proliferation of infected cells. Importantly, many genes linked with the attenuation of the related T. annulata-infected cells were not dysregulated in this study. Furthermore, known T. parva antigens were not dysregulated in attenuated relative to pathogenic cells, indicating that attenuation is not due to enhanced immunogenicity. Overall this study suggests that attenuation is driven by a decrease in proliferation and restoration of the inflammatory profile of T. parva-infected cells. Additionally, it provides a foundation for future mechanistic studies of the attenuation phenotype in Theileria-infected cells.

Toxoplasma gondii rhoptry 16 kinase promotes host resistance to oral infection and intestinal inflammation only in the context of the dense granule protein GRA15.

Toxoplasma gondii transmission between intermediate hosts is dependent on the ingestion of walled cysts formed during the chronic phase of infection. Immediately following consumption, the parasite must ensure survival of the host by preventing adverse inflammatory responses and/or by limiting its own replication. Since the Toxoplasma secreted effectors rhoptry 16 kinase (ROP16) and dense granule 15 (GRA15) activate the JAK-STAT3/6 and NF-κB signaling pathways, respectively, we explored whether a particular combination of these effectors impacted intestinal inflammation and parasite survival in vivo. Here we report that expression of the STAT-activating version of ROP16 in the type II strain (strain II+ROP16I) promotes host resistance to oral infection only in the context of endogenous GRA15 expression. Protection was characterized by a lower intestinal parasite burden and dampened inflammation. Host resistance to the II+ROP16I strain occurred independently of STAT6 and the T cell coinhibitory receptors B7-DC and B7-H1, two receptors that are upregulated by ROP16. In addition, coexpression of ROP16 and GRA15 enhanced parasite susceptibility within tumor necrosis factor alpha/gamma interferon-stimulated macrophages in a STAT3/6-independent manner. Transcriptional profiling of infected STAT3- and STAT6-deficient macrophages and parasitized Peyer's patches from mice orally challenged with strain II+ROP16I suggested that ROP16 activated STAT5 to modulate host gene expression. Consistent with this supposition, the ROP16 kinase induced the sustained phosphorylation and nuclear localization of STAT5 in Toxoplasma-infected cells. In summary, only the combined expression of both GRA15 and ROP16 promoted host resistance to acute oral infection, and Toxoplasma may possibly target the STAT5 signaling pathway to generate protective immunity in the gut.

De novo reconstruction of the Toxoplasma gondii transcriptome improves on the current genome annotation and reveals alternatively spliced transcripts and putative long non-coding RNAs.

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.

Profiling of open chromatin in developing pig (Sus scrofa) muscle to identify regulatory regions.

There is very little information about how the genome is regulated in domestic pigs (Sus scrofa). This lack of knowledge hinders efforts to define and predict the effects of genetic variants in pig breeding programs. To address this knowledge gap, we need to identify regulatory sequences in the pig genome starting with regions of open chromatin. We used the "Improved Protocol for the Assay for Transposase-Accessible Chromatin (Omni-ATAC-Seq)" to identify putative regulatory regions in flash-frozen semitendinosus muscle from 24 male piglets. We collected samples from the smallest-, average-, and largest-sized male piglets from each litter through five developmental time points. Of the 4661 ATAC-Seq peaks identified that represent regions of open chromatin, >50% were within 1 kb of known transcription start sites. Differential read count analysis revealed 377 ATAC-Seq defined genomic regions where chromatin accessibility differed significantly across developmental time points. We found regions of open chromatin associated with downregulation of genes involved in muscle development that were present in small-sized fetal piglets but absent in large-sized fetal piglets at day 90 of gestation. The dataset that we have generated provides a resource for studies of genome regulation in pigs and contributes valuable functional annotation information to filter genetic variants for use in genomic selection in pig breeding programs.

Transcriptional Analyses Identify Genes That Modulate Bovine Macrophage Response to Toxoplasma Infection and Immune Stimulation.

The obligate intracellular parasite, Toxoplasma gondii, is highly prevalent among livestock species. Although cattle are generally resistant to Toxoplasma strains circulating in Europe and North America, the underlying mechanisms are largely unknown. Here, we report that bovine bone marrow-derived macrophage (BMDM) pre-stimulated with interferon gamma (IFNγ) restricts intracellular Toxoplasma growth independently of nitric oxide. While Toxoplasma promoted the expression of genes associated with alternative macrophage activation and lipid metabolism, IFNγ abrogated parasite-induced transcriptional responses and promoted the expression of genes linked to the classical macrophage activation phenotype. Additionally, several chemokines, including CCL22, that are linked to parasite-induced activation of the Wnt/ß-catenin signaling were highly expressed in Toxoplasma-exposed naïve BMDMs. A chemical Wnt/ß-catenin signaling pathway antagonist (IWR-1-endo) significantly reduced intracellular parasite burden in naïve BMDMs, suggesting that Toxoplasma activates this pathway to evade bovine macrophage anti-parasitic responses. Congruently, intracellular burden of a mutant Toxoplasma strain (RHΔASP5) that does not secrete dense granule proteins into the host cell, which is an essential requirement for parasite-induced activation of the Wnt/ß-catenin pathway, was significantly reduced in naïve BMDMs. However, both the Wnt/ß-catenin antagonist and RHASPΔ5 did not abolish parasite burden differences in naïve and IFNγ-stimulated BMDMs. Finally, we observed that parasites infecting IFNγ-stimulated BMDMs largely express genes associated with the slow dividing bradyzoite stage. Overall, this study provides novel insights into bovine macrophage transcriptional response to Toxoplasma. It establishes a foundation for a mechanistic analysis IFNγ-induced bovine anti-Toxoplasma responses and the counteracting Toxoplasma survival strategies.

Single-cell RNA-seq reveals CD16- monocytes as key regulators of human monocyte transcriptional response to Toxoplasma.

Monocytes are among the major myeloid cells that respond to Toxoplasma, a ubiquitous foodborne that infects ≥ 1 billion people worldwide, in human peripheral blood. As such, a molecular understanding of human monocyte-Toxoplasma interactions can expedite the development of novel human toxoplasmosis control strategies. Current molecular studies on monocyte-Toxoplasma interactions are based on average cell or parasite responses across bulk cell populations. Although informative, population-level averages of monocyte responses to Toxoplasma have sometimes produced contradictory results, such as whether CCL2 or IL12 define effective monocyte responses to the parasite. Here, we used single-cell dual RNA sequencing (scDual-Seq) to comprehensively define, for the first time, the monocyte and parasite transcriptional responses that underpin human monocyte-Toxoplasma encounters at the single cell level. We report extreme transcriptional variability between individual monocytes. Furthermore, we report that Toxoplasma-exposed and unexposed monocytes are transcriptionally distinguished by a reactive subset of CD14+CD16- monocytes. Functional cytokine assays on sorted monocyte populations show that the infection-distinguishing monocytes secrete high levels of chemokines, such as CCL2 and CXCL5. These findings uncover the Toxoplasma-induced monocyte transcriptional heterogeneity and shed new light on the cell populations that largely define cytokine and chemokine secretion in human monocytes exposed to Toxoplasma.

Genome-wide screens identify Toxoplasma gondii determinants of parasite fitness in IFNγ-activated murine 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.

Clonal and atypical Toxoplasma strain differences in virulence vary with mouse sub-species.

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.

Genome-Wide Identification of Splicing Quantitative Trait Loci (sQTLs) in Diverse Ecotypes of Arabidopsis thaliana.

Alternative splicing (AS) of pre-mRNAs contributes to transcriptome diversity and enables plants to generate different protein isoforms from a single gene and/or fine-tune gene expression during different development stages and environmental changes. Although AS is pervasive, the genetic basis for differential isoform usage in plants is still emerging. In this study, we performed genome-wide analysis in 666 geographically distributed diverse ecotypes of Arabidopsis thaliana to identify genomic regions [splicing quantitative trait loci (sQTLs)] that may regulate differential AS. These ecotypes belong to different microclimatic conditions and are part of the relict and non-relict populations. Although sQTLs were spread across the genome, we observed enrichment for trans-sQTL (trans-sQTLs hotspots) on chromosome one. Furthermore, we identified several sQTL (911) that co-localized with trait-linked single nucleotide polymorphisms (SNP) identified in the Arabidopsis genome-wide association studies (AraGWAS). Many sQTLs were enriched among circadian clock, flowering, and stress-responsive genes, suggesting a role for differential isoform usage in regulating these important processes in diverse ecotypes of Arabidopsis. In conclusion, the current study provides a deep insight into SNPs affecting isoform ratios/genes and facilitates a better mechanistic understanding of trait-associated SNPs in GWAS studies. To the best of our knowledge, this is the first report of sQTL analysis in a large set of Arabidopsis ecotypes and can be used as a reference to perform sQTL analysis in the Brassicaceae family. Since whole genome and transcriptome datasets are available for these diverse ecotypes, it could serve as a powerful resource for the biological interpretation of trait-associated loci, splice isoform ratios, and their phenotypic consequences to help produce more resilient and high yield crop varieties.

Three Toxoplasma gondii Dense Granule Proteins Are Required for Induction of Lewis Rat Macrophage Pyroptosis.

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.