VEuPathDB: the eukaryotic pathogen, vector and host bioinformatics resource center

Malate dehydrogenase in parasitic protozoans: roles in metabolism and potential therapeutic applications

The role of malate dehydrogenase (MDH) in the metabolism of various medically significant protozoan parasites is reviewed. MDH is an NADH-dependent oxidoreductase that catalyzes interconversion between oxaloacetate and malate, provides metabolic intermediates for both catabolic and anabolic pathways, and can contribute to NAD+/NADH balance in multiple cellular compartments. MDH is present in nearly all organisms; isoforms of MDH from apicomplexans (Plasmodium falciparum, Toxoplasma gondii, Cryptosporidium spp.), trypanosomatids (Trypanosoma brucei, T. cruzi) and anaerobic protozoans (Trichomonas vaginalis, Giardia duodenalis) are presented here. Many parasitic species have complex life cycles and depend on the environment of their hosts for carbon sources and other nutrients. Metabolic plasticity is crucial to parasite transition between host environments; thus, the regulation of metabolic processes is an important area to explore for therapeutic intervention. Common themes in protozoan parasite metabolism include emphasis on glycolytic catabolism, substrate-level phosphorylation, non-traditional uses of common pathways like tricarboxylic acid cycle and adapted or reduced mitochondria-like organelles. We describe the roles of MDH isoforms in these pathways, discuss unusual structural or functional features of these isoforms relevant to activity or drug targeting, and review current studies exploring the therapeutic potential of MDH and related genes. These studies show that MDH activity has important roles in many metabolic pathways, and thus in the metabolic transitions of protozoan parasites needed for success as pathogens.

Identification of the drug/metabolite transporter 1 as a marker of quinine resistance in a NF54×Cam3.II P. falciparum genetic cross

The genetic basis of Plasmodium falciparum resistance to quinine (QN), a drug used to treat severe malaria, has long been enigmatic. To gain further insight, we used FRG-NOD human liver-chimeric mice to conduct a P. falciparum genetic cross between QN-sensitive and QN-resistant parasites, which also differ in their susceptibility to chloroquine (CQ). By applying different selective conditions to progeny pools prior to cloning, we recovered 120 unique recombinant progeny. These progeny were subjected to drug profiling and QTL analyses with QN, CQ, and monodesethyl-CQ (md-CQ, the active metabolite of CQ), which revealed predominant peaks on chromosomes 7 and 12, consistent with a multifactorial mechanism of resistance. A shared chromosome 12 region mapped to resistance to all three antimalarials and was preferentially co-inherited with pfcrt. We identified an ATP-dependent zinc metalloprotease (FtsH1) as one of the top candidates and observed using CRISPR/Cas9 SNP-edited lines that ftsh1 is a potential mediator of QN resistance and a modulator of md-CQ resistance. As expected, CQ and md-CQ resistance mapped to a chromosome 7 region harboring pfcrt. However, for QN, high-grade resistance mapped to a chromosome 7 peak centered 295kb downstream of pfcrt. We identified the drug/metabolite transporter 1 (DMT1) as the top candidate due to its structural similarity to PfCRT and proximity to the peak. Deleting DMT1 in QN-resistant Cam3.II parasites significantly sensitized the parasite to QN but not to the other drugs tested, suggesting that DMT1 mediates QN response specifically. We localized DMT1 to structures associated with vesicular trafficking, as well as the parasitophorous vacuolar membrane, lipid bodies, and the digestive vacuole. We also observed that mutant DMT1 transports more QN than the wild-type isoform in vitro. Our study demonstrates that DMT1 is a novel marker of QN resistance and a new chromosome 12 locus associates with CQ and QN response, with ftsh1 is a potential candidate, suggesting these genes should be genotyped in surveillance and clinical settings.

Identification of potent and orally efficacious phosphodiesterase inhibitors in Cryptosporidium parvum-infected immunocompromised male mice

Cryptosporidium parvum and C. hominis are parasites that cause life-threatening diarrhea in children and immunocompromised people. There is only one approved treatment that is modestly effective for children and ineffective for AIDS patients. Here, screening 278 compounds from the Merck KGaA, Darmstadt, Germany collection and accelerated follow-up enabled by prior investigation of the compounds identifies a series of pyrazolopyrimidine human phosphodiesterase (PDE)-V (hPDE-V) inhibitors with potent anticryptosporidial activity and efficacy following oral administration in C. parvum-infected male mice. The lead compounds affect parasite host cell egress, inhibit both C. parvum and C. hominis, work rapidly, and have minimal off-target effects in a safety screening panel. Interestingly, the hPDE-V inhibitors sildenafil and the 4-aminoquinoline compound 7a do not affect Cryptosporidium. C. parvum expresses one PDE (CpPDE1) continuously during asexual growth, the inhibited life stage. According to homology modeling and docking, the lead compounds interact with CpPDE1. Bulkier amino acids (Val900 and His884) in the CpPDE1 active site replace alanines in hPDE-V and block sildenafil binding. Supporting this, sildenafil kills a CRISPR-engineered Cryptosporidium CpPDE1 V900A mutant. The CpPDE1 mutation also alters parasite susceptibility to pyrazolopyrimidines. CpPDE1 is therefore a validated pyrazolopyrimidine molecular target to exploit for target-based optimization for improved anticryptosporidial development.

Cellulose binding and the timing of expression influence protein targeting to the double-layered cyst wall of Acanthamoeba

The cyst wall of the eye pathogen Acanthamoeba castellanii contains cellulose and has ectocyst and endocyst layers connected by conical ostioles. Cyst walls contain families of lectins that localize to the ectocyst layer (Jonah) or the endocyst layer and ostioles (Luke and Leo). How lectins and an abundant laccase bind cellulose and why proteins go to locations in the wall are not known and are the focus of the studies here. Structural predictions identified β-jelly-roll folds (BJRFs) of Luke and sets of four disulfide knots (4DKs) of Leo, each of which contains linear arrays of aromatic amino acids, also present in carbohydrate-binding modules of bacterial and plant endocellulases. Ala mutations showed that these aromatics are necessary for cellulose binding and proper localization of Luke and Leo in the Acanthamoeba cyst wall. BJRFs of Luke, 4DKs of Leo, a single β-helical fold (BHF) of Jonah, and a copper oxidase domain of the laccase each bind to glycopolymers in both layers of deproteinated cyst walls. Promoter swaps showed that ectocyst localization does not just correlate with but is caused by early encystation-specific expression, while localization in the endocyst layer and ostioles is caused by later expression. Evolutionary studies showed distinct modes of assembly of duplicated domains in Luke, Leo, and Jonah lectins and suggested Jonah BHFs originated from bacteria, Luke BJRFs share common ancestry with slime molds, while 4DKs of Leo are unique to Acanthamoeba.IMPORTANCEAcanthamoebae is the only human parasite with cellulose in its cyst wall and conical ostioles that connect its inner and outer layers. Cyst walls are important virulence factors because they make Acanthamoebae resistant to surface disinfectants, hand sanitizers, contact lens sterilizers, and antibiotics applied to the eye. The goal here was to understand better how proteins are targeted to specific locations in the cyst wall. To this end, we identified three new proteins in the outer layer of the cyst wall, which may be targets for diagnostic antibodies in corneal scrapings. We used structural predictions and mutated proteins to show linear arrays of aromatic amino acids of two unrelated wall proteins are necessary for binding cellulose and proper wall localization. We showed early expression during encystation causes proteins to localize to the outer layer, while later expression causes proteins to localize to the inner layer and the ostioles.

DNA sequence and chromatin differentiate sequence-specific transcription factor binding in the human malaria parasite Plasmodium falciparum

Development of the malaria parasite, Plasmodium falciparum, is regulated by a limited number of sequence-specific transcription factors (TFs). However, the mechanisms by which these TFs recognize genome-wide binding sites is largely unknown. To address TF specificity, we investigated the binding of two TF subsets that either bind CACACA or GTGCAC DNA sequence motifs and further characterized two additional ApiAP2 TFs, PfAP2-G and PfAP2-EXP, which bind unique DNA motifs (GTAC and TGCATGCA). We also interrogated the impact of DNA sequence and chromatin context on P. falciparum TF binding by integrating high-throughput in vitro and in vivo binding assays, DNA shape predictions, epigenetic post-translational modifications, and chromatin accessibility. We found that DNA sequence context minimally impacts binding site selection for paralogous CACACA-binding TFs, while chromatin accessibility, epigenetic patterns, co-factor recruitment, and dimerization correlate with differential binding. In contrast, GTGCAC-binding TFs prefer different DNA sequence context in addition to chromatin dynamics. Finally, we determined that TFs that preferentially bind divergent DNA motifs may bind overlapping genomic regions due to low-affinity binding to other sequence motifs. Our results demonstrate that TF binding site selection relies on a combination of DNA sequence and chromatin features, thereby contributing to the complexity of P. falciparum gene regulatory mechanisms.

