Transcriptomic analysis of albendazole resistance in human diarrheal parasite Giardia duodenalis.

Benzimidazole-2-carbamates (BZ, e.g., albendazole; ALB), which bind ß-tubulin to disrupt microtubule polymerization, are one of two primary compound classes used to treat giardiasis. In most parasitic nematodes and fungi, BZ-resistance is caused by ß-tubulin mutations and its molecular mode of action (MOA) is well studied. In contrast, in Giardia duodenalis BZ MOA or resistance is less well understood, may involve target-specific and broader impacts including cellular damage and oxidative stress, and its underlying cause is not clearly determined. Previously, we identified acquisition of a single nucleotide polymorphism, E198K, in ß-tubulin in ALB-resistant (ALB-R) G. duodenalis WB-1B relative to ALB-sensitive (ALB-S) parental controls. E198K is linked to BZ-resistance in fungi and its allelic frequency correlated with the magnitude of BZ-resistance in G. duodenalis WB-1B. Here, we undertook detailed transcriptomic comparisons of these ALB-S and ALB-R G. duodenalis WB-1B cultures. The primary transcriptional changes with ALB-R in G. duodenalis WB-1B indicated increased protein degradation and turnover, and up-regulation of tubulin, and related genes, associated with the adhesive disc and basal bodies. These findings are consistent with previous observations noting focused disintegration of the disc and associated structures in Giardia duodenalis upon ALB exposure. We also saw transcriptional changes with ALB-R in G. duodenalis WB-1B consistent with prior observations of a shift from glycolysis to arginine metabolism for ATP production and possible changes to aspects of the vesicular trafficking system that require further investigation. Finally, we saw mixed transcriptional changes associated with DNA repair and oxidative stress responses in the G. duodenalis WB-1B line. These changes may be indicative of a role for H2O2 degradation in ALB-R, as has been observed in other G. duodenalis cell cultures. However, they were below the transcriptional fold-change threshold (log2FC > 1) typically employed in transcriptomic analyses and appear to be contradicted in ALB-R G. duodenalis WB-1B by down-regulation of the NAD scavenging and conversion pathways required to support these stress pathways and up-regulation of many highly oxidation sensitive iron-sulphur (FeS) cluster based metabolic enzymes.

Conservation, abundance, glycosylation profile, and localization of the TSP protein family in Cryptosporidium parvum.

Cryptosporidium parvum is a zoonotic apicomplexan parasite and a common cause of diarrheal disease worldwide. The development of vaccines to prevent or limit infection remains an important goal for tackling cryptosporidiosis. At present, the only approved vaccine against any apicomplexan parasite targets a conserved adhesin possessing a thrombospondin repeat domain. C. parvum possesses 12 orthologous thrombospondin repeat domain-containing proteins known as CpTSP1-12, though little is known about these potentially important antigens. Here, we explore the architecture and conservation of the CpTSP protein family, as well as their abundance at the protein level within the sporozoite stage of the life cycle. We examine the glycosylation states of these proteins using a combination of glycopeptide enrichment techniques to demonstrate that these proteins are modified with C-, O-, and N-linked glycans. Using expansion microscopy, and an antibody against the C-linked mannose that is unique to the CpTSP protein family within C. parvum, we show that these proteins are found both on the cell surface and in structures that resemble the secretory pathway of C. parvum sporozoites. Finally, we generated a polyclonal antibody against CpTSP1 to show that it is found at the cell surface and within micronemes, in a pattern reminiscent of other apicomplexan motility-associated adhesins, and is present both in sporozoites and meronts. This work sheds new light on an understudied family of C. parvum proteins that are likely to be important to both parasite biology and the development of vaccines against cryptosporidiosis.

Recent genetic exchanges and admixture shape the genome and population structure of the zoonotic pathogen Cryptosporidium parvum.

Cryptosporidium parvum is a globally distributed zoonotic pathogen and a major cause of diarrhoeal disease in humans and ruminants. The parasite's life cycle comprises an obligatory sexual phase, during which genetic exchanges can occur between previously isolated lineages. Here, we compare 32 whole genome sequences from human- and ruminant-derived parasite isolates collected across Europe, Egypt and China. We identify three strongly supported clusters that comprise a mix of isolates from different host species, geographic origins, and subtypes. We show that: (1) recombination occurs between ruminant isolates into human isolates; (2) these recombinant regions can be passed on to other human subtypes through gene flow and population admixture; (3) there have been multiple genetic exchanges, and most are probably recent; (4) putative virulence genes are significantly enriched within these genetic exchanges, and (5) this results in an increase in their nucleotide diversity. We carefully dissect the phylogenetic sequence of two genetic exchanges, illustrating the long-term evolutionary consequences of these events. Our results suggest that increased globalization and close human-animal contacts increase the opportunity for genetic exchanges between previously isolated parasite lineages, resulting in spillover and spillback events. We discuss how this can provide a novel substrate for natural selection at genes involved in host-parasite interactions, thereby potentially altering the dynamic coevolutionary equilibrium in the Red Queens arms race.

