Temporal gene expression during asexual development of the apicomplexan Sarcocystis neurona.

Asexual replication in the apicomplexan Sarcocystis neurona involves two main developmental stages: the motile extracellular merozoite and the sessile intracellular schizont. Merozoites invade host cells and transform into schizonts that undergo replication via endopolygeny to form multiple (64) daughter merozoites that are invasive to new host cells. Given that the capabilities of the merozoite vary significantly from the schizont, the patterns of transcript levels throughout the asexual lifecycle were determined and compared in this study. RNA-Seq data were generated from extracellular merozoites and four intracellular schizont development time points. Of the 6,938 genes annotated in the S. neurona genome, 6,784 were identified in the transcriptome. Of these, 4,111 genes exhibited significant differential expression between the merozoite and at least one schizont development time point. Transcript levels were significantly higher for 2,338 genes in the merozoite and 1,773 genes in the schizont stages. Included in this list were genes encoding the secretory pathogenesis determinants (SPDs), which encompass the surface antigen and SAG-related sequence (SAG/SRS) and the secretory organelle proteins of the invasive zoite stage (micronemes, rhoptries, and dense granules). As anticipated, many of the S. neurona SPD gene transcripts were abundant in merozoites. However, several SPD transcripts were elevated in intracellular schizonts, suggesting roles unrelated to host cell invasion and the initial establishment of the intracellular niche. The hypothetical genes that are potentially unique to the genus Sarcocystis are of particular interest. Their conserved expression patterns are instructive for future investigations into the possible functions of these putative Sarcocystis-unique genes. IMPORTANCE: The genus Sarcocystis is an expansive clade within the Apicomplexa, with the species S. neurona being an important cause of neurological disease in horses. Research to decipher the biology of S. neurona and its host-pathogen interactions can be enhanced by gene expression data. This study has identified conserved apicomplexan orthologs in S. neurona, putative Sarcocystis-unique genes, and gene transcripts abundant in the merozoite and schizont stages. Importantly, we have identified distinct clusters of genes with transcript levels peaking during different intracellular schizont development time points, reflecting active gene expression changes across endopolygeny. Each cluster also has subsets of transcripts with unknown functions, and investigation of these seemingly Sarcocystis-unique transcripts will provide insights into the interesting biology of this parasite genus.

A new chromosome-level genome assembly and annotation of Cryptosporidium meleagridis.

Cryptosporidium spp. are medically and scientifically relevant protozoan parasites that cause severe diarrheal illness in infants and immunosuppressed populations as well as animals. Although most human Cryptosporidium infections are caused by C. parvum and C. hominis, there are several other human-infecting species including C. meleagridis, which is commonly observed in developing countries. Here, we polished and annotated a long-read genome sequence assembly for C. meleagridis TU1867, a species which infects birds and humans. The genome sequence was generated using a combination of whole genome amplification (WGA) and long-read Oxford Nanopore Technologies sequencing. The assembly was then polished with Illumina data. The chromosome-level genome assembly is 9.2 Mbp with a contig N50 of 1.1 Mb. Annotation revealed 3,923 protein-coding genes. A BUSCO analysis indicates a completeness of 96.6% (n=446), including 430 (96.4%) single-copy and 1 (0.224%) duplicated apicomplexan conserved gene(s). The new C. meleagridis genome assembly is nearly gap-free and provides a valuable new resource for the Cryptosporidium community and future studies on evolution and host-specificity.

"Evaluating the Benefits and Limits of Multiple Displacement Amplification with Whole-Genome Oxford Nanopore Sequencing".

Multiple Displacement Amplification (MDA) outperforms conventional PCR in long fragment and whole genome amplification which makes it attractive to couple with long-read sequencing of samples with limited quantities of DNA to obtain improved genome assemblies. Here, we explore the efficacy and limits of MDA for genome sequence assembly using Oxford Nanopore Technologies (ONT) rapid library preparations and minION sequencing. We successfully generated almost complete genome sequences for all organisms examined, including Cryptosporidium meleagridis, Staphylococcus aureus, Enterococcus faecium, and Escherichia coli, with the ability to generate high-quality data from samples starting with only 0.025 ng of total DNA. Controlled sheared DNA samples exhibited a distinct pattern of size-increase after MDA, which may be associated with the amplification of long, low-abundance fragments present in the assay, as well as generating concatemeric sequences during amplification. To address concatemers, we developed a computational pipeline (CADECT: Concatemer Detection Tool) to identify and remove putative concatemeric sequences. This study highlights the efficacy of MDA in generating high-quality genome assemblies from limited amounts of input DNA. Also, the CADECT pipeline effectively mitigated the impact of concatemeric sequences, enabling the assembly of contiguous sequences even in cases where the input genomic DNA was degraded. These results have significant implications for the study of organisms that are challenging to culture in vitro, such as Cryptosporidium, and for expediting critical results in clinical settings with limited quantities of available genomic DNA.