Mosquito salivary apyrase regulates blood meal hemostasis and facilitates malaria parasite transmission

The evolution of hematophagy involves a series of adaptations that allow blood-feeding insects to access and consume blood efficiently while managing and circumventing the host's hemostatic and immune responses. Mosquito, and other insects, utilize salivary proteins to regulate these responses at the bite site during and after blood feeding. We investigated the function of Anopheles gambiae salivary apyrase (AgApyrase) in regulating hemostasis in the mosquito blood meal and in Plasmodium transmission. Our results demonstrate that salivary apyrase, a known inhibitor of platelet aggregation, interacts with and activates tissue plasminogen activator, facilitating the conversion of plasminogen to plasmin, a human protease that degrades fibrin and facilitates Plasmodium transmission. We show that mosquitoes ingest a substantial amount of apyrase during blood feeding, which reduces coagulation in the blood meal by enhancing fibrin degradation and inhibiting platelet aggregation. AgApyrase significantly enhanced Plasmodium infection in the mosquito midgut, whereas AgApyrase immunization inhibited Plasmodium mosquito infection and sporozoite transmission. This study highlights a pivotal role for mosquito salivary apyrase for regulation of hemostasis in the mosquito blood meal and for Plasmodium transmission to mosquitoes and to the mammalian host, underscoring the potential for strategies to prevent malaria transmission.

Protein disulfide isomerase PDI8 is indispensable for parasite growth and associated with secretory protein processing in Toxoplasma gondii

Protein disulfide isomerase, containing thioredoxin (Trx) domains, serves as a vital enzyme responsible for oxidative protein folding (the formation, reduction, and isomerization of disulfide bonds in newly synthesized proteins) in the endoplasmic reticulum (ER). However, the role of ER-localized PDI proteins in parasite growth and their interaction with secretory proteins remain poorly understood. In this study, we identified two ER-localized PDI proteins, TgPDI8 and TgPDI6, in Toxoplasma gondii. Conditional knockdown of TgPDI8 resulted in a significant reduction in intracellular proliferation and invasion abilities, leading to a complete block in plaque formation on human foreskin fibroblast monolayers, whereas parasites lacking TgPDI6 did not exhibit any apparent fitness defects. The complementation of TgPDI8 with mutant variants highlighted the critical role of the CXXC active site cysteines within its Trx domains for its enzymatic activity. By utilizing TurboID-based proximity labeling, we uncovered a close association between PDI proteins and canonical secretory proteins. Furthermore, parasites lacking TgPDI8 showed a significant reduction in the expression of secretory proteins, especially those from micronemes and dense granules. In summary, our study elucidates the roles of TgPDI8 and sets the stage for future drug discovery studies.

IMPORTANCE

Apicomplexans, a phylum of intracellular parasites, encompass various zoonotic pathogens, including Plasmodium, Cryptosporidium, Toxoplasma, and Babesia, causing a significant economic burden on human populations. These parasites exhibit hypersensitivity to disruptions in endoplasmic reticulum (ER) redox homeostasis, necessitating the presence of ER-localized thioredoxin (Trx) superfamily proteins, particularly protein disulfide isomerase (PDI), for proper oxidative folding. However, the functional characteristics of ER-localized PDI proteins in Toxoplasma gondii remain largely unexplored. In this study, we identified two ER-localized proteins, namely, TgPDI8 and TgPDI6, and demonstrated the indispensable role of TgPDI8 in parasite survival. Through a comprehensive multi-omics analysis, we elucidated the crucial role of TgPDI8 in the processing of secretory proteins in T. gondii. Additionally, we introduced a novel ER-anchored TurboID method to label and identify canonical secretory proteins in T. gondii. This research opens up new avenues for understanding oxidative folding and the secretory pathway in apicomplexan parasites, laying the groundwork for future advancements in antiparasitic drug development.

ATM1, an essential conserved transporter in Apicomplexa, bridges mitochondrial and cytosolic [Fe-S] biogenesis

The Apicomplexa phylum encompasses numerous obligate intracellular parasites, some associated with severe implications for human health, including Plasmodium, Cryptosporidium, and Toxoplasma gondii. The iron-sulfur cluster [Fe-S] biogenesis ISC pathway, localized within the mitochondrion or mitosome of these parasites, is vital for parasite survival and development. Previous work on T. gondii and Plasmodium falciparum provided insights into the mechanisms of [Fe-S] biogenesis within this phylum, while the transporter linking mitochondria-generated [Fe-S] with the cytosolic [Fe-S] assembly (CIA) pathway remained elusive. This critical step is catalyzed by a well-conserved ABC transporter, termed ATM1 in yeast, ATM3 in plants and ABCB7 in mammals. Here, we identify and characterize this transporter in two clinically relevant Apicomplexa. We demonstrate that depletion of TgATM1 does not specifically impair mitochondrial metabolism. Instead, proteomic analyses reveal that TgATM1 expression levels inversely correlate with the abundance of proteins that participate in the transfer of [Fe-S] to cytosolic proteins at the outer mitochondrial membrane. Further insights into the role of TgATM1 are gained through functional complementation with the well-characterized yeast homolog. Biochemical characterization of PfATM1 confirms its role as a functional ABC transporter, modulated by oxidized glutathione (GSSG) and [4Fe-4S].

Cytokinetic abscission in Toxoplasma gondii is governed by protein phosphatase 2A and the daughter cell scaffold complex

Cytokinetic abscission marks the final stage of cell division, during which the daughter cells physically separate through the generation of new barriers, such as the plasma membrane or cell wall. While the contractile ring plays a central role during cytokinesis in bacteria, fungi and animal cells, the process diverges in Apicomplexa. In Toxoplasma gondii, two daughter cells are formed within the mother cell by endodyogeny. The mechanism by which the progeny cells acquire their plasma membrane during the disassembly of the mother cell, allowing daughter cells to emerge, remains unknown. Here we identify and characterize five T. gondii proteins, including three protein phosphatase 2A subunits, which exhibit a distinct and dynamic localization pattern during parasite division. Individual downregulation of these proteins prevents the accumulation of plasma membrane at the division plane, preventing the completion of cellular abscission. Remarkably, the absence of cytokinetic abscission does not hinder the completion of subsequent division cycles. The resulting progeny are able to egress from the infected cells but fail to glide and invade, except in cases of conjoined twin parasites.

Annotation of Functions of Sequences of Concern and Its Relevance to the New Biosecurity Regulatory Framework in the United States

Introduction

Recent regulations from United States Government agencies reshape the screening of synthetic nucleic acids. These take a step away from categorizing hazard on the basis of "bad" taxa and invoke the function of the sequence in pathogenesis or intoxication. Ascertaining functions related to pathogenesis and distinguishing these from other molecular abilities that are unproblematic is not simple. Some have suggested that this information can be readily obtained from existing databases of pathogens.

Objectives

We evaluate how virulence factors are described in current databases of pathogens and their adequacy for biothreat data science. We discuss limitations of how virulence factors have been conceived and propose using the sequence of concern (SoC) term to distinguish sequences with biothreat from those without. We discuss ways in which databases of SoCs might be implemented for research and regulatory purposes. We describe ongoing work improving functional descriptions of SoCs.

Methods

We assess the adequacy of descriptions of virulence factors in pathogen databases following extensive engagement with the literature in microbial pathogenesis.

Results/Conclusions

Descriptions of virulence factors in pathogen databases are inadequate for understanding biothreats. Many are not biothreats and would not be concerning if transferred to another pathogen. New gene ontology terms have been authored, and those specific to pathogenic viral processes are being generalized to make them relevant to other pathogenic taxa. This allows better understanding by humans and better recognition by machines. A database of annotated functions of SoCs could benefit the evolving biosecurity regulatory framework in the United States.

Phosphoproteomic analysis of the response to DNA damage in Trypanosoma brucei

