The parasitic nematode Strongyloides ratti exists predominantly as populations of long-lived asexual lineages.

Nematodes are important parasites of people and animals, and in natural ecosystems they are a major ecological force. Strongyloides ratti is a common parasitic nematode of wild rats and we have investigated its population genetics using single-worm, whole-genome sequencing. We find that S. ratti populations in the UK consist of mixtures of mainly asexual lineages that are widely dispersed across a host population. These parasite lineages are likely very old and may have originated in Asia from where rats originated. Genes that underly the parasitic phase of the parasite's life cycle are hyperdiverse compared with the rest of the genome, and this may allow the parasites to maximise their fitness in a diverse host population. These patterns of parasitic nematode population genetics have not been found before and may also apply to Strongyloides spp. that infect people, which will affect how we should approach their control.

The Human Blood Fluke, Schistosoma mansoni, Harbors Bacteria Throughout the Parasite's Life Cycle.

While symbiotic relationships between invertebrates and bacteria have been extensively described, studies of microbial communities inhabiting parasitic worms remain scarce. Exploring the microbiota associated with helminths responsible for major infectious diseases will inform on parasite biology, host-pathogen interactions, and disease pathophysiology. We investigated the presence of microorganisms inhabiting tissues of the human parasite Schistosoma mansoni. In situ hybridization using a pan-bacterial 16S rRNA gene probe revealed bacteria colonizing key developmental stages that were successfully removed after antibiotic treatment of live parasites. Understanding the composition and function of the S. mansoni-associated microbiota may lead to the development of novel microbiome-targeting control strategies.

Geographic Origin and Vertical Transmission of Leishmania infantum Parasites in Hunting Hounds, United States.

Vertical transmission of leishmaniasis is common but is difficult to study against the background of pervasive vector transmission. We present genomic data from dogs in the United States infected with Leishmania infantum parasites; these infections have persisted in the apparent absence of vector transmission. We demonstrate that these parasites were introduced from the Old World separately and more recently than L. infantum from South America. The parasite population shows unusual genetics consistent with a lack of meiosis: a high level of heterozygous sites shared across all isolates and no decrease in linkage with genomic distance between variants. Our data confirm that this parasite population has been evolving with little or no sexual reproduction. This demonstration of vertical transmission has profound implications for the population genetics of Leishmania parasites. When investigating transmission in complex natural settings, considering vertical transmission alongside vector transmission is vital.

WormBase in 2022-data, processes, and tools for analyzing Caenorhabditis elegans.

WormBase (www.wormbase.org) is the central repository for the genetics and genomics of the nematode Caenorhabditis elegans. We provide the research community with data and tools to facilitate the use of C. elegans and related nematodes as model organisms for studying human health, development, and many aspects of fundamental biology. Throughout our 22-year history, we have continued to evolve to reflect progress and innovation in the science and technologies involved in the study of C. elegans. We strive to incorporate new data types and richer data sets, and to provide integrated displays and services that avail the knowledge generated by the published nematode genetics literature. Here, we provide a broad overview of the current state of WormBase in terms of data type, curation workflows, analysis, and tools, including exciting new advances for analysis of single-cell data, text mining and visualization, and the new community collaboration forum. Concurrently, we continue the integration and harmonization of infrastructure, processes, and tools with the Alliance of Genome Resources, of which WormBase is a founding member.

Genomic and Phenotypic Characterization of Experimentally Selected Resistant Leishmania donovani Reveals a Role for Dynamin-1-Like Protein in the Mechanism of Resistance to a Novel Antileishmanial Compound.