VIVID: A Web Application for Variant Interpretation and Visualization in Multi-dimensional Analyses.

Large-scale comparative genomics- and population genetic studies generate enormous amounts of polymorphism data in the form of DNA variants. Ultimately, the goal of many of these studies is to associate genetic variants to phenotypes or fitness. We introduce VIVID, an interactive, user-friendly web application that integrates a wide range of approaches for encoding genotypic to phenotypic information in any organism or disease, from an individual or population, in three-dimensional (3D) space. It allows mutation mapping and annotation, calculation of interactions and conservation scores, prediction of harmful effects, analysis of diversity and selection, and 3D visualization of genotypic information encoded in Variant Call Format on AlphaFold2 protein models. VIVID enables the rapid assessment of genes of interest in the study of adaptive evolution and the genetic load, and it helps prioritizing targets for experimental validation. We demonstrate the utility of VIVID by exploring the evolutionary genetics of the parasitic protist Plasmodium falciparum, revealing geographic variation in the signature of balancing selection in potential targets of functional antibodies.

Global Population Genomics of Two Subspecies of Cryptosporidium hominis during 500 Years of Evolution.

Cryptosporidiosis is a major global health problem and a primary cause of diarrhea, particularly in young children in low- and middle-income countries (LMICs). The zoonotic Cryptosporidium parvum and anthroponotic Cryptosporidium hominis cause most human infections. Here, we present a comprehensive whole-genome study of C. hominis, comprising 114 isolates from 16 countries within five continents. We detect two lineages with distinct biology and demography, which diverged circa 500 years ago. We consider these lineages two subspecies and propose the names C. hominis hominis and C. hominis aquapotentis (gp60 subtype IbA10G2). In our study, C. h. hominis is almost exclusively represented by isolates from LMICs in Africa and Asia and appears to have undergone recent population contraction. In contrast, C. h. aquapotentis was found in high-income countries, mainly in Europe, North America, and Oceania, and appears to be expanding. Notably, C. h. aquapotentis is associated with high rates of direct human-to-human transmission, which may explain its success in countries with well-developed environmental sanitation infrastructure. Intriguingly, we detected genomic regions of introgression following secondary contact between the subspecies. This resulted in high diversity and divergence in genomic islands of putative virulence genes, including muc5 (CHUDEA2_430) and a hypothetical protein (CHUDEA6_5270). This diversity is maintained by balancing selection, suggesting a co-evolutionary arms race with the host. Finally, we find that recent gene flow from C. h. aquapotentis to C. h. hominis, likely associated with increased human migration, maybe driving the evolution of more virulent C. hominis variants.

In vitro selection of Giardia duodenalis for Albendazole resistance identifies a ß-tubulin mutation at amino acid E198K.

Benzimidazole-2-carbamate (BZ) compounds, including Albendazole (Alb), are one of just two drug classes approved to treat the gastrointestinal protist Giardia duodenalis. Benzimidazoles bind to the tubulin dimer interface overlapping the colchicine binding site (CBS) of ß-tubulin, thereby inhibiting microtubule polymerisation and disrupting microtubule networks. These BZ compounds are widely used as anthelmintic, anti-fungal and anti-giardial drugs. However, in helminths and fungi, BZ-resistance is widespread and caused by specific point mutations primarily occurring at F167, E198 and F200 in ß-tubulin isoform 1. BZ-resistance in Giardia is reported clinically and readily generated in vitro, with significant implications for Giardia control. In Giardia, BZ mode of action (MOA) and resistance mechanisms are presumed but not proven, and no mutations in ß-tubulin have been reported in association with Alb resistance (AlbR). Herein, we undertook detailed in vitro drug-susceptibility screens of 13 BZ compounds and 7 Alb structural analogues in isogenic G. duodenalis isolates selected for AlbR and podophyllotoxin, another ß-tubulin inhibitor, as well as explored cross-resistance to structurally unrelated, metronidazole (Mtz). AlbR lines exhibited co-resistance to many structural variants in the BZ-pharmacophore, and cross-resistance to podophyllotoxin. AlbR lines were not cross-resistant to Mtz, but MtzR lines had enhanced survival in Alb. Lastly, Alb analogues with longer thioether substituents had decreased potency against our AlbR lines. In silico modelling indicated the Alb-ß-tubulin interaction in Giardia partially overlaps the CBS and corresponds to residues associated with BZ-resistance in helminths and fungi (F167, E198, F200). Sequencing of Giardia ß-tubulin identified a single nucleotide polymorphism resulting in a mutation from glutamic acid to lysine at amino acid 198 (E198K). To our knowledge, this is the first ß-tubulin mutation reported for protistan BZ-resistance. This study provides insight into BZ mode of action and resistance in Giardia, and presents a potential avenue for a genetic test for clinically resistance isolates.

