IL-33-binding HpARI family homologues with divergent effects in suppressing or enhancing type 2 immune responses.

HpARI is an immunomodulatory protein secreted by the intestinal nematode Heligmosomoides polygyrus bakeri, which binds and blocks IL-33. Here, we find that the H. polygyrus bakeri genome contains three HpARI family members and that these have different effects on IL-33-dependent responses in vitro and in vivo, with HpARI1+2 suppressing and HpARI3 amplifying these responses. All HpARIs have sub-nanomolar affinity for mouse IL-33; however, HpARI3 does not block IL-33-ST2 interactions. Instead, HpARI3 stabilizes IL-33, increasing the half-life of the cytokine and amplifying responses to it in vivo. Together, these data show that H. polygyrus bakeri secretes a family of HpARI proteins with both overlapping and distinct functions, comprising a complex immunomodulatory arsenal of host-targeted proteins.

Small Intestinal Levels of the Branched Short-Chain Fatty Acid Isovalerate Are Elevated during Infection with Heligmosomoides polygyrus and Can Promote Helminth Fecundity.

Heligmosomoides polygyrus is a helminth which naturally infects mice and is widely used as a laboratory model of chronic small intestinal helminth infection. While it is known that infection with H. polygyrus alters the composition of the host's bacterial microbiota, the functional implications of this alteration are unclear. We investigated the impact of H. polygyrus infection on short-chain fatty acid (SCFA) levels in the mouse intestine and sera. We found that helminth infection resulted in significantly upregulated levels of the branched SCFA isovaleric acid, exclusively in the proximal small intestine, which is the site of H. polygyrus colonization. We next set out to test the hypothesis that elevating local levels of isovaleric acid was a strategy used by H. polygyrus to promote its own fitness within the mammalian host. To test this, we supplemented the drinking water of mice with isovalerate during H. polygyrus infection and examined whether this affected helminth fecundity or chronicity. We did not find that isovaleric acid supplementation affected helminth chronicity; however, we found that it did promote helminth fecundity, as measured by helminth egg output in the feces of mice. Through antibiotic treatment of helminth-infected mice, we found that the bacterial microbiota was required in order to support elevated levels of isovaleric acid in the proximal small intestine during helminth infection. Overall, our data reveal that during H. polygyrus infection there is a microbiota-dependent localized increase in the production of isovaleric acid in the proximal small intestine and that this supports helminth fecundity in the murine host.

Improved enzyme annotation with EC-specific cutoffs using DETECT v2.

Summary: We present DETECT v2-an enzyme annotation tool which considers the effect of sequence diversity when assigning enzymatic function [as an Enzyme Commission (EC) number] to a protein sequence. In addition to capturing more enzyme classes than the previous version, we now provide EC-specific cutoffs that greatly increase precision and recall of assignments and show its performance in the context of pathways. Availability and implementation: https://github.com/ParkinsonLab/DETECT-v2. Supplementary information: Supplementary data are available at Bioinformatics online.

The genomes of four tapeworm species reveal adaptations to parasitism.

Tapeworms (Cestoda) cause neglected diseases that can be fatal and are difficult to treat, owing to inefficient drugs. Here we present an analysis of tapeworm genome sequences using the human-infective species Echinococcus multilocularis, E. granulosus, Taenia solium and the laboratory model Hymenolepis microstoma as examples. The 115- to 141-megabase genomes offer insights into the evolution of parasitism. Synteny is maintained with distantly related blood flukes but we find extreme losses of genes and pathways that are ubiquitous in other animals, including 34 homeobox families and several determinants of stem cell fate. Tapeworms have specialized detoxification pathways, metabolism that is finely tuned to rely on nutrients scavenged from their hosts, and species-specific expansions of non-canonical heat shock proteins and families of known antigens. We identify new potential drug targets, including some on which existing pharmaceuticals may act. The genomes provide a rich resource to underpin the development of urgently needed treatments and control.

