Comparative Genomic Analysis of Mycobacteriaceae Reveals Horizontal Gene Transfer-Mediated Evolution of the CRISPR-Cas System in the Mycobacterium tuberculosis Complex

Pangenome-wide association study reveals the selective absence of CRISPR genes (Rv2816c-19c) in drug-resistant Mycobacterium tuberculosis

The presence of intermittently dispersed insertion sequences and transposases in the Mycobacterium tuberculosis (Mtb) genome makes intra-genome recombination events inevitable. Understanding their effect on the gene repertoires (GR), which may contribute to the development of drug-resistant Mtb, is critical. In this study, publicly available WGS data of clinical Mtb isolates (endemic region n = 2,601; non-endemic region n = 1,130) were de novo assembled, filtered, scaffolded into assemblies, and functionally annotated. Out of 2,601 Mtb WGS data sets from endemic regions, 2,184 (drug resistant/sensitive: 1,386/798) qualified as high quality. We identified 3,784 core genes, 123 softcore genes, 224 shell genes, and 762 cloud genes in the pangenome of Mtb clinical isolates from endemic regions. Sets of 33 and 39 genes showed positive and negative associations (P < 0.01) with drug resistance status, respectively. Gene ontology clustering showed compromised immunity to phages and impaired DNA repair in drug-resistant Mtb clinical isolates compared to the sensitive ones. Multidrug efflux pump repressor genes (Rv3830c and Rv3855c) and CRISPR genes (Rv2816c-19c) were absent in the drug-resistant Mtb. A separate WGS data analysis of drug-resistant Mtb clinical isolates from the Netherlands (n = 1130) also showed the absence of CRISPR genes (Rv2816c-17c). This study highlights the role of CRISPR genes in drug resistance development in Mtb clinical isolates and helps in understanding its evolutionary trajectory and as useful targets for diagnostics development.IMPORTANCEThe results from the present Pan-GWAS study comparing gene sets in drug-resistant and drug-sensitive Mtb clinical isolates revealed intricate presence-absence patterns of genes encoding DNA-binding proteins having gene regulatory as well as DNA modification and DNA repair roles. Apart from the genes with known functions, some uncharacterized and hypothetical genes that seem to have a potential role in drug resistance development in Mtb were identified. We have been able to extrapolate many findings of the present study with the existing literature on the molecular aspects of drug-resistant Mtb, further strengthening the relevance of the results presented in this study.

Characteristics and immune functions of the endogenous CRISPR-Cas systems in myxobacteria

The clustered regularly interspaced short palindromic repeats and their associated proteins (CRISPR-Cas) system widely occurs in prokaryotic organisms to recognize and destruct genetic invaders. Systematic collation and characterization of endogenous CRISPR-Cas systems are conducive to our understanding and potential utilization of this natural genetic machinery. In this study, we screened 39 complete and 692 incomplete genomes of myxobacteria using a combined strategy to dispose of the abridged genome information and revealed at least 19 CRISPR-Cas subtypes, which were distributed with a taxonomic difference and often lost stochastically in intraspecies strains. The cas genes in each subtype were evolutionarily clustered but deeply separated, while most of the CRISPRs were divided into four types based on the motif characteristics of repeat sequences. The spacers recorded in myxobacterial CRISPRs were in high G+C content, matching lots of phages, tiny amounts of plasmids, and, surprisingly, massive organismic genomes. We experimentally demonstrated the immune and self-target immune activities of three endogenous systems in Myxococcus xanthus DK1622 against artificial genetic invaders and revealed the microhomology-mediated end-joining mechanism for the immunity-induced DNA repair but not homology-directed repair. The panoramic view and immune activities imply potential omnipotent immune functions and applications of the endogenous CRISPR-Cas machinery.

IMPORTANCE

Serving as an adaptive immune system, clustered regularly interspaced short palindromic repeats and their associated proteins (CRISPR-Cas) empower prokaryotes to fend off the intrusion of external genetic materials. Myxobacteria are a collective of swarming Gram-stain-negative predatory bacteria distinguished by intricate multicellular social behavior. An in-depth analysis of their intrinsic CRISPR-Cas systems is beneficial for our understanding of the survival strategies employed by host cells within their environmental niches. Moreover, the experimental findings presented in this study not only suggest the robust immune functions of CRISPR-Cas in myxobacteria but also their potential applications.

Evolution and emergence of Mycobacterium tuberculosis