Damage to the genetic material of the cell poses a universal threat to all forms of life. The DNA damage response is a coordinated cellular response to a DNA break, key to which is the phosphorylation signaling cascade. Identifying which proteins are phosphorylated is therefore crucial to understanding the mechanisms that underlie it. We have used stable isotopic labeling of amino acids in cell culture-based quantitative phosphoproteomics to profile changes in phosphorylation site abundance following double stranded DNA breaks, at two distinct loci in the genome of the single cell eukaryote Trypanosoma brucei. Here, we report on the T. brucei phosphoproteome following a single double-strand break at either a chromosome internal or subtelomeric locus, specifically the bloodstream form expression site. We detected >6500 phosphorylation sites, of which 211 form a core set of double-strand break responsive phosphorylation sites. Along with phosphorylation of canonical DNA damage factors, we have identified two novel phosphorylation events on histone H2A and found that in response to a chromosome internal break, proteins are predominantly phosphorylated, while a greater proportion of proteins dephosphorylated following a DNA break at a subtelomeric bloodstream form expression site. Our data represent the first DNA damage phosphoproteome and provides novel insights into repair at distinct chromosomal contexts in T. brucei.

DNA damage drives antigen diversification through mosaic Variant Surface Glycoprotein (VSG) formation in Trypanosoma brucei

Antigenic variation, using large genomic repertoires of antigen-encoding genes, allows pathogens to evade host antibody. Many pathogens, including the African trypanosome Trypanosoma brucei, extend their antigenic repertoire through genomic diversification. While evidence suggests that T. brucei depends on the generation of new variant surface glycoprotein (VSG) genes to maintain a chronic infection, a lack of experimentally tractable tools for studying this process has obscured its underlying mechanisms. Here, we present a highly sensitive targeted sequencing approach for measuring VSG diversification. Using this method, we demonstrate that a Cas9-induced DNA double-strand break within the VSG coding sequence can induce VSG recombination with patterns identical to those observed during infection. These newly generated VSGs are antigenically distinct from parental clones and thus capable of facilitating immune evasion. Together, these results provide insight into the mechanisms of VSG diversification and an experimental framework for studying the evolution of antigen repertoires in pathogenic microbes.

The Entamoeba histolytica Vps26 (EhVps26) retromeric protein is involved in phagocytosis: Bioinformatic and experimental approaches

The retromer is a cellular structure that recruits and recycles proteins inside the cell. In mammalian and yeast, the retromer components have been widely studied, but very little in parasites. In yeast, it is formed by a SNX-BAR membrane remodeling heterodimer and the cargo selecting complex (CSC), composed by three proteins. One of them, the Vps26 protein, possesses a flexible and intrinsically disordered region (IDR), that facilitates interactions with other proteins and contributes to the retromer binding to the endosomal membrane. In Entamoeba histolytica, the protozoan parasite responsible for human amoebiasis, the retromer actively participates during the high mobility and phagocytosis of trophozoites, but the molecular details in these events, are almost unknown. Here, we studied the EhVps26 role in phagocytosis. Bioinformatic analyses of EhVps26 revealed a typical arrestin folding structure of the protein, and a long and charged IDR, as described in other systems. EhVps26 molecular dynamics simulations (MDS) allowed us to predict binding pockets for EhVps35, EhSNX3, and a PX domain-containing protein; these pockets were disorganized in a EhVps26 truncated version lacking the IDR. The AlphaFold2 software predicted the interaction of EhVps26 with EhVps35, EhVps29 and EhSNX3, in a model similar to the reported mammalian crystals. By confocal and transmission electron microscopy, EhVps26 was found in the trophozoites plasma membrane, cytosol, endosomes, and Golgi-like apparatus. During phagocytosis, it followed the erythrocytes pathway, probably participating in cargoes selection and recycling. Ehvps26 gene knocking down evidenced that the EhVps26 protein is necessary for efficient phagocytosis.

Novel metagenomics analysis of stony coral tissue loss disease

Stony coral tissue loss disease (SCTLD) has devastated coral reefs off the coast of Florida and continues to spread throughout the Caribbean. Although a number of bacterial taxa have consistently been associated with SCTLD, no pathogen has been definitively implicated in the etiology of SCTLD. Previous studies have predominantly focused on the prokaryotic community through 16S rRNA sequencing of healthy and affected tissues. Here, we provide a different analytical approach by applying a bioinformatics pipeline to publicly available metagenomic sequencing samples of SCTLD lesions and healthy tissues from 4 stony coral species. To compensate for the lack of coral reference genomes, we used data from apparently healthy coral samples to approximate a host genome and healthy microbiome reference. These reads were then used as a reference to which we matched and removed reads from diseased lesion tissue samples, and the remaining reads associated only with disease lesions were taxonomically classified at the DNA and protein levels. For DNA classifications, we used a pathogen identification protocol originally designed to identify pathogens in human tissue samples, and for protein classifications, we used a fast protein sequence aligner. To assess the utility of our pipeline, a species-level analysis of a candidate genus, Vibrio, was used to demonstrate the pipeline's effectiveness. Our approach revealed both complementary and unique coral microbiome members compared with a prior metagenome analysis of the same dataset.

LeishGEM Team. LeishGEM: genome-wide deletion mutant fitness and protein localisations in Leishmania

LeishGEM is a genome-wide functional annotation community resource for Leishmania mexicana, where deletion mutant growth in vitro and in vivo is measured and protein localisation is determined by endogenous tagging and LOPIT-DC (localisation of organelle proteins by isotope tagging with differential centrifugation) spatial proteomics. Data are being made available pre-publication via http://leishgem.org which allows data-driven identification of the mechanisms for Leishmania parasitism.

Pathogen-associated molecular patterns (PAMPs) derived from Leishmania and bacteria increase gene expression of antimicrobial peptides and gut surface proteins in sand flies

The interaction between pathogens and vectors' physiology can impact parasite transmission. Studying this interaction at the molecular level can help in developing control strategies. We study leishmaniases, diseases caused by Leishmania parasites transmitted by sand fly vectors, posing a significant global public health concern. Lipophosphoglycan (LPG), the major surface glycoconjugate of Leishmania, has been described to have several roles throughout the parasite's life cycle, both in the insect and vertebrate hosts. In addition, the sand fly midgut possesses a rich microbiota expressing lipopolysaccharides (LPS). However, the effect of LPG and LPS on the gene expression of sand fly midgut proteins or immunity effectors has not yet been documented. We experimentally fed Lutzomyia longipalpis and Phlebotomus papatasi sand flies with blood containing purified LPG from Leishmania infantum, Leishmania major, or LPS from Escherichia coli. The effect on the expression of genes encoding gut proteins galectin and mucin, digestive enzymes trypsin and chymotrypsin, and antimicrobial peptides (AMPs) attacin and defensins was assessed by quantitative PCR (qPCR). The gene expression of a mucin-like protein in L. longipalpis was increased by L. infantum LPG and E. coli LPS. The gene expression of a galectin was increased in L. longipalpis by L. major LPG, and in P. papatasi by E. coli LPS. Nevertheless, the gene expression of trypsins and chymotrypsins did not significantly change. On the other hand, both L. infantum and L. major LPG significantly enhanced expression of the AMP attacin in both sand fly species and defensin in L. longipalpis. In addition, E. coli LPS increased the expression of attacin and defensin in L. longipalpis. Our study showed that Leishmania LPG and E. coli LPS differentially modulate the expression of sand fly genes involved in gut maintenance and defence. This suggests that the glycoconjugates from microbiota or Leishmania may increase the vector's immune response and the gene expression of a gut coating protein in a permissive vector.

The influence of amoeba metal homeostasis on antifungal activity against Cryptococcus gattii

Free-living amoebas are natural predators of fungi, including human pathogens of the Cryptococcus genus. To survive and proliferate inside phagocytes, cryptococcal cells must acquire several nutrients. Zinc is fundamental for all life forms and develops a crucial role in the virulence of fungal pathogens, phagocytes reduce the availability of this metal to reduce the development of infection. The Acanthamoeba castellanii ACA1_271600 gene codes a metal transporter that is possibly associated with such antifungal strategy. Here, we evaluated the impact of A. castellanii metal homeostasis on C. gattii survival. Gene silencing of ACA1_271600 was performed and the interaction outcome of amoeba cells with both WT and zinc homeostasis-impaired mutant cryptococcal cells was evaluated. Decreased levels of ACA1_271600 in silenced amoeba cells led to higher proliferation of such cryptococcal strains. This effect was more pronounced in the zip1 mutant of C. gattii, suggesting that ACA1_271600 gene product modulates metal availability in Cryptococcus-infected amoebae. In addition, a systems biology analysis allowed us to infer that ACA1_271600 may also be involved in other biological processes that could compromise amoebae activity over cryptococcal cells. These results support the hypothesis that A. castellanii can apply nutritional immunity to hamper cryptococcal survival.

zDB: bacterial comparative genomics made easy