Does the arthropod microbiota impact the establishment of vector-borne diseases in mammalian hosts?

The impact of the microbiota on the immune status of its host is a source of intense research and publicity. In comparison, the effect of arthropod microbiota on vector-borne infectious diseases has received little attention. A better understanding of the vector microbiota in relation to mammalian host immune responses is vital, as it can lead to strategies that affect transmission and improve vaccine design in a field of research where few vaccines exist and effective treatment is rare. Recent demonstrations of how microbiota decrease pathogen development in arthropods, and thus alter vector permissiveness to vector-borne diseases (VBDs), have led to renewed interest. However, hypotheses on the interactions between the arthropod-derived microbiota and the mammalian hosts have yet to be addressed. Advances in DNA sequencing technology, increased yield and falling costs, mean that these studies are now feasible for many microbiologists and entomologists. Here, we distill current knowledge and put forward key questions and experimental designs to shed light on this burgeoning research topic.

Expanding the view on the evolution of the nematode dauer signalling pathways: refinement through gene gain and pathway co-option.

BACKGROUND: Signalling pathways underlie development, behaviour and pathology. To understand patterns in the evolution of signalling pathways, we undertook a comprehensive investigation of the pathways that control the switch between growth and developmentally quiescent dauer in 24 species of nematodes spanning the phylum. RESULTS: Our analysis of 47 genes across these species indicates that the pathways and their interactions are not conserved throughout the Nematoda. For example, the TGF-ß pathway was co-opted into dauer control relatively late in a lineage that led to the model species Caenorhabditis elegans. We show molecular adaptations described in C. elegans that are restricted to its genus or even just to the species. Similarly, our analyses both identify species where particular genes have been lost and situations where apparently incorrect orthologues have been identified. CONCLUSIONS: Our analysis also highlights the difficulties of working with genome sequences from non-model species as reliance on the published gene models would have significantly restricted our understanding of how signalling pathways evolve. Our approach therefore offers a robust standard operating procedure for genomic comparisons.

Kablammo: an interactive, web-based BLAST results visualizer.

MOTIVATION: Kablammo is a web-based application that produces interactive, vector-based visualizations of sequence alignments generated by BLAST. These visualizations can illustrate many features, including shared protein domains, chromosome structural modifications and genome misassembly. AVAILABILITY AND IMPLEMENTATION: Kablammo can be used at http://kablammo.wasmuthlab.org. For a local installation, the source code and instructions are available under the MIT license at http://github.com/jwintersinger/kablammo. CONTACT: jeff@wintersinger.org.

Figmop: a profile HMM to identify genes and bypass troublesome gene models in draft genomes.

MOTIVATION: Gene models from draft genome assemblies of metazoan species are often incorrect, missing exons or entire genes, particularly for large gene families. Consequently, labour-intensive manual curation is often necessary. We present Figmop (Finding Genes using Motif Patterns) to help with the manual curation of gene families in draft genome assemblies. The program uses a pattern of short sequence motifs to identify putative genes directly from the genome sequence. Using a large gene family as a test case, Figmop was found to be more sensitive and specific than a BLAST-based approach. The visualization used allows the validation of potential genes to be carried out quickly and easily, saving hours if not days from an analysis. AVAILABILITY AND IMPLEMENTATION: Source code of Figmop is freely available for download at https://github.com/dave-the-scientist, implemented in C and Python and is supported on Linux, Unix and MacOSX. CONTACT: curran.dave.m@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comparative genomics of the apicomplexan parasites Toxoplasma gondii and Neospora caninum: Coccidia differing in host range and transmission strategy.

