Harnessing the power of new genetic tools to illuminate Giardia biology and pathogenesis.

Giardia is a prevalent single-celled microaerophilic intestinal parasite causing diarrheal disease and significantly impacting global health. Double diploid (essentially tetraploid) Giardia trophozoites have presented a formidable challenge to the development of molecular genetic tools to interrogate gene function. High sequence divergence and the high percentage of hypothetical proteins lacking homology to proteins in other eukaryotes have limited our understanding of Giardia protein function, slowing drug target validation and development. For more than 25 years, Giardia A and B assemblages have been readily amenable to transfection with plasmids or linear DNA templates. Here, we highlight the utility and power of genetic approaches developed to assess protein function in Giardia, with particular emphasis on the more recent clustered regularly interspaced palindromic repeats/Cas9-based methods for knockdowns and knockouts. Robust and reliable molecular genetic approaches are fundamental toward the interrogation of Giardia protein function and evaluation of druggable targets. New genetic approaches tailored for the double diploid Giardia are imperative for understanding Giardia's unique biology and pathogenesis.

Characterization of a unique attachment organelle: Single-cell force spectroscopy of Giardia duodenalis trophozoites.

The unicellular parasite Giardia duodenalis is the causative agent of giardiasis, a gastrointestinal disease with global spread. In its trophozoite form, G. duodenalis can adhere to the human intestinal epithelium and a variety of other, artificial surfaces. Its attachment is facilitated by a unique microtubule-based attachment organelle, the so-called ventral disc. The mechanical function of the ventral disc, however, is still debated. Earlier studies postulated that a dynamic negative pressure under the ventral disc, generated by persistently beating flagella, mediates the attachment. Later studies suggested a suction model based on structural changes of the ventral discs, substrate clutching or grasping, or unspecific contact forces. In this study, we aim to contribute to the understanding of G. duodenalis attachment by investigating detachment characteristics and determining adhesion forces of single trophozoites on a smooth glass surface (RMS = 1.1 ± 0.2 nm) by fluidic force microscopy (FluidFM)-based single-cell force spectroscopy (SCFS). Briefly, viable adherent trophozoites were approached with a FluidFM micropipette, immobilized to the micropipette aperture by negative pressure, and detached from the surface by micropipette retraction while retract force curves were recorded. These force curves displayed novel and so far undescribed characteristics for a microorganism, namely, gradual force increase on the pulled trophozoite, with localization of adhesion force shortly before cell detachment length. Respective adhesion forces reached 7.7 ± 4.2 nN at 1 µm s-1 pulling speed. Importantly, this unique force pattern was different from that of other eukaryotic cells such as Candida albicans or oral keratinocytes, considered for comparison in this study. The latter both displayed a force pattern with force peaks of different values or force plateaus (for keratinocytes) indicative of breakage of molecular bonds of cell-anchored classes of adhesion molecules or membrane components. Furthermore, the attachment mode of G. duodenalis trophozoites was mechanically resilient to tensile forces, when the pulling speeds were raised up to 10 µm s-1 and adhesion forces increased to 28.7 ± 10.5 nN. Taken together, comparative SCSF revealed novel and unique retract force curve characteristics for attached G. duodenalis, suggesting a ligand-independent suction mechanism, that differ from those of other well described eukaryotes.

ß-Variational autoencoders and transformers for reduced-order modelling of fluid flows.

Variational autoencoder architectures have the potential to develop reduced-order models for chaotic fluid flows. We propose a method for learning compact and near-orthogonal reduced-order models using a combination of a ß-variational autoencoder and a transformer, tested on numerical data from a two-dimensional viscous flow in both periodic and chaotic regimes. The ß-variational autoencoder is trained to learn a compact latent representation of the flow velocity, and the transformer is trained to predict the temporal dynamics in latent-space. Using the ß-variational autoencoder to learn disentangled representations in latent-space, we obtain a more interpretable flow model with features that resemble those observed in the proper orthogonal decomposition, but with a more efficient representation. Using Poincaré maps, the results show that our method can capture the underlying dynamics of the flow outperforming other prediction models. The proposed method has potential applications in other fields such as weather forecasting, structural dynamics or biomedical engineering.

The domed architecture of Giardias ventral disc is necessary for attachment and host pathogenesis.

