NextGen sequencing reveals short double crossovers contribute disproportionately to genetic diversity in Toxoplasma gondii

Brain -cyst-driven genes expression in Toxoplasma Gondii Tehran strain: a parasitic-immunogenicity assessment by dint of RNA-Seq

Toxoplasma gondii (T. gondii) is an obligate intracellular parasite of warm-blooded vertebrates. At present, High-throughput RNA sequencing analysis have made it possible to determine the role of effective genes in host immune response. The aim of the present study is to global transcriptome analysis of the brain of mice infected with T. gondii Tehran strain for the first time and also to evaluate the expression of effective genes in the chronic form of infection. RNA was extracted from the samples and the library was prepared and sequenced using the IlluminaNovaSeq 6000 system. After analyzing gene expression changes, the results were confirmed by real-time method. We found 125 genes that were significantly differentially expressed between infected and non-infected samples (p < 0.0005). Gene ontology analysis revealed that the expression of many genes is critical for pathways such as T cell receptor signaling pathway, Natural Killer cell mediated cytotoxicity, Lysosome and Apoptosis of the host. As infection with Tehran strain leads to chronic infection in mice, therefore, we investigated the genes effective in creating the chronic form of Toxoplasma infection. The comparative analysis of genes showed increases in the expression of genes ctla4, ccl4, cd3e, c3, lcn2, gbp5, usp18, cyba, tap1 and samhd1 in the in the infected sample, which highlights their role in causing chronic infection. RNA-seq provides a valuable tool for analyzing host transcriptomes, better understanding the parasite-host interaction, and developing future drug and vaccine targets.

Mendelian segregation and high recombination rates facilitate genetic analyses in Cryptosporidium parvum

Very little is known about the process of meiosis in the apicomplexan parasite Cryptosporidium despite the essentiality of sex in its life cycle. Most cell lines only support asexual growth of Cryptosporidium parvum (C. parvum), but stem cell derived intestinal epithelial cells grown under air-liquid interface (ALI) conditions support the sexual cycle. To examine chromosomal dynamics during meiosis in C. parvum, we generated two transgenic lines of parasites that were fluorescently tagged with mCherry or GFP on chromosomes 1 or 5, respectively. Infection of ALI cultures or Ifngr1-/- mice with mCherry and GFP parasites resulted in cross-fertilization and the formation of "yellow" oocysts, which contain 4 haploid sporozoites that are the product of meiosis. Recombinant oocysts from the F1 generation were purified and used to infect HCT-8 cultures, and phenotypes of the progeny were observed by microscopy. All possible phenotypes predicted by independent segregation were represented equally (~25%) in the population, indicating that C. parvum chromosomes exhibit a Mendelian inheritance pattern. The most common pattern observed from the outgrowth of single oocysts included all possible parental and recombinant phenotypes derived from a single meiotic event, suggesting a high rate of crossover. To estimate the frequency of crossover, additional loci on chromosomes 1 and 5 were tagged and used to monitor intrachromosomal crosses in Ifngr1-/- mice. Both chromosomes showed a high frequency of crossover compared to other apicomplexans with map distances (i.e., 1% recombination) of 3-12 kb. Overall, a high recombination rate may explain many unique characteristics observed in Cryptosporidium spp. such as high rates of speciation, wide variation in host range, and rapid evolution of host-specific virulence factors.

Mendelian segregation and high recombination rates facilitate genetic analyses in Cryptosporidium parvum

Very little is known about the process of meiosis in the apicomplexan parasite Cryptosporidium despite the essentiality of sex in its life cycle. Most cell lines only support asexual growth of Cryptosporidium parvum (C. parvum), but stem cell derived intestinal epithelial cells grown under air-liquid interface (ALI) conditions support the sexual cycle. To examine chromosomal dynamics during meiosis in C. parvum, we generated two transgenic lines of parasites that were fluorescently tagged with mCherry or GFP on chromosomes 1 or 5, respectively. Infection of ALI cultures or Ifngr1-/- mice with mCherry and GFP parasites produced "yellow" oocysts generated by cross-fertilization. Outcrossed oocysts from the F1 generation were purified and used to infect HCT-8 cultures, and phenotypes of the progeny were observed by microscopy. All possible phenotypes predicted by independent segregation were represented equally (~25%) in the population, indicating that C. parvum chromosomes exhibit a Mendelian inheritance pattern. Unexpectedly, the most common pattern observed from the outgrowth of single oocysts included all possible parental and recombinant phenotypes derived from a single meiotic event, suggesting a high rate of crossover. To estimate the frequency of crossover, additional loci on chromosomes 1 and 5 were tagged and used to monitor intrachromosomal crosses in Ifngr1-/- mice. Both chromosomes showed a high frequency of crossover compared to other apicomplexans with map distances (i.e., 1% recombination) of 3-12 kb. Overall, a high recombination rate may explain many unique characteristics observed in Cryptosporidium spp. such as high rates of speciation, wide variation in host range, and rapid evolution of host-specific virulence factors.