Toxoplasma gondii is a zoonotic protozoan parasite which infects nearly one third of the human population and is found in an extraordinary range of vertebrate hosts. Its epidemiology depends heavily on horizontal transmission, especially between rodents and its definitive host, the cat. Neospora caninum is a recently discovered close relative of Toxoplasma, whose definitive host is the dog. Both species are tissue-dwelling Coccidia and members of the phylum Apicomplexa; they share many common features, but Neospora neither infects humans nor shares the same wide host range as Toxoplasma, rather it shows a striking preference for highly efficient vertical transmission in cattle. These species therefore provide a remarkable opportunity to investigate mechanisms of host restriction, transmission strategies, virulence and zoonotic potential. We sequenced the genome of N. caninum and transcriptomes of the invasive stage of both species, undertaking an extensive comparative genomics and transcriptomics analysis. We estimate that these organisms diverged from their common ancestor around 28 million years ago and find that both genomes and gene expression are remarkably conserved. However, in N. caninum we identified an unexpected expansion of surface antigen gene families and the divergence of secreted virulence factors, including rhoptry kinases. Specifically we show that the rhoptry kinase ROP18 is pseudogenised in N. caninum and that, as a possible consequence, Neospora is unable to phosphorylate host immunity-related GTPases, as Toxoplasma does. This defense strategy is thought to be key to virulence in Toxoplasma. We conclude that the ecological niches occupied by these species are influenced by a relatively small number of gene products which operate at the host-parasite interface and that the dominance of vertical transmission in N. caninum may be associated with the evolution of reduced virulence in this species.

The VPS33B-binding protein VPS16B is required in megakaryocyte and platelet α-granule biogenesis.

Patients with platelet α or dense δ-granule defects have bleeding problems. Although several proteins are known to be required for δ-granule development, less is known about α-granule biogenesis. Our previous work showed that the BEACH protein NBEAL2 and the Sec1/Munc18 protein VPS33B are required for α-granule biogenesis. Using a yeast two-hybrid screen, mass spectrometry, coimmunoprecipitation, and bioinformatics studies, we identified VPS16B as a VPS33B-binding protein. Immunoblotting confirmed VPS16B expression in various human tissues and cells including megakaryocytes and platelets, and also in megakaryocytic Dami cells. Characterization of platelets from a patient with arthrogryposis, renal dysfunction, and cholestasis (ARC) syndrome containing mutations in C14orf133 encoding VPS16B revealed pale-appearing platelets in blood films and electron microscopy revealed a complete absence of α-granules, whereas δ-granules were observed. Soluble and membrane-bound α-granule proteins were reduced or undetectable, suggesting that both releasable and membrane-bound α-granule constituents were absent. Immunofluorescence microscopy of Dami cells stably expressing GFP-VPS16B revealed that similar to VPS33B, GFP-VPS16B colocalized with markers of the trans-Golgi network, late endosomes and α-granules. We conclude that VPS16B, similar to its binding partner VPS33B, is essential for megakaryocyte and platelet α-granule biogenesis.

The impact of FASTQ and alignment read order on structural variant calling from long-read sequencing data.

Background: Structural variant (SV) calling from DNA sequencing data has been challenging due to several factors, including the ambiguity of short-read alignments, multiple complex SVs in the same genomic region, and the lack of "truth" datasets for benchmarking. Additionally, caller choice, parameter settings, and alignment method are known to affect SV calling. However, the impact of FASTQ read order on SV calling has not been explored for long-read data. Results: Here, we used PacBio DNA sequencing data from 15 Caenorhabditis elegans strains and four Arabidopsis thaliana ecotypes to evaluate the sensitivity of different SV callers on FASTQ read order. Comparisons of variant call format files generated from the original and permutated FASTQ files demonstrated that the order of input data affected the SVs predicted by each caller. In particular, pbsv was highly sensitive to the order of the input data, especially at the highest depths where over 70% of the SV calls generated from pairs of differently ordered FASTQ files were in disagreement. These demonstrate that read order sensitivity is a complex, multifactorial process, as the differences observed both within and between species varied considerably according to the specific combination of aligner, SV caller, and sequencing depth. In addition to the SV callers being sensitive to the input data order, the SAMtools alignment sorting algorithm was identified as a source of variability following read order randomization. Conclusion: The results of this study highlight the sensitivity of SV calling on the order of reads encoded in FASTQ files, which has not been recognized in long-read approaches. These findings have implications for the replication of SV studies and the development of consistent SV calling protocols. Our study suggests that researchers should pay attention to the input order sensitivity of read alignment sorting methods when analyzing long-read sequencing data for SV calling, as mitigating a source of variability could facilitate future replication work. These results also raise important questions surrounding the relationship between SV caller read order sensitivity and tool performance. Therefore, tool developers should also consider input order sensitivity as a potential source of variability during the development and benchmarking of new and improved methods for SV calling.

