Multiomics analysis reveals B. MO1 as a distinct Babesia species and provides insights into its evolution and virulence.

Babesiosis, caused by protozoan parasites of the genus Babesia , is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of various Babesia species underscores the ongoing risk of new zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental shifts impacting the distribution and transmission dynamics of parasites, their vectors, and reservoir hosts. One such species, Babesia MO1, previously implicated in severe cases of human babesiosis in the midwestern United States, was initially considered closely related to B. divergens , the predominant agent of human babesiosis in Europe. Yet, uncertainties persist regarding whether these pathogens represent distinct variants of the same species or are entirely separate species. We show that although both B. MO1 and B. divergens share similar genome sizes, comprising three nuclear chromosomes, one linear mitochondrial chromosome, and one circular apicoplast chromosome, major differences exist in terms of genomic sequence divergence, gene functions, transcription profiles, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens , and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

Inferring ligand-receptor cellular networks from bulk and spatial transcriptomic datasets with BulkSignalR.

The study of cellular networks mediated by ligand-receptor interactions has attracted much attention recently owing to single-cell omics. However, rich collections of bulk data accompanied with clinical information exists and continue to be generated with no equivalent in single-cell so far. In parallel, spatial transcriptomic (ST) analyses represent a revolutionary tool in biology. A large number of ST projects rely on multicellular resolution, for instance the Visium™ platform, where several cells are analyzed at each location, thus producing localized bulk data. Here, we describe BulkSignalR, a R package to infer ligand-receptor networks from bulk data. BulkSignalR integrates ligand-receptor interactions with downstream pathways to estimate statistical significance. A range of visualization methods complement the statistics, including functions dedicated to spatial data. We demonstrate BulkSignalR relevance using different datasets, including new Visium liver metastasis ST data, with experimental validation of protein colocalization. A comparison with other ST packages shows the significantly higher quality of BulkSignalR inferences. BulkSignalR can be applied to any species thanks to its built-in generic ortholog mapping functionality.

Babesia duncani multi-omics identifies virulence factors and drug targets.

Babesiosis is a malaria-like disease in humans and animals that is caused by Babesia species, which are tick-transmitted apicomplexan pathogens. Babesia duncani causes severe to lethal infection in humans, but despite the risk that this parasite poses as an emerging pathogen, little is known about its biology, metabolic requirements or pathogenesis. Unlike other apicomplexan parasites that infect red blood cells, B. duncani can be continuously cultured in vitro in human erythrocytes and can infect mice resulting in fulminant babesiosis and death. We report comprehensive, detailed molecular, genomic, transcriptomic and epigenetic analyses to gain insights into the biology of B. duncani. We completed the assembly, 3D structure and annotation of its nuclear genome, and analysed its transcriptomic and epigenetics profiles during its asexual life cycle stages in human erythrocytes. We used RNA-seq data to produce an atlas of parasite metabolism during its intraerythrocytic life cycle. Characterization of the B. duncani genome, epigenome and transcriptome identified classes of candidate virulence factors, antigens for diagnosis of active infection and several attractive drug targets. Furthermore, metabolic reconstitutions from genome annotation and in vitro efficacy studies identified antifolates, pyrimethamine and WR-99210 as potent inhibitors of B. duncani to establish a pipeline of small molecules that could be developed as effective therapies for the treatment of human babesiosis.

FT-GPI, a highly sensitive and accurate predictor of GPI-anchored proteins, reveals the composition and evolution of the GPI proteome in Plasmodium species.

BACKGROUND: Protozoan parasites are known to attach specific and diverse group of proteins to their plasma membrane via a GPI anchor. In malaria parasites, GPI-anchored proteins (GPI-APs) have been shown to play an important role in host-pathogen interactions and a key function in host cell invasion and immune evasion. Because of their immunogenic properties, some of these proteins have been considered as malaria vaccine candidates. However, identification of all possible GPI-APs encoded by these parasites remains challenging due to their sequence diversity and limitations of the tools used for their characterization. METHODS: The FT-GPI software was developed to detect GPI-APs based on the presence of a hydrophobic helix at both ends of the premature peptide. FT-GPI was implemented in C ++and applied to study the GPI-proteome of 46 isolates of the order Haemosporida. Using the GPI proteome of Plasmodium falciparum strain 3D7 and Plasmodium vivax strain Sal-1, a heuristic method was defined to select the most sensitive and specific FT-GPI software parameters. RESULTS: FT-GPI enabled revision of the GPI-proteome of P. falciparum and P. vivax, including the identification of novel GPI-APs. Orthology- and synteny-based analyses showed that 19 of the 37 GPI-APs found in the order Haemosporida are conserved among Plasmodium species. Our analyses suggest that gene duplication and deletion events may have contributed significantly to the evolution of the GPI proteome, and its composition correlates with speciation. CONCLUSION: FT-GPI-based prediction is a useful tool for mining GPI-APs and gaining further insights into their evolution and sequence diversity. This resource may also help identify new protein candidates for the development of vaccines for malaria and other parasitic diseases.