Cefiderocol heteroresistance associated with mutations in TonB-dependent receptor genes in Pseudomonas aeruginosa of clinical origin.

The siderophore-cephalosporin cefiderocol(FDC) presents a promising treatment option for carbapenem-resistant (CR) P. aeruginosa (PA). FDC circumvents traditional porin and efflux mediated resistance by utilizing TonB-dependent receptors (TBDRs) to access the periplasmic space. Emerging FDC resistance has been associated with loss of function mutations within TBDR genes or the regulatory genes controlling TBDR expression. Further, difficulties with antimicrobial susceptibility testing (AST) and unexpected negative clinical treatment outcomes have prompted concerns for heteroresistance, where a single lineage isolate contains resistant subpopulations not detectable by standard AST. This study aimed to evaluate the prevalence of TBDR mutations among clinical isolates of P. aeruginosa and the phenotypic effect on FDC susceptibility and heteroresistance. We evaluated the sequence of pirR , pirS , pirA , piuA or piuD from 498 unique isolates collected before the introduction of FDC from 4 clinical sites in Portland, OR (1), Houston, TX (2), and Santiago, Chile (1). At some clinical sites, TBDR mutations were seen in up to 25% of isolates, and insertion, deletion, or frameshift mutations were predicted to impair protein function were seen in 3% of all isolates (n=15). Using population analysis profile testing, we found that P. aeruginosa with major TBDR mutations were enriched for a heteroresistant phenotype and undergo a shift in the susceptibility distribution of the population as compared to susceptible strains with wild type TBDR genes. Our results indicate that mutations in TBDR genes predate the clinical introduction of FDC, and these mutations may predispose to the emergence of FDC resistance.

Genomic and virulence analysis of in vitro cultured Cryptosporidium parvum.

Recent advances in the in vitro cultivation of Cryptosporidium parvum using hollow fiber bioreactor technology (HFB) have permitted continuous growth of parasites that complete all life cycle stages. The method provides access to all stages of the parasite and provides a method for non-animal production of oocysts for use in clinical trials. Here we examined the effect of long-term (>20 months) in vitro culture on virulence-factors, genome conservation, and in vivo pathogenicity of the host by in vitro cultured parasites. We find low-level sequence variation that is consistent with that observed in calf-passaged parasites. Further using a calf model infection, oocysts obtained from the HFB caused diarrhea of the same volume, duration and oocyst shedding intensity as in vivo passaged parasites.

Epidemiologic investigation and genetic characterization of canine respiratory coronavirus in the Southeastern United States.

Canine respiratory coronavirus (CRCoV) is one of the main causative agents of canine infectious respiratory disease (CIRD), an illness whose epidemiology is poorly understood. We assessed the prevalence, risk factors, and genetic characterization of CRCoV in privately owned dogs in the Southeastern United States. We PCR-screened 189 nasal swabs from dogs with and without CIRD clinical signs for 9 CIRD-related pathogens, including CRCoV; 14% of dogs, all diagnosed with CIRD, were positive for CRCoV, with a significantly higher rate of cases in younger dogs and during warmer weather. Notably, the presence of CRCoV, alone or in coinfection with other CIRD pathogens, was statistically associated with a worse prognosis. We estimated a CRCoV seroprevalence of 23.7% retrospectively from 540 serum samples, with no statistical association to dog age, sex, or season, but with a significantly higher presence in urban counties. Additionally, the genomes of 6 CRCoVs were obtained from positive samples using an in-house developed targeted amplicon-based approach specific to CRCoV. Subsequent phylogeny clustered their genomes in 2 distinct genomic groups, with most isolates sharing a higher similarity with CRCoVs from Sweden and only 1 more closely related to CRCoVs from Asia. We provide new insights into CIRD and CRCoV epidemiology in the Southeastern United States and further support the association of CRCoV with more severe cases of CIRD. Additionally, we developed and successfully tested a new amplicon-based approach for whole-genome sequencing of CRCoV that can be used to further investigate the genetic diversity within CRCoVs.

