Building pangenome graphs.

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder, a pipeline for constructing pangenome graphs without bias or exclusion. The PanGenome Graph Builder uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events and infer phylogenetic relationships.

Comprehensive blueprint of Salmonella genomic plasticity identifies hotspots for pathogenicity genes.

Understanding the dynamic evolution of Salmonella is vital for effective bacterial infection management. This study explores the role of the flexible genome, organised in regions of genomic plasticity (RGP), in shaping the pathogenicity of Salmonella lineages. Through comprehensive genomic analysis of 12,244 Salmonella spp. genomes covering 2 species, 6 subspecies, and 46 serovars, we uncover distinct integration patterns of pathogenicity-related gene clusters into RGP, challenging traditional views of gene distribution. These RGP exhibit distinct preferences for specific genomic spots, and the presence or absence of such spots across Salmonella lineages profoundly shapes strain pathogenicity. RGP preferences are guided by conserved flanking genes surrounding integration spots, implicating their involvement in regulatory networks and functional synergies with integrated gene clusters. Additionally, we emphasise the multifaceted contributions of plasmids and prophages to the pathogenicity of diverse Salmonella lineages. Overall, this study provides a comprehensive blueprint of the pathogenicity potential of Salmonella. This unique insight identifies genomic spots in nonpathogenic lineages that hold the potential for harbouring pathogenicity genes, providing a foundation for predicting future adaptations and developing targeted strategies against emerging human pathogenic strains.

Distribution and molecular evolution of the anti-CRISPR family AcrIF7.

Anti-clustered regularly interspaced short palindromic repeats (CRISPRs) are proteins capable of blocking CRISPR-Cas systems and typically their genes are located on mobile genetic elements. Since their discovery, numerous anti-CRISPR families have been identified. However, little is known about the distribution and sequence diversity of members within a family, nor how these traits influence the anti-CRISPR's function and evolution. Here, we use AcrIF7 to explore the dissemination and molecular evolution of an anti-CRISPR family. We uncovered 5 subclusters and prevalent anti-CRISPR variants within the group. Remarkably, AcrIF7 homologs display high similarity despite their broad geographical, ecological, and temporal distribution. Although mainly associated with Pseudomonas aeruginosa, AcrIF7 was identified in distinct genetic backgrounds indicating horizontal dissemination, primarily by phages. Using mutagenesis, we recreated variation observed in databases but also extended the sequence diversity of the group. Characterisation of the variants identified residues key for the anti-CRISPR function and other contributing to its mutational tolerance. Moreover, molecular docking revealed that variants with affected function lose key interactions with its CRISPR-Cas target. Analysis of publicly available data and the generated variants suggests that the dominant AcrIF7 variant corresponds to the minimal and optimal anti-CRISPR selected in the family. Our study provides a blueprint to investigate the molecular evolution of anti-CRISPR families.

Diversity and phage sensitivity to phages of porcine enterotoxigenic Escherichia coli.

