Impact of CRISPR immunity on the emergence and virulence of bacterial pathogens

The role of noncoding RNAs in bacterial immunity

The evolutionary arms race between bacteria and phages has driven the development of diverse anti-phage defense mechanisms. Recent studies have identified noncoding RNAs (ncRNAs) as key players in bacteria-phage conflicts, including CRISPR-Cas, toxin-antitoxin (TA), and reverse transcriptase (RT)-based defenses; however, our understanding of their roles in immunity is still emerging. In this review, we explore the multifaceted roles of ncRNAs in bacterial immunity, offering insights into their contributions to defense and anti-defense mechanisms, their influence on immune regulatory networks, and potential biotechnological applications. Finally, we highlight key outstanding questions in the field to spark future research directions.

Role of CRISPR-cas system on virulence traits and carbapenem resistance in clinical Klebsiella pneumoniae isolates

BACKGROUND AND OBJECTIVES

The bacterial adaptive immune system known as CRISPR-Cas (clustered regularly interspersed short palindromic repeats-CRISPR-associated protein) is engaged in defense against various mobile genetic elements (MGEs) such as plasmids and bacteriophages. The purpose of this study was to characterize the CRISPR-Cas systems in carbapenem-resistant Klebsiella pneumoniae isolates and assess any possible correlation between these systems with antibiotic susceptibility, biofilm formation, and bacterial virulence.

MATERIALS AND METHODS

A total of 156 CRKP isolates were collected from different specimens of the inpatients. Biofilm formation and antibiotic susceptibility testing were evaluated using standard methods. Furthermore, the CRISPR-Cas system subtype genes, 11 carbapenemase genes, and 17 virulence genes were identified using separate standard PCR reactions. The diversity of the isolates was determined by random amplified polymorphic DNA (RAPD)-PCR.

RESULTS

The development of biofilms and antibiotic susceptibility of several CRKP isolates were significantly correlated with the absence or presence of the CRISPR-Cas system. PCR analysis of carbapenemase genes revealed that the frequency of the blaNDM-1 gene was significantly higher in the isolates with the subtype I-E CRISPR-Cas system. Moreover, the isolates with the subtype I-E CRISPR-Cas system exhibited a propensity to possess more virulence genes such as allS, k2A, wcaG, aerobactin, rmpA, iroN, magA, rmpA2, kfu, iutA, iucB, ybtS, repA, and terW.

CONCLUSION

CRISPR-Cas systems could affect the antibiotic susceptibility, capacity for biofilm formation, and virulence of Klebsiella pneumoniae. Our findings showed that the isolates containing the CRISPR-Cas system were moderate or strong biofilm producers and had a higher frequency of virulence genes.

Prevalence of the CRISPR-cas system and its association with antibiotic resistance in clinical Klebsiella pneumoniae isolates

BACKGROUND AND OBJECTIVE(S)

CRISPR-Cas is a prokaryotic adaptive immune system that protects bacteria and archaea against mobile genetic elements (MGEs) such as bacteriophages plasmids, and transposons. In this study, we aimed to assess the prevalence of the CRISPR-Cas systems and their association with antibiotic resistance in one of the most challenging bacterial pathogens, Klebsiella pneumoniae.

MATERIALS AND METHODS

A total of 105 K. pneumoniae isolates were collected from various clinical infections. Extended-spectrum β-lactamases (ESBLs) phenotypically were detected and the presence of ESBL, aminoglycoside-modifying enzymes (AME), and CRISPR-Cas system subtype genes were identified using PCR. Moreover, the diversity of the isolates was determined by enterobacterial repetitive intergenic consensus (ERIC)-PCR.

RESULTS

Phenotypically, 41.9% (44/105) of the isolates were found to be ESBL producers. A significant inverse correlation existed between the subtype I-E CRISPR-Cas system's presence and ESBL production in K. pneumoniae isolates. Additionally, the frequency of the ESBL genes blaCTX-M1 (3%), blaCTX-M9 (12.1%), blaSHV (51.5%), and blaTEM (33.3%), as well as some AME genes such as aac(3)-Iva (21.2%) and ant(2'')-Ia (3%) was significantly lower in the isolates with the subtype I-E CRISPR-Cas system in comparison to CRISPR-negative isolates. There was a significant inverse correlation between the presence of ESBL and some AME genes with subtype I-E CRISPR-Cas system.

CONCLUSION

The presence of the subtype I-E CRISPR-Cas system was correlated with the antibiotic-resistant gene (ARGs). The isolates with subtype I-E CRISPR-Cas system had a lower frequency of ESBL genes and some AME genes than CRISPR-negative isolates.

Exploring the comparative genome of rice pathogen Burkholderia plantarii: unveiling virulence, fitness traits, and a potential type III secretion system effector

This study presents a comprehensive genomic analysis of Burkholderia plantarii, a rice pathogen that causes blight and grain rot in seedlings. The entire genome of B. plantarii KACC 18964 was sequenced, followed by a comparative genomic analysis with other available genomes to gain insights into its virulence, fitness, and interactions with rice. Multiple secondary metabolite gene clusters were identified. Among these, 12 demonstrated varying similarity levels to known clusters linked to bioactive compounds, whereas eight exhibited no similarity, indicating B. plantarii as a source of potentially novel secondary metabolites. Notably, the genes responsible for tropolone and quorum sensing were conserved across the examined genomes. Additionally, B. plantarii was observed to possess three complete CRISPR systems and a range of secretion systems, exhibiting minor variations among the analyzed genomes. Genomic islands were analyzed across the four genomes, and a detailed study of the B. plantarii KACC 18964 genome revealed 59 unique islands. These islands were thoroughly investigated for their gene contents and potential roles in virulence. Particular attention has been devoted to the Type III secretion system (T3SS), a crucial virulence factor. An in silico analysis of potential T3SS effectors identified a conserved gene, aroA. Further mutational studies, in planta and in vitro analyses validated the association between aroA and virulence in rice. Overall, this study enriches our understanding of the genomic basis of B. plantarii pathogenicity and emphasizes the potential role of aroA in virulence. This understanding may guide the development of effective disease management strategies.

Analysis of twelve genomes of the bacterium Kerstersia gyiorum from brown-throated sloths (Bradypus variegatus), the first from a non-human host

Kerstersia gyiorum is a Gram-negative bacterium found in various animals, including humans, where it has been associated with various infections. Knowledge of the basic biology of K. gyiorum is essential to understand the evolutionary strategies of niche adaptation and how this organism contributes to infectious diseases; however, genomic data about K. gyiorum is very limited, especially from non-human hosts. In this work, we sequenced 12 K. gyiorum genomes isolated from healthy free-living brown-throated sloths (Bradypus variegatus) in the Parque Estadual das Fontes do Ipiranga (São Paulo, Brazil), and compared them with genomes from isolates of human origin, in order to gain insights into genomic diversity, phylogeny, and host specialization of this species. Phylogenetic analysis revealed that these K. gyiorum strains are structured according to host. Despite the fact that sloth isolates were sampled from a single geographic location, the intra-sloth K. gyiorum diversity was divided into three clusters, with differences of more than 1,000 single nucleotide polymorphisms between them, suggesting the circulation of various K. gyiorum lineages in sloths. Genes involved in mobilome and defense mechanisms against mobile genetic elements were the main source of gene content variation between isolates from different hosts. Sloth-specific K. gyiorum genome features include an IncN2 plasmid, a phage sequence, and a CRISPR-Cas system. The broad diversity of defense elements in K. gyiorum (14 systems) may prevent further mobile element flow and explain the low amount of mobile genetic elements in K. gyiorum genomes. Gene content variation may be important for the adaptation of K. gyiorum to different host niches. This study furthers our understanding of diversity, host adaptation, and evolution of K. gyiorum, by presenting and analyzing the first genomes of non-human isolates.