Tuberculosis (TB) remains one of the deadliest infectious diseases in human history, prevailing even in the 21st century. The causative agents of TB are represented by a group of closely related bacteria belonging to the Mycobacterium tuberculosis complex (MTBC), which can be subdivided into several lineages of human- and animal-adapted strains, thought to have shared a last common ancestor emerged by clonal expansion from a pool of recombinogenic Mycobacterium canettii-like tubercle bacilli. A better understanding of how MTBC populations evolved from less virulent mycobacteria may allow for discovering improved TB control strategies and future epidemiologic trends. In this review, we highlight new insights into the evolution of mycobacteria at the genus level, describing different milestones in the evolution of mycobacteria, with a focus on the genomic events that have likely enabled the emergence and the dominance of the MTBC. We also review the recent literature describing the various MTBC lineages and highlight their particularities and differences with a focus on host preferences and geographic distribution. Finally, we discuss on putative mechanisms driving the evolution of tubercle bacilli and mycobacteria in general, by taking the mycobacteria-specific distributive conjugal transfer as an example.

Identifying ARG-carrying bacteriophages in a lake replenished by reclaimed water using deep learning techniques

As important mobile genetic elements, phages support the spread of antibiotic resistance genes (ARGs). Previous analyses of metaviromes or metagenome-assembled genomes (MAGs) failed to assess the extent of ARGs transferred by phages, particularly in the generation of antibiotic pathogens. Therefore, we have developed a bioinformatic pipeline that utilizes deep learning techniques to identify ARG-carrying phages and predict their hosts, with a special focus on pathogens. Using this method, we discovered that the predominant types of ARGs carried by temperate phages in a typical landscape lake, which is fully replenished by reclaimed water, were related to multidrug resistance and β-lactam antibiotics. MAGs containing virulent factors (VFs) were predicted to serve as hosts for these ARG-carrying phages, which suggests that the phages may have the potential to transfer ARGs. In silico analysis showed a significant positive correlation between temperate phages and host pathogens (R = 0.503, p < 0.001), which was later confirmed by qPCR. Interestingly, these MAGs were found to be more abundant than those containing both ARGs and VFs, especially in December and March. Seasonal variations were observed in the abundance of phages harboring ARGs (from 5.62 % to 21.02 %) and chromosomes harboring ARGs (from 18.01 % to 30.94 %). In contrast, the abundance of plasmids harboring ARGs remained unchanged. In summary, this study leverages deep learning to analyze phage-transferred ARGs and demonstrates an alternative method to track the production of potential antibiotic-resistant pathogens by metagenomics that can be extended to microbiological risk assessment.

Horizontal Gene Transfer and Drug Resistance Involving Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) acquires drug resistance at a rate comparable to that of bacterial pathogens that replicate much faster and have a higher mutation rate. One explanation for this rapid acquisition of drug resistance in Mtb is that drug resistance may evolve in other fast-replicating mycobacteria and then be transferred to Mtb through horizontal gene transfer (HGT). This paper aims to address three questions. First, does HGT occur between Mtb and other mycobacterial species? Second, what genes after HGT tend to survive in the recipient genome? Third, does HGT contribute to antibiotic resistance in Mtb? I present a conceptual framework for detecting HGT and analyze 39 ribosomal protein genes, 23S and 16S ribosomal RNA genes, as well as several genes targeted by antibiotics against Mtb, from 43 genomes representing all major groups within Mycobacterium. I also included mgtC and the insertion sequence IS6110 that were previously reported to be involved in HGT. The insertion sequence IS6110 shows clearly that the Mtb complex participates in HGT. However, the horizontal transferability of genes depends on gene function, as was previously hypothesized. HGT is not observed in functionally important genes such as ribosomal protein genes, rRNA genes, and other genes chosen as drug targets. This pattern can be explained by differential selection against functionally important and unimportant genes after HGT. Functionally unimportant genes such as IS6110 are not strongly selected against, so HGT events involving such genes are visible. For functionally important genes, a horizontally transferred diverged homologue from a different species may not work as well as the native counterpart, so the HGT event involving such genes is strongly selected against and eliminated, rendering them invisible to us. In short, while HGT involving the Mtb complex occurs, antibiotic resistance in the Mtb complex arose from mutations in those drug-targeted genes within the Mtb complex and was not gained through HGT.

Cold Cas: reevaluating the occurrence of CRISPR/Cas systems in Mycobacteriaceae