The analysis and comparison of genomes rely on different tools for tasks such as annotation, orthology prediction, and phylogenetic inference. Most tools are specialized for a single task, and additional efforts are necessary to integrate and visualize the results. To fill this gap, we developed zDB, an application integrating a Nextflow analysis pipeline and a Python visualization platform built on the Django framework. The application is available on GitHub (https://github.com/metagenlab/zDB) and from the bioconda channel. Starting from annotated Genbank files, zDB identifies orthologs and infers a phylogeny for each orthogroup. A species phylogeny is also constructed from shared single-copy orthologs. The results can be enriched with Pfam protein domain prediction, Cluster of Orthologs Genes and Kyoto Encyclopedia of Genes and Genomes annotations, and Swissprot homologs. The web application allows searching for specific genes or annotations, running Blast queries, and comparing genomic regions and whole genomes. The metabolic capacities of organisms can be compared at either the module or pathway levels. Finally, users can run queries to examine the conservation of specific genes or annotations across a chosen subset of genomes and display the results as a list of genes, Venn diagram, or heatmaps. Those features make zDB useful for both bioinformaticians and researchers more accustomed to laboratory research.IMPORTANCEGenome comparison and analysis rely on many independent tools, leaving to scientists the burden to integrate and visualize their results for interpretation. To alleviate this burden, we have built zDB, a comparative genomics tool that includes both an analysis pipeline and a visualization platform. The analysis pipeline automates gene annotation, orthology prediction, and phylogenetic inference, while the visualization platform allows scientists to easily explore the results in a web browser. Among other features, the interface allows users to visually compare whole genomes and targeted regions, assess the conservation of genes or metabolic pathways, perform Blast searches, or look for specific annotations. Altogether, this tool will be useful for a broad range of applications in comparative studies between two and hundred genomes. Furthermore, it is designed to allow sharing of data sets easily at a local or international scale, thereby supporting exploratory analyses for non-bioinformaticians on the genome of their favorite organisms.

Establishing a Male-Positive Genetic Sexing Strain in the Asian Malaria Vector Anopheles stephensi

Genetic biocontrol interventions targeting mosquito-borne diseases require the release of male mosquitoes exclusively, as only females consume blood and transmit human pathogens. This reduces the risk of spreading pathogens while enabling effective population control. Robust sex sorting methods to enable early larval sorting in mosquitoes need to be developed to allow for scalable sex sorting for genetic biocontrol interventions. This study applies the SEPARATOR (Sexing Element Produced by Alternative RNA-splicing of A Transgenic Observable Reporter) system, previously developed for Aedes aegypti, to the Asian malaria vector Anopheles stephensi. We hypothesized that the intron from the doublesex gene in Anopheles gambiae would function in An. stephensi due to evolutionary conservation. Our results confirm that the splicing module from An. gambiae operates effectively in An. stephensi, demonstrating evolutionary conservation in sex-specific splicing events between these species. This system enables reliable positive male selection from first instar larval to pupal stages. RT-PCR analysis demonstrates that male-specific EGFP expression is dependent on doublesex sex-specific splicing events. The SEPARATOR system's independence from sex-chromosome linkage confers resistance to meiotic recombination and chromosomal rearrangements. This approach may facilitate the mass release of males, and the cross-species portability of SEPARATOR establishes it as a valuable tool for genetic biocontrol interventions across various pest species.

Entamoeba histolytica: EhADH, an Alix Protein, Participates in Several Virulence Events through Its Different Domains

Entamoeba histolytica is the protozoan causative of human amoebiasis. The EhADH adhesin (687 aa) is a protein involved in tissue invasion, phagocytosis and host-cell lysis. EhADH adheres to the prey and follows its arrival to the multivesicular bodies. It is an accessory protein of the endosomal sorting complexes required for transport (ESCRT) machinery. Here, to study the role of different parts of EhADH during virulence events, we produced trophozoites overexpressing the three domains of EhADH, Bro1 (1-400 aa), Linker (246-446 aa) and Adh (444-687 aa) to evaluate their role in virulence. The TrophozBro11-400 slightly increased adherence and phagocytosis, but these trophozoites showed a higher ability to destroy cell monolayers, augment the permeability of cultured epithelial cells and mouse colon, and produce more damage to hamster livers. The TrophozLinker226-446 also increased the virulence properties, but with lower effect than the TrophozBro11-400. In addition, this fragment participates in cholesterol transport and GTPase binding. Interestingly, the TrophozAdh444-687 produced the highest effect on adherence and phagocytosis, but it poorly influenced the monolayers destruction; nevertheless, they augmented the colon and liver damage. To identify the protein partners of each domain, we used recombinant peptides. Pull-down assays and mass spectrometry showed that Bro1 domain interplays with EhADH, Gal/GalNAc lectin, EhCPs, ESCRT machinery components and cytoskeleton proteins. While EhADH, ubiquitin, EhRabB, EhNPC1 and EhHSP70 were associated to the Linker domain, and EhADH, EhHSP70, EhPrx and metabolic enzymes interacted to the Adh domain. The diverse protein association confirms that EhADH is a versatile molecule with multiple functions probably given by its capacity to form distinct molecular complexes.

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.

Imaging Assays to Detect DNA Damage in Trypanosome Parasites Using γH2A

Diseases caused by trypanosomatid parasites remain a significant unmet medical need for millions of people globally. Trypanosomatid parasites such as Trypanosoma cruzi and subspecies of Trypanosoma brucei cause Chagas disease and human African trypanosomiasis (HAT), respectively. Although efforts to find novel treatments have been successful for HAT, Chagas disease is still treated with decades-old therapies that suffer from long treatment durations and severe safety concerns. We recently described the identification and characterization of the cyanotriazole compound class that kills trypanosomes, in vitro and in vivo, by selective inhibition of the trypanosome nuclear topoisomerase II enzyme. To evaluate whether inhibition of the topoisomerase II enzyme led to parasite death due to lethal double-strand DNA breaks, we developed assays for detecting DNA damage in both intracellular amastigotes of T. cruzi and bloodstream-form T. brucei by using the canonical DNA damage marker γH2A. Herein, this article describes the protocols for detecting DNA damage using an immunofluorescence assessment of γH2A by microscopy in trypanosome parasites. Key features • Immunofluorescence-based assay to detect the γH2A response in T. brucei and T. cruzi parasites. • Robust DNA damage pathway-based cellular assays to evaluate topoisomerase II poisons' ability to cause DNA damage. • A 384-well plate-based T. cruzi protocol allows high-resolution and high-throughput evaluation of compounds that cause DNA damage by measuring γH2A in intracellular parasites. • This assay could be modifiable for evaluation of DNA damage responses in various intracellular and extracellular eukaryotic pathogens.

Endosymbiosis in trypanosomatids: The bacterium division depends on microtubule dynamism

Microtubules are components of the cytoskeleton that perform essential functions in eukaryotes, such as those related to shape change, motility and cell division. In this context some characteristics of these filaments are essential, such as polarity and dynamic instability. In trypanosomatids, microtubules are integral to ultrastructure organization, intracellular transport and mitotic processes. Some species of trypanosomatids co-evolve with a symbiotic bacterium in a mutualistic association that is marked by extensive metabolic exchanges and a coordinated division of the symbiont with other cellular structures, such as the nucleus and the kinetoplast. It is already established that the bacterium division is microtubule-dependent, so in this work, it was investigated whether the dynamism and remodeling of these filaments is capable of affecting the prokaryote division. To this purpose, Angomonas deanei was treated with Trichostatin A (TSA), a deacetylase inhibitor, and mutant cells for histone deacetylase 6 (HDAC6) were obtained by CRISPR-Cas9. A decrease in proliferation, an enhancement in tubulin acetylation, as well as morphological and ultrastructural changes, were observed in TSA-treated protozoa and mutant cells. In both cases, symbiont filamentation occurred, indicating that prokaryote cell division is dependent on microtubule dynamism.

BCC0 collaborates with IMC32 and IMC43 to form the Toxoplasma gondii essential daughter bud assembly complex

Toxoplasma gondii divides by endodyogeny, in which two daughter buds are formed within the cytoplasm of the maternal cell using the inner membrane complex (IMC) as a scaffold. During endodyogeny, components of the IMC are synthesized and added sequentially to the nascent daughter buds in a tightly regulated manner. We previously showed that the early recruiting proteins IMC32 and IMC43 form an essential daughter bud assembly complex which lays the foundation of the daughter cell scaffold in T. gondii. In this study, we identify the essential, early recruiting IMC protein BCC0 as a third member of this complex by using IMC32 as bait in both proximity labeling and yeast two-hybrid screens. We demonstrate that BCC0's localization to daughter buds depends on the presence of both IMC32 and IMC43. Deletion analyses and functional complementation studies reveal that residues 701-877 of BCC0 are essential for both its localization and function and that residues 1-899 are sufficient for function despite minor mislocalization. Pairwise yeast two-hybrid assays additionally demonstrate that BCC0's essential domain binds to the coiled-coil region of IMC32 and that BCC0 and IMC43 do not directly interact. This data supports a model for complex assembly in which an IMC32-BCC0 subcomplex initially recruits to nascent buds via palmitoylation of IMC32 and is locked into the scaffold once bud elongation begins by IMC32 binding to IMC43. Together, this study dissects the organization and function of a complex of three early recruiting daughter proteins which are essential for the proper assembly of the IMC during endodyogeny.

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.

Plasmodium infection induces phenotypic, clonal, and spatial diversity among differentiating CD4+ T cells

Naive CD4+ T cells must differentiate in order to orchestrate immunity to Plasmodium, yet understanding of their emerging phenotypes, clonality, spatial distributions, and cellular interactions remains incomplete. Here, we observe that splenic polyclonal CD4+ T cells differentiate toward T helper 1 (Th1) and T follicular helper (Tfh)-like states and exhibit rarer phenotypes not elicited among T cell receptor (TCR) transgenic counterparts. TCR clones present at higher frequencies exhibit Th1 skewing, suggesting that variation in major histocompatibility complex class II (MHC-II) interaction influences proliferation and Th1 differentiation. To characterize CD4+ T cell interactions, we map splenic microarchitecture, cellular locations, and molecular interactions using spatial transcriptomics at near single-cell resolution. Tfh-like cells co-locate with stromal cells in B cell follicles, while Th1 cells in red pulp co-locate with activated monocytes expressing multiple chemokines and MHC-II. Spatial mapping of individual transcriptomes suggests that proximity to chemokine-expressing monocytes correlates with stronger effector phenotypes in Th1 cells. Finally, CRISPR-Cas9 gene disruption reveals a role for CCR5 in promoting clonal expansion and Th1 differentiation. A database of cellular locations and interactions is presented: https://haquelab.mdhs.unimelb.edu.au/spatial_gui/.

Novel metagenomics analysis of stony coral tissue loss disease

Stony coral tissue loss disease (SCTLD) has devastated coral reefs off the coast of Florida and continues to spread throughout the Caribbean. Although a number of bacterial taxa have consistently been associated with SCTLD, no pathogen has been definitively implicated in the etiology of SCTLD. Previous studies have predominantly focused on the prokaryotic community through 16S rRNA sequencing of healthy and affected tissues. Here, we provide a different analytical approach by applying a bioinformatics pipeline to publicly available metagenomic sequencing samples of SCTLD lesions and healthy tissues from four stony coral species. To compensate for the lack of coral reference genomes, we used data from apparently healthy coral samples to approximate a host genome and healthy microbiome reference. These reads were then used as a reference to which we matched and removed reads from diseased lesion tissue samples, and the remaining reads associated only with disease lesions were taxonomically classified at the DNA and protein levels. For DNA classifications, we used a pathogen identification protocol originally designed to identify pathogens in human tissue samples, and for protein classifications, we used a fast protein sequence aligner. To assess the utility of our pipeline, a species-level analysis of a candidate genus, Vibrio, was used to demonstrate the pipeline's effectiveness. Our approach revealed both complementary and unique coral microbiome members compared to a prior metagenome analysis of the same dataset.

Plasmodium LCCL domain-containing modular proteins have their origins in the ancestral alveolate

Plasmodium species encode a unique set of six modular proteins named LCCL lectin domain adhesive-like proteins (LAPs) that operate as a complex and that are essential for malaria parasite transmission from mosquito to vertebrate. LAPs possess complex architectures obtained through unique assemblies of conserved domains associated with lipid, protein and carbohydrate interactions, including the name-defining LCCL domain. Here, we assessed the prevalence of Plasmodium LAP orthologues across eukaryotic life. Our findings show orthologous conservation in all apicomplexans, with lineage-specific repertoires acquired through differential lap gene loss and duplication. Besides Apicomplexa, LAPs are found in their closest relatives: the photosynthetic chromerids, which encode the broadest repertoire including a novel membrane-bound LCCL protein. LAPs are notably absent from other alveolate lineages (dinoflagellates, perkinsids and ciliates), but are encoded by predatory colponemids, a sister group to the alveolates. These results reveal that the LAPs are much older than previously thought and pre-date not only the Apicomplexa but the Alveolata altogether.

Properties and predicted functions of large genes and proteins of apicomplexan parasites

Evolutionary constraints greatly favor compact genomes that efficiently encode proteins. However, several eukaryotic organisms, including apicomplexan parasites such as Toxoplasma gondii, Plasmodium falciparum and Babesia duncani, the causative agents of toxoplasmosis, malaria and babesiosis, respectively, encode very large proteins, exceeding 20 times their average protein size. Although these large proteins represent <1% of the total protein pool and are generally expressed at low levels, their persistence throughout evolution raises important questions about their functions and possible evolutionary pressures to maintain them. In this study, we examined the trends in gene and protein size, function and expression patterns within seven apicomplexan pathogens. Our analysis revealed that certain large proteins in apicomplexan parasites harbor domains potentially important for functions such as antigenic variation, erythrocyte invasion and immune evasion. However, these domains are not limited to or strictly conserved within large proteins. While some of these proteins are predicted to engage in conventional metabolic pathways within these parasites, others fulfill specialized functions for pathogen-host interactions, nutrient acquisition and overall survival.

The many roles of sulfur in the fungal-host interaction

Sulfur is an essential macronutrient for life, and consequently, all living organisms must acquire it from external sources to thrive and grow. Sulfur is a constituent of a multitude of crucial molecules, such as the S-containing proteinogenic amino acids cysteine and methionine; cofactors and prosthetic groups, such as coenzyme-A and iron-sulfur (Fe-S) clusters; and other essential organic molecules, such as glutathione or S-adenosylmethionine. Additionally, sulfur in cysteine thiols is an active redox group that plays paramount roles in protein stability, enzyme catalysis, and redox homeostasis. Furthermore, H2S is gaining more attention as a crucial signaling molecule that influences metabolism and physiological functions. Given its importance, it is not surprising that sulfur plays key roles in the host-pathogen interaction. However, in contrast to its well-recognized involvement in the plant-pathogen interaction, the specific contributions of sulfur to the human-fungal interaction are much less understood. In this short review, I highlight some of the most important known mechanisms and propose directions for further research.

Elucidating the spatio-temporal dynamics of the Plasmodium falciparum basal complex

Asexual replication of Plasmodium falciparum occurs via schizogony, wherein 16-36 daughter cells are produced within the parasite during one semi-synchronized cytokinetic event. Schizogony requires a divergent contractile ring structure known as the basal complex. Our lab has previously identified PfMyoJ (PF3D7_1229800) and PfSLACR (PF3D7_0214700) as basal complex proteins recruited midway through segmentation. Using ultrastructure expansion microscopy, we localized both proteins to a novel basal complex subcompartment. While both colocalize with the basal complex protein PfCINCH upon recruitment, they form a separate, more basal subcompartment termed the posterior cup during contraction. We also show that PfSLACR is recruited to the basal complex prior to PfMyoJ, and that both proteins are removed unevenly as segmentation concludes. Using live-cell microscopy, we show that actin dynamics are dispensable for basal complex formation, expansion, and contraction. We then show that EF-hand containing P. falciparum Centrin 2 partially localizes to this posterior cup of the basal complex and that it is essential for growth and replication, with variable defects in basal complex contraction and synchrony. Finally, we demonstrate that free intracellular calcium is necessary but not sufficient for basal complex contraction in P. falciparum. Thus, we demonstrate dynamic spatial compartmentalization of the Plasmodium falciparum basal complex, identify an additional basal complex protein, and begin to elucidate the unique mechanism of contraction utilized by P. falciparum, opening the door for further exploration of Apicomplexan cellular division.

Lack of paternal silencing and ecotype-specific expression in head and body lice hybrids

Paternal genome elimination (PGE) is a non-Mendelian inheritance system, described in numerous arthropod species, in which males develop from fertilized eggs, but their paternally inherited chromosomes are eliminated before or during spermatogenesis. Therefore, PGE males only transmit their maternally inherited set of chromosomes to their offspring. In addition to the elimination of paternal chromosomes, diverse PGE species have also repeatedly evolved the transcriptional silencing of the paternal genome, making males effectively haploid. However, it is unclear if this paternal chromosome silencing is mechanistically linked to the chromosome elimination or has evolved at a later stage, and if so, what drives the haploidization of males under PGE. In order to understand these questions, here we study the human louse, Pediculus humanus, which represents an ideal model system, as it appears to be the only instance of PGE where males eliminate, but not silence their paternal chromosomes, although the latter remains to be shown conclusively. In this study, we analyzed parent-of-origin allele-specific expression patterns in male offspring of crosses between head and body lice ecotypes. We show that hybrid adult males of P. humanus display biparental gene expression, which constitutes the first case of a species with PGE in which genetic activity of paternal chromosomes in the soma is not affected by embryonic silencing or (partial or complete) elimination. We did however also identify a small number of maternally biased genes (potentially imprinted genes), which may be involved in the elimination of paternal chromosomes during spermatogenesis. Finally, we have identified genes that show ecotype-specific expression bias. Given the low genetic diversity between ecotypes, this is suggestive for a role of epigenetic processes in ecotype differences.

scRNA-seq reveals transcriptional dynamics of Encephalitozoon intestinalis parasites in human macrophages

Microsporidia are single-celled intracellular parasites that cause opportunistic diseases in humans. Encephalitozoon intestinalis is a prevalent human-infecting species that invades the small intestine. Dissemination to other organ systems is also observed, and is potentially facilitated by macrophages. The macrophage response to infection and the developmental trajectory of the parasite are not well studied. Here we use single cell RNA sequencing to investigate transcriptional changes in both the host and parasite during infection. While a small population of infected macrophages mount a response, most remain transcriptionally unchanged, suggesting that the majority of parasites may avoid host detection. The parasite transcriptome reveals large transcriptional changes throughout the life cycle, providing a blueprint for parasite development. The stealthy microsporidian lifestyle likely allows these parasites to harness macrophages for replication and dissemination. Together, our data provide insights into the host response in primary human macrophages and the E. intestinalis developmental program.

Transgene removal using an in cis programmed homing endonuclease via single-strand annealing in the mosquito Aedes aegypti

While gene drive strategies have been proposed to aid in the control of mosquito-borne diseases, additional genome engineering technologies may be required to establish a defined end-of-product-life timeline. We previously demonstrated that single-strand annealing (SSA) was sufficient to program the scarless elimination of a transgene while restoring a disrupted gene in the disease vector mosquito Aedes aegypti. Here, we extend these findings by establishing that complete transgene removal (four gene cassettes comprising ~8-kb) can be programmed in cis. Reducing the length of the direct repeat from 700-bp to 200-bp reduces, but does not eliminate, SSA activity. In contrast, increasing direct repeat length to 1.5-kb does not increase SSA rates, suggesting diminishing returns above a certain threshold size. Finally, we show that while the homing endonuclease Y2-I-AniI triggered both SSA and NHEJ at significantly higher rates than I-SceI at one genomic locus (P5-EGFP), repair events are heavily skewed towards NHEJ at another locus (kmo), suggesting the nuclease used and the genomic region targeted have a substantial influence on repair outcomes. Taken together, this work establishes the feasibility of engineering temporary transgenes in disease vector mosquitoes, while providing critical details concerning important operational parameters.

Characterization and functional analysis of Toxoplasma Golgi-associated proteins identified by proximity labelling

Toxoplasma gondii possesses a highly polarized secretory pathway that contains both broadly conserved eukaryotic organelles and unique apicomplexan organelles which play essential roles in the parasite's lytic cycle. As in other eukaryotes, the T. gondii Golgi apparatus sorts and modifies proteins prior to their distribution to downstream organelles. Many of the typical trafficking factors found involved in these processes are missing from apicomplexan genomes, suggesting that these parasites have evolved unique proteins to fill these roles. Here we identify a novel Golgi-localizing protein (ULP1) which contains structural homology to the eukaryotic trafficking factor p115/Uso1. We demonstrate that depletion of ULP1 leads to a dramatic reduction in parasite fitness and replicative ability. Using ULP1 as bait for TurboID proximity labelling and immunoprecipitation, we identify eleven more novel Golgi-associated proteins and demonstrate that ULP1 interacts with the T. gondii COG complex. These proteins include both conserved trafficking factors and parasite-specific proteins. Using a conditional knockdown approach, we assess the effect of each of these eleven proteins on parasite fitness. Together, this work reveals a diverse set of novel T. gondii Golgi-associated proteins that play distinct roles in the secretory pathway. As several of these proteins are absent outside of the Apicomplexa, they represent potential targets for the development of novel therapeutics against these parasites.

Importance

Apicomplexan parasites such as Toxoplasma gondii infect a large percentage of the world's population and cause substantial human disease. These widespread pathogens use specialized secretory organelles to infect their host cells, modulate host cell functions, and cause disease. While the functions of the secretory organelles are now better understood, the Golgi apparatus of the parasite remains largely unexplored, particularly regarding parasite-specific innovations that may help direct traffic intracellularly. In this work, we characterize ULP1, a protein that is unique to parasites but shares structural similarity to the eukaryotic trafficking factor p115/Uso1. We show that ULP1 plays an important role in parasite replication and demonstrate that it interacts with the conserved oligomeric Golgi (COG) complex. We then use ULP1 proximity labelling to identify eleven additional Golgi-associated proteins which we functionally analyze via conditional knockdown. This work expands our knowledge of the Toxoplasma Golgi apparatus and identifies potential targets for therapeutic intervention.

Insights into copper sensing and tolerance in Pneumocystis species

Introduction

Pneumocystis species are pathogenic fungi known to cause pneumonia in immunocompromised mammals. They are obligate to their host, replicate extracellularly in lung alveoli and thrive in the copper-enriched environment of mammalian lungs. In this study, we investigated the proteome of Pneumocystis murina, a model organism that infects mice, in the context of its copper sensing and tolerance.

Methods and results

The query for copper-associated annotations in FungiDB followed by a manual curation identified only 21 genes in P. murina, significantly fewer compared to other clinically relevant fungal pathogens or phylogenetically similar free-living fungi. We then employed instrumental analyses, including Size-Exclusion Chromatography Inductively Coupled Plasma Mass Spectrometry (SEC-ICP-MS), Immobilized Metal Affinity Chromatography (IMAC), and Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS), to isolate and identify copper-binding proteins from freshly extracted organisms, revealing 29 distinct cuproproteins. The RNA sequencing (RNA-seq) analysis of P. murina exposed to various CuSO4 concentrations at three temporal intervals (0.5, 2, and 5 h) indicated that significant gene expression changes occurred only under the highest CuSO4 concentration probed (100 μM) and the longest exposure duration (5 h). This stimulus led to the upregulation of 43 genes and downregulation of 27 genes compared to untreated controls. Quantitative PCR (qPCR) confirmed the expression of four out of eight selected upregulated genes, including three assumed transcription factors (PNEG_01236, PNEG_01675, and PNEG_01730) and a putative copper transporter (PNEG_02609). Notably, the three applied methodologies - homology-based annotation, SEC-ICP-MS/IMAC/LC-MS/MS, and RNA-seq - yielded largely distinct findings, with only four genes (PNEG_02587, PNEG_03319, PNEG_02584, and PNEG_02989) identified by both instrumental methods.

