Annotation and visualization of parasite, fungi and arthropod genomes with Companion.

As sequencing genomes has become increasingly popular, the need for annotation of the resulting assemblies is growing. Structural and functional annotation is still challenging as it includes finding the correct gene sequences, annotating other elements such as RNA and being able to submit those data to databases to share it with the community. Compared to de novo assembly where contiguous chromosomes are a sign of high quality, it is difficult to visualize and assess the quality of annotation. We developed the Companion web server to allow non-experts to annotate their genome using a reference-based method, enabling them to assess the output before submitting to public databases. In this update paper, we describe how we have included novel methods for gene finding and made the Companion server more efficient for annotation of genomes of up to 1 Gb in size. The reference set was increased to include genomes of interest for human and animal health from the fungi and arthropod kingdoms. We show that Companion outperforms existing comparable tools where closely related references are available.

A novel computational pipeline for var gene expression augments the discovery of changes in the Plasmodium falciparum transcriptome during transition from in vivo to short-term in vitro culture.

The pathogenesis of severe Plasmodium falciparum malaria involves cytoadhesive microvascular sequestration of infected erythrocytes, mediated by P. falciparum erythrocyte membrane protein 1 (PfEMP1). PfEMP1 variants are encoded by the highly polymorphic family of var genes, the sequences of which are largely unknown in clinical samples. Previously, we published new approaches for var gene profiling and classification of predicted binding phenotypes in clinical P. falciparum isolates (Wichers et al., 2021), which represented a major technical advance. Building on this, we report here a novel method for var gene assembly and multidimensional quantification from RNA-sequencing that outperforms the earlier approach of Wichers et al., 2021, on both laboratory and clinical isolates across a combination of metrics. Importantly, the tool can interrogate the var transcriptome in context with the rest of the transcriptome and can be applied to enhance our understanding of the role of var genes in malaria pathogenesis. We applied this new method to investigate changes in var gene expression through early transition of parasite isolates to in vitro culture, using paired sets of ex vivo samples from our previous study, cultured for up to three generations. In parallel, changes in non-polymorphic core gene expression were investigated. Modest but unpredictable var gene switching and convergence towards var2csa were observed in culture, along with differential expression of 19% of the core transcriptome between paired ex vivo and generation 1 samples. Our results cast doubt on the validity of the common practice of using short-term cultured parasites to make inferences about in vivo phenotype and behaviour.

Sterile protection against P. vivax malaria by repeated blood stage infection in the Aotus monkey model.

The malaria parasite Plasmodium vivax remains a major global public health challenge, and no vaccine is approved for use in humans. Here, we assessed whether P. vivax strain-transcendent immunity can be achieved by repeated infection in Aotus monkeys. Sterile immunity was achieved after two homologous infections, whereas subsequent heterologous challenge provided only partial protection. IgG levels based on P. vivax lysate ELISA and protein microarray increased with repeated infections and correlated with the level of homologous protection. Parasite transcriptional profiles provided no evidence of major antigenic switching upon homologous or heterologous challenge. However, we observed significant sequence diversity and transcriptional differences in the P. vivax core gene repertoire between the two strains used in the study, suggesting that partial protection upon heterologous challenge is due to molecular differences between strains rather than immune evasion by antigenic switching. Our study demonstrates that sterile immunity against P. vivax can be achieved by repeated homologous blood stage infection in Aotus monkeys, thus providing a benchmark to test the efficacy of candidate blood stage P. vivax malaria vaccines.

MalariaSED: a deep learning framework to decipher the regulatory contributions of noncoding variants in malaria parasites.

Malaria remains one of the deadliest infectious diseases. Transcriptional regulation effects of noncoding variants in this unusual genome of malaria parasites remain elusive. We developed a sequence-based, ab initio deep learning framework, MalariaSED, for predicting chromatin profiles in malaria parasites. The MalariaSED performance was validated by published ChIP-qPCR and TF motifs results. Applying MalariaSED to ~ 1.3 million variants shows that geographically differentiated noncoding variants are associated with parasite invasion and drug resistance. Further analysis reveals chromatin accessibility changes at Plasmodium falciparum rings are partly associated with artemisinin resistance. MalariaSED illuminates the potential functional roles of noncoding variants in malaria parasites.

