Spacer-Based Macroarrays for CRISPR Genotyping

Bioarchaeological and molecular evidence of tuberculosis in human skeletal remains from 18th-19th century orthodox cemeteries in Irkutsk, Eastern Siberia

In this study, we tested the skeletal human remains from the 18th - early 19th century Orthodox cemeteries in Irkutsk, Eastern Siberia, for tuberculosis-associated morphological alterations and Mycobacterium tuberculosis DNA. The morphologically studied bone collection included 591 individuals of mainly Caucasian origin. The molecular methods (IS6110-PCR and spoligotyping) suggested that at least four individuals (out of 15 TB-suspected, DNA-tested) were positive for the presence of M. tuberculosis DNA. All of them were males (3 maturus, 1 maturus senilis). Two of them date back to the second and third quarters of the 18th century, another to the last quarter of the 18th century, and the last one to the second half of the 19th century. The combined molecular analysis cautiously suggested presence of different strains and at least some of them represented not the currently predominant in Siberia Beijing genotype (M. tuberculosis East-Asian lineage) but strains of European origin. In conclusion, this study presented bioarchaeological and molecular evidence of tuberculosis in human skeletal remains from 18th-19th century Orthodox cemeteries in Irkutsk, Eastern Siberia. The samples are not M. bovis and represent human M. tuberculosis sensu stricto. Their precise phylogenetic identity is elusive but evokes the European/Russian origin of at least some isolates.

Insight into pathogenomics and phylogeography of hypervirulent and highly-lethal Mycobacterium tuberculosis strain cluster

BACKGROUND

. The Mycobacterium tuberculosis Beijing genotype is globally spread lineage with important medical properties that however vary among its subtypes. M. tuberculosis Beijing 14717-15-cluster was recently discovered as both multidrug-resistant, hypervirulent, and highly-lethal strain circulating in the Far Eastern region of Russia. Here, we aimed to analyze its pathogenomic features and phylogeographic pattern.

RESULTS

. The study collection included M. tuberculosis DNA collected between 1996 and 2020 in different world regions. The bacterial DNA was subjected to genotyping and whole genome sequencing followed by bioinformatics and phylogenetic analysis. The PCR-based assay to detect specific SNPs of the Beijing 14717-15-cluster was developed and used for its screening in the global collections. Phylogenomic and phylogeographic analysis confirmed endemic prevalence of the Beijing 14717-15-cluster in the Asian part of Russia, and distant common ancestor with isolates from Korea (> 115 SNPs). The Beijing 14717-15-cluster isolates had two common resistance mutations RpsL Lys88Arg and KatG Ser315Thr and belonged to spoligotype SIT269. The Russian isolates of this cluster were from the Asian Russia while 4 isolates were from the Netherlands and Spain. The cluster-specific SNPs that significantly affect the protein function were identified in silico in genes within different categories (lipid metabolism, regulatory proteins, intermediary metabolism and respiration, PE/PPE, cell wall and cell processes).

CONCLUSIONS

. We developed a simple method based on real-time PCR to detect clinically significant MDR and hypervirulent Beijing 14717-15-cluster. Most of the identified cluster-specific mutations were previously unreported and could potentially be associated with increased pathogenic properties of this hypervirulent M. tuberculosis strain. Further experimental study to assess the pathobiological role of these mutations is warranted.

CrisprVi: a software for visualizing and analyzing CRISPR sequences of prokaryotes

Background

Clustered regularly interspaced short palindromic repeats (CRISPR) and their spacers are important components of prokaryotic CRISPR-Cas systems. In order to analyze the CRISPR loci of multiple genomes more intuitively and comparatively, here we propose a visualization analysis tool named CrisprVi.

Results

CrisprVi is a Python package consisting of a graphic user interface (GUI) for visualization, a module for commands parsing and data transmission, local SQLite and BLAST databases for data storage and a functions layer for data processing. CrisprVi can not only visually present information of CRISPR direct repeats (DRs) and spacers, such as their orders on the genome, IDs, start and end coordinates, but also provide interactive operation for users to display, label and align the CRISPR sequences, which help researchers investigate the locations, orders and components of the CRISPR sequences in a global view. In comparison to other CRISPR visualization tools such as CRISPRviz and CRISPRStudio, CrisprVi not only improves the interactivity and effects of the visualization, but also provides basic statistics of the CRISPR sequences, and the consensus sequences of DRs/spacers across the input strains can be inspected from a clustering heatmap based on the BLAST results of the CRISPR sequences hitting against the genomes.