Targeting the splicing isoforms of spleen tyrosine kinase affects the viability of colorectal cancer cells.

Spleen tyrosine kinase (Syk) expression have been both positively and negatively associated with tumorigenesis. Our goal was to evaluate the contribution of Syk and its two splice variants, full length Syk (L) and short isoform Syk (S), in the tumor biology of colorectal cancer cells (CRC). The analysis of Syk expression in primary human colorectal tumors, as well as the analysis of TCGA database, revealed a high Syk mRNA expression score in colorectal cancer tumors, suggesting a tumor promotor role of Syk in CRC. Our analysis showed that Syk (L) isoform is highly expressed in the majority of the tumor tissues and that it remains expressed in tumors in which global Syk expression is downregulated, suggesting the dependence of tumors to Syk (L) isoform. We also identified a small cluster of tumor tissues, which express a high proportion of Syk (S) isoform. This specific cluster is associated with overexpressed genes related to translation and mitochondria, and down regulated genes implicated in the progression of mitosis. For our functional studies, we used short hairpin RNA tools to target the expression of Syk in CRC cells bearing the activating K-Ras (G13D) mutation. Our results showed that while global Syk knock down increases cell proliferation and cell motility, Syk (L) expression silencing affects the viability and induces the apoptosis of the cells, confirming the dependence of cells on Syk (L) isoform for their survival. Finally, we report the promising potential of compound C-13, an original non-enzymatic inhibitor of Syk isolated in our group. In vitro studies showed that C-13 exerts cytotoxic effects on Syk-positive CRC cells by inhibiting their proliferation and their motility, and by inducing their apoptosis, while Syk-negative cell lines viability was not affected. Moreover, the oral and intraperitoneal administration of C-13 reduced the tumor growth of CRC DLD-1 cells xenografts in Nude mice in vivo.

Chemoradiation triggers antitumor Th1 and tissue resident memory-polarized immune responses to improve immune checkpoint inhibitors therapy.

BACKGROUND: Multiple synergistic combination approaches with cancer drugs are developed to overcome primary resistance to immunotherapy; however, the mechanistic rationale to combine chemoradiotherapy (CRT) with immune checkpoint inhibitors remains elusive. METHODS: This study described the immunological landscape of tumor microenvironment (TME) exposed to CRT. Tumor samples from patients with rectal cancer (n=43) treated with neoadjuvant CRT or radiotherapy were analyzed by nanostring and immunohistochemistry. Studies in mice were performed using three syngeneic tumors (TC1, CT26 and MC38). Tumor-bearing mice were treated either with platinum-based CRT, radiotherapy or chemotherapy. Anti-CTLA-4 and/or anti-Programmed Cell Death Receptor-1 (PD-1) therapy was used in combination with CRT. The therapy-exposed TME was screened by RNA sequencing and flow cytometry and tumor-infiltrating T lymphocyte functionality was evaluated by interferon (IFN)-γ ELIspot and intracellular cytokine staining. RESULTS: Front-to-front comparison analysis revealed the synergistic effect of CRT to establish a highly inflamed and Th1-polarized immune signature in the TME of patients and mice. In both settings, CRT-exposed TMEs were highly enriched in newly-infiltrated tumor-specific CD8+ T cells as well as tissue resident memory CD103+CD8+ T cells. In mice, CD8 T cells were involved in the antitumor response mediated by CRT and were primed by CRT-activated CD103+ dendritic cells. In the three tumor models, we showed that concurrent combination of CRT with a dual CTLA-4 and PD-1 blockade was required to achieve an optimal antitumor effect and to establish a broad and long-lasting protective antitumor T cell immunity. CONCLUSIONS: Our results highlight the ability of CRT to stimulate strong antitumor T-cell-mediated immunity and tissue resident memory T activation in TME, to foster immune checkpoint inhibitors action. These findings have implications in clinic for the design clinical trials combining chemoradiation with immunotherapy.