Genotyping of European Toxoplasma gondii strains by a new high-resolution next-generation sequencing-based method

PURPOSE

A new high-resolution next-generation sequencing (NGS)-based method was established to type closely related European type II Toxoplasma gondii strains.

METHODS

T. gondii field isolates were collected from different parts of Europe and assessed by whole genome sequencing (WGS). In comparison to ME49 (a type II reference strain), highly polymorphic regions (HPRs) were identified, showing a considerable number of single nucleotide polymorphisms (SNPs). After confirmation by Sanger sequencing, 18 HPRs were used to design a primer panel for multiplex PCR to establish a multilocus Ion AmpliSeq typing method. Toxoplasma gondii isolates and T. gondii present in clinical samples were typed with the new method. The sensitivity of the method was tested with serially diluted reference DNA samples.

RESULTS

Among type II specimens, the method could differentiate the same number of haplotypes as the reference standard, microsatellite (MS) typing. Passages of the same isolates and specimens originating from abortion outbreaks were identified as identical. In addition, seven different genotypes, two atypical and two recombinant specimens were clearly distinguished from each other by the method. Furthermore, almost all SNPs detected by the Ion AmpliSeq method corresponded to those expected based on WGS. By testing serially diluted DNA samples, the method exhibited a similar analytical sensitivity as MS typing.

CONCLUSION

The new method can distinguish different T. gondii genotypes and detect intra-genotype variability among European type II T. gondii strains. Furthermore, with WGS data additional target regions can be added to the method to potentially increase typing resolution.

Late Embryogenesis Abundant Proteins Contribute to the Resistance of Toxoplasma gondii Oocysts against Environmental Stresses

Toxoplasma gondii oocysts, which are shed in large quantities in the feces from infected felines, are very stable in the environment, resistant to most inactivation procedures, and highly infectious. The oocyst wall provides an important physical barrier for sporozoites contained inside oocysts, protecting them from many chemical and physical stressors, including most inactivation procedures. Furthermore, sporozoites can withstand large temperature changes, even freeze-thawing, as well as desiccation, high salinity, and other environmental insults; however, the genetic basis for this environmental resistance is unknown. Here, we show that a cluster of four genes encoding Late Embryogenesis Abundant (LEA)-related proteins are required to provide Toxoplasma sporozoites resistance to environmental stresses. Toxoplasma LEA-like genes (TgLEAs) exhibit the characteristic features of intrinsically disordered proteins, explaining some of their properties. Our in vitro biochemical experiments using recombinant TgLEA proteins show that they have cryoprotective effects on the oocyst-resident lactate dehydrogenase enzyme and that induced expression in E. coli of two of them leads to better survival after cold stress. Oocysts from a strain in which the four LEA genes were knocked out en bloc were significantly more susceptible to high salinity, freezing, and desiccation compared to wild-type oocysts. We discuss the evolutionary acquisition of LEA-like genes in Toxoplasma and other oocyst-producing apicomplexan parasites of the Sarcocystidae family and discuss how this has likely contributed to the ability of sporozoites within oocysts to survive outside the host for extended periods. Collectively, our data provide a first molecular detailed view on a mechanism that contributes to the remarkable resilience of oocysts against environmental stresses. IMPORTANCE Toxoplasma gondii oocysts are highly infectious and may survive in the environment for years. Their resistance against disinfectants and irradiation has been attributed to the oocyst and sporocyst walls by acting as physical and permeability barriers. However, the genetic basis for their resistance against stressors like changes in temperature, salinity, or humidity, is unknown. We show that a cluster of four genes encoding Toxoplasma Late Embryogenesis Abundant (TgLEA)-related proteins are important for this resistance to environmental stresses. TgLEAs have features of intrinsically disordered proteins, explaining some of their properties. Recombinant TgLEA proteins show cryoprotective effects on the parasite's lactate dehydrogenase, an abundant enzyme in oocysts, and expression in E. coli of two TgLEAs has a beneficial effect on growth after cold stress. Moreover, oocysts from a strain lacking all four TgLEA genes were more susceptible to high salinity, freezing, and desiccation compared to wild-type oocysts, highlighting the importance of the four TgLEAs for oocyst resilience.