The implementation of prospective drug resistance (DR) studies in the research-and-development (R&D) pipeline is a common practice for many infectious diseases but not for neglected tropical diseases (NTDs). Here, we explored and demonstrated the importance of this approach using as paradigms Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GlaxoSmithKline (GSK) "Leishbox" to treat VL. We experimentally selected resistance to TCMDC-143345 in vitro and characterized resistant parasites at the genomic and phenotypic levels. We found that it took more time to develop resistance to TCMDC-143345 than to other drugs in clinical use and that there was no cross-resistance to these drugs, suggesting a new and unique mechanism. By whole-genome sequencing, we found two mutations in the gene encoding the L. donovani dynamin-1-like protein (LdoDLP1) that were fixed at the highest drug pressure. Through phylogenetic analysis, we identified LdoDLP1 as a family member of the dynamin-related proteins, a group of proteins that impacts the shapes of biological membranes by mediating fusion and fission events, with a putative role in mitochondrial fission. We found that L. donovani lines genetically engineered to harbor the two identified LdoDLP1 mutations were resistant to TCMDC-143345 and displayed altered mitochondrial properties. By homology modeling, we showed how the two LdoDLP1 mutations may influence protein structure and function. Taken together, our data reveal a clear involvement of LdoDLP1 in the adaptation/reduced susceptibility of L. donovani to TCMDC-143345. IMPORTANCE Humans and their pathogens are continuously locked in a molecular arms race during which the eventual emergence of pathogen drug resistance (DR) seems inevitable. For neglected tropical diseases (NTDs), DR is generally studied retrospectively once it has already been established in clinical settings. We previously recommended to keep one step ahead in the host-pathogen arms race and implement prospective DR studies in the R&D pipeline, a common practice for many infectious diseases but not for NTDs. Here, using Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GSK Leishbox to treat VL, as paradigms, we experimentally selected resistance to the compound and proceeded to genomic and phenotypic characterization of DR parasites. The results gathered in the present study suggest a new DR mechanism involving the L. donovani dynamin-1-like protein (LdoDLP1) and demonstrate the practical relevance of prospective DR studies.

Mapping Rora expression in resting and activated CD4+ T cells.

The transcription factor Rora has been shown to be important for the development of ILC2 and the regulation of ILC3, macrophages and Treg cells. Here we investigate the role of Rora across CD4+ T cells in general, but with an emphasis on Th2 cells, both in vitro as well as in the context of several in vivo type 2 infection models. We dissect the function of Rora using overexpression and a CD4-conditional Rora-knockout mouse, as well as a RORA-reporter mouse. We establish the importance of Rora in CD4+ T cells for controlling lung inflammation induced by Nippostrongylus brasiliensis infection, and have measured the effect on downstream genes using RNA-seq. Using a systematic stimulation screen of CD4+ T cells, coupled with RNA-seq, we identify upstream regulators of Rora, most importantly IL-33 and CCL7. Our data suggest that Rora is a negative regulator of the immune system, possibly through several downstream pathways, and is under control of the local microenvironment.

Individual-level variations in malaria susceptibility and acquisition of clinical protection.

After decades of research, our understanding of when and why individuals infected with Plasmodium falciparum develop clinical malaria is still limited. Correlates of immune protection are often sought through prospective cohort studies, where measured host factors are correlated against the incidence of clinical disease over a set period of time. However, robustly inferring individual-level protection from these population-level findings has proved difficult due to small effect sizes and high levels of variance underlying such data. In order to better understand the nature of these inter-individual variations, we analysed the long-term malaria epidemiology of children ≤12 years old growing up under seasonal exposure to the parasite in the sub-location of Junju, Kenya. Despite the cohort's limited geographic expanse (ca. 3km x 10km), our data reveal a high degree of spatial and temporal variability in malaria prevalence and incidence rates, causing individuals to experience varying levels of exposure to the parasite at different times during their life. Analysing individual-level infection histories further reveal an unexpectedly high variability in the rate at which children experience clinical malaria episodes. Besides exposure to the parasite, measured as disease prevalence in the surrounding area, we find that the birth time of year has an independent effect on the individual's risk of experiencing a clinical episode. Furthermore, our analyses reveal that those children with a history of an above average number of episodes are more likely to experience further episodes during the upcoming transmission season. These findings are indicative of phenotypic differences in the rates by which children acquire clinical protection to malaria and offer important insights into the natural variability underlying malaria epidemiology.

Complete representation of a tapeworm genome reveals chromosomes capped by centromeres, necessitating a dual role in segregation and protection.

BACKGROUND: Chromosome-level assemblies are indispensable for accurate gene prediction, synteny assessment, and understanding higher-order genome architecture. Reference and draft genomes of key helminth species have been published, but little is yet known about the biology of their chromosomes. Here, we present the complete genome of the tapeworm Hymenolepis microstoma, providing a reference quality, end-to-end assembly that represents the first fully assembled genome of a spiralian/lophotrochozoan, revealing new insights into chromosome evolution. RESULTS: Long-read sequencing and optical mapping data were added to previous short-read data enabling complete re-assembly into six chromosomes, consistent with karyology. Small genome size (169 Mb) and lack of haploid variation (1 SNP/3.2 Mb) contributed to exceptionally high contiguity with only 85 gaps remaining in regions of low complexity sequence. Resolution of repeat regions reveals novel gene expansions, micro-exon genes, and spliced leader trans-splicing, and illuminates the landscape of transposable elements, explaining observed length differences in sister chromatids. Syntenic comparison with other parasitic flatworms shows conserved ancestral linkage groups indicating that the H. microstoma karyotype evolved through fusion events. Strikingly, the assembly reveals that the chromosomes terminate in centromeric arrays, indicating that these motifs play a role not only in segregation, but also in protecting the linear integrity and full lengths of chromosomes. CONCLUSIONS: Despite strong conservation of canonical telomeres, our results show that they can be substituted by more complex, species-specific sequences, as represented by centromeres. The assembly provides a robust platform for investigations that require complete genome representation.

