Single-cell topographical profiling of the immune synapse reveals a biomechanical signature of cytotoxicity.

Immune cells have intensely physical lifestyles characterized by structural plasticity and force exertion. To investigate whether specific immune functions require stereotyped mechanical outputs, we used super-resolution traction force microscopy to compare the immune synapses formed by cytotoxic T cells with contacts formed by other T cell subsets and by macrophages. T cell synapses were globally compressive, which was fundamentally different from the pulling and pinching associated with macrophage phagocytosis. Spectral decomposition of force exertion patterns from each cell type linked cytotoxicity to compressive strength, local protrusiveness, and the induction of complex, asymmetric topography. These features were validated as cytotoxic drivers by genetic disruption of cytoskeletal regulators, live imaging of synaptic secretion, and in silico analysis of interfacial distortion. Synapse architecture and force exertion were sensitive to target stiffness and size, suggesting that the mechanical potentiation of killing is biophysically adaptive. We conclude that cellular cytotoxicity and, by implication, other effector responses are supported by specialized patterns of efferent force.

The Toxoplasma gondii F-Box Protein L2 Functions as a Repressor of Stage Specific Gene Expression.

Toxoplasma gondii is a foodborne pathogen that can cause severe and life-threatening infections in fetuses and immunocompromised patients. Felids are its only definitive hosts, and a wide range of animals, including humans, serve as intermediate hosts. When the transmissible bradyzoite stage is orally ingested by felids, they transform into merozoites that expand asexually, ultimately generating millions of gametes for the parasite sexual cycle. However, bradyzoites in intermediate hosts differentiate exclusively to disease-causing tachyzoites, which rapidly disseminate throughout the host. Though tachyzoites are well-studied, the molecular mechanisms governing transitioning between developmental stages are poorly understood. Each parasite stage can be distinguished by a characteristic transcriptional signature, with one signature being repressed during the other stages. Switching between stages require substantial changes in the proteome, which is achieved in part by ubiquitination. F-box proteins mediate protein poly-ubiquitination by recruiting substrates to SKP1, Cullin-1, F-Box protein E3 ubiquitin ligase (SCF-E3) complexes. We have identified an F-box protein named Toxoplasma gondii F-Box Protein L2 (TgFBXL2), which localizes to distinct perinucleolar sites. TgFBXL2 is stably engaged in an SCF-E3 complex that is surprisingly also associated with a COP9 signalosome complex that negatively regulates SCF-E3 function. At the cellular level, TgFBXL2-depleted parasites are severely defective in centrosome replication and daughter cell development. Most remarkable, RNAseq data show that TgFBXL2 conditional depletion induces the expression of stage-specific genes including a large cohort of genes necessary for sexual commitment. Together, these data suggest that TgFBXL2 is a latent guardian of stage specific gene expression in Toxoplasma and poised to remove conflicting proteins in response to an unknown trigger of development.

Shaping membrane vesicles by adsorption of hinge-like nanoparticles.

The adsorption of particles onto fluid membranes can lead to membrane-mediated interactions between particles that promote their self-assembly and lead to changes in membrane morphology. However, in contrast with rigid particles, relatively little is known about deformable particles, which introduce additional complexities due to the mutual deformability of the particles and the membrane. Here, we use Monte Carlo simulations and umbrella sampling to investigate the equilibrium properties of hinge-like particles adsorbed on membrane vesicles by means of anisotropic, attractive interactions. We vary the hinge stiffness, adhesive area fraction, patterning of adhesive regions, and number of adsorbed particles. Depending on their properties, isolated particles can conform to the vesicle, induce invaginations of the membrane, or exhibit multistable behavior in which they sample distinct classes of configurations due to the interplay of particle and membrane deformations. With two adsorbed particles, the properties of the particles can be used to promote aggregation, bias the particles to different parts of the vesicle, or stabilize the coexistence of both cases. With multiple adsorbed particles, the number and type control their organization and collective impact on the vesicle, which can adopt shapes ranging from roughly spherical to dumbbell-like and multi-lobed. Our results highlight how modifying the mechanical properties and patterned adhesion of deformable particles, which is possible with DNA nanotechnology, influences their self-assembly and the resulting shapes of both the particles and vesicles.