Third-generation sequencing revises the molecular karyotype for Toxoplasma gondii and identifies emerging copy number variants in sexual recombinants

Toxoplasma gondii is a useful model for intracellular parasitism given its ease of culture in the laboratory and genomic resources. However, as for many other eukaryotes, the T. gondii genome contains hundreds of sequence gaps owing to repetitive and/or unclonable sequences that disrupt the assembly process. Here, we use the Oxford Nanopore Minion platform to generate near-complete de novo genome assemblies for multiple strains of T. gondii and its near relative, N. caninum. We significantly improved T. gondii genome contiguity (average N50 of ∼6.6 Mb) and added ∼2 Mb of newly assembled sequence. For all of the T. gondii strains that we sequenced (RH, ME49, CTG, II×III progeny clones CL13, S27, S21, S26, and D3X1), the largest contig ranged in size between 11.9 and 12.1 Mb in size, which is larger than any previously reported T. gondii chromosome, and found to be due to a consistent fusion of Chromosomes VIIb and VIII. These data were validated by mapping existing T. gondii ME49 Hi-C data to our assembly, providing parallel lines of evidence that the T. gondii karyotype consists of 13, rather than 14, chromosomes. By using this technology, we also resolved hundreds of tandem repeats of varying lengths, including in well-known host-targeting effector loci like rhoptry protein 5 (ROP5) and ROP38. Finally, when we compared T. gondii with N. caninum, we found that although the 13-chromosome karyotype was conserved, extensive, previously unappreciated chromosome-scale rearrangements had occurred in T. gondii and N. caninum since their most recent common ancestry.

QTL Mapping and CRISPR/Cas9 Editing to Identify a Drug Resistance Gene in Toxoplasma gondii

Scientific knowledge is intrinsically linked to available technologies and methods. This article will present two methods that allowed for the identification and verification of a drug resistance gene in the Apicomplexan parasite Toxoplasma gondii, the method of Quantitative Trait Locus (QTL) mapping using a Whole Genome Sequence (WGS) -based genetic map and the method of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 -based gene editing. The approach of QTL mapping allows one to test if there is a correlation between a genomic region(s) and a phenotype. Two datasets are required to run a QTL scan, a genetic map based on the progeny of a recombinant cross and a quantifiable phenotype assessed in each of the progeny of that cross. These datasets are then formatted to be compatible with R/qtl software that generates a QTL scan to identify significant loci correlated with the phenotype. Although this can greatly narrow the search window of possible candidates, QTLs span regions containing a number of genes from which the causal gene needs to be identified. Having WGS of the progeny was critical to identify the causal drug resistance mutation at the gene level. Once identified, the candidate mutation can be verified by genetic manipulation of drug sensitive parasites. The most facile and efficient method to genetically modify T. gondii is the CRISPR/Cas9 system. This system comprised of just 2 components both encoded on a single plasmid, a single guide RNA (gRNA) containing a 20 bp sequence complementary to the genomic target and the Cas9 endonuclease that generates a double-strand DNA break (DSB) at the target, repair of which allows for insertion or deletion of sequences around the break site. This article provides detailed protocols to use CRISPR/Cas9 based genome editing tools to verify the gene responsible for sinefungin resistance and to construct transgenic parasites.

Genetic Mapping of Pathogenesis Determinants in Toxoplasma gondii

Toxoplasma gondii is a widespread parasite of warm-blooded vertebrates that also causes opportunistic infections in humans. Rodents are a natural host for asexually replicating forms, whereas cats serve as the definitive host for sexual development. The laboratory mouse provides a model to study pathogenesis. Strains of T. gondii are globally diverse, with more than 16 distinct haplogroups clustered into 6 major clades. Forward genetic analysis of genetic crosses between different lineages has been used to define the molecular basis of acute virulence in the mouse. These studies have identified a family of secretory serine/threonine rhoptry kinases that target innate immune pathways to protect intracellular parasites from destruction. Rhoptry kinases target immunity-related GTPases, a family of immune effectors that is expanded in rodents. Similar forward genetic studies may be useful to define the basis of pathogenesis in other hosts, including humans, where infections of different strains present with variable clinical severity.

The Knowns Unknowns: Exploring the Homologous Recombination Repair Pathway in Toxoplasma gondii