Discussion

The insights contribute to the broader knowledge of Pneumocystis copper homeostasis and provide novel facets of host-pathogen interactions for extracellular pathogens. We suggest that future studies of Pneumocystis pathogenicity and copper stress survival should consider the entire spectrum of identified genes.

Identification of 30 transition fibre proteins in Trypanosoma brucei reveals a complex and dynamic structure

Transition fibres and distal appendages surround the distal end of mature basal bodies and are essential for ciliogenesis, but only a few of the proteins involved have been identified and functionally characterised. Here, through genome-wide analysis, we have identified 30 transition fibre proteins (TFPs) and mapped their arrangement in the flagellated eukaryote Trypanosoma brucei. We discovered that TFPs are recruited to the mature basal body before and after basal body duplication, with differential expression of five TFPs observed at the assembling new flagellum compared to the existing fixed-length old flagellum. RNAi-mediated depletion of 17 TFPs revealed six TFPs that are necessary for ciliogenesis and a further three TFPs that are necessary for normal flagellum length. We identified nine TFPs that had a detectable orthologue in at least one basal body-forming eukaryotic organism outside of the kinetoplastid parasites. Our work has tripled the number of known transition fibre components, demonstrating that transition fibres are complex and dynamic in their composition throughout the cell cycle, which relates to their essential roles in ciliogenesis and flagellum length regulation.