Pediococcus inopinatus with a well-developed CRISPR-Cas system dominates in long-term fermented kimchi, Mukeunji

Kimchi is produced through a low-temperature fermentation without pre-sterilization, resulting in a heterogeneous microbial community. As fermentation progresses, dominant lactic acid bacteria (LAB) species emerge and undergo a transition process. In this study, LAB were isolated from Mukeunji, a long-term fermented kimchi that is in the final stage of kimchi fermentation process. It was confirmed, through culture-dependent and independent analysis, as well as metagenome analysis, that Pediococcus inopinatus are generally dominant in long-term fermented kimchi. Comparative analysis of the de novo assembled whole genome of P. inopinatus with other kimchi LAB revealed that this species has a well-developed clustered regularly interspaced short palindromic repeats (CRISPR) system. The CRISPR system of P. inopinatus has an additional copy of the csa3 gene, a transcription factor for cas genes. Indeed, this species not only highly expresses cas1 and cas2, which induce spacer acquisition, but also has many diverse spacers that are actively expressed. These findings indicate that the well-developed CRISPR-Cas system is enabling P. inopinatus to dominate in long-fermented kimchi. Overall, this study revealed that LAB with a robust defense system dominate in the final stage of kimchi fermentation and presented a model for the succession mechanism of kimchi LAB.

Exploring the Interplay of the CRISPR-CAS System with Antibiotic Resistance in Staphylococcus aureus: A Poultry Meat Study from Lahore, Pakistan

Staphylococcus aureus is one of the major pathogens responsible for causing food poisoning worldwide. The emergence of antibiotic resistance in this bacterium is influenced by various factors. Among them, bacterial acquired defense systems described as clustered regularly interspaced short palindromic repeats (CRISPR)-cas system might be involved in antibiotic resistance development in bacteria. The current study was designed to assess the prevalence of S. aureus and its antibiotic resistance profile and identify the relationship of the CRISPR-cas system with antimicrobial resistance, followed by phylogenetic analysis. Total samples (n = 188) of poultry meat were collected from the poultry bird market of Lahore, Punjab, Pakistan. We used both phenotypic (antibiotic disc diffusion) and genotypic methods (PCR) to identify multi-drug resistant (MDR) strains of S. aureus. Additionally, the role of the CRISPR-Cas system in the isolated MDR S. aureus was also assessed. In addition, real-time quantitative PCR (qRT-PCR) was used to evaluate the association of the CRISPR-cas system with antimicrobial resistance. All of the S. aureus isolates showed 100% resistance against erythromycin, 97.5% were resistant to tetracycline, and 75% were resistant to methicillin. Eleven isolates were MDR in the current study. The CRISPR system was found in all MDR isolates, and fifteen spacers were identified within the CRISPR locus. Furthermore, MDR S. aureus isolates and the standard strain showed higher expression levels of CRISPR-associated genes. The correlation of said system with MDR isolates points to foreign gene acquisition by horizontal transfer. Current knowledge could be utilized to tackle antibiotic-resistant bacteria, mainly S. aureus.

Advances in CRISPR/Cas-Based Strategies on Zoonosis

Clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated (Cas) has emerged as the predominant technique for gene editing technique due to its high efficiency and low cost. In the area of zoonosis, CRISPR/Cas is also widely used in different research areas. This paper reviewed the principles of CRISPR/Cas technique and its applications in zoonosis. Moreover, we analyze the shortcomings and weaknesses that currently limit its use, highlight its direction for improvement, and foresee its application prospects in the prevention and treatment of zoonosis. For the purpose of preventing and controlling zoonosis, we need to develop diagnostic method with high sensitivity and specificity, highly protective vaccines, and also better understanding of the pathogenesis. Our review aimed to promote the application and improvement of CRISPR/Cas technique in the above-mentioned areas, and provide brief and comprehensive references for CRISPR/Cas-based research. Through reviewing the advances in CRISPR/Cas-based strategies on zoonosis, we believe that CRISPR/Cas technique will provide more powerful assistance for the prevention and control of zoonosis.

Correlation of CRISPR/Cas and Antimicrobial Resistance in Klebsiella pneumoniae Clinical Isolates Recovered from Patients in Egypt Compared to Global Strains

The CRISPR/Cas system has been long known to interfere with the acquisition of foreign genetic elements and was recommended as a tool for fighting antimicrobial resistance. The current study aimed to explore the prevalence of the CRISPR/Cas system in Klebsiella pneumoniae isolates recovered from patients in Egypt in comparison to global strains and correlate the CRISPR/Cas to susceptibility to antimicrobial agents. A total of 181 clinical isolates were PCR-screened for cas and selected antimicrobial resistance genes (ARGs). In parallel, 888 complete genome sequences were retrieved from the NCBI database for in silico analysis. CRISPR/Cas was found in 46 (25.4%) isolates, comprising 18.8% type I-E and 6.6% type I-E*. Multidrug resistance (MDR) and extensive drug resistance (XDR) were found in 73.5% and 25.4% of the isolates, respectively. More than 95% of the CRISPR/Cas-bearing isolates were MDR (65.2%) or XDR (32.6%). No significant difference was found in the susceptibility to the tested antimicrobial agents among the CRISPR/Cas-positive and -negative isolates. The same finding was obtained for the majority of the screened ARGs. Among the published genomes, 23.2% carried CRISPR/Cas, with a higher share of I-E* (12.8%). They were confined to specific sequence types (STs), most commonly ST147, ST23, ST15, and ST14. More plasmids and ARGs were carried by the CRISPR/Cas-negative group than others, but their distribution in the two groups was not significantly different. The prevalence of some ARGs, such as blaKPC, blaTEM, and rmtB, was significantly higher among the genomes of the CRISPR/Cas-negative strains. A weak, nonsignificant positive correlation was found between the number of spacers and the number of resistance plasmids and ARGs. In conclusion, the correlation between CRISPR/Cas and susceptibility to antimicrobial agents or bearing resistance plasmids and ARGs was found to be nonsignificant. Plasmid-targeting spacers might not be naturally captured by CRISPR/Cas. Spacer match analysis is recommended to provide a clearer image of the exact behavior of CRISPR/Cas towards resistance plasmids.

An expectation-maximization algorithm for estimating proportions of deletions among bacterial populations with application to study antibiotic resistance gene transfer in Enterococcus faecalis

The emergence of antibiotic resistance in bacteria limits the availability of antibiotic choices for treatment and infection control, thereby representing a major threat to human health. The de novo mutation of bacterial genomes is an essential mechanism by which bacteria acquire antibiotic resistance. Previously, deletion mutations within bacterial immune systems, ranging from dozens to thousands of base pairs (bps) in length, have been associated with the spread of antibiotic resistance. Most current methods for evaluating genomic structural variations (SVs) have concentrated on detecting them, rather than estimating the proportions of populations that carry distinct SVs. A better understanding of the distribution of mutations and subpopulations dynamics in bacterial populations is needed to appreciate antibiotic resistance evolution and movement of resistance genes through populations. Here, we propose a statistical model to estimate the proportions of genomic deletions in a mixed population based on Expectation-Maximization (EM) algorithms and next-generation sequencing (NGS) data. The method integrates both insert size and split-read mapping information to iteratively update estimated distributions. The proposed method was evaluated with three simulations that demonstrated the production of accurate estimations. The proposed method was then applied to investigate the horizontal transfers of antibiotic resistance genes in concert with changes in the CRISPR-Cas system of E. faecalis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00144-z.