Toxoplasma gondii is an apicomplexan parasite of medical and veterinary importance which causes toxoplasmosis in humans. Great effort is currently being devoted toward the identification of novel drugs capable of targeting such illness. In this context, we believe that the thorough understanding of the life cycle of this model parasite will facilitate the identification of new druggable targets in T. gondii. It is important to exploit the available knowledge of pathways which could modulate the sensitivity of the parasite to DNA damaging agents. The homologous recombination repair (HRR) pathway may be of particular interest in this regard as its inactivation sensitizes other cellular models such as human cancer to targeted therapy. Herein we discuss the information available on T. gondii's HRR pathway from the perspective of its conservation with respect to yeast and humans. Special attention was devoted to BRCT domain-containing and end-resection associated proteins in T. gondii as in other experimental models such proteins have crucial roles in early/late steps or HRR and in the pathway choice for double strand break resolution. We conclude that T. gondii HRR pathway is a source of several lines of investigation that allow to to comprehend the extent of diversification of HRR in T. gondii. Such an effort will serve to determine if HRR could represent a potential targer for the treatment of toxoplasmosis.

The Past, Present, and Future of Genetic Manipulation in Toxoplasma gondii

Toxoplasma gondii is a classic model for studying obligate intracellular microorganisms as various genetic manipulation tools have been developed in T. gondii over the past 20 years. Here we summarize the major strategies for T. gondii genetic manipulation including genetic crosses, insertional mutagenesis, chemical mutagenesis, homologous gene replacement, conditional knockdown techniques, and the recently developed clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system. We evaluate the advantages and limitations of each of these tools in a historical perspective. We also discuss additional applications of modified CRISPR-Cas9 systems for use in T. gondii, such as regulation of gene expression, labeling of specific genomic loci, and epigenetic modifications. These approaches have the potential to revolutionize the analysis of T. gondii biology and help us to better develop new drugs and vaccines.

A new coccidian, Isospora parnaitatiaiensis n. sp. (Apicomplexa, Eimeriidae), from the white-shouldered fire-eye Pyriglena leucoptera (Passeriformes, Thamnophilidae) from South America

A new coccidian species (Protozoa: Apicomplexa: Isospora) parasitizing the white-shouldered fire-eye Pyriglena leucoptera (Vieillot, 1818) is described in the Parque Nacional do Itatiaia. This park is a protected area in southeastern Brazil with a high degree of vulnerability, representing a "conservation island" of biodiversity. Isospora parnaitatiaiensis n. sp. has oocysts that are ellipsoidal, 23.8 × 19.4 μm, with smooth, bilayered wall, ~1.1 μm thick. Micropyle and oocyst residuum are absent, but one or two polar granules are present. Sporocysts are ellipsoidal, 14.6 × 9.3 μm. The Stieda body is nipple- to knob-like and sub-Stieda body rounded to rectangular. Sporocyst residuum is present, usually as a cluster of numerous granules. Sporozoites are vermiform with two refractile bodies and a nucleus. This is the second isosporoid coccidian described from antbirds (Thamnophilidae).

Rhoptry Proteins ROP5 and ROP18 Are Major Murine Virulence Factors in Genetically Divergent South American Strains of Toxoplasma gondii

Toxoplasma gondii has evolved a number of strategies to evade immune responses in its many hosts. Previous genetic mapping of crosses between clonal type 1, 2, and 3 strains of T. gondii, which are prevalent in Europe and North America, identified two rhoptry proteins, ROP5 and ROP18, that function together to block innate immune mechanisms activated by interferon gamma (IFNg) in murine hosts. However, the contribution of these and other virulence factors in more genetically divergent South American strains is unknown. Here we utilized a cross between the intermediately virulent North American type 2 ME49 strain and the highly virulent South American type 10 VAND strain to map the genetic basis for differences in virulence in the mouse. Quantitative trait locus (QTL) analysis of this new cross identified one peak that spanned the ROP5 locus on chromosome XII. CRISPR-Cas9 mediated deletion of all copies of ROP5 in the VAND strain rendered it avirulent and complementation confirmed that ROP5 is the major virulence factor accounting for differences between type 2 and type 10 strains. To extend these observations to other virulent South American strains representing distinct genetic populations, we knocked out ROP5 in type 8 TgCtBr5 and type 4 TgCtBr18 strains, resulting in complete loss of virulence in both backgrounds. Consistent with this, polymorphisms that show strong signatures of positive selection in ROP5 were shown to correspond to regions known to interface with host immunity factors. Because ROP5 and ROP18 function together to resist innate immune mechanisms, and a significant interaction between them was identified in a two-locus scan, we also assessed the role of ROP18 in the virulence of South American strains. Deletion of ROP18 in South American type 4, 8, and 10 strains resulted in complete attenuation in contrast to a partial loss of virulence seen for ROP18 knockouts in previously described type 1 parasites. These data show that ROP5 and ROP18 are conserved virulence factors in genetically diverse strains from North and South America, suggesting they evolved to resist innate immune defenses in ancestral T. gondii strains, and they have subsequently diversified under positive selection.