Corrigendum: Trickle infection with Heligmosomoides polygyrus results in decreased worm burdens but increased intestinal inflammation and scarring.

[This corrects the article DOI: 10.3389/fimmu.2022.1020056.].

Identification of potential molecular mimicry in pathogen-host interactions.

Pathogens have evolved sophisticated strategies to manipulate host signaling pathways, including the phenomenon of molecular mimicry, where pathogen-derived biomolecules imitate host biomolecules. In this study, we resurrected, updated, and optimized a sequence-based bioinformatics pipeline to identify potential molecular mimicry candidates between humans and 32 pathogenic species whose proteomes' 3D structure predictions were available at the start of this study. We observed considerable variation in the number of mimicry candidates across pathogenic species, with pathogenic bacteria exhibiting fewer candidates compared to fungi and protozoans. Further analysis revealed that the candidate mimicry regions were enriched in solvent-accessible regions, highlighting their potential functional relevance. We identified a total of 1,878 mimicked regions in 1,439 human proteins, and clustering analysis indicated diverse target proteins across pathogen species. The human proteins containing mimicked regions revealed significant associations between these proteins and various biological processes, with an emphasis on host extracellular matrix organization and cytoskeletal processes. However, immune-related proteins were underrepresented as targets of mimicry. Our findings provide insights into the broad range of host-pathogen interactions mediated by molecular mimicry and highlight potential targets for further investigation. This comprehensive analysis contributes to our understanding of the complex mechanisms employed by pathogens to subvert host defenses and we provide a resource to assist researchers in the development of novel therapeutic strategies.

Transcriptional patterns of sexual dimorphism and in host developmental programs in the model parasitic nematode Heligmosomoides bakeri.

BACKGROUND: Heligmosomoides bakeri (often mistaken for Heligmosomoides polygyrus) is a promising model for parasitic nematodes with the key advantage of being amenable to study and manipulation within a controlled laboratory environment. While draft genome sequences are available for this worm, which allow for comparative genomic analyses between nematodes, there is a notable lack of information on its gene expression. METHODS: We generated biologically replicated RNA-seq datasets from samples taken throughout the parasitic life of H. bakeri. RNA from tissue-dwelling and lumen-dwelling worms, collected under a dissection microscope, was sequenced on an Illumina platform. RESULTS: We find extensive transcriptional sexual dimorphism throughout the fourth larval and adult stages of this parasite and identify alternative splicing, glycosylation, and ubiquitination as particularly important processes for establishing and/or maintaining sex-specific gene expression in this species. We find sex-linked differences in transcription related to aging and oxidative and osmotic stress responses. We observe a starvation-like signature among transcripts whose expression is consistently upregulated in males, which may reflect a higher energy expenditure by male worms. We detect evidence of increased importance for anaerobic respiration among the adult worms, which coincides with the parasite's migration into the physiologically hypoxic environment of the intestinal lumen. Furthermore, we hypothesize that oxygen concentration may be an important driver of the worms encysting in the intestinal mucosa as larvae, which not only fully exposes the worms to their host's immune system but also shapes many of the interactions between the host and parasite. We find stage- and sex-specific variation in the expression of immunomodulatory genes and in anthelmintic targets. CONCLUSIONS: We examine how different the male and female worms are at the molecular level and describe major developmental events that occur in the worm, which extend our understanding of the interactions between this parasite and its host. In addition to generating new hypotheses for follow-up experiments into the worm's behavior, physiology, and metabolism, our datasets enable future more in-depth comparisons between nematodes to better define the utility of H. bakeri as a model for parasitic nematodes in general.

