Induction and Elimination of Prophages Using CRISPR Interference

Prophages are widely spread among bacterial genomes, and they can have positive or negative effects on their hosts. A key aspect in the study of prophages is the discovery of their induction signals. Prophage induction can occur by inactivating a phage transcriptional repressor, which is responsible for maintaining the lysogenic state. This repressor can be inactivated through the bacterial SOS response. However, the induction signals for numerous prophages do not involve the SOS system, and therefore significant efforts are needed to identify these conditions. Similarly, curing bacterial strains of inducible prophages is a tedious process, requiring the screening of several colonies. Here, we investigated whether transcriptional silencing of a prophage repressor using CRISPR interference (CRISPRi) would lead to prophage induction. Using Escherichia coli phages λ and P2 as models, we demonstrated the efficiency of CRISPRi for prophage induction and for curing lysogenic strains of their prophages.

Streamlining CRISPR spacer-based bacterial host predictions to decipher the viral dark matter

Thousands of new phages have recently been discovered thanks to viral metagenomics. These phages are extremely diverse and their genome sequences often do not resemble any known phages. To appreciate their ecological impact, it is important to determine their bacterial hosts. CRISPR spacers can be used to predict hosts of unknown phages, as spacers represent biological records of past phage–bacteria interactions. However, no guidelines have been established to standardize host prediction based on CRISPR spacers. Additionally, there are no tools that use spacers to perform host predictions on large viral datasets. Here, we developed a set of tools that includes all the necessary steps for predicting the hosts of uncharacterized phages. We created a database of >11 million spacers and a program to execute host predictions on large viral datasets. Our host prediction approach uses biological criteria inspired by how CRISPR–Cas naturally work as adaptive immune systems, which make the results easy to interpret. We evaluated the performance using 9484 phages with known hosts and obtained a recall of 49% and a precision of 69%. We also found that this host prediction method yielded higher performance for phages that infect gut-associated bacteria, suggesting it is well suited for gut-virome characterization.

Genomic diversity and CRISPR-Cas systems in the cyanobacterium Nostoc in the High Arctic

Nostoc (Nostocales, Cyanobacteria) has a global distribution in the Polar Regions. However, the genomic diversity of Nostoc is little known and there are no genomes available for polar Nostoc. Here we carried out the first genomic analysis of the Nostoc commune morphotype with a recent sample from the High Arctic and a herbarium specimen collected during the British Arctic Expedition (1875-76). Comparisons of the polar genomes with 26 present-day non-polar members of the Nostocales family highlighted that there are pronounced genetic variations among Nostoc strains and species. Osmoprotection and other stress genes were found in all Nostoc strains, but the two Arctic strains had markedly higher numbers of biosynthetic gene clusters for uncharacterised non-ribosomal peptide synthetases, suggesting a high diversity of secondary metabolites. Since viral-host interactions contribute to microbial diversity, we analysed the CRISPR-Cas systems in the Arctic and two temperate Nostoc species. There were a large number of unique repeat-spacer arrays in each genome, indicating diverse histories of viral attack. All Nostoc strains had a subtype I-D system, but the polar specimens also showed evidence of a subtype I-B system that has not been previously reported in cyanobacteria, suggesting diverse cyanobacteria-virus interactions in the Arctic.

Functional Study of the Type II-A CRISPR-Cas System of Streptococcus agalactiae Hypervirulent Strains

Nearly all strains of Streptococcus agalactiae, the leading cause of invasive infections in neonates, encode a type II-A clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system. Interestingly, S. agalactiae strains belonging to the hypervirulent Sequence Type 17 (ST17) contain significantly fewer spacers in their CRISPR locus than other lineages, which could be the result of a less functional CRISPR-Cas system. Here, we revealed one large deletion in the ST17 cas promoter region and we evaluated its impact on the transcription of cas genes as well as the functionalities of the CRISPR-Cas system. We demonstrated that Cas9 interference is functional and that the CRISPR-Cas system of ST17 strains can still acquire new spacers, despite the absence of a regular cas promoter. We demonstrated that a promoter sequence upstream of srn036, a small RNA partially overlapping the antisense tracrRNA, is responsible for the ST17 CRISPR-Cas adaptation and interference activities.