Stereotactic Body Radiation Therapy (SBRT) for Hepatocellular Carcinoma (HCC) With Single Photon Emission Computed Tomography (SPECT) Functional Treatment Planning in Patients With Advanced Hepatic Cirrhosis.

Purpose: We report on the feasibility and outcomes of liver stereotactic body radiation therapy (SBRT) for hepatocellular carcinoma (HCC) with single-photon emission computed tomography (SPECT) functional treatment planning in patients with Child-Pugh (CP) B/C cirrhosis. Methods and Materials: Liver SPECT with 99mTc-sulfur colloid was coregistered to treatment planning computed tomography (CT) for the guided avoidance of functional hepatic parenchyma during SBRT. Functional liver volumes (FLVs) obtained from SPECT were compared with anatomic liver volumes defined on the planning CT. Radiation dose constraints were adapted exclusively to FLV. Local control, toxicity, and survival were reported with at least 6 months of radiographic follow-up. Pre- and posttransplant outcomes were analyzed in a subset of patients who completed SBRT as a bridge to liver transplant. Model of End-Stage Liver Disease was used to score hepatic function before and after SBRT completion. Results: With a median follow-up of 32 months, 45 patients (58 lesions) with HCC and CP-B/C cirrhosis received SBRT to a median dose of 45 Gy (3-5 fractions). FLV loss (34%, P < .001) was observed in all patients, and the functional and anatomic liver volumes matched well in a control group of noncirrhotic/non-HCC patients. Despite marked functional parenchyma retraction, the amount of FLV on SPECT exposed to the threshold irradiation was significantly less than the CT liver volumes (P < .001) because of the optimized beam placement during dosimetry planning. Twenty-three patients (51%) successfully completed orthotopic liver transplant, with a median time to transplant of 9.2 months. With 91% in-field local control, the overall 2-year survival was 65% (90% after the orthotopic liver transplant), with no incidence of radiation-induced liver disease observed within 3 to 4 months or accelerated CP class migration from B to C within the first 6 months post-SBRT. Mean Model of End-Stage Liver Disease-Na score was not significantly elevated at 3-month intervals after SBRT completion. Conclusions: Functional treatment planning with 99mTc sulfur colloid SPECT/CT allows identification and avoidance of functional hepatic parenchyma in patients with CP-B/C cirrhosis, leading to low toxicity and satisfactory transplant outcomes.

Compartmental exchange regulates steady states and stochastic switching of a phosphorylation network.

The phosphoregulation of proteins with multiple phosphorylation sites is governed by biochemical reaction networks that can exhibit multistable behavior. However, the behavior of such networks is typically studied in a single reaction volume, while cells are spatially organized into compartments that can exchange proteins. In this work, we use stochastic simulations to study the impact of compartmentalization on a two-site phosphorylation network. We characterize steady states and fluctuation-driven transitions between them as a function of the rate of protein exchange between two compartments. Surprisingly, the average time spent in a state before stochastically switching to another depends nonmonotonically on the protein exchange rate, with the most frequent switching occurring at intermediate exchange rates. At sufficiently small exchange rates, the state of the system and mean switching time are controlled largely by fluctuations in the balance of enzymes in each compartment. This leads to negatively correlated states in the compartments. For large exchange rates, the two compartments behave as a single effective compartment. However, when the compartmental volumes are unequal, the behavior differs from a single compartment with the same total volume. These results demonstrate that exchange of proteins between distinct compartments can regulate the emergent behavior of a common signaling motif.

Proteome-Wide Identification of RNA-dependent proteins and an emerging role for RNAs in Plasmodium falciparum protein complexes.