The immune contexture of primary central nervous system diffuse large B cell lymphoma associates with patient survival and specific cell signaling.

Rationale: Primary central nervous system diffuse large B-cell lymphoma (PCNSL) is a rare and aggressive entity that resides in an immune-privileged site. The tumor microenvironment (TME) and the disruption of the immune surveillance influence lymphoma pathogenesis and immunotherapy resistance. Despite growing knowledge on heterogeneous therapeutic responses, no comprehensive description of the PCNSL TME is available. We hence investigated the immune subtypes of PCNSL and their association with molecular signaling and survival. Methods: Analysis of PCNSL transcriptomes (sequencing, n = 20; microarrays, n = 34). Integrated correlation analysis and signaling pathway topology enabled us to infer intercellular interactions. Immunohistopathology and digital imaging were used to validate bioinformatic results. Results: Transcriptomics revealed three immune subtypes: immune-rich, poor, and intermediate. The immune-rich subtype was associated to better survival and characterized by hyper-activation of STAT3 signaling and inflammatory signaling, e.g., IFNγ and TNF-α, resembling the hot subtype described in primary testicular lymphoma and solid cancer. WNT/ß-catenin, HIPPO, and NOTCH signaling were hyper-activated in the immune-poor subtype. HLA down-modulation was clearly associated with a low or intermediate immune infiltration and the absence of T-cell activation. Moreover, HLA class I down-regulation was also correlated with worse survival with implications on immune-intermediate PCNSL that frequently feature reduced HLA expression. A ligand-receptor intercellular network revealed high expression of two immune checkpoints, i.e., CTLA-4/CD86 and TIM-3/LAGLS9. TIM-3 and galectin-9 proteins were clearly upregulated in PCNSL. Conclusion: Altogether, our study reveals that patient stratification according to immune subtypes, HLA status, and immune checkpoint molecule quantification should be considered prior to immune checkpoint inhibitor therapy. Moreover, TIM-3 protein should be considered an axis for future therapeutic development.

The molecular landscape and microenvironment of salivary duct carcinoma reveal new therapeutic opportunities.

Purpose: Salivary duct carcinoma (SDC) is a rare and aggressive salivary gland cancer subtype with poor prognosis. The mutational landscape of SDC has already been the object of several studies, however little is known regarding the functional genomics and the tumor microenvironment despite their importance in oncology. Our investigation aimed at describing both the functional genomics of SDC and the SDC microenvironment, along with their clinical relevance. Methods: RNA-sequencing (24 tumors), proteomics (17 tumors), immunohistochemistry (22 tumors), and multiplexed immunofluorescence (3 tumors) data were obtained from three different patient cohorts and analyzed by digital imaging and bioinformatics. Adjacent non-tumoral tissue from patients in two cohorts were used in transcriptomic and proteomic analyses. Results: Transcriptomic and proteomic data revealed the importance of Notch, TGF-ß, and interferon-γ signaling for all SDCs. We confirmed an overall strong desmoplastic reaction by measuring α-SMA abundance, the level of which was associated with recurrence-free survival (RFS). Two distinct immune phenotypes were observed: immune-poor SDCs (36%) and immune-infiltrated SDCs (64%). Advanced bioinformatics analysis of the transcriptomic data suggested 72 ligand-receptor interactions occurred in the microenvironment and correlated with the immune phenotype. Among these interactions, three immune checkpoints were validated by immunofluorescence, including CTLA-4/DC86 and TIM-3/galectin-9 interactions, previously unidentified in SDC. Immunofluorescence analysis also confirmed an important immunosuppressive role of macrophages and NK cells, also supported by the transcriptomic data. Conclusions: Together our data significantly increase the understanding of SDC biology and open new perspectives for SDC tumor treatment. Before applying immunotherapy, patient stratification according to the immune infiltrate should be taken into account. Immune-infiltrated SDC could benefit from immune checkpoint-targeting therapy, with novel options such as anti-CTLA-4. Macrophages or NK cells could also be targeted. The dense stroma, i.e., fibroblasts or hyaluronic acid, may also be the focus for immune-poor SDC therapies, e.g. in combination with Notch or TGF-ß inhibitors, or molecules targeting SDC mutations.