Conclusions

CrisprVi is a convenient tool for visualizing and analyzing the CRISPR sequences and it would be helpful for users to inspect novel CRISPR-Cas systems of prokaryotes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04716-9.

Practical approach to detection and surveillance of emerging highly resistant Mycobacterium tuberculosis Beijing 1071-32-cluster

Ancient sublineage of the Mycobacterium tuberculosis Beijing genotype is endemic and prevalent in East Asia and rare in other world regions. While these strains are mainly drug susceptible, we recently identified a novel clonal group Beijing 1071-32 within this sublineage emerging in Siberia, Russia and present in other Russian regions. This cluster included only multi/extensive drug resistant (MDR/XDR) isolates. Based on the phylogenetic analysis of the available WGS data, we identified three synonymous SNPs in the genes Rv0144, Rv0373c, and Rv0334 that were specific for the Beijing 1071-32-cluster and developed a real-time PCR assay for their detection. Analysis of the 2375 genetically diverse M. tuberculosis isolates collected between 1996 and 2020 in different locations (European and Asian parts of Russia, former Soviet Union countries, Albania, Greece, China, Vietnam, Japan and Brazil), confirmed 100% specificity and sensitivity of this real-time PCR assay. Moreover, the epidemiological importance of this strain and the newly developed screening assay is further stressed by the fact that all identified Beijing 1071-32 isolates were found to exhibit MDR genotypic profiles with concomitant resistance to additional first-line drugs due to a characteristic signature of six mutations in rpoB450, rpoC485, katG315, katG335, rpsL43 and embB497. In conclusion, this study provides a set of three concordant SNPs for the detection and screening of Beijing 1071-32 isolates along with a validated real-time PCR assay easily deployable across multiple settings for the epidemiological tracking of this significant MDR cluster.

Genomic signatures of drug resistance in highly resistant Mycobacterium tuberculosis strains of the early ancient sublineage of Beijing genotype in Russia

The Mycobacterium tuberculosis Beijing genotype is a clinically and epidemiologically important lineage that is subdivided into ancient/ancestral and modern strains. In our previous study in western Siberia, we identified variable number of tandem repeats (VNTR)-based clusters within the early ancient sublineage of the Beijing genotype characterized by an unexpectedly high rate of extensive drug resistance (XDR). In the current study, next generation sequencing data were analysed to gain insight into genomic signatures underlying drug resistance of these strains. A total of 184 genomes of the Beijing early ancient sublineage from Russia (16), China (15), Japan (36), Korea (25), Vietnam (18), Thailand (73), and the USA (1) were used for phylogenetic analysis. The drug-resistant profile was deduced genotypically. The Russian isolates were distributed into two clusters and were all drug resistant, mainly pre-XDR and XDR. The largest of these clusters included only Russian isolates from remote locations in both Asian and European parts of the country. All its isolates had a quadruple drug resistance (to isoniazid, rifampin, ethambutol and streptomycin) due to the 6-mutation signature (KatG Ser315Thr, KatG Ile335Val, RpoB Ser450Leu, RpoC Asp485Asn, EmbB Gln497Arg, and RpsL Lys43Arg). In most samples, it was complemented with additional and different pncA, gyrA and rrs mutations leading to the pre-XDR/XDR genotype. Phylogenomic analysis indicates a distant origin of this Russian resistant cluster in the early 1970s but location and circumstances are yet to be clarified.

First insight into phylogeography of Mycobacterium bovis and M. caprae from cattle in Bulgaria

Bovine tuberculosis (bTB) represents a significant economic burden to the agriculture. In spite of decades of the control program, Mycobacterium bovis infection levels in cattle in Bulgaria continued to rise over recent years. In order to gain a better understanding of the M. bovis diversity, we used spoligotyping for strain differentiation and the data were compared to the international databases Mbovis.org and SITVIT2 for shared type and clade assignment. Study sample included 30 M. tuberculosis complex isolates from cattle originating from different regions of Bulgaria. The isolates were subdivided by spoligotyping into 4 spoligotypes: 2 types shared by 20 and 8 isolates and 2 singletons. SITVIT2-defined types SIT645 and SIT647 belonged to the common and classical bovine ecotype M. bovis (9 isolates) while types SIT120 and SIT339 belonged to the M. caprae ecotype (21 isolates). A certain phylogeographic gradient of the spoligotypes and clades at the within-country level was observed: M. caprae was prevalent in the central/southwestern, while classical M. bovis in the northeastern Bulgaria. Whereas all four types have global or European circulation, neither was described in the neighboring Balkan countries. M. caprae isolates identified in this study mostly belong to the Central/Eastern European cluster. In summary, this study provided a first insight into phylogeography of M. bovis in Bulgaria and described, for the first time, M. caprae as an important infectious agent of bTB in this country.