Ribonucleoprotein complexes are composed of RNA, RNA-dependent proteins (RDPs) and RNA-binding proteins (RBPs), and play fundamental roles in RNA regulation. However, in the human malaria parasite, Plasmodium falciparum, identification and characterization of these proteins are particularly limited. In this study, we use an unbiased proteome-wide approach, called R-DeeP, a method based on sucrose density gradient ultracentrifugation, to identify RDPs. Quantitative analysis by mass spectrometry identifies 898 RDPs, including 545 proteins not yet associated with RNA. Results are further validated using a combination of computational and molecular approaches. Overall, this method provides the first snapshot of the Plasmodium protein-protein interaction network in the presence and absence of RNA. R-DeeP also helps to reconstruct Plasmodium multiprotein complexes based on co-segregation and deciphers their RNA-dependence. One RDP candidate, PF3D7_0823200, is functionally characterized and validated as a true RBP. Using enhanced crosslinking and immunoprecipitation followed by high-throughput sequencing (eCLIP-seq), we demonstrate that this protein interacts with various Plasmodium non-coding transcripts, including the var genes and ap2 transcription factors.

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.

Particle Deformability Enables Control of Interactions between Membrane-Anchored Nanoparticles.

Nanoparticles adsorbed on a membrane can induce deformations of the membrane that give rise to effective interactions between the particles. Previous studies have focused primarily on rigid nanoparticles with fixed shapes. However, DNA origami technology has enabled the creation of deformable nanostructures with controllable shapes and mechanical properties, presenting new opportunities to modulate interactions between particles adsorbed on deformable surfaces. Here we use coarse-grained molecular dynamics simulations to investigate deformable, hinge-like nanostructures anchored to lipid membranes via cholesterol anchors. We characterize deformations of the particles and membrane as a function of the hinge stiffness. Flexible particles adopt open configurations to conform to a flat membrane, whereas stiffer particles induce deformations of the membrane. We further show that particles spontaneously aggregate and that cooperative effects lead to changes in their shape when they are close together. Using umbrella sampling methods, we quantify the effective interaction between two particles and show that stiffer hinge-like particles experience stronger and longer-ranged attraction. Our results demonstrate that interactions between deformable, membrane-anchored nanoparticles can be controlled by modifying mechanical properties of the particles, suggesting new ways to modulate the self-assembly of particles on deformable surfaces.

Plant Development in the Garden Pea as Revealed by Mutations in the Crd/PsYUC1 Gene.

In common with other plant species, the garden pea (Pisum sativum) produces the auxin indole-3-acetic acid (IAA) from tryptophan via a single intermediate, indole-3-pyruvic acid (IPyA). IPyA is converted to IAA by PsYUC1, also known as Crispoid (Crd). Here, we extend our understanding of the developmental processes affected by the Crd gene by examining the phenotypic effects of crd gene mutations on leaves, flowers, and roots. We show that in pea, Crd/PsYUC1 is important for the initiation and identity of leaflets and tendrils, stamens, and lateral roots. We also report on aspects of auxin deactivation in pea.

Propulsive cell entry diverts pathogens from immune degradation by remodeling the phagocytic synapse.

Phagocytosis is a critical immune function for infection control and tissue homeostasis. During phagocytosis, pathogens are internalized and degraded in phagolysosomes. For pathogens that evade immune degradation, the prevailing view is that virulence factors are required to disrupt the biogenesis of phagolysosomes. In contrast, we present here that physical forces from motile pathogens during cell entry divert them away from the canonical degradative pathway. This altered fate begins with the force-induced remodeling of the phagocytic synapse formation. We used the parasite Toxoplasma gondii as a model because live Toxoplasma actively invades host cells using gliding motility. To differentiate the effects of physical forces from virulence factors in phagocytosis, we employed magnetic forces to induce propulsive entry of inactivated Toxoplasma into macrophages. Experiments and computer simulations show that large propulsive forces hinder productive activation of receptors by preventing their spatial segregation from phosphatases at the phagocytic synapse. Consequently, the inactivated parasites are engulfed into vacuoles that fail to mature into degradative units, similar to the live motile parasite's intracellular pathway. Using yeast cells and opsonized beads, we confirmed that this mechanism is general, not specific to the parasite used. These results reveal new aspects of immune evasion by demonstrating how physical forces during active cell entry, independent of virulence factors, enable pathogens to circumvent phagolysosomal degradation.