Multimodal regulation of encystation in Giardia duodenalis revealed by deep proteomics.

Cyst formation in the parasitic protist Giardia duodenalis is critical to its transmission. Existing proteomic data quantifies only 17% of coding genes transcribed during encystation and does not cover the complete process from trophozoite to mature cyst. Using high-resolution mass spectrometry, we have quantified proteomic changes across encystation and compared this with published transcriptomic data. We reproducibly identified 3863 (64.5% of Giardia proteins) and quantified 3382 proteins (56.5% of Giardia proteins) over standard trophozoite growth (TY), during low-bile encystation priming (LB), 16 h into encystation (EC), and at cyst maturation (C). This work provides the first known expanded observation of encystation at the proteomic level and triples the coverage of previous encystation proteomes. One-third (1169 proteins) of the quantified proteome is differentially expressed in the mature cyst relative to the trophozoite, including proteasomal machinery, metabolic pathways, and secretory proteins. Changes in lipid metabolism indicated a shift in lipid species dependency during encystation. Consistent with this, we identified the first, putative lipid transporters in this species, representing the steroidogenic acute regulatory protein-related lipid transfer (StARkin), oxysterol binding protein related protein (ORP/Osh) and glycosphingolipid transfer protein (GLTP) families, and follow their differential expression over cyst formation. Lastly, we undertook correlation analyses of the transcriptome and proteome of trophozoites and cysts, and found evidence of post-transcriptional regulation of key protein classes (RNA binding proteins) and stage-specific genes (encystation markers) implicating translation-repression in encystation. We provide the most extensive proteomic analysis of encystation in Giardia to date and the first known exploration across its complete duration. This work identifies encystation as highly coordinated, involving major changes in proteostasis, metabolism and membrane dynamics, and indicates a potential role for post-transcriptional regulation, mediated through RNA-binding proteins. Together our work provides a valuable resource for Giardia research and the development of transmission-blocking anti-giardials.

Comparative genomics revealed adaptive admixture in Cryptosporidium hominis in Africa.

Cryptosporidiosis is a major cause of diarrhoeal illness among African children, and is associated with childhood mortality, malnutrition, cognitive development and growth retardation. Cryptosporidium hominis is the dominant pathogen in Africa, and genotyping at the glycoprotein 60 (gp60) gene has revealed a complex distribution of different subtypes across this continent. However, a comprehensive exploration of the metapopulation structure and evolution based on whole-genome data has yet to be performed. Here, we sequenced and analysed the genomes of 26 C. hominis isolates, representing different gp60 subtypes, collected at rural sites in Gabon, Ghana, Madagascar and Tanzania. Phylogenetic and cluster analyses based on single-nucleotide polymorphisms showed that isolates predominantly clustered by their country of origin, irrespective of their gp60 subtype. We found a significant isolation-by-distance signature that shows the importance of local transmission, but we also detected evidence of hybridization between isolates of different geographical regions. We identified 37 outlier genes with exceptionally high nucleotide diversity, and this group is significantly enriched for genes encoding extracellular proteins and signal peptides. Furthermore, these genes are found more often than expected in recombinant regions, and they show a distinct signature of positive or balancing selection. We conclude that: (1) the metapopulation structure of C. hominis can only be accurately captured by whole-genome analyses; (2) local anthroponotic transmission underpins the spread of this pathogen in Africa; (3) hybridization occurs between distinct geographical lineages; and (4) genetic introgression provides novel substrate for positive or balancing selection in genes involved in host-parasite coevolution.

Eukaryote-Conserved Methylarginine Is Absent in Diplomonads and Functionally Compensated in Giardia.

