First insight into the whole-genome sequence variations in Mycobacterium bovis BCG-1 (Russia) vaccine seed lots and their progeny clinical isolates from children with BCG-induced adverse events

Background The only licensed live Bacille Calmette-Guérin (BCG) vaccine used to prevent severe childhood tuberculosis comprises genetically divergent strains with variable protective efficacy and rates of BCG-induced adverse events. The whole-genome sequencing (WGS) allowed evaluating the genome stability of BCG strains and the impact of spontaneous heterogeneity in seed and commercial lots on the efficacy of BCG-vaccines in different countries. Our study aimed to assess sequence variations and their putative effects on genes and protein functions in the BCG-1 (Russia) seed lots compared to their progeny isolates available from immunocompetent children with BCG-induced disease (mainly, osteitis). Results Based on the WGS data, we analyzed the links between seed lots 361, 367, and 368 used for vaccine manufacture in Russia in different periods, and their nine progeny isolates recovered from immunocompetent children with BCG-induced disease. The complete catalog of variants in genes relative to the reference genome (GenBank: CP013741) included 4 synonymous and 8 nonsynonymous single nucleotide polymorphisms, and 3 frameshift deletions. Seed lot 361 shared variants with 2 of 6 descendant isolates that had higher proportions of such polymorphisms in several genes, including ppsC, eccD5, and eccA5 involved in metabolism and cell wall processes and reportedly associated with virulence in mycobacteria. One isolate preserved variants of its parent seed lot 361 without gain of further changes in the sequence profile within 14 years. Conclusions The background genomic information allowed us for the first time to follow the BCG diversity starting from the freeze-dried seed lots to descendant clinical isolates. Sequence variations in several genes of seed lot 361 did not alter the genomic stability and viability of the vaccine and appeared accumulated in isolates during the survival in the human organism. The impact of the observed variations in the context of association with the development of BCG-induced disease should be evaluated in parallel with the immune status and host genetics. Comparative genomic studies of BCG seed lots and their descendant clinical isolates represent a beneficial approach to better understand the molecular bases of efficacy and adverse events during the long-term survival of BCG in the host organism.