DNA-binding protein PfAP2-P regulates parasite pathogenesis during malaria parasite blood stages.

Malaria-associated pathogenesis such as parasite invasion, egress, host cell remodelling and antigenic variation requires concerted action by many proteins, but the molecular regulation is poorly understood. Here we have characterized an essential Plasmodium-specific Apicomplexan AP2 transcription factor in Plasmodium falciparum (PfAP2-P; pathogenesis) during the blood-stage development with two peaks of expression. An inducible knockout of gene function showed that PfAP2-P is essential for trophozoite development, and critical for var gene regulation, merozoite development and parasite egress. Chromatin immunoprecipitation sequencing data collected at timepoints matching the two peaks of pfap2-p expression demonstrate PfAP2-P binding to promoters of genes controlling trophozoite development, host cell remodelling, antigenic variation and pathogenicity. Single-cell RNA sequencing and fluorescence-activated cell sorting revealed de-repression of most var genes in Δpfap2-p parasites. Δpfap2-p parasites also overexpress early gametocyte marker genes, indicating a regulatory role in sexual stage conversion. We conclude that PfAP2-P is an essential upstream transcriptional regulator at two distinct stages of the intra-erythrocytic development cycle.

A Study in Scaffold Hopping: Discovery and Optimization of Thiazolopyridines as Potent Herbicides That Inhibit Acyl-ACP Thioesterase.

In the search for new chemical entities that can control resistant weeds by addressing novel modes of action (MoAs), we were interested in further exploring a compound class that contained a 1,8-naphthyridine core. By leveraging scaffold hopping methodologies, we were able to discover the new thiazolopyridine compound class that act as potent herbicidal molecules. Further biochemical investigations allowed us to identify that the thiazolopyridines inhibit acyl-acyl carrier protein (ACP) thioesterase (FAT), with this being further confirmed via an X-ray cocrystal structure. Greenhouse trials revealed that the thiazolopyridines display excellent control of grass weed species in pre-emergence application coupled with dose response windows that enable partial selectivity in certain crops.

Plasmodium ARK2 and EB1 drive unconventional spindle dynamics, during chromosome segregation in sexual transmission stages.

The Aurora family of kinases orchestrates chromosome segregation and cytokinesis during cell division, with precise spatiotemporal regulation of its catalytic activities by distinct protein scaffolds. Plasmodium spp., the causative agents of malaria, are unicellular eukaryotes with three unique and highly divergent aurora-related kinases (ARK1-3) that are essential for asexual cellular proliferation but lack most canonical scaffolds/activators. Here we investigate the role of ARK2 during sexual proliferation of the rodent malaria Plasmodium berghei, using a combination of super-resolution microscopy, mass spectrometry, and live-cell fluorescence imaging. We find that ARK2 is primarily located at spindle microtubules in the vicinity of kinetochores during both mitosis and meiosis. Interactomic and co-localisation studies reveal several putative ARK2-associated interactors including the microtubule-interacting protein EB1, together with MISFIT and Myosin-K, but no conserved eukaryotic scaffold proteins. Gene function studies indicate that ARK2 and EB1 are complementary in driving endomitotic division and thereby parasite transmission through the mosquito. This discovery underlines the flexibility of molecular networks to rewire and drive unconventional mechanisms of chromosome segregation in the malaria parasite.

Novel insights into the role of long non-coding RNA in the human malaria parasite, Plasmodium falciparum.

