Building pyramids against the evolutionary emergence of pathogens.

Mutations allowing pathogens to escape host immunity promote the spread of infectious diseases in heterogeneous host populations and can lead to major epidemics. Understanding the conditions that slow down this evolution is key for the development of durable control strategies against pathogens. Here, we use theory and experiments to compare the efficacy of three strategies for the deployment of resistance: (i) a mixing strategy where the host population contains two single-resistant genotypes, (ii) a pyramiding strategy where the host carries a double-resistant genotype, (iii) a combining strategy where the host population is a mix of a single-resistant genotype and a double-resistant genotype. First, we use evolutionary epidemiology theory to clarify the interplay between demographic stochasticity and evolutionary dynamics to show that the pyramiding strategy always yields lower probability of evolutionary emergence. Second, we test experimentally these predictions with the introduction of bacteriophages into bacterial populations where we manipulated the diversity and the depth of immunity using a Clustered Regularly Interspaced Short Palindromic Repeats-CRISPR associated (CRISPR-Cas) system. These biological assays confirm that pyramiding multiple defences into the same host genotype and avoiding combination with single-defence genotypes is a robust way to reduce pathogen evolutionary emergence. The experimental validation of these theoretical recommendations has practical implications in various areas, including for the optimal deployment of resistance varieties in agriculture and for the design of durable vaccination strategies.

Longitudinal Study of Lactococcus Phages in a Canadian Cheese Factory.

The presence of virulent phages is closely monitored during cheese manufacturing, as these bacterial viruses can significantly slow down the milk fermentation process and lead to low-quality cheeses. From 2001 to 2020, whey samples from cheddar cheese production in a Canadian factory were monitored for the presence of virulent phages capable of infecting proprietary strains of Lactococcus cremoris and Lactococcus lactis used in starter cultures. Phages were successfully isolated from 932 whey samples using standard plaque assays and several industrial Lactococcus strains as hosts. A multiplex PCR assay assigned 97% of these phage isolates to the Skunavirus genus, 2% to the P335 group, and 1% to the Ceduovirus genus. DNA restriction profiles and a multilocus sequence typing (MLST) scheme distinguished at least 241 unique lactococcal phages from these isolates. While most phages were isolated only once, 93 of them (out of 241, 39%) were isolated multiple times. Phage GL7 was isolated 132 times from 2006 to 2020, demonstrating that phages can persist in a cheese factory for long periods of time. Phylogenetic analysis of MLST sequences showed that phages could be clustered based on their bacterial hosts rather than their year of isolation. Host range analysis showed that Skunavirus phages exhibited a very narrow host range, whereas some Ceduovirus and P335 phages had a broader host range. Overall, the host range information was useful in improving the starter culture rotation by identifying phage-unrelated strains and helped mitigating the risk of fermentation failure due to virulent phages. IMPORTANCE Although lactococcal phages have been observed in cheese production settings for almost a century, few longitudinal studies have been performed. This 20-year study describes the close monitoring of dairy lactococcal phages in a cheddar cheese factory. Routine monitoring was conducted by factory staff, and when whey samples were found to inhibit industrial starter cultures under laboratory conditions, they were sent to an academic research laboratory for phage isolation and characterization. This led to a collection of at least 241 unique lactococcal phages, which were characterized through PCR typing and MLST profiling. Phages of the Skunavirus genus were by far the most dominant. Most phages lysed a small subset of the Lactococcus strains. These findings guided the industrial partner in adapting the starter culture schedule by using phage-unrelated strains in starter cultures and removing some strains from the starter rotation. This phage control strategy could be adapted for other large-scale bacterial fermentation processes.

DNA tandem repeats contribute to the genetic diversity of Brevibacterium aurantiacum phages.

This report presents the characterization of the first virulent phages infecting Brevibacterium aurantiacum, a bacterial species used during the manufacture of surface-ripened cheeses. These phages were also responsible for flavour and colour defects in surface-ripened cheeses. Sixteen phages (out of 62 isolates) were selected for genome sequencing and comparative analyses. These cos-type phages with a long non-contractile tail currently belong to the Siphoviridae family (Caudovirales order). Their genome sizes vary from 35,637 to 36,825 bp and, similar to their host, have a high GC content (~61%). Genes encoding for an immunity repressor, an excisionase and a truncated integrase were found, suggesting that these virulent phages may be derived from a prophage. Their genomic organization is highly conserved, with most of the diversity coming from the presence of long (198 bp) DNA tandem repeats (TRs) within an open reading frame coding for a protein of unknown function. We categorized these phages into seven genomic groups according to their number of TR, which ranged from two to eight. Moreover, we showed that TRs are widespread in phage genomes, found in more than 85% of the genomes available in public databases.