The emergence of cefiderocol resistance in Pseudomonas aeruginosa from a heteroresistant isolate during prolonged therapy.

Cefiderocol is a siderophore cephalosporin designed to target multi-drug-resistant Gram-negative bacteria. Previously, the emergence of cefiderocol non-susceptibility has been associated with mutations in the chromosomal cephalosporinase (PDC) along with mutations in the PirA and PiuA/D TonB-dependent receptor pathways. Here, we report a clinical case of cefiderocol-resistant P. aeruginosa that emerged in a patient during treatment. This resistance was associated with mutations not previously reported, suggesting potential novel pathways to cefiderocol resistance.

Inherently Reduced Expression of ASC Restricts Caspase-1 Processing in Hepatocytes and Promotes Plasmodium Infection.

Inflammasome-mediated caspase-1 activation facilitates innate immune control of Plasmodium in the liver, thereby limiting the incidence and severity of clinical malaria. However, caspase-1 processing occurs incompletely in both mouse and human hepatocytes and precludes the generation of mature IL-1ß or IL-18, unlike in other cells. Why this is so or how it impacts Plasmodium control in the liver has remained unknown. We show that an inherently reduced expression of the inflammasome adaptor molecule apoptosis-associated specklike protein containing CARD (ASC) is responsible for the incomplete proteolytic processing of caspase-1 in murine hepatocytes. Transgenically enhancing ASC expression in hepatocytes enabled complete caspase-1 processing, enhanced pyroptotic cell death, maturation of the proinflammatory cytokines IL-1ß and IL-18 that was otherwise absent, and better overall control of Plasmodium infection in the liver of mice. This, however, impeded the protection offered by live attenuated antimalarial vaccination. Tempering ASC expression in mouse macrophages, on the other hand, resulted in incomplete processing of caspase-1. Our work shows how caspase-1 activation and function in host cells are fundamentally defined by ASC expression and offers a potential new pathway to create better disease and vaccination outcomes by modifying the latter.

Genetic crosses within and between species of Cryptosporidium.

Parasites and their hosts are engaged in reciprocal coevolution that balances competing mechanisms of virulence, resistance, and evasion. This often leads to host specificity, but genomic reassortment between different strains can enable parasites to jump host barriers and conquer new niches. In the apicomplexan parasite Cryptosporidium, genetic exchange has been hypothesized to play a prominent role in adaptation to humans. The sexual lifecycle of the parasite provides a potential mechanism for such exchange; however, the boundaries of Cryptosporidium sex are currently undefined. To explore this experimentally, we established a model for genetic crosses. Drug resistance was engineered using a mutated phenylalanyl tRNA synthetase gene and marking strains with this and the previously used Neo transgene enabled selection of recombinant progeny. This is highly efficient, and genomic recombination is evident and can be continuously monitored in real time by drug resistance, flow cytometry, and PCR mapping. Using this approach, multiple loci can now be modified with ease. We demonstrate that essential genes can be ablated by crossing a Cre recombinase driver strain with floxed strains. We further find that genetic crosses are also feasible between species. Crossing Cryptosporidium parvum, a parasite of cattle and humans, and Cryptosporidium tyzzeri a mouse parasite resulted in progeny with a recombinant genome derived from both species that continues to vigorously replicate sexually. These experiments have important fundamental and translational implications for the evolution of Cryptosporidium and open the door to reverse- and forward-genetic analysis of parasite biology and host specificity.

Genetic crosses within and between species of Cryptosporidium.

Parasites and their hosts are engaged in rapid coevolution that balances competing mechanisms of virulence, resistance, and evasion. This often leads to host specificity, but genomic reassortment between different strains can enable parasites to jump host barriers and conquer new niches. In the apicomplexan parasite Cryptosporidium genetic exchange has been hypothesized to play a prominent role in adaptation to humans. The sexual lifecycle of the parasite provides a potential mechanism for such exchange; however, the boundaries of Cryptosporidium sex are currently undefined. To explore this experimentally, we established a model for genetic crosses. Drug resistance was engineered using a mutated phenylalanyl tRNA synthetase gene and marking strains with this and the previously used Neo transgene enabled selection of recombinant progeny. This is highly efficient, and genomic recombination is evident and can be continuously monitored in real time by drug resistance, flow cytometry, and PCR mapping. Using this approach multiple loci can now be modified with ease. We demonstrate that essential genes can be ablated by crossing a Cre recombinase driver strain with floxed strains. We further find that genetic crosses are also feasible between species. Crossing C. parvum, a parasite of cattle and humans, and C. tyzzeri a mouse parasite resulted in progeny with a recombinant genome derived from both species that continues to vigorously replicate sexually. These experiments have important fundamental and translational implications for the evolution of Cryptosporidium and open the door to reverse- and forward- genetic analysis of parasite biology and host specificity.