The complex life cycle of Plasmodium falciparum requires coordinated gene expression regulation to allow host cell invasion, transmission, and immune evasion. Increasing evidence now suggests a major role for epigenetic mechanisms in gene expression in the parasite. In eukaryotes, many lncRNAs have been identified to be pivotal regulators of genome structure and gene expression. To investigate the regulatory roles of lncRNAs in P. falciparum we explore the intergenic lncRNA distribution in nuclear and cytoplasmic subcellular locations. Using nascent RNA expression profiles, we identify a total of 1768 lncRNAs, of which 718 (~41%) are novels in P. falciparum. The subcellular localization and stage-specific expression of several putative lncRNAs are validated using RNA-FISH. Additionally, the genome-wide occupancy of several candidate nuclear lncRNAs is explored using ChIRP. The results reveal that lncRNA occupancy sites are focal and sequence-specific with a particular enrichment for several parasite-specific gene families, including those involved in pathogenesis and sexual differentiation. Genomic and phenotypic analysis of one specific lncRNA demonstrate its importance in sexual differentiation and reproduction. Our findings bring a new level of insight into the role of lncRNAs in pathogenicity, gene regulation and sexual differentiation, opening new avenues for targeted therapeutic strategies against the deadly malaria parasite.

Small RNA sequencing of field Culex mosquitoes identifies patterns of viral infection and the mosquito immune response.

Mosquito-borne disease remains a significant burden on global health. In the United States, the major threat posed by mosquitoes is transmission of arboviruses, including West Nile virus by mosquitoes of the Culex genus. Virus metagenomic analysis of mosquito small RNA using deep sequencing and advanced bioinformatic tools enables the rapid detection of viruses and other infecting organisms, both pathogenic and non-pathogenic to humans, without any precedent knowledge. In this study, we sequenced small RNA samples from over 60 pools of Culex mosquitoes from two major areas of Southern California from 2017 to 2019 to elucidate the virome and immune responses of Culex. Our results demonstrated that small RNAs not only allowed the detection of viruses but also revealed distinct patterns of viral infection based on location, Culex species, and time. We also identified miRNAs that are most likely involved in Culex immune responses to viruses and Wolbachia bacteria, and show the utility of using small RNA to detect antiviral immune pathways including piRNAs against some pathogens. Collectively, these findings show that deep sequencing of small RNA can be used for virus discovery and surveillance. One could also conceive that such work could be accomplished in various locations across the world and over time to better understand patterns of mosquito infection and immune response to many vector-borne diseases in field samples.

Inspired by Nature: Isostere Concepts in Plant Hormone Chemistry.

Chemical concepts such as isosteres and scaffold hopping have proven to be powerful tools in agrochemical innovation processes. They offer opportunities to modify known molecular lead structures with the aim to improve a range of parameters, including biological efficacy and spectrum, physicochemical properties, stability, and toxicity. While recent biochemical insights into plant-specific receptors and signaling pathways trigger the discovery of the first lead structures, the disclosure of such a new chemical structure sparks a broad range of synthesis activities giving rise to diverse chemical innovation and often a considerable boost in biological activity. Herein, recent examples of isostere concepts in plant-hormone chemistry will be discussed, outlining how synthetic creativity can broaden the scope of natural product chemistry and giving rise to new opportunities in research fields such as abiotic stress tolerance and growth promotion.

PfAP2-MRP DNA-binding protein is a master regulator of parasite pathogenesis during malaria parasite blood stages.