Methylation is a common posttranslational modification of arginine and lysine in eukaryotic proteins. Methylproteomes are best characterized for higher eukaryotes, where they are functionally expanded and evolved complex regulation. However, this is not the case for protist species evolved from the earliest eukaryotic lineages. Here, we integrated bioinformatic, proteomic, and drug-screening data sets to comprehensively explore the methylproteome of Giardia duodenalis-a deeply branching parasitic protist. We demonstrate that Giardia and related diplomonads lack arginine-methyltransferases and have remodeled conserved RGG/RG motifs targeted by these enzymes. We also provide experimental evidence for methylarginine absence in proteomes of Giardia but readily detect methyllysine. We bioinformatically infer 11 lysine-methyltransferases in Giardia, including highly diverged Su(var)3-9, Enhancer-of-zeste and Trithorax proteins with reduced domain architectures, and novel annotations demonstrating conserved methyllysine regulation of eukaryotic elongation factor 1 alpha. Using mass spectrometry, we identify more than 200 methyllysine sites in Giardia, including in species-specific gene families involved in cytoskeletal regulation, enriched in coiled-coil features. Finally, we use known methylation inhibitors to show that methylation plays key roles in replication and cyst formation in this parasite. This study highlights reduced methylation enzymes, sites, and functions early in eukaryote evolution, including absent methylarginine networks in the Diplomonadida. These results challenge the view that arginine methylation is eukaryote conserved and demonstrate that functional compensation of methylarginine was possible preceding expansion and diversification of these key networks in higher eukaryotes.

Disulphide bond restrains the C-terminal region of thermostable direct hemolysin during folding to promote oligomerization.

Pore-forming toxins (PFTs) are typically produced as water-soluble monomers, which upon interacting with target cells assemble into transmembrane oligomeric pores. Vibrio parahaemolyticus thermostable direct hemolysin (TDH) is an atypical PFT that exists as a tetramer in solution, prior to membrane binding. The TDH structure highlights a core ß-sandwich domain similar to those found in the eukaryotic actinoporin family of PFTs. However, the TDH structure harbors an extended C-terminal region (CTR) that is not documented in the actinoporins. This CTR remains tethered to the ß-sandwich domain through an intra-molecular disulphide bond. Part of the CTR is positioned at the inter-protomer interface in the TDH tetramer. Here we show that the truncation, as well as mutation, of the CTR compromise tetrameric assembly, and the membrane-damaging activity of TDH. Our study also reveals that intra-protomer disulphide bond formation during the folding/assembly process of TDH restrains the CTR to mediate its participation in the formation of inter-protomer contact, thus facilitating TDH oligomerization. However, once tetramerization is achieved, disruption of the disulphide bond does not affect oligomeric assembly. Our study provides critical insights regarding the regulation of the oligomerization mechanism of TDH, which has not been previously documented in the PFT family.

Effect of STAT3 inhibitor in chronic myeloid leukemia associated signaling pathway: a mathematical modeling, simulation and systems biology study.

Chronic myeloid leukemia (CML) is a hematopoietic stem-cell disorder which proliferates due to abnormal growth of basophil cells. Several proangiogenic molecules have been reported to be associated in CML progression, including the hepatocyte growth factor (HGF). However, detail mechanism about the cellular distribution and function of HGF in CML is yet to be revealed. The proliferation of hematopoietic cells are regulated by some of the growth factors like interleukin 3 (IL-3), IL-6, erythropoietin, thrombopoietin, etc. In this study IL-6 pathways have been taken into consideration which induces JAK/STAT and MAPK pathways to decipher the CML progression stages. An attempt has been made to model these pathways with the help of ordinary differential equations (ODEs) and estimating unknown parameters through fminsearch optimization algorithm. Some of the specific component like STAT3, of the pathway has been analyzed in detail and their role in CML progression has been elucidated. The roles of STAT3 inhibitors into the treatment of CML have been thoroughly studied and optimum concentration of the inhibitors have been predicted.

Parameters Involved in Autophosphorylation in Chronic Myeloid Leukemia: a Systems Biology Approach.

BACKGROUND: Chronic myeloid leukemia (CML) is a stem cell disorder characterized by the fusion of two oncogenes namely BCR and ABL with their aberrant expression. Autophosphorylation of BCR-ABL oncogenes results in proliferation of CML. The study deals with estimation of rate constant involved in each step of the cellular autophosphorylation process, which are consequently playing important roles in the proliferation of cancerous cells. MATERIALS AND METHODS: A mathematical model was proposed for autophosphorylation of BCR-ABL oncogenes utilizing ordinary differential equations to enumerate the rate of change of each responsible system component. The major difficulty to model this process is the lack of experimental data, which are needed to estimate unknown model parameters. Initial concentration data of each substrate and product for BCR-ABL systems were collected from the reported literature. All parameters were optimized through time interval simulation using the fminsearch algorithm. RESULTS: The rate of change versus time was estimated to indicate the role of each state variable that are crucial for the systems. The time wise change in concentration of substrate shows the convergence of each parameter in autophosphorylation process. CONCLUSIONS: The role of each constituent parameter and their relative time dependent variations in autophosphorylation process could be inferred.