After ingestion of dormant cysts, the widespread protozoan parasite Giardia lamblia colonizes the host gastrointestinal tract via direct and reversible attachment using a novel microtubule organelle, the ventral disc. Extracellular attachment to the host allows the parasite to resist peristaltic flow, facilitates colonization and is proposed to cause damage to the microvilli of host enterocytes as well as disrupt host barrier integrity. The 9 um in diameter ventral disc is defined by a highly complex architecture of unique protein complexes scaffolded onto a spiral microtubule (MT) array of one hundred parallel, uniformly spaced MT polymers that bend approximately one and a quarter turns to form a domed structure. To investigate the role of disc-mediated attachment in causing epithelial cell damage, we used a new approach to rapidly create a stable quadruple knockout of Giardia of an essential ventral disc protein, MBP, using a new method of CRISPR-mediated gene disruption with multiple positive selectable markers. MBP quadruple KO mutant discs lack the characteristic domed architecture and possess a flattened crescent or horseshoe-shaped conformation that lacks the overlapping region, with severe defects in the microribbon-crossbridge (MR-CB) complex structure. MBP KO mutants are also unable to resist fluid flow required for attachment to inert surfaces. Importantly, MBP KO mutants have 100% penetrance off positive selection, which is essential for quantification of in vivo impacts of disc and attachment mutants with host cells. Using a new gastrointestinal organoid model of pathogenesis, we found that MBP KO infections had a significantly reduced ability to cause the barrier breakdown characteristic of wild-type infections. Overall, this work provides direct evidence of the role of MBP in creating the domed disc, as well as the first direct evidence that parasite attachment is necessary for host pathology, specifically epithelial barrier breakdown.

Highly contiguous genomes of human clinical isolates of Giardia duodenalis reveal assemblage- and sub-assemblage-specific presence-absence variation in protein-coding genes.

Giardia duodenalis (syn. G. intestinalis, G. lamblia) is a widespread gastrointestinal protozoan parasite with debated taxonomic status. Currently, eight distinct genetic sub-groups, termed assemblages A-H, are defined based on a few genetic markers. Assemblages A and B may represent distinct species and are both of human public health relevance. Genomic studies are scarce and the few reference genomes available, in particular for assemblage B, are insufficient for adequate comparative genomics. Here, by combining long- and short-read sequences generated by PacBio and Illumina sequencing technologies, we provide nine annotated genome sequences for reference from new clinical isolates (four assemblage A and five assemblage B parasite isolates). Isolates chosen represent the currently accepted classification of sub-assemblages AI, AII, BIII and BIV. Synteny over the whole genome was generally high, but we report chromosome-level translocations as a feature that distinguishes assemblage A from B parasites. Orthologue gene group analysis was used to define gene content differences between assemblage A and B and to contribute a gene-set-based operational definition of respective taxonomic units. Giardia is tetraploid, and high allelic sequence heterogeneity (ASH) for assemblage B vs. assemblage A has been observed so far. Noteworthy, here we report an extremely low ASH (0.002%) for one of the assemblage B isolates (a value even lower than the reference assemblage A isolate WB-C6). This challenges the view of low ASH being a notable feature that distinguishes assemblage A from B parasites, and low ASH allowed assembly of the most contiguous assemblage B genome currently available for reference. In conclusion, the description of nine highly contiguous genome assemblies of new isolates of G. duodenalis assemblage A and B adds to our understanding of the genomics and species population structure of this widespread zoonotic parasite.

Nucleophilic Thiol Proteins Bind Covalently to Abasic Sites in DNA.

In the course of studies on the enhancement of 1,2-dibromoethane-induced DNA base pair mutations by O6-alkylguanine-DNA alkyltransferase (AGT, MGMT), we discovered the facile reaction of AGT with an abasic site in DNA, leading to covalent cross-linking. The binding of AGT differs from the mechanism reported for the protein HMCES; instead it appears to involve formation of a stable thioglycoside. Facile cross-linking was also observed with the protease papain, which like AGT has a low pKa cysteine, and the tripeptide glutathione.

Efficient CRISPR/Cas9-mediated gene disruption in the tetraploid protist Giardia intestinalis.

CRISPR/Cas9-mediated genome editing has become an extremely powerful technique used to modify gene expression in many organisms, including parasitic protists. Giardia intestinalis, a protist parasite that infects approximately 280 million people around the world each year, has been eluding the use of CRISPR/Cas9 to generate knockout cell lines due to its tetraploid genome. In this work, we show the ability of the in vitro assembled CRISPR/Cas9 components to successfully edit the genome of G. intestinalis. The cell line that stably expresses Cas9 in both nuclei of G. intestinalis showed effective recombination of the cassette containing the transcription units for the gRNA and the resistance marker. This highly efficient process led to the removal of all gene copies at once for three independent experimental genes, mem, cwp1 and mlf1. The method was also applicable to incomplete disruption of the essential gene, as evidenced by significantly reduced expression of tom40. Finally, testing the efficiency of Cas9-induced recombination revealed that homologous arms as short as 150 bp can be sufficient to establish a complete knockout cell line in G. intestinalis.