Malaria pathogenicity results from the parasite's ability to invade, multiply within and then egress from the host red blood cell (RBC). Infected RBCs are remodeled, expressing antigenic variant proteins (such as PfEMP1, coded by the var gene family) for immune evasion and survival. These processes require the concerted actions of many proteins, but the molecular regulation is poorly understood. We have characterized an essential Plasmodium specific Apicomplexan AP2 (ApiAP2) transcription factor in Plasmodium falciparum (PfAP2-MRP; Master Regulator of Pathogenesis) during the intraerythrocytic developmental cycle (IDC). An inducible gene knockout approach showed that PfAP2-MRP is essential for development during the trophozoite stage, and critical for var gene regulation, merozoite development and parasite egress. ChIP-seq experiments performed at 16 hour post invasion (h.p.i.) and 40 h.p.i. matching the two peaks of PfAP2-MRP expression, demonstrate binding of PfAP2-MRP to the promoters of genes controlling trophozoite development and host cell remodeling at 16 h.p.i. and antigenic variation and pathogenicity at 40 h.p.i. Using single-cell RNA-seq and fluorescence-activated cell sorting, we show de-repression of most var genes in Δpfap2-mrp parasites that express multiple PfEMP1 proteins on the surface of infected RBCs. In addition, the Δpfap2-mrp parasites overexpress several early gametocyte marker genes at both 16 and 40 h.p.i., indicating a regulatory role in the sexual stage conversion. Using the Chromosomes Conformation Capture experiment (Hi-C), we demonstrate that deletion of PfAP2-MRP results in significant reduction of both intra-chromosomal and inter-chromosomal interactions in heterochromatin clusters. We conclude that PfAP2-MRP is a vital upstream transcriptional regulator controlling essential processes in two distinct developmental stages during the IDC that include parasite growth, chromatin structure and var gene expression.

Propulsive cell entry diverts pathogens from immune degradation by remodeling the phagocytic synapse.

Phagocytosis is a critical immune function for infection control and tissue homeostasis. This process is typically described as non-moving pathogens being internalized and degraded in phagolysosomes. For pathogens that evade immune degradation, the prevailing view is that virulence factors that biochemically disrupt the biogenesis of phagoslysosomes are required. In contrast, here we report that physical forces exerted by pathogens during cell entry divert them away from the canonical phagolysosomal degradation pathway, and this altered intracellular fate is determined at the time of phagocytic synapse formation. We used the eukaryotic parasite Toxoplasma gondii as a model because live Toxoplasma uses gliding motility to actively invade into host cells. To differentiate the effect of physical forces from that of virulence factors in phagocytosis, we developed a strategy that used magnetic forces to induce propulsive entry of inactivated Toxoplasma into macrophage cells. Experiments and computer simulations collectively reveal that large propulsive forces suppress productive activation of receptors by hindering their spatial segregation from phosphatases at the phagocytic synapse. Consequently, the inactivated parasites, instead of being degraded in phagolysosomes, are engulfed into vacuoles that fail to mature into degradative units, following an intracellular pathway strikingly similar to that of the live motile parasite. Using opsonized beads, we further confirmed that this mechanism is general, not specific to the parasite used. These results reveal previously unknown aspects of immune evasion by demonstrating how physical forces exerted during active cell entry, independent of virulence factors, can help pathogens circumvent phagolysosomal degradation.

Topographical analysis of immune cell interactions reveals a biomechanical signature for immune cytolysis.

Immune cells live intensely physical lifestyles characterized by structural plasticity, mechanosensitivity, and force exertion. Whether specific immune functions require stereotyped patterns of mechanical output, however, is largely unknown. To address this question, we used super-resolution traction force microscopy to compare cytotoxic T cell immune synapses with contacts formed by other T cell subsets and macrophages. T cell synapses were globally and locally protrusive, which was fundamentally different from the coupled pinching and pulling of macrophage phagocytosis. By spectrally decomposing the force exertion patterns of each cell type, we associated cytotoxicity with compressive strength, local protrusiveness, and the induction of complex, asymmetric interfacial topographies. These features were further validated as cytotoxic drivers by genetic disruption of cytoskeletal regulators, direct imaging of synaptic secretory events, and in silico analysis of interfacial distortion. We conclude that T cell-mediated killing and, by implication, other effector responses are supported by specialized patterns of efferent force.

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.