Double-negative-2 B cells are the major synovial plasma cell precursor in rheumatoid arthritis.

B cells are key pathogenic drivers of chronic inflammation in rheumatoid arthritis (RA). There is limited understanding of the relationship between synovial B cell subsets and pathogenic antibody secreting cells (ASCs). This knowledge is crucial for the development of more targeted B-cell depleting therapies. While CD11c+ double-negative 2 (DN2) B cells have been suggested as an ASC precursor in lupus, to date there is no proven link between the two subsets in RA. We have used both single-cell gene expression and BCR sequencing to study synovial B cells from patients with established RA, in addition to flow cytometry of circulating B cells. To better understand the differentiation patterns within the diseased tissue, a combination of RNA-based trajectory inference and clonal lineage analysis of BCR relationships were used. Both forms of analysis indicated that DN2 B cells serve as a major precursors to synovial ASCs. This study advances our understanding of B cells in RA and reveals the origin of pathogenic ASCs in the RA synovium. Given the significant role of DN2 B cells as a progenitor to pathogenic B cells in RA, it is important to conduct additional research to investigate the origins of DN2 B cells in RA and explore their potential as therapeutic targets in place of the less specific pan-B cells depletion therapies currently in use.

From contigs towards chromosomes: automatic improvement of long read assemblies (ILRA).

Recent advances in long read technologies not only enable large consortia to aim to sequence all eukaryotes on Earth, but they also allow individual laboratories to sequence their species of interest with relatively low investment. Long read technologies embody the promise of overcoming scaffolding problems associated with repeats and low complexity sequences, but the number of contigs often far exceeds the number of chromosomes and they may contain many insertion and deletion errors around homopolymer tracts. To overcome these issues, we have implemented the ILRA pipeline to correct long read-based assemblies. Contigs are first reordered, renamed, merged, circularized, or filtered if erroneous or contaminated. Illumina short reads are used subsequently to correct homopolymer errors. We successfully tested our approach by improving the genome sequences of Homo sapiens, Trypanosoma brucei, and Leptosphaeria spp., and by generating four novel Plasmodium falciparum assemblies from field samples. We found that correcting homopolymer tracts reduced the number of genes incorrectly annotated as pseudogenes, but an iterative approach seems to be required to correct more sequencing errors. In summary, we describe and benchmark the performance of our new tool, which improved the quality of novel long read assemblies up to 1 Gbp. The pipeline is available at GitHub: https://github.com/ThomasDOtto/ILRA.

The exception that proves the rule: Virulence gene expression at the onset of Plasmodium falciparum blood stage infections.

Controlled human malaria infections (CHMI) are a valuable tool to study parasite gene expression in vivo under defined conditions. In previous studies, virulence gene expression was analyzed in samples from volunteers infected with the Plasmodium falciparum (Pf) NF54 isolate, which is of African origin. Here, we provide an in-depth investigation of parasite virulence gene expression in malaria-naïve European volunteers undergoing CHMI with the genetically distinct Pf 7G8 clone, originating in Brazil. Differential expression of var genes, encoding major virulence factors of Pf, PfEMP1s, was assessed in ex vivo parasite samples as well as in parasites from the in vitro cell bank culture that was used to generate the sporozoites (SPZ) for CHMI (Sanaria PfSPZ Challenge (7G8)). We report broad activation of mainly B-type subtelomeric located var genes at the onset of a 7G8 blood stage infection in naïve volunteers, mirroring the NF54 expression study and suggesting that the expression of virulence-associated genes is generally reset during transmission from the mosquito to the human host. However, in 7G8 parasites, we additionally detected a continuously expressed single C-type variant, Pf7G8_040025600, that was most highly expressed in both pre-mosquito cell bank and volunteer samples, suggesting that 7G8, unlike NF54, maintains expression of some previously expressed var variants during transmission. This suggests that in a new host, the parasite may preferentially express the variants that previously allowed successful infection and transmission. Trial registration: ClinicalTrials.gov - NCT02704533; 2018-004523-36.