Bacterial CRISPR/Cas systems target foreign genetic elements such as phages and regulate gene expression by some pathogens, even in the host. The system is a marker for evolutionary history and has been used for inferences in Mycobacterium tuberculosis for 30 years. However, knowledge about mycobacterial CRISPR/Cas systems remains limited. It is believed that Type III-A Cas systems are exclusive to Mycobacterium canettii and the M. tuberculosis complex (MTBC) of organisms and that very few of the >200 diverse species of non-tuberculous mycobacteria (NTM) possess any CRISPR/Cas system. This study sought unreported CRISPR/Cas loci across NTM to better understand mycobacterial evolution, particularly in species phylogenetically near the MTBC. An analysis of available mycobacterial genomes revealed that Cas systems are widespread across Mycobacteriaceae and that some species contain multiple types. The phylogeny of Cas loci shows scattered presence in many NTM, with variation even within species, suggesting gains/losses of these loci occur frequently. Cas Type III-A systems were identified in pathogenic Mycobacterium heckeshornense and the geological environmental isolate Mycobacterium SM1. In summary, mycobacterial CRISPR/Cas systems are numerous, Type III-A systems are unreliable as markers for MTBC evolution, and mycobacterial horizontal gene transfer appears to be a frequent source of genetic variation.

Comparative Analysis of CRISPR-Cas Systems in Pseudomonas Genomes

Pseudomonas is a bacterial genus with some saprophytic species from land and others associated with opportunistic infections in humans and animals. Factors such as pathogenicity or metabolic aspects have been related to CRISPR-Cas, and in silico studies into it have focused more on the clinical and non-environmental setting. This work aimed to perform an in silico analysis of the CRISPR-Cas systems present in Pseudomonas genomes. It analyzed 275 complete genomic sequences of Pseudomonas taken from the NCBI database. CRISPR loci were obtained from CRISPRdb. The genes associated with CRISPR (cas) and CAS proteins, and the origin and diversity of spacer sequences, were identified and compared by BLAST. The presence of self-targeting sequences, PAMs, and the conservation of DRs were visualized using WebLogo 3.6. The CRISPR-like RNA secondary structure prediction was analyzed using RNAFold and MFold. CRISPR structures were identified in 19.6% of Pseudomonas species. In all, 113 typical CRISPR arrays with 18 putative cas were found, as were 2050 spacers, of which 52% showed homology to bacteriophages, 26% to chromosomes, and 22% to plasmids. No potential self-targeting was detected within the CRISPR array. All the found DRs can form thermodynamically stable secondary RNA structures. The comparison of the CRISPR/Cas system can help understand the environmental adaptability of each evolutionary lineage of clinically and environmentally relevant species, providing data support for bacterial typing, traceability, analysis, and exploration of unconventional CRISPR.

Characteristics of subtype III-A CRISPR-Cas system in Mycobacterium tuberculosis: An overview

CRISPR-Cas systems are the only RNA- guided adaptive immunity pathways that trigger the detection and destruction of invasive phages and plasmids in bacteria and archaea. Due to its prevalence and mystery, the Class 1 CRISPR-Cas system has lately been the subject of several studies. This review highlights the specificity of CRISPR-Cas system III-A in Mycobacterium tuberculosis, the tuberculosis-causing pathogen, for over twenty years. We discuss the difference between the several subtypes of Type III and their defence mechanisms. The anti-CRISPRs (Acrs) recently described, the critical role of Reverse transcriptase (RT) and housekeeping nuclease for type III CRISPR-Cas systems, and the use of this cutting-edge technology, its impact on the search for novel anti-tuberculosis drugs.

Revealing bacteria-phage interactions in human microbiome through the CRISPR-Cas immune systems

The human gut microbiome is composed of a diverse consortium of microorganisms. Relatively little is known about the diversity of the bacteriophage population and their interactions with microbial organisms in the human microbiome. Due to the persistent rivalry between microbial organisms (hosts) and phages (invaders), genetic traces of phages are found in the hosts' CRISPR-Cas adaptive immune system. Mobile genetic elements (MGEs) found in bacteria include genetic material from phage and plasmids, often resultant from invasion events. We developed a computational pipeline (BacMGEnet), which can be used for inference and exploratory analysis of putative interactions between microbial organisms and MGEs (phages and plasmids) and their interaction network. Given a collection of genomes as the input, BacMGEnet utilizes computational tools we have previously developed to characterize CRISPR-Cas systems in the genomes, which are then used to identify putative invaders from publicly available collections of phage/prophage sequences. In addition, BacMGEnet uses a greedy algorithm to summarize identified putative interactions to produce a bacteria-MGE network in a standard network format. Inferred networks can be utilized to assist further examination of the putative interactions and for discovery of interaction patterns. Here we apply the BacMGEnet pipeline to a few collections of genomic/metagenomic datasets to demonstrate its utilities. BacMGEnet revealed a complex interaction network of the Phocaeicola vulgatus pangenome with its phage invaders, and the modularity analysis of the resulted network suggested differential activities of the different P. vulgatus' CRISPR-Cas systems (Type I-C and Type II-C) against some phages. Analysis of the phage-bacteria interaction network of human gut microbiome revealed a mixture of phages with a broad host range (resulting in large modules with many bacteria and phages), and phages with narrow host range. We also showed that BacMGEnet can be used to infer phages that invade bacteria and their interactions in wound microbiome. We anticipate that BacMGEnet will become an important tool for studying the interactions between bacteria and their invaders for microbiome research.