Parasites and the hosts they infect are in constant struggle with each other for survival. On the one hand, the host needs to control parasite growth, while the parasite needs to evade the host response long enough to allow for efficient transmission. The parasite Toxoplasma gondii has evolved virulence factors ROP5 and ROP18 to evade innate immune mechanisms of its natural intermediate host, small rodents. These genes were initially identified in clonal parasite types isolated from Europe and North America, but the factors that contribute to virulence in genetically divergent South American strains have not been tested. Here we used forward and reverse genetic analyses to show that ROP5 and ROP18 are also major virulence factors in genetically distinct virulent South American strains. Given that ROP5 and ROP18 function as virulence factors in strains from North America, Europe, and South America they likely acquired their functions before Toxoplasma gondii radiated into its present global population structure.

REDHORSE-REcombination and Double crossover detection in Haploid Organisms using next-geneRation SEquencing data

Background

Next-generation sequencing technology provides a means to study genetic exchange at a higher resolution than was possible using earlier technologies. However, this improvement presents challenges as the alignments of next generation sequence data to a reference genome cannot be directly used as input to existing detection algorithms, which instead typically use multiple sequence alignments as input. We therefore designed a software suite called REDHORSE that uses genomic alignments, extracts genetic markers, and generates multiple sequence alignments that can be used as input to existing recombination detection algorithms. In addition, REDHORSE implements a custom recombination detection algorithm that makes use of sequence information and genomic positions to accurately detect crossovers. REDHORSE is a portable and platform independent suite that provides efficient analysis of genetic crosses based on Next-generation sequencing data.

Results

We demonstrated the utility of REDHORSE using simulated data and real Next-generation sequencing data. The simulated dataset mimicked recombination between two known haploid parental strains and allowed comparison of detected break points against known true break points to assess performance of recombination detection algorithms. A newly generated NGS dataset from a genetic cross of Toxoplasma gondii allowed us to demonstrate our pipeline. REDHORSE successfully extracted the relevant genetic markers and was able to transform the read alignments from NGS to the genome to generate multiple sequence alignments. Recombination detection algorithm in REDHORSE was able to detect conventional crossovers and double crossovers typically associated with gene conversions whilst filtering out artifacts that might have been introduced during sequencing or alignment. REDHORSE outperformed other commonly used recombination detection algorithms in finding conventional crossovers. In addition, REDHORSE was the only algorithm that was able to detect double crossovers.

Conclusion

REDHORSE is an efficient analytical pipeline that serves as a bridge between genomic alignments and existing recombination detection algorithms. Moreover, REDHORSE is equipped with a recombination detection algorithm specifically designed for Next-generation sequencing data. REDHORSE is portable, platform independent Java based utility that provides efficient analysis of genetic crosses based on Next-generation sequencing data. REDHORSE is available at http://redhorse.sourceforge.net/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1309-7) contains supplementary material, which is available to authorized users.

Genetic mapping reveals that sinefungin resistance in Toxoplasma gondii is controlled by a putative amino acid transporter locus that can be used as a negative selectable marker

Quantitative trait locus (QTL) mapping studies have been integral in identifying and understanding virulence mechanisms in the parasite Toxoplasma gondii. In this study, we interrogated a different phenotype by mapping sinefungin (SNF) resistance in the genetic cross between type 2 ME49-FUDR(r) and type 10 VAND-SNF(r). The genetic map of this cross was generated by whole-genome sequencing of the progeny and subsequent identification of single nucleotide polymorphisms (SNPs) inherited from the parents. Based on this high-density genetic map, we were able to pinpoint the sinefungin resistance phenotype to one significant locus on chromosome IX. Within this locus, a single nonsynonymous SNP (nsSNP) resulting in an early stop codon in the TGVAND_290860 gene was identified, occurring only in the sinefungin-resistant progeny. Using CRISPR/CAS9, we were able to confirm that targeted disruption of TGVAND_290860 renders parasites sinefungin resistant. Because disruption of the SNR1 gene confers resistance, we also show that it can be used as a negative selectable marker to insert either a positive drug selection cassette or a heterologous reporter. These data demonstrate the power of combining classical genetic mapping, whole-genome sequencing, and CRISPR-mediated gene disruption for combined forward and reverse genetic strategies in T. gondii.