In Vivo Large-Scale Mapping of Protein Turnover in Human Cerebrospinal Fluid.

The extraction of accurate physiological parameters from clinical samples provides a unique perspective to understand disease etiology and evolution, including under therapy. We introduce a new methodologic framework to map patient proteome dynamics in vivo, either proteome-wide or in large targeted panels. We applied it to ventricular cerebrospinal fluid (CSF) and could determine the turnover parameters of almost 200 proteins, whereas a handful were known previously. We covered a large number of neuron biology- and immune system-related proteins, including many biomarkers and drug targets. This first large data set unraveled a significant relationship between turnover and protein origin that relates to our ability to investigate organ physiology with protein-labeling strategy specifics. Our data constitute the first draft of CSF proteome dynamics as well as a repertoire of peptides for the community to design new analyses. The disclosed methods apply to other fluids or tissues provided sequential sample collection can be performed. We show that the proposed mathematical modeling applies to other analytical methods in the field.

Babesia divergens glycosylphosphatidylinositols modulate blood coagulation and induce Th2-biased cytokine profiles in antigen presenting cells.

Glycosylphosphatidylinositols (GPIs) are glycolipids described as toxins of protozoan parasites due to their inflammatory properties in mammalian hosts characterized by the production of interleukin (IL)-1, IL-12 and tumor necrosis factor (TNF)-α. In the present work, we studied the cytokines produced by antigen presenting cells in response to ten different GPI species extracted from Babesia divergens, responsible for babesiosis. Interestingly, B. divergens GPIs induced the production of anti-inflammatory cytokines (IL-2, IL-5) and of the regulatory cytokine IL-10 by macrophages and dendritic cells. In contrast to all protozoan GPIs studied until now, GPIs from B. divergens did not stimulate the production of TNF-α and IL-12, leading to a unique Th1/Th2 profile. Analysis of the carbohydrate composition of the B. divergens GPIs indicated that the di-mannose structure was different from the evolutionary conserved tri-mannose structure, which might explain the particular cytokine profile they induce. Expression of major histocompatibility complex (MHC) molecules on dendritic cells and apoptosis of mouse peritoneal cells were also analysed. B. divergens GPIs did not change expression of MHC class I, but decreased expression of MHC class II at the cell surface, while GPIs slightly increased the percentages of apoptotic cells. During pathogenesis of babesiosis, the inflammation-coagulation auto-amplification loop can lead to thrombosis and the effect of GPIs on coagulation parameters was investigated. Incubation of B. divergens GPIs with rat plasma ex vivo led to increase of fibrinogen levels and to prolonged activated partial thromboplastin time, suggesting a direct modulation of the extrinsic coagulation pathway by GPIs.

Case report of the patient source of the Babesia microti R1 reference strain and implications for travelers.

Background: In 2002, a previously healthy 69-year-old man travelled to France from the United States and presented to our hospital with a febrile illness that subsequently was determined to be babesiosis. The blood isolated from this patient served as a source for propagation of the Babesia microti R1 strain with subsequent sequencing and annotation of the parasite genome. Methods: Upon admission, we obtained a medical history, performed a physical examination, and examined his blood for the presence of a blood borne pathogen by microscopy, PCR and indirect immunofluorescence antibody testing. Once the diagnosis of babesiosis was made, we reviewed the literature to assess the distribution of B. microti-associated babesiosis cases in immunocompetent patients from outside the USA. Results: The patient recalled a tick bite during the previous month on Cape Cod, Massachusetts. The diagnosis was confirmed by identification of Babesia-infected red blood cells on blood smears, amplification of B. microti DNA in blood by PCR and the presence of B. microti antibody in the serum. This strain was the first isolate of B. microti to be fully sequenced and its annotated genome serves as a reference for molecular and cell biology studies aimed at understanding B. microti pathophysiology and developing diagnostic tests and therapies. A review of babesiosis cases demonstrates a worldwide distribution of B. microti and identifies potential emerging endemic areas where travelers may be at risk of contracting B. microti infection. Conclusion: This case provides clinical information about the patient infected with the R1 isolate and a review of travel risk, diagnosis and treatment of babesiosis in endemic and non-endemic areas.