Evaluating the efficacy of stool sample on Xpert MTB/RIF Ultra and its comparison with other sample types by meta-analysis for TB diagnostics

Precise and timely detection of tuberculosis (TB) is crucial to reduce transmission. This study aims to assess the accuracy of Xpert MTB/RIF Ultra on stool samples and systematically review the performance of Xpert MTB/RIF Ultra with different sample types by meta-analysis. Stool samples of smear-negative pulmonary TB (PTB), cervical lymph node TB, and abdominal TB patients were tested on the Xpert MTB/RIF Ultra system. Meta-analysis was performed on a set of 44 studies. Data were grouped by sample type, and the pooled sensitivity and specificity of Xpert MTB/RIF Ultra were calculated. The sensitivity of Xpert MTB/RIF Ultra with stool samples was 100% for smear-negative PTB, 27.27% for cervical lymph node TB, and 50% for abdominal TB patients, with 100% specificity for all included TB groups. The summary estimate for all PTB samples showed 84.2% sensitivity and 94.5% specificity, and EPTB samples showed 88.6% sensitivity and 96.4% specificity. Among all sample types included in our meta-analysis, urine showed the best performance for EPTB diagnosis. This pilot study supports the use of stool as an alternative non-invasive sample on Xpert MTB/RIF Ultra for rapid testing, suitable for both PTB and EPTB diagnosis.

The Involvement of Mycobacterium Type III-A CRISPR-Cas System in Oxidative Stress

Type I and type II CRISPR-Cas systems are employed to evade host immunity by targeting interference of bacteria’s own genes. Although Mycobacterium tuberculosis (M. tuberculosis), the causative agent of tuberculosis, possesses integrated type III-A CRISPR-Cas system, its role in mycobacteria remains obscure. Here, we observed that seven cas genes (csm2∼5, cas10, cas6) were upregulated in Mycobacterium bovis BCG under oxidative stress treatment, indicating the role of type III-A CRISPR-Cas system in oxidative stress. To explore the functional role of type III-A CRISPR-Cas system, TCC (Type III-A CRISPR-Cas system, including cas6, cas10, and csm2-6) mutant was generated. Deletion of TCC results in increased sensitivity in response to hydrogen peroxide and reduced cell envelope integrity. Analysis of RNA-seq dataset revealed that TCC impacted on the oxidation-reduction process and the composition of cell wall which is essential for mycobacterial envelop integrity. Moreover, disrupting TCC led to poor intracellular survival in vivo and in vitro. Finally, we showed for the first time that TCC contributed to the regulation of regulatory T cell population, supporting a role of TCC in modulating host immunity. Our finding reveals the important role of TCC in cell envelop homeostasis. Our work also highlights type III-A CRISPR-Cas system as an important factor for intracellular survival and host immunoregulation in mycobacteria, thus may be a potential target for therapy.

Complete Genome Sequence of Weissella confusa LM1 and Comparative Genomic Analysis

The genus Weissella is attracting an increasing amount of attention because of its multiple functions and probiotic potential. In particular, the species Weissella confusa is known to have great potential in industrial applications and exhibits numerous biological functions. However, the knowledge on this bacterium in insects is not investigated. Here, we isolated and identified W. confusa as the dominant lactic acid bacteria in the gut of the migratory locust. We named this strain W. confusa LM1, which is the first genome of an insect-derived W. confusa strain with one complete chromosome and one complete plasmid. Among all W. confusa strains, W. confusa LM1 had the largest genome. Its genome was the closest to that of W. confusa 1001271B_151109_G12, a strain from human feces. Our results provided accurate evolutionary relationships of known Weissella species and W. confusa strains. Based on genomic analysis, the pan-genome of W. confusa is in an open state. Most strains of W. confusa had the unique genes, indicating that these strains can adapt to different ecological niches and organisms. However, the variation of strain-specific genes did represent significant correlations with their hosts and ecological niches. These strains were predicted to have low potential to produce secondary metabolites. Furthermore, no antibiotic resistance genes were identified. At the same time, virulence factors associated with toxin production and secretion system were not found, indicating that W. confusa strains were not sufficient to perform virulence. Our study facilitated the discovery of the functions of W. confusa LM1 in locust biology and their potential application to locust management.