New T2T assembly of Cryptosporidium parvum IOWA annotated with reference genome gene identifiers.

Cryptosporidium parvum is a significant pathogen causing gastrointestinal infections in humans and animals, that is spread through the ingestion of contaminated food and water. Despite its global impact on public health, generating a C. parvum genome sequence has always been challenging due to a lack of in vitro cultivation systems and challenging sub-telomeric gene families. A gapless telomere to telomere genome assembly has been created for Cryptosporidium parvum IOWA obtained from Bunch Grass Farms, named here as CpBGF. There are 8 chromosomes that total 9,259,183 bp. The new hybrid assembly which was generated with Illumina and Oxford Nanopore resolves complex sub-telomeric regions of chromosomes 1, 7 and 8. To facilitate ease of use and consistency with the literature, whenever possible, chromosomes have been oriented and genes in this annotation have been given the same gene IDs used in the current reference genome sequence generated in 2004. The annotation of this assembly utilized considerable RNA expression evidence, thus, untranslated regions, long noncoding RNAs and antisense RNAs are annotated. The CpBGF genome assembly serves as a valuable resource for understanding the biology, pathogenesis, and transmission of C. parvum, and it facilitates the development of diagnostics, drugs, and vaccines against cryptosporidiosis.

Membrane Phenotypic, Metabolic and Genotypic Adaptations of Streptococcus oralis Strains Destined to Rapidly Develop Stable, High-Level Daptomycin Resistance during Daptomycin Exposures.

The Streptococcus mitis-oralis subgroup of viridans group streptococci are important human pathogens. We previously showed that a substantial portion of S. mitis-oralis strains (>25%) are 'destined' to develop rapid, high-level, and stable daptomycin (DAP) resistance (DAP-R) during DAP exposures in vitro. Such DAP-R is often accompanied by perturbations in distinct membrane phenotypes and metabolic pathways. The current study evaluated two S. oralis bloodstream isolates, 73 and 205. Strain 73 developed stable, high-level DAP-R (minimum inhibitory concentration [MIC] > 256 µg/mL) within 2 days of in vitro DAP passage ("high level" DAP-R [HLDR]). In contrast, strain 205 evolved low-level and unstable DAP-R (MIC = 8 µg/mL) under the same exposure conditions in vitro ("non-HLDR"). Comparing the parental 73 vs. 73-D2 (HLDR) strain-pair, we observed the 73-D2 had the following major differences: (i) altered cell membrane (CM) phospholipid profiles, featuring the disappearance of phosphatidylglycerol (PG) and cardiolipin (CL), with accumulation of the PG-CL pathway precursor, phosphatidic acid (PA); (ii) enhanced CM fluidity; (iii) increased DAP surface binding; (iv) reduced growth rates; (v) decreased glucose utilization and lactate accumulation; and (vi) increased enzymatic activity within the glycolytic (i.e., lactate dehydrogenase [LDH]) and lipid biosynthetic (glycerol-3-phosphate dehydrogenase [GPDH]) pathways. In contrast, the 205 (non-HLDR) strain-pair did not show these same phenotypic or metabolic changes over the 2-day DAP exposure. WGS analyses confirmed the presence of mutations in genes involved in the above glycolytic and phospholipid biosynthetic pathways in the 73-D2 passage variant. These data suggest that S. oralis strains which are 'destined' to rapidly develop HLDR do so via a conserved cadre of genotypic, membrane phenotypic, and metabolic adaptations.

An X-Domain Phosphoinositide Phospholipase C (PI-PLC-like) of Trypanosoma brucei Has a Surface Localization and Is Essential for Proliferation.