Identification of domains in Plasmodium falciparum proteins of unknown function using DALI search on AlphaFold predictions

Plasmodium falciparum, the causative agent of malaria, poses a significant global health challenge, yet much of its biology remains elusive. A third of the genes in the P. falciparum genome lack annotations regarding their function, impeding our understanding of the parasite's biology. In this study, we employ structure predictions and the DALI search algorithm to analyse proteins encoded by uncharacterized genes in the reference strain 3D7 of P. falciparum. By comparing AlphaFold predictions to experimentally determined protein structures in the Protein Data Bank, we found similarities to known domains in 353 proteins of unknown function, shedding light on their potential functions. The lowest-scoring 5% of similarities were additionally validated using the size-independent TM-align algorithm, confirming the detected similarities in 88% of the cases. Notably, in over 70 P. falciparum proteins the presence of domains resembling heptatricopeptide repeats, which are typically involvement in RNA binding and processing, was detected. This suggests this family, which is important in transcription in mitochondria and apicoplasts, is much larger in Plasmodium parasites than previously thought. The results of this domain search provide a resource to the malaria research community that is expected to inform and enable experimental studies.

Shared Features Underlying Compact Genomes and Extreme Habitat Use in Chironomid Midges

Nonbiting midges (family Chironomidae) are found throughout the world in a diverse array of aquatic and terrestrial habitats, can often tolerate harsh conditions such as hypoxia or desiccation, and have consistently compact genomes. Yet we know little about the shared molecular basis for these attributes and how they have evolved across the family. Here, we address these questions by first creating high-quality, annotated reference assemblies for Tanytarsus gracilentus (subfamily Chironominae, tribe Tanytarsini) and Parochlus steinenii (subfamily Podonominae). Using these and other publicly available assemblies, we created a time-calibrated phylogenomic tree for family Chironomidae with outgroups from order Diptera. We used this phylogeny to test for features associated with compact genomes, as well as examining patterns of gene family evolution and positive selection that may underlie chironomid habitat tolerances. Our results suggest that compact genomes evolved in the common ancestor of Chironomidae and Ceratopogonidae and that this occurred mainly through reductions in noncoding regions (introns, intergenic sequences, and repeat elements). Significantly expanded gene families in Chironomidae included biological processes that may relate to tolerance of stressful environments, such as temperature homeostasis, carbohydrate transport, melanization defense response, and trehalose transport. We identified several positively selected genes in Chironomidae, notably sulfonylurea receptor, CREB-binding protein, and protein kinase D. Our results improve our understanding of the evolution of small genomes and extreme habitat use in this widely distributed group.