Novel Genus of Phages Infecting Streptococcus thermophilus: Genomic and Morphological Characterization.

Streptococcus thermophilus is a lactic acid bacterium commonly used for the manufacture of yogurt and specialty cheeses. Virulent phages represent a major risk for milk fermentation processes worldwide, as they can inactivate the added starter bacterial cells, leading to low-quality fermented dairy products. To date, four genetically distinct groups of phages infecting S. thermophilus have been described. Here, we describe a fifth group. Phages P738 and D4446 are virulent siphophages that infect a few industrial strains of S. thermophilus The genomes of phages P738 and D4446 were sequenced and found to contain 34,037 and 33,656 bp as well as 48 and 46 open reading frames, respectively. Comparative genomic analyses revealed that the two phages are closely related to each other but display very limited similarities to other S. thermophilus phages. In fact, these two novel S. thermophilus phages share similarities with streptococcal phages of nondairy origin, suggesting that they emerged recently in the dairy environment.IMPORTANCE Despite decades of research and adapted antiphage strategies such as CRISPR-Cas systems, virulent phages are still a persistent risk for the milk fermentation industry worldwide, as they can cause manufacturing failures and alter product quality. Phages P738 and D4446 are novel virulent phages that infect the food-grade Gram-positive bacterial species Streptococcus thermophilus These two related viruses represent a fifth group of S. thermophilus phages, as they are significantly distinct from other known S. thermophilus phages. Both phages share similarities with phages infecting nondairy streptococci, suggesting their recent emergence and probable coexistence in dairy environments. These findings highlight the necessity of phage surveillance programs as the phage population evolves in response to the application of antiphage strategies.

Lactococcus lactis type III-A CRISPR-Cas system cleaves bacteriophage RNA.

CRISPR-Cas defends microbial cells against invading nucleic acids including viral genomes. Recent studies have shown that type III-A CRISPR-Cas systems target both RNA and DNA in a transcription-dependent manner. We previously found a type III-A system on a conjugative plasmid in Lactococcus lactis which provided resistance against virulent phages of the Siphoviridae family. Its naturally occurring spacers are oriented to generate crRNAs complementary to target phage mRNA, suggesting transcription-dependent targeting. Here, we show that only constructs whose spacers produce crRNAs complementary to the phage mRNA confer phage resistance in L. lactis. In vivo nucleic acid cleavage assays showed that cleavage of phage dsDNA genome was not detected within phage-infected L. lactis cells. On the other hand, Northern blots indicated that the lactococcal CRISPR-Cas cleaves phage mRNA in vivo. These results cannot exclude that single-stranded phage DNA is not being targeted, but phage DNA replication has been shown to be impaired.

Comparison of advanced whole genome sequence-based methods to distinguish strains of Salmonella enterica serovar Heidelberg involved in foodborne outbreaks in Québec.

Salmonella enterica serovar Heidelberg (S. Heidelberg) is one of the top serovars causing human salmonellosis. This serovar ranks second and third among serovars that cause human infections in Québec and Canada, respectively, and has been associated with severe infections. Traditional typing methods such as PFGE do not display adequate discrimination required to resolve outbreak investigations due to the low level of genetic diversity of isolates belonging to this serovar. This study evaluates the ability of four whole genome sequence (WGS)-based typing methods to differentiate among 145 S. Heidelberg strains involved in four distinct outbreak events and sporadic cases of salmonellosis that occurred in Québec between 2007 and 2016. Isolates from all outbreaks were indistinguishable by PFGE. The core genome single nucleotide variant (SNV), core genome multilocus sequence typing (MLST) and whole genome MLST approaches were highly discriminatory and separated outbreak strains into four distinct phylogenetic clusters that were concordant with the epidemiological data. The clustered regularly interspaced short palindromic repeats (CRISPR) typing method was less discriminatory. However, CRISPR typing may be used as a secondary method to differentiate isolates of S. Heidelberg that are genetically similar but epidemiologically unrelated to outbreak events. WGS-based typing methods provide a highly discriminatory alternative to PFGE for the laboratory investigation of foodborne outbreaks.

Genomic Diversity of Phages Infecting Probiotic Strains of Lactobacillus paracasei.