Phenotypic Screens Identify Genetic Factors Associated with Gametocyte Development in the Human Malaria Parasite Plasmodium falciparum.

Transmission of the deadly malaria parasite Plasmodium falciparum from humans to mosquitoes is achieved by specialized intraerythrocytic sexual forms called gametocytes. Though the crucial regulatory mechanisms leading to gametocyte commitment have recently come to light, networks of genes that control sexual development remain to be elucidated. Here, we report a pooled-mutant screen to identify genes associated with gametocyte development in P. falciparum. Our results categorized genes that modulate gametocyte progression as hypoproducers or hyperproducers of gametocytes, and the in-depth analysis of individual clones confirmed phenotypes in sexual commitment rates and putative functions in gametocyte development. We present a new set of genes that have not been implicated in gametocytogenesis before and demonstrate the potential of forward genetic screens in isolating genes impacting parasite sexual biology, an exciting step toward the discovery of new antimalarials for a globally significant pathogen. IMPORTANCE Blocking human-to-vector transmission is an essential step toward malaria elimination. Gametocytes are solely responsible for achieving this transmission and represent an opportunity for therapeutic intervention. While these falciform-shaped parasite stages were first discovered in the 1880s, our understanding of the genetic determinants responsible for their formation and molecular mechanisms that drive their development is limited. In this work, we developed a scalable screening methodology with piggyBac mutants to identify genes that influence the development of gametocytes in the most lethal human malaria parasite, P. falciparum. By doing so, we lay the foundation for large-scale functional genomic studies specifically designed to address remaining questions about sexual commitment, maturation, and mosquito infection in P. falciparum. Such functional genetic screens will serve to expedite the identification of essential pathways and processes for the development of novel transmission-blocking agents.

Profiling the bloodstream form and procyclic form Trypanosoma brucei cell cycle using single-cell transcriptomics.

African trypanosomes proliferate as bloodstream forms (BSFs) and procyclic forms in the mammal and tsetse fly midgut, respectively. This allows them to colonise the host environment upon infection and ensure life cycle progression. Yet, understanding of the mechanisms that regulate and drive the cell replication cycle of these forms is limited. Using single-cell transcriptomics on unsynchronised cell populations, we have obtained high resolution cell cycle regulated (CCR) transcriptomes of both procyclic and slender BSF Trypanosoma brucei without prior cell sorting or synchronisation. Additionally, we describe an efficient freeze-thawing protocol that allows single-cell transcriptomic analysis of cryopreserved T. brucei. Computational reconstruction of the cell cycle using periodic pseudotime inference allowed the dynamic expression patterns of cycling genes to be profiled for both life cycle forms. Comparative analyses identify a core cycling transcriptome highly conserved between forms, as well as several genes where transcript levels dynamics are form specific. Comparing transcript expression patterns with protein abundance revealed that the majority of genes with periodic cycling transcript and protein levels exhibit a relative delay between peak transcript and protein expression. This work reveals novel detail of the CCR transcriptomes of both forms, which are available for further interrogation via an interactive webtool.

A novel Modulator of Ring Stage Translation (MRST) gene alters artemisinin sensitivity in Plasmodium falciparum.