Different structural variant prediction tools yield considerably different results in Caenorhabditis elegans.

The accurate characterization of structural variation is crucial for our understanding of how large chromosomal alterations affect phenotypic differences and contribute to genome evolution. Whole-genome sequencing is a popular approach for identifying structural variants, but the accuracy of popular tools remains unclear due to the limitations of existing benchmarks. Moreover, the performance of these tools for predicting variants in non-human genomes is less certain, as most tools were developed and benchmarked using data from the human genome. To evaluate the use of long-read data for the validation of short-read structural variant calls, the agreement between predictions from a short-read ensemble learning method and long-read tools were compared using real and simulated data from Caenorhabditis elegans. The results obtained from simulated data indicate that the best performing tool is contingent on the type and size of the variant, as well as the sequencing depth of coverage. These results also highlight the need for reference datasets generated from real data that can be used as 'ground truth' in benchmarks.

Draft Hybrid Genome Assembly of a Canadian Cyclospora cayetanensis Isolate.

The apicomplexan parasite Cyclospora cayetanensis causes foodborne gastrointestinal disease in humans. Here, we report the first hybrid assembly for C. cayetanensis, which uses both Illumina MiSeq and Oxford Nanopore Technologies MinION platforms to generate genomic sequence data. The final genome assembly consists of 44,586,677 bases represented in 313 contigs.

Genotyping Canadian Cyclospora cayetanensis Isolates to Supplement Cyclosporiasis Outbreak Investigations.

Cyclospora cayetanensis is an emerging foodborne parasite that causes cyclosporiasis, an enteric disease of humans. Domestically acquired outbreaks have been reported in Canada every spring or summer since 2013. To date, investigations into the potential sources of infection have relied solely on epidemiological data. To supplement the epidemiological data with genetic information, we genotyped 169 Canadian cyclosporiasis cases from stool specimens collected from 2010 to 2021 using an existing eight-marker targeted amplicon deep (TADS) scheme specific to C. cayetanensis as previously described by the US Centers for Disease Control and Prevention (CDC). This is the first study to genotype Canadian Cyclospora cayetanensis isolates, and it focuses on evaluating the genotyping performance and genetic clustering. Genotyping information was successfully collected with at least part of one of the markers in the TADS assay for 97.9% of specimens, and 81.1% of cyclosporiasis cases met the minimum requirements to genetically cluster into 20 groups. The performance of the scheme suggests that examining cyclosporiasis cases genetically will be a valuable tool for supplementing epidemiological outbreak investigations and to minimize further infections. Further research is required to expand the number of discriminatory markers to improve genetic clustering.

Genomic signatures of selection associated with benzimidazole drug treatments in Haemonchus contortus field populations.

Genome-wide methods offer a powerful approach to detect signatures of drug selection. However, limited availability of suitable reference genomes and the difficulty of obtaining field populations with well-defined, distinct drug treatment histories mean there is little information on the signatures of selection in parasitic nematodes and on how best to detect them. This study addresses these knowledge gaps by using field populations of Haemonchus contortus with well-defined benzimidazole treatment histories, leveraging a recently completed chromosomal-scale reference genome assembly. We generated a panel of 49,393 genomic markers to genotype 20 individual adult worms from each of four H. contortus populations: two from closed sheep flocks with an approximate 20 year history of frequent benzimidazole treatment, and two populations with a history of little or no treatment. Sampling occurred in the same geographical region to limit genetic differentiation and maximise the detection sensitivity. A clear signature of selection was detected on chromosome I, centred on the isotype-1 ß-tubulin gene. Two additional, but weaker, signatures of selection were detected; one near the middle of chromosome I spanning 3.75 Mbp and 259 annotated genes, and one on chromosome II spanning a region of 3.3 Mbp and 206 annotated genes, including the isotype-2 ß-tubulin locus. We also assessed how sensitivity was impacted by sequencing depth, worm number, and pooled versus individual worm sequence data. This study provides the first known direct genome-wide evidence for any parasitic nematode, that the isotype-1 ß-tubulin gene is quantitatively the single most important benzimidazole resistance locus. It also identified two additional genomic regions that likely contain benzimidazole resistance loci of secondary importance. This study provides an experimental framework to maximise the power of genome-wide approaches to detect signatures of selection driven by anthelmintic drug treatments in field populations of parasitic nematodes.