Strains of the Lactobacillus casei group have been extensively studied because some are used as probiotics in foods. Conversely, their phages have received much less attention. We analyzed the complete genome sequences of five L. paracasei temperate phages: CL1, CL2, iLp84, iLp1308, and iA2. Only phage iA2 could not replicate in an indicator strain. The genome lengths ranged from 34,155 bp (iA2) to 39,474 bp (CL1). Phages iA2 and iLp1308 (34,176 bp) possess the smallest genomes reported, thus far, for phages of the L. casei group. The GC contents of the five phage genomes ranged from 44.8 to 45.6%. As observed with many other phages, their genomes were organized as follows: genes coding for DNA packaging, morphogenesis, lysis, lysogeny, and replication. Phages CL1, CL2, and iLp1308 are highly related to each other. Phage iLp84 was also related to these three phages, but the similarities were limited to gene products involved in DNA packaging and structural proteins. Genomic fragments of phages CL1, CL2, iLp1308, and iLp84 were found in several genomes of L. casei strains. Prophage iA2 is unrelated to these four phages, but almost all of its genome was found in at least four L. casei strains. Overall, these phages are distinct from previously characterized Lactobacillus phages. Our results highlight the diversity of L. casei phages and indicate frequent DNA exchanges between phages and their hosts.

A virulent phage infecting Lactococcus garvieae, with homology to Lactococcus lactis phages.

A new virulent phage belonging to the Siphoviridae family and able to infect Lactococcus garvieae strains was isolated from compost soil. Phage GE1 has a prolate capsid (56 by 38 nm) and a long noncontractile tail (123 nm). It had a burst size of 139 and a latent period of 31 min. Its host range was limited to only two L. garvieae strains out of 73 tested. Phage GE1 has a double-stranded DNA genome of 24,847 bp containing 48 predicted open reading frames (ORFs). Putative functions could be assigned to only 14 ORFs, and significant matches in public databases were found for only 17 ORFs, indicating that GE1 is a novel phage and its genome contains several new viral genes and encodes several new viral proteins. Of these 17 ORFs, 16 were homologous to deduced proteins of virulent phages infecting the dairy bacterium Lactococcus lactis, including previously characterized prolate-headed phages. Comparative genome analysis confirmed the relatedness of L. garvieae phage GE1 to L. lactis phages c2 (22,172 bp) and Q54 (26,537 bp), although its genome organization was closer to that of phage c2. Phage GE1 did not infect any of the 58 L. lactis strains tested. This study suggests that phages infecting different lactococcal species may have a common ancestor.

A new Microviridae phage isolated from a failed biotechnological process driven by Escherichia coli.

Bacteriophages are present in every environment that supports bacterial growth, including man made ecological niches. Virulent phages may even slow or, in more severe cases, interrupt bioprocesses driven by bacteria. Escherichia coli is one of the most widely used bacteria for large-scale bioprocesses; however, literature describing phage-host interactions in this industrial context is sparse. Here, we describe phage MED1 isolated from a failed industrial process. Phage MED1 (Microviridae family, with a single-stranded DNA [ssDNA] genome) is highly similar to the archetypal phage phiX174, sharing >95% identity between their genomic sequences. Whole-genome phylogenetic analysis of 52 microvirus genomes from public databases revealed three genotypes (alpha3, G4, and phiX174). Phage MED1 belongs to the phiX174 group. We analyzed the distribution of single nucleotide variants in MED1 and 18 other phiX174-like genomes and found that there are more missense mutations in genes G, B, and E than in the other genes of these genomes. Gene G encodes the spike protein, involved in host attachment. The evolution of this protein likely results from the selective pressure on phages to rapidly adapt to the molecular diversity found at the surface of their hosts.

Complete genome of Escherichia phages Carena and JoYop.