Affinity-purified Plasmodium tubulin provides a key reagent for antimalarial drug development.

Hirst et al. used a TOG-domain-based affinity-purification approach to reconstitute and define the in vitro dynamics of blood-stage Plasmodium falciparum αß-tubulin. This provides a key reagent for defining parasite microtubule (MT) dynamics and for evaluating the efficacy of anti-MT drugs throughout the complex parasite life cycle.

Integrating multi-fidelity blood flow data with reduced-order data assimilation.

High-fidelity patient-specific modeling of cardiovascular flows and hemodynamics is challenging. Direct blood flow measurement inside the body with in-vivo measurement modalities such as 4D flow magnetic resonance imaging (4D flow MRI) suffer from low resolution and acquisition noise. In-vitro experimental modeling and patient-specific computational fluid dynamics (CFD) models are subject to uncertainty in patient-specific boundary conditions and model parameters. Furthermore, collecting blood flow data in the near-wall region (e.g., wall shear stress) with experimental measurement modalities poses additional challenges. In this study, a computationally efficient data assimilation method called reduced-order modeling Kalman filter (ROM-KF) was proposed, which combined a sequential Kalman filter with reduced-order modeling using a linear model provided by dynamic mode decomposition (DMD). The goal of ROM-KF was to overcome low resolution and noise in experimental and uncertainty in CFD modeling of cardiovascular flows. The accuracy of the method was assessed with 1D Womersley flow, 2D idealized aneurysm, and 3D patient-specific cerebral aneurysm models. Synthetic experimental data were used to enable direct quantification of errors using benchmark datasets. The accuracy of ROM-KF in reconstructing near-wall hemodynamics was assessed by applying the method to problems where near-wall blood flow data were missing in the experimental dataset. The ROM-KF method provided blood flow data that were more accurate than the computational and synthetic experimental datasets and improved near-wall hemodynamics quantification.

Data-driven cardiovascular flow modelling: examples and opportunities.

High-fidelity blood flow modelling is crucial for enhancing our understanding of cardiovascular disease. Despite significant advances in computational and experimental characterization of blood flow, the knowledge that we can acquire from such investigations remains limited by the presence of uncertainty in parameters, low resolution, and measurement noise. Additionally, extracting useful information from these datasets is challenging. Data-driven modelling techniques have the potential to overcome these challenges and transform cardiovascular flow modelling. Here, we review several data-driven modelling techniques, highlight the common ideas and principles that emerge across numerous such techniques, and provide illustrative examples of how they could be used in the context of cardiovascular fluid mechanics. In particular, we discuss principal component analysis (PCA), robust PCA, compressed sensing, the Kalman filter for data assimilation, low-rank data recovery, and several additional methods for reduced-order modelling of cardiovascular flows, including the dynamic mode decomposition and the sparse identification of nonlinear dynamics. All techniques are presented in the context of cardiovascular flows with simple examples. These data-driven modelling techniques have the potential to transform computational and experimental cardiovascular research, and we discuss challenges and opportunities in applying these techniques in the field, looking ultimately towards data-driven patient-specific blood flow modelling.

Draft Genome Sequence of the Free-Living, Iridescent Bacterium Tenacibaculum mesophilum Strain ECR.

Here, we report the genome sequence of Tenacibaculum mesophilum strain ECR, which was isolated from the river/ocean interface at Trunk River in Falmouth, Massachusetts. The isolation and sequencing were performed as part of the 2016 and 2018 Microbial Diversity courses at the Marine Biological Laboratory in Woods Hole, Massachusetts.

The Tail of Kinesin-14a in Giardia Is a Dual Regulator of Motility.

Kinesin-14s are microtubule-based motor proteins that play important roles in mitotic spindle assembly [1]. Ncd-type kinesin-14s are a subset of kinesin-14 motors that exist as homodimers with an N-terminal microtubule-binding tail, a coiled-coil central stalk (central stalk), a neck, and two identical C-terminal motor domains. To date, no Ncd-type kinesin-14 has been found to naturally exhibit long-distance minus-end-directed processive motility on single microtubules as individual homodimers. Here, we show that GiKIN14a from Giardia intestinalis [2] is an unconventional Ncd-type kinesin-14 that uses its N-terminal microtubule-binding tail to achieve minus-end-directed processivity on single microtubules over micrometer distances as a homodimer. We further find that although truncation of the N-terminal tail greatly reduces GiKIN14a processivity, the resulting tailless construct GiKIN14a-Δtail is still a minimally processive motor and moves its center of mass via discrete 8-nm steps on the microtubule. In addition, full-length GiKIN14a has significantly higher stepping and ATP hydrolysis rates than does GiKIN14a-Δtail. Inserting a flexible polypeptide linker into the central stalk of full-length GiKIN14a nearly reduces its ATP hydrolysis rate to that of GiKIN14a-Δtail. Collectively, our results reveal that the N-terminal tail of GiKIN14a is a de facto dual regulator of motility and reinforce the notion of the central stalk as a key mechanical determinant of kinesin-14 motility [3].