Predicting and visualizing features of CRISPR-Cas systems

Pervasive application of CRISPR-Cas systems in genome editing has prompted an increase in both interest and necessity to further elucidate existing systems as well as discover putative novel systems. The ubiquity and power of current computational platforms have made in silico approaches to CRISPR-Cas identification and characterization accessible to a wider audience and increasingly amenable for processing extensive data sets. Here, we describe in silico methods for predicting and visualizing notable features of CRISPR-Cas systems, including Cas domain determination, CRISPR array visualization, and inference of the protospacer-adjacent motif. The efficiency of these tools enables rapid exploration of CRISPR-Cas diversity across prokaryotic genomes and supports scalable analysis of large genomic data sets.

Mycobacterial biomaterials and resources for researchers

There are many resources available to mycobacterial researchers, including culture collections around the world that distribute biomaterials to the general scientific community, genomic and clinical databases, and powerful bioinformatics tools. However, many of these resources may be unknown to the research community. This review article aims to summarize and publicize many of these resources, thus strengthening the quality and reproducibility of mycobacterial research by providing the scientific community access to authenticated and quality-controlled biomaterials and a wealth of information, analytical tools and research opportunities.

CRISPR Visualizer: rapid identification and visualization of CRISPR loci via an automated high-throughput processing pipeline

A CRISPR locus, defined by an array of repeat and spacer elements, constitutes a genetic record of the ceaseless battle between bacteria and viruses, showcasing the genomic integration of spacers acquired from invasive DNA. In particular, iterative spacer acquisitions represent unique evolutionary histories and are often useful for high-resolution bacterial genotyping, including comparative analysis of closely related organisms, clonal lineages, and clinical isolates. Current spacer visualization methods are typically tedious and can require manual data manipulation and curation, including spacer extraction at each CRISPR locus from genomes of interest. Here, we constructed a high-throughput extraction pipeline coupled with a local web-based visualization tool which enables CRISPR spacer and repeat extraction, rapid visualization, graphical comparison, and progressive multiple sequence alignment. We present the bioinformatic pipeline and investigate the loci of reference CRISPR-Cas systems and model organisms in 4 well-characterized subtypes. We illustrate how this analysis uncovers the evolutionary tracks and homology shared between various organisms through visual comparison of CRISPR spacers and repeats, driven through progressive alignments. Due to the ability to process unannotated genome files with minimal preparation and curation, this pipeline can be implemented promptly. Overall, this efficient high-throughput solution supports accelerated analysis of genomic data sets and enables and expedites genotyping efforts based on CRISPR loci.

Evolutionary History and Ongoing Transmission of Phylogenetic Sublineages of Mycobacterium tuberculosis Beijing Genotype in China

Mycobacterium tuberculosis Beijing genotype originated in China and has undergone a dramatic population growth and global spread in the last century. Here, a collection of M. tuberculosis Beijing family isolates from different provinces across all China was genotyped by high-resolution (24-MIRU-VNTR) and low-resolution, high-rank (modern and ancient sublineages) markers. The molecular profiles and global and local phylogenies were compared to the strain phenotype and patient data. The phylogeographic patterns observed in the studied collection demonstrate that large-scale (but not middle/small-scale) distance remains one of the decisive factors of the genetic divergence of M. tuberculosis populations. Analysis of diversity and network topology of the local collections appears to corroborate a recent intriguing hypothesis about Beijing genotype originating in South China. Placing our results within the Eurasian context suggested that important Russian B0/W148 and Asian/Russian A0/94-32 epidemic clones of the Beijing genotype could trace their origins to the northeastern and northwestern regions of China, respectively. The higher clustering of the modern isolates in children and lack of increased MDR rate in any sublineage suggest that not association with drug resistance but other (e.g., speculatively, virulence-related) properties underlie an enhanced dissemination of the evolutionarily recent, modern sublineage of the Beijing genotype in China.