Escherichia phages Carena and JoYop were isolated from water samples in Abidjan (Cote d'Ivoire). Their genomes comprise 39,283 and 169,193 bp, encoding 44 and 246 predicted genes, respectively. Carena shares 93.4% nucleotide identity with Escherichia podophage CarlSpitteler (Berlinvirus), and JoYop shows 95.6% identity with Escherichia myophage ADUt (Tequatrovirus).

Identification of a new P335 subgroup through molecular analysis of lactococcal phages Q33 and BM13.

Lactococcal dairy starter strains are under constant threat from phages in dairy fermentation facilities, especially by members of the so-called 936, P335, and c2 species. Among these three phage groups, members of the P335 species are the most genetically diverse. Here, we present the complete genome sequences of two P335-type phages, Q33 and BM13, isolated in North America and representing a novel lineage within this phage group. The Q33 and BM13 genomes exhibit homology, not only to P335-type, but also to elements of the 936-type phage sequences. The two phage genomes also have close relatedness to phages infecting Enterococcus and Clostridium, a heretofore unknown feature among lactococcal P335 phages. The Q33 and BM13 genomes are organized in functionally related clusters with genes encoding functions such as DNA replication and packaging, morphogenesis, and host cell lysis. Electron micrographic analysis of the two phages highlights the presence of a baseplate more reminiscent of the baseplate of 936 phages than that of the majority of members of the P335 group, with the exception of r1t and LC3.

Involvement of the major capsid protein and two early-expressed phage genes in the activity of the lactococcal abortive infection mechanism AbiT.

The dairy industry uses the mesophilic, Gram-positive, lactic acid bacterium (LAB) Lactococcus lactis to produce an array of fermented milk products. Milk fermentation processes are susceptible to contamination by virulent phages, but a plethora of phage control strategies are available. One of the most efficient is to use LAB strains carrying phage resistance systems such as abortive infection (Abi) mechanisms. Yet, the mode of action of most Abi systems remains poorly documented. Here, we shed further light on the antiviral activity of the lactococcal AbiT system. Twenty-eight AbiT-resistant phage mutants derived from the wild-type AbiT-sensitive lactococcal phages p2, bIL170, and P008 were isolated and characterized. Comparative genomic analyses identified three different genes that were mutated in these virulent AbiT-insensitive phage derivatives: e14 (bIL170 [e14(bIL170)]), orf41 (P008 [orf41(P008)]), and orf6 (p2 [orf6(p2)] and P008 [orf6(P008)]). The genes e14(bIL170) and orf41(P008) are part of the early-expressed genomic region, but bioinformatic analyses did not identify their putative function. orf6 is found in the phage morphogenesis module. Antibodies were raised against purified recombinant ORF6, and immunoelectron microscopy revealed that it is the major capsid protein (MCP). Coexpression in L. lactis of ORF6(p2) and ORF5(p2), a protease, led to the formation of procapsids. To our knowledge, AbiT is the first Abi system involving distinct phage genes.

A truncated anti-CRISPR protein prevents spacer acquisition but not interference.

CRISPR-Cas systems in prokaryotic cells provide an adaptive immunity against invading nucleic acids. For example, phage infection leads to addition of new immunity (spacer acquisition) and DNA cleavage (interference) in the bacterial model species Streptococcus thermophilus, which primarily relies on Cas9-containing CRISPR-Cas systems. Phages can counteract this defense system through mutations in the targeted protospacers or by encoding anti-CRISPR proteins (ACRs) that block Cas9 interference activity. Here, we show that S. thermophilus can block ACR-containing phages when the CRISPR immunity specifically targets the acr gene. This in turn selects for phage mutants carrying a deletion within the acr gene. Remarkably, a truncated acrIIA allele, found in a wild-type virulent streptococcal phage, does not block the interference activity of Cas9 but still prevents the acquisition of new immunities, thereby providing an example of an ACR specifically inhibiting spacer acquisition.

Genome Sequence of SN1, a Bacteriophage That Infects Sphaerotilus natans and Pseudomonas aeruginosa.

Phage SN1 infects Sphaerotilus natans and Pseudomonas aeruginosa strains. Its genome consists of 61,858 bp (64.3% GC) and 89 genes, including 32 with predicted functions. SN1 genome is very similar to Pseudomonas phage M6, which contains hypermodified thymidines. Genome analyses revealed similar base-modifying genes as those found in M6.

In through the Out Door: A Functional Virulence Factor Secretion System Is Necessary for Phage Infection in Ralstonia solanacearum.

Bacteriophages put intense selective pressure on microbes, which must evolve diverse resistance mechanisms to survive continuous phage attacks. We used a library of spontaneous Bacteriophage Insensitive Mutants (BIMs) to learn how the plant pathogen Ralstonia solanacearum resists the virulent lytic podophage phiAP1. Phenotypic and genetic characterization of many BIMs suggested that the R. solanacearum Type II Secretion System (T2SS) plays a key role in phiAP1 infection. Using precision engineered mutations that permit T2SS assembly but either inactivate the T2SS GspE ATPase or sterically block the secretion portal, we demonstrated that phiAP1 needs a functional T2SS to infect R. solanacearum. This distinction between the static presence of T2SS components, which is necessary but not sufficient for phage sensitivity, and the energized and functional T2SS, which is sufficient, implies that binding interactions alone cannot explain the role of the T2SS in phiAP1 infection. Rather, our results imply that some aspect of the resetting of the T2SS, such as disassembly of the pseudopilus, is required. Because R. solanacearum secretes multiple virulence factors via the T2SS, acquiring resistance to phiAP1 also dramatically reduced R. solanacearum virulence on tomato plants. This acute fitness trade-off suggests this group of phages may be a sustainable control strategy for an important crop disease. IMPORTANCE Ralstonia solanacearum is a destructive plant pathogen that causes lethal bacterial wilt disease in hundreds of diverse plant hosts, including many economically important crops. Phages that kill R. solanacearum could offer effective and environmentally friendly wilt disease control, but only if the bacterium cannot easily evolve resistance. Encouragingly, most R. solanacearum mutants resistant to the virulent lytic phage phiAP1 no longer secreted multiple virulence factors and had much reduced fitness and virulence on tomato plants. Further analysis revealed that phage phiAP1 needs a functional type II secretion system to infect R. solanacearum, suggesting this podophage uses a novel infection mechanism.

Staphylococcus epidermidis bacteriophages from the anterior nares of humans.

The role of virulent bacteriophages in staphylococcal colonization of the human anterior nares is not known. This report of lytic bacteriophages against Staphylococcus epidermidis in the anterior nares of 5.5% of human subjects (n = 202) suggests their potential role in modulating staphylococcal colonization in this ecological niche.

Complete Genome Sequences of 10 Lactococcal Skunavirus Phages Isolated from Cheddar Cheese Whey Samples in Canada.

We report the complete genome sequences of 10 virulent phages of the Skunavirus genus (Siphoviridae) that infect Lactococcus lactis strains used for cheddar cheese production in Canada. Their linear genomes range from 28,969 bp to 31,042 bp with GC contents of 34.1 to 35.1% and 55 to 60 predicted open reading frames (ORFs).

A Lactococcal Phage Protein Promotes Viral Propagation and Alters the Host Proteomic Response During Infection.

The lactococcal virulent phage p2 is a model for studying the Skunavirus genus, the most prevalent group of phages causing milk fermentation failures in cheese factories worldwide. This siphophage infects Lactococcus lactis MG1363, a model strain used to study Gram-positive lactic acid bacteria. The structural proteins of phage p2 have been thoroughly described, while most of its non-structural proteins remain uncharacterized. Here, we developed an integrative approach, making use of structural biology, genomics, physiology, and proteomics to provide insights into the function of ORF47, the most conserved non-structural protein of unknown function among the Skunavirus genus. This small phage protein, which is composed of three α-helices, was found to have a major impact on the bacterial proteome during phage infection and to significantly reduce the emergence of bacteriophage-insensitive mutants.

Source Tracking Based on Core Genome SNV and CRISPR Typing of Salmonella enterica Serovar Heidelberg Isolates Involved in Foodborne Outbreaks in Québec, 2012.

Whole-genome sequencing (WGS) is the method of choice for bacterial subtyping and it is rapidly replacing the more traditional methods such as pulsed-field gel electrophoresis (PFGE). Here we used the high-resolution core genome single nucleotide variant (cgSNV) typing method to characterize clinical and food from Salmonella enterica serovar Heidelberg isolates in the context of source attribution. Additionally, clustered regularly interspaced short palindromic repeats (CRISPR) analysis was included to further support this method. Our results revealed that cgSNV was highly discriminatory and separated the outbreak isolates into distinct clusters (0-4 SNVs). CRISPR analysis was also able to distinguish outbreak strains from epidemiologically unrelated isolates. Specifically, our data clearly demonstrated the strength of these two methods to determine the probable source(s) of a 2012 epidemiologically characterized outbreak of S. Heidelberg. Using molecular cut-off of 0-10 SNVs, the cgSNV analysis of 246 clinical and food isolates of S. Heidelberg collected in Québec, in the same year of the outbreak event, revealed that retail and abattoir chicken isolates likely represent an important source of human infection to S. Heidelberg. Interestingly, the isolates genetically related by cgSNV also harbored the same CRISPR as outbreak isolates and clusters. This indicates that CRISPR profiles can be useful as a complementary approach to determine source attribution in foodborne outbreaks. Use of the genomic analysis also allowed to identify a large number of cases that were missed by PFGE, indicating that most outbreaks are probably underestimated. Although epidemiological information must still support WGS-based results, cgSNV method is a highly discriminatory method for the resolution of outbreak events and the attribution of these events to their respective sources. CRISPR typing can serve as a complimentary tool to this analysis during source tracking.

Author Correction: A mutation in the methionine aminopeptidase gene provides phage resistance in Streptococcus thermophilus.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.