Ectopic Spacer Acquisition in Streptococcus thermophilus CRISPR3 Array

Streptococcus thermophilus relies heavily on two type II-A CRISPR-Cas systems, CRISPR1 and CRISPR3, to resist siphophage infections. One hallmark of these systems is the integration of a new spacer at the 5' end of the CRISPR arrays following phage infection. However, we have previously shown that ectopic acquisition of spacers can occur within the CRISPR1 array. Here, we present evidence of the acquisition of new spacers within the array of CRISPR3 of S. thermophilus. The analysis of randomly selected bacteriophage-insensitive mutants of the strain Uy01 obtained after phage infection, as well as the comparison with other S. thermophilus strains with similar CRISPR3 content, showed that a specific spacer within the array could be responsible for misguiding the adaptation complex. These results also indicate that while the vast majority of new spacers are added at the 5' end of the CRISPR array, ectopic spacer acquisition is a common feature of both CRISPR1 and CRISPR3 systems in S. thermophilus, and it can still provide phage resistance. Ectopic spacer acquisition also appears to have occurred naturally in some strains of Streptococcus pyogenes, suggesting that it is a general phenomenon, at least in type II-A systems.

The endless battle between phages and CRISPR-Cas systems in Streptococcus thermophilus

This review describes the contribution of basic research on phage-bacteria interactions to the understanding of CRISPR-Cas systems and their various applications. It focuses on the natural function of CRISPR-Cas systems as adaptive defense mechanisms against mobile genetic elements such as bacteriophage genomes and plasmids. Some of the advances in the characterization of the type II-A CRISPR-Cas system of Streptococcus thermophilus and Streptococcus pyogenes led to the development of the CRISPR-Cas9 genome-editing technology. We mostly discuss the 3 stages of the CRISPR-Cas system in S. thermophilus, namely the adaptation stage, which is unique to this resistance mechanism; the CRISPR RNA biogenesis; and the DNA-cutting activity in the interference stage to protect bacteria against phages. Finally, we look into applications of CRISPR-Cas in microbiology, including overcoming limitations in genome editing.

Delivery of CRISPR-Cas systems using phage-based vectors

Antimicrobial resistance has spread quickly on a worldwide scale, reducing therapeutic options for bacterial infections. CRISPR-Cas is an adaptive immune system found in many prokaryotes that can be designed to target bacterial genomes, leading to cell death. Repurposing the CRISPR-Cas system as a therapeutic strategy offers an attractive way to overcome antimicrobial resistance. However, this strategy requires efficient vectors for the CRISPR-Cas system to reach the bacterial genomes. Engineered phages offer an attractive option as cargo delivery vectors. In this review, we discuss the production of phage-based vectors and the relevance of using repurposed CRISPR-Cas systems as antimicrobials. We also discuss recent progress in phage engineering that can potentially overcome the limitations and increase the efficiency of CRISPR-Cas delivery.

Detection of preQ0 deazaguanine modifications in bacteriophage CAjan DNA using Nanopore sequencing reveals same hypermodification at two distinct DNA motifs

In the constant evolutionary battle against mobile genetic elements (MGEs), bacteria have developed several defense mechanisms, some of which target the incoming, foreign nucleic acids e.g. restriction-modification (R-M) or CRISPR-Cas systems. Some of these MGEs, including bacteriophages, have in turn evolved different strategies to evade these hurdles. It was recently shown that the siphophage CAjan and 180 other viruses use 7-deazaguanine modifications in their DNA to evade bacterial R-M systems. Among others, phage CAjan genome contains a gene coding for a DNA-modifying homolog of a tRNA-deazapurine modification enzyme, together with four 7-cyano-7-deazaguanine synthesis genes. Using the CRISPR-Cas9 genome editing tool combined with the Nanopore Sequencing (ONT) we showed that the 7-deazaguanine modification in the CAjan genome is dependent on phage-encoded genes. The modification is also site-specific and is found mainly in two separate DNA sequence contexts: GA and GGC. Homology modeling of the modifying enzyme DpdA provides insight into its probable DNA binding surface and general mode of DNA recognition.