Systematic characterization of all Toxoplasma gondii TBC domain-containing proteins identifies an essential regulator of Rab2 in the secretory pathway

Toxoplasma gondii resides in its intracellular niche by employing a series of specialized secretory organelles that play roles in invasion, host cell manipulation, and parasite replication. Rab GTPases are major regulators of the parasite's secretory traffic that function as nucleotide-dependent molecular switches to control vesicle trafficking. While many of the Rab proteins have been characterized in T. gondii, precisely how these Rabs are regulated remains poorly understood. To better understand the parasite's secretory traffic, we investigated the entire family of Tre2-Bub2-Cdc16 (TBC) domain-containing proteins, which are known to be involved in vesicle fusion and secretory protein trafficking. We first determined the localization of all 18 TBC domain-containing proteins to discrete regions of the secretory pathway or other vesicles in the parasite. Second, we use an auxin-inducible degron approach to demonstrate that the protozoan-specific TgTBC9 protein, which localizes to the endoplasmic reticulum (ER), is essential for parasite survival. Knockdown of TgTBC9 results in parasite growth arrest and affects the organization of the ER and mitochondrial morphology. TgTBC9 knockdown also results in the formation of large lipid droplets (LDs) and multi-membranous structures surrounded by ER membranes, further indicating a disruption of ER functions. We show that the conserved dual-finger active site in the TBC domain of the protein is critical for its GTPase-activating protein (GAP) function and that the Plasmodium falciparum orthologue of TgTBC9 can rescue the lethal knockdown. We additionally show by immunoprecipitation and yeast 2 hybrid analyses that TgTBC9 preferentially binds Rab2, indicating that the TBC9-Rab2 pair controls ER morphology and vesicular trafficking in the parasite. Together, these studies identify the first essential TBC protein described in any protozoan and provide new insight into intracellular vesicle trafficking in T. gondii.

Two Nested Inversions in the X Chromosome Differentiate the Dominant Malaria Vectors in Europe, Anopheles atroparvus and Anopheles messeae

The Maculipennis subgroup of malaria mosquitoes includes both dominant malaria vectors and non-vectors in Eurasia. Understanding the genetic factors, particularly chromosomal inversions, that differentiate Anopheles species can provide valuable insights for vector control strategies. Although autosomal inversions between the species in this subgroup have been characterized based on the chromosomal banding patterns, the number and positions of rearrangements in the X chromosome remain unclear due to the divergent banding patterns. Here, we identified two large X chromosomal inversions, approximately 13 Mb and 10 Mb in size, using fluorescence in situ hybridization. The inversion breakpoint regions were mapped by hybridizing 53 gene markers with polytene chromosomes of An. messeae. The DNA probes were designed based on gene sequences from the annotated An. atroparvus genome. The two nested inversions resulted in five syntenic blocks. Only two small syntenic blocks, which encompass 181 annotated genes in the An. atroparvus genome, changed their position and orientation in the X chromosome. The analysis of the An. atroparvus genome revealed an enrichment of gene ontology terms associated with immune system and mating behavior in the rearranged syntenic blocks. Additionally, the enrichment of DNA transposons was found in sequences homologous to three of the four breakpoint regions. This study demonstrates the successful application of the physical genome mapping approach to identify rearrangements that differentiate species in insects with polytene chromosomes.

Comparative transcriptional analysis identifies genes associated with the attenuation of Theileria parva infected cells after long-term in vitro culture

Autologous administration of attenuated Theileria parva-infected cells induces immunity to T. parva in cattle. The mechanism of attenuation, however, is largely unknown. Here, we used RNA sequencing of pathogenic and attenuated T. parva-infected T-cells to elucidate the transcriptional changes underpinning attenuation. We observed differential expression of several host genes, including TRAIL, PD-1, TGF-β and granzymes that are known to regulate inflammation and proliferation of infected cells. Importantly, many genes linked with the attenuation of the related T. annulata-infected cells were not dysregulated in this study. Furthermore, known T. parva antigens were not dysregulated in attenuated relative to pathogenic cells, indicating that attenuation is not due to enhanced immunogenicity. Overall this study suggests that attenuation is driven by a decrease in proliferation and restoration of the inflammatory profile of T. parva-infected cells. Additionally, it provides a foundation for future mechanistic studies of the attenuation phenotype in Theileria-infected cells.

Cryogenic electron tomography reveals novel structures in the apical complex of Plasmodium falciparum

Intracellular infectious agents, like the malaria parasite, Plasmodium falciparum, face the daunting challenge of how to invade a host cell. This problem may be even harder when the host cell in question is the enucleated red blood cell, which lacks the host machinery co-opted by many pathogens for internalization. Evolution has provided P. falciparum and related single-celled parasites within the phylum Apicomplexa with a collection of organelles at their apical end that mediate invasion. This apical complex includes at least two sets of secretory organelles, micronemes and rhoptries, and several structural features like apical rings and a putative pore through which proteins may be introduced into the host cell during invasion. We perform cryogenic electron tomography (cryo-ET) equipped with Volta Phase Plate on isolated and vitrified merozoites to visualize the apical machinery. Through tomographic reconstruction of cellular compartments, we see new details of known structures like the rhoptry tip interacting directly with a rosette resembling the recently described rhoptry secretory apparatus (RSA), or with an apical vesicle docked beneath the RSA. Subtomogram averaging reveals that the apical rings have a fixed number of repeating units, each of which is similar in overall size and shape to the units in the apical rings of tachyzoites of Toxoplasma gondii. Comparison of these polar rings in Plasmodium and Toxoplasma parasites also reveals them to have a structurally conserved assembly pattern. These results provide new insight into the essential and structurally conserved features of this remarkable machinery used by apicomplexan parasites to invade their respective host cells.

IMPORTANCE

Malaria is an infectious disease caused by parasites of the genus Plasmodium and is a leading cause of morbidity and mortality globally. Upon infection, Plasmodium parasites invade and replicate in red blood cells, where they are largely protected from the immune system. To enter host cells, the parasites employ a specialized apparatus at their anterior end. In this study, advanced imaging techniques like cryogenic electron tomography (cryo-ET) and Volta Phase Plate enable unprecedented visualization of whole Plasmodium falciparum merozoites, revealing previously unknown structural details of their invasion machinery. Key findings include new insights into the structural conservation of apical rings shared between Plasmodium and its apicomplexan cousin, Toxoplasma. These discoveries shed light on the essential and conserved elements of the invasion machinery used by these pathogens. Moreover, the research provides a foundation for understanding the molecular mechanisms underlying parasite-host interactions, potentially informing strategies for combating diseases caused by apicomplexan parasites.