Trypanosoma brucei is the causative agent of African trypanosomiasis, a deadly disease that affects humans and cattle. There are very few drugs to treat it, and there is evidence of mounting resistance, raising the need for new drug development. Here, we report the presence of a phosphoinositide phospholipase C (TbPI-PLC-like), containing an X and a PDZ domain, that is similar to the previously characterized TbPI-PLC1. TbPI-PLC-like only possesses the X catalytic domain and does not have the EF-hand, Y, and C2 domains, having instead a PDZ domain. Recombinant TbPI-PLC-like does not hydrolyze phosphatidylinositol 4,5-bisphosphate (PIP2) and does not modulate TbPI-PLC1 activity in vitro. TbPI-PLC-like shows a plasma membrane and intracellular localization in permeabilized cells and a surface localization in non-permeabilized cells. Surprisingly, knockdown of TbPI-PLC-like expression by RNAi significantly affected proliferation of both procyclic and bloodstream trypomastigotes. This is in contrast with the lack of effect of downregulation of expression of TbPI-PLC1.

AIM2 sensors mediate immunity to Plasmodium infection in hepatocytes.

Malaria, caused by Plasmodium parasites is a severe disease affecting millions of people around the world. Plasmodium undergoes obligatory development and replication in the hepatocytes, before initiating the life-threatening blood-stage of malaria. Although the natural immune responses impeding Plasmodium infection and development in the liver are key to controlling clinical malaria and transmission, those remain relatively unknown. Here we demonstrate that the DNA of Plasmodium parasites is sensed by cytosolic AIM2 (absent in melanoma 2) receptors in the infected hepatocytes, resulting in Caspase-1 activation. Remarkably, Caspase-1 was observed to undergo unconventional proteolytic processing in hepatocytes, resulting in the activation of the membrane pore-forming protein, Gasdermin D, but not inflammasome-associated proinflammatory cytokines. Nevertheless, this resulted in the elimination of Plasmodium-infected hepatocytes and the control of malaria infection in the liver. Our study uncovers a pathway of natural immunity critical for the control of malaria in the liver.

Accessing the Variability of Multicopy Genes in Complex Genomes using Unassembled Next-Generation Sequencing Reads: The Case of Trypanosoma cruzi Multigene Families.

Repetitive elements cause assembly fragmentation in complex eukaryotic genomes, limiting the study of their variability. The genome of Trypanosoma cruzi, the parasite that causes Chagas disease, has a high repetitive content, including multigene families. Although many T. cruzi multigene families encode surface proteins that play pivotal roles in host-parasite interactions, their variability is currently underestimated, as their high repetitive content results in collapsed gene variants. To estimate sequence variability and copy number variation of multigene families, we developed a read-based approach that is independent of gene-specific read mapping and de novo assembly. This methodology was used to estimate the copy number and variability of MASP, TcMUC, and Trans-Sialidase (TS), the three largest T. cruzi multigene families, in 36 strains, including members of all six parasite discrete typing units (DTUs). We found that these three families present a specific pattern of variability and copy number among the distinct parasite DTUs. Inter-DTU hybrid strains presented a higher variability of these families, suggesting that maintaining a larger content of their members could be advantageous. In addition, in a chronic murine model and chronic Chagasic human patients, the immune response was focused on TS antigens, suggesting that targeting TS conserved sequences could be a potential avenue to improve diagnosis and vaccine design against Chagas disease. Finally, the proposed approach can be applied to study multicopy genes in any organism, opening new avenues to access sequence variability in complex genomes. IMPORTANCE Sequences that have several copies in a genome, such as multicopy-gene families, mobile elements, and microsatellites, are among the most challenging genomic segments to study. They are frequently underestimated in genome assemblies, hampering the correct assessment of these important players in genome evolution and adaptation. Here, we developed a new methodology to estimate variability and copy numbers of repetitive genomic regions and employed it to characterize the T. cruzi multigene families MASP, TcMUC, and transsialidase (TS), which are important virulence factors in this parasite. We showed that multigene families vary in sequence and content among the parasite's lineages, whereas hybrid strains have a higher sequence variability that could be advantageous to the parasite's survivability. By identifying conserved sequences within multigene families, we showed that the mammalian host immune response toward these multigene families is usually focused on the TS multigene family. These TS conserved and immunogenic peptides can be explored in future works as diagnostic targets or vaccine candidates for Chagas disease. Finally, this methodology can be easily applied to any organism of interest, which will aid in our understanding of complex genomic regions.

Three Distinct Annotation Platforms Differ in Detection of Antimicrobial Resistance Genes in Long-Read, Short-Read, and Hybrid Sequences Derived from Total Genomic DNA or from Purified Plasmid DNA.