Correction: Human babesiosis: Indication of a molecular mimicry between thrombospondin domains from a novel Babesia microti BmP53 protein and host platelets molecules.

[This corrects the article DOI: 10.1371/journal.pone.0185372.].

Functional analysis of Plasmodium falciparum subpopulations associated with artemisinin resistance in Cambodia.

BACKGROUND: Plasmodium falciparum malaria is one of the most widespread parasitic infections in humans and remains a leading global health concern. Malaria elimination efforts are threatened by the emergence and spread of resistance to artemisinin-based combination therapy, the first-line treatment of malaria. Promising molecular markers and pathways associated with artemisinin drug resistance have been identified, but the underlying molecular mechanisms of resistance remains unknown. The genomic data from early period of emergence of artemisinin resistance (2008-2011) was evaluated, with aim to define k13 associated genetic background in Cambodia, the country identified as epicentre of anti-malarial drug resistance, through characterization of 167 parasite isolates using a panel of 21,257 SNPs. RESULTS: Eight subpopulations were identified suggesting a process of acquisition of artemisinin resistance consistent with an emergence-selection-diffusion model, supported by the shifting balance theory. Identification of population specific mutations facilitated the characterization of a core set of 57 background genes associated with artemisinin resistance and associated pathways. The analysis indicates that the background of artemisinin resistance was not acquired after drug pressure, rather is the result of fixation followed by selection on the daughter subpopulations derived from the ancestral population. CONCLUSIONS: Functional analysis of artemisinin resistance subpopulations illustrates the strong interplay between ubiquitination and cell division or differentiation in artemisinin resistant parasites. The relationship of these pathways with the P. falciparum resistant subpopulation and presence of drug resistance markers in addition to k13, highlights the major role of admixed parasite population in the diffusion of artemisinin resistant background. The diffusion of resistant genes in the Cambodian admixed population after selection resulted from mating of gametocytes of sensitive and resistant parasite populations.

Incidence of dengue and chikungunya viruses in mosquitoes and human patients in border provinces of Vietnam.

BACKGROUND: Dengue virus remains a major threat in Vietnam, while chikungunya virus is expected to become one. Surveillance was conducted from 2012 to 2014 in Vietnam to assess the presence of dengue and chikungunya viruses in patients hospitalized with acute fever in five Vietnam provinces neighboring Lao PDR and Cambodia. Surveillance was extended to mosquitoes present in the vicinity of the patients' households. RESULTS: A total 558 human serum samples were collected along with 1104 adult mosquitoes and 12,041 larvae from 2250 households. Dengue virus was found in 17 (3%) human serum samples and in 9 (0.8%) adult mosquitoes. Chikungunya virus was detected in 2 adult mosquitoes (0.18%) while no chikungunya virus was detected in humans. Differing densities of mosquito populations were found, with the highest in the Long An Province border with Cambodia. Long An Province also displayed the lowest rate of infection, despite a very high Breteau Index, high human population density and presence of the main cross border road system. The highest incidence was found in Dac Nong Province, where the Breteau and Container indices were the second lowest. Dengue virus was detected in five Aedes albopictus, three Aedes aegypti and one Culex vishnui. Chikungunya virus was detected in two Ae. aegypti. All infected mosquitoes belonged to haplotypes described in other parts of the world and a number of novel haplotypes were found among uninfected mosquitoes. CONCLUSIONS: Dengue is considered to be regularly introduced to Vietnam from Cambodia, mostly through human movement. The data reported here provides a complementary picture. Due to intensive international trade, long-distance transportation of mosquito populations may play a role in the regular importation of dengue in Vietnam through Ho Chi Minh City. It is important to decipher the movement of mosquitoes in Vietnam, not only at the Lao PDR and Cambodia borders but also through international trade routes. Mosquito surveillance programs should address and follow mosquito populations instead of mosquito species.

Human babesiosis: Indication of a molecular mimicry between thrombospondin domains from a novel Babesia microti BmP53 protein and host platelets molecules.

