Genomic analysis of the B12 biosynthetic pathway
Sequences of each of the 16 cob genes present in M. tuberculosis H37Rv 1 were aligned using the NCBI Nucleotide BLAST (BLASTn) tool, against the rest of the MTBC genomes available. Mutations were assigned to specific lineages based on inspection of literature and/or description of the sequenced strains. Experimental confirmations of specific mutations were performed in clinical isolates from our laboratory collection previously assigned to MTBC lineages. After DNA extraction from each isolate, the 16 cob genes were PCR amplified and sequenced by Sanger methodology to detect polymorphisms.
Bacterial strains and growth conditions
Mycobacteria were cultured at 37ºC in Middlebrook 7H9 liquid medium (Difco) containing 0.05% (vol/vol) Tween-80 (Sigma) and supplemented with 10% (vol/vol) of ADC (0,5% bovine serum albumin, 0.2% dextrose, 0.085% NaCL and 0.0003% beef catalase) growth supplement (Middlebrook). For solid media, Middlebrook 7H10 broth containing 10% (vol/vol) ADC was used. For B12 uptake in vitro experiments AdoCbl (Sigma-Aldrich) or CNCbl (Sigma-Aldrich) were added to the liquid medium at a concentration of 10 mg/mL. In in vitro growth assays, both liquid and solid were supplemented with AdoCbl and L-Methionine (Sigma-Aldrich) when required at concentrations of 10 mg/mL and 25 mg/mL, respectively. When necessary, 50 mg/mL hygromycin (Hyg) or 20 mg/mL kanamycin (Kan) were added to the media to select mutant strains. Bacterial suspensions of M. tuberculosis H37Rv, M. canettii C59, and their derivatives used for intranasal infections of mice were properly diluted in PBS from previously quantitated glycerol stock solutions. Escherichia coli was cultured in liquid media Luria-Bertani (LB) broth or in solid media LB agar at 37ºC, or at 30ºC for the recombineering event in the case of pKD46 temperature-sensitive plasmid containing strains. When required, media were supplemented with 100 mg/mL ampicillin (Amp), 12.5 mg/mL chloramphenicol (Cm) or 20 mg/mL Kan. A detailed description of each strain used in the study is provided in (Table 1).
Measures of endogenous B12 production and exogenous B12 uptake by bacteria
For B12 production assays, strains from M. canettii, MTBC, environmental, and oportunistic mycobacteria were cultured at 37ºC until log phase (OD600=0,6) in 10 ml 7H9-0.05% Tween-80-ADC. Cells were harvested by centrifugation at 4000 x g for 10 min and extensively washed three times with sterile PBS 1x. Pellets were resuspended in PBS containing 1% Triton X-100 and transferred into tubes containing glass beads (MP Biomedicals). Suspensions were disrupted by Fast-Prep (6.5 m/s, 45 s) twice and samples were cooled on ice between the cycles. Supernatants containing soluble fractions were filtered through a 0.22 µm-pore-size filter after cold centrifugation at 14000 x g for 5 min. Protein content in the whole-cell extracts were quantified using QuantiPro BCA assay (Sigma Aldrich) for subsequent normalization of B12 measurements. Bacterial whole-cell extracts were used for the quantification of B12 levels in a Cobas analyser from the haematology section of the University Clinical Hospital “Lozano Blesa” (Zaragoza) based on a competitive electrochemiluminescent binding assay. The values of total B12 obtained in pg/mL were normalized against the total protein values in mg/mL. For B12 uptake assays M. canettii and MTBC strains were cultured at 37ºC in 30 ml 7H9-0.05% Tween-80-ADC supplemented when required with AdoCbl (10 mg/mL) and CNCbl (10 mg/mL). Ten mL of culture were taken at day 2, day 7 and 2 months, and bacteria were processed as reported for B12 production studies.
RNA extraction and quantitative Reverse Transcription-PCR (qRT-PCR)
These protocols were previously described in 2. Sequences of primers used are specified in (Table 2). Melting curves were performed for each amplicon to verify specificity. Four replicates of each gene CT value were obtained and normalized to the CT of the sigA housekeeping gene (amplified from the same samples), obtaining a ΔCT = CT, j –CT, sigA, where j is a gene different from sigA. We calculated a ΔΔCT value specific for the B21-supplementation condition (ΔΔCT(+ B12)), subtracting the ΔCT mean value for each gene obtained from non-B12-supplemented condition from the ΔCT obtained for each gene from + B12 condition. Finally, change of expression, denoted as Relative quantity, was calculated with the equation \({2}^{-\varDelta \varDelta {C}_{T}(+B12)}\).
RNA-seq analysis
RNA-seq analysis were performed by STAB-vida (Portugal). First, the library construction of cDNA molecules was carried out using a Ribosomal Depletion Library Preparation Kit. The generated DNA fragments (DNA library) were sequenced in the lllumina Novaseq platform, using 150bp paired-end sequencing reads, and the analysis of the generated raw sequence data was carried out using CLC Genomics Workbench 12.0.3. The bioinformatics analysis started with trimming of raw sequences to generate high quality data only. The high-quality sequencing reads were mapped against the reference genomes, M. tuberculosis H37Rv (NC_ 000962.3) or M. bovis AF2122/97 (NC_002945.3) when appropriate, using the following parameters: length fraction = 0.8 and similarity fraction = 0.8. The result of mapping served to determine the gene expression levels based on the TPM (Trancripts per Million). To analyse variation and bring out strong patterns in the dataset Principal Component Analysis was performed. Differential expression analysis among samples from each strain and replicate was carried out. Fold changes were calculated from the Generalized Linear Model GLM, which corrects for differences in library size between the samples and the effects of confounding factors. The differentially expressed genes were filtered using the condition "Fold change ≥ 2 or ≤-2” and the genes that fulfilled that condition were listed (Tables S1 to S3).
Targeted proteomics by SRM/MS
The SRM/MS approach was used for the quantification of specific PrpC, PrpD and MetE peptides in whole-cell extracts of M. tuberculosis H37Rv and GC1237. Cultures were grown until logarithmic growth in 35 ml 7H9-0.05% Tween-80-ADC supplemented with or without AdoCbl. Cells were washed three times with sterile PBS 1x and pelleted. Pellets were resuspended in a volume equivalent to 10% of the initial culture volume with PBS containing 1% Triton X-100 and transferred into tubes containing glass beads (MP Biomedicals). Suspensions were disrupted by Fast-Prep (6.5 m/s, 45 s) twice and samples were cooled on ice between the cycles. Supernatants containing soluble proteins were filtered through a 0.22 µm-pore-size filter after cold centrifugation at 14000 x g for 5 min. 1.5 mL of each bacterial extract was taken and a volume equivalent to 10% of the initial volume of TCA. After overnight cold precipitation of the samples and a long cold centrifugation at 14000 x g, the supernatants were discarded, and the pellets washed with acetone. A second centrifugation was performed for 10 min and the supernatants were discarded again. The pellets were left at room temperature until completely dry. Finally, they were resuspended in 1 mL of 100 mM Tris solution and the total proteins present in the final samples were quantified using QuantiPro BCA assay (Sigma Aldrich). The SRM/MS was performed as previously described 3. The identified and quantified peptides were as follows: for PrpC, GLAGVVVDTTAISK, GELPTDAELALFSQR and VVPQTNSLTYR; for PrpD, IIDNAAVSAASMVR, FTELADGVVEPVEQQR and VPLPAPGEPK; and for MetE, SWLAFGAEK, VPSAEEMADSLR, IEAIVASGAHR and NVDEVTASLHNMVAAAR.
Construction of mutant strains in M. canettii and M. tuberculosis
Deletion of cobM and cobK genes in M. canettii C59 was achieved using the BAC-recombineering (BAC-rec) strategy 4. Briefly, the thermosensitive plasmid pKD46 containing the red recombinase from lambda phage 5 was co-transformed into the E. coli DH10B clone carrying BAC Rv412 (containing cobMK genes) 6. DH10B Rv412 pKD46 transformants incubated with arabinose 0.15% were subsequently transformed with a PCR product obtained using “KO BAC cobMK Mcan-P1 Fw” and “KO BAC cobMK Mcan P1-Rv” primers. This PCR product consists on a Kan resistance cassette (KanR), with FRT sites from pKD4, flanked by 40 bp identity arms to the genes of interest. Recombinants were selected on LB agar containing Kan incubated at 30ºC overnight, and genes deletion in the BAC was confirmed by PCR amplification using “Conf-KO BAC cobMK Mcan-Fw” and “P1 inv” primers for the 5’ end and “P2 inv long” and “Conf-KO BAC cobMK Mcan-Rv” for the 3’ end. Allelic exchange substrates (AES) containing the KanR flanked by approximately 1kb identity arms for site-specific recombination were obtained by high fidelity PCR using BAC Rv412-ΔcobM,K::KanR as template and “Conf-KO BAC cobMK-Fw”/“Conf-KO BAC cobMK-Rv” primers. AESs sequences are shown in (Supplementary File 1). AES were transformed by electroporation in M. canettii C59 carrying pJV53H recombineering plasmid 7 and cultured in the presence of 0.2% acetamide. The BAC-rec strategy was also used in M. tuberculosis H37Rv for construction of a ΔmetH::KanR mutant. We followed the procedure described above but starting from BAC Rv73 containing the gene Rv2124c (metH). The AES (Supplementary File 1) were obtained using the BAC-knockouts as templates for a PCR using primers described in (Table 2). Then, M. tuberculosis H37Rv carrying the pJV53H was electroporated with the AES and recombinant colonies were selected and confirmed. The M. tuberculosis H37Rv ΔmetE::KanR mutant was constructed using a synthetic AES (GenScript) (Supplementary File 1) introduced in H37Rv-pJV53H by electroporation. To favour the recombination event and recovery of recombinants, the transformation mixtures were incubated overnight in liquid medium without antibiotic and in the presence of AdoCbl, since MetH requires B12 as cofactor. The AESs used for the construction of H37Rv ΔmetE::KanR and ΔmetH::KanR mutants were also electroporated in the M. canettii C59 strain containing pJV53H in order to obtain metE and metH mutants in this strain following the same procedure explained above for M. tuberculosis. Finally, to construct the double recombinant M. tuberculosis H37Rv ΔmetH::KanR PrAg85a_metE (HygR), a pMV361H integrative plasmid containing the metE gene regulated by the promoter of fbpA (Ag85a) (Supplementary File 1) was electroporated in H37Rv ΔmetH::KanR.
For the final steps of mutant construction and confirmation, recombinant colonies were selected by plating on 7H10-ADC with the appropriate antibiotic/supplement and confirmed by PCR using specific primers (Table 2 and Figure S12). All plasmids used in this work are listed in (Table 3).
Construction of mutant strains in M. smegmatis
To construct a metE knockout in M. smegmatis mc2155 and to mimic the B12 phenotype of M. tuberculosis, B12 synthesis in mc2155 was first abrogated by deletion of cobLMK genes. M. smegmatis mc2155 carrying pJV53H 7 were electroporated with AES containing the KanR flanked by approximately 50 bp identity arms for site-specific recombination. AES (Supplementary File 1) were synthesized by PCR amplification of the KanR cassette of pKD4 using primers containing 50 bp identity arms (Table 2). Transformants were plated on 7H10-ADC-Kan plates and recombinants were confirmed by PCR using specific primers (Table 2 and Figure S12)
Once constructed the M. smegmatis mc2155 DcobLMK strain, and after confirmation loss of pJV53H, this strain was electroporated with the pRES-FLP-Mtb plasmid to resolve the KanR resistance cassette. The resolved strain was transformed again with the pJV53H to subsequently generate double mutants. For the construction of the metE mutant, the AES -which contained a KanR cassette and the GFP gene optimized for its expression in mycobacteria (eGFP) (Supplementary File 1)- was synthesized by GenScript. Construction was confirmed by PCR (Figure S12) of transformant colonies grown on 7H10-ADC-Kan-L-Met plates.
Growth of methionine synthesis mutants in laboratory cultures
The L-Met synthesis mutants were first grown in 7H9-Tween-ADC liquid medium, both with and without B12 (AdoCbl, 10 µg/mL), to assess their growth in the presence and absence of the vitamin and thus validate their phenotypes. To evaluate growth on solid media, in all cases -with the exception of the ΔmetE mutant, which only grew in the presence of B12 and, therefore, plating with and without B12 was performed from liquid cultures supplemented with B12- the liquid cultures grown without B12 were plated in solid 7H10-ADC plates supplemented with or without AdoCbl (10 µg/mL) and/or L-Met (25 µg/mL).
Riboswitch structure prediction
Prediction of the secondary structure of the metE B12-riboswitch was based on previous studies 8,9. Using these studies, we identified invariant residues across approximately 200 B12-riboswitches, and M. tuberculosis-specific residues. We also identified a conserved B12-box.
Mouse model of B12 deficiency
Female C57BL/6 and SCID mice were fed with a diet similar to control group normal diet (Teklad 2014S, ENVIGO), but with no added B12 and with 5% Pectin (Tekland custom diet TD.170206, ENVIGO). Note that pectin has been shown to bind vitamin B-12 in the intestine making it less bioavailable 10. Moreover, cages were renewed more frequent than usual to avoid coprophagia, as it could be an additional source of B12. After 8 weeks of feeding, concentrations of vitamin B12 were measured in plasma. For this, mice were euthanized, and blood was extracted by cardiac puncture. The blood extractions were incubated at 4ºC for 24h and after centrifugation for 5 min at 4000 x g, the upper phase containing the serum was extracted. Serum samples were analyzed by electrochemiluminescence at the University Clinical Hospital “Lozano Blesa” to quantify the B12-levels. For B12 supplementation experiments, we used the previously mentioned B12-deficient diet (Diet #1 in the text; Tekland custom diet TD.170206, ENVIGO), and a second independent diet reproducing the composition of the standard 2014S diet which was supplemented with a vitamin mix without B12 (Vitamin A, Vitamin D3, Vitamin E, Vitamin K3, Vitamin B1, Vitamin B2, Niacin, Vitamin B6, Pantothenic Acid, Biotin, Folate and Choline) to ensure nutritional adequacy after autoclaving (Diet #2 in the text). Subcutaneous B12 supplementation of SCID mice subsequent to M. tuberculosis infection, in the corresponding group (Diet #1), consisted on 2.5 µg CNCbl (Optovite B12, Normon laboratories) per week and per mouse, which translates the treatment prescribed for patients with B12 deficits into the animal model.
Mouse infection experiments
All mice were kept under controlled conditions and observed for any sign of disease. Immunocompromised SCID and immunocompetent C57BL/6 female mice were infected with a low dose (≈ 200 CFUs) of the different M. tuberculosis wild type and mutant strains and with a high dose (≈ 50000 CFUs) of the M. canettii C59 wild type and the ΔcobMK mutant strain by the intranasal route. The animals were anesthetized by inhalation route with Isofluorane (Isboa Vet) using a vaporizer and intranasal administration was performed with two instillations of 20 µl of the bacterial suspension prepared in PBS. Bacterial suspensions for infection were plated in solid agar medium to confirm the CFUs used for in vivo challenges. For survival experiments, SCID mice were monitored daily for the development of clinical signs of disease and examined in case any abnormality of behavior was observed. Weight of mice was followed during the experiment. The euthanasia endpoint was defined at the point that the loss of weight was more than 20%. Four weeks postinoculation, the bacterial burden was evaluated in the lungs and spleen of C57BL/6 mice. For this, the organs were aseptically removed and homogenized in 1 ml of H2O using a GentleMacs dissociator (Miltenyi Biotec). CFUs were determined by plating serial dilutions onto 7H10-ADC plates supplemented with AdoCbl when required.
Ethics statement
Experimental animal studies were performed in agreement with European and national directives for the protection of animals for experimental purposes. All procedures were carried out under Project License PI 31/21, approved by the Ethics Committee for Animal Experiments from the University of Zaragoza.
Statistics
Mice were randomly distributed in groups of 6 animals per cage prior to start experimental procedures. Results were not blinded for analysis. No statistical method was used to calculate sample size in animal experiments. GraphPrism software was used for statistical analysis. Statistical tests used for each experiment are indicated in the figure legends. All statistical tests were two-tailed and differences were considered significant at p < 0.05.
METHODS AND TABLES REFERENCES
1 Cole, S. T. et al. Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature 393, 537-544, doi:10.1038/31159 (1998).
2 Perez, I. et al. Live attenuated TB vaccines representing the three modern Mycobacterium tuberculosis lineages reveal that the Euro-American genetic background confers optimal vaccine potential. EBioMedicine 55, 102761, doi:10.1016/j.ebiom.2020.102761 (2020).
3 Perez, I. et al. The Mycobacterium tuberculosis PhoPR virulence system regulates expression of the universal second messenger c-di-AMP and impacts vaccine safety and efficacy. Mol Ther Nucleic Acids 27, 1235-1248, doi:10.1016/j.omtn.2022.02.011 (2022).
4 Aguilo, N. et al. Reactogenicity to major tuberculosis antigens absent in BCG is linked to improved protection against Mycobacterium tuberculosis. Nat Commun 8, 16085, doi:10.1038/ncomms16085 (2017).
5 Datsenko, K. A. & Wanner, B. L. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci U S A 97, 6640-6645, doi:10.1073/pnas.120163297 (2000).
6 Brosch, R. et al. Use of a Mycobacterium tuberculosis H37Rv bacterial artificial chromosome library for genome mapping, sequencing, and comparative genomics. Infect Immun 66, 2221-2229, doi:10.1128/iai.66.5.2221-2229.1998 (1998).
7 van Kessel, J. C. & Hatfull, G. F. Recombineering in Mycobacterium tuberculosis. Nat Methods 4, 147-152, doi:10.1038/nmeth996 (2007).
8 Vitreschak, A. G., Rodionov, D. A., Mironov, A. A. & Gelfand, M. S. Regulation of the vitamin B12 metabolism and transport in bacteria by a conserved RNA structural element. RNA 9, 1084-1097, doi:10.1261/rna.5710303 (2003).
9 Warner, D. F., Savvi, S., Mizrahi, V. & Dawes, S. S. A riboswitch regulates expression of the coenzyme B12-independent methionine synthase in Mycobacterium tuberculosis: implications for differential methionine synthase function in strains H37Rv and CDC1551. J Bacteriol 189, 3655-3659, doi:10.1128/JB.00040-07 (2007).
10 Cullen, R. W. & Oace, S. M. Dietary pectin shortens the biologic half-life of vitamin B-12 in rats by increasing fecal and urinary losses. J Nutr 119, 1121-1127, doi:10.1093/jn/119.8.1121 (1989).
11 Snapper, S. B., Melton, R. E., Mustafa, S., Kieser, T. & Jacobs, W. R., Jr. Isolation and characterization of efficient plasmid transformation mutants of Mycobacterium smegmatis. Mol Microbiol 4, 1911-1919, doi:10.1111/j.1365-2958.1990.tb02040.x (1990).
12 Supply, P. et al. Genomic analysis of smooth tubercle bacilli provides insights into ancestry and pathoadaptation of Mycobacterium tuberculosis. Nat Genet 45, 172-179, doi:10.1038/ng.2517 (2013).
13 Garnier, T. et al. The complete genome sequence of Mycobacterium bovis. Proc Natl Acad Sci U S A 100, 7877-7882, doi:10.1073/pnas.1130426100 (2003).
14 Keating, L. A. et al. The pyruvate requirement of some members of the Mycobacterium tuberculosis complex is due to an inactive pyruvate kinase: implications for in vivo growth. Mol Microbiol 56, 163-174, doi:10.1111/j.1365-2958.2005.04524.x (2005).
15 Caminero, J. A. et al. Epidemiological evidence of the spread of a Mycobacterium tuberculosis strain of the Beijing genotype on Gran Canaria Island. Am J Respir Crit Care Med 164, 1165-1170, doi:10.1164/ajrccm.164.7.2101031 (2001).
16 Jackson, M. et al. Inactivation of the antigen 85C gene profoundly affects the mycolate content and alters the permeability of the Mycobacterium tuberculosis cell envelope. Mol Microbiol 31, 1573-1587, doi:10.1046/j.1365-2958.1999.01310.x (1999).
17 Gonzalo-Asensio, J. et al. PhoP: a missing piece in the intricate puzzle of Mycobacterium tuberculosis virulence. PLoS One 3, e3496, doi:10.1371/journal.pone.0003496 (2008).
18 Broset, E. et al. MTBVAC-Based TB-HIV Vaccine Is Safe, Elicits HIV-T Cell Responses, and Protects against Mycobacterium tuberculosis in Mice. Mol Ther Methods Clin Dev 13, 253-264, doi:10.1016/j.omtm.2019.01.014 (2019).
19 Broset, E. et al. Engineering a new vaccine platform for heterologous antigen delivery in live-attenuated Mycobacterium tuberculosis. Comput Struct Biotechnol J 19, 4273-4283, doi:10.1016/j.csbj.2021.07.035 (2021).
20 Song, H. & Niederweis, M. Functional expression of the Flp recombinase in Mycobacterium bovis BCG. Gene 399, 112-119, doi:10.1016/j.gene.2007.05.005 (2007).