Enterotoxigenic Escherichia coli (ETEC) is a diverse and poorly characterized E. coli pathotype that causes diarrhea in humans and animals. Phages have been proposed for the veterinary biocontrol of ETEC, but effective solutions require understanding of porcine ETEC diversity that affects phage infection. Here, we sequenced and analyzed the genomes of the PHAGEBio ETEC collection, gathering 79 diverse ETEC strains isolated from European pigs with post-weaning diarrhea (PWD). We identified the virulence factors characterizing the pathotype and several antibiotic resistance genes on plasmids, while phage resistance genes and other virulence factors were mostly chromosome encoded. We experienced that ETEC strains were highly resistant to Enterobacteriaceae phage infection. It was only by enrichment of numerous diverse samples with different media and conditions, using the 41 ETEC strains of our collection as hosts, that we could isolate two lytic phages that could infect a large part of our diverse ETEC collection: vB_EcoP_ETEP21B and vB_EcoS_ETEP102. Based on genome and host range analyses, we discussed the infection strategies of the two phages and identified components of lipopolysaccharides ( LPS) as receptors for the two phages. Our detailed computational structural analysis highlights several loops and pockets in the tail fibers that may allow recognition and binding of ETEC strains, also in the presence of O-antigens. Despite the importance of receptor recognition, the diversity of the ETEC strains remains a significant challenge for isolating ETEC phages and developing sustainable phage-based products to address ETEC-induced PWD.IMPORTANCEEnterotoxigenic Escherichia coli (ETEC)-induced post-weaning diarrhea is a severe disease in piglets that leads to weight loss and potentially death, with high economic and animal welfare costs worldwide. Phage-based approaches have been proposed, but available data are insufficient to ensure efficacy. Genome analysis of an extensive collection of ETEC strains revealed that phage defense mechanisms were mostly chromosome encoded, suggesting a lower chance of spread and selection by phage exposure. The difficulty in isolating lytic phages and the molecular and structural analyses of two ETEC phages point toward a multifactorial resistance of ETEC to phage infection and the importance of extensive phage screenings specifically against clinically relevant strains. The PHAGEBio ETEC collection and these two phages are valuable tools for the scientific community to expand our knowledge on the most studied, but still enigmatic, bacterial species-E. coli.

Identification and classification of antiviral defence systems in bacteria and archaea with PADLOC reveals new system types.

