F1000Research F1000Research 2046-1402 F1000 Research Limited London, UK 10.12688/f1000research.14239.2 Research Article Articles High monocyte to lymphocyte ratio is associated with impaired protection after subcutaneous administration of BCG in a mouse model of tuberculosis

[version 2; peer review: 2 approved]

 Zelmer Andrea Formal Analysis Funding Acquisition Investigation Methodology Supervision Visualization Writing – Original Draft Preparation Writing – Review & Editing https://orcid.org/0000-0001-9621-4406 a 1 Stockdale Lisa Investigation Writing – Review & Editing 1 Prabowo Satria A. Investigation Writing – Review & Editing https://orcid.org/0000-0002-3808-0638 1 Cia Felipe Investigation Writing – Review & Editing 1 Spink Natasha Investigation Writing – Review & Editing 1 Gibb Matthew Investigation Writing – Review & Editing https://orcid.org/0000-0002-6774-9943 1 2 Eddaoudi Ayad Investigation Methodology Writing – Review & Editing 3 Fletcher Helen A. Conceptualization Funding Acquisition Methodology Supervision Writing – Review & Editing https://orcid.org/0000-0002-0435-1006 1 London School of Hygiene and Tropical Medicine, Department of Immunology and Infection, Keppel Street, London, WC1E 7HT, UK Baylor Institute for Immunology Research, 3434 Live Oak Street, Dallas, Texas, 75204, USA UCL Great Ormond Street Institute of Child Health, 30 Guilford Street, London, WC1N 1EH, UK andrea.zelmer@lshtm.ac.uk

No competing interests were disclosed.

 27 6 2018 2018 7 296 20 6 2018 Copyright: © 2018 Zelmer A et al. 2018 This is an open access article distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Background: The only available tuberculosis (TB) vaccine, Bacillus Calmette-Guérin (BCG), has variable efficacy. New vaccines are therefore urgently needed. Why BCG fails is incompletely understood, and the tools used for early assessment of new vaccine candidates do not account for BCG variability. Taking correlates of risk of TB disease observed in human studies and back-translating them into mice to create models of BCG variability should allow novel vaccine candidates to be tested early in animal models that are more representative of the human populations most at risk. Furthermore, this could help to elucidate the immunological mechanisms leading to BCG failure. We have chosen the monocyte to lymphocyte (ML) ratio as a correlate of risk of TB disease and have back-translated this into a mouse model.

Methods: Four commercially available, inbred mouse strains were chosen. We investigated their baseline ML ratio by flow cytometry; extent of BCG-mediated protection from Mycobacterium tuberculosis infection by experimental challenge; vaccine-induced interferon gamma (IFNγ) response by ELISPOT assay; and tissue distribution of BCG by plating tissue homogenates.

Results: The ML ratio varied significantly between A/J, DBA/2, C57Bl/6 and 129S2 mice. A/J mice showed the highest BCG-mediated protection and lowest ML ratio, while 129S2 mice showed the lowest protection and higher ML ratio. We also found that A/J mice had a lower antigen specific IFNγ response than 129S2 mice. BCG tissue distribution appeared higher in A/J mice, although this was not statistically significant.

Conclusions: These results suggest that the ML ratio has an impact on BCG-mediated protection in mice, in alignment with observations from clinical studies. A/J and 129S2 mice may therefore be useful models of BCG vaccine variability for early TB vaccine testing. We speculate that failure of BCG to protect from TB disease is linked to poor tissue distribution in a ML high immune environment.

 Tuberculosis animal models BCG vaccine ML ratio mice London School of Hygiene and Tropical Medicine AthenaSWANCareerRe-EntryFellowship European Commission 643381 The work in this article was funded by a grant to HF from the European Commission HORIZON2020 program (TBVAC2020 grant no. 643381) and an Athena SWAN Career Re-entry award funded by the London School of Hygiene and Tropical Medicine to AZ. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Amendments from Version 1

This version has been revised to address the reviewers' comments. We have now included an analysis of the CD11b intermediate, CD11c+ cell population in the lung which should be enriched in alveolar macrophages (Fig 1E). We have also included a discussion about the main genetic differences between the mouse strains, and provided some clarification of experimental and technical details.

 Introduction

Tuberculosis (TB), caused by Mycobacterium tuberculosis ( Mtb) is the leading cause of death from a single infectious agent. Multidrug-resistant TB remains a public health crisis, and only one vaccine (the M. bovis-derived Bacillus Calmette-Guérin, BCG) is currently licensed for clinical use. To meet the sustainable development goal of ending the TB epidemic by 2030, new treatments and vaccines are both urgently needed.

BCG efficacy is highly variable 1, 2 . The reasons why BCG protects when it does and why it fails when it doesn’t are incompletely understood, but are crucial to the successful design and testing of new vaccines 3 . Furthermore, to be able to accurately assess vaccine candidates in models where BCG both protects and does not protect very early on in the vaccine development pipeline would be highly advantageous and de-risk failure in later stage clinical trials. There is however a lack of a broad range of tools that can collectively predict with some confidence whether a vaccine will be protective.

One of the tools for very early testing of TB vaccine candidates is a mouse model. However, in order to obtain meaningful information, new vaccine candidates should be tested in models that are clinically relevant and reflect the breadth and heterogeneity of immune environments and varying BCG efficacy found in human populations. We propose that this could be achieved by back-translating observations from clinical studies, such as correlates of risk of TB disease, into animal models.

A number of correlates of risk of TB disease have recently been identified, including transcriptomic mRNA signatures in blood, T cell activation, and monocyte to lymphocyte (ML) ratio 4– 8 , while others are being investigated (reviewed in 9). We have chosen the ML ratio to provide proof of principle that correlates of risk can be back-translated into the mouse to develop a model for vaccine testing that better reflects the populations most at risk of BCG failure.

In this study, we show that inbred, commercially available mouse strains have differing ML ratios, and that a high ML ratio is associated with lower BCG-mediated protection from experimental Mtb challenge. We further suggest that lack of BCG dissemination and/or persistence in ML high mice impairs protection.

 Methods Ethics statement

All animal work was carried out in accordance with the Animals (Scientific Procedures) Act 1986 under a license granted by the UK Home Office (PPL 70/8043), and approved locally by the London School of Hygiene and Tropical Medicine Animal Welfare and Ethics Review Body.

 Animals and experiment design

The following mouse strains were used for the experiments reported here (abbreviated names used throughout the manuscript are given in brackets): A/JOlaHsd (A/J); DBA/2OlaHsd (DBA/2); C57BL/6JOlaHsd (C57Bl/6); 129S2/SvHsd (129S2). Female mice were acquired from Envigo UK at 5–7 weeks of age. Animals were housed in specific pathogen-free individually vented cages with environmental enrichment (play tunnel and tapvei block), with 12 hours light / 12 hours dark cycles, at temperatures between 19° – 23°C and relative humidity of 45 – 65%. Mice were fed sterilized diet RM1 and filtered water ad libitum, and were allowed to acclimatize for at least 5 days before the start of any experimental procedure. Mice were allocated to cages as groups of 5 by technical staff not involved in experimental procedures or data analysis, and mice of the same strain were housed together. Each cage was allocated to a treatment in no particular order, but without formal randomisation. Animal welfare was assessed twice every day before and during the study.