Prevalence and analysis of CRISPR/Cas systems in Pseudomonas aeruginosa isolates from Greece

The objective of this study was to investigate the prevalence of CRISPR/Cas systems in P. aeruginosa, isolated from a Greek hospital. Additionally, we aimed to determine the origin of the sequenced spacers. A collection of 100 nonrepetitive P. aeruginosa was analyzed. Isolates were typed by MLST. The presence of CRISPR/Cas systems, as well as amplification of CRISPR arrays, was examined by PCR using specific primers. CRISPR/Cas systems were detected in 36 isolates, of which 27 isolates exhibited resistance to carbapenems, with 10 of the later isolates producing a VIM-type MβL. The majority (n = 19) of CRISPR/Cas-positive isolates harbored a type I-F system, while I-C and I-E systems were found in 9 and 8 isolates, respectively. Based on MLST, isolates carrying I-E and I-F systems belonged to different STs and included CRISPR arrays with diverse number of spacers. Isolates with I-C systems belonged to clonal complex 235 and exhibited identical CRISPR arrays. Among 425 unique spacers, identified during this study, BLASTn search showed that they matched with P. aeruginosa chromosomal sequences (47.0%), phages (31.9%), plasmids, PAGIs, and an ICE. 16.3% of the spacers exhibited no significant similarity with sequences submitted to GenBank database. In conclusion, we observed the presence of type I-C, I-E and I-F CRISPR/Cas systems in P. aeruginosa of clinical origin. CRISPR/Cas were also observed among isolates carrying the carbapenemase-encoding bla_VIM gene, which is usually associated with integrons, questioning the defense role against mobile elements. Therefore, further experimental characterization is needed to clarify their functional role.

Genomic Insights Into the Interspecific Diversity and Evolution of Mobiluncus, a Pathogen Associated With Bacterial Vaginosis

Bacterial vaginosis (BV) is a common vaginal infection and has been associated with increased risk for a wide array of health issues. BV is linked with a variety of heterogeneous pathogenic anaerobic bacteria, among which Mobiluncus is strongly associated with BV diagnosis. However, their genetic features, pathogenicity, interspecific diversity, and evolutionary characters have not been illustrated at genomic level. The current study performed phylogenomic and comparative genomic analyses of Mobiluncus. Phylogenomic analyses revealed remarkable phylogenetic distinctions among different species. Compared with M. curtisii, M. mulieris had a larger genome and pangenome size with more insertion sequences but less CRISPR-Cas systems. In addition, these two species were diverse in profile of virulence factors, but harbored similar antibiotic resistance genes. Statistically different functional genome profiles between strains from the two species were determined, as well as correlations of some functional genes/pathways with putative pathogenicity. We also showed that high levels of horizontal gene transfer might be an important strategy for species diversification and pathogenicity. Collectively, this study provides the first genome sequence level description of Mobiluncus, and may shed light on its virulence/pathogenicity, functional diversification, and evolutionary dynamics. Our study could facilitate the further investigations of this important pathogen, and might improve the future treatment of BV.

CRISPR/Cas9 a simple, inexpensive and effective technique for gene editing

In recent years, the number of tools and techniques that enable genetic material to be added, removed or altered at specific locations in the genome has increased significantly. The objective is to know the structure of genomes, the function of genes and improve gene therapy.In this work we intend to explain the functioning of the CRISPR/Cas9 (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR associated protein 9) and the advantages that this technique may have compared to previously developed techniques, such as RNA interference (RNAi), Zinc Finger Nucleases (ZFNs) and transcription activator-like effector nucleases (TALENs) in gene and genome editing.We will start with the story of the discovery, then its biological function in the adaptive immune system of bacteria against bacteriophage attack, and ending with a description of the mechanism of action and its use in gene editing. We will also discuss other Cas enzymes with great potential for use in genome editing as an alternative to Cas9.CRISPR/Cas9 is a simple, inexpensive, and effective technique for gene editing with multiple applications from the development of functional genomics and epigenetics. This technique will, in the near future, have great applications in the development of cell models for use in medical and pharmaceutical processes, in targeted therapy, and improvement of agricultural and environmental species.

Comparative Genomics Provides Insights Into Genetic Diversity of Clostridium tyrobutyricum and Potential Implications for Late Blowing Defects in Cheese

Clostridium tyrobutyricum has been recognized as the main cause of late blowing defects (LBD) in cheese leading to considerable economic losses for the dairy industry. Although differences in spoilage ability among strains of this species have been acknowledged, potential links to the genetic diversity and functional traits remain unknown. In the present study, we aimed to investigate and characterize genomic variation, pan-genomic diversity and key traits of C. tyrobutyricum by comparing the genomes of 28 strains. A comparative genomics analysis revealed an “open” pangenome comprising 9,748 genes and a core genome of 1,179 genes shared by all test strains. Among those core genes, the majority of genes encode proteins related to translation, ribosomal structure and biogenesis, energy production and conversion, and amino acid metabolism. A large part of the accessory genome is composed of sets of unique, strain-specific genes ranging from about 5 to more than 980 genes. Furthermore, functional analysis revealed several strain-specific genes related to replication, recombination and repair, cell wall, membrane and envelope biogenesis, and defense mechanisms that might facilitate survival under stressful environmental conditions. Phylogenomic analysis divided strains into two clades: clade I contained human, mud, and silage isolates, whereas clade II comprised cheese and milk isolates. Notably, these two groups of isolates showed differences in certain hypothetical proteins, transcriptional regulators and ABC transporters involved in resistance to oxidative stress. To the best of our knowledge, this is the first study to provide comparative genomics of C. tyrobutyricum strains related to LBD. Importantly, the findings presented in this study highlight the broad genetic diversity of C. tyrobutyricum, which might help us understand the diversity in spoilage potential of C. tyrobutyricum in cheese and provide some clues for further exploring the gene modules responsible for the spoilage ability of this species.

Distribution of CRISPR in Escherichia coli Isolated from Bulk Tank Milk and Its Potential Relationship with Virulence

Escherichia coli is one of the most common causes of mastitis on dairy farms around the world, but its clinical severity is determined by a combination of virulence factors. Recently, clustered regularly interspaced short palindromic repeat (CRISPR) arrays have been reported as a novel typing method because of their usefulness in discriminating pathogenic bacterial isolates. Therefore, this study aimed to investigate the virulence potential of E. coli isolated from bulk tank milk, not from mastitis, and to analyze its pathogenic characterization using the CRISPR typing method. In total, 164 (89.6%) out of 183 E. coli isolated from the bulk tank milk of 290 farms carried one or more of eighteen virulence genes. The most prevalent virulence gene was fimH (80.9%), followed by iss (38.3%), traT (26.8%), ompT (25.7%), afa/draBC (24.0%), and univcnf (21.9%). Moreover, the phylogenetic group with the highest prevalence was B1 (64.0%), followed by A (20.1%), D (8.5%), and C (7.3%) (p < 0.05). Among the four CRISPR loci, only two, CRISPR 1 and CRISPR 2, were found. Interestingly, the distribution of CRISPR 1 was significantly higher in groups A and B1 compared to that of CRISPR 2 (p < 0.05), but there were no significant differences in groups C and D. The prevalence of CRISPR 1 by virulence gene ranged from 91.8% to 100%, whereas that of CRISPR 2 ranged from 57.5% to 93.9%. The distribution of CRISPR 1 was significantly higher in fimH, ompT, afa/draBC, and univcnf genes than that of CRISPR 2 (p < 0.05). The most prevalent E. coli sequence types (EST) among 26 ESTs was EST 22 (45.1%), followed by EST 4 (23.2%), EST 16 (20.1%), EST 25 (19.5%), and EST 24 (18.3%). Interestingly, four genes, fimH, ompT, afa/draBC, and univcnf, had a significantly higher prevalence in both EST 4 and EST 22 (p < 0.05). Among the seven protospacers derived from CRISPR 1, protospacer 163 had the highest prevalence (20.4%), and it only existed in EST 4 and EST 22. This study suggests that the CRISPR sequence-typing approach can help to clarify and trace virulence potential, although the E. coli isolates were from normal bulk tank milk and not from mastitis.