Background
Treatment of tuberculosis (TB) and vaccination remain inevitable health-care interventions to prevent new Mycobacterium tuberculosis infections and their progression to disease. The live Bacille Calmette-Guérin (BCG) vaccine is the only licensed vaccine that is successfully used for nearly a century to prevent childhood TB, although it shows a modest protective effect in adults [1,2]. The history of BCG development and spreading was thoroughly reviewed [3,4]. Briefly, the ancestral form of BCG was derived from wild-type Mycobacterium bovis progenitor strain, while attenuation by numerous in vitro passages resulted in world distribution of parental BCG strain in 1924. During the next decades, maintenance of the BCG strain progenies by further passages using different culture media and transfer schedules in different countries contributed to the in vitro evolution of BCG. The production of seed stocks without lyophilization until the 1960s resulted in the generation of BCG daughter strains (substrains) with variable morphological, culture, antigenic characteristics, and residual virulence due to adaptation in vitro [3][4][5][6][7].
The up-to-date genealogy of widely used BCG vaccine strains is based on genomic variations that occurred due to large deletions, tandem duplications, insertion sequences, and the subsequent series of insertions/deletions and single nucleotide polymorphisms (SNPs). Accordingly, BCG strains were classified into "early" strains represented by BCG Japan, Moreau, Russia (and the descendant BCG Sofia) strains that show fewer chromosomal deletions than "late" strains -Danish, Pasteur, Glaxo, and others [8][9][10][11]. Altogether, these variations underpin the molecular basis of attenuation, different levels of gene expression of immunodominant proteins leading to variable protective efficacy, and different rates of adverse events of currently used genetically divergent BCG strains otherwise considered different BCG vaccines [3,8].
Historically, the original M. bovis BCG strain was transferred to Russia in 1924, and since 1925 it was referred to as M. bovis BCG-1 (also spelled as BCG-I). In the 1940s, M. bovis BCG-1 was lyophilized, and in 1954 the seed lot system was adopted in Russia. The first seed lot used in 1963 for the manufacturing of freeze-dried BCG vaccines in the former Soviet Union was seed lot 352 ("ch"). Since 2006, the 368 ("shch") generation of M. bovis BCG-1 (Russia) strain no. 700001 is in use for the production of the BCG vaccine in Russia [12,13]. The same BCG generation 368 is also adopted by two manufacturers in India for vaccine production and distribution worldwide under the name "BCG Russia 368" [14].
The attenuated BCG vaccines are believed safe and effective mainly to prevent the most severe meningitis and miliary forms of childhood TB associated with high mortality in infants and young children, especially in high TB-burden countries [1,2]. However, adverse events following immunization with BCG in 1-10% of vaccinees are likely to be substantially underreported since assays to confirm isolates as BCG are still relatively rare [15][16][17][18][19][20] despite the advent of molecular techniques and commercially available kits for BCG identification [21]. In Russia, the estimated frequency of complications per 100,000 BCGimmunized decreased from 35.3 in 2005 to 11.2 in 2013 [22].
Vaccine-related complications, i.e., BCG-induced disease (BCG-ID), may occur later than 12 months after vaccination as local site reactions or distal to the site of inoculation, and they appear to differ from cutaneous lesions, lymphadenopathy/lymphadenitis, and hypersensitivity reactions to osteitis and disseminated BCG infections (almost exclusively in individuals with compromised cellular immunity) [1,16,19,[23][24][25].
The BCG-induced adverse reactions are believed to be due to the use of genetically distant vaccine strains with variable residual virulence, the culture technique, the route of administration and the dose of viable cells in the BCG vaccine delivered, along with the age, immune status, and host genetics of the vaccinees [4,8,9,16,17,26,27].
To standardize the vaccine production and to prevent adverse reactions related to BCG vaccination, the WHO has established the international requirements for the manufacture and quality control including genetic characterization of seeds and final commercial lots of the reference strains, i.e., BCG Danish 1331, Tokyo 172-1, and Russia BCG-1 used for vaccine manufacture and distribution worldwide [11,28,29].
Given recommendations, the identity and genome stability of seven seed lots of M. bovis BCG-1 (Russia) used for immunization against TB in Russia since 1948 till present, and 32 commercial lots of the vaccine produced by Russian manufacturers were proved by multiplex polymerase chain reaction [12].
During the last decades, the whole-genome sequencing (WGS) was applied to assess the genome stability of BCG vaccine strains from different culture collections and to evaluate the impact of minor sequence variations occurring in seed and commercial lots on the efficacy of BCG-vaccines used in different countries [30][31][32][33][34]. Recently, sequence variations were identified in clinical isolates from two infants with BCGinduced disease compared with the sub-cultured M. bovis BCG Moscow strain (Serum Institute of India, India) in South Africa [15].
The complete genome sequences of the current M. bovis BCG-1 (Russia) 368 ("shch") freeze-dried seed lot and its working seed lot used for vaccine production by one of the Russian manufacturers are available in Gen-Bank (NCBI) under accession numbers CP011455 [13], CP009243 [35], and CP013741 [36]. The genome constancy of M. bovis BCG-1 (Russia) 368 within the entire vaccine production flow (from working seed lot, the last production passage to final vaccine product) was confirmed by spoligotyping, 24 loci VNTR-typing, and WGS [36,37]. However, so far, there is a lack of data on the contribution of BCG genomic diversity to the viability, fitness, and residual virulence of vaccine populations in vivo in the context of the development of BCGinduced adverse events. Therefore, our study aimed to assess sequence variations and their putative effects on genes and protein functions in the BCG-1 (Russia) seed lots compared to their progeny clinical isolates available from immunocompetent children with BCG-induced osteitis and lymphadenitis.

Results
Sequence variations identified in BCG-1 (Russia) parent seed lots and their progeny clinical isolates Raw reads of seed lots 361, 367, 368, and descendant clinical isolates were mapped to the reference genome of BCG-1 (Russia) (GenBank: CP013741). Variant calling allowed to produce the complete catalog of 15 differences, including 12 SNPs (four synonymous and eight nonsynonymous) and three disruptive deletions (frameshift variants) in CDS of the seed lots and clinical isolates relative to the reference genome. Overview of sequence variation patterns identified in BCG seed lots and their progeny clinical isolates are presented in Table 1, Additional Files 1 and 2.
His57Asp mutation in pncA (Rv2043c) gene conferring high intrinsic resistance to pyrazinamide in M. bovis was detected in all BCG-1 (Russia) progeny isolates. A single-nucleotide insertion in the recA gene was not detected in any of three seed lots, nor their descendant clinical isolates.
Comparative analysis of sequence variants in BCG-1 (Russia) seed lots and their progeny clinical isolates We examined a set of nine BCG-1 (Russia) clinical isolates recovered from duly BCG-vaccinated immunocompetent children (except one, immunized with delay for 4 months after birth) with a culture-confirmed diagnosis of BCG-ID. For each case, analysis of medical records allowed to follow the history of primary BCG vaccination, including information on commercial vaccine products used for intradermal inoculation, its manufacture, and finally, the parent seed lot (Tables 2 and 3).
BCG clinical isolates shared the majority of genomic features with their parent seed lots. However, the proportions of reads with alternative alleles varied in BCG seed lots, and their progeny isolates. These findings were summarized in Table 1 and Additional file 1. Based on the comparative analyses of the data obtained, we built a schematic view of the relationships between BCG seed lots and their descendant clinical isolates (Fig. 1).
The population of BCG-1 (Russia) seed lot 361 was polymorphic with variable proportions of SNPs and deletions in reads spanning seven CDS (Table 1, Fig. 1, Additional file 1).
Synonymous SNP in the ornithine aminotransferase gene (rocD) was detected in almost 1/3 reads of seed lot 361, but it was absent in the descendant isolates. Alternatively, seed lot 361 shared nonsynonymous (missense) SNPs in type VII secretion system ESX-5 (eccD5 and eccA5), dolichol-phosphate mannosyltransferase (ppm1) and hypothetical protease genes with its progeny clinical isolates 2925 and 5448 that had higher proportions (up to 100%) of alternative variants. However, sequences of these isolates differed in amino acid transporter (cycA) and polyketide synthase (ppsC) genes. A nonsynonymous SNP in the PE-PGRS family protein PE_PGRS7 gene was found in 38% reads of isolate 5340, the progeny of seed lot 361.
Isolates 1986, 4159, and 5075, progenies of the seed lot 361, and seed lot 367 did not show sequence polymorphisms compared to the reference sequence, i.e., BCG-1 (Russia) current seed lot 368 (Table 1, Additional file 1).  Isolate 3363, a single progeny of the seed lot 367, carried a missense variant in amidase (amiD) gene and a frameshift deletion in the thioesterase (pks13) gene (Table 1, Fig. 1, Additional file 1).
The sequenced seed lot 368 shared a single base deletion in methionine aminopeptidase (mapA) with descendant isolate 577, also demonstrating a polymorphism (43% reads) in the ATP synthase subunit beta gene (atpD). Polymorphism in the N-acetylglucosaminyldiphospho-decaprenol L-rhamnosyltransferase gene (wbbL1) and a frameshift deletion in the universal stress protein UspA gene were identified in clinical isolate 1032 (18 and 28% reads, respectively), the progeny of seed lot 368 (Table 1, Fig. 1; Additional file 1 and 3).

Discussion
All nine BCG clinical isolates recovered from tissue samples taken during the surgery were susceptible to first-line antibiotics: streptomycin, isoniazid, rifampin, ethambutol. However, they were pyrazinamide-resistant due to the typical His57Asp mutation in pncA (Rv2043c) gene conferring high intrinsic resistance in M. bovis, and consequently, all BCG vaccine strains [38].
Previously, spontaneous heterogeneity of BCG seed lots and commercial vaccines during vaccine production determined by variant calling was analyzed in the Tokyo-172 vaccine strain [32]. In line with the published data, we observed a synonymous A2139G single nucleotide base substitution (position 3,175,529, "-" chain) in BOVR_15,090, i.e., the glnD (ortholog Rv2918c in M. tuberculosis H37Rv) gene, without Ala replacement in position 713 of 808 amino acids in the conservative domain of the PII uridylyltransferase protein in all genomes studied. The G variant was present in 100% reads of sequenced seed lots 361 and 367, and eight of nine clinical isolates. However, it was detected only in 25% reads of current seed lot 368 (sequenced in this study), and 46% reads of its progeny isolate 577.
These findings inspired further examination of A/G polymorphisms in the glnD in other sequenced BCG genomes.
The G variant was previously reported in the glnD in BCG-1 (Russia) current seed lot 368 and some earlier   vaccine products [35], though alternative variant A was detected in two other BCG-1 (Russia) 368 complete genomes [13,36]. The stability of variant A in the glnD gene of the BCG-1 (Russia) 368 working seed lot (Gen-Bank: CP013741) and consistent vaccine batches within 2 years was confirmed by WGS [36,37]. Noteworthy, all three complete genomes of BCG-1 (Russia) were sequenced independently from the same generation 368 ("shch"), which is still in use since 2008 for the BCG vaccine production by two Russian manufactures: Medgamal (Moscow) and Microgen (Stavropol). Interestingly, allele G was recorded in the genomes of BCG str. Tokyo 172 which represents a well-defined Type I subpopulation of BCG Tokyo-172, the closest to BCG Russia (GenBank: NZ_CUWO01000001) and Russia ATCC 35740, members of presumably monophyletic "early" substrain group of BCG [10,30,32]. The G variant was also common in the reference sequences of some other BCG strains, e.g., Pasteur 1173P2, Danish 1331, ATCC 35743 (Tice), Moreau RDJ, used worldwide for vaccine production and research purposes.
Altogether, these data suggest that intrinsic polymorphism in the glnD gene of BCG-1 (Russia) strain represented by different seed lots and descendant clinical isolates occurs due to the selection of variants during laboratory and vaccine production manipulations, rather than mutations in patients' organisms. Otherwise, the observed heterogeneity can be partly explained as a result of processing nucleotide variants among the mapped reads using different bioinformatics tools and quality criteria.
Another interesting finding: a notorious singlenucleotide insertion in the recA gene leading to nonfunctional truncated recombinase A reported previously in a vaccine production lot of M. bovis BCG Russia (corresponding to ATCC 35740) [39] was not detected in seed lots studied, nor their descendant clinical isolates. De facto, there were no alterations in recA in the  Table 1 for more details previously sequenced complete genomes of BCG-1 (Russia) [13,35,36] along with other hyperconserved BCG strains [33].
Based on the comparative study of WGS data, we attempted to analyze the links between BCG-1 (Russia) seed lots used for vaccine manufacture and their progeny clinical isolates associated with specific adverse events, e.g., BCG-osteitis and an abscess at the point of injection.
The population of BCG-1 (Russia) seed lot 361 appeared heterogeneous with variable proportions of SNPs and deletions in reads spanning seven CDS. Seed lot 361 shared SNPs in eccD5, eccA5, ppm1 genes, and hypothetical protease gene with its progeny clinical isolates 2925 and 5448. The latter typically had higher proportions (reaching up to 100%) of alternative alleles than the parent seed lot 361, thus suggesting the likely accumulation of such changes under selective pressure in the host organism. The similar SNPs we observed in corresponding genomic positions in M. tuberculosis H37Rv, BCG strains Tokyo-172 type I and Pasteur 1173P2. These polymorphisms affected type VII specialized bacterial secretion ESX-5 (integral membrane protein EccD and secretion AAA-ATPase), reportedly associated with virulence in mycobacteria, and dolichol-phosphate mannosyltransferase and protease genes involved in the vital cell processes (protein modification, lipoprotein biosynthesis), intermediary metabolism, and respiration [8,40,41]. The data obtained confirmed the close genetic relatedness between seed lot 361 and both progeny isolates 2925 and 5448 that otherwise differed from each other in amino acid transporter (cycA) and polyketide synthase (ppsC) sequences.
A unique missense variant in the PE-PGRS family protein PE_PGRS7 gene with an unclear function in antigenic variation and interactions of bacteria with the host immune system was found in 38% reads of isolate 5340, the progeny of seed lot 361.
The sequences of seed lot 367 (used for vaccine production in 2001-2008) also did not show variations in 14 of 15 CDS compared to the reference genome and thus appeared the most conservative in our study. However, its single progeny, isolate 3363, carried a missense variant in amidase (amiD) gene essential for protein translation and required for mycobacterial survival during infection (virulence factor) and a frameshift deletion in the thioesterase region of polyketide synthase 13 (pks13) gene involved in the final assembly of mycolic acids, major and specific lipid components of the envelope essential for the survival of a mycobacterial cell [40,42].
The population of the sequenced seed lot 368 was almost homogenous containing a single base deletion (69% reads) in methionine aminopeptidase (mapA) gene compared to the previously reported complete genomes of BCG-1 (Russia) current-generation 368 [13,35,36] and BCG Tokyo 172 type 1 [32]. This frameshift leading to protein truncation was inherited in the descendant isolate 577 (47% reads), also bearing an isolate-specific allele (43% reads) in the ATP synthase subunit beta gene (atpD). Both proteins encoded by mapA and atpD genes are involved in the processes of intermediary metabolism and respiration [8].
A synonymous SNP in the N-acetylglucosaminyldiphospho-decaprenol L-rhamnosyltransferase gene (wbbL1) not altering the enzyme structure and function in the cell wall and cell processes was identified in clinical isolate 1032 (18% reads). Another finding in isolate 1032 (28% reads) was a frameshift due to a single-base cytosine deletion in codon 630 of the universal stress protein UspA gene (orthologue Rv1996 in M. tuberculosis H37Rv), one of eight Usp paralogues essential for the intracellular survival of pathogenic mycobacteria in hypoxia, low levels of nitric oxide and carbon monoxide environment. This deletion completely altering 27 amino acids from the position 212 onwards resulted in a gained opal stop codon (TGA). The latter was leading to a premature protein truncation in the residue 239 yielding to the loss of 3 ligand-binding loci in the putative ATP binding motif of the UspA domain 2 [43].
It appears that the accumulation of shared SNPs and unique variations occurred in clinical isolates 577 and 1032 (both progenies of the sequenced seed lot 368), 3363 (progeny of seed lot 367), and 5340 (progeny of seed lot 361) did not affect the vaccine survival in the organisms of vaccinees. Moreover, none of the frameshifts, which presumably could alter essential cell functions, were crucial for the microbial population in vitro, nor in vivo, thus emphasizing the long-term viability of BCG-1 (Russia) strain in the human organism.
Adverse events following immunization with BCG are considered underreported since assays to confirm isolates as BCG are still relatively rare. Moreover, based on the routine laboratory diagnostics and even molecular techniques, it is difficult to associate the clinical BCG isolates with the particular vaccine lot administered [15][16][17][18][19][20]. Therefore, the contribution of BCG genomic heterogeneity to the viability and residual virulence of vaccine populations in vivo remains unclear in terms of the development of BCG-related disease that otherwise can be associated with immune responses and human genetics.
Unfortunately, the administered commercial vaccines were unavailable for WGS, which limited the comparative analysis of sequence variations through the whole vaccine production chain. Nevertheless, bioinformatics tools provided a possibility to identify sequence variants and their putative effects on genes and protein functions in BCG-1(Russia) seed lots 361 and 368 used for vaccine manufacture in Russia in different periods relative to their progeny clinical isolates from patients with confirmed BCG-related disease and the genomes of internationally recognized BCG reference strains. BCG osteitis is a rare but severe adverse event that typically presents in children before the age of 5 years, and there are few reports concerning the late onset of the disease in teenagers [17,19]. In our setting, we analyzed confirmed cases of BCG-osteitis typically developed at the age under 5 years (median 22 months) among immunocompetent children immunized with BCG-1 (Russia) vaccine after, or soon after birth. The latest onset of the osteitis affecting the left eighth rib was diagnosed in the patient aged 14 years who once received as a newborn the BCG vaccine derived from the seed lot 361. The isolate 2925 recovered from tissue samples taken during surgery accumulated four nonsynonymous and one synonymous SNPs assigned to its parent seed lot 361 but did not gain any further sequence variants in the patient's organism. Noteworthy, such a prolonged period of survival in a human organism of the life BCG vaccine without significant changes in its genome appears remarkable.

Conclusions
To our knowledge, we reported for the first time the results of a WGS-based comparative analysis of the BCG-1 (Russia) seed lots used for the vaccine production in different periods and their progeny clinical isolates recovered from immunocompetent children with BCG-induced adverse reactions (mainly, BCG osteitis). Our study provided the background genomic information, which allowed us to follow the BCG diversity starting from the freeze-dried seed lots to descendant clinical isolates presenting the results of the long-term survival of vaccine populations in the human organism.
Two of three BCG-1 (Russia) seed lots were heterogeneous with various proportions of sequence variants in several genes, including ppsC, eccD5, and eccA5 involved in metabolism and cell wall processes and reportedly associated with virulence in mycobacteria. However, such polymorphisms did not alter the genomic stability and viability of BCG-1 (Russia)-derived vaccines in vivo. Sequence variants occurred in polymorphic CDS of seed lot 361 used for vaccine manufacture before 2006 appeared accumulated in two of six progeny isolates, hypothetically, under selective pressure in the human organism. Although polymorphisms were identified in two unrelated clinical isolates 2925 and 5448 from children immunized in different years with different vaccine lots produced from the BCG-1 (Russia) seed lot 361 by the same manufacturer, the impact of such changes in the context of association with the development of BCG-induced disease remains questionable due to the lack of detailed information on the immune status and host genetics of presumably immunocompetent patients with BCG-osteitis.
Despite promising results in the development of a new generation of BCG vaccines, there is no reliable alternative to those which are still in use in the nearest future.
The contributions of genomic diversity and variations in gene expression profiles during the BCG vaccine manufacture flow and adaptations in vivo are yet to be uncovered in parallel with the molecular mechanisms of innate and adaptive immune responses in the human organism. Thereby, the comparative studies of genomic variations in production strains, their seeds, administered vaccine lots, and the descendant clinical isolates (if available) represent a beneficial approach to better understand the bases of vaccine efficacy and adverse reactions of present and future BCG-based vaccines at genomic, transcriptomic, and proteomic levels.

Study design
We retrospectively reviewed 65 medical records at the clinics of St. Petersburg Research Institute of Phthisiopulmonology selected as codes A18.0 (Tuberculosis of bones and joints, acute or chronic osteomyelitis) and A18.2 (Tuberculosis peripheral lymphadenopathy) according to the International Statistical Classification of Diseases and Related Health Problems (ICD-10). The probable BCG-ID was suspected in nine immunocompetent children (eight patients aged under 5 years and the 14-year old patient) with developed osteitis (eight patients), or cutaneous lesion near the BCG inoculation site (one patient), welldefined BCG inoculation records, and without a history of TB contacts admitted for treatment (including surgical intervention) in 2006-2018 (Table 2).

Inclusion criteria
The definition for the diagnosis of BCG-ID ("BCG-itis") included 1) the histology showing granulation tissue with caseous necrosis in biopsy samples taken while surgery and 2) mycobacteria cultures recovered from tissue samples on the Löwenstein-Jensen (L-J) medium if they were identified as M. bovis BCG using GenoType MTBC (Hain Lifescience, Germany). The M. bovis BCG cultures were stored at − 80°C for further examination.