The implementation of artemisinin (ART) combination therapies (ACTs) has greatly decreased deaths caused by Plasmodium falciparum malaria, but increasing ACT resistance in Southeast Asia and Africa could reverse this progress. Parasite population genetic studies have identified numerous genes, single-nucleotide polymorphisms (SNPs), and transcriptional signatures associated with altered artemisinin activity with SNPs in the Kelch13 (K13) gene being the most well-characterized artemisinin resistance marker. However, there is an increasing evidence that resistance to artemisinin in P. falciparum is not related only to K13 SNPs, prompting the need to characterize other novel genes that can alter ART responses in P. falciparum. In our previous analyses of P. falciparum piggyBac mutants, several genes of unknown function exhibited increased sensitivity to artemisinin that was similar to a mutant of K13. Further analysis of these genes and their gene co-expression networks indicated that the ART sensitivity cluster was functionally linked to DNA replication and repair, stress responses, and maintenance of homeostatic nuclear activity. In this study, we have characterized PF3D7_1136600, another member of the ART sensitivity cluster. Previously annotated as a conserved Plasmodium gene of unknown function, we now provide putative annotation of this gene as a Modulator of Ring Stage Translation (MRST). Our findings reveal that the mutagenesis of MRST affects gene expression of multiple translation-associated pathways during the early ring stage of asexual development via putative ribosome assembly and maturation activity, suggesting an essential role of MRST in protein biosynthesis and another novel mechanism of altering the parasite's ART drug response.IMPORTANCEPlasmodium falciparum malaria killed more than 600,000 people in 2021, though ACTs have been critical in reducing malaria mortality as a first-line treatment for infection. However, ACT resistance in Southeast Asia and emerging resistance in Africa are detrimental to this progress. Mutations to Kelch13 (K13) have been identified to confer increased artemisinin tolerance in field isolates, however, genes other than K13 are implicated in altering how the parasite responds to artemisinin prompts additional analysis. Therefore, in this study we have characterized a P. falciparum mutant clone with altered sensitivity to artemisinin and identified a novel gene (PF3D7_1136600) that is associated with alterations to parasite translational metabolism during critical timepoints for artemisinin drug response. Many genes of the P. falciparum genome remain unannotated, posing a challenge for drug-gene characterizations in the parasite. Therefore, through this study, we have putatively annotated PF3D7_1136600 as a novel MRST gene and have identified a potential link between MRST and parasite stress response mechanisms.

Chemogenomic Profiling of a Plasmodium falciparum Transposon Mutant Library Reveals Shared Effects of Dihydroartemisinin and Bortezomib on Lipid Metabolism and Exported Proteins.

The antimalarial activity of the frontline drug artemisinin involves generation of reactive oxygen species (ROS) leading to oxidative damage of parasite proteins. To achieve homeostasis and maintain protein quality control in the overwhelmed parasite, the ubiquitin-proteasome system kicks in. Even though molecular markers for artemisinin resistance like pfkelch13 have been identified, the intricate network of mechanisms driving resistance remains to be elucidated. Here, we report a forward genetic screening strategy that enables a broader identification of genetic factors responsible for altering sensitivity to dihydroartemisinin (DHA) and a proteasome inhibitor, bortezomib (BTZ). Using a library of isogenic piggyBac mutants in P. falciparum, we defined phenotype-genotype associations influencing drug responses and highlighted shared mechanisms between the two processes, which mainly included proteasome-mediated degradation and the lipid metabolism genes. Additional transcriptomic analysis of a DHA/BTZ-sensitive piggyBac mutant showed it is possible to find differences between the two response mechanisms on the specific components for regulation of the exportome. Our results provide further insight into the molecular mechanisms of antimalarial drug resistance. IMPORTANCE Malaria control is seriously threatened by the emergence and spread of Plasmodium falciparum resistance to the leading antimalarial, artemisinin. The potent killing activity of artemisinin results from oxidative damage unleashed by free heme activation released by hemoglobin digestion. Although the ubiquitin-proteasome system is considered critical for parasite survival of this toxicity, the diverse genetic changes linked to artemisinin resistance are complex and, so far, have not included the ubiquitin-proteasome system. In this study, we use a systematic forward genetic approach by screening a library of P. falciparum random piggyBac mutants to decipher the genetic factors driving malaria parasite responses to the oxidative stress caused by antimalarial drugs. This study compares phenotype-genotype associations influencing dihydroartemisinin responses with the proteasome inhibitor bortezomib to delineate the role of ubiquitin-proteasome system. Our study highlights shared and unique pathways from the complex array of molecular processes critical for P. falciparum survival resulting from the oxidative damage of artemisinin.

peaks2utr: a robust Python tool for the annotation of 3' UTRs.