Endogenous CRISPR-Cas Systems in Group I Clostridium botulinum and Clostridium sporogenes Do Not Directly Target the Botulinum Neurotoxin Gene Cluster

Most strains of proteolytic group I Clostridium botulinum (G1 C. botulinum) and some strains of Clostridium sporogenes possess genes encoding botulinum neurotoxin (BoNT), a potent neuroparalytic agent. Within G1 C. botulinum, conserved bont gene clusters of three major toxin serotypes (bont/A/B/F) can be found on conjugative plasmids and/or within chromosomal pathogenicity islands. CRISPR-Cas systems enable site-specific targeting of previously encountered mobile genetic elements (MGE) such as plasmids and bacteriophage through the creation of a spacer library complementary to protospacers within the MGEs. To examine whether endogenous CRISPR-Cas systems restrict the transfer of bont gene clusters across strains we conducted a bioinformatic analysis profiling endogenous CRISPR-Cas systems from 241 G1 C. botulinum and C. sporogenes strains. Approximately 6,200 CRISPR spacers were identified across the strains and Type I-B, III-A/B/D cas genes and CRISPR array features were identified in 83% of the strains. Mapping the predicted spacers against the masked strain and RefSeq plasmid dataset identified 56,000 spacer-protospacer matches. While spacers mapped heavily to targets within bont(+) plasmids, no protospacers were identified within the bont gene clusters. These results indicate the toxin is not a direct target of CRISPR-Cas but the plasmids predominantly responsible for its mobilization are. Finally, while the presence of a CRISPR-Cas system did not reliably indicate the presence or absence of a bont gene cluster, comparative genomics across strains indicates they often occupy the same hypervariable loci common to both species, potentially suggesting similar mechanisms are involved in the acquisition and curation of both genomic features.

CRISPR-Cas is associated with fewer antibiotic resistance genes in bacterial pathogens

The acquisition of antibiotic resistance (ABR) genes via horizontal gene transfer (HGT) is a key driver of the rise in multidrug resistance amongst bacterial pathogens. Bacterial defence systems per definition restrict the influx of foreign genetic material, and may therefore limit the acquisition of ABR. CRISPR-Cas adaptive immune systems are one of the most prevalent defences in bacteria, found in roughly half of bacterial genomes, but it has remained unclear if and how much they contribute to restricting the spread of ABR. We analysed approximately 40 000 whole genomes comprising the full RefSeq dataset for 11 species of clinically important genera of human pathogens, including Enterococcus, Staphylococcus, Acinetobacter and Pseudomonas. We modelled the association between CRISPR-Cas and indicators of HGT, and found that pathogens with a CRISPR-Cas system were less likely to carry ABR genes than those lacking this defence system. Analysis of the mobile genetic elements (MGEs) targeted by CRISPR-Cas supports a model where this host defence system blocks important vectors of ABR. These results suggest a potential 'immunocompromised' state for multidrug-resistant strains that may be exploited in tailored interventions that rely on MGEs, such as phages or phagemids, to treat infections caused by bacterial pathogens. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives

Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design.

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Coevolution between bacterial CRISPR-Cas systems and their bacteriophages

CRISPR-Cas systems provide bacteria and archaea with adaptive, heritable immunity against their viruses (bacteriophages and phages) and other parasitic genetic elements. CRISPR-Cas systems are highly diverse, and we are only beginning to understand their relative importance in phage defense. In this review, we will discuss when and why CRISPR-Cas immunity against phages evolves, and how this, in turn, selects for the evolution of immune evasion by phages. Finally, we will discuss our current understanding of if, and when, we observe coevolution between CRISPR-Cas systems and phages, and how this may be influenced by the mechanism of CRISPR-Cas immunity.

Identification of Natural CRISPR Systems and Targets in the Human Microbiome

Many bacteria resist invasive DNA by incorporating sequences into CRISPR loci, which enable sequence-specific degradation. CRISPR systems have been well studied from isolate genomes, but culture-independent metagenomics provide a new window into their diversity. We profiled CRISPR loci and cas genes in the body-wide human microbiome using 2,355 metagenomes, yielding functional and taxonomic profiles for 2.9 million spacers by aligning the spacer content to each sample's metagenome and corresponding gene families. Spacer and repeat profiles agree qualitatively with those from isolate genomes but expand their diversity by approximately 13-fold, with the highest spacer load present in the oral microbiome. The taxonomy of spacer sequences parallels that of their source community, with functional targets enriched for viral elements. When coupled with cas gene systems, CRISPR-Cas subtypes are highly site and taxon specific. Our analysis provides a comprehensive collection of natural CRISPR-cas loci and targets in the human microbiome.

Unravelling the Features of Success of VIM-Producing ST111 and ST235 Pseudomonas aeruginosa in a Greek Hospital

The objective of this study was to analyze the characteristics that contribute to the successful dissemination of VIM-producing Pseudomonas aeruginosa (P. aeruginosa), belonging to ST111 and ST235, in a Greek hospital. A total of 120 non-repetitive P. aeruginosa, which had meropenem minimal inhibitory concentrations (MICs) greater than 2 mg/L, were studied. VIM-encoding genes were amplified and sequenced within their integrons. Isolates were typed by multilocus sequence typing (MLST). Six VIM-producers, representative of different integron structures and sequence types (STs), were completely sequenced using Illumina platform. Sixty-one P. aeruginosa were confirmed to produce VIM-type carbapenemases. ST111 dominated (n = 34) among VIM-producers, while 15 VIM-producers belonged to ST235. The blaVIM-like genes were located in three integron types, including In59, In595 and In1760, which were integrated into P. aeruginosa chromosomes. Whole genome sequencing (WGS) data demonstrated that ST111 and ST235 MBL producers carried several resistance and virulence genes. Additionally, the presence of type I-C and type I-E clustered regularly interspaced short palindromic repeats (CRISPR)/Cas locus was observed in ST235 and ST395 isolates, respectively. In conclusion, our findings confirmed the clonal spread of ST111 P. aeruginosa, carrying the VIM-2-encoding integron In59, in the University Hospital of Larissa (UHL). In addition, they highlighted the important role of high-risk clones, ST111 and ST235, in the successful dissemination and establishment into hospital settings of clinically important pathogens carrying resistance determinants.

Association of CRISPR/Cas System with the Drug Resistance in Klebsiella pneumoniae

BACKGROUND

Klebsiella pneumoniae is a common opportunistic pathogen and its production of extended-spectrum β-lactamases (ESBL) and carbapenemases leads to drug resistance. Clustered regularly interspaced short palindromic repeats (CRISPRs) and their associated genes (Cas) are widespread in the genome of many bacteria and are a defense mechanism against foreign invaders such as plasmids and viruses.

PURPOSE

To investigate the prevalence of the CRISPR/Cas system in wild type strains of K. pneumoniae in the hospital and its association with drug resistance.

MATERIALS AND METHODS

A total of 136 strains were collected and characterized their susceptibility to antimicrobial agents. The prevalence of CRISPR/Cas system was detected by PCR and DNA sequencing was analyzed by CRISPRFinder. The statistical analysis of the results was performed by SPSS.

RESULTS