GeneDB and Wikidata.

Publishing authoritative genomic annotation data, keeping it up to date, linking it to related information, and allowing community annotation is difficult and hard to support with limited resources. Here, we show how importing GeneDB annotation data into Wikidata allows for leveraging existing resources, integrating volunteer and scientific communities, and enriching the original information.

Identification and expression profiling of microRNAs in Hymenolepis.

Tapeworms (cestodes) of the genus Hymenolepis are the causative agents of hymenolepiasis, a neglected zoonotic disease. Hymenolepis nana is the most prevalent human tapeworm, especially affecting children. The genomes of Hymenolepis microstoma and H. nana have been recently sequenced and assembled. MicroRNAs (miRNAs), a class of small non-coding RNAs, are principle regulators of gene expression at the post-transcriptional level and are involved in many different biological processes. In previous work, we experimentally identified miRNA genes in the cestodes Echinococcus, Taenia and Mesocestoides. However, current knowledge about miRNAs in Hymenolepis is limited. In this work we described for the first known time the expression profile of the miRNA complement in H. microstoma, and discovered miRNAs in H. nana. We found a reduced complement of 37 evolutionarily conserved miRNAs, putatively reflecting their low morphological complexity and parasitic lifestyle. We found high expression of a few miRNAs in the larval stage of H. microstoma that are conserved in other cestodes, suggesting that these miRNAs may have important roles in development, survival and for host-parasite interplay. We performed a comparative analysis of the identified miRNAs across the Cestoda and showed that most of the miRNAs in Hymenolepis are located in intergenic regions, implying that they are independently transcribed. We found a Hymenolepis-specific cluster composed of three members of the mir-36 family. Also, we found that one of the neighboring genes of mir-10 was a Hox gene as in most bilaterial species. This study provides a valuable resource for further experimental research in cestode biology that might lead to improved detection and control of these neglected parasites. The comprehensive identification and expression analysis of Hymenolepis miRNAs can help to identify novel biomarkers for diagnosis and/or novel therapeutic targets for the control of hymenolepiasis.

Evolutionary analysis of the most polymorphic gene family in falciparum malaria.

The var gene family of the human malaria parasite Plasmodium falciparum encode proteins that are crucial determinants of both pathogenesis and immune evasion and are highly polymorphic. Here we have assembled nearly complete var gene repertoires from 2398 field isolates and analysed a normalised set of 714 from across 12 countries. This therefore represents the first large scale attempt to catalogue the worldwide distribution of var gene sequences We confirm the extreme polymorphism of this gene family but also demonstrate an unexpected level of sequence sharing both within and between continents. We show that this is likely due to both the remnants of selective sweeps as well as a worrying degree of recent gene flow across continents with implications for the spread of drug resistance. We also address the evolution of the var repertoire with respect to the ancestral genes within the Laverania and show that diversity generated by recombination is concentrated in a number of hotspots. An analysis of the subdomain structure indicates that some existing definitions may need to be revised From the analysis of this data, we can now understand the way in which the family has evolved and how the diversity is continuously being generated. Finally, we demonstrate that because the genes are distributed across the genome, sequence sharing between genotypes acts as a useful population genetic marker.

Genome-wide transcriptome profiling and spatial expression analyses identify signals and switches of development in tapeworms.