Two independent experiments with separate primary outcomes were carried out to obtain the data described in this report. In Experiment 1, 5 naïve mice of each strain (20 mice total) were culled by anaesthetic overdose, and cardiac blood, lungs and spleens were collected for determination of the ML ratio by flow cytometric analysis (see below for details). The primary outcome for this experiment was the ML ratio ( Figure 1). In Experiment 2, parallel groups of mice (n=5 per group) were BCG-immunised or left untreated. These were allocated at the start of the study to either an immunogenicity group (IFNγ response and BCG dissemination; 10 mice per strain; 40 mice total; Figure 3 and Figure 4), or an Mtb challenge group (bacterial burden in lung; 10 mice per strain; 40 mice total; Figure 2). Six weeks after immunisation, mice were infected with Mtb or culled for isolation of cells from lung and spleen. In all, a total of 100 mice (25 per strain) was used to obtain all data presented here.

The ML ratio differs significantly in four different inbred mouse strains.

A Gating strategy for flow cytometric analysis. Cells from naïve animals were fixed, stained and data acquired as described in Materials and Methods. Cell debris was gated out by use of a FSC-SSC gate, followed by gating on single cells (FSC-H and FSC-W). A sequential gating strategy was then applied to determine the frequency of T cells (CD3 +), B cells (B220 +), neutrophils (CD11b + Ly6G +), monocytes/macrophages (CD11b+ CD11c low-int, red gate) and CD11b int CD11c + cells (blue gate) as a percentage of single cells. Plots shown are from a sample of a C57Bl/6 spleen. B– D The ML ratio was calculated by dividing the percentage of monocytes/macrophages by the sum of the percentages of B and T cells. ML ratio was analysed in blood ( B), lung ( C) and spleen ( D) of four different mouse strains. E Percentage of CD11b int CD11c + cells in the lung of four different mouse strains, likely to be enriched in alveolar macrophages. Each symbol represents one animal; box plots represent the median (middle line), 25 th to 75 th percentile (box) and minimum to maximum value (error bars). Data sets are presented in order of decreasing protection. p values were determined using ordinary ANOVA with Holm-Sidak test for multiple comparisons. Multiplicity adjusted p values are reported. A p value <0.05 was considered statistically significant.

 Protection from Mtb infection varies after BCG vaccination in ML high and low mouse strains.

A Mice of each strain were immunised s.c. with 2x10e5 CFU BCG Pasteur (open circles) or left untreated (filled circle), and infected i.n. with 30 CFU Mtb Erdmann 6 weeks later. Bacterial burden in the lungs of all animals was enumerated 7 weeks after challenge. Each symbol represents one animal. Bacterial numbers are given as log10 CFU per whole organ. No bacteria were detected in one of the C57Bl/6 samples; this value was set to zero. Δ indicates the difference in bacterial burden between naïve and BCG immunised mice. p values were determined by Kruskal-Wallis test with Dunn’s post-test for multiple comparisons between naïve groups, and multiplicity adjusted p values are reported. Individual Mann-Whitney tests to compare naïve with BCG immunised groups of each mouse strain. p<0.05 was considered statistically significant. Error bars represent the median and interquartile range. LoD: limit of detection. B ML ratio in comparison to protection in the lung of ML high and low mice. The median ML ratio in lung (red; as in Figure 1) is plotted together with protection expressed as the median % decrease in CFU (blue; as in A).

 Antigen-specific interferon gamma response varies in splenocytes of ML high and ML low mouse strains.

Splenocytes were isolated from BCG-immunised or control mice at the time of Mtb challenge (6 weeks after immunisation), and restimulated with PPD. The number of IFNγ producing cells were enumerated using an ELISPOT assay (presented as spot forming units [SFU] per 10e6 cells). Non-specific background measured in unstimulated duplicate wells was removed. Data sets are presented in order of decreasing protection. Each symbol represents one animal. p values were determined by using one-way ANOVA and Dunn’s post test for multiple comparisons. Error bars represent the median and interquartile range. Multiplicity adjusted p values are reported. p<0.05 was considered statistically significant.

 BCG distribution varies between lung and spleen after vaccination of ML high and ML low mice.

Approximately half of each lung ( A) and spleen ( B) was homogenised at the time of Mtb challenge (6 weeks after immunisation) and plated on 7H11 agar plates to determine viable BCG bacteria. Total CFU per organ are reported. Data sets are presented in order of decreasing protection. Each symbol represents one animal. In some instances, bacteria could not be detected and values were set to zero. p values were determined by using Kruskal-Wallis test and Dunn’s post test for multiple comparisons. Error bars represent the median and interquartile range. Multiplicity adjusted p values are reported. p<0.05 was considered statistically significant.

The study was not blinded.

 Vaccination

The BCG Pasteur strain was obtained from Aeras (Rockville, MD, USA) as frozen aliquots. These were stored at -80°C until needed. BCG was then thawed at room temperature and diluted to a final concentration of 2×10 6 CFU/ml in physiological saline solution for irrigation (Baxter Healthcare, Newbury, UK). Each animal received a subcutaneous injection of 100 μl BCG containing 2×10 5 CFU BCG (vaccinated groups). BCG dose was confirmed by plating of an aliquot of the prepared vaccine suspension on 7H11 agar plates. Colonies were counted after 12–14 days of incubation at 37°C. Animals were then rested for 6 weeks before either infection with Mtb, or sacrifice for cell isolation.

 BCG enumeration in tissues

Approximately half of each spleen and the right-hand side lobes of the lungs from mice in the immunogenicity group were used to determine the number of viable BCG bacteria in each organ six weeks after vaccination. Tissues were removed aseptically in a microbiological safety cabinet and placed in sterile 2 ml screwcap vials containing 500 μl PBS + 0.05% Tween80 and Precellys 1.4mm ceramic beads (CK14; Peqlab, Sarisbury Green, UK). A Precellys 24 homogeniser (Peqlab) was used to homogenise tissues for 15 s at 5000 rpm before plating. Each entire homogenate was plated onto two 7H11 agar plates containing 10% OADC supplement (Yorlab, York, UK) and 0.5% glycerol. Colonies were counted after 3 weeks of incubation at 37°C.

 Infection with <italic toggle="yes">M. tuberculosis</italic>

Mice were infected intranasally with M. tuberculosis Erdman (BEI Resources, Manassas, VA, USA) 6 weeks after BCG immunization and kept in isolators under CL-3 containment. Frozen aliquots of Mtb Erdman were thawed at room temperature, and diluted in saline. Mice were anaesthetized by an intraperitoneal injection of a combination of Ketamine (50 mg/kg; Ketalar Pfizer Itd, Kent, UK) and Xylazine (10 mg/kg; Rompun; Berkshire, UK) in saline. Each animal then received 50 μl of the inoculum, estimated to contain 30 CFU. The number of bacteria in the inoculum was confirmed by plating aliquots on 7H11 agar plates containing 10 % OADC and 0.5 % glycerol.

Seven weeks after infection, animals were killed by cervical dislocation. Lungs and spleens were removed aseptically and homogenized by mechanical disruption in sterile PBS, using the plunger of a 5 ml syringe and a 100 μm cell strainer. A series of 10-fold dilutions of tissue homogenates in PBS with 0.05 % Tween 80 were plated onto 7H11 agar plates with 10 % OADC supplement and 0.5 % glycerol. Plates were incubated at 37 °C and colonies counted after 3 weeks.