We found that 50/136 (37%) isolates produced ESBL and 30/136 (22%) isolates were resistant to carbapenems. These isolates were liable to be multidrug resistant against β-lactams, quinolones, and aminoglycosides. Among the carbapenem-resistant isolates, blaKPC was the main drug resistance-associated gene and different types of ESBL and AmpC genes were present. Resistance to β-lactams, quinolones, aminoglycosides, tetracyclines, and β-lactams/enzyme inhibitor were higher in absence of the CRISPR/Cas system. Eighteen spacers within the CRISPR arrays matched with the genomes of plasmids or phages, some of which carried drug resistance genes.

CONCLUSION

ESBL-producing and carbapenem-resistant K. pneumoniae are more likely to develop multidrug resistance and show an inverse correlation between drug resistance and CRISPR/Cas system. Absence of CRISPR/Cas modules allow for the acquisition of external drug resistance genes.

Genomic investigation of emerging zoonotic pathogen Shewanella xiamenensis

Objective

Shewanella xiamenensis is an emerging zoonotic pathogen commonly found in aquatic ecosystem. Clustered regularly interspaced short palindromic repeats (CRISPR) and (CRISPR)-associated gene systems act as adaptive immune system of prokaryotes. Recently, growing evidence suggested their role in bacterial virulence and resistance. Despite its medical importance, little is known about the genomic characteristics of S. xiamenensis.

Materials and Methods

Strain ZYW6 was isolated from Epinephelus awoara. We sequenced the 16S rRNA gene and blast against the GenBank bacterial database. Antibiotic susceptibility tests and interpretation were performed by automatic VITEK 2 system. We extracted the genomic DNA with QIAGEN Genomic-tip 100/G kit and QIAGEN Genomic DNA Buffer Set. Whole-genome shotgun sequencing was performed using the Illumina MiSeq sequencer. To identify the CRISPR-Cas System in the genome of S. xiamenensis ZYW6, the Integrated Microbial Genomes and Microbiomes and CRISPRFinder were used.

Results

We characterized the genome of a S. xiamenensis strain. The genome is 4,765,190 bp in length and encodes 4262 open-reading frames. Type I CRISPR-Cas system and serine biosynthesis genes were identified.

Conclusion

Our results demonstrate the genetic structure of CRISPR-Cas system, l-serine synthesis, and oxacillinase in S. xiamenensis. The report of antibiotics resistance genes in the study might indicate a possible reservoir of antimicrobial drug resistance determinants in food animal, resulting in potential infection source. The findings provide insights into the structure and composition of CRISPR-Cas system in S. xiamenensis and foundation for future biological validation.

Heterogeneity of Streptococcus anginosus ß-hemolysis in relation to CRISPR/Cas

Streptococcus anginosus is a commensal of the oral mucosa that can cause severe invasive infections. A considerable proportion of Streptococcus anginosus strains are ß-hemolytic due to the presence of an SLS-like gene cluster. However, the majority of strains do not display ß-hemolysis. To investigate ß-hemolysin heterogeneity in S. anginosus, we determined the presence of sag genes and correlated it with the presence of CRISPR/Cas genes in a collection of ß-hemolytic and non-ß-hemolytic strains. All of the ß-hemolytic strains carried the sag gene cluster. In contrast to other streptococci, clinical S. anginosus strains that do not display ß-hemolysis do not harbor sag genes. Phylogenetic analysis of the ß-hemolytic strains revealed that they belong to two previously defined clusters within S. anginosus. Correlation with CRISPR/Cas genes showed a significant difference for the presence of CRISPR/Cas in ß-hemolytic versus non-ß-hemolytic isolates. The presence of the CRISPR/Cas type IIA or type IIC locus is associated with the absence of sag genes; in 65% of the non-ß-hemolytic strains a CRISPR/Cas locus was found, while only 24% of ß-hemolytic strains carry CRISPR/Cas genes. Further analysis of the spacer content of the CRISPR systems revealed the presence of multiple self-targeting sequences directed against S. anginosus genes. These results support the hypothesis that horizontal gene transfer is involved in the acquisition of ß-hemolysin genes and that CRISPR/Cas may limit DNA uptake in S. anginosus.

CRISPR-Cas9 Genome Editing in Human Cell Lines with Donor Vector Made by Gibson Assembly

CRISPR Cas9 genome editing allows researchers to modify genes in a multitude of ways including to obtain deletions, epitope-tagged loci, and knock-in mutations. Within 6 years of its initial application, CRISPR-Cas9 genome editing has been widely employed, but disadvantages to this method, such as low modification efficiencies and off-target effects, need careful consideration. Obtaining custom donor vectors can also be expensive and time-consuming. This chapter details strategies to overcome barriers to CRISPR-Cas9 genome editing as well as recent developments in employing this technique.

Analysis of the structures of confirmed and questionable CRISPR loci in 325 Staphylococcus genomes

The CRISPR-Cas (clustered regular interspaced short palindromic repeats and CRISPR-associated proteins) system is a newly discovered immune defense system in the genome of prokaryotes, which can resist the invasion of foreign genetic elements, such as plasmids or phage. In this study, 154 strains of Staphylococcus published in the CRISPRDatabase and 171 strains included in NCBI were downloaded, the confirmed and questionable CRISPR loci of which were analyzed by bioinformatics methods, including their distribution, characteristics of the structure (including the direct repeats, spacers and cas genes), and the relationship between the presence of CRISPR and the mecA gene. Meanwhile, a comprehensive analysis of orphan CRISPR arrays was performed on this basis. A total of 196 confirmed and 1757 questionable CRISPR loci were found in 325 Staphylococcus genomes. Only 25 strains contained cas genes, which were classified into III-A (48.1%) and II-C (51.9%). The difference between the presence of the cas gene and the carrying rate of mecA was statistically significant, and they were negatively correlated. A total of 137 confirmed and 1755 questionable CRISPR loci were assumed to be false-CRISPR. The present study also analyzed the questionable CRISPR array for the first time while analyzing the confirmed CRISPR array in the Staphylococcal genome and screened the false-CRISPR elements in the orphan CRISPR array.

CRISPR-cas system: biological function in microbes and its use to treat antimicrobial resistant pathogens

The development of antibiotic resistance in bacteria is a major public health threat. Infection rates of resistant pathogens continue to rise against nearly all antimicrobials, which has led to development of different strategies to combat the antimicrobial resistance. In this review, we discuss how the newly popular CRISPR-cas system has been applied to combat antibiotic resistance in both extracellular and intracellular pathogens. We also review a recently developed method in which nano-size CRISPR complex was used without any phage to target the mecA gene. However, there is still challenge to practice these methods in field against emerging antimicrobial resistant pathogens.

A phylogenetic test of the role of CRISPR-Cas in limiting plasmid acquisition and prophage integration in bacteria

CRISPR-Cas is a prokaryotic defense system capable of protecting the cell from damaging foreign genetic elements. However, some genetic elements can be beneficial, which suggests the hypothesis that bacteria with CRISPR-Cas incur a cost of reduced intake of mutualistic plasmids and prophage. Here we present the first robust test of this hypothesis that controls for phylogenic and ecological biases in the distribution of CRISPR-Cas. We filtered the available genomic data (~7000 strains from ~2100 species) by first selecting all pairs of conspecific strains, one with and one without CRISPR-Cas (controlling ecological bias), and second by retaining only one such pair per bacterial family (controlling phylogenetic bias), resulting in pairs representing 38 bacterial families. Analysis of these pairs of bacterial strains showed that on average the CRISPR-Cas strain of each pair contained significantly fewer plasmids than its CRISPR-Cas negative partner (0.86 vs. 1.86). It also showed that the CRISPR-Cas positive strains had 31% fewer intact prophage (1.17 vs. 1.75), but the effect was highly variable and not significant. These results support the hypothesis that CRISPR-Cas reduces the rate of plasmid-mediated HGT and, given the abundant evidence of beneficial genes carried by plasmids, provide a clear example of a cost associated with the CRISPR-Cas system.