Structural Insights into Lactococcal Siphophage p2 Baseplate Activation Mechanism

Virulent phages infecting L. lactis, an industry-relevant bacterium, pose a significant risk to the quality of the fermented milk products. Phages of the Skunavirus genus are by far the most isolated lactococcal phages in the cheese environments and phage p2 is the model siphophage for this viral genus. The baseplate of phage p2, which is used to recognize its host, was previously shown to display two conformations by X-ray crystallography, a rested state and an activated state ready to bind to the host. The baseplate became only activated and opened in the presence of Ca²⁺. However, such an activated state was not previously observed in the virion. Here, using nanobodies binding to the baseplate, we report on the negative staining electron microscopy structure of the activated form of the baseplate directly observed in the p2 virion, that is compatible with the activated baseplate crystal structure. Analyses of this new structure also established the presence of a second distal tail (Dit) hexamer as a component of the baseplate, the topology of which differs largely from the first one. We also observed an uncoupling between the baseplate activation and the tail tip protein (Tal) opening, suggesting an infection mechanism more complex than previously expected.

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.

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.

The EcoChip 2: An Autonomous Sensor Platform for Multimodal Bio-environmental Monitoring of the Northern Habitat

This paper presents the EcoChip 2, an autonomous multimodal bio-environmental sensor platform for the monitoring of microorganisms in the northern habitat. The EcoChip 2 prototype includes an array of 96-wells for the continuous monitoring of microbiological growth through a multichannel electrochemical impedance analyzer circuit. In addition, the platform includes luminosity, humidity, temperature sensors and monitoring. The developed electronic board uses an ultra-low-power microcontroller unit, a custom power management unit, a low-power wireless ISM-2.45 GHz transceiver, and a flash memory to accumulate and store the sensor data over extended monitoring periods. When a wireless base station is placed within the transmission range of the EcoChip 2, an embedded low-power wireless transceiver transmits the 96-wells impedance data and the other sensor data stored in the flash memory to the user interface. We present the measured performance of the prototype, along with laboratory test results of bacterial growth measurements inside the 96 wells in parallel. We show that the EcoChip 2 can successfully measure the impedances associated with bacterial growth over several hours using an excitation frequency of 2 kHz with power consumption of 114.6 mW under operating mode.

A short overview of the CRISPR-Cas adaptation stage

CRISPR research began over 30 years ago with the incidental discovery of an unusual nucleotide arrangement in the Escherichia coli genome. It took 20 years to find the main function of CRISPR-Cas systems as an adaptive defence mechanism against invading nucleic acids, and our knowledge of their biology has steadily increased ever since. In parallel, the number of applications derived from CRISPR-Cas systems has risen spectacularly. The CRISPR-based genome editing tool is arguably the most exciting application in both basic and applied research. Lately, CRISPR-Cas research has partially shifted to the least understood aspect of its biology: the ability of CRISPR-Cas systems to acquire new immunities during the so-called adaptation step. To date, the most efficient natural system to readily acquire new spacers is the type II-A system of the gram-positive dairy bacterium Streptococcus thermophilus. The discovery of additional systems able to acquire new spacers will hopefully draw more attention to this step of CRISPR-Cas biology. This review focuses on the breakthroughs that have helped to unravel the adaptation phase and on questions that remain to be answered.

Comparative genomic analysis of 142 bacteriophages infecting Salmonella enterica subsp. enterica

BACKGROUND