BACKGROUND: Tapeworms are agents of neglected tropical diseases responsible for significant health problems and economic loss. They also exhibit adaptations to a parasitic lifestyle that confound comparisons of their development with other animals. Identifying the genetic factors regulating their complex ontogeny is essential to understanding unique aspects of their biology and for advancing novel therapeutics. Here we use RNA sequencing to identify up-regulated signalling components, transcription factors and post-transcriptional/translational regulators (genes of interest, GOI) in the transcriptomes of Larvae and different regions of segmented worms in the tapeworm Hymenolepis microstoma and combine this with spatial gene expression analyses of a selection of genes. RESULTS: RNA-seq reads collectively mapped to 90% of the > 12,000 gene models in the H. microstoma v.2 genome assembly, demonstrating that the transcriptome profiles captured a high percentage of predicted genes. Contrasts made between the transcriptomes of Larvae and whole, adult worms, and between the Scolex-Neck, mature strobila and gravid strobila, resulted in 4.5-30% of the genes determined to be differentially expressed. Among these, we identified 190 unique GOI up-regulated in one or more contrasts, including a large range of zinc finger, homeobox and other transcription factors, components of Wnt, Notch, Hedgehog and TGF-ß/BMP signalling, and post-transcriptional regulators (e.g. Boule, Pumilio). Heatmap clusterings based on overall expression and on select groups of genes representing 'signals' and 'switches' showed that expression in the Scolex-Neck region is more similar to that of Larvae than to the mature or gravid regions of the adult worm, which was further reflected in large overlap of up-regulated GOI. CONCLUSIONS: Spatial expression analyses in Larvae and adult worms corroborated inferences made from quantitative RNA-seq data and in most cases indicated consistency with canonical roles of the genes in other animals, including free-living flatworms. Recapitulation of developmental factors up-regulated during larval metamorphosis suggests that strobilar growth involves many of the same underlying gene regulatory networks despite the significant disparity in developmental outcomes. The majority of genes identified were investigated in tapeworms for the first time, setting the stage for advancing our understanding of developmental genetics in an important group of flatworm parasites.

Cyclin-dependent kinase 12 is a drug target for visceral leishmaniasis.

Visceral leishmaniasis causes considerable mortality and morbidity in many parts of the world. There is an urgent need for the development of new, effective treatments for this disease. Here we describe the development of an anti-leishmanial drug-like chemical series based on a pyrazolopyrimidine scaffold. The leading compound from this series (7, DDD853651/GSK3186899) is efficacious in a mouse model of visceral leishmaniasis, has suitable physicochemical, pharmacokinetic and toxicological properties for further development, and has been declared a preclinical candidate. Detailed mode-of-action studies indicate that compounds from this series act principally by inhibiting the parasite cdc-2-related kinase 12 (CRK12), thus defining a druggable target for visceral leishmaniasis.

Long read assemblies of geographically dispersed Plasmodium falciparum isolates reveal highly structured subtelomeres.

Background: Although thousands of clinical isolates of Plasmodium falciparum are being sequenced and analysed by short read technology, the data do not resolve the highly variable subtelomeric regions of the genomes that contain polymorphic gene families involved in immune evasion and pathogenesis. There is also no current standard definition of the boundaries of these variable subtelomeric regions. Methods: Using long-read sequence data (Pacific Biosciences SMRT technology), we assembled and annotated the genomes of 15 P. falciparum isolates, ten of which are newly cultured clinical isolates. We performed comparative analysis of the entire genome with particular emphasis on the subtelomeric regions and the internal var genes clusters.   Results: The nearly complete sequence of these 15 isolates has enabled us to define a highly conserved core genome, to delineate the boundaries of the subtelomeric regions, and to compare these across isolates. We found highly structured variable regions in the genome. Some exported gene families purportedly involved in release of merozoites show copy number variation. As an example of ongoing genome evolution, we found a novel CLAG gene in six isolates.  We also found a novel gene that was relatively enriched in the South East Asian isolates compared to those from Africa. Conclusions: These 15 manually curated new reference genome sequences with their nearly complete subtelomeric regions and fully assembled genes are an important new resource for the malaria research community. We report the overall conserved structure and pattern of important gene families and the more clearly defined subtelomeric regions.

Histone methylation changes are required for life cycle progression in the human parasite Schistosoma mansoni.