SUMMARY: Annotation of nonmodel organisms is an open problem, especially the detection of untranslated regions (UTRs). Correct annotation of UTRs is crucial in transcriptomic analysis to accurately capture the expression of each gene yet is mostly overlooked in annotation pipelines. Here we present peaks2utr, an easy-to-use Python command line tool that uses the UTR enrichment of single-cell technologies, such as 10× Chromium, to accurately annotate 3' UTRs for a given canonical annotation. AVAILABILITY AND IMPLEMENTATION: peaks2utr is implemented in Python 3 (≥3.8). It is available via PyPI at https://pypi.org/project/peaks2utr and GitHub at https://github.com/haessar/peaks2utr. It is licensed under GNU GPLv3.

New technologies to study helminth development and host-parasite interactions.

How parasites develop and survive, and how they stimulate or modulate host immune responses are important in understanding disease pathology and for the design of new control strategies. Microarray analysis and bulk RNA sequencing have provided a wealth of data on gene expression as parasites develop through different life-cycle stages and on host cell responses to infection. These techniques have enabled gene expression in the whole organism or host tissue to be detailed, but do not take account of the heterogeneity between cells of different types or developmental stages, nor the spatial organisation of these cells. Single-cell RNA-seq (scRNA-seq) adds a new dimension to studying parasite biology and host immunity by enabling gene profiling at the individual cell level. Here we review the application of scRNA-seq to establish gene expression cell atlases for multicellular helminths and to explore the expansion and molecular profile of individual host cell types involved in parasite immunity and tissue repair. Studying host-parasite interactions in vivo is challenging and we conclude this review by briefly discussing the applications of organoids (stem-cell derived mini-tissues) to examine host-parasite interactions at the local level, and as a potential system to study parasite development in vitro. Organoid technology and its applications have developed rapidly, and the elegant studies performed to date support the use of organoids as an alternative in vitro system for research on helminth parasites.

Protein KIC5 is a novel regulator of artemisinin stress response in the malaria parasite Plasmodium falciparum.

Artemisinin combination therapies (ACTs) have led to a significant decrease in Plasmodium falciparum malaria mortality. This progress is now threatened by emerging artemisinin resistance (ART-R) linked originally in SE Asia to polymorphisms in the Kelch propeller protein (K13) and more recently to several other seemingly unrelated genetic mutations. To better understand the parasite response to ART, we are characterizing a P. falciparum mutant with altered sensitivity to ART that was created via piggyBac transposon mutagenesis. The transposon inserted near the putative transcription start site of a gene defined as a "Plasmodium-conserved gene of unknown function," now functionally linked to K13 as the Kelch13 Interacting Candidate 5 protein (KIC5). Phenotype analysis of the KIC5 mutant during intraerythrocytic asexual development identified transcriptional changes associated with DNA stress response and altered mitochondrial metabolism, linking dysregulation of the KIC5 gene to the parasite's ability to respond to ART exposure. Through characterization of the KIC5 transcriptome, we hypothesize that this gene may be essential under ART exposure to manage gene expression of the wild-type stress response at early ring stage, thereby providing a better understanding of the parasite's processes that can alter ART sensitivity.

Single cell and spatial transcriptomic analyses reveal microglia-plasma cell crosstalk in the brain during Trypanosoma brucei infection.

Human African trypanosomiasis, or sleeping sickness, is caused by the protozoan parasite Trypanosoma brucei and induces profound reactivity of glial cells and neuroinflammation when the parasites colonise the central nervous system. However, the transcriptional and functional responses of the brain to chronic T. brucei infection remain poorly understood. By integrating single cell and spatial transcriptomics of the mouse brain, we identify that glial responses triggered by infection are readily detected in the proximity to the circumventricular organs, including the lateral and 3rd ventricle. This coincides with the spatial localisation of both slender and stumpy forms of T. brucei. Furthermore, in silico predictions and functional validations led us to identify a previously unknown crosstalk between homeostatic microglia and Cd138+ plasma cells mediated by IL-10 and B cell activating factor (BAFF) signalling. This study provides important insights and resources to improve understanding of the molecular and cellular responses in the brain during infection with African trypanosomes.

Using topic modeling to detect cellular crosstalk in scRNA-seq.