CRISPR-Cas in Streptococcus pyogenes

The discovery and characterization of the prokaryotic CRISPR-Cas immune system has led to a revolution in genome editing and engineering technologies. Despite the fact that most applications emerged after the discovery of the type II-A CRISPR-Cas9 system of Streptococcus pyogenes, its biological importance in this organism has received little attention. Here, we provide a comprehensive overview of the current knowledge about CRISPR-Cas systems from S. pyogenes. We discuss how the interplay between CRISPR-mediated immunity and horizontal gene transfer might have modeled the evolution of this pathogen. We review the current literature about the CRISPR-Cas systems present in S. pyogenes (types I-C and II-A), and describe their distinctive biochemical and functional features. Finally, we summarize the main biotechnological applications that have arisen from the discovery of the CRISPR-Cas9 system in S. pyogenes.

Insertion sequences in the CRISPR-Cas system regulate horizontal antimicrobial resistance gene transfer in Shigella strains

Multidrug-resistant (MDR) Shigella strains are an enormous threat to public health. Antimicrobial resistance genes are frequently located on plasmids, phages and integrons, which enter bacterial cells by horizontal gene transfer (HGT). CRISPR-Cas systems are adaptive prokaryotic immune systems in bacteria that confer resistance to foreign genetic material such as phages and other mobile genetic elements. However, this may come at a cost of inhibiting the acquisition of other beneficial genes through HGT. This study investigated how Shigella strains regulate the activity of the CRISPR-Cas system spontaneously when they require an exogenous gene necessary for survival. Insertion sequence (IS) elements were identified in cas genes, such as IS600 in cse2, ISSfl2 in cas6e and IS629 in cse1-cas3. The number of spacers in CRISPR-Cas arrays in strains containing an IS was less than that for strains with no IS. Interestingly, fewer spacers were also found in MDR Shigella isolates. Furthermore, an antimicrobial-resistant strain was constructed by electrotransformation of a resistance plasmid in order to detect changes in the CRISPR-Cas system. It was found that the cse2 gene had a new IS (IS600) in the antimicrobial-resistant strain. Bioinformatics analyses showed that the IS600 insertion hotspot was TGC-GGC in the cse2 gene, and the tertiary structure of the Cse2 protein was different with IS600. IS600 caused a five-order of magnitude decrease in relative expression of the cse2 gene. This study sheds mechanistic light on CRISPR-Cas-mediated HGT of antimicrobial resistance genes in Shigella spp. isolates.

Detection and variability analyses of CRISPR-like loci in the H. pylori genome

Helicobacter pylori is a human pathogenic bacterium with a high genomic plasticity. Although the functional CRISPR-Cas system has not been found in its genome, CRISPR-like loci have been recently identified. In this work, 53 genomes from different geographical areas are analyzed for the search and analysis of variability of this type of structure. We confirm the presence of a locus that was previously described in the VlpC gene in al lgenomes, and we characterize new CRISPR-like loci in other genomic locations. By studying the variability and gene location of these loci, the evolution and the possible roles of these sequences are discussed. Additionally, the usefulness of this type of sequences as a phylogenetic marker has been demonstrated, associating the different strains by geographical area.

In silico Analysis Suggests Common Appearance of scaRNAs in Type II Systems and Their Association With Bacterial Virulence

In addition to its well-established defense function, CRISPR/Cas can also exhibit crucial non-canonical activity through endogenous gene expression regulation, which was found to mainly affect bacterial virulence. These non-canonical functions depend on scaRNA, which is a small RNA encoded outside of CRISPR array, that is typically flanked by a transcription start site (TSS) and a terminator, and is in part complementary to another small CRISPR/Cas-associated RNA (tracrRNAs). Identification of scaRNAs is however largely complicated by the scarcity of RNA-Seq data across different bacteria, so that they were identified only in a relatively rare CRISPR/Cas subtype (IIB), and the possibility of finding them in other Type II systems is currently unclear. This study presents the first effort toward systematic detection of small CRISPR/Cas-associated regulatory RNAs, where obtained predictions can guide future experiments. The core of our approach is ab initio detection of small RNAs from bacterial genome, which is based on jointly predicting transcription signals - TSS and terminators - and homology to CRISPR array repeat. Particularly, we employ our improved approach for detecting bacterial TSS, since accurate TSS detection is the main limiting factor for accurate small RNA prediction. We also explore how our predictions match to available RNA-Seq data and analyze their conservation across related bacterial species. In Type IIB systems, our predictions are consistent with experimental data, and we systematically identify scaRNAs throughout this subtype. Furthermore, we identify scaRNA:tracrRNA pairs in a number of IIA/IIC systems, where the appearance of scaRNAs co-occurs with the strains being pathogenic. RNA-Seq and conservation analysis show that our method is well suited for predicting CRISPR/Cas-associated small RNAs. We also find possible existence of a modified mechanism of CRISPR-associated small RNA action, which, interestingly, closely resembles the setup employed in biotechnological applications. Overall, our findings indicate that scaRNA:tracrRNA pairs are present in all subtypes of Type II systems, and point to an underlying connection with bacterial virulence. In addition to formulating these hypotheses, careful manual curation that we performed, makes an important first step toward fully automated predictor of CRISPR/Cas-associated small RNAs, which will allow their large scale analysis across diverse bacterial genomes.

Comparative genomic analysis of Myroides odoratimimus isolates

Myroides odoratimimus is an important nosocomial pathogen. Management of M. odoratimimus infection is difficult owing to the multidrug resistance and the unknown pathogenesis mechanisms. Based on our previous genomic sequencing data of M. odoratimimus PR63039 (isolated from a patient with the urinary tract infection), in this study, we further performed comparative genomic analysis for 10 selected Myroides strains. Our results showed that these Myroides genome contexts were very similar and phylogenetically related. Various prophages were identified in the four clinical isolate genomes, which possibly contributed to the genome evolution among the Myroides strains. CRISPR elements were only detected in the two clinical (PR63039 and CCUG10230) isolates and two environmental (CCUG12700 and H1bi) strains. With more stringent cutoff parameters in CARD analysis, the four clinical M. odoratimimus contained roughly equal antibiotic resistance genes, indicating their similar antibiotic resistance profiles. The three clinical (CCUG10230, CCUG12901, CIP101113) and three environmental (CCUG12700, L41, H1bi) M. odoratimimus strains were speculated to carry the indistinguishable virulent factors (VFs), which may involve in the similar pathogenesis mechanism. Moreover, some VFs might confer to the high capacity of dissemination, attacking tissue cells and induction of autoimmune complications. Our results facilitate the research of antibiotic resistance and the development of therapeutic regimens for the M. odoratimimus infections.

Reconstituting the History of Cronobacter Evolution Driven by Differentiated CRISPR Activity