Trickle infection with Heligmosomoides polygyrus results in decreased worm burdens but increased intestinal inflammation and scarring.

Introduction: Intestinal roundworms cause chronic debilitating disease in animals, including humans. Traditional experimental models of these types of infection use a large single-dose infection. However, in natural settings, hosts are exposed to parasites on a regular basis and when mice are exposed to frequent, smaller doses of Heligmosomoides polygyrus, the parasites are cleared more quickly. Whether this more effective host response has any negative consequences for the host is not known. Results: Using a trickle model of infection, we found that worm clearance was associated with known resistance-related host responses: increased granuloma and tuft cell numbers, increased levels of granuloma IgG and decreased intestinal transit time, as well as higher serum IgE levels. However, we found that the improved worm clearance was also associated with an inflammatory phenotype in and around the granuloma, increased smooth muscle hypertrophy/hyperplasia, and elevated levels of Adamts gene expression. Discussion: To our knowledge, we are the first to identify the involvement of this protein family of matrix metalloproteinases (MMPs) in host responses to helminth infections. Our results highlight the delicate balance between parasite clearance and host tissue damage, which both contribute to host pathology. When continually exposed to parasitic worms, improved clearance comes at a cost.

Characterization of the ß-tubulin gene family in Ascaris lumbricoides and Ascaris suum and its implication for the molecular detection of benzimidazole resistance.

BACKGROUND: The treatment coverage of control programs providing benzimidazole (BZ) drugs to eliminate the morbidity caused by soil-transmitted helminths (STHs) is unprecedently high. This high drug pressure may result in the development of BZ resistance in STHs and so there is an urgent need for surveillance systems detecting molecular markers associated with BZ resistance. A critical prerequisite to develop such systems is an understanding of the gene family encoding ß-tubulin proteins, the principal targets of BZ drugs. METHODOLOGY AND PRINCIPAL FINDINGS: First, the ß-tubulin gene families of Ascaris lumbricoides and Ascaris suum were characterized through the analysis of published genomes. Second, RNA-seq and RT-PCR analyses on cDNA were applied to determine the transcription profiles of the different gene family members. The results revealed that Ascaris species have at least seven different ß-tubulin genes of which two are highly expressed during the entire lifecycle. Third, deep amplicon sequencing was performed on these two genes in more than 200 adult A. lumbricoides (Ethiopia and Tanzania) and A. suum (Belgium) worms, to investigate the intra- and inter-species genetic diversity and the presence of single nucleotide polymorphisms (SNPs) that are associated with BZ resistance in other helminth species; F167Y (TTC>TAC or TTT>TAT), E198A (GAA>GCA or GAG>GCG), E198L (GAA>TTA) and F200Y (TTC>TAC or TTT>TAT). These particular SNPs were absent in the two investigated genes in all three Ascaris populations. SIGNIFICANCE: This study demonstrated the presence of at least seven ß-tubulin genes in Ascaris worms. A new nomenclature was proposed and prioritization of genes for future BZ resistance research was discussed. This is the first comprehensive description of the ß-tubulin gene family in Ascaris and provides a framework to investigate the prevalence and potential role of ß-tubulin sequence polymorphisms in BZ resistance in a more systematic manner than previously possible.