To provide protection against viral infection and limit the uptake of mobile genetic elements, bacteria and archaea have evolved many diverse defence systems. The discovery and application of CRISPR-Cas adaptive immune systems has spurred recent interest in the identification and classification of new types of defence systems. Many new defence systems have recently been reported but there is a lack of accessible tools available to identify homologs of these systems in different genomes. Here, we report the Prokaryotic Antiviral Defence LOCator (PADLOC), a flexible and scalable open-source tool for defence system identification. With PADLOC, defence system genes are identified using HMM-based homologue searches, followed by validation of system completeness using gene presence/absence and synteny criteria specified by customisable system classifications. We show that PADLOC identifies defence systems with high accuracy and sensitivity. Our modular approach to organising the HMMs and system classifications allows additional defence systems to be easily integrated into the PADLOC database. To demonstrate application of PADLOC to biological questions, we used PADLOC to identify six new subtypes of known defence systems and a putative novel defence system comprised of a helicase, methylase and ATPase. PADLOC is available as a standalone package (https://github.com/padlocbio/padloc) and as a webserver (https://padloc.otago.ac.nz).

Prophages are associated with extensive CRISPR-Cas auto-immunity.

CRISPR-Cas systems require discriminating self from non-self DNA during adaptation and interference. Yet, multiple cases have been reported of bacteria containing self-targeting spacers (STS), i.e. CRISPR spacers targeting protospacers on the same genome. STS has been suggested to reflect potential auto-immunity as an unwanted side effect of CRISPR-Cas defense, or a regulatory mechanism for gene expression. Here we investigated the incidence, distribution, and evasion of STS in over 100 000 bacterial genomes. We found STS in all CRISPR-Cas types and in one fifth of all CRISPR-carrying bacteria. Notably, up to 40% of I-B and I-F CRISPR-Cas systems contained STS. We observed that STS-containing genomes almost always carry a prophage and that STS map to prophage regions in more than half of the cases. Despite carrying STS, genetic deterioration of CRISPR-Cas systems appears to be rare, suggesting a level of escape from the potentially deleterious effects of STS by other mechanisms such as anti-CRISPR proteins and CRISPR target mutations. We propose a scenario where it is common to acquire an STS against a prophage, and this may trigger more extensive STS buildup by primed spacer acquisition in type I systems, without detrimental autoimmunity effects as mechanisms of auto-immunity evasion create tolerance to STS-targeted prophages.

CRISPR-based DNA and RNA detection with liquid-liquid phase separation.

The ability to detect specific nucleic acid sequences allows for a wide range of applications such as the identification of pathogens, clinical diagnostics, and genotyping. CRISPR-Cas proteins Cas12a and Cas13a are RNA-guided endonucleases that bind and cleave specific DNA and RNA sequences, respectively. After recognition of a target sequence, both enzymes activate indiscriminate nucleic acid cleavage, which has been exploited for sequence-specific molecular diagnostics of nucleic acids. Here, we present a label-free detection approach that uses a readout based on solution turbidity caused by liquid-liquid phase separation (LLPS). Our approach relies on the fact that the LLPS of oppositely charged polymers requires polymers to be longer than a critical length. This length dependence is predicted by the Voorn-Overbeek model, which we describe in detail and validate experimentally in mixtures of polynucleotides and polycations. We show that the turbidity resulting from LLPS can be used to detect the presence of specific nucleic acid sequences by employing the programmable CRISPR-nucleases Cas12a and Cas13a. Because LLPS of polynucleotides and polycations causes solutions to become turbid, the detection of specific nucleic acid sequences can be observed with the naked eye. We furthermore demonstrate that there is an optimal polynucleotide concentration for detection. Finally, we provide a theoretical prediction that hints towards possible improvements of an LLPS-based detection assay. The deployment of LLPS complements CRISPR-based molecular diagnostic applications and facilitates easy and low-cost nucleotide sequence detection.

Development of Styrene Maleic Acid Lipid Particles as a Tool for Studies of Phage-Host Interactions.

The infection of a bacterium by a phage starts with attachment to a receptor molecule on the host cell surface by the phage. Since receptor-phage interactions are crucial to successful infections, they are major determinants of phage host range and, by extension, of the broader effects that phages have on bacterial communities. Many receptor molecules, particularly membrane proteins, are difficult to isolate because their stability is supported by their native membrane environments. Styrene maleic acid lipid particles (SMALPs), a recent advance in membrane protein studies, are the result of membrane solubilizations by styrene maleic acid (SMA) copolymer chains. SMALPs thereby allow for isolation of membrane proteins while maintaining their native environment. Here, we explore SMALPs as a tool to isolate and study phage-receptor interactions. We show that SMALPs produced from taxonomically distant bacterial membranes allow for receptor-specific decrease of viable phage counts of several model phages that span the three largest phage families. After characterizing the effects of incubation time and SMALP concentration on the activity of three distinct phages, we present evidence that the interaction between two model phages and SMALPs is specific to bacterial species and the phage receptor molecule. These interactions additionally lead to DNA ejection by nearly all particles at high phage titers. We conclude that SMALPs are a potentially highly useful tool for phage-host interaction studies.IMPORTANCE Bacteriophages (viruses that infect bacteria or phages) impact every microbial community. All phage infections start with the binding of the viral particle to a specific receptor molecule on the host cell surface. Due to its importance in phage infections, this first step is of interest to many phage-related research and applications. However, many phage receptors are difficult to isolate. Styrene maleic acid lipid particles (SMALPs) are a recently developed approach to isolate membrane proteins in their native environment. In this study, we explore SMALPs as a tool to study phage-receptor interactions. We find that different phage species bind to SMALPs, while maintaining specificity to their receptor. We then characterize the time and concentration dependence of phage-SMALP interactions and furthermore show that they lead to genome ejection by the phage. The results presented here show that SMALPs are a useful tool for future studies of phage-receptor interactions.

An educational guide for nanopore sequencing in the classroom.

The last decade has witnessed a remarkable increase in our ability to measure genetic information. Advancements of sequencing technologies are challenging the existing methods of data storage and analysis. While methods to cope with the data deluge are progressing, many biologists have lagged behind due to the fast pace of computational advancements and tools available to address their scientific questions. Future generations of biologists must be more computationally aware and capable. This means they should be trained to give them the computational skills to keep pace with technological developments. Here, we propose a model that bridges experimental and bioinformatics concepts using the Oxford Nanopore Technologies (ONT) sequencing platform. We provide both a guide to begin to empower the new generation of educators, scientists, and students in performing long-read assembly of bacterial and bacteriophage genomes and a standalone virtual machine containing all the required software and learning materials for the course.

Cas4-Cas1 fusions drive efficient PAM selection and control CRISPR adaptation.

Microbes have the unique ability to acquire immunological memories from mobile genetic invaders to protect themselves from predation. To confer CRISPR resistance, new spacers need to be compatible with a targeting requirement in the invader's DNA called the protospacer adjacent motif (PAM). Many CRISPR systems encode Cas4 proteins to ensure new spacers are integrated that meet this targeting prerequisite. Here we report that a gene fusion between cas4 and cas1 from the Geobacter sulfurreducens I-U CRISPR-Cas system is capable of introducing functional spacers carrying interference proficient TTN PAM sequences at much higher frequencies than unfused Cas4 adaptation modules. Mutations of Cas4-domain catalytic residues resulted in dramatically decreased naïve and primed spacer acquisition, and a loss of PAM selectivity showing that the Cas4 domain controls Cas1 activity. We propose the fusion gene evolved to drive the acquisition of only PAM-compatible spacers to optimize CRISPR interference.

Screening and characterization of novel specific peptides targeting MDA-MB-231 claudin-low breast carcinoma by computer-aided phage display methodologies.

BACKGROUND: Claudin-low breast carcinoma represents 19% of all breast cancer cases and is characterized by an aggressive progression with metastatic nature and high rates of relapse. Due to a lack of known specific molecular biomarkers for this breast cancer subtype, there are no targeted therapies available, which results in the worst prognosis of all breast cancer subtypes. Hence, the identification of novel biomarkers for this type of breast cancer is highly relevant for an early diagnosis. Additionally, claudin-low breast carcinoma peptide ligands can be used to design powerful drug delivery systems that specifically target this type of breast cancer. METHODS: In this work, we propose the identification of peptides for the specific recognition of MDA-MB-231, a cell line representative of claudin-low breast cancers, using phage display (both conventional panning and BRASIL). Binding assays, such as phage forming units and ELISA, were performed to select the most interesting peptides (i.e., specific to the target cells) and bioinformatics approaches were applied to putatively identify the biomarkers to which these peptides bind. RESULTS: Two peptides were selected using this methodology specifically targeting MDA-MB-231 cells, as demonstrated by a 4 to 9 log higher affinity as compared to control cells. The use of bioinformatics approaches provided relevant insights into possible cell surface targets for each peptide identified. CONCLUSIONS: The peptides herein identified may contribute to an earlier detection of claudin-low breast carcinomas and possibly to develop more individualized therapies.

Molecular aspects and comparative genomics of bacteriophage endolysins.

Phages are recognized as the most abundant and diverse entities on the planet. Their diversity is determined predominantly by their dynamic adaptation capacities when confronted with different selective pressures in an endless cycle of coevolution with a widespread group of bacterial hosts. At the end of the infection cycle, progeny virions are confronted with a rigid cell wall that hinders their release into the environment and the opportunity to start a new infection cycle. Consequently, phages encode hydrolytic enzymes, called endolysins, to digest the peptidoglycan. In this work, we bring to light all phage endolysins found in completely sequenced double-stranded nucleic acid phage genomes and uncover clues that explain the phage-endolysin-host ecology that led phages to recruit unique and specialized endolysins.

The Role of Flagellum and Flagellum-Based Motility on Salmonella Enteritidis and Escherichia coli Biofilm Formation.

Flagellum-mediated motility has been suggested to contribute to virulence by allowing bacteria to colonize and spread to new surfaces. In Salmonella enterica and Escherichia coli species, mutants affected by their flagellar motility have shown a reduced ability to form biofilms. While it is known that some species might act as co-aggregation factors for bacterial adhesion, studies of food-related biofilms have been limited to single-species biofilms and short biofilm formation periods. To assess the contribution of flagella and flagellum-based motility to adhesion and biofilm formation, two Salmonella and E. coli mutants with different flagellar phenotypes were produced: the fliC mutants, which do not produce flagella, and the motAB mutants, which are non-motile. The ability of wild-type and mutant strains to form biofilms was compared, and their relative fitness was determined in two-species biofilms with other foodborne pathogens. Our results showed a defective and significant behavior of E. coli in initial surface colonization (p < 0.05), which delayed single-species biofilm formation. Salmonella mutants were not affected by the ability to form biofilm (p > 0.05). Regarding the effect of motility/flagellum absence on bacterial fitness, none of the mutant strains seems to have their relative fitness affected in the presence of a competing species. Although the absence of motility may eventually delay initial colonization, this study suggests that motility is not essential for biofilm formation and does not have a strong impact on bacteria's fitness when a competing species is present.

Bacterial defense systems exhibit synergistic anti-phage activity.

Bacterial defense against phage predation involves diverse defense systems acting individually and concurrently, yet their interactions remain poorly understood. We investigated >100 defense systems in 42,925 bacterial genomes and identified numerous instances of their non-random co-occurrence and negative association. For several pairs of defense systems significantly co-occurring in Escherichia coli strains, we demonstrate synergistic anti-phage activity. Notably, Zorya II synergizes with Druantia III and ietAS defense systems, while tmn exhibits synergy with co-occurring systems Gabija, Septu I, and PrrC. For Gabija, tmn co-opts the sensory switch ATPase domain, enhancing anti-phage activity. Some defense system pairs that are negatively associated in E. coli show synergy and significantly co-occur in other taxa, demonstrating that bacterial immune repertoires are largely shaped by selection for resistance against host-specific phages rather than negative epistasis. Collectively, these findings demonstrate compatibility and synergy between defense systems, allowing bacteria to adopt flexible strategies for phage defense.

Accumulation of defense systems in phage-resistant strains of Pseudomonas aeruginosa.

Prokaryotes encode multiple distinct anti-phage defense systems in their genomes. However, the impact of carrying a multitude of defense systems on phage resistance remains unclear, especially in a clinical context. Using a collection of antibiotic-resistant clinical strains of Pseudomonas aeruginosa and a broad panel of phages, we demonstrate that defense systems contribute substantially to defining phage host range and that overall phage resistance scales with the number of defense systems in the bacterial genome. We show that many individual defense systems target specific phage genera and that defense systems with complementary phage specificities co-occur in P. aeruginosa genomes likely to provide benefits in phage-diverse environments. Overall, we show that phage-resistant phenotypes of P. aeruginosa with at least 19 phage defense systems exist in the populations of clinical, antibiotic-resistant P. aeruginosa strains.

Building pangenome graphs.

Pangenome graphs can represent all variation between multiple genomes, but existing methods for constructing them are biased due to reference-guided approaches. In response, we have developed PanGenome Graph Builder (PGGB), a reference-free pipeline for constructing unbi-ased pangenome graphs. PGGB uses all-to-all whole-genome alignments and learned graph embeddings to build and iteratively refine a model in which we can identify variation, measure conservation, detect recombination events, and infer phylogenetic relationships.

Phage-encoded carbohydrate-interacting proteins in the human gut.

In the human gastrointestinal tract, the gut mucosa and the bacterial component of the microbiota interact and modulate each other to accomplish a variety of critical functions. These include digestion aid, maintenance of the mucosal barrier, immune regulation, and production of vitamins, hormones, and other metabolites that are important for our health. The mucus lining of the gut is primarily composed of mucins, large glycosylated proteins with glycosylation patterns that vary depending on factors including location in the digestive tract and the local microbial population. Many gut bacteria have evolved to reside within the mucus layer and thus encode mucus-adhering and -degrading proteins. By doing so, they can influence the integrity of the mucus barrier and therefore promote either health maintenance or the onset and progression of some diseases. The viral members of the gut - mostly composed of bacteriophages - have also been shown to have mucus-interacting capabilities, but their mechanisms and effects remain largely unexplored. In this review, we discuss the role of bacteriophages in influencing mucosal integrity, indirectly via interactions with other members of the gut microbiota, or directly with the gut mucus via phage-encoded carbohydrate-interacting proteins. We additionally discuss how these phage-mucus interactions may influence health and disease states.

Structural basis for broad anti-phage immunity by DISARM.

In the evolutionary arms race against phage, bacteria have assembled a diverse arsenal of antiviral immune strategies. While the recently discovered DISARM (Defense Island System Associated with Restriction-Modification) systems can provide protection against a wide range of phage, the molecular mechanisms that underpin broad antiviral targeting but avoiding autoimmunity remain enigmatic. Here, we report cryo-EM structures of the core DISARM complex, DrmAB, both alone and in complex with an unmethylated phage DNA mimetic. These structures reveal that DrmAB core complex is autoinhibited by a trigger loop (TL) within DrmA and binding to DNA substrates containing a 5' overhang dislodges the TL, initiating a long-range structural rearrangement for DrmAB activation. Together with structure-guided in vivo studies, our work provides insights into the mechanism of phage DNA recognition and specific activation of this widespread antiviral defense system.

Gut virome profiling identifies a widespread bacteriophage family associated with metabolic syndrome.

There is significant interest in altering the course of cardiometabolic disease development via gut microbiomes. Nevertheless, the highly abundant phage members of the complex gut ecosystem -which impact gut bacteria- remain understudied. Here, we show gut virome changes associated with metabolic syndrome (MetS), a highly prevalent clinical condition preceding cardiometabolic disease, in 196 participants by combined sequencing of bulk whole genome and virus like particle communities. MetS gut viromes exhibit decreased richness and diversity. They are enriched in phages infecting Streptococcaceae and Bacteroidaceae and depleted in those infecting Bifidobacteriaceae. Differential abundance analysis identifies eighteen viral clusters (VCs) as significantly associated with either MetS or healthy viromes. Among these are a MetS-associated Roseburia VC that is related to healthy control-associated Faecalibacterium and Oscillibacter VCs. Further analysis of these VCs revealed the Candidatus Heliusviridae, a highly widespread gut phage lineage found in 90+% of participants. The identification of the temperate Ca. Heliusviridae provides a starting point to studies of phage effects on gut bacteria and the role that this plays in MetS.

Genomic characterization of four novel bacteriophages infecting the clinical pathogen Klebsiella pneumoniae.

Bacteriophages are an invaluable source of novel genetic diversity. Sequencing of phage genomes can reveal new proteins with potential uses as biotechnological and medical tools, and help unravel the diversity of biological mechanisms employed by phages to take over the host during viral infection. Aiming to expand the available collection of phage genomes, we have isolated, sequenced, and assembled the genome sequences of four phages that infect the clinical pathogen Klebsiella pneumoniae: vB_KpnP_FBKp16, vB_KpnP_FBKp27, vB_KpnM_FBKp34, and Jumbo phage vB_KpnM_FBKp24. The four phages show very low (0-13%) identity to genomic phage sequences deposited in the GenBank database. Three of the four phages encode tRNAs and have a GC content very dissimilar to that of the host. Importantly, the genome sequences of the phages reveal potentially novel DNA packaging mechanisms as well as distinct clades of tubulin spindle and nucleus shell proteins that some phages use to compartmentalize viral replication. Overall, this study contributes to uncovering previously unknown virus diversity, and provides novel candidates for phage therapy applications against antibiotic-resistant K. pneumoniae infections.