Bacteriophages are bacterial parasites and are considered the most abundant and diverse biological entities on the planet. Previously we identified 154 prophages from 151 serovars of Salmonella enterica subsp. enterica. A detailed analysis of Salmonella prophage genomics is required given the influence of phages on their bacterial hosts and should provide a broader understanding of Salmonella biology and virulence and contribute to the practical applications of phages as vectors and antibacterial agents.

RESULTS

Here we provide a comparative analysis of the full genome sequences of 142 prophages of Salmonella enterica subsp. enterica which is the full complement of the prophages that could be retrieved from public databases. We discovered extensive variation in genome sizes (ranging from 6.4 to 358.7 kb) and guanine plus cytosine (GC) content (ranging from 35.5 to 65.4%) and observed a linear correlation between the genome size and the number of open reading frames (ORFs). We used three approaches to compare the phage genomes. The NUCmer/MUMmer genome alignment tool was used to evaluate linkages and correlations based on nucleotide identity between genomes. Multiple sequence alignment was performed to calculate genome average nucleotide identity using the Kalgin program. Finally, genome synteny was explored using dot plot analysis. We found that 90 phage genome sequences grouped into 17 distinct clusters while the remaining 52 genomes showed no close relationships with the other phage genomes and are identified as singletons. We generated genome maps using nucleotide and amino acid sequences which allowed protein-coding genes to be sorted into phamilies (phams) using the Phamerator software. Out of 5796 total assigned phamilies, one phamily was observed to be dominant and was found in 49 prophages, or 34.5% of the 142 phages in our collection. A majority of the phamilies, 4330 out of 5796 (74.7%), occurred in just one prophage underscoring the high degree of diversity among Salmonella bacteriophages.

CONCLUSIONS

Based on nucleotide and amino acid sequences, a high diversity was found among Salmonella bacteriophages which validate the use of prophage sequence analysis as a highly discriminatory subtyping tool for Salmonella. Thorough understanding of the conservation and variation of prophage genomic characteristics will facilitate their rational design and use as tools for bacterial strain construction, vector development and as anti-bacterial agents.

Variability in the durability of CRISPR-Cas immunity

The durability of host resistance is challenged by the ability of pathogens to escape the defence of their hosts. Understanding the variability in the durability of host resistance is of paramount importance for designing more effective control strategies against infectious diseases. Here, we study the durability of various clustered regularly interspaced short palindromic repeats-Cas (CRISPR-Cas) alleles of the bacteria Streptococcus thermophilus against lytic phages. We found substantial variability in durability among different resistant bacteria. Since the escape of the phage is driven by a mutation in the phage sequence targeted by CRISPR-Cas, we explored the fitness costs associated with these escape mutations. We found that, on average, escape mutations decrease the fitness of the phage. Yet, the magnitude of this fitness cost does not predict the durability of CRISPR-Cas immunity. We contend that this variability in the durability of resistance may be because of variations in phage mutation rate or in the proportion of lethal mutations across the phage genome. These results have important implications on the coevolutionary dynamics between bacteria and phages and for the optimal deployment of resistance strategies against pathogens and pests. Understanding the durability of CRISPR-Cas immunity may also help develop more effective gene-drive strategies based on CRISPR-Cas9 technology. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Phage diversity, genomics and phylogeny

Recent advances in viral metagenomics have enabled the rapid discovery of an unprecedented catalogue of phages in numerous environments, from the human gut to the deep ocean. Although these advances have expanded our understanding of phage genomic diversity, they also revealed that we have only scratched the surface in the discovery of novel viruses. Yet, despite the remarkable diversity of phages at the nucleotide sequence level, the structural proteins that form viral particles show strong similarities and conservation. Phages are uniquely interconnected from an evolutionary perspective and undergo multiple events of genetic exchange in response to the selective pressure of their hosts, which drives their diversity. In this Review, we explore phage diversity at the structural, genomic and community levels as well as the complex evolutionary relationships between phages, moulded by the mosaicity of their genomes.