Human babesiosis is caused by the apicomplexan parasite Babesia microti, which is of major public health concern in the United States and elsewhere, resulting in malaise and fatigue, followed by a fever and hemolytic anemia. In this paper we focus on the characterization of a novel B. microti thrombospondin domain (TSP1)-containing protein (BmP53) from the new annotation of the B. microti genome (locus 'BmR1_04g09041'). This novel protein (BmP53) had a single TSP1 and a transmembrane domain, with a short cytoplasmic tail containing a sub-terminal glutamine residue, but no signal peptide and Von Willebrand factor type A domains (VWA), which are found in classical thrombospondin-related adhesive proteins (TRAP). Co-localization assays of BmP53 and Babesia microti secreted antigen 1 (BmSA1) suggested that BmP53 might be a non-secretory membranous protein. Molecular mimicry between the TSP1 domain from BmP53 and host platelets molecules was indicated through different measures of sequence homology, phylogenetic analysis, 3D structure and shared epitopes. Indeed, hamster isolated platelets cross-reacted with mouse anti-BmP53-TSP1. Molecular mimicry are used to help parasites to escape immune defenses, resulting in immune evasion or autoimmunity. Furthermore, specific host reactivity was also detected against the TSP1-free part of BmP53 in infected hamster sera. In conclusion, the TSP1 domain mimicry might help in studying the mechanisms of parasite-induced thrombocytopenia, with the TSP1-free truncate of the protein representing a potential safe candidate for future vaccine studies.

Valine/isoleucine variants drive selective pressure in the VP1 sequence of EV-A71 enteroviruses.

BACKGROUND: In 2011-2012, Northern Vietnam experienced its first large scale hand foot and mouth disease (HFMD) epidemic. In 2011, a major HFMD epidemic was also reported in South Vietnam with fatal cases. This 2011-2012 outbreak was the first one to occur in North Vietnam providing grounds to study the etiology, origin and dynamic of the disease. We report here the analysis of the VP1 gene of strains isolated throughout North Vietnam during the 2011-2012 outbreak and before. METHODS: The VP1 gene of 106 EV-A71 isolates from North Vietnam and 2 from Central Vietnam were sequenced. Sequence alignments were analyzed at the nucleic acid and protein level. Gene polymorphism was also analyzed. A Factorial Correspondence Analysis was performed to correlate amino acid mutations with clinical parameters. RESULTS: The sequences were distributed into four phylogenetic clusters. Three clusters corresponded to the subgenogroup C4 and the last one corresponded to the subgenogroup C5. Each cluster displayed different polymorphism characteristics. Proteins were highly conserved but three sites bearing only Isoleucine (I) or Valine (V) were characterized. The isoleucine/valine variability matched the clusters. Spatiotemporal analysis of the I/V variants showed that all variants which emerged in 2011 and then in 2012 were not the same but were all present in the region prior to the 2011-2012 outbreak. Some correlation was found between certain I/V variants and ethnicity and severity. CONCLUSIONS: The 2011-2012 outbreak was not caused by an exogenous strain coming from South Vietnam or elsewhere but by strains already present and circulating at low level in North Vietnam. However, what triggered the outbreak remains unclear. A selective pressure is applied on I/V variants which matches the genetic clusters. I/V variants were shown on other viruses to correlate with pathogenicity. This should be investigated in EV-A71. I/V variants are an easy and efficient way to survey and identify circulating EV-A71 strains.

Babesia microti from humans and ticks hold a genomic signature of strong population structure in the United States.

BACKGROUND: Babesia microti is an emerging tick-borne apicomplexan parasite with increasing geographic range and incidence in the United States. The rapid expansion of B. microti into its current distribution in the northeastern USA has been due to the range expansion of the tick vector, Ixodes scapularis, upon which the causative agent is dependent for transmission to humans. RESULTS: To reconstruct the history of B. microti in the continental USA and clarify the evolutionary origin of human strains, we used multiplexed hybrid capture of 25 B. microti isolates obtained from I. scapularis and human blood. Despite low genomic variation compared with other Apicomplexa, B. microti was strongly structured into three highly differentiated genetic clusters in the northeastern USA. Bayesian analyses of the apicoplast genomes suggest that the origin of the current diversity of B. microti in northeastern USA dates back 46 thousand years with a signature of recent population expansion in the last 1000 years. Human-derived samples belonged to two rarely intermixing clusters, raising the possibility of highly divergent infectious phenotypes in humans. CONCLUSIONS: Our results validate the multiplexed hybrid capture strategy for characterizing genome-wide diversity and relatedness of B. microti from ticks and humans. We find strong population structure in B. microti samples from the Northeast indicating potential barriers to gene flow.