Cronobacter strains harboring the CRISPR-Cas system are important foodborne pathogens causing serious neonatal infections. However, the specific role of the CRISPR-Cas system in bacterial evolution remains relatively unexplored. In this study, we investigated the impact of the CRISPR-Cas system on Cronobacter evolution and obtained 137 new whole-genome Cronobacter sequences by next-generation sequencing technology. Among the strains examined (n = 240), 90.6% (193/213) of prevalent species Cronobacter sakazakii, Cronobactermalonaticus, and Cronobacterdublinensis strains had intact CRISPR-Cas systems. Two rare species, Cronobactercondimenti (n = 2) and Cronobacteruniversalis (n = 6), lacked and preserved the CRISPR-Cas system at a low frequency (1/6), respectively. These results suggest that the presence of one CRISPR-Cas system is important for a Cronobacter species to maintain genome homeostasis for survival. The Cronobacter ancestral strain is likely to have harbored both subtype I-E and I-F CRISPR-Cas systems; during the long evolutionary process, subtype I-E was retained while subtype I-F selectively degenerated in Cronobacter species and was even lost by the major Cronobacter pathovars. Moreover, significantly higher CRISPR activity was observed in the plant-associated species Cdublinensis than in the virulence-related species C. sakazakii and Cmalonaticus Similar spacers of CRISPR arrays were rarely found among species, suggesting intensive change through adaptive acquisition and loss. Differentiated CRISPR activity appears to be the product of environmental selective pressure and might contribute to the bidirectional divergence and speciation of CronobacterIMPORTANCE This study reports the evolutionary history of Cronobacter under the selective pressure of the CRISPR-Cas system. One CRISPR-Cas system in Cronobacter is important for maintaining genome homeostasis, whereas two types of systems may be redundant and not conducive to acquiring beneficial DNA for environmental adaptation and pathogenicity. Differentiated CRISPR activity has contributed to the bidirectional divergence and genetic diversity of Cronobacter This perspective makes a significant contribution to the literature by providing new insights into CRISPR-Cas systems in general, while further expanding the roles of CRISPR beyond conferring adaptive immunity and demonstrating a link to adaptation and species divergence in a genus. Moreover, our study provides new insights into the balance between genome homeostasis and the uptake of beneficial DNA related to CRISPR-based activity in the evolution of Cronobacter.

The physicist's guide to one of biotechnology's hottest new topics: CRISPR-Cas

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated proteins (Cas) constitute a multi-functional, constantly evolving immune system in bacteria and archaea cells. A heritable, molecular memory is generated of phage, plasmids, or other mobile genetic elements that attempt to attack the cell. This memory is used to recognize and interfere with subsequent invasions from the same genetic elements. This versatile prokaryotic tool has also been used to advance applications in biotechnology. Here we review a large body of CRISPR-Cas research to explore themes of evolution and selection, population dynamics, horizontal gene transfer, specific and cross-reactive interactions, cost and regulation, non-immunological CRISPR functions that boost host cell robustness, as well as applicable mechanisms for efficient and specific genetic engineering. We offer future directions that can be addressed by the physics community. Physical understanding of the CRISPR-Cas system will advance uses in biotechnology, such as developing cell lines and animal models, cell labeling and information storage, combatting antibiotic resistance, and human therapeutics.

The truncated Rv2820c of Mycobacterium tuberculosis Beijing family augments intracellular survival of M. smegmatis by altering cytokine profile and inhibiting NO generation

Genetic variations among genes of Mycobacterium tuberculosis may be associated with antigenic variation and immune evasion, which complicates the pathogenesis of M. tuberculosis. The hyper-virulent M. tuberculosis Beijing strains harbored several large sequence deletions, among which RD207 attributed to the deletion of CRISPR loci and several Cas genes. RD207 also gave rise to a truncated gene Rv2820c-Bj with 60% deletion in length at the 3'-end and a new 3'-end of five amino acid mutations. It has been reported that Rv2820c-Bj correlated with enhanced intracellular survival of M. smegmatis in macrophages when compared to its full-length counterpart Rv2820c in M. tuberculosis, however, the respective contribution of the truncation and the new 3'-end of Rv2820c-Bj to this enhancement was unclear. Here, by infecting THP-1 macrophages with Ms_Rv2820c-Bj, Ms_Rv2820c and MS_Rv2820c-Tr (expressing the truncated Rv2820c without five amino acid mutations at 3'-end), we found only Ms_Rv2820c-Bj was responsible for the enhancement of survival of M. smegmatis in macrophages. Furthermore, we detected that Ms_Rv2820c-Tr and Ms_Rv2820c-Bj induced similar cytokine profile and NO production after infection of macrophages, which was distinctly different from Ms_Rv2820c. However, Ms_Rv2820c-Bj evoked higher levels of interleukin-10 (IL-10) and lower levels of interleukin- 6 (IL-6), interleukin-1β (IL-1β) and interleukin-12 (IL-12) in infected THP-1 macrophages than Ms_Rv2820c-Tr. Accordingly, we concluded that the new 3'-end of Rv2820c-Bj was important to dampen host defense and enhance the intracellular survival of M. smegmatis.

Genome sequence and comparative analysis of Jiangella alba YIM 61503T isolated from a medicinal plant Maytenus austroyunnanensis

A draft genome sequence of Jiangella alba YIM 61503^T revealed a genome size of 7,664,864 bp arranged in 33 scaffolds. The genome was predicted to contain 7196 predicted genes, including 51 coding for RNA. Phylogenetic and comparative analyses of the draft genome of J. alba YIM 61503^T with the available genomes of other Jiangella species suggested a proximal similarity between strains J. alba YIM 61503^T and J. muralis DSM 45357^T, while indicating a high divergence between J. gansuensis YIM 002^T and other Jiangella species. The genome of J. alba YIM 61503^T also revealed genes involved in indole-3-acetic acid biosynthesis and an alkylresorcinols gene cluster. Further, detection of phosphotransferase genes in the genome of all Jiangella species indicated that they can uptake and phosphorylate sugars. The presences of TreX-Z, TreS and OtsA-OtsB genes in some of the Jiangella strains also indicated a possible mechanism for their tolerance of high salinity. Besides providing new insights into its genetic features, our results suggested that J. alba YIM 61503^T could be a potential strain for further genome mining studies. The release of this genome may, therefore, provide a better prospect for understanding "evolutionary taxonomy" about this genus in future.

Strategies for Editing Virulent Staphylococcal Phages Using CRISPR-Cas10

Staphylococci are prevalent skin-dwelling bacteria that are also leading causes of antibiotic-resistant infections. Viruses that infect and lyse these organisms (virulent staphylococcal phages) can be used as alternatives to conventional antibiotics and represent promising tools to eliminate or manipulate specific species in the microbiome. However, since over half their genes have unknown functions, virulent staphylococcal phages carry inherent risk to cause unknown downstream side effects. Further, their swift and destructive reproductive cycle make them intractable by current genetic engineering techniques. CRISPR-Cas10 is an elaborate prokaryotic immune system that employs small RNAs and a multisubunit protein complex to detect and destroy phages and other foreign nucleic acids. Some staphylococci naturally possess CRISPR-Cas10 systems, thus providing an attractive tool already installed in the host chromosome to harness for phage genome engineering. However, the efficiency of CRISPR-Cas10 immunity against virulent staphylococcal phages and corresponding utility as a tool to facilitate their genome editing has not been explored. Here, we show that the CRISPR-Cas10 system native to Staphylococcus epidermidis exhibits robust immunity against diverse virulent staphylococcal phages. On the basis of this activity, a general two-step approach was developed to edit these phages that relies upon homologous recombination machinery encoded in the host. Variations of this approach to edit toxic phage genes and access phages that infect CRISPR-less staphylococci are also presented. This versatile set of genetic tools enables the systematic study of phage genes of unknown functions and the design of genetically defined phage-based antimicrobials that can eliminate or manipulate specific Staphylococcus species.

CRISPR-Cas: biology, mechanisms and relevance