Virulent coliphages in 1-year-old children fecal samples are fewer, but more infectious than temperate coliphages

Bacteriophages constitute an important part of the human gut microbiota, but their impact on this community is largely unknown. Here, we cultivate temperate phages produced by 900 E. coli strains isolated from 648 fecal samples from 1-year-old children and obtain coliphages directly from the viral fraction of the same fecal samples. We find that 63% of strains hosted phages, while 24% of the viromes contain phages targeting E. coli. 150 of these phages, half recovered from strain supernatants, half from virome (73% temperate and 27% virulent) were tested for their host range on 75 E. coli strains isolated from the same cohort. Temperate phages barely infected the gut strains, whereas virulent phages killed up to 68% of them. We conclude that in fecal samples from children, temperate coliphages dominate, while virulent ones have greater infectivity and broader host range, likely playing a role in gut microbiota dynamics.

The impact of bacteriophages in the human gut microbiome remains poorly understood. Here, the authors characterize coliphages isolated from a large cohort of 1-year-old infants and show that temperate coliphages dominate, while virulent ones have greater infectivity and broader host range.

Versatile and robust genome editing with Streptococcus thermophilus CRISPR1-Cas9

Targeting definite genomic locations using CRISPR-Cas systems requires a set of enzymes with unique protospacer adjacent motif (PAM) compatibilities. To expand this repertoire, we engineered nucleases, cytosine base editors, and adenine base editors from the archetypal Streptococcus thermophilus CRISPR1-Cas9 (St1Cas9) system. We found that St1Cas9 strain variants enable targeting to five distinct A-rich PAMs and provide a structural basis for their specificities. The small size of this ortholog enables expression of the holoenzyme from a single adeno-associated viral vector for in vivo editing applications. Delivery of St1Cas9 to the neonatal liver efficiently rewired metabolic pathways, leading to phenotypic rescue in a mouse model of hereditary tyrosinemia. These robust enzymes expand and complement current editing platforms available for tailoring mammalian genomes.

7-Deazaguanine modifications protect phage DNA from host restriction systems

Genome modifications are central components of the continuous arms race between viruses and their hosts. The archaeosine base (G+), which was thought to be found only in archaeal tRNAs, was recently detected in genomic DNA of Enterobacteria phage 9g and was proposed to protect phage DNA from a wide variety of restriction enzymes. In this study, we identify three additional 2'-deoxy-7-deazaguanine modifications, which are all intermediates of the same pathway, in viruses: 2'-deoxy-7-amido-7-deazaguanine (dADG), 2'-deoxy-7-cyano-7-deazaguanine (dPreQ0) and 2'-deoxy-7- aminomethyl-7-deazaguanine (dPreQ1). We identify 180 phages or archaeal viruses that encode at least one of the enzymes of this pathway with an overrepresentation (60%) of viruses potentially infecting pathogenic microbial hosts. Genetic studies with the Escherichia phage CAjan show that DpdA is essential to insert the 7-deazaguanine base in phage genomic DNA and that 2'-deoxy-7-deazaguanine modifications protect phage DNA from host restriction enzymes.

A mutation in the methionine aminopeptidase gene provides phage resistance in Streptococcus thermophilus