Protection is expressed by the percent reduction of median CFU/lung in the BCG group compared to the control group for each mouse strain.

 Flow cytometry

Single cell suspensions from spleens were prepared in RPMI-1640 media (Sigma-Aldrich, Dorset, UK) containing 10% heat-inactivated FBS (Labtech International Ltd, Uckfield, UK) and 2 mM L-Glutamine (Fisher Scientific, Loughborough, UK) as soon as possible after sacrifice. Spleens were mechanically disrupted by mashing through a 100μm cell strainer using the rubber end of the plunger from a 5ml syringe. Lungs were collected into RPMI-1640 media without FBS and cut into small pieces of approx. 2mm 3 before incubation with 0.5mg/ml Liberase TL (Sigma-Aldrich) and 10μg/ml DNAse (Sigma-Aldrich). The enzyme reaction was stopped by adding an equal amount of media containing 10% FBS and the tissue was mashed as above to obtain single cells. Cells were fixed and red blood cells lysed by adding lyse-fix solution (PhosFlow; Becton Dickinson, Oxford, UK). After fixing, cell suspensions were made up in PBS + 1% FBS.

Approximately 10 6 cells were stained with the following antibody cocktail in BD Brilliant Stain buffer as per manufacturer’s instructions (Becton Dickinson): CD3-APC/Cy7 (clone 17A2, 1:80), B220-BV510 (clone RA3-6B2, 1:40), Ly6G-BV711 (clone 1A8, 1:40), CD11b-BV650 (clone M1/70, 1:60), CD11c-BV605 (clone N418, 1:40). All antibodies were purchased from Biolegend (via Fisher Scientific).

Cells from each tissue from one mouse per strain were used as fluorescence minus one (FMO) controls. These were stained with the antibody cocktail as described above, but without one of the antibodies. This was done for each antibody. FMO controls were used to guide gating.

OneComp beads (eBioscience via Fisher Scientific, Loughborough, UK) were stained with single antibodies as per manufacturer’s instructions and used to calculate compensation using the automatic function in FlowJo version 10.4. Compensation matrices were manually checked and adjusted where necessary.

Analysis was carried out using FlowJo version 10.4.

The ML ratio was calculated by dividing the percentage of monocytes/macrophages (Mono_Mac) (B220 - CD3 - Ly6G - CD11b + CD11c lo-int) relative to single cells by the percentage of B cells (B220 +) and T cells (CD3 +): Mono_Mac (% of single cells) / [B cells (% of single cells) + T cells (% of single cells)]

 IFN-γ ELISPOT

To quantify IFN-γ secreting antigen-specific splenocytes, single cell suspensions were prepared by mechanical disruption of the remaining spleen samples from the immunogenicity group through a 100μm cell strainer as soon as possible after sacrifice. After lysis of red blood cells, single cell suspensions were made up in RPMI-1640 media containing 10% heat-inactivated FBS and 2 mM L-Glutamine. 96-well microtiter ELISPOT plates (MAIPS4510, Millipore, Watford, UK) were coated with 10 µg/ml rat anti-mouse IFN-γ (clone AN18, Mabtech, Nacka Strand, Sweden). Free binding sites were blocked with RMPI-1640 supplemented with 10% heat-inactivated FBS and 2 mM L-Glutamine. 2×10 5 of total splenocytes were added and incubated in duplicate with PPD (10 µg/ml), supplemented RPMI as a negative control, or Phorbol myristate acetate (PMA) (0.1 µg/ml, Sigma-Aldrich) and Phytohemagglutinin (PHA) (1 µg/ml, Sigma-Aldrich) as a positive control. After overnight incubation at 37°C in 5% CO 2, IFN-γ was detected with 1 µg/ml biotin labelled rat anti-mouse antibody (clone R4-6A2, Mabtech) and 1 µg/ml alkaline phosphatase-conjugated streptavidin (Mabtech). The enzyme reaction was developed with BCIP/NBT substrate (5-Bromo-4-chloro-3-indolyl phosphate/Nitro blue tetrazolium) (MP Biochemicals, UK) and stopped by washing the plates with tap water when individual spots could be visually detected (up to 5min). ELISPOT plates were analysed using an automatic plate reader. IFN-γ-specific cells are expressed as number of spot-forming units (SFU) per million spleen cells after non-specific background was subtracted using negative control wells.

 Statistical analysis

Groups of animals were compared, and a p value of <0.05 was considered statistically significant. Statistical analysis was carried out using GraphPad Prism version 6. The specific test used for each analysis is described in the figure legends.

 Results The ML ratio differs significantly in four inbred mouse strains