Latin-American-Mediterranean lineage of Mycobacterium tuberculosis: Human traces across pathogen's phylogeography

Currently, Mycobacterium tuberculosis isolates of Latin-American Mediterranean (LAM) family may be detected far beyond the geographic areas that coined its name 15years ago. Here, we established the framework phylogeny of this geographically intriguing and pathobiologically important mycobacterial lineage and hypothesized how human demographics and migration influenced its phylogeography. Phylogenetic analysis of LAM isolates from all continents based on 24 variable number of tandem repeats (VNTR) loci and other markers identified three global sublineages with certain geographic affinities and defined by large deletions RD115, RD174, and by spoligotype SIT33. One minor sublineage (spoligotype SIT388) appears endemic in Japan. One-locus VNTR signatures were established for sublineages and served for their search in published literature and geographic mapping. We suggest that the LAM family originated in the Western Mediterranean region. The most widespread RD115 sublineage seems the most ancient and encompasses genetically and geographically distant branches, including extremely drug resistant KZN in South Africa and LAM-RUS recently widespread across Northern Eurasia. The RD174 sublineage likely started its active spread in Brazil; its earlier branch is relatively dominated by isolates from South America and the derived one is dominated by Portuguese and South/Southeastern African isolates. The relatively most recent SIT33-sublineage is marked with enigmatic gaps and peaks across the Americas and includes South African clade F11/RD761, which likely emerged within the SIT33 subpopulation after its arrival to Africa. In addition to SIT388-sublineage, other deeply rooted, endemic LAM sublineages may exist that remain to be discovered. As a general conclusion, human mass migration appears to be the major factor that shaped the M. tuberculosis phylogeography over large time-spans.

Trends in molecular epidemiology of drug-resistant tuberculosis in Republic of Karelia, Russian Federation

Background

Russian Republic of Karelia is located at the Russian-Finnish border. It contains most of the historical Karelia land inhabited with autochthonous Karels and more recently migrated Russians. Although tuberculosis (TB) incidence in Karelia is decreasing, it remains high (45.8/100 000 in 2014) with the rate of multi-drug resistance (MDR) among newly diagnosed TB patients reaching 46.5 %. The study aimed to genetically characterize Mycobacterium tuberculosis isolates obtained at different time points from TB patients from Karelia to gain insight into the phylogeographic specificity of the circulating genotypes and to assess trends in evolution of drug resistant subpopulations.

Methods

The sample included 150 M. tuberculosis isolates: 78 isolated in 2013–2014 (“new” collection) and 72 isolated in 2006 (“old” collection). Drug susceptibility testing was done by the method of absolute concentrations. Spoligotyping was used to test genotype-specific markers of a Latin-American-Mediterranean (LAM) family and its sublineages as well as a Beijing B0/W148-cluster.

Results

The largest spoligotypes were SIT1 (Beijing family, n = 42) and SIT40 (T family, n = 5). Beijing family was the largest (n = 43) followed by T (n = 11), Ural (n = 10) and LAM (n = 8). Successful Russian clone, Beijing В0/W148, was identified in 15 (34.9 %) of 43 Beijing isolates; all В0/W148 isolates were drug-resistant. Seven of 8 LAM isolates belonged to the RD115/LAM-RUS branch, 1 - to the LAM RD174/RD-Rio sublineage. MDR was found in Beijing (32/43), Ural (3/10), and LAM (3/8). In contrast, all T isolates were pansusceptible. Comparison of drug resistant subgroups of the new and old collections showed an increasing prevalence of the B0/W148 clonal cluster, from 18.0 % (mainly polyresistant) in 2006 to 32.6 % in 2014 (mainly MDR and pre-XDR). The West–east increasing gradient is observed for the Ural genotype that may be defined a ‘Russian’ strain. In contrast, the spoligotype SIT40 of the T family appears to be a historical Karelian strain.

Conclusions

Circulation of the MDR M. tuberculosis isolates of the Beijing genotype and its B0/W148 cluster continues to critically influence the current situation with the MDR-TB control in northwestern Russia including the Republic of Karelia. Revealed phylogeographic patterns of some genotypes reflect a complex demographic history of Karelia within the course of the 20^th century.