Integrating high-throughput analysis to create an atlas of replication origins in Trypanosoma cruzi in the context of genome structure and variability

Trypanosoma cruzi is the etiologic agent of the most prevalent human parasitic disease in Latin America, Chagas disease. Its genome is rich in multigenic families that code for virulent antigens and are present in the rapidly evolving genomic compartment named Disruptive. DNA replication is a meticulous biological process in which flaws can generate mutations and changes in chromosomal and gene copy numbers. Here, integrating high-throughput and single-molecule analyses, we were able to identify Predominant, Flexible, and Dormant Orc1Cdc6-dependent origins as well as Orc1Cdc6-independent origins. Orc1Cdc6-dependent origins were found in multigenic family loci, while independent origins were found in the Core compartment that contains conserved and hypothetical protein-coding genes, in addition to multigenic families. In addition, we found that Orc1Cdc6 density is related to the firing of origins and that Orc1Cdc6-binding sites within fired origins are depleted of a specific class of nucleosomes that we previously categorized as dynamic. Together, these data suggest that Orc1Cdc6-dependent origins may contribute to the rapid evolution of the Disruptive compartment and, therefore, to the success of T. cruzi infection and that the local epigenome landscape is also involved in this process.IMPORTANCETrypanosoma cruzi, responsible for Chagas disease, affects millions globally, particularly in Latin America. Lack of vaccine or treatment underscores the need for research. Parasite's genome, with virulent antigen-coding multigenic families, resides in the rapidly evolving Disruptive compartment. Study sheds light on the parasite's dynamic DNA replication, discussing the evolution of the Disruptive compartment. Therefore, the findings represent a significant stride in comprehending T. cruzi's biology and the molecular bases that contribute to the success of infection caused by this parasite.

Characterization of the First Secreted Sorting Nexin Identified in the Leishmania Protists

Proteins of the sorting nexin (SNX) family present a modular structural architecture with a phox homology (PX) phosphoinositide (PI)-binding domain and additional PX structural domains, conferring to them a wide variety of vital eukaryotic cell's functions, from signal transduction to membrane deformation and cargo binding. Although SNXs are well studied in human and yeasts, they are poorly investigated in protists. Herein, is presented the characterization of the first SNX identified in Leishmania protozoan parasites encoded by the LdBPK_352470 gene. In silico secondary and tertiary structure prediction revealed a PX domain on the N-terminal half and a Bin/amphiphysin/Rvs (BAR) domain on the C-terminal half of this protein, with these features classifying it in the SNX-BAR subfamily of SNXs. We named the LdBPK_352470.1 gene product LdSNXi, as it is the first SNX identified in Leishmania (L.) donovani. Its expression was confirmed in L. donovani promastigotes under different cell cycle phases, and it was shown to be secreted in the extracellular medium. Using an in vitro lipid binding assay, it was demonstrated that recombinant (r) LdSNXi (rGST-LdSNXi) tagged with glutathione-S-transferase (GST) binds to the PtdIns3P and PtdIns4P PIs. Using a specific a-LdSNXi antibody and immunofluorescence confocal microscopy, the intracellular localization of endogenous LdSNXi was analyzed in L. donovani promastigotes and axenic amastigotes. Additionally, rLdSNXi tagged with enhanced green fluorescent protein (rLdSNXi-EGFP) was heterologously expressed in transfected HeLa cells and its localization was examined. All observed localizations suggest functions compatible with the postulated SNX identity of LdSNXi. Sequence, structure, and evolutionary analysis revealed high homology between LdSNXi and the human SNX2, while the investigation of protein-protein interactions based on STRING (v.11.5) predicted putative molecular partners of LdSNXi in Leishmania.

Myosin F controls actin organization and dynamics in Toxoplasma gondii

Intracellular cargo transport is a ubiquitous cellular process in all eukaryotes. In many cell types, membrane bound cargo is associated with molecular motors which transport cargo along microtubule and actin tracks. In Toxoplasma gondii (T. gondii), an obligate intracellular parasite in the phylum Apicomplexa, organization of the endomembrane pathway depends on actin and an unconventional myosin motor, myosin F (MyoF). Loss of MyoF and actin disrupts vesicle transport, organelle positioning, and division of the apicoplast, a nonphotosynthetic plastid organelle. How this actomyosin system contributes to these cellular functions is still unclear. Using live-cell imaging, we observed that MyoF-EmeraldFP (MyoF-EmFP) displayed a dynamic and filamentous-like organization in the parasite cytosol, reminiscent of cytosolic actin filament dynamics. MyoF was not associated with the Golgi, apicoplast or dense granule surfaces, suggesting that it does not function using the canonical cargo transport mechanism. Instead, we found that loss of MyoF resulted in a dramatic rearrangement of the actin cytoskeleton in interphase parasites accompanied by significantly reduced actin dynamics. However, actin organization during parasite replication and motility was unaffected by the loss of MyoF. These findings revealed that MyoF is an actin organizing protein in Toxoplasma and facilitates cargo movement using an unconventional transport mechanism.

Endoplasmic reticulum localized TMEM33 domain-containing protein is crucial for all life cycle stages of the malaria parasite

Endoplasmic reticulum (ER) plays a pivotal role in the regulation of stress responses in multiple eukaryotic cells. However, little is known about the effector mechanisms that regulate stress responses in ER of the malaria parasite. Herein, we aimed to identify the importance of a transmembrane protein 33 (TMEM33)-domain-containing protein in life cycle of the rodent malaria parasite Plasmodium berghei. TMEM33 is an ER membrane-resident protein that is involved in regulating stress responses in various eukaryotic cells. A C-terminal tagged TMEM33 was localized in the ER throughout the blood and mosquito stages of development. Targeted deletion of TMEM33 confirmed its importance for asexual blood stages and ookinete development, in addition to its essential role for sporozoite infectivity in the mammalian host. Pilot scale analysis shows that the loss of TMEM33 results in the initiation of ER stress response and induction of autophagy. Our findings conclude an important role of TMEM33 in the development of all life cycle stages of the malaria parasite, which indicates its potential as an antimalarial target.

FAIR Header Reference genome: a TRUSTworthy standard

The lack of interoperable data standards among reference genome data-sharing platforms inhibits cross-platform analysis while increasing the risk of data provenance loss. Here, we describe the FAIR bioHeaders Reference genome (FHR), a metadata standard guided by the principles of Findability, Accessibility, Interoperability and Reuse (FAIR) in addition to the principles of Transparency, Responsibility, User focus, Sustainability and Technology. The objective of FHR is to provide an extensive set of data serialisation methods and minimum data field requirements while still maintaining extensibility, flexibility and expressivity in an increasingly decentralised genomic data ecosystem. The effort needed to implement FHR is low; FHR's design philosophy ensures easy implementation while retaining the benefits gained from recording both machine and human-readable provenance.

Ten common issues with reference sequence databases and how to mitigate them

Metagenomic sequencing has revolutionized our understanding of microbiology. While metagenomic tools and approaches have been extensively evaluated and benchmarked, far less attention has been given to the reference sequence database used in metagenomic classification. Issues with reference sequence databases are pervasive. Database contamination is the most recognized issue in the literature; however, it remains relatively unmitigated in most analyses. Other common issues with reference sequence databases include taxonomic errors, inappropriate inclusion and exclusion criteria, and sequence content errors. This review covers ten common issues with reference sequence databases and the potential downstream consequences of these issues. Mitigation measures are discussed for each issue, including bioinformatic tools and database curation strategies. Together, these strategies present a path towards more accurate, reproducible and translatable metagenomic sequencing.

Genomic diversity of the pathogenic fungus Aspergillus fumigatus in Japan reveals the complex genomic basis of azole resistance

Aspergillus fumigatus is a pathogenic fungus with a global distribution. The emergence of azole-resistant A. fumigatus (ARAf) other than the TR-mutants is a problem in Japan. Additionally, the genetic diversity of A. fumigatus strains in Japan remains relatively unknown. Here we show the diversity in the A. fumigatus strains isolated in Japan as well as the complexity in the global distribution of the pathogenic strains. First, we analyzed the genome sequences of 171 strains from Japan as well as the antifungal susceptibility of these strains. Next, we conducted a population analysis of 876 strains by combining the available genomic data for strains isolated worldwide, which were grouped in six clusters. Finally, a genome-wide association study identified the genomic loci associated with ARAf strains, but not the TR-mutants. These results highlight the complexity of the genomic mechanism underlying the emergence of ARAf strains other than the TR-mutants.