Disc-associated proteins mediate the unusual hyperstability of the ventral disc in Giardia lamblia.

Giardia lamblia, a widespread parasitic protozoan, attaches to the host gastrointestinal epithelium by using the ventral disc, a complex microtubule (MT) organelle. The 'cup-like' disc is formed by a spiral MT array that scaffolds numerous disc-associated proteins (DAPs) and higher-order protein complexes. In interphase, the disc is hyperstable and has limited MT dynamics; however, it remains unclear how DAPs confer these properties. To investigate mechanisms of hyperstability, we confirmed the disc-specific localization of over 50 new DAPs identified by using both a disc proteome and an ongoing GFP localization screen. DAPs localize to specific disc regions and many lack similarity to known proteins. By screening 14 CRISPRi-mediated DAP knockdown (KD) strains for defects in hyperstability and MT dynamics, we identified two strains - DAP5188KD and DAP6751KD -with discs that dissociate following high-salt fractionation. Discs in the DAP5188KD strain were also sensitive to treatment with the MT-polymerization inhibitor nocodazole. Thus, we confirm here that at least two of the 87 known DAPs confer hyperstable properties to the disc MTs, and we anticipate that other DAPs contribute to disc MT stability, nucleation and assembly.

Microtubule organelles in Giardia.

Giardia lamblia is a widespread parasitic protist with a complex MT cytoskeleton that is critical for motility, attachment, mitosis and cell division, and transitions between its two life cycle stages-the infectious cyst and flagellated trophozoite. Giardia trophozoites have both highly dynamic and highly stable MT organelles, including the ventral disc, eight flagella, the median body and the funis. The ventral disc, an elaborate MT organelle, is essential for the parasite's attachment to the intestinal villi to avoid peristalsis. Giardia's four flagellar pairs enable swimming motility and may also promote attachment. They are maintained at different equilibrium lengths and are distinguished by their long cytoplasmic regions and novel extra-axonemal structures. The functions of the median body and funis, MT organelles unique to Giardia, remain less understood. In addition to conserved MT-associated proteins, the genome is enriched in ankyrins, NEKs, and novel hypothetical proteins that also associate with the MT cytoskeleton. High-resolution ultrastructural imaging and a current inventory of more than 300 proteins associated with Giardia's MT cytoskeleton lay the groundwork for future mechanistic analyses of parasite attachment to the host, motility, cell division, and encystation/excystation. Giardia's unique MT organelles exemplify the capacity of MT polymers to generate intricate structures that are diverse in both form and function. Thus, beyond its relevance to pathogenesis, the study of Giardia's MT cytoskeleton informs basic cytoskeletal biology and cellular evolution. With the availability of new molecular genetic tools to disrupt gene function, we anticipate a new era of cytoskeletal discovery in Giardia.

Recent advances in functional research in Giardia intestinalis.

This review considers current advances in tools to investigate the functional biology of Giardia, it's coding and non-coding genes, features and cellular and molecular biology. We consider major gaps in current knowledge of the parasite and discuss the present state-of-the-art in its in vivo and in vitro cultivation. Advances in in silico tools, including for the modelling non-coding RNAs and genomic elements, as well as detailed exploration of coding genes through inferred homology to model organisms, have provided significant, primary level insight. Improved methods to model the three-dimensional structure of proteins offer new insights into their function, and binding interactions with ligands, other proteins or precursor drugs, and offer substantial opportunities to prioritise proteins for further study and experimentation. These approaches can be supplemented by the growing and highly accessible arsenal of systems-based methods now being applied to Giardia, led by genomic, transcriptomic and proteomic methods, but rapidly incorporating advanced tools for detection of real-time transcription, evaluation of chromatin states and direct measurement of macromolecular complexes. Methods to directly interrogate and perturb gene function have made major leaps in recent years, with CRISPr-interference now available. These approaches, coupled with protein over-expression, fluorescent labelling and in vitro and in vivo imaging, are set to revolutionize the field and herald an exciting time during which the field may finally realise Giardia's long proposed potential as a model parasite and eukaryote.