Drug susceptibility testing
Drug susceptibility to first-line antibiotics (streptomycin, isoniazid, rifampin, ethambutol) and pyrazinamide of the M. bovis BCG cultures was assessed by using conventional absolute concentration method and/or BACTEC MGIT 960 System (Becton-Dickinson, Sparks, MD, USA).

BCG-1 (Russia) seed lots and clinical isolates used for WGS
For genomic analysis, we used genomic DNA extracted as described [44] from M. bovis BCG-1 (Russia) freezedried seed lots and nine BCG clinical isolates cultured on L-J medium, each recovered from a child patient with confirmed BCG-ID.
The lyophilized seed lots of three generations 361 ("sh"), 367 ("shch"), and 368 ("shch") of M. bovis BCG-1 (Russia) substrain No. 700001 in ampoules were obtained from the National State Collection of Pathogenic Microorganisms (NSCPM) of Scientific Center for Expertise of Medical Application Products of the Ministry of Health, Russia along with the characteristics (the origin and manufacture) of corresponding BCG vaccine commercial lots used for immunization (Table 3).
The pncA gene was searched for mutations conferring resistance to pyrazinamide in BCG clinical isolates [38].
Additional file 2: Figure S1. Comparative circular genome diagram of M. bovis BCG 1 (Russia) vaccine seeds and clinical isolates (circos plot was created with Circa (http://omgenomics.com/circa). Tracks from inside: scale in megabases; G.C. content histogram (yellow fill with black stroke) calculated from the ratio (G + C)/(A + T + G + C) using a 1 kilobase (kb) non-overlapping sliding window; G.C. bias plots calculated from the ratio (G C)/(G + C) using a 1 kb non-overlapping window (blue plot) and a 10 kb non-overlapping window (translucent black plot); strain number labeled tracks (alternating grey and white) each depicting strain-specific insertions/deletions (blue/red transverse marks) and SNPs (dots*) identified relative to the reference genome of M. bovis BCG-1 (Russia) (GenBank accession number CP013741); tracks depicting ORFs (blue transverse marks) alongside affected CDSs (black, orange and red transverse marks † ) on forward and reverse strands; names of affected genes colored according to an impact type of the most severe variant affecting a particular gene. Notes. *Sizes and colors determine a particular SNP type: small blacksynonymous SNP, large orangemissense SNP; † color determines an impact type of a particular variant on a particular CDS assigned according to the following concept: low impact (black mark)synonymous SNPs, moderate impact (orange mark)missense SNPs and conservative insertions/deletions, high impact (red mark)nonsense SNPs and frameshift