Prokaryotes have evolved several defence mechanisms to protect themselves from viral predators. Clustered regularly interspaced short palindromic repeats (CRISPR) and their associated proteins (Cas) display a prokaryotic adaptive immune system that memorizes previous infections by integrating short sequences of invading genomes—termed spacers—into the CRISPR locus. The spacers interspaced with repeats are expressed as small guide CRISPR RNAs (crRNAs) that are employed by Cas proteins to target invaders sequence-specifically upon a reoccurring infection. The ability of the minimal CRISPR-Cas9 system to target DNA sequences using programmable RNAs has opened new avenues in genome editing in a broad range of cells and organisms with high potential in therapeutical applications. While numerous scientific studies have shed light on the biochemical processes behind CRISPR-Cas systems, several aspects of the immunity steps, however, still lack sufficient understanding. This review summarizes major discoveries in the CRISPR-Cas field, discusses the role of CRISPR-Cas in prokaryotic immunity and other physiological properties, and describes applications of the system as a DNA editing technology and antimicrobial agent.

This article is part of the themed issue ‘The new bacteriology’.

Gaussian Accelerated Molecular Dynamics: Theory, Implementation, and Applications

A novel Gaussian Accelerated Molecular Dynamics (GaMD) method has been developed for simultaneous unconstrained enhanced sampling and free energy calculation of biomolecules. Without the need to set predefined reaction coordinates, GaMD enables unconstrained enhanced sampling of the biomolecules. Furthermore, by constructing a boost potential that follows a Gaussian distribution, accurate reweighting of GaMD simulations is achieved via cumulant expansion to the second order. The free energy profiles obtained from GaMD simulations allow us to identify distinct low energy states of the biomolecules and characterize biomolecular structural dynamics quantitatively. In this chapter, we present the theory of GaMD, its implementation in the widely used molecular dynamics software packages (AMBER and NAMD), and applications to the alanine dipeptide biomolecular model system, protein folding, biomolecular large-scale conformational transitions and biomolecular recognition.

Comparative Analysis of the Flavobacterium columnare Genomovar I and II Genomes

Columnaris disease caused by Gram-negative rod Flavobacterium columnare is one of the most common diseases of catfish. F. columnare is also a common problem in other cultured fish species worldwide. F. columnare has three major genomovars; we have sequenced a representative strain from genomovar I (ATCC 49512, which is avirulent in catfish) and genomovar II (94-081, which is highly pathogenic in catfish). Here, we present a comparative analysis of the two genomes. Interestingly, F. columnare ATCC 49512 and 94-081 meet criteria to be considered different species based on the Average Nucleotide Identity (90.71% similar) and DNA–DNA Hybridization (42.6% similar). Genome alignment indicated the two genomes have a large number of rearrangements. However, function-based comparative genomics analysis indicated that the two strains have similar functional capabilities with 2,263 conserved orthologous clusters; strain ATCC 49512 has 290 unique orthologous clusters while strain 94-081 has 391. Both strains carry type I secretion system, type VI secretion system, and type IX secretion system. The two genomes also have similar CRISPR capacities. The F. columnare strain ATCC 49512 genome contains a higher number of insertion sequence families and phage regions, while the F. columnare strain 94-081 genome has more genomic islands and more regulatory gene capacity. Transposon mutagenesis using Tn4351 in pathogenic strain 94-081 yielded six mutants, and experimental infections of fish showed hemolysin and glycine cleavage protein mutants had 15 and 10% mortalities, respectively, while the wild-type strain caused 100% mortalities. Our comparative and mutational analysis yielded important information on classification of genomovars I and II F. columnare as well as potential virulence genes in F. columnare strain 94-081.

Genome editing approaches: manipulating of lovastatin and taxol synthesis of filamentous fungi by CRISPR/Cas9 system

Filamentous fungi are prolific repertoire of structurally diverse secondary metabolites of remarkable biological activities such as lovastatin and paclitaxel that have been approved by FDA as drugs for hypercholesterolemia and cancer treatment. The clusters of genes encoding lovastatin and paclitaxel are cryptic at standard laboratory cultural conditions (Kennedy et al. Science 284:1368-1372, 1999; Bergmann et al. Nature Chem Biol 3:213-217, 2007). The expression of these genes might be triggered in response to nutritional and physical conditions; nevertheless, the overall yield of these metabolites does not match the global need. Consequently, overexpression of the downstream limiting enzymes and/or blocking the competing metabolic pathways of these metabolites could be the most successful technologies to enhance their yield. This is the first review summarizing the different strategies implemented for fungal genome editing, molecular regulatory mechanisms, and prospective of clustered regulatory interspaced short palindromic repeat/Cas9 system in metabolic engineering of fungi to improve their yield of lovastatin and taxol to industrial scale. Thus, elucidating the putative metabolic pathways in fungi for overproduction of lovastatin and taxol was the ultimate objective of this review.

CRISPR-Cas Defense System and Potential Prophages in Cyanobacteria Associated with the Coral Black Band Disease

Understanding how pathogens maintain their virulence is critical to developing tools to mitigate disease in animal populations. We sequenced and assembled the first draft genome of Roseofilum reptotaenium AO1, the dominant cyanobacterium underlying pathogenicity of the virulent coral black band disease (BBD), and analyzed parts of the BBD-associated Geitlerinema sp. BBD_1991 genome in silico. Both cyanobacteria are equipped with an adaptive, heritable clustered regularly interspaced short palindromic repeats (CRISPR)-Cas defense system type I-D and have potential virulence genes located within several prophage regions. The defense system helps to prevent infection by viruses and mobile genetic elements via identification of short fingerprints of the intruding DNA, which are stored as templates in the bacterial genome, in so-called "CRISPRs." Analysis of CRISPR target sequences (protospacers) revealed an unusually high number of self-targeting spacers in R. reptotaenium AO1 and extraordinary long CRIPSR arrays of up to 260 spacers in Geitlerinema sp. BBD_1991. The self-targeting spacers are unlikely to be a form of autoimmunity; instead these target an incomplete lysogenic bacteriophage. Lysogenic virus induction experiments with mitomycin C and UV light did not reveal an actively replicating virus population in R. reptotaenium AO1 cultures, suggesting that phage functionality is compromised or excision could be blocked by the CRISPR-Cas system. Potential prophages were identified in three regions of R. reptotaenium AO1 and five regions of Geitlerinema sp. BBD_1991, containing putative BBD relevant virulence genes, such as an NAD-dependent epimerase/dehydratase (a homolog in terms of functionality to the third and fourth most expressed gene in BBD), lysozyme/metalloendopeptidases and other lipopolysaccharide modification genes. To date, viruses have not been considered to be a component of the BBD consortium or a contributor to the virulence of R. reptotaenium AO1 and Geitlerinema sp. BBD_1991. We suggest that the presence of virulence genes in potential prophage regions, and the CRISPR-Cas defense systems are evidence of an arms race between the respective cyanobacteria and their bacteriophage predators. The presence of such a defense system likely reduces the number of successful bacteriophage infections and mortality in the cyanobacteria, facilitating the progress of BBD.

Survival and Evolution of CRISPR-Cas System in Prokaryotes and Its Applications

Prokaryotes have developed numerous innate immune mechanisms in order to fend off bacteriophage or plasmid attack. One of these immune systems is clustered regularly interspaced short palindromic repeats (CRISPR). CRISPR-associated proteins play a key role in survival of prokaryotes against invaders, as these systems cleave DNA of foreign genetic elements. Beyond providing immunity, these systems have significant impact in altering the bacterial physiology in term of its virulence and pathogenicity, as well as evolution. Also, due to their diverse nature of functionality, cas9 endoribonuclease can be easily reprogrammed with the help of guide RNAs, showing unprecedented potential and significance for gene editing in treating genetic diseases. Here, we also discuss the use of NgAgo–gDNA system in genome editing of human cells.