Microbial community dynamics and coexistence in a sulfide-driven phototrophic bloom.

BACKGROUND: Lagoons are common along coastlines worldwide and are important for biogeochemical element cycling, coastal biodiversity, coastal erosion protection and blue carbon sequestration. These ecosystems are frequently disturbed by weather, tides, and human activities. Here, we investigated a shallow lagoon in New England. The brackish ecosystem releases hydrogen sulfide particularly upon physical disturbance, causing blooms of anoxygenic sulfur-oxidizing phototrophs. To study the habitat, microbial community structure, assembly and function we carried out in situ experiments investigating the bloom dynamics over time. RESULTS: Phototrophic microbial mats and permanently or seasonally stratified water columns commonly contain multiple phototrophic lineages that coexist based on their light, oxygen and nutrient preferences. We describe similar coexistence patterns and ecological niches in estuarine planktonic blooms of phototrophs. The water column showed steep gradients of oxygen, pH, sulfate, sulfide, and salinity. The upper part of the bloom was dominated by aerobic phototrophic Cyanobacteria, the middle and lower parts by anoxygenic purple sulfur bacteria (Chromatiales) and green sulfur bacteria (Chlorobiales), respectively. We show stable coexistence of phototrophic lineages from five bacterial phyla and present metagenome-assembled genomes (MAGs) of two uncultured Chlorobaculum and Prosthecochloris species. In addition to genes involved in sulfur oxidation and photopigment biosynthesis the MAGs contained complete operons encoding for terminal oxidases. The metagenomes also contained numerous contigs affiliating with Microviridae viruses, potentially affecting Chlorobi. Our data suggest a short sulfur cycle within the bloom in which elemental sulfur produced by sulfide-oxidizing phototrophs is most likely reduced back to sulfide by Desulfuromonas sp. CONCLUSIONS: The release of sulfide creates a habitat selecting for anoxygenic sulfur-oxidizing phototrophs, which in turn create a niche for sulfur reducers. Strong syntrophism between these guilds apparently drives a short sulfur cycle that may explain the rapid development of the bloom. The fast growth and high biomass yield of Chlorobi-affiliated organisms implies that the studied lineages of green sulfur bacteria can thrive in hypoxic habitats. This oxygen tolerance is corroborated by oxidases found in MAGs of uncultured Chlorobi. The findings improve our understanding of the ecology and ecophysiology of anoxygenic phototrophs and their impact on the coupled biogeochemical cycles of sulfur and carbon.

Length-dependent disassembly maintains four different flagellar lengths in Giardia.

With eight flagella of four different lengths, the parasitic protist Giardia is an ideal model to evaluate flagellar assembly and length regulation. To determine how four different flagellar lengths are maintained, we used live-cell quantitative imaging and mathematical modeling of conserved components of intraflagellar transport (IFT)-mediated assembly and kinesin-13-mediated disassembly in different flagellar pairs. Each axoneme has a long cytoplasmic region extending from the basal body, and transitions to a canonical membrane-bound flagellum at the 'flagellar pore'. We determined that each flagellar pore is the site of IFT accumulation and injection, defining a diffusion barrier functionally analogous to the transition zone. IFT-mediated assembly is length-independent, as train size, speed, and injection frequencies are similar for all flagella. We demonstrate that kinesin-13 localization to the flagellar tips is inversely correlated to flagellar length. Therefore, we propose a model where a length-dependent disassembly mechanism controls multiple flagellar lengths within the same cell.

Draft Genome Sequence of the Caffeine-Degrading Methylotroph Methylorubrum populi Pinkel.

A pink-pigmented facultative methylotroph, Methylorubrum populi Pinkel, was isolated from compost by selective enrichment with caffeine (3,5,7-trimethylxanthine) as the sole carbon, nitrogen, and energy source. We report here its high-quality draft genome sequence, assembled in 35 contigs totaling 5,630,907 bp. We identified 5,681 protein-coding sequences, including those putatively involved in caffeine degradation.

Draft Genome Sequence of the Iridescent Marine Bacterium Tenacibaculum discolor Strain IMLK18.

We report here the draft genome sequence of a strain of Tenacibaculum discolor (Bacteroidetes) that was isolated from the river-ocean interface at Trunk River in Falmouth, Massachusetts. The isolation and genomic sequencing were performed during the 2016 and 2018 Microbial Diversity summer programs at the Marine Biological Laboratory in Woods Hole, Massachusetts.