Genome-wide diversity and gene expression profiling of Babesia microti isolates identify polymorphic genes that mediate host-pathogen interactions.

Babesia microti, a tick-transmitted, intraerythrocytic protozoan parasite circulating mainly among small mammals, is the primary cause of human babesiosis. While most cases are transmitted by Ixodes ticks, the disease may also be transmitted through blood transfusion and perinatally. A comprehensive analysis of genome composition, genetic diversity, and gene expression profiling of seven B. microti isolates revealed that genetic variation in isolates from the Northeast United States is almost exclusively associated with genes encoding the surface proteome and secretome of the parasite. Furthermore, we found that polymorphism is restricted to a small number of genes, which are highly expressed during infection. In order to identify pathogen-encoded factors involved in host-parasite interactions, we screened a proteome array comprised of 174 B. microti proteins, including several predicted members of the parasite secretome. Using this immuno-proteomic approach we identified several novel antigens that trigger strong host immune responses during the onset of infection. The genomic and immunological data presented herein provide the first insights into the determinants of B. microti interaction with its mammalian hosts and their relevance for understanding the selective pressures acting on parasite evolution.

Genetic Characterization of Plasmodium Putative Pantothenate Kinase Genes Reveals Their Essential Role in Malaria Parasite Transmission to the Mosquito.

The metabolic machinery for the biosynthesis of Coenzyme A (CoA) from exogenous pantothenic acid (Vitamin B5) has long been considered as an excellent target for the development of selective antimicrobials. Earlier studies in the human malaria parasite Plasmodium falciparum have shown that pantothenate analogs interfere with pantothenate phosphorylation and block asexual blood stage development. Although two eukaryotic-type putative pantothenate kinase genes (PanK1 and PanK2) have been identified in all malaria parasite species, their role in the development of Plasmodium life cycle stages remains unknown. Here we report on the genetic characterization of PanK1 and PanK2 in P. yoelii. We show that P. yoelii parasites lacking either PanK1 or PanK2 undergo normal asexual stages development and sexual stages differentiation, however they are severely deficient in ookinete, oocyst and sporozoite formation inside the mosquito vector. Quantitative transcriptional analyses in wild-type and knockout parasites demonstrate an important role for these genes in the regulation of expression of other CoA biosynthesis genes. Together, our data provide the first genetic evidence for the importance of the early steps of pantothenate utilization in the regulation of CoA biosynthesis and malaria parasite transmission to Anopheles mosquitoes.

Plasmodium falciparum parasite population structure and gene flow associated to anti-malarial drugs resistance in Cambodia.

BACKGROUND: Western Cambodia is recognized as the epicentre of emergence of Plasmodium falciparum multi-drug resistance. The emergence of artemisinin resistance has been observed in this area since 2008-2009 and molecular signatures associated to artemisinin resistance have been characterized in k13 gene. At present, one of the major threats faced, is the possible spread of Asian artemisinin resistant parasites over the world threatening millions of people and jeopardizing malaria elimination programme efforts. To anticipate the diffusion of artemisinin resistance, the identification of the P. falciparum population structure and the gene flow among the parasite population in Cambodia are essential. METHODS: To this end, a mid-throughput PCR-LDR-FMA approach based on LUMINEX technology was developed to screen for genetic barcode in 533 blood samples collected in 2010-2011 from 16 health centres in malaria endemics areas in Cambodia. RESULTS: Based on successful typing of 282 samples, subpopulations were characterized along the borders of the country. Each 11-loci barcode provides evidence supporting allele distribution gradient related to subpopulations and gene flow. The 11-loci barcode successfully identifies recently emerging parasite subpopulations in western Cambodia that are associated with the C580Y dominant allele for artemisinin resistance in k13 gene. A subpopulation was identified in northern Cambodia that was associated to artemisinin (R539T resistant allele of k13 gene) and mefloquine resistance. CONCLUSIONS: The gene flow between these subpopulations might have driven the spread of artemisinin resistance over Cambodia.