Cell-cell interactions are vital for numerous biological processes including development, differentiation, and response to inflammation. Currently, most methods for studying interactions on scRNA-seq level are based on curated databases of ligands and receptors. While those methods are useful, they are limited to our current biological knowledge. Recent advances in single cell protocols have allowed for physically interacting cells to be captured, and as such we have the potential to study interactions in a complemantary way without relying on prior knowledge. We introduce a new method based on Latent Dirichlet Allocation (LDA) for detecting genes that change as a result of interaction. We apply our method to synthetic datasets to demonstrate its ability to detect genes that change in an interacting population compared to a reference population. Next, we apply our approach to two datasets of physically interacting cells to identify the genes that change as a result of interaction, examples include adhesion and co-stimulatory molecules which confirm physical interaction between cells. For each dataset we produce a ranking of genes that are changing in subpopulations of the interacting cells. In addition to the genes discussed in the original publications, we highlight further candidates for interaction in the top 100 and 300 ranked genes. Lastly, we apply our method to a dataset generated by a standard droplet-based protocol not designed to capture interacting cells, and discuss its suitability for analysing interactions. We present a method that streamlines detection of interactions and does not require prior clustering and generation of synthetic reference profiles to detect changes in expression.

Host cell maturation modulates parasite invasion and sexual differentiation in Plasmodium berghei.

Malaria remains a global health problem causing more than 400,000 deaths annually. Plasmodium parasites, the causative agents of malaria, replicate asexually in red blood cells (RBCs) of their vertebrate host, while a subset differentiates into sexual stages (gametocytes) for mosquito transmission. Parasite replication and gametocyte maturation in the erythropoietic niches of the bone marrow and spleen contribute to pathogenesis and drive transmission, but the mechanisms underlying this organ enrichment remain unknown. Here, we performed a comprehensive analysis of rodent P. berghei infection by flow cytometry and single-cell RNA sequencing. We identified CD71 as a host receptor for reticulocyte invasion and found that parasites metabolically adapt to the host cell environment. Transcriptional analysis and functional assays further revealed a nutrient-dependent tropism for gametocyte formation in reticulocytes. Together, we provide a thorough characterization of host-parasite interactions in erythropoietic niches and define host cell maturation state as the key driver of parasite adaptation.

Varia: a tool for prediction, analysis and visualisation of variable genes.

BACKGROUND: Parasites use polymorphic gene families to evade the immune system or interact with the host. Assessing the diversity and expression of such gene families in pathogens can inform on the repertoire or host interaction phenotypes of clinical relevance. However, obtaining the sequences and quantifying their expression is a challenge. In Plasmodium falciparum, the highly polymorphic var genes encode the major virulence protein, PfEMP1, which bind a range of human receptors through varying combinations of DBL and CIDR domains. Here we present a tool, Varia, to predict near full-length gene sequences and domain compositions of query genes from database genes sharing short sequence tags. Varia generates output through two complementary pipelines. Varia_VIP returns all putative gene sequences and domain compositions of the query gene from any partial sequence provided, thereby enabling experimental validation of specific genes of interest and detailed assessment of their putative domain structure. Varia_GEM accommodates rapid profiling of var gene expression in complex patient samples from DBLα expression sequence tags (EST), by computing a sample overall transcript profile stratified by PfEMP1 domain types. RESULTS: Varia_VIP was tested querying sequence tags from all DBL domain types using different search criteria. On average 92% of query tags had one or more 99% identical database hits, resulting in the full-length query gene sequence being identified (> 99% identical DNA > 80% of query gene) among the five most prominent database hits, for ~ 33% of the query genes. Optimized Varia_GEM settings allowed correct prediction of > 90% of domains placed among the four most N-terminal domains, including the DBLα domain, and > 70% of C-terminal domains. With this accuracy, N-terminal domains could be predicted for > 80% of queries, whereas prediction rates of C-terminal domains dropped with the distance from the DBLα from 70 to 40%. CONCLUSION: Prediction of var sequence and domain composition is possible from short sequence tags. Varia can be used to guide experimental validation of PfEMP1 sequences of interest and conduct high-throughput analysis of var type expression in patient samples.