Streptococcus thermophilus is a lactic acid bacterium widely used by the dairy industry for the manufacture of yogurt and specialty cheeses. It is also a Gram-positive bacterial model to study phage-host interactions. CRISPR-Cas systems are one of the most prevalent phage resistance mechanisms in S. thermophilus. Little information is available about other host factors involved in phage replication in this food-grade streptococcal species. We used the model strain S. thermophilus SMQ-301 and its virulent phage DT1, harboring the anti-CRISPR protein AcrIIA6, to show that a host gene coding for a methionine aminopeptidase (metAP) is necessary for phage DT1 to complete its lytic cycle. A single mutation in metAP provides S. thermophilus SMQ-301 with strong resistance against phage DT1. The mutation impedes a late step of the lytic cycle since phage adsorption, DNA replication, and protein expression were not affected. When the mutated strain was complemented with the wild-type version of the gene, the phage sensitivity phenotype was restored. When this mutation was introduced into other S. thermophilus strains it provided resistance against cos-type (Sfi21dt1virus genus) phages but replication of pac-type (Sfi11virus genus) phages was not affected. The mutation in the gene coding for the MetAP induces amino acid change in a catalytic domain conserved across many bacterial species. Introducing the same mutation in Streptococcus mutans also provided a phage resistance phenotype, suggesting the wide-ranging importance of the host methionine aminopeptidase in phage replication.

Diversity and Host Specificity Revealed by Biological Characterization and Whole Genome Sequencing of Bacteriophages Infecting Salmonella enterica

Phages infecting members of the opportunistic human pathogen, Salmonella enterica, are widespread in natural environments and offer a potential source of agents that could be used for controlling populations of this bacterium; yet, relatively little is known about these phages. Here we describe the isolation and characterization of 45 phages of Salmonella enterica from disparate geographic locations within British Columbia, Canada. Host-range profiling revealed host-specific patterns of susceptibility and resistance, with several phages identified that have a broad-host range (i.e., able to lyse >40% of bacterial hosts tested). One phage in particular, SE13, is able to lyse 51 out of the 61 Salmonella strains tested. Comparative genomic analyses also revealed an abundance of sequence diversity in the sequenced phages. Alignment of the genomes grouped the phages into 12 clusters with three singletons. Phages within certain clusters exhibited extraordinarily high genome homology (>98% nucleotide identity), yet between clusters, genomes exhibited a span of diversity (<50% nucleotide identity). Alignment of the major capsid protein also supported the clustering pattern observed with alignment of the whole genomes. We further observed associations between genomic relatedness and the site of isolation, as well as genetic elements related to DNA metabolism and host virulence. Our data support the knowledge framework for phage diversity and phage-host interactions that are required for developing phage-based applications for various sectors, including biocontrol, detection and typing.

Mobilome of Brevibacterium aurantiacum Sheds Light on Its Genetic Diversity and Its Adaptation to Smear-Ripened Cheeses

Brevibacterium aurantiacum is an actinobacterium that confers key organoleptic properties to washed-rind cheeses during the ripening process. Although this industrially relevant species has been gaining an increasing attention in the past years, its genome plasticity is still understudied due to the unavailability of complete genomic sequences. To add insights on the mobilome of this group, we sequenced the complete genomes of five dairy Brevibacterium strains and one non-dairy strain using PacBio RSII. We performed phylogenetic and pan-genome analyses, including comparisons with other publicly available Brevibacterium genomic sequences. Our phylogenetic analysis revealed that these five dairy strains, previously identified as Brevibacterium linens, belong instead to the B. aurantiacum species. A high number of transposases and integrases were observed in the Brevibacterium spp. strains. In addition, we identified 14 and 12 new insertion sequences (IS) in B. aurantiacum and B. linens genomes, respectively. Several stretches of homologous DNA sequences were also found between B. aurantiacum and other cheese rind actinobacteria, suggesting horizontal gene transfer (HGT). A HGT region from an iRon Uptake/Siderophore Transport Island (RUSTI) and an iron uptake composite transposon were found in five B. aurantiacum genomes. These findings suggest that low iron availability in milk is a driving force in the adaptation of this bacterial species to this niche. Moreover, the exchange of iron uptake systems suggests cooperative evolution between cheese rind actinobacteria. We also demonstrated that the integrative and conjugative element BreLI (Brevibacterium Lanthipeptide Island) can excise from B. aurantiacum SMQ-1417 chromosome. Our comparative genomic analysis suggests that mobile genetic elements played an important role into the adaptation of B. aurantiacum to cheese ecosystems.