Recent advances and lower costs in rapid high-throughput sequencing have engendered hope that whole genome sequencing (WGS) might afford complete resistome characterization in bacterial isolates. WGS is particularly useful for the clinical characterization of fastidious and slow-growing bacteria. Despite its potential, several challenges should be addressed before adopting WGS to detect antimicrobial resistance (AMR) genes in the clinical laboratory. Here, with three distinct ESKAPE bacteria (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.), different approaches were compared to identify best practices for detecting AMR genes, including: total genomic DNA and plasmid DNA extractions, the solo assembly of Illumina short-reads and of Oxford Nanopore Technologies (ONT) long-reads, two hybrid assembly pipelines, and three in silico AMR databases. We also determined the susceptibility of each strain to 21 antimicrobials. We found that all AMR genes detected in pure plasmid DNA were also detectable in total genomic DNA, indicating that, at least in these three enterobacterial genera, the purification of plasmid DNA was not necessary to detect plasmid-borne AMR genes. Illumina short-reads used with ONT long-reads in either hybrid or polished assemblies of total genomic DNA enhanced the sensitivity and accuracy of AMR gene detection. Phenotypic susceptibility closely corresponded with genotypes identified by sequencing; however, the three AMR databases differed significantly in distinguishing mobile dedicated AMR genes from non-mobile chromosomal housekeeping genes in which rare spontaneous resistance mutations might occur. This study indicates that each method employed in a WGS workflow has an impact on the detection of AMR genes. A combination of short- and long-reads, followed by at least three different AMR databases, should be used for the consistent detection of such genes. Further, an additional step for plasmid DNA purification and sequencing may not be necessary. This study reveals the need for standardized biochemical and informatic procedures and database resources for consistent, reliable AMR genotyping to take full advantage of WGS in order to expedite patient treatment and track AMR genes within the hospital and community.

Replacement of Leishmania (Leishmania) infantum Populations in an Endemic Focus of Visceral Leishmaniasis in Brazil.

Visceral leishmaniasis is an important global health problem with an estimated of 50,000 to 90,000 new cases per year. VL is the most serious form of leishmaniasis as it can be fatal in 95% of the cases if it remains untreated. VL is a particularly acute problem in Brazil which contributed with 97% of all cases reported in 2020 in the Americas. In this country, VL affects mainly the poorest people in both urban and rural areas and continues to have a high mortality rate estimated around 8.15%. Here, we performed a temporal parasite population study using whole genome sequence data from a set of 34 canine isolates sampled in 2008, 2012 and 2015 from a re-emergent focus in Southeastern Brazil. Our study found the presence of two distinct sexual subpopulations that corresponded to two isolation periods. These subpopulations diverged hundreds of years ago with no apparent gene flow between them suggesting a process of rapid replacement during a two-year period. Sequence comparisons and analysis of nucleotide diversity also showed evidence of balancing selection acting on transport-related genes and antigenic families. To our knowledge this is the first population genomic study showing a turn-over of parasite populations in an endemic region for leishmaniasis. The complexity and rapid adaptability of these parasites pose new challenges to control activities and demand more integrated approaches to understand this disease in New World foci.

Small and intermediate size structural RNAs in the unicellular parasite Cryptosporidium parvum as revealed by sRNA-seq and comparative genomics.

Small and intermediate-size noncoding RNAs (sRNAs and is-ncRNAs) have been shown to play important regulatory roles in the development of several eukaryotic organisms. However, they have not been thoroughly explored in Cryptosporidium parvum, an obligate zoonotic protist parasite responsible for the diarrhoeal disease cryptosporidiosis. Using Illumina sequencing of a small RNA library, a systematic identification of novel small and is-ncRNAs was performed in C. parvum excysted sporozoites. A total of 79 novel is-ncRNA candidates, including antisense, intergenic and intronic is-ncRNAs, were identified, including 7 new small nucleolar RNAs (snoRNAs). Expression of select novel is-ncRNAs was confirmed by RT-PCR. Phylogenetic conservation was analysed using covariance models (CMs) in related Cryptosporidium and apicomplexan parasite genome sequences. A potential new type of small ncRNA derived from tRNA fragments was observed. Overall, a deep profiling analysis of novel is-ncRNAs in C. parvum and related species revealed structural features and conservation of these novel is-ncRNAs. Covariance models can be used to detect is-ncRNA genes in other closely related parasites. These findings provide important new sequences for additional functional characterization of novel is-ncRNAs in the protist pathogen C. parvum.