Four commercially available inbred mouse strains were chosen to investigate the impact of the baseline ML ratio on BCG vaccine efficacy. A/J, DBA/2, C57Bl/6 and 129S2 mice were selected to represent varying monocyte frequencies and a range of BCG-mediated protection based on available data (Jax Phenome Database https://phenome.jax.org; 10, 11 ). To allow the direct comparison of the ML ratio between strains, animals were age and sex matched and all samples were processed at the same time, stained with aliquots of the same antibody cocktail, and data acquired as a batch. All animals used in this experiment were included in the analysis.

We found significant differences between mouse strains in their baseline ML ratio in blood, spleen, and lung ( Figure 1). Across all tissues, A/J mice showed the lowest ML ratio ( Figure 1 B–D). In spleen and blood, the ML ratio of DBA/2, C57Bl/6 and 129S2 mice was higher than A/J, but only minor differences were found between those three strains ( Figure 1 B and D). In the lung, differences between all strains were more apparent ( Figure 1C). 129S2 mice had the highest ML ratio, with DBA/2 and C57Bl/6 showing intermediate levels. These differences in ML ratio were mostly driven by differing frequencies in monocytes/macrophages, with minor differences in B and T cell frequencies between mouse strains ( Figure. S1–S3).

This file contains the data underlying the analysis of the ML ratio and cell subsets in blood, spleen, and lung shown in Figure 1, Figure 2B, Figure S1, Figure S2, and Figure S3.

To assess the impact of these differences in ML ratio on BCG efficacy in our model, we infected mice of all strains with virulent Mtb Erdmann.

 Mouse strains with varying ML ratios are differentially protected from <italic toggle="yes">Mtb</italic> infection by BCG

Mice of each strain were vaccinated s.c. with BCG, or left untreated, and infected with Mtb Erdmann six weeks later. The extent of protection by BCG was measured by determining the bacterial burden in lung and spleen of vaccinated and control mice ( Figure 2). All mice were immunised with the same inoculum of BCG and infected using the same inoculum of Mtb, so that results of the different strains are directly comparable. All animals were included in the analysis and no unexpected adverse events were observed. No more than 10% weight loss was recorded over the duration of the experiment. We found varying degrees of protection ( Figure 2A). A/J mice were the most protected with a difference in median CFU (Δ) between vaccinated and control groups of 1.41 log 10, followed by DBA/2 (Δ0.81 log 10) and C57Bl/6 (Δ0.54 log 10), and 129S2 mice were the least protected (Δ0.37 log 10). Interestingly, the mouse strain with the highest ML ratio in the lung showed the lowest protection (129S2), and the strain with the lowest ML ratio showed the highest protection (A/J; Figure 2B). Some differences were also observed in the innate susceptibility between strains, most notably C57Bl/6 mice were significantly less susceptible than DBA/2 (Δ1.16 log 10 CFU) or A/J mice (Δ1.41 log 10 CFU). Alveolar macrophages play an important role in the defence against respiratory pathogens, however they are not included in the gate used here to quantify monocytes/macrophages ( Figure 2A). While we did not include specific markers for this cell population, they carry a CD11b intermediate (CD11bint), CD11c+ phenotype. We found that the frequency of CD11bint CD11c+ cells was significantly higher in C57Bl/6 mice compared to the other mouse strains ( Figure 2E).

Next, to investigate whether differential BCG-mediated protection from Mtb infection was associated with a differential immune response, we carried out an IFNγ ELISPOT assay.

This files contains the raw data as colony forming units (CFU) of the graphs in Figure 2.

 A differential antigen-specific IFNγ response is associated with differential protection

To assess the immune response to mycobacterial antigens in the differentially protected mouse strains, splenocytes from mice vaccinated with BCG were stimulated with PPD, and the number of IFNγ-producing cells was measured by ELISPOT ( Figure 3). IFNγ is a cytokine associated with control of infection in mice 12 and reduced risk of TB disease in humans 6 . We found that there were significantly more antigen-specific IFNγ producing splenocytes in 129S2 mice than in A/J or DBA/2 mice, both of which were better protected by BCG from Mtb challenge than 129S2 mice.

Since BCG is a live mycobacterium, and monocytes and macrophages are the natural host cells for these organisms, it was plausible that a high ML ratio and high IFNγ production in 129S2 mice impair BCG survival and dissemination to and/or persistence in tissues, and thus prevent the formation of a protective immune response in the lung.

This file contains the raw data for the graph shown in Figure 3

BCG dissemination to lung and spleen is impaired in ML high mice

To determine the extent to which BCG is found in lungs and spleens after s.c. administration, we plated organ homogenates onto 7H11 agar plates and enumerated viable BCG after 3 weeks of culture ( Figure 4). We generally found low numbers of bacteria per organ, especially in the lung, where the bacterial count did not exceed 80 CFU/lung in any sample. Moderate differences were detected between mouse strains. Most notably, the most protected ML low A/J mice had the highest number of BCG in the lung, while no BCG was detectable in the lung of the least protected ML high 129S2 mice.

These data led us to speculate that certain host immune-environments, such as a high ML ratio, can influence BCG dissemination and/or persistence, which in turn could impact BCG efficacy.

The data presented here are from one experiment as detailed in the Methods section. To ensure reproducibility, we have endeavoured to include as much experimental detail as possible to allow others to repeat the study and confirm the results. Groups of animals of all strains received the same treatments at the same time to minimise technical variation. We are therefore confident that the major differences between strains are true biological differences.

This file contains the raw data as colony forming units (CFU) for the graph shown in Figure 4.

 Discussion

The varying efficacy of BCG has long been a concern for the control of the TB epidemic. Numerous factors have been discussed as contributors to this, such as population differences, different BCG strains, vaccination schedules, exposure to environmental mycobacteria, co-infections with viruses and/or parasites, and geographical location 1, 3, 13– 16 . While BCG works in some populations and protects children from pulmonary and extra-pulmonary disease, new vaccines are urgently needed for the populations where BCG fails to protect, such as adolescents and adults in endemic areas 1, 2 . The incomplete understanding of the mechanisms behind BCG failure, together with limited tools for early assessment of vaccine efficacy, are hampering the development of new vaccines. In order to allow early vaccine testing in a model that better represents the populations most at risk of TB disease, we have taken a correlate of risk observed in human studies and back-translated it into a mouse model. The monocyte to lymphocyte (ML) ratio is a non-specific marker of inflammation and has been shown to be associated with risk of TB disease in several populations, including pregnant HIV infected women, people starting anti-retroviral therapy, BCG vaccinated infants and latently Mtb infected adolescents 7, 8, 17– 19 . Scriba et al. show a cascade of inflammatory events as adolescents progress towards disease, which appears to start with increased Type I/II IFN signalling, followed by an increase in monocytes and a decrease in lymphocytes 19 . An increased risk of TB disease in BCG vaccinated individuals with a high ML ratio indicates that BCG is failing in those populations.

In order to provide proof of principle that the ML ratio, and observations from human studies more generally, can be back-translated to develop more clinically relevant mouse models, we chose to use commercially available, genetically tractable inbred mouse strains. We could show that these have highly varying ML ratios in blood, spleen, and lung, and that these differences, particularly in the lung, were associated with varying BCG-mediated protection from Mtb infection. Using mouse strains with genetically different backgrounds means that there are likely other confounding factors that could impact vaccine efficacy. This is exemplified by the DBA/2 mice, which are well protected despite having a relatively high ML ratio across all tissues. This could be due to the fact that they have a higher neutrophil frequency than some of the other strains (data not shown). Although confounding factors may be at play here, this would also be the case in a clinical setting since an altered ML ratio can have different underlying causes. An increased ML ratio is an indicator for inflammation, likely driven by monocytosis during inflammation. Some factors that are thought to influence BCG efficacy also drive inflammation, such as co-infection with viruses or parasites, malnutrition, or exposure to environmental mycobacteria. However, as manipulation of ML ratio in vitro impacts on ability to control mycobacterial growth, it is likely that ML ratio itself is contributing to TB risk, independent of the factor driving inflammation 20 . It is also possible that host genetic factors determine baseline inflammation and ML ratio 21 . If this is the case in TB endemic populations, inbred mouse strains with naturally varying ML ratios may be a useful model for BCG vaccine variability.

There are genetic differences between the inbred mouse strains used in this study, some of which have been well described, while others are likely to be unknown. Among the most well described genetic factors is the mutation of the interleukin-3 (IL-3) receptor alpha subunit gene in A/J mice, which leads to impaired IL-3 signalling and reduced proliferation of haematopoietic stem cells and their differentiation into myeloid precursors. This is likely to play a role in the low ML ratio observed in these mice 22, 23 . A/J mice additionally carry a defect in the gene for Naip5, an intracellular pattern recognition receptor involved in control of proliferation of intracellular pathogens 24, 25 . It may thus have an effect on BCG persistence and dissemination. C57Bl/6 mice carry a polymorphism in the slc11a1 (formerly nramp1) gene, which has been shown to determine susceptibility or resistance to intravenous infection with BCG as defined by BCG load in the spleen 26 . It does not seem to determine susceptibility to primary Mtb infection in mice, and the effect on BCG efficacy is not well studied 27 . In addition, different inbred mouse strains carry different MHC haplotypes. For the strains used in this study, they are as follows: 129S2, H2 b; C57Bl/6, H2 b; A/J, H2 a; DBA/2, H2 d. Each one confers the ability to recognise a slightly different range of epitopes, although there is overlap. This may also impact BCG- mediated immunity.

It is currently unclear what the relative contributions are of genetic and environmental factors to a change in ML ratio 21 . Does a naturally high ML ratio in an individual make them more prone to TB disease, or do other factors including Mtb itself drive up the ML ratio, causing BCG to fail? Does a naturally high ML ratio in an individual exacerbate an inflammatory response initiated by environmental factors? There is striking heterogeneity in the immune response to BCG in infants vaccinated at birth, including in the ML ratio and cytokine responses 17 . This points towards an involvement of host factors, as the influence of environmental factors would be limited so early in life. On the other hand, Scriba et al. observe increases in inflammation and ML ratio over time in adolescents as an individual progresses towards active TB disease and diagnosis 19 . This may point towards the pathogen itself as the cause of an altered immune environment.

Interestingly, we found a higher number of antigen-specific IFNγ-producing splenocytes in unprotected 129S2 mice compared to protected A/J mice at the time of Mtb challenge. This seemed counterintuitive at first, but is in line with recent findings showing that a sub-population of BCG-vaccinated infants with an increased ML ratio showed increased T cell-mediated cytokine production, and this was associated with increased risk of TB disease 17 . High monocyte count and pro-inflammatory cytokines are indicators for inflammation. It is possible that BCG efficacy is impaired in inflammatory immune environments, which may be caused by a variety of factors, such as chronic viral infection, age, malnutrition, or prolonged contact with environmental mycobacteria, all factors thought to impact BCG efficacy. Ultimately, BCG is administered into host immune environments with a heterogeneity that is not fully characterised or understood, yet may have implications for its efficacy.

The exact mechanisms behind an altered BCG vaccine efficacy in a ML high immune environment are currently unknown. Our data indicate that impaired dissemination and/or persistence of BCG may play a role in initiating an effective local immune response. A study by Kaveh and colleagues has shown that BCG persists for up to 16 months in lymph nodes and spleen of Balb/c mice and that persistence of live bacilli contributes to optimal protection from M. bovis challenge 28 . Direct investigation of this phenomenon in humans is not possible; there are however reports that indicate that BCG can disseminate and persist for years in a variety of organs in humans 29, 30 . It is therefore plausible that an immune environment that impairs survival of BCG also impairs vaccine efficacy, although we did not determine in this study whether the absence of live bacilli in lungs of 129S2 mice is due to lack of dissemination or lack of persistence.

The absolute numbers of BCG might be under-estimated in this study as we only enumerated bacilli that were readily growing on standard 7H11 agar plates. Use of PCR methods or culture with resuscitation promoting factors may increase yield and allow for detection of dormant bacilli and their contribution to protection 31, 32 .

If impaired dissemination and/or persistence in certain immune environments is indeed reducing BCG vaccine efficacy, it is reasonable to assume that other live mycobacterial vaccines would have similar limitations. Variability in the ML ratio and more general heterogeneity in the host immune environment should thus be taken into account when designing new live TB vaccines or vaccines based on boosting BCG.

To circumvent potentially impaired dissemination of parenterally administered live vaccines, mucosal administration of live vaccines could be considered. Recent evidence shows that mucosal administration of BCG confers increased protection in animal models, for example by inducing homing of T cells into protected mucosal niches 33– 35 . This agrees with the notion that BCG needs to be present in the lungs to initiate a protective immune response in this organ. Furthermore, subunit vaccines should not show decreased efficacy in ML high immune environments if our hypothesis is correct. Both these possibilities will need careful investigation to shed light on the mechanisms behind BCG failure.

We did not investigate the phenotype of monocytes/macrophages or T cells in this study, although it is likely that this plays a role. In particular, inflammatory monocytes might lead to killing of BCG. It will be interesting in future to investigate this, as well as comparing the ML ratio and immune cell phenotypes before and after BCG vaccination in the different mouse strains and locations.

We found a significantly higher proportion of CD11b int CD11c + cells in the lungs of C57Bl/6 mice, which are likely to be enriched in alveolar macrophages. Interestingly, C57Bl/6 mice are also the strain least susceptible to Mtb infection without BCG vaccination. While the frequency of this population does not seem to be associated with BCG-mediated protection after subcutaneous administration of the vaccine, it could play a role after intranasal vaccination. Alveolar macrophages are located on the luminal side of the epithelium, therefore encountering antigens and pathogens taken up from environmental air. The monocyte/macrophage population included in our ML ratio on the other hand is likely to be of interstitial nature, residing on the basal side of the epithelium and encountering antigens from the blood stream.

While our data on the susceptibility of naïve mouse strains to Mtb infection agrees with other studies 36 , we found some discrepancies between the data presented here and previously published data on BCG-mediated protection. A study using DBA/2 mice has shown that these mice are not protected by s.c. administration of BCG, but control infection with Mtb better when vaccinated i.n. 11 This is in contrast to our finding that DBA/2 mice are protected from Mtb challenge after s.c. vaccination. However, there are some differences between these two studies that may explain the different findings: a different source of mice was used; protection was assessed at 8 weeks vs 6 weeks after vaccination; the BCG strain and Mtb strain differed; and the Mtb dose used for challenge was three times higher in the study by Aguilo and colleagues compared to ours. Another report detailing the susceptibility of different mouse strains to Mtb challenge shows that A/J mice are notably less protected than C57Bl/6 mice 10 , while in our hands the opposite is the case. There are important differences between the two studies, including the age and sex of mice used, the time point after vaccination at which protection was measured, and the use of different BCG and Mtb strains. In combination, these factors may lead to a differing extent of protection. With particular relevance to this study, the immune cell frequency between male and female animals differs in A/J mice. For example, dataset Donahue5 on the Jackson Laboratory Mouse Phenome Database indicates a higher percentage of monocytes in blood for males compared to females, while T cells are less frequent in males, leading to a higher ML ratio in males 37 .

In conclusion, the back-translation of a correlate of risk of TB disease, the ML ratio, into an animal model is a step forward for better tools to study the immune mechanisms behind BCG vaccine failure and to test vaccines in a model more relevant to populations most at risk. Hopefully other risk factors can also be back-translated in the future to obtain a more complete picture of vaccine efficacy at very early stages of testing. We would encourage the use of more diverse mouse models in TB vaccine testing and development to reflect the heterogeneity of immune environments found in human populations.

 Data availability

The data referenced by this article are under copyright with the following copyright statement: Copyright: � 2018 Zelmer A et al.

Data associated with the article are available under the terms of the Creative Commons Zero "No rights reserved" data waiver (CC0 1.0 Public domain dedication).

Dataset 1: Flow cytometry raw data – This file contains the data underlying the analysis of the ML ratio and cell subsets in blood, spleen, and lung shown in Figure 1, Figure 2B, Figure S1, Figure S2, and Figure S3. 10.5256/f1000research.14239.d207942 38

Dataset 2: Mtb challenge CFU raw data – This files contains the raw data as colony forming units (CFU) of the graphs in Figure 2. 10.5256/f1000research.14239.d197118 39

Dataset 3: ELISPOT raw data – This file contains the raw data for the graph shown in Figure 3. 10.5256/f1000research.14239.d197136 40

Dataset 4: BCG distribution raw data – This file contains the raw data as colony forming units (CFU) for the graph shown in Figure 4. 10.5256/f1000research.14239.d197137 41

Acknowledgements

We would like to thank Christopher Sassetti for helpful discussions and selecting the mouse strains, and Stephanie Canning for expert technical assistance with flow cytometry studies.

 Supplementary material

Figure S1: The frequency of T cells varies in tissues of ML high and ML low mouse strains. T cell (CD3 +) frequencies were determined in blood ( A), lung ( B) and spleen ( C) of four different mouse strains using the gating strategy described in Figure 2A. The ML ratio was calculated using these data. Each symbol represents one animal; box plots represent the median (middle line), 25 th to 75 th percentile (box) and minimum to maximum value (error bars). Data sets are presented in order of decreasing protection. p values were determined using ordinary ANOVA with Holm-Sidak test for multiple comparisons. Multiplicity adjusted p values are reported. A p value <0.05 was considered statistically significant.

Click here to access the data.

Figure S2: The frequency of B cells varies in tissues of ML high and ML low mouse strains. B cell frequencies (B220 +) were determined in blood ( A), lung ( B) and spleen ( C) of four different mouse strains using the gating strategy described in Figure 2A. The ML ratio was calculated using these data. Each symbol represents one animal; box plots represent the median (middle line), 25 th to 75 th percentile (box) and minimum to maximum value (error bars). Data sets are presented in order of decreasing protection. p values were determined using ordinary ANOVA with Holm-Sidak test for multiple comparisons. Multiplicity adjusted p values are reported. A p value <0.05 was considered statistically significant.

Click here to access the data.

Figure S3: The frequency of monocytes/macrophages varies in tissues of ML high and ML low mouse strains. Monocyte/macrophage (CD11b + CD11c low-int) frequencies were determined in blood ( A), lung ( B) and spleen ( C) of four different mouse strains using the gating strategy described in Figure 2A. The ML ratio was calculated using these data. Each symbol represents one animal; box plots represent the median (middle line), 25 th to 75 th percentile (box) and minimum to maximum value (error bars). Data sets are presented in order of decreasing protection. p values were determined using ordinary ANOVA with Holm-Sidak test for multiple comparisons. Multiplicity adjusted p values are reported. A p value <0.05 was considered statistically significant.

Click here to access the data.

 Mangtani P Abubakar I Ariti C : Protection by BCG vaccine against tuberculosis: a systematic review of randomized controlled trials. Clin Infect Dis. 2014;58(4):470–80. 24336911 10.1093/cid/cit790 Colditz GA Brewer TF Berkey CS : Efficacy of BCG vaccine in the prevention of tuberculosis. Meta-analysis of the published literature. JAMA. 1994;271(9):698–702. 8309034 10.1001/jama.1994.03510330076038 Dockrell HM Smith SG : What Have We Learnt about BCG Vaccination in the Last 20 Years? Front Immunol. 2017;8:1134. 28955344 10.3389/fimmu.2017.01134 5601272 Zak DE Penn-Nicholson A Scriba TJ : A blood RNA signature for tuberculosis disease risk: a prospective cohort study. Lancet. 2016;387(10035):2312–22. 27017310 10.1016/S0140-6736(15)01316-1 5392204 Berry MP Graham CM McNab FW : An interferon-inducible neutrophil-driven blood transcriptional signature in human tuberculosis. Nature. 2010;466(7309):973–7. 20725040 10.1038/nature09247 3492754 Fletcher HA Snowden MA Landry B : T-cell activation is an immune correlate of risk in BCG vaccinated infants. Nat Commun. 2016;7: 11290. 27068708 10.1038/ncomms11290 4832066 Naranbhai V Hill AV Abdool Karim SS : Ratio of monocytes to lymphocytes in peripheral blood identifies adults at risk of incident tuberculosis among HIV-infected adults initiating antiretroviral therapy. J Infect Dis. 2014;209(4):500–9. 24041796 10.1093/infdis/jit494 3903371 Naranbhai V Moodley D Chipato T : The association between the ratio of monocytes: lymphocytes and risk of tuberculosis among HIV-infected postpartum women. J Acquir Immune Defic Syndr. 2014;67(5):573–5. 25247435 10.1097/QAI.0000000000000353 4229408 Petruccioli E Scriba TJ Petrone L : Correlates of tuberculosis risk: predictive biomarkers for progression to active tuberculosis. Eur Respir J. 2016;48(6):1751–63. 27836953 10.1183/13993003.01012-2016 5898936 Smith CM Proulx MK Olive AJ : Tuberculosis Susceptibility and Vaccine Protection Are Independently Controlled by Host Genotype. mBio. 2016;7(5):pii: e01516–6. 27651361 10.1128/mBio.01516-16 5030360 Aguilo N Alvarez-Arguedas S Uranga S : Pulmonary but Not Subcutaneous Delivery of BCG Vaccine Confers Protection to Tuberculosis-Susceptible Mice by an Interleukin 17-Dependent Mechanism. J Infect Dis. 2016;213(5):831–9. 26494773 10.1093/infdis/jiv503 Cooper AM Dalton DK Stewart TA : Disseminated tuberculosis in interferon gamma gene-disrupted mice. J Exp Med. 1993;178(6):2243–7. 8245795 10.1084/jem.178.6.2243 2191280 Brosch R Gordon SV Garnier T : Genome plasticity of BCG and impact on vaccine efficacy. Proc Natl Acad Sci U S A. 2007;104(13):5596–601. 17372194 10.1073/pnas.0700869104 1838518 Muller J Matsumiya M Snowden MA : Cytomegalovirus infection is a risk factor for TB disease in Infants. bioRxiv. 2017;222646. 10.1101/222646 Lalor MK Ben-Smith A Gorak-Stolinska P : Population differences in immune responses to Bacille Calmette-Guérin vaccination in infancy. J Infect Dis. 2009;199(6):795–800. 19434928 10.1086/597069 3276835 Chatterjee S Natman TB : Helminth-induced immune regulation: implications for immune responses to tuberculosis. PLoS Pathog. 2015;11(1):e1004582. 25632943 10.1371/journal.ppat.1004582 4310592 Fletcher HA Filali-Mouhim A Nemes E : Human newborn bacille Calmette-Guérin vaccination and risk of tuberculosis disease: a case-control study. BMC Med. 2016;14(1):76. 27183822 10.1186/s12916-016-0617-3 4869393 Naranbhai V Kim S Fletcher H : The association between the ratio of monocytes:lymphocytes at age 3 months and risk of tuberculosis (TB) in the first two years of life. BMC Med. 2014;12(1):120. 25034889 10.1186/s12916-014-0120-7 4223414 Scriba TJ Penn-Nicholson A Shankar S : Sequential inflammatory processes define human progression from M. tuberculosis infection to tuberculosis disease. PLoS Pathog. 2017;13(11):e1006687. 29145483 10.1371/journal.ppat.1006687 5689825 Naranbhai V Fletcher HA Tanner R : Distinct Transcriptional and Anti-Mycobacterial Profiles of Peripheral Blood Monocytes Dependent on the Ratio of Monocytes: Lymphocytes. EBioMedicine. 2015;2(11):1619–26. 28958495 10.1016/j.ebiom.2015.09.027 Lin BD Willemsen G Fedko IO : Heritability and GWAS Studies for Monocyte-Lymphocyte Ratio. Twin Res Hum Genet. 2017;20(2):97–107. 28193307 10.1017/thg.2017.3 Ichihara M Hara T Takagi M : Impaired interleukin-3 (IL-3) response of the A/J mouse is caused by a branch point deletion in the IL-3 receptor alpha subunit gene. EMBO J. 1995;14(5):939–950. 7889941 398166 Morris CF Salisbury J Kobayashi M : Interleukin 3 alone does not support the proliferation of bone marrow cells from A/J mice: a novel system for studying the synergistic activities of IL-3. Br J Haematol. 1990;74(2):131–137. 2138495 10.1111/j.1365-2141.1990.tb02555.x Kang TJ Lee GS Kim SK : Comparison of two mice strains, A/J and C57BL/6, in caspase-1 activity and IL-1beta secretion of macrophage to Mycobacterium leprae infection. Mediators Inflamm. 2010;2010:1–5, 708713. 20671924 10.1155/2010/708713 2910506 Lamkanfi M Amer A Kanneganti TD : The Nod-like receptor family member Naip5/Birc1e restricts Legionella pneumophila growth independently of caspase-1 activation. J Immunol. 2007;178(12):8022–8027. 17548639 10.4049/jimmunol.178.12.8022 Sellers RS Clifford CB Treuting PM : Immunological variation between inbred laboratory mouse strains: points to consider in phenotyping genetically immunomodified mice. Vet Pathol. 2012;49(1):32–43. 22135019 10.1177/0300985811429314 North RJ LaCourse R Ryan L : Consequence of Nramp1 deletion to Mycobacterium tuberculosis infection in mice. Infect Immun. 1999;67(11):5811–5814. 10531233 96959 Kaveh DA Garcia-Pelayo MC Hogarth PJ : Persistent BCG bacilli perpetuate CD4 T effector memory and optimal protection against tuberculosis. Vaccine. 2014;32(51):6911–8. 25444816 10.1016/j.vaccine.2014.10.041 Tajima Y Takagi R Nakajima T : An infant with asymptomatic hepatic granuloma probably caused by bacillus Calmette-Guérin (BCG) vaccination found incidentally at autopsy: a case report. Cases J. 2008;1(1):337. 19019255 10.1186/1757-1626-1-337 2615765 Gormsen H : On The Occurrence Of Epithelioid Cell Granulomas In The Organs Of Bcg-Vaccinated Human Beings. Acta Pathol Microbiol Scand Suppl. 1956;38(Suppl 111):117–20. 13372286 10.1111/j.1600-0463.1956.tb06761.x Loraine J Pu F Turapov O : Development of an In Vitro Assay for Detection of Drug-Induced Resuscitation-Promoting-Factor-Dependent Mycobacteria. Antimicrob Agents Chemother. 2016;60(10):6227–33. 27503641 10.1128/AAC.00518-16 5038329 Jin TH Qu T Raina A : Identification of Growth Promoting Effect of rBCG/BCG Culture Supernatant and Its Potential Applications. WJV. 2013;2013(02):32–8. 10.4236/wjv.2013.32006 Lai R Afkhami S Haddadi S : Mucosal immunity and novel tuberculosis vaccine strategies: route of immunisation-determined T-cell homing to restricted lung mucosal compartments. Eur Respir Rev. 2015;24(136):356–60. 26028646 10.1183/16000617.00002515 Perdomo C Zedler U Kühl AA : Mucosal BCG Vaccination Induces Protective Lung-Resident Memory T Cell Populations against Tuberculosis. mBio. 2016;7(6):pii: e01686–16. 27879332 10.1128/mBio.01686-16 5120139 Jeyanathan M Afkhami S Khera A : CXCR3 Signaling Is Required for Restricted Homing of Parenteral Tuberculosis Vaccine-Induced T Cells to Both the Lung Parenchyma and Airway. J Immunol. 2017;199(7):2555–69. 28827285 10.4049/jimmunol.1700382 Chackerian AA Behar SM : Susceptibility to Mycobacterium tuberculosis: lessons from inbred strains of mice. Tuberculosis (Edinb). 2003;83(5):279–85. 12972341 10.1016/S1472-9792(03)00017-9 Donahue L Morgan J : Donahue 5. Mouse Phenome Database. (Accessed: 9 February 2018). Reference Source Zelmer A Stockdale L Prabowo SA : Dataset 1 in: High monocyte to lymphocyte ratio is associated with impaired protection after subcutaneous administration of BCG in a mouse model of tuberculosis. F1000Research. 2018. Data Source Zelmer A Stockdale L Prabowo SA : Dataset 2 in: High monocyte to lymphocyte ratio is associated with impaired protection after subcutaneous administration of BCG in a mouse model of tuberculosis. F1000Research. 2018. Data Source Zelmer A Stockdale L Prabowo SA : Dataset 3 in: High monocyte to lymphocyte ratio is associated with impaired protection after subcutaneous administration of BCG in a mouse model of tuberculosis. F1000Research. 2018. Data Source Zelmer A Stockdale L Prabowo SA : Dataset 4 in: High monocyte to lymphocyte ratio is associated with impaired protection after subcutaneous administration of BCG in a mouse model of tuberculosis. F1000Research. 2018. Data Source 10.5256/f1000research.16467.r35522 Reviewer response for version 2 Aguilo Nacho 1 Referee Department of Microbiology, Preventive Medicine, and Public Health, University of Zaragoza, Zaragoza, Spain

Competing interests: No competing interests were disclosed.

 4 7 2018 Copyright: © 2018 Aguilo N 2018 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve

I am satisfied with the answers given to my comments and therefore I recommend the indexing of the study.

Is the work clearly and accurately presented and does it cite the current literature?

Yes

If applicable, is the statistical analysis and its interpretation appropriate?

Partly

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Yes

Are the conclusions drawn adequately supported by the results?

Partly

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

tuberculosis vaccines immunology

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

 10.5256/f1000research.15495.r32475 Reviewer response for version 1 Gutierrez Maximiliano G. 1 Referee Host-Pathogen Interactions in Tuberculosis Laboratory, Francis Crick Institute, London, USA

Competing interests: No competing interests were disclosed.

 20 4 2018 Copyright: © 2018 Gutierrez MG 2018 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve

The work from Zelmer and collaborators is very interesting, timely, well executed and clearly written. I anticipate a great interest in the tuberculosis vaccine field.

I have only minor comments and suggestions:

Regarding introduction and/or discussion, it would be important the authors discuss what are the main genetic differences between the 4 mouse strains selected for this study.

In Figure 4, it is right to state that colonies vary in tissues but only 2 tissues were analysed here so it would be more appropriate to state "between lung and spleen". Moreover, would it be possible to get the data in lymph nodes as well? Lymph nodes are critical for the responses to vaccination.

The variation in CFU/organ seem to be higher in animals that have high levels of protection in this study (e.g. A/J mice). Could the authors comment on this? Is it possible to know the ML ratio in the 2 mice showing high CFU/lung from Fig 4A?

Not sure if the authors have the samples available, but it would have been interesting to see the profiles of Type I/II IFN in tissue at least by qPCR. This information could be of great interest for the tuberculosis field.

The authors should discuss the limitations of a relatively low number of biological experiments in this work.

Is the work clearly and accurately presented and does it cite the current literature?

Yes

If applicable, is the statistical analysis and its interpretation appropriate?

I cannot comment. A qualified statistician is required.

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Yes

Are the conclusions drawn adequately supported by the results?

Yes

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

Zelmer Andrea LSHTM

Competing interests: No competing interests were disclosed.

 19 6 2018

Dear Dr Gutierrez

We would like to thank you for taking the time to review our manuscript, and for the thoughtful feedback provided. We address your individual comments below. 

The work from Zelmer and collaborators is very interesting, timely, well executed and clearly written. I anticipate a great interest in the tuberculosis vaccine field. I have only minor comments and suggestions:

Regarding introduction and/or discussion, it would be important the authors discuss what are the main genetic differences between the 4 mouse strainsselected for this study.

We have now included a paragraph in the Discussion to address this.

In Figure 4, it is right to state that colonies vary in tissues but only 2 tissues were analysed here so it would be more appropriate to state "between lung and spleen". Moreover, would it be possible to get the data in lymph nodes as well? Lymph nodes are critical for the responses to vaccination.

We have now specified “lung and spleen” in the legend to Fig. 4 and also in the sub-heading in the results section. 

It would be very interesting to look at tissue distribution to lymph nodes, and this is a point both reviewers have brought up. Unfortunately, we do not have the samples available and no additional funding is available at the moment to repeat this study. Furthermore, we believe it would be difficult to justify the use of further animals solely for this purpose; however, we will investigate this in future experiments and thank you for the suggestion.

The variation in CFU/organ seem to be higher in animals that have high levels of protection in this study (e.g. A/J mice). Could the authors comment on this ? Is it possible to know the ML ratio in the 2 mice showing high CFU/lung from Fig 4A?

Since a different set of animals was used for flow cytometry and BCG distribution studies (outlined in the Methods section), it is not possible to directly compare the ML ratio and BCG CFU of individual animals.

With regards to variation in Fig. 4, we believe this is likely to be a technical effect – we detected very few bacilli in general, and especially in the lung, we did not detect any bacteria in some samples from lesser protected animals. This does not mean there are non – it just means that we were not able to detect them with the method used. As such, it is difficult to interpret the real variation within those mouse strains.

Not sure if the authors have the samples available, but it would have been interesting to see the profiles of Type I/II IFN in tissue at least by qPCR. This information could be of great interest for the tuberculosis field.

Unfortunately we do not have those samples available. It is an important point, and would indeed be very relevant. We plan to investigate this in future studies.

The authors should discuss the limitations of a relatively low number of biological experiments in this work.

We have added a comment to this effect at the end of the results section.

 10.5256/f1000research.15495.r32536 Reviewer response for version 1 Aguilo Nacho 1 Referee Department of Microbiology, Preventive Medicine, and Public Health, University of Zaragoza, Zaragoza, Spain

Competing interests: No competing interests were disclosed.

 9 4 2018 Copyright: © 2018 Aguilo N 2018 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve-with-reservations

The study by Zelmer et al. evaluates in different mouse genetic backgrounds whether the ratio of monocytes/macrophages vs T lymphocytes could affect BCG-induced protection. This study results are highly interesting as it evaluates in an animal model data derived from a previous vaccine efficacy trial. In addition, it highlights a critical question in the field of vaccine evaluation in mice, which is the genetic heterogeneity among inbred mouse strains. This issue is of particular importance when we are assessing vaccines which aim to be moved to humans. During the reading of the manuscript I have found some minor questions that I think should be amended to improve the quality of the study.

1. I think that the gating criteria used to enclose the monocyte/macrophage population (CD11b +CD11c int cells) is excluding myeloid populations that could be highly informative for the outcome of the study. This is particularly important in the case of the lungs, where naïve alveolar macrophages have a phenotype CD11b intCD11c +. Alveolar macrophages represent a first-line defense against respiratory pathogens, and therefore results could be very informative to know which is the status of these cells in the different strains used. More even if we consider that authors compare directly data from BCG-induced protection with the ML ratio in the lungs (Figure 2B).

2. I suggest to make the comparison between ML ratio and protection using a linear regression, so statistical significance of the shown data could be assessed.

3. I find highly interesting the data showing the differences in BCG persistence among the four strains, more even considering the influence of vaccine persistence in protective efficacy. This result evidences that live vaccine persistence does not only depends on the vaccine itself, but also on the host, and this is a question that should be considered in the future to study novel tuberculosis vaccines. My question is whether the authors could provide some data of BCG persistence in the local draining lymph nodes (in addition to lungs and spleen), as they are the primary sites of BCG biodistribution following subcutaneous/intradermal vaccination.

Is the work clearly and accurately presented and does it cite the current literature?

Yes

If applicable, is the statistical analysis and its interpretation appropriate?

Partly

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Yes

Are the conclusions drawn adequately supported by the results?

Partly

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

tuberculosis vaccines immunology

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.

Zelmer Andrea LSHTM

Competing interests: No competing interests were disclosed.

 19 6 2018

Dear Dr Aguilo

Thank you for reviewing our manuscript, and for providing timely and helpful comments and feedback. We discuss your individual points below.

The study by Zelmer  et al. evaluates in different mouse genetic backgrounds whether the ratio of monocytes/macrophages vs T lymphocytes could affect BCG-induced protection. This study results are highly interesting as it evaluates in an animal model data derived from a previous vaccine efficacy trial. In addition, it highlights a critical question in the field of vaccine evaluation in mice, which is the genetic heterogeneity among inbred mouse strains. This issue is of particular importance when we are assessing vaccines which aim to be moved to humans. During the reading of the manuscript I have found some minor questions that I think should be amended to improve the quality of the study.

1. I think that the gating criteria used to enclose the monocyte/macrophage population (CD11b +CD11c int cells) is excluding myeloid populations that could be highly informative for the outcome of the study. This is particularly important in the case of the lungs, where naïve alveolar macrophages have a phenotype CD11b intCD11c +. Alveolar macrophages represent a first-line defense against respiratory pathogens, and therefore results could be very informative to know which is the status of these cells in the different strains used. More even if we consider that authors compare directly data from BCG-induced protection with the ML ratio in the lungs (Figure 2B).

We have analysed the frequency of CD11b intCD11c +cells in the lungs of the four different mouse strains, and have added Fig. 1E, as well as described the data in the Results section and discussed it in the Discussion. In brief, the frequency of this population is significantly higher in C57Bl/6 mice, and seems to be associated with susceptibility of naïve mice, but not BCG-mediated protection.

2. I suggest to make the comparison between ML ratio and protection using a linear regression, so statistical significance of the shown data could be assessed.

This would indeed be an informative analysis; however, our experimental design does not allow it here. Since the data shown from flow cytometry and protection studies were acquired from two different sets of mice (see Methods section), the data are not paired, and it would not be appropriate to link them in a linear regression.

3. I find highly interesting the data showing the differences in BCG persistence among the four strains, more even considering the influence of vaccine persistence in protective efficacy. This result evidences that live vaccine persistence does not only depends on the vaccine itself, but also on the host, and this is a question that should be considered in the future to study novel tuberculosis vaccines. My question is whether the authors could provide some data of BCG persistence in the local draining lymph nodes (in addition to lungs and spleen), as they are the primary sites of BCG biodistribution following subcutaneous/intradermal vaccination.

This is an important point, and was also raised by Dr Gutierrez, the other reviewer. As you will see from our response to Dr Gutierrez’ comments, we unfortunately do not have the samples available and no additional funding is available at the moment to repeat this study. Furthermore, we believe it would be difficult to justify the use of further animals solely for this purpose; however, we are very interested to investigate the role of BCG distribution in vaccine efficacy in more detail in the future, and thank you for the suggestion.