Epigenetic mechanisms and chromatin structure play an important role in development. Their impact is therefore expected to be strong in parasites with complex life cycles and multiple, strikingly different, developmental stages, i.e. developmental plasticity. Some studies have already described how the chromatin structure, through histone modifications, varies from a developmental stage to another in a few unicellular parasites. While H3K4me3 profiles remain relatively constant, H3K27 trimethylation and bivalent methylation show strong variation. Inhibitors (A366 and GSK343) of H3K27 histone methyltransferase activity in S. mansoni efficiently blocked miracidium to sporocyst transition indicating that H3K27 trimethylation is required for life cycle progression. As S. mansoni is a multicellular parasite that significantly affects both the health and economy of endemic areas, a better understanding of fluke developmental processes within the definitive host will likely highlight novel disease control strategies. Towards this goal, we also studied H4K20me1 in female cercariae and adults. In particular, we found that bivalent trimethylation of H3K4 and H3K27 at the transcription start site of genes is a landmark of the cercarial stage. In cercariae, H3K27me3 presence and strong enrichment in H4K20me1 over long regions (10-100 kb) is associated with development related genes. Here, we provide a broad overview of the chromatin structure of a metazoan parasite throughout its most important lifecycle stages. The five developmental stages studied here present distinct chromatin structures, indicating that histone methylation plays an important role during development. Hence, components of the histone methylation (and demethylation) machinery may provide suitable Schistosomiasis control targets.

An improved Plasmodium cynomolgi genome assembly reveals an unexpected methyltransferase gene expansion.

BACKGROUND: Plasmodium cynomolgi, a non-human primate malaria parasite species, has been an important model parasite since its discovery in 1907. Similarities in the biology of P. cynomolgi to the closely related, but less tractable, human malaria parasite P. vivax make it the model parasite of choice for liver biology and vaccine studies pertinent to P. vivax malaria. Molecular and genome-scale studies of P. cynomolgi have relied on the current reference genome sequence, which remains highly fragmented with 1,649 unassigned scaffolds and little representation of the subtelomeres.  Methods: Using long-read sequence data (Pacific Biosciences SMRT technology), we assembled and annotated a new reference genome sequence, PcyM, sourced from an Indian rhesus monkey. We compare the newly assembled genome sequence with those of several other Plasmodium species, including a re-annotated P. coatneyi assembly. RESULTS: The new PcyM genome assembly is of significantly higher quality than the existing reference, comprising only 56 pieces, no gaps and an improved average gene length. Detailed manual curation has ensured a comprehensive annotation of the genome with 6,632 genes, nearly 1,000 more than previously attributed to P. cynomolgi. The new assembly also has an improved representation of the subtelomeric regions, which account for nearly 40% of the sequence. Within the subtelomeres, we identified more than 1300 Plasmodium interspersed repeat ( pir) genes, as well as a striking expansion of 36 methyltransferase pseudogenes that originated from a single copy on chromosome 9. CONCLUSIONS: The manually curated PcyM reference genome sequence is an important new resource for the malaria research community. The high quality and contiguity of the data have enabled the discovery of a novel expansion of methyltransferase in the subtelomeres, and illustrates the new comparative genomics capabilities that are being unlocked by complete reference genomes.

Population genomics reveals the origin and asexual evolution of human infective trypanosomes.

Evolutionary theory predicts that the lack of recombination and chromosomal re-assortment in strictly asexual organisms results in homologous chromosomes irreversibly accumulating mutations and thus evolving independently of each other, a phenomenon termed the Meselson effect. We apply a population genomics approach to examine this effect in an important human pathogen, Trypanosoma brucei gambiense. We determine that T.b. gambiense is evolving strictly asexually and is derived from a single progenitor, which emerged within the last 10,000 years. We demonstrate the Meselson effect for the first time at the genome-wide level in any organism and show large regions of loss of heterozygosity, which we hypothesise to be a short-term compensatory mechanism for counteracting deleterious mutations. Our study sheds new light on the genomic and evolutionary consequences of strict asexuality, which this pathogen uses as it exploits a new biological niche, the human population.

What helminth genomes have taught us about parasite evolution.

The genomes of more than 20 helminths have now been sequenced. Here we perform a meta-analysis of all sequenced genomes of nematodes and Platyhelminthes, and attempt to address the question of what are the defining characteristics of helminth genomes. We find that parasitic worms lack systems for surface antigenic variation, instead maintaining infections using their surfaces as the first line of defence against the host immune system, with several expanded gene families of genes associated with the surface and tegument. Parasite excretory/secretory products evolve rapidly, and proteases even more so, with each parasite exhibiting unique modifications of its protease repertoire. Endoparasitic flatworms show striking losses of metabolic capabilities, not matched by nematodes. All helminths do however exhibit an overall reduction in auxiliary metabolism (biogenesis of co-factors and vitamins). Overall, the prevailing pattern is that there are few commonalities between the genomes of independently evolved parasitic worms, with each parasite having undergone specific adaptations for their particular niche.