Antibiotic Resistance Protocols, Fourth Edition Springer, 2024
Antibiotic Resistance Protocols, Fourth Edition Springer, 2024
Antibiotic Resistance Protocols, Fourth Edition Springer, 2024
Antibiotic
Resistance
Protocols
Fourth Edition
METHODS IN MOLECULAR BIOLOGY
Series Editor
John M. Walker
School of Life and Medical Sciences
University of Hertfordshire
Hatfield, Hertfordshire, UK
Edited by
Stephen H. Gillespie
Division of Infection and Global Health, School of Medicine, University of St Andrews, St Andrews, Fife, UK
Editor
Stephen H. Gillespie
Division of Infection and Global
Health, School of Medicine
University of St Andrews
St Andrews, Fife, UK
This Humana imprint is published by the registered company Springer Science+Business Media, LLC, part of Springer
Nature.
The registered company address is: 1 New York Plaza, New York, NY 10004, U.S.A.
Much has changed in the world of antibiotics resistance research since the first volume of this
protocol series was published in 2001. The most important factor is that the tide of
antibiotic resistance is rising steadily [1]. Fortunately, recognition of the scale of the
problem has also occurred with the publication of national and international reports that
set out the significant danger to human health posed by this problem [2]. International
organizations to fight antibiotic resistance, such as the Global Antibiotic Research and
Development Partnership (GARDP), have been established with the aim of developing
five new treatments by 2025 [3]. Similarly, academic groups and consortia have been
established to address antibiotic resistance from fundamental science through to implemen-
tation science.
When the first edition of this book was published, it took perhaps as much as a year to
sequence and assemble a genome. Now with the availability of high-throughput platforms,
sequencing of multiple strain has become commonplace, and it is now starting to take its
place in routine microbiology practice [4]. There remain many questions about how samples
are collected, strains isolated and sequenced, and the results integrated in the microbiolo-
gical workflow. Several of our chapters address these issues. Sequencing is also a powerful
tool to understand the development, emergence, and spread of resistance, as is discussed in
this volume. The challenge of integrating sequencing methodologies is not trivial, and this
subject is addressed in the current edition.
We now have molecular methods to detect and quantify the number of live bacteria that
will aid in the monitoring of treatment response in difficult infections caused by slow-
growing organisms such as Mycobacterium tuberculosis [5]. The sensitivity and specificity
of such a technique can be harnessed to make a diagnosis from specimens other than
sputum, such as stool [6], that may make it easier to diagnose tuberculosis in children.
Despite the increasing role of sequencing, growing microorganisms will probably
remain a major function of microbiology laboratories. The question is, however, can we
do it more rapidly and more cleverly? We present a number of techniques in this volume that
unequivocally answer that question in the affirmative with novel methods to test resistance
and interactions between antibiotics, physiological conditions or using innovative tools like
the hollow fiber, which can mimic human pharmacology well, or Raman spectroscopy that
enables changes in the bacteria to be measured in real time.
Sometimes the problem of antibiotic resistance can seem very large and only really be
encompassed numerically. In this volume, we present mathematical models that can describe
resistance within host [7].
v
vi Preface
The main purpose of this volume is to present a range of different techniques that the
reader can use to address important questions in antibiotic resistance. In the coming years, I
believe several of the techniques presented here will play an important role in answering
questions as to how to control antibiotic resistance, how to develop new agents, and how to
address the problems posed by microbes that have become resistant to our antibiotics.
References
1. Collaborators AR, Murray CJ, Ikuta KS, et al (2022) Global burden of bacterial antimicrobial
resistance in 2019: a systematic analysis. Lancet Lond Engl 399:629–655
2. O’Neill J (2016) Tackling drug resistant infections globally. Final report and recomendations, HM
Government, London
3. Global Antibiotic Research and Development Partnership. https://gardp.org/
4. Parcell BJ, Gillespie SH, Pettigrew KA, et al (2020) Clinical perspectives in integrating whole genome
sequencing into the investigation of healthcare and public health outbreaks – hype or help? J Hosp
Infect 109:1–9
5. Sabiiti W, Azam K, Farmer ECW, et al (2020) Tuberculosis bacillary load, an early marker of disease
severity: the utility of tuberculosis molecular bacterial load assay. Thorax 75:606–608
6. Musisi E, Sessolo A, Kaswabuli S, et al (2022) High mycobacterium tuberculosis bacillary loads
detected by tuberculosis molecular bacterial load assay in patient stool: a potential alternative for
nonsputum diagnosis and treatment response monitoring of tuberculosis. Microbiol Spectr 10:
e02100-21
7. Bowness R, Chaplain MAJ, Powathil GG, et al (2018) Modelling the effects of bacterial cell state and
spatial location on tuberculosis treatment: Insights from a hybrid multiscale cellular automaton model.
J Theor Biol 446:87 100
Contents
Preface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . v
Contributors. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ix
1 Using the Zebrafish Larval Model of Infection to Investigate Antibiotic
Efficacy and Combination Treatments Against Staphylococcus aureus. . . . . . . . . . . 1
Amy K. Tooke and Simon J. Foster
2 Three-Dimensional Rotary Culture to Model Mycobacterial
Biofilms in Low-Shear Detergent-Free Liquid Suspension . . . . . . . . . . . . . . . . . . . 11
Daire Cantillon and Simon J. Waddell
3 Rapid Gene Silencing Followed by Antimicrobial Susceptibility
Testing for Target Validation in Antibiotic Discovery . . . . . . . . . . . . . . . . . . . . . . . . 23
Chris Daniel, Sam Willcocks, and Sanjib Bhakta
4 Modified HT-SPOTi: An Antimicrobial Susceptibility Testing
to Evaluate Formulated Therapeutic Combinations
Against Bacterial Growth and Viability. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
Anushandan Navaratnarajah, Chris Daniel, and Sanjib Bhakta
5 Investigating Combination Therapy as a Means
to Enhance Activity and Repurpose Antimicrobials . . . . . . . . . . . . . . . . . . . . . . . . . . 43
Robert J. H. Hammond
6 Investigating Photoactive Antimicrobials as Alternatives (or Adjuncts)
to Traditional Therapy. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
Robert J. H. Hammond and Marianna Leite De Avellar
7 Using Hollow Fiber to Model Treatment
of Antimicrobial-Resistant Organisms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
Robert J. H. Hammond
8 A Microtiter Plate Assay at Acidic pH to Identify Potentiators
that Enhance Pyrazinamide Activity Against Mycobacterium tuberculosis. . . . . . . . 65
Christopher William Moon, Eleanor Porges,
Stephen Charles Taylor, and Joanna Bacon
9 Within-Host Mathematical Models of Antibiotic Resistance . . . . . . . . . . . . . . . . . . 79
Aminat Yetunde Saula, Gwenan Knight, and Ruth Bowness
10 Use of Individual-Based Mathematical Modelling
to Understand More About Antibiotic Resistance Within-Host. . . . . . . . . . . . . . . 93
Aminat Yetunde Saula, Christopher Rowlatt, and Ruth Bowness
11 Monitoring Live Mycobacteria in Real-Time Using
a Microfluidic Acoustic-Raman Platform . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
Mingzhou Chen, Vincent Baron, Björn Hammarström,
Robert J. H. Hammond, Peter Glynne-Jones,
Stephen H. Gillespie, and Kishan Dholakia
12 Phylogenetic Survival Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 121
Arturo Torres Ortiz and Louis Grandjean
vii
viii Contents
Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 227
Contributors
JOANNA BACON • Discovery Group, UK Health Security Agency, Porton Down, Salisbury, UK
IVAN BARILAR • Molecular and Experimental Mycobacteriology, Research Center Borstel,
Borstel, Germany
VINCENT BARON • Division of Infection and Global Health, School of Medicine, University of
St Andrews, St Andrews, UK
SANJIB BHAKTA • Mycobacteria Research Laboratory, Institute of Structural and Molecular
Biology, School of Natural Sciences, Birkbeck, University of London, London, UK
RUTH BOWNESS • Department of Mathematical Sciences, University of Bath, Bath, UK
DAIRE CANTILLON • Department of Global Health and Infection, Brighton and Sussex
Medical School, University of Sussex, Brighton, UK; Department of Tropical Biology,
LiverpoolLiverpool, UK
MARIE ANNE CHATTAWAY • Gastrointestinal Bacteria Unit, United Kingdom Health Security
Agency, London, UK
MINGZHOU CHEN • School of Physics, University of St Andrews, St Andrews, UK
CHRIS DANIEL • Mycobacteria Research Laboratory, Institute of Structural and Molecular
Biology, School of Natural Sciences, Birkbeck, University of London, London, UK
MARIANNA LEITE DE AVELLAR • School of Physics, University of St Andrews, St Andrews, UK
KISHAN DHOLAKIA • School of Physics, University of St Andrews, St Andrews, UK; Centre of
Light for Life, University of Adelaide, Adelaide, SA, Australia
MARGO DIRICKS • Molecular and Experimental Mycobacteriology, Research Center Borstel,
Borstel, Germany
VIOLA DREYER • Molecular and Experimental Mycobacteriology, Research Center Borstel,
Borstel, Germany
SIMON J. FOSTER • School of Biosciences, Krebs Institute, and Florey Institute, University of
Sheffield, Sheffield, UK; Newcastle University Biosciences Institute, Newcastle upon Tyne,
UK
STEPHEN H. GILLESPIE • Division of Infection and Global Health, School of Medicine,
University of St Andrews, St Andrews, UK
PETER GLYNNE-JONES • Faculty of Engineering and Physical Sciences, University of
Southampton, Southampton, UK
LOUIS GRANDJEAN • Department of Infection, Immunity and Inflammation, Institute of
Child Health, University College London, London, UK
BJÖRN HAMMARSTRÖM • Biophysics, KTH Royal Institute of Technology, Stockholm, Sweden
ROBERT J. H. HAMMOND • Division of Infection and Global Health, School of Medicine,
University of St Andrews, St Andrews, UK
GWENAN KNIGHT • Department of Infectious Disease Epidemiology, London School of Hygiene
and Tropical Medicine, London, UK
THOMAS A. KOHL • Molecular and Experimental Mycobacteriology, Research Center Borstel,
Borstel, Germany
MATTHIAS MERKER • Evolution of the Resistome, Research Center Borstel, Borstel, Germany
VANESSA MOHR • Molecular and Experimental Mycobacteriology, Research Center Borstel,
Borstel, Germany
ix
x Contributors
CHRISTOPHER WILLIAM MOON • Discovery Group, UK Health Security Agency, Porton Down,
Salisbury, UK
BARIKI MTAFYA • National Institute for Medical Research – Mbeya Medical Research Centre,
Mbeya, Tanzania
EMMANUEL MUSISI • Infectious Diseases Research Collaboration, Kampala, Uganda; Adroit
Biomedical & Bioentreprenuership Research Services, Kampala, Uganda; Division of
Infection and Global Health, School of Medicine, University of St Andrews, St Andrews, UK
ANUSHANDAN NAVARATNARAJAH • Mycobacteria Research Laboratory, Institute of Structural
and Molecular Biology, School of Natural Sciences, Birkbeck, University of London, London,
UK
STEFAN NIEMANN • Molecular and Experimental Mycobacteriology, Research Center Borstel,
Borstel, Germany
NYANDA ELIAS NTINGINYA • National Institute for Medical Research –, MbeyaMbeya,
Tanzania
KATARINA ORAVCOVA • School of Biodiversity, One Health and Veterinary Medicine,
University of Glasgow, Glasgow, UK
BENJAMIN J. PARCELL • Clinical Senior Lecturer and Honorary Consultant, Population
Health and Genomics, School of Medicine, University of Dundee, Dundee, UK
KERRY A. PETTIGREW • School of Biomedical Science, School of Health and Life Sciences,
Teesside University, Middlesbrough, UK
ELEANOR PORGES • Discovery Group, UK Health Security Agency, Porton Down, Salisbury,
UK
PASCAL QWARAY • National Institute for Medical Research –, MbeyaMbeya, Tanzania
CHRISTOPHER ROWLATT • Department of Mathematical Sciences, University of Bath, Bath,
UK
WILBER SABIITI • Division of Infection and Global Health, School of Medicine, University of St
Andrews, St Andrews, UK
AMINAT YETUNDE SAULA • Department of Mathematical Sciences, University of Bath, Bath,
UK
ABDULWAHAB SESSOLO • Infectious Diseases Research Collaboration, Kampala, Uganda
EMANUEL SICHONE • National Institute for Medical Research – Mbeya Medical Research
Centre, Mbeya, Tanzania
LINDSAY SONNENKALB • Molecular and Experimental Mycobacteriology, Research Center
Borstel, Borstel, Germany
WILLY SSENGOOBA • School of Biomedical Sciences, College of Health Sciences, Makerere
University, Kampala, Uganda
STEPHEN CHARLES TAYLOR • Pathogen Immunology Group, UK Health Security Agency,
Porton Down, Salisbury, UK
AMY K. TOOKE • School of Biosciences, University of Sheffield, Sheffield, UK; Newcastle
University Biosciences Institute, Newcastle upon Tyne, UK
ARTURO TORRES ORTIZ • Department of Infection, Immunity and Inflammation, Institute
of Child Health, University College London, London, UK
CHRISTIAN UTPATEL • Molecular and Experimental Mycobacteriology, Research Center
Borstel, Borstel, Germany
SIMON J. WADDELL • Department of Global Health and Infection, Brighton and Sussex
Medical School, University of Sussex, Brighton, UK
NATASHA WALBAUM • Division of Infection and Global Health, School of Medicine, University
of St Andrews, St Andrews, UK
Contributors xi
WILLIAM SAAVA WAMBI • School of Biomedical Sciences, College of Health Sciences, Makerere
University, Kampala, Uganda
SAM WILLCOCKS • Mycobacteria Research Laboratory, Institute of Structural and Molecular
Biology, School of Natural Sciences, Birkbeck, University of London, London, UK; Division
of Biosciences, College of Health, Medicine and Life Sciences, Brunel University London,
Uxbridge, UK
JOSEPHINE ZAWEDDE • Infectious Diseases Research Collaboration, Kampala, Uganda
Chapter 1
Abstract
There is an increasing need for new treatment regimens to combat antibiotic-resistant strains of bacteria.
Staphylococcus aureus is a clinically important, opportunist pathogen that has developed resistance to a range
of antibiotics. The zebrafish larval model of systemic disease has been increasingly utilized to elucidate
S. aureus virulence mechanisms and host-pathogen interactions. Here, we outline how this model can be
used to investigate the effects of different antibiotics alone and in combination against S. aureus.
Key words Zebrafish, Antibiotics, Antimicrobials, Drug screening, Staphylococcus aureus, In vivo
model, Infection, Combination treatment
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_1,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
1
2 Amy K. Tooke and Simon J. Foster
2 Materials
Fig. 1 A range of outcomes can occur after exposing zebrafish larvae to antibiotics. If the antibiotic is toxic
against uninfected larvae at a particular concentration it cannot be used in assays against infected larvae. If
the antibiotic is non-toxic, it can be either curative or non-curative. Non-curative doses do not prevent
bacterial growth and larval survival is not increased in the treated group. Curative doses can have a
bactericidal or bacteriostatic effect, which is established by determining the bacterial numbers in sampled
larvae before and after antibiotic treatment is given
3 Methods
3.1 Choice of Dug An in vitro minimum inhibitory concentration (MIC) assay carried
Concentrations Based out over 24 h of growth in a sterile 96 well plate is necessary to
on In Vitro Minimum determine the range of concentrations that should be tested against
Inhibitory infection in the zebrafish larvae. This method can give 3 biological
Concentration repeats for 2 different strains in one plate (see Note 4).
1. Grow an overnight culture of the strain(s) of S. aureus of
interest in 10 mL TSB in a 50 mL centrifuge tube at 37 °C,
shaking at 200 revolutions per minute (rpm) with antibiotics
for selection where appropriate. Grow one overnight culture
per biological repeat.
2. In a sterile 96-well plate (with lid) set up serial 2× dilutions of
your antibiotic(s) of interest in TSB, leaving columns 11 and
12 with TSB with no antibiotic, to use as controls (see Fig. 2).
3. Inoculate each well (columns 1–11) with 10 μL overnight
bacterial culture back diluted to OD600 0.01, leaving column
12 as a blank.
4. Seal the lid onto the plate with parafilm (see Note 5).
5. Incubate the plate overnight at 37 °C (see Note 6).
6. After 24 h of growth, read the OD600 of each well in an
automated plate reader, using column 12 as a blank. Column
11 should have a high bacterial density, which decreases
moving across the plate back to column 1. There should be
no growth in column 12; if there is, the media is contaminated.
7. The MIC can be interpreted as the lowest concentration of
antibiotic at which there is no bacterial growth.
Antibiotic Efficacy Testing in Zebrafish 5
Fig. 2 Layout of the MIC assay. Serial dilutions of antibiotics of interest are carried out across columns 1–10,
with column 11 a positive control with bacteria added with no antibiotic, and column 12 a blank with no
antibiotic or bacteria added (just TSB). Grey wells can be left empty. This method gives 3 biological repeats of
2 strains each per 96-well plate
3.2 Toxicity Assays Before an antibiotic is tested against infected larvae, it must be
and Choice of Solvent determined whether the drug itself is toxic to the zebrafish. This
is carried out by titrating a range of concentrations of the test
compound diluted in E3, which the larvae are then immersed
in. Survival is monitored over 92 h compared to an untreated
control. The highest concentration in the range tested where sur-
vival is 100% is the maximum of what can be used to treat infected
embryos, as toxic concentrations negate any potential antibiotic
curative effect.
1. Dechorionate larvae at approximately 24–28 hpf.
2. At 31 hpf, prepare a range of dilutions of the compounds of
interest in E3 buffer. Prepare them at 6× the final concentration
and add 50 μL antibiotic to each well containing 250 μL with a
larva in it (otherwise the concentration will be diluted when the
larva is transferred from its dish of E3 (see Note 7)).
3. Repeat step 2 with equivalent concentrations of solvent that
each antibiotic is dissolved in. For example, if adding 50 μL of
antibiotic in 50% v/v ethanol, add 50 μL of 50% v/v ethanol
minus the antibiotic (see Notes 8 and 9).
6 Amy K. Tooke and Simon J. Foster
3.3 Antibiotic To determine whether the antibiotic of interest degrades over the
Stability Assay time course of an experiment (e.g., up to 92 h) in E3, compounds
at the concentration used in in vivo experiments with the zebrafish
larvae are incubated in E3 at 28.5 °C for a range of timepoints
corresponding to the duration of the infection studies (e.g.,
0–92 hpi). Then bacteria are inoculated into the E3, incubated
overnight and live cells enumerated by plating. If cells grow this
means that the antibiotics have degraded throughout the time
course, whereas if growth is inhibited, this confirms the stability
of the compounds under the experimental conditions.
For a range of timepoints, incubate antibiotics in E3 for up to
92 h at the concentration to be tested with zebrafish larvae (see
Table 1). Additionally set up the same conditions with PBS as a
control to rule out any interactions between the E3 components
and the antibiotic.
1. Set up an overnight culture of your strain of interest in the
afternoon/evening of day 4 by inoculating a single colony into
10 mL of TSB (with antibiotics for selection where necessary)
and incubating overnight at 37 °C with shaking at 200 rpm.
2. On day 5, dilute the overnight culture to OD600 0.01 and add
15 μL to 135 μL of each antibiotic in E3 in a 96-well plate.
Additionally, dilute each antibiotic E3 mix 10 and 100 times
and inoculate these as well (see Fig. 3).
Table 1
Timepoints for setting up each condition in the antibiotic stability assay
Fig. 3 Layout of the 96-well plates and corresponding TSA plates for antibiotic stability assay. Bacteria are
incubated with antibiotics that have been incubated in PBS or E3 for predetermined amounts of time overnight
before plating out spots from each well on a TSA plate to establish whether the antibiotics have degraded or
are still effective against the bacteria
3.5 Determination of Larvae can be sampled from a population over the course of an
Live Bacterial experiment to determine the live bacterial numbers in the fish at
Numbers With/Without that given timepoint. This allows monitoring of whether antibiotic
Antibiotics treatments that appear to “cure” fish by increasing survival
8 Amy K. Tooke and Simon J. Foster
4 Notes
References
1. O’Neill J (2016) Tackling drug-resistant infec- 5. Torraca V, Mostowy S (2018) Zebrafish infec-
tions globally: final report and recommenda- tion: from pathogenesis to cell biology. Trends
tions. Review on Antimicrobial Resistance. Cell Biol 28(2):143–156. https://doi.org/10.
Wellcome Trust and HM Government. 1016/j.tcb.2017.10.002
https://wellcomecollection.org/works/ 6. Renshaw SA, Trede NS (2012) A model 450 mil-
thvwsuba lion years in the making: zebrafish and vertebrate
2. Murray CJ, Ikuta KS, Sharara F et al (2022) immunity. Dis Model Mech 5:38–47. https://
Global burden of bacterial antimicrobial resis- doi.org/10.1242/dmm.007138
tance in 2019: a systematic analysis. Lancet 7. Herbomel P, Thisse B, Thisse C (1999) Ontog-
399:629–655. https://doi.org/10.1016/ eny and behaviour of early macrophages in the
S0140-6736(21)02724-0 zebrafish embryo. Development 126(17):
3. Righi E, Scudeller L, Chiamenti M et al (2020) 3735–3745. https://doi.org/10.1242/dev.
In vivo studies on antibiotic combination for the 126.17.3735
treatment of carbapenem-resistant Gram-nega- 8. Prajsnar TK, McVicker G, Williams A et al
tive bacteria: a systematic review and meta- (2018) Use of larval zebrafish model to study
analysis protocol. BMJ Open Sci 4:e100055. within-host infection dynamics. Methods Mol
h t t p s : // d o i . o r g / 1 0 . 1 1 3 6 / b m j o s - Biol 1736:147–156. https://doi.org/10.
2019-100055 1007/978-1-4939-7638-6_14
4. Klug DM, Idiris FIM, Blaskovich MAT et al 9. McVicker G, Prajsnar TK, Williams A et al
(2021) There is no market for new antibiotics: (2014) Clonal expansion during Staphylococcus
this allows an open approach to research and aureus infection dynamics reveals the effect of
development. Wellcome Open Res 6:146. antibiotic intervention. PLoS Pathog 10:
https://doi.org/10.12688/wellcomeopenres. e1003959. https://doi.org/10.1371/journal.
16847.1 ppat.1003959
Chapter 2
Abstract
In vitro biofilm models have allowed researchers to investigate the role biofilms play in the pathogenesis,
virulence, and antimicrobial drug susceptibility of a wide range of bacterial pathogens. Rotary cell culture
systems create three-dimensional cellular structures, primarily applied to eukaryotic organoids, that better
capture characteristics of the cells in vivo. Here, we describe how to apply a low-shear, detergent-free rotary
cell culture system to generate biofilms of Mycobacterium bovis BCG. The three-dimensional biofilm model
forms mycobacterial cell aggregates in suspension as surface-detached biomass, without severe nutrient
starvation or environmental stress, that can be harvested for downstream experiments. Mycobacterium bovis
BCG derived from cell clusters display antimicrobial drug tolerance, presence of an extracellular matrix, and
evidence of cell wall remodeling, all features of biofilm-associated bacteria that may be relevant to the
treatment of tuberculosis.
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_2,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
11
12 Daire Cantillon and Simon J. Waddell
Fig. 1 Overview of the Rotary cell culture system. (a) RCCS 10 mL High Aspect Ratio Vessel (HARV). (i) Large
outlet port for harvesting biofilms; (ii) fill port butterfly valve; (iii) exit port butterfly valve; (iv) fill port screwcap
for addition of media; (v) exit port screwcap for removal of air bubbles; (vi) RCCS attachment point to the
vertical support. The semi-permeable membrane is visible on the back of the vessel. (b) RCCS-4DQ vertical
support and control panel, where the rotation of each vessel is independently controlled. RCCS vertical support
platform with 4 10 mL HARV mounted. (c) Centrifugal force (Fcf) is controlled through rotation of the vessels,
enabling biomass to accumulate in a low-shear liquid suspension. Gravity (g) collects biomass to the bottom of
non-rotating RCCS vessels
Generation of Mycobacterial Biofilms in Suspension 13
2 Materials
3 Methods
3.1 Preparation 1. For 1 L Middlebrook 7H9 ADC with 0.05% Tween 80 liquid
of Media culture media—Add 4.7 g Middlebrook 7H9 broth powder to
a 1 L Duran bottle with 900 mL ultrapure water and a mag-
netic stirring bar. Then add 500 μL Tween 80 and stir on a
magnetic plate until dissolved. Autoclave at 121 C for 15 min,
then cool to 52 C in a water bath. Add 100 mL Middlebrook
ADC growth supplement and mix well (see Note 1). Store at
4 C, warm to 37 C before use.
2. For 1 L Middlebrook 7H10 OADC with 0.5% glycerol agar
plates (approximately 25 plates)—Add 19.5 g Middlebrook
7H10 agar powder to a 1 L Duran bottle with 900 mL ultra-
pure water and a magnetic stirring bar. Add 5 mL glycerol, then
stir until dissolved. Autoclave at 121 C for 15 min, then cool
to 52 C in a water bath. Add 100 mL OADC and mix well (see
Note 2). Pour 30 mL per 90 mm Petri dish, allow to set, and
then store at 4 C. Warm to 37 C before use.
3. For 1 L Sauton media with 0.05% Tween 80—Add approxi-
mately 800 mL ultrapure water to a 1 L glass beaker with a
stirring bar. With the bar spinning, add 0.5 g (monobasic)
potassium phosphate, 0.5 g magnesium sulfate, 4.0 g L-aspar-
agine, 0.05 g ferric ammonium citrate, and 2.0 g citric acid.
Generation of Mycobacterial Biofilms in Suspension 15
3.3 Inoculation of the 1. RCCS vessels are sterile, single-use bioreactors that do not
RCCS Culture System need to be autoclaved (see Note 8). The day before inoculation
of the RCCS, pre-treat each of the RCCS vessels, adding 10 mL
sterile PBS through the fill port of the vessels. Attach to the
RCCS vertical support platform and place inside a humidified
tissue culture incubator at 37 C overnight (see Note 9).
2. The following day, remove the PBS from the vessels using the
exit port with a 5 mL syringe and discard the liquid (see
Note 10).
3. Using a 20 mL syringe, add 10 mL Sauton no Tween 80 myco-
bacterial inoculum through the fill port while simultaneously
removing air with a second 20 mL syringe from the exit port.
Ensure no air bubbles remain by drawing excess inoculum
through the fill and exit ports to remove air from the system.
Clean the necks of the fill and exit ports after closing with a
swab soaked in 70% ethanol to surface disinfect.
16 Daire Cantillon and Simon J. Waddell
4. Once all four vessels have been inoculated, attach to the RCCS
vertical support platform and place inside the humidified tissue
culture incubator at 37 C. Run the two-meter flat RCCS cable
from the incubator to the external RCCS interface controls (see
Note 11). The rotation of each vessel is controlled indepen-
dently; set two vessels to rotate at 15 rpm, the other two vessels
at 0 rpm as static controls (see Note 12).
3.4 Establishing and 1. Check the cultures daily for air bubbles in the vessels. Air
Maintaining Cell bubbles can disrupt the formation of biofilms and should be
Aggregates in the removed. To maintain the humidified environment of the incu-
RCCS System bator and reduce evaporation from the vessels, only open the
incubator once per day (see Note 13).
2. Remove air bubbles by adding Sauton no Tween 80 media
through the fill port with a 5 mL syringe, while drawing air
from the exit port with a second 5 mL syringe. Remove/add
the same culture volume to all vessels to maintain nutrient
balance across replicates or conditions (see Note 14).
3. Clean the necks of the fill and exit ports after closing with a
swab soaked in 70% ethanol (see Note 15) to maintain sterility.
4. Visible clusters of M. bovis BCG that are free-falling in suspen-
sion will be macroscopically visible by day 7. As the cellular
aggregates form and increase in size, decrease the speed of
rotation of the RCCS vessels to 10 rpm (see Note 16).
5. Media constituents, cells, macromolecules, chemicals, and anti-
microbial drugs may be added to the vessels through the fill
port to measure the impact on bacterial biofilm formation as
required, or aggregates may be allowed to form without inter-
vention for analysis in downstream applications.
3.5 Harvesting 1. Harvest mature M. bovis BCG biofilms after 21 days (see Fig. 2
Biofilms from the and Note 17). Carefully transfer the RCCS vessel to a class II
RCCS biosafety cabinet. Surface disinfect the external surfaces of the
vessel with a swab soaked in 70% ethanol to prevent contami-
nation being introduced when removing the mycobacterial
biomass.
2. Remove the cap from the large outlet-sample removal port and
loosen the fill port to prevent the vessel becoming air-locked.
3. Use a 25 mL serological pipette to gently remove the entire
vessel contents to sterile containers for downstream applica-
tions, such as drug sensitivity testing or biochemical analyses
(see Note 18).
Generation of Mycobacterial Biofilms in Suspension 17
Fig. 2 Representative images of M. bovis BCG cell aggregates cultured in the RCCS system. Mycobacterial
biofilm formation in the RCCS system with rotation at day 7 (a), day 14 (b), and day 21 (c). Replicate vessels
without rotation as control cultures at day 7 (d), day 14 (e), and day 21 (f). Control vessels (d, e, f) were
agitated briefly for the purposes of imaging the biomass settled at the bottom of the vessels. Black arrows
mark mycobacterial biomass grown in these 10 mL HARV RCCS vessels that are 6.5 cm in diameter
18 Daire Cantillon and Simon J. Waddell
4 Notes
Acknowledgments
References
1. Gonzalez JF, Hahn MM, Gunn JS (2018) 8. Cantillon D, Wroblewska J, Cooper I, New-
Chronic biofilm-based infections: skewing of port MJ, Waddell SJ (2021)
the immune response. Pathog Dis 76(3): Three-dimensional low shear culture of Myco-
fty023 bacterium bovis BCG induces biofilm forma-
2. Sharma D, Misba L, Khan AU (2019) Antibio- tion and antimicrobial drug tolerance. NPJ
tics versus biofilm: an emerging battleground Biofilms Microbiomes 7(1):12
in microbial communities. Antimicrob Resist 9. Ojha AK, Baughn AD, Sambandan D, Hsu T,
Infect Control 8:76 Trivelli X, Guerardel Y et al (2008) Growth of
3. Ciofu O, Moser C, Jensen PO, Hoiby N Mycobacterium tuberculosis biofilms containing
(2022) Tolerance and resistance of microbial free mycolic acids and harbouring drug-
biofilms. Nat Rev Microbiol 20(10):621–635 tolerant bacteria. Mol Microbiol 69(1):
4. Canetti G (1955) Tubercle bacillus in the pul- 164–174
monary lesion of man: histobacteriology and 10. Flores-Valdez MA, Aceves-Sanchez MJ, Peter-
its bearing on the therapy of pulmonary tuber- son EJR, Baliga N, Bravo-Madrigal J, De la
culosis, vol 231. Springer, New York, p 480 Cruz-Villegas MA et al (2020) Transcriptional
5. Chakraborty P, Bajeli S, Kaushal D, Radotra portrait of M. bovis BCG during biofilm pro-
BD, Kumar A (2021) Biofilm formation in duction shows genes differentially expressed
the lung contributes to virulence and drug tol- during intercellular aggregation and substrate
erance of Mycobacterium tuberculosis. Nat attachment. Sci Rep 10(1):12578
Commun 12(1):1606 11. Trivedi A, Mavi PS, Bhatt D, Kumar A (2016)
6. Basaraba RJ, Ojha AK (2017) Mycobacterial Thiol reductive stress induces cellulose-
biofilms: revisiting tuberculosis bacilli in extra- anchored biofilm formation in Mycobacterium
cellular necrotizing lesions. Microbiol Spectr tuberculosis. Nat Commun 7:11392
5(3):10.1128/microbiolspec.TBTB2- 12. Ackart DF, Hascall-Dove L, Caceres SM, Kirk
0024-2016. NM, Podell BK, Melander C et al (2014)
7. Lenaerts AJ, Hoff D, Aly S, Ehlers S, Expression of antimicrobial drug tolerance by
Andries K, Cantarero L et al (2007) Location attached communities of Mycobacterium tuber-
of persisting mycobacteria in a Guinea pig culosis. Pathog Dis 70(3):359–369
model of tuberculosis revealed by r207910. 13. Keating T, Lethbridge S, Allnutt JC, Hendon-
Antimicrob Agents Chemother 51(9): Dunn CL, Thomas SR, Alderwick LJ et al
3338–3345 (2021) Mycobacterium tuberculosis modifies
Generation of Mycobacterial Biofilms in Suspension 21
Abstract
Mycobacterium tuberculosis is the main causative agent of tuberculosis (TB)—an ancient yet widespread
global infectious disease to which 1.6 million people lost their lives in 2021. Antimicrobial resistance
(AMR) has been an ongoing crisis for decades; 4.95 million deaths were associated with antibiotic resistance
in 2019. While AMR is a multi-faceted problem, drug discovery is an urgent part of the solution and is at
the forefront of modern research.
The landscape of drug discovery for TB has undoubtedly been transformed by the development of high-
throughput gene-silencing techniques that enable interrogation of every gene in the genome, and their
relative contribution to fitness, virulence, and AMR. A recent advance in this area is CRISPR interference
(CRISPRi). The application of this technique to antimicrobial susceptibility testing (AST) is the subject of
ongoing research in basic science.
CRISPRi technology can be used in conjunction with the high-throughput SPOT-culture growth
inhibition assay (HT-SPOTi) to rapidly evaluate and assess gene essentiality including non-essential,
conditionally essential (by using appropriate culture conditions), and essential genes. In addition, the
HT-SPOTi method can develop drug susceptibility and drug resistance profiles.
This technology is further useful for drug discovery groups who have designed target-based inhibitors
rationally and wish to validate the primary mechanisms of their novel compounds’ antibiotic action against
the proposed target.
Key words Mycobacteria, CRISPRi, Gene essentiality, Antimicrobial susceptibility testing, Target
validation
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_3,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
23
24 Chris Daniel et al.
2 Materials
2.1 Golden Gate 1. Golden Gate Cloning (NEBridge® Golden Gate Assembly Kit,
Assembly [16] BsmBI-v2, NEB).
2. pIRL2 plasmid (Addgene #163631).
3. Oligonucleotides.
4. PCR buffer (containing MgCl2).
5. Molecular biology grade water, nuclease free (e.g.,
HyClone™).
6. Competent E. coli cells (e.g., E. coli DH5α).
Rapid Gene Silencing Followed by Antimicrobial Susceptibility Testing. . . 25
2.4 Plasmid DNA 1. Plasmid DNA extraction kit (e.g., Qiagen QIAprep Spin
Extraction and Miniprep Kit).
Sequencing 2. 1.5 mL Microcentrifuge tubes.
3 Methods
3.1 CRISPRi Mutant 1. Select the gene of interest for silencing. Identify a PAM
Generation Using sequence on the non-coding strand using sequences recog-
Golden Gate Assembly nized and identified by the Sth1dCas9 [4, 6]. There are differ-
ent strategies for selecting and designing sgRNA using either
the coding or the non-coding strand [4], the method we
describe has been optimized for this specific dCas9 system (see
Note 1).
2. Immediately upstream of the PAM, select 20–25 bp to be the
sgRNA. inhA from M. tuberculosis H37Rv is used as an exam-
ple below [8].
inhA (Rv1484) PAM 5′ CCAGAAT 3′
sgRNA 5′ GTCGGTGATGATTCCGCTAA 3′
5′ AACGGATTCTGGTTAGCGGAATCATCACCGACTC 3′
3′ TTGCCTAAGACCAATCGCCTTAGTAGTGGCTGAG 5′
Component Volume
HotStar Taq master mix 2x 25 μL
Molecular-grade water 15 μL
2 μM of sgRNA and pIRL2 primer 5 μL
Bacterial suspension 5 μL
Rapid Gene Silencing Followed by Antimicrobial Susceptibility Testing. . . 29
3.2 Agarose Gel 1. Prepare a 1% mix of agarose and 1 × TAE buffer for the gel.
Electrophoresis Dissolve 1 g of agarose in 100 mL of 1 × TAE in a microwave
for 2 min, gently swirling every 30 s. Let the solution cool, add
1 μL of GelRed® (10,000 × concentrated) and pour into a gel
tray. Add an appropriate comb ensuring no bubbles are in the
gel, leave to cool and completely set.
2. Once cooled, carefully remove the comb and the gel from the
casting tray. Place the gel tray into electrophoresis unit and
cover with 1 × TAE.
3. Add the appropriate volume (see recommendation by manu-
facturer) of ladder to the first well, mix 2 μL of loading dye to
the PCR product and add 10 μL of the final solution to
each well.
4. Connect the power supply in the correct orientation (negative
to positive), and run the gel for approximately 40 min at 100 V.
5. Once run, disconnect the power supply and visualize the gel
using a gel imaging system. If the cloning is successful, based
on the reverse pIRL2 primer provided a band should be visible
at approximately 160 bp.
3.3 Plasmid DNA 1. When the cultures have grown, make glycerol stocks using
Extraction and 500 μL of sterile 50% glycerol and 500 μL of culture.
Sequencing 2. Using the remaining culture, extract the plasmid DNA accord-
ing to the manufacturer’s instructions using the Qiagen QIA-
prep Spin Miniprep Kit.
3. Once extracted, measure the concentration, A260/280,
A260/230 using a nanodrop. Sequence the purified plasmid
to verify correct cloning of the sgRNA.
4 Notes
Rock et al
(PMID:
Feature 28165460) Latest scores
type Feature value Coefficient Standard Error scores (v4)
References
1. WHO (2022) Global Tuberculosis 11. Faulkner V, Cox AA, Goh S et al (2021) Re-
Report:2022 sensitization of Mycobacterium smegmatis to
2. Murray CJ, Ikuta KS, Sharara F et al (2022) Rifampicin Using CRISPR Interference
Global burden of bacterial antimicrobial resis- Demonstrates Its Utility for the Study of
tance in 2019: a systematic analysis. Lancet Non-essential Drug Resistance Traits. Front
6736:02724-0 Microbiol 11:619427
3. Andries K, Villellas C, Coeck N et al (2014) 12. Evangelopoulos D, Bhakta S (2010) Rapid
Acquired resistance of mycobacterium tuber- methods for testing inhibitors of mycobacterial
culosis to bedaquiline. PLoS One 9:e102135 growth. In: Gillespie S, McHugh T (eds) Anti-
4. Rock JM, Hopkins FF, Chavez A et al (2017) biotic resistance protocols. Humana Press, pp
Programmable transcriptional repression in 193–201
mycobacteria using an orthogonal CRISPR 13. Danquah CA, Maitra A, Gibbons S et al (2016)
interference platform. Nat Microbiol 2:16274 HT-SPOTi: a rapid drug susceptibility test
5. Rock J (2019) Tuberculosis drug discovery in (DST) to evaluate antibiotic resistance profiles
the CRISPR era. PLoS Pathog 15:1–10 and novel chemicals for anti-infective drug dis-
covery. Curr Protoc Microbiol 40:17.8.1–
6. Bosch B, DeJesus MA, Poulton NC et al 17.8.12
(2021) Genome-wide gene expression tuning
reveals diverse vulnerabilities of M. tuberculo- 14. Abrahams GL, Kumar A, Savvi S et al (2012)
sis. Cell 184:4579–4592.e24 Pathway-selective sensitization of mycobacte-
rium tuberculosis for target-based whole-cell
7. Li S, Poulton NC, Chang JS et al (2022) CRIS- screening. Chem Biol 19:844–854
PRi chemical genetics and comparative geno-
mics identify genes mediating drug potency in 15. Kurepina N, Chen L, Composto K et al (2022)
Mycobacterium tuberculosis. Nat Microbiol CRISPR inhibition of essential peptidoglycan
7:766–779 biosynthesis genes in mycobacterium abscessus
and its impact on b-lactam susceptibility. Anti-
8. McNeil MB, Keighley LM, Cook JR et al microb Agents Chemother 66:e0009322
(2021) CRISPR interference identifies vulner-
able cellular pathways with bactericidal pheno- 16. Bird JE, Marles-Wright J, Giachino A (2022) A
types in Mycobacterium tuberculosis. Mol user’s guide to golden gate cloning methods
Microbiol 116:1033–1043 and standards. ACS Synth Biol 11:3551–3563
9. Rock JM, Hopkins F, Chavez A et al (2016) 17. Parish T, Brown AC (2009) Mycobacteria Pro-
Programmable transcriptional repression in tocols Second Edition Humana Press Totowa
mycobacteria using an orthogonal CRISPR NJ
interference platform. Nat Microbiol 2:16274 18. Gupta R, Rohde KH (2023) Implementation
10. El Bakali J, Blaszczyk M, Evans JC et al (2023) of a mycobacterial CRISPRi platform in Myco-
Chemical validation of mycobacterium tuber- bacterium abscessus and demonstration of the
culosis phosphopantetheine adenylyltransferase essentiality of ftsZ. Tuberculosis 138:102292.
using fragment linking and CRISPR interfer- https://doi.org/10.1016/j.tube.2022.
ence. Angew Chem Int Ed Engl 62: 102292
e202300221
Chapter 4
Abstract
Antimicrobial resistance (AMR) poses a serious threat to global health, potentially causing 10 million
deaths per year globally by 2050. To tackle AMR, researchers from all around the world have generated a
selection of various formulated (viz. nanoparticulate, liposomal) therapeutic combinations to be evaluated
for new antimicrobial drug discovery. To meet the urgent need for accelerating new antibacterial drug
development, we need rapid but reliable whole-cell assay methods and models to test formulated therapeu-
tic combinations against several pathogens in different in vitro conditions as models of actual infections.
Over the past two decades, high-throughput spot-culture growth inhibition assay (HT-SPOTi) has been
demonstrated to be a gold-standard drug susceptibility method for evaluating novel chemotherapeutic
entities and existing drugs against various microbes of global concern. Our modified HT-SPOTi method
serves the purpose of evaluating drug combinations against Gram-positive/negative microorganisms as well
as acid-fast bacilli. The newly developed and modified HT-SPOTi assay builds upon the limitations of our
previously published method to incorporate antimicrobial susceptibility testing with formulated therapeutic
combinations. The modified HT-SPOTi is compared with a range of other antimicrobial susceptibility
testing methods and validated using a library of existing antibiotics as well as formulated therapeutic
combinations. The modified HT-SPOTi assay can serve as an efficient and reliable high-throughput drug
screening platform to discover new potential antimicrobial molecules, including as part of therapeutic
formulations.
This chapter describes the generation of drug susceptibility profile for formulated therapeutic combina-
tions using modified HT-SPOTi in a semi-automated system.
Key words Antimicrobial resistance (AMR), ESKAPE pathogens, Formulated therapeutic combina-
tions, Antimicrobial susceptibility testing (AST), Combination therapy
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_4,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
35
36 Anushandan Navaratnarajah et al.
Fig. 1 Workflow of modified HT-SPOTi. (Created with BioRender.com): On the 96-well plate, different test
compounds or formulations were added into row A to G. From column 1 to column 11 show two-fold
decreasing concentrations of the respective compound/formulation. A dark circle around the well is to indicate
the minimum inhibitory concentration of the compound/formulation tested. Abbreviations: ZN Ziehl–Neelsen,
DMSO Dimethyl sulfoxide, MIC Minimum inhibitory concentration)
38 Anushandan Navaratnarajah et al.
2 Materials
3 Methods
3.3 Agar Plate 1. After autoclaving the Middlebrook 7H10 (M7H10) agar
Preparation medium supplemented with 0.5% glycerol, keep the agar hot
(~50 °C) in a steam bath and add 10% (v/v) OADC to it. Place
it in a water bath kept at a temperature of 55° to 60 °C to
prevent it from solidifying.
2. Sterilize the dispensing cassette and tubing by running warm
70% ethanol through followed by priming it with hot (~65 °C)
sterile water (see Note 5).
3. Program the dispenser to suit the plate type (96-well, flat
bottom) and the volume of agar to be dispensed (e.g.,
200 μL) (see Note 6).
40 Anushandan Navaratnarajah et al.
4. Dispense the agar medium and leave the plates in the safety
cabinet undisturbed to solidify (~5 min). Once set, cover the
plates with lid (see Note 7).
3.5 Spotting 1. Sterilize the dispensing cassette and tubing by running warm
Bacterial Cultures 70% ethanol through followed by priming it with hot (~65 °C)
sterile water.
2. Program the dispenser to suit the plate type (96-well, F bot-
tom), volume of bacteria to be dispensed (2 μL), and the wells
in column 1–11 (see Note 10).
3. Dispense the bacteria on growth medium. Leave the plates in
the safety cabinet undisturbed for 5 min.
4. Seal the plates with parafilm, and leave in a static incubator at
37 °C and record spots observed daily.
4 Notes
Acknowledgments
References
1. Murray CJ, Ikuta KS, Sharara F et al (2022) 4. Evangelopoulos D, Bhakta S (2010) Rapid
Global burden of bacterial antimicrobial resis- methods for testing inhibitors of mycobacterial
tance in 2019: a systematic analysis. Lancet growth. Methods Mol Biol 642:193–201
399:629–655 5. Daraee H, Etemadi A, Kouhi M et al (2016)
2. Saifullah B, Chrzastek A, Maitra A et al (2017) Application of liposomes in medicine and drug
Novel anti-tuberculosis nano delivery formula- delivery. Artif Cells Nanomedicine Biotechnol
tion of ethambutol with graphene oxide. Mole- 44:381–391
cules 22(10):1560 6. Danquah CA, Maitra A, Gibbons S, et al
3. Wang Y, Grainger DW (2022) Regulatory con- (2016) HT-SPOTi: a rapid drug susceptibility
siderations specific to liposome drug develop- test (DST) to evaluate antibiotic resistance pro-
ment as complex drug products. Front Drug files and novel chemicals for anti-infective drug
Deliv 2:901281 discovery. 17.8.1–17.8.12
42 Anushandan Navaratnarajah et al.
Abstract
Current clinical practice assumes that a single antibiotic given as a bolus or as a course will successfully treat
most infections. In modern medicine, this is becoming less and less true with drug-resistant, multi-drug-
resistant, extensively drug-resistant, and untreatable infections becoming more common. Where single-
drug therapy (monotherapy) fails, we will turn to multi-drug therapy. Alternatively, combination therapy
could be useful to prevent the emergence of resistance. Multi-drug therapy is already standard for some
multi-drug resistant infections and is the standard for the treatment of some pathogens such as Mycobacte-
rium tuberculosis.
The use of combination therapy for everyday infections could be a clear course out of the current AMR
crisis we are facing. With every additional drug added to a combination (n + 1) the likelihood of the
pathogen evolving resistance drops exponentially.
Many generic antibiotics are cheap to manufacture as they have fallen out of patent protection but are less
effective at pharmacologically effective doses due to overuse in the past. Combination therapy can combine
these generic compounds into cocktails that can not only treat susceptible and resistant infections but can
also reduce the risk of new resistances arising and can resuscitate the use of antimicrobials once thought
defunct.
In this chapter, we will summarize theory behind combination therapy and standard in vitro
methodologies used.
Key words Antibiotic resistance, Combination therapy, Therapeutics, Chequerboard analysis, Syn-
ergy, Fractional inhibition concentration
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_5,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
43
44 Robert J. H. Hammond
1.1 Definition and Combination therapy refers to the use of two or more antimicrobial
Types of Combination agents in a single treatment regimen. These agents can be from
Therapy similar or different classes with similar or highly dissimilar mechan-
isms of action. Indeed, as we will see, where combinations include
agents with dissimilar but complementary mechanisms of action
the overall effect can be enhanced [12]. Synergistic, additive, indif-
ferent, and antagonistic interactions can occur between agents in a
combination therapy regimen as set out below [11, 13–15].
the FICI falls above the crucial value of 0.5 other language is
necessary to define the reaction. An FICI of 0.5–2 will indicate an
additive reaction, 2–4 is defined as an indifferent reaction and an
FICI of >4 is an antagonistic reaction [21, 22] where the antibiotics
are actually inhibiting each other’s activity and causing a poorer
outcome than monotherapy.
In an in vitro case study presented below a combination of only
2 agents is used. This is the standard use-case of synergy in medi-
cine, however there are exceptions. The standard intensive phase
therapy for tuberculosis treatment includes 4 antimicrobials used in
tandem and cases of 3 or more agents used in MDR infections are
not uncommon [23]. Where this occurs, the mathematics for
defining synergy no longer work and will need to be redefined in
the future.
2 Materials
2.1 Bacterial Growth 1. Mueller Hinton cation adjusted broth (Sigma-Aldrich, UK)
(see Note 1).
2. Resistant organism, in this case a KPC Klebsiella pneumonia
(ATCC BAA 1705) (see Note 2).
3. Antimicrobials of choice, in this instance colistin and Zidovu-
dine (Sigma-Aldrich, UK) (see Notes 3 and 4).
4. Phosphate Buffered Saline (PBS) (Sigma-Aldrich, UK).
5. Sterile water.
6. 96- or 384-well plates (Nunc, UK) (see Note 5).
7. Appropriate multichannel pipettor (Gilson, or similar) (see
Notes 6 and 7).
8. Incubator (Genlab (Fisher Scientific)).
9. Spectrophotometer (BioChrom WPA Co8000, UK).
10. 96- or 384-well plate reader with the ability to detect at
600 nm absorbance (Clariostar Pro).
11. Virkon solution (2%) for decontamination (Sigma-
Aldrich, UK).
12. Microcentrifuge tubes (Eppendorf tubes®, (Sigma-
Aldrich, UK)).
3 Methods
3.2 Checkerboard A similar assay was conducted to establish if any synergistic effect
existed when CSS and AZT were combined. This checkerboard
assay can be thought of as a two-dimensional BMD.
1. Dilute the experiment drugs in the same way as noted above in
Subheading 3.1, step 2.
2. Add the drug solution to the 96-well plate in the manner
displayed in Fig. 1. Add sterile media and bacteria as noted
above to make a final volume of 200μL and bacterial concen-
tration of ~5 105 cfu/mL
3. Measure the OD for all experimental case results using the
appropriate plate reader at 600 nm.
CSS mg/L
4 2 1 0.5 0.25 0.125 0.063 0.031 0.016 0.008 0
64 324 194 129 97 80 72 68 66 65 65 64
32 164 98 65 49 40 36 34 33 33 32 32
16 84 50 33 25 20 18 17 17 16 16 16
AZT 8 44 26 17 13 10 9 9 8 8 8 8
mg/L 4 24 14 9 7 5 5 4 4 4 4 4
2 14 8 5 4 3 2 2 2 2 2 2
1 9 5 3 2 2 1 1 1 1 1 1
0 4 2 1 1 0 0 0 0 0 0 0
Fig. 1 Schematic representation of the layout of checkerboard assay for AZT + CSS
Evaluating Antibiotic Combinations 47
4 Notes
References
1. Gajdacs M, Abrok M, Lazar A, Burian K 9. Soothill G, Hu Y, Coates A (2013) Can
(2020) Increasing relevance of gram-positive we prevent antimicrobial resistance by using
cocci in urinary tract infections: a 10-year anal- antimicrobials better? Pathogens 2(2):
ysis of their prevalence and resistance trends. 4 2 2 – 4 3 5 . h t t p s : // d o i . o r g / 1 0 . 3 3 9 0 /
Sci Rep 10(1):17658. https://doi.org/10. pathogens2020422
1038/s41598-020-74834-y 10. Siriyong T, Murray RM, Bidgood LE, Young
2. Vasala A, Hytonen VP, Laitinen OH (2020) SA, Wright F, Parcell BJ et al (2019) Dual
Modern tools for rapid diagnostics of antimi- β-lactam combination therapy for multi-drug
crobial resistance. Front Cell Infect Microbiol resistant Pseudomonas aeruginosa infection:
10:308. https://doi.org/10.3389/fcimb. enhanced efficacy in vivo and comparison with
2020.00308 monotherapies of penicillin-binding protein
3. Galata V, Laczny CC, Backes C, Hemmrich- inhibition. Sci Rep 9(1):9098. https://doi.
Stanisak G, Schmolke S, Franke A, Meese E, org/10.1038/s41598-019-45550-z
Herrmann M, von Müller L, Plum A, Müller R, 11. Coates ARM, Hu Y, Holt J, Yeh P (2020)
Stähler C, Posch AE, Keller A et al (2019) Antibiotic combination therapy against resis-
Integrating culture-based antibiotic resistance tant bacterial infections: synergy, rejuvenation
profiles with whole-genome sequencing data and resistance reduction. Expert Rev Anti-
for 11,087 clinical isolates. Genomics Proteo- Infect Ther 18(1):5–15. https://doi.org/10.
mics Bioinformatics 17(2):169. https://doi. 1080/14787210.2020.1705155
org/10.1016/j.gpb.2018.11.002 12. Hu Y, Liu Y, Coates A (2019) Azidothymidine
4. O’Neill J (2016) Tackling drug-resistant infec- produces synergistic activity in combination
tions globally: final report and recommenda- with Colistin against antibiotic-resistant Enter-
tions. Government of the United Kingdom obacteriaceae. Antimicrob Agents Chemother
5. Raymond B (2019) Five rules for resistance 63(1). https://d oi.o rg/10.1 128 /aac.
management in the antibiotic apocalypse, a 01630-18
road map for integrated microbial manage- 13. Deshpande D, Srivastava S, Nuermberger E,
ment. Evol Appl 12(6):1079–1091. https:// Pasipanodya JG, Swaminathan S, Gumbo T
doi.org/10.1111/eva.12808 (2016) Concentration-dependent synergy and
6. Wilkinson A, Ebata A, MacGregor H (2018) antagonism of linezolid and Moxifloxacin in
Interventions to reduce antibiotic prescribing the treatment of childhood tuberculosis: the
in LMICs: a scoping review of evidence from dynamic duo. Clin Infect Dis 32016:S88–S94
human and animal health systems. Antibiotics 14. Dutta NK, Pinn ML, Karakousis PC (2014)
(Basel) 8(1). https://doi.org/10.3390/ Reduced emergence of isoniazid resistance
antibiotics8010002 with concurrent use of thioridazine against
7. Manyi-Loh C, Mamphweli S, Meyer E, Okoh acute murine tuberculosis. Antimicrob Agents
A (2018) Antibiotic use in agriculture and its Chemother 58(7):4048–4053. https://doi.
consequential resistance in environmental org/10.1128/AAC.02981-14
sources: potential public health implications. 15. Tripathi RP, Tewari N, Dwivedi N, Tiwari VK
Molecules 23(4). https://doi.org/10.3390/ (2005) Fighting tuberculosis: an old disease
molecules23040795 with new challenges. Med Res Rev 25(1):
8. Bebell LM, Muiru AN (2014) Antibiotic use 93–131. https://doi.org/10.1002/med.
and emerging resistance: how can resource- 20017
limited countries turn the tide? Glob Heart 16. Bonapace CR, White RL, Friedrich LV, Bosso
9(3):347–358. https://doi.org/10.1016/j. JA (2000) Evaluation of antibiotic synergy
gheart.2014.08.009 against Acinetobacter baumannii: a comparison
Evaluating Antibiotic Combinations 49
with Etest, time-kill, and checkerboard meth- and dominance. Dominance and antibiotic
ods. Diagn Microbiol Infect Dis 38(1):43–50. combination]. Pathol Biol (Paris) 35(5):
https://doi.org/10.1016/s0732-8893(00) 495–500
00163-2 22. Rakholiya KD, Kaneria MJ, Chanda SV (2013)
17. Long H, Miller SF, Strauss C, Zhao C, Chapter 11 – medicinal plants as alternative
Cheng L, Ye Z et al (2016) Antibiotic treat- sources of therapeutics against multidrug-
ment enhances the genome-wide mutation rate resistant pathogenic microorganisms based on
of target cells. Proc Natl Acad Sci U S A 113: their antimicrobial potential and synergistic
E2498. https://doi.org/10.1073/pnas. properties. In: Rai MK, Kon KV (eds) Fighting
1601208113 multidrug resistance with herbal extracts,
18. Kloprogge F, Hammond R, Copas A, Gillespie essential oils and their components. Academic,
SH, Della PO (2019) Can phenotypic data San Diego, pp 165–179
complement our understanding of antimyco- 23. Umemura T, Kato H, Hagihara M, Hirai J,
bacterial effects for drug combinations? J Anti- Yamagishi Y, Mikamo H (2022) Efficacy of
microb Chemother 74:3530. https://doi.org/ combination therapies for the treatment of
10.1093/jac/dkz369 multi-drug resistant gram-negative bacterial
19. Kloprogge F, Hammond R, Copas A, Gillespie infections based on meta-analyses. Antibiotics
SH, Della PO (2019) Can phenotypic data (Basel) 11(4). https://doi.org/10.3390/
complement our understanding of antimyco- antibiotics11040524
bacterial effects for drug combinations? J Anti- 24. Koeth LM, King A, Knight H, May J, Miller
microb Chemother 74(12):3530–3536. LA, Phillips I et al (2000) Comparison of
https://doi.org/10.1093/jac/dkz369 cation-adjusted Mueller–Hinton broth with
20. Lewis K (2008) Multidrug tolerance of bio- Iso-Sensitest broth for the NCCLS broth
films and persister cells. Curr Top Microbiol microdilution method. J Antimicrob Che-
Immunol 322:107–131 mother 46(3):369–376. https://doi.org/10.
21. Drugeon HB, Dellamonica P, Caillon J (1987) 1093/jac/46.3.369
[Synergy, addition, indifference, antagonism
Chapter 6
Abstract
Photodynamic therapy (PDT) is an established therapy used for the treatment of cutaneous skin cancers and
other non-infective ailments. There has been recent interest in the opportunity to use aPDT (antimicrobial
PDT) to treat skin and soft tissue infections. PDT utilizes photosensitizers that infiltrate all cells and
“sensitize” them to a given wavelength of light. The photosensitizer is simply highly absorbent to a given
wavelength of light and when excited will produce, in the presence of oxygen, damaging oxygen radicals
and singlet oxygen. Bacterial cells are comparatively poor at combatting oxidative stress when compared
with human cells therefore a degree of selective toxicity can be achieved with aPDT.
In this chapter, we outline methodologies for testing aPDT in vitro using standard lab equipment.
Key words Antibiotic combination, Photodynamic therapy, Drug development, MDR therapy,
Photosensitizers, Susceptibility
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_6,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
51
52 Robert J. H. Hammond and Marianna Leite De Avellar
2 Materials
3 Methods
3.1 Bacterial Strains 1. Bacterial culture S. aureus (ATCC 25923), S. epidermidis, and
and Growth Conditions (MRSA/S. pyogenes) from glycerol stock (50%) grown in
Mueller-Hinton cation adjusted broth (MHB2) (Sigma-
Aldrich) at 37 C for 16 h.
2. The optical density of the culture should be measured on a Cell
density meter (Biochrome WPA CO8000, UK) at wavelength
600 nm (OD600) and diluted to be used as final OD600 equal
to 0.001 (approximately 105 cfu/mL) or 0.1 (approximately
107 cfu/mL).
4 Notes
References
1. Alamon-Reig F, Marti-Marti I, Loughlin CR, 3. Boltes Cecatto R, Siqueira de Magalhães L,
Garcia A, Carrera C, Aguilera-Peiro P (2022) Fernanda Setúbal Destro Rodrigues M,
Successful treatment of facial cutaneous leish- Pavani C, Lino-dos-Santos-Franco A, Teixeira
maniasis with photodynamic therapy. Indian J Gomes M et al (2020) Methylene blue
Dermatol Venereol Leprol 88(5):667–670. mediated antimicrobial photodynamic therapy
https://doi.org/10.25259/IJDVL_1175_ in clinical human studies: the state of the art.
2021 Photodiagn Photodyn Ther 2020(31):
2. Benov L (2015) Photodynamic therapy: cur- 101828. https://doi.org/10.1016/j.pdpdt.
rent status and future directions. Med Princ 2020.101828
Pract 24(Suppl 1):14–28. https://doi.org/ 4. Coates ARM, Hu Y, Holt J, Yeh P (2020)
10.1159/000362416 Antibiotic combination therapy against
56 Robert J. H. Hammond and Marianna Leite De Avellar
resistant bacterial infections: synergy, rejuve- (2006) Clinical and research applications of
nation and resistance reduction. Expert Rev photodynamic therapy in dermatology: experi-
Anti-Infect Ther 18(1):5–15. https://doi. ence of the Scottish PDT Centre. Lasers Surg
org/10.1080/14787210.2020.1705155 Med 38(5):403–416. https://doi.org/10.
5. Dai T, Huang YY, Hamblin MR (2009) Pho- 1002/lsm.20369
todynamic therapy for localized infections-- 14. O’Neill J (2016) Tackling drug-resistant infec-
state of the art. Photodiagn Photodyn Ther tions globally: final report and recommenda-
6(3–4):170–188. https://doi.org/10.1016/j. tions. Government of the United Kingdom.
pdpdt.2009.10.008 15. Ozlem-Caliskan S, Ilikci-Sagkan R, Karakas H,
6. Fu XJ, Fang Y, Yao M (2013) Antimicrobial Sever S, Yildirim C, Balikci M et al (2022)
photodynamic therapy for methicillin-resistant Efficacy of malachite green mediated photody-
Staphylococcus aureus infection. Biomed Res namic therapy on treatment of cutaneous
Int 2013:159157. https://doi.org/10.1155/ Leishmaniasis: in vitro study. Photodiagn
2013/159157 Photodyn Ther 40:103111. https://doi.org/
7. Gajdacs M, Abrok M, Lazar A, Burian K 10.1016/j.pdpdt.2022.103111
(2020) Increasing relevance of gram-positive 16. Raetz C, Guan Z, Ingram B, Six D, Song F,
cocci in urinary tract infections: a 10-year anal- Wang X et al (2009) Discovery of new biosyn-
ysis of their prevalence and resistance trends. thetic pathways: the lipid a story. J Lipid Res 50
Sci Rep 10(1):17658. https://doi.org/10. (Suppl):S103–S108. https://doi.org/10.
1038/s41598-020-74834-y 1194/jlr.R800060-JLR200. R800060-
8. Koeth LM, King A, Knight H, May J, Miller JLR200 [pii]
LA, Phillips I et al (2000) Comparison of 17. Stevens DL, Bisno AL, Chambers HF, Everett
cation-adjusted Mueller–Hinton broth with ED, Dellinger P, Goldstein EJC et al (2005)
Iso-Sensitest broth for the NCCLS broth Practice guidelines for the diagnosis and Man-
microdilution method. J Antimicrob Che- agement of Skin and Soft-Tissue Infections.
mother 46(3):369–376. https://doi.org/10. Clin Infect Dis 41(10):1373–1406. https://
1093/jac/46.3.369 doi.org/10.1086/497143
9. Kofler B, Romani A, Pritz C, Steinbichler TB, 18. Galata V, Laczny CC, Backes C, Hemmrich-
Schartinger VH, Riechelmann H et al (2018) Stanisak G, Schmolke S, Franke A, Meese E,
Photodynamic effect of methylene blue and Herrmann M, von Müller L, Plum A, Müller R,
low level laser radiation in head and neck squa- Stähler C, Posch AE, Keller A et al (2019)
mous cell carcinoma cell lines. Int J Mol Sci Integrating culture-based antibiotic resistance
1 9 ( 4 ) . h t t p s : // d o i . o r g / 1 0 . 3 3 9 0 / profiles with whole-genome sequencing data
ijms19041107 for 11,087 clinical isolates. Genom Proteom
10. Laxminarayan R, Duse A, Wattal C, Zaidi AK, Bioinform 17(2). https://doi.org/10.1016/j.
Wertheim HF, Sumpradit N et al (2013) Anti- gpb.2018.11.002
biotic resistance-the need for global solutions. 19. Varzandeh M, Mohammadinejad R,
Lancet Infect Dis 13(12):1057–1098. https:// Esmaeilzadeh-Salestani K, Dehshahri A,
doi.org/10.1016/s1473-3099(13)70318-9 Zarrabi A, Aghaei-Afshar A (2021) Photody-
11. Mantravadi PK, Kalesh KA, Dobson RCJ, namic therapy for leishmaniasis: recent
Hudson AO, Parthasarathy A (2019) The advances and future trends. Photodiagn Photo-
quest for novel antimicrobial compounds: dyn Ther 36:102609. https://doi.org/10.
emerging trends in research, development, 1016/j.pdpdt.2021.102609
and technologies. Antibiotics (Basel) 8(1). 20. Vasala A, Hytonen VP, Laitinen OH (2020)
h t t p s : // d o i . o r g / 1 0 . 3 3 9 0 / Modern tools for rapid diagnostics of antimi-
antibiotics8010008 crobial resistance. Front Cell Infect Microbiol
12. Moseley H, Allen JW, Ibbotson S, Lesar A, 10:308. https://doi.org/10.3389/fcimb.
McNeill A, Camacho-Lopez MA et al (2006) 2020.00308
Ambulatory photodynamic therapy: a new con- 21. Vital-Fujii DG, Baptista MS (2021) Progress in
cept in delivering photodynamic therapy. Br J the photodynamic therapy treatment of Leish-
Dermatol 154(4):747–750. https://doi.org/ maniasis. Braz J Med Biol Res 54(12):e11570.
10.1111/j.1365-2133.2006.07145.x h t t p s : // d o i . o r g / 1 0 . 1 5 9 0 / 1 4 1 4 -
13. Moseley H, Ibbotson S, Woods J, 431X2021e11570
Brancaleon L, Lesar A, Goodman C et al
Chapter 7
Abstract
The use of animal models is still widespread in science but there is a movement away from this manner of
experimentation. One option approved by the FDA for human-like studies is the hollow fiber bioreactor
(HFS). HFSs are highly controllable, self-contained systems that allow for the modeling of individual
tissues and disease phenotypes. Oxygen, drug concentration & half-life, and immune cell invasion are all
scalable to human and veterinary conditions using a HFS. There are drawbacks to the systems including cost
and contamination so the use of these systems must be carefully managed.
With these limitations in mind, the scope of the technology is great. Antimicrobial susceptibility testing
(AST) is possible with greater accuracy and clinical validity than classical in vitro techniques making minimal
inhibitory concentration (MIC) data generated on the bench more translatable to the clinic.
In this chapter, we will outline the background of the HFS and some typical uses.
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_7,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
57
58 Robert J. H. Hammond
B
1
D
E
F
A
Fig. 1 Schematic diagram of the hollow fiber system fully constructed. In our hollow fiber experiments, media
was supplied via input A under motive force from pump 1 (Masterflex C/L). The FiberCell Duet pump (2) then
draws the media along tube B to tube C and supplies it to the cartridge. The media passes through the
cartridge and flows down tube D to re-join the media in the central reservoir. At the same time pump 1 is also
drawing media out of the central reservoir (at the same rate as it is being introduced from tube A) along tube E
and into tube F and finally into the elimination reservoir where it is discarded
allowing for the removal of waste products [6, 7]. The fibers are
connected to a fluid source, such as a culture medium, which
provides a suitable environment for the growth of microorganisms.
1.2 Advantages of The principal advantage of hollow fiber bioreactors is that the
Using Hollow Fiber experimenter can control many of the critical parameters. Hollow
Bioreactors fiber bioreactors allow for precise control of the culture environ-
ment, including temperature, pH, drug delivery, and oxygen levels.
This means that the conditions can be reproduced accurately
throughout the course of the experiment and between
experiments.
No less important is the scalability of the design. Bioreactor
design allows for the cultivation of large quantities of microorgan-
isms, making it a suitable model for large-scale studies. 1011 cfu/
mL of some bacteria is possible in a volume of 20 mL of culture
medium, far outstripping conventional high-density
methodologies.
Using Hollow Fiber to Model Treatment of Antimicrobial-Resistant Organisms 59
1.3 Applications of Hollow fiber bioreactors have myriad uses and some examples are
Hollow Fiber detailed below.
Bioreactors in In a world where AMR is becoming increasingly important and
Antimicrobial research into it is picking up pace novel or updated methods of
Resistance Studies antibiotic susceptibility testing are crucial. Hollow fiber bioreactors
can be used to determine the susceptibility of microorganisms to
various antibiotics, providing valuable information for the selection
of appropriate treatments. In addition, they can be used to establish
treatment efficacy. Researchers can use the bioreactor to study the
efficacy of various antibiotics and combination therapies on resis-
tant microorganisms (see Chap. 5).
While considering AMR it is vital to understand how a patho-
gen resists the antibiotic in question. Understanding the mechan-
isms of resistance allows researchers to gain insights into potential
strategies to overcome pathogen resistance or develop new
therapies.
1.4 Limitations and Hollow fiber bioreactors are relatively complex. At their fullest
Considerations extent they can contain meters of tubing and several bottles of
media. This complexity of design means that understanding and
successfully running a hollow fiber bioreactor requires a certain
level of technical expertise to set up and use. It also means that
there is considerable opportunity for contamination if mitigation
strategies are not put into place before the experiment begins.
The expense of the hollow fiber can also be prohibitive. The
consumable element is the hollow fiber cartridge itself and these
can cost hundreds of pounds each. The other aspects of the set up
are more long-lasting but are still costly. This can be considered a
significant restrictive factor as unfortunately its accessibility will be
limited to some researchers.
While the FDA has approved the hollow fiber bioreactor for
in vivo-like studies it may not fully capture the complexities of the
human body limiting its scope as a research tool. This, in turn, leads
to limitations in the conclusions that can be drawn from any results
generated when compared with in vivo data and especially when
comparing with first-in-human trials.
60 Robert J. H. Hammond
1.5 Initial The use of hollow fiber bioreactors provides a valuable tool for
Conclusions researchers studying the treatment of antimicrobial-resistant organ-
isms. The controlled environment and scalability of the bioreactor
make it a useful model for studies on antibiotic susceptibility,
treatment efficacy, and mechanisms of resistance. However, it is
important to keep in mind the limitations and considerations
when interpreting the results of studies using this model.
In the below example of an AST testing method, the hollow
fiber assay was used to simulate in vivo conditions without the use
of human or animal test subjects. As a highly manageable test
system, it allows control of O2 concentration, drug half-life (T1/
2), maximal concentration of drug (Cmax), and time over which
maximal concentration of drug is reached (Tmax) as well as other
factors. This makes it ideal for testing novel antibiotics or antibiotic
combination therapies.
2 Materials
3 Methods
3.1 Hollow Fiber The overall shape of the rig used is influenced by the experiment
Bioreactor being undertaken. In the below example an AST test is being done
wherein components required are as in the material section above.
Briefly: C3008 cartridge, media “feed” bottle, empty “waste” bot-
tle, connectors, drug infusing pump (syringe driver), 3-way stop-
cock, silicone tubing.
1. Prepare Hollow fiber cartridges as per the manufacturer’s
instructions and attached to the overall rig, aseptically, as is
shown in Fig. 1 (see Note 7).
2. The attachment is best undertaken in a sterile cabinet with all
tubing ends that are going to meet being sprayed liberally with
70% ethanol.
3. Introduce bacteria into the cartridge via the sampling ports,
indicated in Fig. 1, and allowed to settle for 1 h (see Note 8).
4. After this time, activate the Duet pump to ensure flow through
the cartridge is smooth and to eliminate bubbles.
5. Antibiotics are introduced into the system at the red X (see
Fig. 1) using a 3-way stopcock and a syringe driver (Aladdin
Programmable Syringe Pump AL 1000-220).
6. Set the pump to deliver drug such that the final concentration
at the end of dosing was equal to the Cmax of the drug over the
time Tmax.
7. Stop pump 1 for the duration of drug delivery at a predeter-
mined time (Tmax).
Calculate the specific Cmax and Tmax parameters for each
drug and drug combination used based on the drug’s half-life
in vivo, the flow rate of the system, and other factors [7].
62 Robert J. H. Hammond
3.2 Sampling 1. Sample the bacteria at regular intervals throughout the experi-
mental period (7 days) to address the pharmacology of the
studied drug.
2. Select the sampling schedule that best assesses changes in bac-
terial numbers in response to the administered drugs and mea-
sure the associated drug concentration.
3. In order to sample the bacteria from the cartridge Switch flow
pump (FiberCell Duet (2)) is switched off temporarily, 5 mL of
medium is removed as described by the manufacturer
(FiberCell).
4. Remove the filled syringe from the assembly. Spray the Luer
lock with ethanol and add a new sterile syringe to the vacated
Luer lock.
5. Restart the flow pump and continue the experiment until the
next sampling interval.
3.3 Sample Depending on the initial conditions for the experiment and the
Examination reason for the experiment, the samples produced and the generated
data will vary.
In the exemplar above an AST test was attempted against a
resistant organism. As this is an in vivo-like system the results
generated here might be different to those found in a static assay
like a broth microdilution, disc diffusion, or static time kill assay.
Examining these samples requires understanding of the patho-
gen used and an understanding of common culture techniques.
Where sampling has taken place throughout a dosing regimen
bacteria might be in different growth phases, stress responses, or
stages of death. These are valuable information and can be collected
in a variety of ways; growth kinetic experiments, RNA sequencing,
and viability assays to name a few.
Generally the most compelling piece of information, and the
most common, is the simplest; are the cells alive or not? This data
can, again, be extracted in a variety of ways such as viability plating,
colorimetric assays, and DNA sequencing.
Using Hollow Fiber to Model Treatment of Antimicrobial-Resistant Organisms 63
4 Notes
References
1. Gajdacs M, Abrok M, Lazar A, Burian K 6. Kloprogge F, Hammond R, Copas A, Gillespie
(2020) Increasing relevance of gram-positive SH, Della PO (2019) Can phenotypic data
cocci in urinary tract infections: a 10-year anal- complement our understanding of antimyco-
ysis of their prevalence and resistance trends. bacterial effects for drug combinations? J Anti-
Sci Rep 10(1):17658. https://doi.org/10. microb Chemother 2019:3530. https://doi.
1038/s41598-020-74834-y org/10.1093/jac/dkz369
2. O’Neill J (2016) Tackling drug-resistant infec- 7. Kloprogge F, Hammond R, Kipper K, Gillespie
tions globally: final report and recommenda- SH, Della PO (2019) Mimicking in-vivo expo-
tions. Government of the United Kingdom sures to drug combinations in-vitro: anti-
3. Galata V, Laczny CC, Backes C, Hemmrich- tuberculosis drugs in lung lesions and the hol-
Stanisak G, Schmolke S, Franke A, Meese E, low fiber model of infection. Sci Rep 9(1):
Herrmann M, von Müller L, Plum A, Müller R 13228. https://doi.org/10.1038/s41598-
(2019) Integrating culture-based antibiotic 019-49556-5
resistance profiles with whole-genome 8. Storm MP, Sorrell I, Shipley R, Regan S,
sequencing data for 11,087 clinical isolates. Luetchford KA, Sathish J et al (2016) Hollow
Genom Proteom Bioinform 17(2):169–182. fiber bioreactors for in vivo-like mammalian
https://doi.org/10.1016/j.gpb.2018.11.002 tissue culture. JoVE 111:e53431
4. Vasala A, Hytonen VP, Laitinen OH (2020) 9. Romero K, Clay R, Hanna D (2015) Strategic
Modern tools for rapid diagnostics of antimi- regulatory evaluation and endorsement of the
crobial resistance. Front Cell Infect Microbiol hollow fiber tuberculosis system as a novel
10:308. https://doi.org/10.3389/fcimb. drug development tool. Clin Infect Dis 61
2020.00308 (suppl_1):S5–S9
5. Deshpande D, Srivastava S, Nuermberger E, 10. Koeth LM, King A, Knight H, May J, Miller
Pasipanodya JG, Swaminathan S, Gumbo T LA, Phillips I et al (2000) Comparison of
(2016) Concentration-dependent synergy and cation-adjusted Mueller–Hinton broth with
antagonism of linezolid and moxifloxacin in Iso-Sensitest broth for the NCCLS broth
the treatment of childhood tuberculosis: the microdilution method. J Antimicrob Che-
dynamic duo. Clin Infect Dis 32016:S88–S94 mother 46(3):369–376. https://doi.org/10.
1093/jac/46.3.369
Chapter 8
Abstract
Pyrazinamide (PZA) is a key component of chemotherapy for the treatment of drug-susceptible tuberculo-
sis (TB) and is likely to continue to be included in new drug combinations. Potentiation of PZA could be
used to reduce the emergence of resistance, shorten treatment times, and lead to a reduction in the quantity
of PZA consumed by patients, thereby reducing the toxic effects. Acidified medium is required for the
activity of PZA against Mycobacterium tuberculosis. In vitro assessments of pyrazinamide activity are often
avoided because of the lack of standardization, which has led to a lack of effective in vitro tools for assessing
and/or enhancing PZA activity.
We have developed and optimized a novel, robust, and reproducible, microtiter plate assay, that centers
around acidity levels that are low enough for PZA activity. The assay can be applied to the evaluation of
novel compounds for the identification of potentiators that enhance PZA activity. In this assay, potentiation
of PZA is demonstrated to be statistically significant with the addition of rifampicin (RIF), which can,
therefore, be used as a positive control. Conversely, norfloxacin demonstrates no potentiating activity with
PZA and can be used as a negative control. The method, and the associated considerations, described here,
can be adapted in the search for potentiators of other antimicrobials.
Key words Mycobacterium tuberculosis, Pyrazinamide, Potentiation, Microtiter plate, Acidic pH,
Chequerboard, Rifampicin, Norfloxacin, Highest single agent
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_8,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
65
66 Christopher William Moon et al.
upon have often been avoided because of the complexities and lack
of standardization of in vitro assay conditions, which has led to a
lack of effective in vitro tools for assessing/enhancing PZA activity.
The optimal pH for the activity of PZA in vitro has previously been
shown to be pH 5.6 [14]. Although acidification of the growth
medium is considered to be important for PZA activity, previous
studies have demonstrated that in certain situations, PZA is active
at a neutral pH [15]. Further to this, metabolic potentiation of
PZA activity has been discovered independently of environmental
pH [16]. However, a key component of the assay, described here, is
the growth of M. tuberculosis in a defined, suitably buffered
medium, at pH 5.5, followed by exposure to drug, under the
same conditions in the microtiter plate. This is designed to reflect
the acidic conditions that the organism will be encountering intra-
cellularly in the macrophage and in the granuloma, during disease;
inflammatory lung lesions have shown to be pH 5.5—pH 6.0
[17]. Variations in in vitro parameters such as inoculum size and
medium components can also alter the activity of PZA [7] and
these have been optimized in our method. In this assay, potentia-
tion of PZA is demonstrated to be statistically significant with the
addition of rifampicin (RIF), which is used as a positive control.
Conversely, norfloxacin (NOR) demonstrates no potentiating
activity with PZA and serves as a negative control.
In summary, robust and reproducible in vitro assay conditions
are required for the discovery of potentiators that can boost PZA
activity. We have successfully addressed these issues and have devel-
oped an in vitro microtiter plate chequerboard assay that can be
used to find compounds that will potentiate PZA activity against
M. tuberculosis. The method, and the associated considerations,
described here, can be adapted in the search for potentiators of
other antimicrobials.
2 Materials
2.1 Growth Medium 1. Middlebrook 7H9 medium supplemented with oleic acid, dex-
trose, and catalase (OADC) enrichment (referred to as “Mid-
dlebrook 7H9”). The ingredients required for 1 L of
Middlebrook 7H9: Middlebrook 7H9 broth powder (4.7 g),
glycerol (5 mL), Tween-80 (2 mL), and oleic albumin dextrose
catalase (OADC; (Beckman Dickinson, UK—212240))
(100 mL). Add the Middlebrook 7H9 broth powder to
500 mL of distilled water, and then add the glycerol and
Tween-80. Stir the solution to dissolve all the ingredients and
make the volume up to 900 mL using distilled water. Autoclave
this medium at 121 °C, 15 psi for 15 min. Once cooled to
room temperature, add 100 mL of OADC. Store at room
temperature. For best results, use freshly made medium (see
Note 1).
68 Christopher William Moon et al.
2.4 Equipment 1. Erlenmeyer cell culture vented flasks, 125 mL, sterile (Corning,
USA—431143).
2. Shaking incubator, 37 °C, 200 rpm.
3. Screw cap, O-ring, 2 mL tubes.
4. UV/vis spectrophotometer.
5. Cuvettes.
6. Reagent reservoirs.
7. Nunc 96-well, black plates (Thermo Fisher Scientific, UK—
237105) (see Note 3).
8. SealPlate film (Excel Scientific).
9. MicroAmp Optical Adhesive film (Thermo Fisher Scientific,
UK—4311971).
Enhancing Pyrazinamide Activity Against Mycobacterium tuberculosis 69
3 Methods
3.1 Preparation of M. 1. Derive mycobacterial stocks, used for the inoculations, from
tuberculosis Inoculum the growth of M. tuberculosis H37Rv, in a defined and standar-
dized way. The preferred method is by continuous culture in a
chemostat as described by Bacon et al. in 2004 [18] and Bacon
and Hatch in 2009 [19] (see Note 6).
2. Store culture stocks as 1 mL volumes at -70 °C, until use (see
Note 7).
3.3 Growth of M. 1. Sample the cultures, described in Subheading 3.2, for optical
tuberculosis in density measurements. Fix the cells by the addition of 4%
Acidified Middlebrook formaldehyde (final concentration v/v) to a tenfold dilution
7H9 for 7 Days of the bacterial culture, in PBS pH 7.4, in a screw cap tube with
an O-ring. Invert the tube and leave for 30 min, to fix the
bacteria (see Note 9).
2. Place the formaldehyde-fixed mycobacterial suspension into a
plastic cuvette.
3. Read the optical density at 540 nm (OD540nm) against water as
the blank.
3.4 Growth of M. 1. Prepare six 96-well microtiter plates; three biological repeat
tuberculosis in plates for the positive control drug (RIF) and three biological
Acidified Middlebrook repeat plates for the negative control drug (NOR).
7H9 for 10 Days 2. Add 100 μL of acidified Middlebrook 7H9 to all wells of rows
(Microtiter Plates) in A-G and wells H11 and H12 of each 96-well plate.
the Presence of 3. Add 190.44 μL of acidified Middlebrook 7H9 to wells
Antibiotic H1-H10 of each 96-well plate.
Combinations
4. Add 9.56 μL PZA to wells H1-H10 in all plates from a stock
3.4.1 Preparation of PZA- concentration of 20 mg/mL prepared in acidified Middleb-
RIF and PZA-NOR rook 7H9 to achieve a concentration of 1600 μg/mL in a final
Combinations in 96-Well assay volume of 110 μL per well.
Plates (Chequerboard) 5. Perform a twofold serial dilution of PZA across each plate, by
transferring and mixing 100 μL from wells H1-H10 through to
wells B1-10 in all plates, to achieve concentrations of
1600–25 μg/mL in rows B-H and columns 1–10 (see Fig. 2).
6. To three of the 96-well plates create a RIF twofold dilution
series from a stock concentration of 0.5 mg/mL prepared in
dimethyl sulfoxide (DMSO) to achieve concentrations of
50–0.195 ng/mL in rows A-H and columns 2–10 in a micro-
titer plate (see Fig. 2a).
7. To the remaining three 96-well plates create a norfloxacin
twofold dilution series from a stock concentration of 2 mg/
mL prepared in DMSO to achieve concentrations of
32–0.125 μg/mL in rows A-H and columns 2–10 in a microti-
ter plate (see Fig. 2b) (see Note 12).
(a) Addition of PZA and RIF, or NOR on opposing axis in the
microtiter plate creates a dose matrix of PZA and RIF or
NOR with eight concentrations of PZA and ten concen-
trations of RIF or NOR. No antibiotics are added to well
A1.
8. To columns 1–10, rows A-H add sufficient DMSO to raise and
normalize the DMSO concentration of each well to 3% (v/v)
(see Note 13).
9. To column 11, rows A-H add 3.4 μL DMSO to achieve a final
concentration of 3% (v/v) (see Note 14).
10. No antibiotic or DMSO additions are made to column
12 (medium-only control) (see Note 15).
Enhancing Pyrazinamide Activity Against Mycobacterium tuberculosis 71
Fig. 2 Example plate layout for creating a 10 × 8 dose matrix of rifampicin and
pyrazinamide (a), and norfloxacin and pyrazinamide (b). Concentrations of
rifampicin are 50, 25, 12.5, 6.25, 3.13, 1.56, 0.781, 0.391, and 0.195 ng/mL,
concentrations of norfloxacin are 32, 16, 8, 4, 2, 1, 0.5, 0.25, and 0.125 μg/mL,
concentrations of pyrazinamide are 1600, 800, 400, 200, 100, 50, and 25 μg/mL
72 Christopher William Moon et al.
3.4.2 Inoculation of 1. Sample the cultures that have been growing for 7 days in
Antibiotic Dose Matrix acidified medium (described in Subheading 3.3) for optical
Microtiter Plates with M. density measurements. Fix the cells by the addition of 4%
tuberculosis formaldehyde (final concentration v/v) to a tenfold dilution
of the bacterial culture, in PBS pH 7.4, in a screw cap tube with
an O-ring. Invert the tube and leave for 30 min, to fix the
bacteria.
2. Place cell suspension into a plastic cuvette.
3. Read the optical density at 540 nm (OD540nm) against water.
4. If cultures meet growth requirements of a minimum OD540nm
of 1.0 (see Note 16), dilute the three cultures to an OD540nm of
0.55, using acidified Middlebrook 7H9 (see Notes 17 and 18).
5. Add 10 μL of the prepared M. tuberculosis culture (at 0.55
OD540nm) to columns 1–11 of all six 96-well plates (see
Notes 15 and 19).
6. Add 10 μL of acidified Middlebrook 7H9 to column 12 of the
six 96-well plates (see Note 15).
7. Apply a plate seal (SealPlate film, Excel Scientific) to all plates
(see Note 20).
8. Incubate the six 96-well plates in a contained shaking incubator
for 10 days, at 37 °C, orbitally shaking at 200 rpm (see Note 8)
(see Fig. 1).
3.5 Endpoint 1. Remove the plate seals from the microtiter plates and discard.
Detection/Final 2. Add 10 μL of resazurin (that has been adjusted to pH 8.2 as
Viability Assessment described in Subheading 2.2) to all wells of the microtiter
plates and incubate at room temperature, for 6 h (see Note 21).
3. Add a MicroAmp Optical Adhesive film (Thermo Fisher Scien-
tific) plate seal, to each plate, ensuring a tight, complete seal (see
Notes 22 and 23).
4. Image the microtiter plates using a microplate reader with fluo-
rescent detection (Ex: 525 nm Em: 580–640 nm) (see Note 24).
3.6 Data Analysis 1. Biological repeat averaging—Calculate the average of the resor-
ufin fluorescence data collected from each biological replicate
in each well of the microtiter plates.
2. Baseline correction—Calculate the average of column
12 (medium-only control). Subtract this value from all other
values (columns 1–11), to generate baseline-corrected data.
3. Normalization—Calculate the average of column
11 (M. tuberculosis with 3% DMSO addition). Divide all
baseline-corrected values by this average value for column
11, then multiply by 100 (columns 1–10), to generate normal-
ized data compared to M. tuberculosis growth with no antibi-
otic addition (see Note 25).
Enhancing Pyrazinamide Activity Against Mycobacterium tuberculosis 73
Fig. 3 Example baseline corrected and normalized pyrazinamide and rifampicin dose matrix data (a) and
pyrazinamide and norfloxacin dose matrix data (c). Data values are colored so that growth of 100% that of the
3% DMSO control M. tuberculosis growth (column 11 average) is green and 0% is red. Example of highest
single agent determination of antibiotic combination synergy of pyrazinamide and rifampicin (b) and pyrazi-
namide and norfloxacin (d). Synergistic effects (higher positive values) are colored blue and antagonistic
effects (higher negative values) are colored red
4 Notes
each dose matrix point using Eq. (1). HSA values greater than
10 represent synergistic combinations, values between 10 and
-10 represent additive combinations and values less than -10
represent antagonistic combinations. HSA is only defined on a
dose-by-dose basis with no model fit, it is sensitive to the dose
range selected.
Acknowledgments
References
1. British T (1976) Short-course chemotherapy 9. Ibrahim M et al (2007) Synergistic activity of
in pulmonary tuberculosis. Lancet (London, R207910 combined with pyrazinamide against
England) 2:1102–1104 murine tuberculosis. Antimicrob Agents Che-
2. Mitchison DA (1985) The action of antituber- mother 51:1011–1015
culosis drugs in short-course chemotherapy. 10. Shi W et al (2011) Pyrazinamide inhibits trans-
Tubercle 66:219. https://doi.org/10.1016/ translation in Mycobacterium tuberculosis. Sci-
0041-3879(85)90040-6 ence 333:1630–1632
3. Wade MM, Zhang Y (2006) Effects of weak 11. Keiler KC (2008) Biology of trans-translation.
acids, UV and proton motive force inhibitors Ann Rev Microbiol 62:133–151
on pyrazinamide activity against Mycobacte- 12. Rodionova IA et al (2014) Metabolic and bac-
rium tuberculosis in vitro. J Antimicrob Che- tericidal effects of targeted suppression of
mother 58:936–941 NadD and NadE enzymes in mycobacteria.
4. Zhang Y, Shi W, Zhang W, Mitchison D MBio 5:10-1128
(2014) Mechanisms of pyrazinamide action 13. Cook GM, Greening C, Hards K, Berney M
and resistance. Microbiol Spectr 2:10–1128 (2014) Energetics of pathogenic bacteria and
5. Wade MM, Zhang Y (2004) Anaerobic incuba- opportunities for drug development. Adv
tion conditions enhance pyrazinamide activity Microb Physiol 65:1–62
against Mycobacterium tuberculosis. J Med 14. Gumbo T, Dona CSWS, Meek C, Leff R
Microbiol 53:769–773 (2009) Pharmacokinetics-pharmacodynamics
6. Zhang Y, Wade MM, Scorpio A, Zhang H, Sun of pyrazinamide in a novel in vitro model of
Z (2003) Mode of action of pyrazinamide: tuberculosis for sterilizing effect: a paradigm
disruption of Mycobacterium tuberculosis mem- for faster assessment of new antituberculosis
brane transport and energetics by pyrazinoic drugs. Antimicrob Agents Chemother 53:
acid. J Antimicrob Chemother 52:790–795 3197–3204
7. Zhang Y, Permar S, Sun Z (2002) Conditions 15. Thiede JM et al (2022) Pyrazinamide suscepti-
that may affect the results of susceptibility test- bility is driven by activation of the SigE-
ing of Mycobacterium tuberculosis to pyrazina- dependent cell envelope stress response in
mide. J Med Microbiol 51:42–49 Mycobacterium tuberculosis. MBio 13:
8. Sbarbaro JA, Iseman MD, Crowle AJ (1992) e0043921
The combined effect of rifampin and pyrazina- 16. den Hertog AL et al (2016) Pyrazinamide is
mide within the human macrophage. Am Rev active against Mycobacterium tuberculosis cul-
Respir Dis 146:1448–1451 tures at neutral pH and low temperature.
Enhancing Pyrazinamide Activity Against Mycobacterium tuberculosis 77
Antimicrob Agents Chemother 60:4956. 18. Bacon J et al (2004) The influence of reduced
https://doi.org/10.1128/aac.00654-16 oxygen availability on pathogenicity and gene
17. Lanoix J-P et al (2016) Selective inactivity of expression in Mycobacterium tuberculosis.
pyrazinamide against tuberculosis in C3HeB/ Tuberculosis 84:205–217
FeJ mice is best explained by neutral pH of 19. Bacon J, Hatch KA (2009) Continuous culture
Caseum. Antimicrob Agents Chemother 60: of mycobacteria. Methods Mol Biol 465:153–
735–743 171
Chapter 9
Abstract
Mathematical models have been used to study the spread of infectious diseases from person to person. More
recently studies are developing within-host modeling which provides an understanding of how pathogens—
bacteria, fungi, parasites, or viruses—develop, spread, and evolve inside a single individual and their
interaction with the host’s immune system.
Such models have the potential to provide a more detailed and complete description of the pathogenesis
of diseases within-host and identify other influencing factors that may not be detected otherwise. Mathe-
matical models can be used to aid understanding of the global antibiotic resistance (ABR) crisis and identify
new ways of combating this threat.
ABR occurs when bacteria respond to random or selective pressures and adapt to new environments
through the acquisition of new genetic traits. This is usually through the acquisition of a piece of DNA from
other bacteria, a process called horizontal gene transfer (HGT), the modification of a piece of DNA within a
bacterium, or through. Bacteria have evolved mechanisms that enable them to respond to environmental
threats by mutation, and horizontal gene transfer (HGT): conjugation; transduction; and transformation. A
frequent mechanism of HGT responsible for spreading antibiotic resistance on the global scale is conjuga-
tion, as it allows the direct transfer of mobile genetic elements (MGEs). Although there are several MGEs,
the most important MGEs which promote the development and rapid spread of antimicrobial resistance
genes in bacterial populations are plasmids and transposons. Each of the resistance-spread-mechanisms
mentioned above can be modeled allowing us to understand the process better and to define strategies to
reduce resistance.
Key words Mathematical modeling, Antibiotic resistance, Within-host modeling, Differential equa-
tions, Prediction
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_9,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
79
80 Aminat Yetunde Saula et al.
2 Within-Host Modeling
Fig. 1 Schematic illustration of a model of resistance to one antibiotic. This shows the evolution of bacteria
population following a logistic growth from the initial state of susceptibility to an antibiotic to the decrease in
the population of susceptible bacteria attributed to the actions of the immune cells (γ), the bactericidal effect
of the antibiotic (α) on the susceptible population, and the host’s carrying capacity K indicating scarcity of
resources which limits growth rate (i.e., vN/K). The subsequent growth of resistant bacteria follows fitness
cost, with a decrease in the population of resistant bacteria caused by γ and the carrying capacity K (i.e., v1N/
K ). The antibiotic does not affect resistant bacteria. The movement from susceptible to resistant compartment
(highlighted with a blue arrow) describes the proportion q of resistant bacteria to the antibiotic, that emerge
during the reproduction of susceptible bacteria
Fig. 2 Schematic illustration of a model of resistance to one antibiotic. This shows the evolution of bacteria
population following a logistic growth from the initial state of susceptibility to antibiotics to the decrease in the
population of susceptible bacteria due to the actions of the immune cells (γ), the bactericidal effect of both
antibiotics (αi, i = 1, 2) on the susceptible population and the host’s carrying capacity K (i.e., vN K ). The
subsequent growth of resistant bacteria populations Ri, (i = 1, 2, 3) indicates the acquisition of resistance to
drug 1, drug 2, and both drugs respectively with growth rates vi, (i = 1, 2, 3) following fitness cost. The
reduction in the populations of resistant bacteria is due to γ, the bactericidal effect of drug 2 on R1, the
bactericidal effect of drug 1 on R2, and the carrying capacity K (i.e., viN/K ). Both antibiotics do not affect R3
population due to acquired resistance to both drugs. The movement between compartments (highlighted with
blue arrows) describes the proportion qi, (i = 1, 2) of resistant bacteria to drug 1, drug 2, or both drugs that
emerge during the reproduction of susceptible bacteria
References
1. Awofisayo-Okuyelu A, Pratt A, McCarthy N, 16. Huddleston JR (2014) Horizontal gene trans-
Hall I (2019) Within-host mathematical mod- fer in the human gastrointestinal tract: poten-
elling of the incubation period of salmonella tial spread of antibiotic resistance genes. Infect
typhi. R Soc Open Sci 6:182143 Drug Resist 7:167
2. Martinez J, Baquero F (2000) Mutation fre- 17. McCarthy AJ, Loeffler A, Witney AA, Gould
quencies and antibiotic resistance. Antimicrob KA, Lloyd DH, Lindsay JA (2014) Extensive
Agents Chemother 44:1771–1777 horizontal gene transfer during staphylococcus
3. Boto L, Mart´ınez JL. (2011) Ecological and aureus co-colonization in vivo. Genome Biol
temporal constraints in the evolution of bacte- Evol 6:2697–2708
rial genomes. Genes 2:804–828 18. Lipsitch M, Levin BR (1997) The population
4. Pang Z, Raudonis R, Glick BR, Lin TJ, Cheng dynamics of antimicrobial chemotherapy. Anti-
Z (2019) Antibiotic resistance in pseudomonas microb Agents Chemother 41:363–373
aeruginosa: mechanisms and alternative thera- 19. Holt RD, Barfield M (2006) Within-host path-
peutic strategies. Biotechnol Adv 37:177–192 ogen dynamics: some ecological and evolution-
5. Chevallereau A, Meaden S, van Houte S, Wes- ary consequences of transients, dispersal mode,
tra ER, Rollie C (2019) The effect of bacterial and within-host spatial heterogeneity. In: Dis-
mutation rate on the evolution of crispr-cas ease evolution: models, concepts, and data ana-
adaptive immunity. Philos Trans R Soc B 374: lyses, vol 71, pp 45–66
20180094 20. Drlica K (2003) The mutant selection window
6. Foster PL (2005) Stress responses and genetic and antimicrobial resistance. J Antimicrob
variation in bacteria. Mutat Res/Fundam Mol Chemother 52:11–17
Mech Mutagen 569:3–11 21. Thompson JN (1998) Rapid evolution as an
7. Burmeister AR (2015) Horizontal gene trans- ecological process. Trends Ecol Evol 13:329–
fer. Evol Med Public Health 2015:193 332
8. Bushman F (2002) Lateral DNA transfer. Cold 22. Thomas CM, Nielsen KM (2005) Mechanisms
Spring Harbor Laboratory Press of, and barriers to, horizontal gene transfer
9. Von Wintersdorff CJ, Penders J, Van Niekerk between bacteria. Nat Rev Microbiol 3:711–
JM, Mills ND, Majumder S, Van Alphen LB 721
et al (2016) Dissemination of antimicrobial 23. Culyba MJ, Van Tyne D (2021) Bacterial evo-
resistance in microbial ecosystems through lution during human infection: Adapt and live
horizontal gene transfer. Front Microbiol 7: or adapt and die. PLoS Pathog 17:e1009872
173 24. Alizon S, Luciani F, Regoes RR (2011) Epide-
10. Andam CP, Fournier GP, Gogarten JP (2011) miological and clinical consequences of within-
Multilevel populations and the evolution of host evolution. Trends Microbiol 19:24–32
antibiotic resistance through horizontal gene 25. Kaplan G, Post FA, Moreira AL,
transfer. FEMS Microbiol Rev 35:756–767 Wainwright H, Kreiswirth BN, Tanverdi M
11. Clewell DB (2013) Bacterial conjugation. et al (2003) Mycobacterium tuberculosis
Springer Science & Business Media growth at the cavity surface: a microenviron-
12. Munita JM, Arias CA (2016) Mechanisms of ment with failedimmunity. Infect Immun 71:
antibiotic resistance. Microbiol Spectr 4:4–2 7099–7108
13. Coombs D, Gilchrist MA, Ball CL (2007) 26. Read AF, Huijben S (2009) Perspective: evolu-
Evaluating the importance of within-and tionary biology and the avoidance of antimi-
between-host selection pressures on the evolu- crobial resistance. Evol Appl 2:40–51
tion of chronic pathogens. Theor Popul Biol 27. Boni MF, Feldman MW (2005) Evolution o
72:576–591 antibiotic resistance by human and bacterial
14. Andersson DI, Patin SM, Nilsson AI, Kugel- neche construction. Evolution 59:477–491
berg E (2007) The biological cost of antibiotic 28. Cobey S, Baskerville EB, Colijn C, Hanage W,
resistance. In: Enzyme-mediated resistance to Fraser C, Lipsitch M (2017) Host population
antibiotics: mechanisms, dissemination, and structure and treatment frequency maintain
prospects for inhibition, pp 339–348 balancing selection on drug resistance. J R
15. SchulzzurWiesch P, Engelstadter J, Bonhoeffer Soc Interface 14:20170295
S (2010) Compensation of fitness costs and 29. Blanquart F, Lehtinen S, Lipsitch M, Fraser C
reversibility of antibiotic resistance mutations. (2018) The evolution of antibiotic resistance in
Antimicrob Agents Chemother 54:2085–2095 a structured host population. J R Soc Interface
15:20180040
90 Aminat Yetunde Saula et al.
55. Ganguli S, Gammack D, Kirschner DE (2005) 57. Hao W, Schlesinger LS, Friedman A (2016)
A metapopulation model of granuloma forma- Modeling granulomas in response to infection
tion in the lung during infection with myco- in the lung. PLoS One 11:e0148738
bacterium tuberculosis. Math Biosci Eng 2:535 58. Pitcher MJ, Bowness R, Dobson S, Eftimie R,
56. Golubchik T, Batty EM, Miller RR, Farr H, Gillespie SH (2020) Modelling the effects of
Young BC, Larner-Svensson H et al (2013) environmental heterogeneity within the lung
Within- host evolution of staphylococcus on the tuberculosis life-cycle. J Theor Biol
aureus during asymptomatic carriage. PLoS 506:110381
One 8:e61319
Chapter 10
Abstract
To model complex systems, individual-based models (IBMs), sometimes called “agent-based models”
(ABMs), describe a simplification of the system through an adequate representation of the elements.
IBMs simulate the actions and interaction of discrete individuals/agents within a system in order to discover
the pattern of behavior that comes from these interactions. Examples of individuals/agents in biological
systems are individual immune cells and bacteria that act independently with their own unique attributes
defined by behavioral rules. In IBMs, each of these agents resides in a spatial environment and interactions
are guided by predefined rules. These rules are often simple and can be easily implemented. It is expected
that following the interaction guided by these rules we will have a better understanding of agent-agent
interaction as well as agent-environment interaction. Stochasticity described by probability distributions
must be accounted for. Events that seldom occur such as the accumulation of rare mutations can be easily
modeled.
Thus, IBMs are able to track the behavior of each individual/agent within the model while also obtaining
information on the results of their collective behaviors. The influence of impact of one agent with another
can be captured, thus allowing a full representation of both direct and indirect causation on the aggregate
results. This means that important new insights can be gained and hypotheses tested.
Key words Mathematical modeling, Antibiotic resistance, Agent-based modeling, Differential equa-
tions, Prediction
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_10,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
93
94 Aminat Yetunde Saula et al.
2 Constructing IBMs
2.1 Defining the The first step in IBM construction is to clearly define the individuals
Environment and the that make up the system, the environment which serves as a habitat
Individual Elements for the individuals, and the spatial and temporal extent of the
model. Most IBMs, especially when modeling within-host antibi-
otic resistance, are spatially explicit. This implies that they capture
behavioral shifts of individuals under varying conditions in their
local environment (such as a section of lung tissue, a section of
intestinal epithelium, etc.) at different spatial and temporal scales.
In the following steps, more information will be provided about the
behavioral changes that could occur. For now, we shall focus on the
many questions that revolve around our definition of the environ-
ment and the processes or actions that could be carried out by the
individuals within our defined environment.
When defining an environment, different types of space can be
adopted during the construction of an IBM. This could be through
the use of arbitrary graphs, discrete or continuous spaces. A discrete
(grid-like) space is commonly utilized for within-host IBMs due to
its simplicity, visual appeal, and its easy tracking of the position of
cells over time. With a grid-like space, it is easy to observe the
removal, replication/recruitment of agents distributed over cells,
as well as the events that occur at each time step—which is an
important attribute when modeling the evolution of bacteria.
Thus, we present some of these questions as guidelines for consid-
eration when defining the individuals that make up the system and
98 Aminat Yetunde Saula et al.
Fig. 1 Simple illustration of possible environment and individuals for a within-host model of antibiotic
resistance
Fig. 2 An example of possible initial state (and traits) variables. Green box represents the individuals, blue
indicates the state variables, and red indicates the trait variables
2.3 Defining We have defined the initial states and traits of the individuals in the
Individual Activities/ second step, as a basis for what contributes to state updates. The
Actions and traits of the individual entities and their behavioral rules organize
Environmental their actions and how they interact. In step 2, we have not defined
Attributes Influencing the activities that could occur and the outcome of interactions
State (and Traits) between individuals. However, it is important to note that indivi-
Updates duals have different objectives. Hence, we need to carefully think
through the actions that are important to address these objectives,
again, ignoring irrelevant ones. Hence, for each type of individual, a
rule must be defined that allows the implementation of possible
activities that could be carried out in a particular situation. Usually,
these activities influence a change in the state of the individuals as
well as the environment. Hence, the actions to be incorporated and
their implications must be clearly stated.
In some cases, individuals may make decisions in response to
some form of stimuli, (e.g., the initiation of cell recruitment
through the production of cytokines in response to bacterial infec-
tion, movement of immune cells toward the site of infection,
induced by the high concentration of chemokine molecules)
and/or other criteria (e.g., random/pseudo-random migration).
These activities or actions usually contribute to the overall outcome
observed and may take different forms. Some of these activities
occur randomly, in a hierarchy, concurrently, or could be fixed.
Hence, it becomes important to assign rates or probabilities,
where appropriate, to the activities or actions carried out. Examples
of these actions or decisions taken by individuals include the
recruitment of T cells at a particular density of bacteria or the supply
of antibiotics at a particular time step to inhibit bacterial growth or
induce mortality ( fixed); the death of bacteria after the movement
of active macrophages toward the infection site (hierarchy); the
migration of T cells toward infected macrophages to induce the
mortality of macrophages and internalized bacteria at the same time
(concurrent); and the phagocytosis of bacteria upon contact with
inactive macrophages (random). In the examples given, following
these actions, a state and trait update will occur for bacteria and
macrophages. While the interaction between individuals affects the
environment, environmental attributes such as ventilation, perfu-
sion, and drainage could also influence a state update of individuals,
thus, affecting the overall outcome of our result [29]. A good
example is the response of bacteria to surplus or deprivation of
oxygen and/or different concentrations of antibiotics in an envi-
ronment, where an oxygen deficit environment could reduce the
replication rate of bacteria.
Following clear definitions of the variables characterizing the
state (and trait) of the individuals in the system as well as the
activities that could be carried out, we begin to ask questions of
the form “where, when, and how?” For instance, in modeling
Use of Individual-Based Mathematical Modelling to Understand More About. . . 101
2.4 Declaration of In step 3, we mentioned that various activities could be carried out
Assumptions About by each individual element in a system. However, we are yet to
Behavioral Rules of discuss how individual characteristics or actions could change in an
Individuals environment over time and how the behavior of an individual
might change as a consequence of the behavior exhibited by other
individual(s). It is important to have a clear thought on the possi-
bility that some of the individuals in our system might exhibit
adaptive behavior which could influence their form or actions,
and their interaction with other individuals and the environment
in which they reside. A good example is a probability that some
drug-susceptible bacteria in a particular environment can develop
mechanisms against antibiotics, thereby improving their fitness,
and bringing about resistant bacteria. In this case, how the drug
affects susceptible bacteria will be different from its effect on
mutants. Hence, we need to declare how the adaptability of indivi-
duals to their environment might impact the changes to an indivi-
dual’s behavior as well as the overall outcome of the research
problem. The integration of this behavior into the IBM con-
structed will not pose a problem, as one of the classical properties
of IBMs is the ability to integrate dynamics and feedback within the
model [30].
Further, the interactions between individuals within a spatial
environment might occur randomly. However, we are yet to incor-
porate the possibility that some individuals in an environment
could detect changes in their environment, which could ultimately
influence their decisions. An example is the migration of immune
cells through space toward the site of infection invoked by a signal-
ing protein such as chemokine molecules. Also, a clear definition of
the outcome of the interaction between individuals of different
types should be made. Thus, it is important to state some underly-
ing assumptions of behavioral rules and attributes that could
emerge in our model. This could be stochasticity, the ability to
detect or sense signals, emergence (i.e., the behavior/outcome that
emerges only after the interaction with others), and other attributes
we consider to be significant for our model. Again, we present some
guidelines to the questions that could guide our model construc-
tion when declaring attributes that might occur over time in our
model.
(p) Do we have individual(s) in our model that are able to adapt to
the environment?
(q) How does adaptive behavior influence the actions, and fitness
of individuals, as well as their interaction with other indivi-
duals and the environment?
(r) What level of randomness (stochasticity) is to be included in
the behavioral rules? For example, how much variation is there
in the replication rate of particular bacteria?
Use of Individual-Based Mathematical Modelling to Understand More About. . . 103
(s) Do we have individual(s) that can detect signals, and how does
this attribute contribute to the model to be constructed?
(t) Are there any assumed interactions among the individuals in
our system, and what should the impact of these
interactions be? E.g., should we assume that the interaction
between an active T cell and bacteria will always result in the
death of the bacteria?
In Algorithm 1, we present an example of a step-wise represen-
tation that could be adapted when modeling the evolution of
antibiotic resistance using the steps described above. This includes
some possible activities/actions of the individuals, the order and
situation in which they could be carried out, the subsequent effect
on their state (and trait) updates, as well as how a behavior change
could occur. For clarity, we have emphasized the state (and traits)
updates in red and detailed the processes/actions carried out as
comments. Here, we have included an option for more than one
phenotype of bacteria. Where a single phenotype is modeled, the
process of state change of bacteria should be ignored. Although it is
not explicitly represented, it is assumed that regulated cytokines
production is involved in the recruitment of both innate and adap-
tive immune cells. Figure 3b gives a visualization to Algorithm
1 (found in Fig. 3a) in the form of a flowchart. It is important to
note that the immune system consists of a complex network of cells,
with many processes involved in their ability to perform various
functions that aid in the containment of bacterial infection. Thus,
the flowchart presented should be taken as an example with
assumptions/simplifications to the response of the immune cells
and the treatment of bacterial infection using antibiotic(s).
2.5 Defining the For IBM simulations, one must think of the initial values to assign
Initialization of State to some (or all) of the state and (traits) variables of the individuals in
(and Trait Variables) the model when starting the simulation. In within-host modeling
of antibiotic resistance; this could be an initial deposit of a user-
defined number of bacteria and/or tissue-resident immune cells
having the state (and traits) defined in step 2. To ensure our model
is fully representative of possibilities that could arise, simulations
should be carried out many times. So, we must also decide if we
intend to vary the initial values or use the same values for all
simulations.
2.6 Defining The use of appropriate parameter values during the simulation of
Parameter Input the IBM is as significant as the model design and construction. This
is because the overall outcome of the bacterial infection under
investigation will be in response to the parameter input of the
model constructed. Hence, the use of accurate parameter values
ensures we simulate the model with biologically relevant outcomes.
The choice of parameter inputs is often driven by experimental
Algorithm 1: An example of a model design with particular assumptions and simplifications.
end if
end while
recruitment of resting macrophages (Mr ) through blood vessels l> At a specific rate
end if
Fig. 3 (a): An example of an algorithm that could be adapted when modeling the evolution of antibiotic
resistance. (b): An example of a model design with particular assumptions and simplifications made. The
dashed lines with bullet points indicate activities/actions being carried out. The solid lines with arrows show
the outcome following the activity and the result of the interactions between individuals after some specific
actions
Use of Individual-Based Mathematical Modelling to Understand More About. . . 105
end if
end if
end while
time
end if
end while
end while
end while
end if
Fig. 3 (continued)
106 Aminat Yetunde Saula et al.
Yes No
Is oxygen concentration
surplus?
Yes
Update macrophages state Phagocytosis
Is chemotactic signal = 0?
Replication of bacteria from resting to active
persist (at user-defined
rates)? No
Yes Activation Update macrophages trait
Active macrophages move
from healthy to
chemotactically and infected/chronically infected
Active T cells move
interact with bacteria (kills
randomly and interact
bacteria at a specific rate)
with resting (healthy)
Recruitment macrophages
Antibiotic Concentration >= MBC Inhibit bacteria growth Update bacteria behaviour from
drug-susceptible to resistant
Kill drug-susceptible bacteria
Toxicity (at different rates for slow-and Update bacteria state to
fast-growing) dead
Fig. 3 (continued)
Use of Individual-Based Mathematical Modelling to Understand More About. . . 107
2.7 Decision In an IBM simulation, the length of the simulation may be long
on Result Display with short time steps, hence we need to decide the interval at which
the output should be captured. With this, it will be possible to see
how the system evolves (e.g., the emergence of resistant strains
from an original population of drug-susceptible bacteria), the activ-
ities carried out by individuals over time (i.e., the process outputs
such as replication, growth inhibition, death), and the emergence
that arise from the behavior of individuals as well as the environ-
mental attributes. In addition to displaying the snapshots of the
outcome of the simulation over space and time, due to the inherent
stochasticity in biological processes, to get an accurate representa-
tion of the system, the simulations should be run many times, after
which the mean values and the variances could be plotted to see the
summary outcomes.
Following these steps and providing suitable and well-thought
answers to the questions posed is a basis for the development of a
good IBM, which will make implementation and execution easier.
In addition, if these steps are followed appropriately, we expect that
the location of the individuals in our model, the individual and
environmental attributes, and the spatial relationships that occur
over time to be captured.
Although there is a wide variety of software environments for
IBMs such as Netlogo (http://ccl.northwestern.edu/netlogo),
MASON (https://cs.gmu.edu/~eclab/projects/mason/), Agents.jl
(https://juliapackages.com/p/agents), however, we would encour-
age the use of object-oriented programming (OOP) languages such
as Julia, C++, Python, and Java. There are numerous advantages of
OOP, including the facilitation of code maintenance and reusability,
easier troubleshooting and collaborative development, flexibility
through polymorphism, and many more.
References
1. Macal CM (2016) Everything you need to 4. Bowness R, Chaplain MA, Powathil GG, Gille-
know about agent-based modelling and simu- spie SH (2108) Modelling the effects of bacte-
lation. J Simul 10:144–156 rial cell state and spatial location on
2. Bozic I, Antal T, Ohtsuki H, Carter H, Kim D, tuberculosis treatment: insights from a hybrid
Chen S et al (2010) Accumulation of driver and multiscale cellular automaton model. J Theor
passenger mutations during tumor progres- Biol 446:87–100
sion. Proc Natl Acad Sci 107:18545–18550 5. Strydom N, Gupta SV, Fox WS, Via LE,
3. Wilensky U, Rand W (2015) An introduction Bang H, Lee M et al (2019) Tuberculosis
to agent-based modeling: modeling natural, drugs’ distribution and emergence of resistance
social, and engineered complex systems with in patient’s lung lesions: a mechanistic model
NetLogo. MIT Press and tool for regimen and dose optimization.
PLoS Med 16:e1002773
108 Aminat Yetunde Saula et al.
6. Hult C, Mattila JT, Gideon HP, Linderman JJ, 19. Dancik GM, Jones DE, Dorman KS (2010) An
Kirschner DE (2021) Neutrophil dynamics agent-based model for leishmania major
affect mycobacterium tuberculosis granuloma infection. In: Unifying themes in complex sys-
outcomes and dissemination. Front Immunol tems. Springer, pp 243–250
12:712457 20. Whitlock AO, Juliano JJ, Mideo N (2021)
7. Joslyn LR, Linderman JJ, Kirschner DE (2022) Immune selection suppresses the emergence
A virtual host model of mycobacterium tuber- of drug resistance in malaria parasites but facil-
culosis infection identifies early immune events itates its spread. PLoS Comput Biol 17:
as predictive of infection outcomes. J Theor e1008577
Biol 539:111042 21. Rowlatt CF, Chaplain MA, Hughes DJ, Gille-
8. Feder AF, Pennings PS, Hermisson J, Petrov spie SH, Dockrell DH, Johannessen I et al
DA (2019) Evolutionary dynamics in (2022) Modelling the within-host spread of
structured populations under strong popula- sars-cov-2 infection, and the subsequent
tion genetic forces. G3: genes. Genom Genet immune response, using a hybrid, multiscale,
9:3395–3407 individual-based model. part i: Macrophages.
9. Monaco H, Liu KS, Sereno T, Deforet M, Tay- bioRxiv
lor BP, Chen Y et al (2022) Spatial-temporal 22. Kirschner DE, Hunt CA, Marino S, Fallahi-
dynamics of a microbial cooperative behavior Sichani M, Linderman JJ (2014) Tuneable res-
resistant to cheating. Nat Commun 13:1–11 olution as a systems biology approach for
10. Railsback SF, Grimm V (2019) Agent-based multi-scale, multi-compartment computa-
and individual-based modeling: a practical tional models. Wiley Interdiscip Rev Syst Biol
introduction. Princeton University Press Med 6:289–309
11. Anderson AR, Quaranta V (2008) Integrative 23. Radchuk V, Kramer-Schadt S, Berger U,
mathematical oncology. Nat Rev Cancer 8: Scherer C, Backmann P, Grimm V (2021) Indi-
227–234 vidual- based models. Demographic methods
12. West J, Robertson-Tessi M, Anderson AR across the tree of life. Oxford University Press,
(2022) Agent-based methods facilitate integra- p 213
tive science in cancer. Trends Cell Biol 33:300– 24. Breckling B (2002) Individual-based model-
311 ling potentials and limitations. Sci World J 2:
13. Segovia-Juarez JL, Ganguli S, Kirschner D 1044–1062
(2004) Identifying control mechanisms of 25. Marino S, Hogue IB, Ray CJ, Kirschner DE
granuloma formation during m. tuberculosis (2008) A methodology for performing global
infection using an agent-based model. J Theor uncertainty and sensitivity analysis in systems
Biol 231:357–376 biology. J Theor Biol 254:178–196
14. Wong EA, Joslyn L, Grant NL, Klein E, Lin 26. Renardy M, Hult C, Evans S, Linderman JJ,
PL, Kirschner DE et al (2018) Low levels of Kirschner DE (2019) Global sensitivity analysis
t cell exhaustion in tuberculous lung granulo- of biological multiscale models. Curr Opin
mas. Infect Immun 86:e00426–e00418 Biomed Eng 11:109–116
15. Joslyn LR, Renardy M, Weissman C, Grant 27. Hamis S, Stratiev S, Powathil GG (2021)
NL, Flynn JL, Butler JR et al (2021) Temporal Uncertainty and sensitivity analyses methods
and spatial analyses of tb granulomas to predict for agent- based mathematical models: an
long-term outcomes. In: Complex systems and introductory review. In: The physics of cancer:
computational biology approaches to acute research advances, pp 1–37
inflammation. Springer, pp 273–291 28. Lipworth S, Hammond RJH, Baron VO,
16. Joseph IM, Kirschner D (2004) A model for Hu Y, Coates A, Gillespie SH (2016) Defining
the study of helicobacter pylori interaction with dormancy in mycobacterial disease. Tuberculo-
human gastric acid secretion. J Theor Biol 228: sis 99:131–142
55–80 29. Pitcher MJ, Bowness R, Dobson S, Eftimie R,
17. Weathered C, Pennington K, Escalante P, Pie- Gillespie SH (2020) Modelling the effects of
naar E (2022) The role of biofilms, bacterial environmental heterogeneity within the lung
phenotypes, and innate immune response in on the tuberculosis life-cycle. J Theor Biol
mycobacterium avium colonization to infec- 506:110381
tion. J Theor Biol 534:110949 30. Auchincloss AH, Garcia LMT (2015) Brief
18. Zarrabi N, Mancini E, Tay J, Shahand S, Sloot introductory guide to agent-based modeling
PM (2010) Modeling hiv-1 intracellular repli- and an illustration from urban health research.
cation: two simulation approaches. Procedia Cad Saude Publica 31:65–78
Comput Sci 1:555–564
Chapter 11
Abstract
Tuberculosis (TB) is the most common cause of death from an infectious disease. Although treatment has
been available for more than 70 years, it still takes too long and many patients default risking relapse and the
emergence of resistance. It is known that lipid-rich, phenotypically antibiotic-tolerant, bacteria are more
resistant to antibiotics and may be responsible for relapse necessitating extended therapy. Using a micro-
fluidic system that acoustically traps live mycobacteria, M. smegmatis, a model organism for M. tuberculosis
we can perform optical analysis in the form of wavelength-modulated Raman spectroscopy (WMRS) on the
trapped organisms. This system can allow observations of the mycobacteria for up to 8 h. By adding
antibiotics, it is possible to study the effect of antibiotics in real-time by comparing the Raman fingerprints
in comparison to the unstressed condition. This microfluidic platform may be used to study any microor-
ganism and to dynamically monitor its response to many conditions including antibiotic stress, and changes
in the growth media. This opens the possibility of understanding better the stimuli that trigger the lipid-
rich downregulated and phenotypically antibiotic-resistant cell state.
Key words Raman spectroscopy, Tuberculosis, Dormancy, Single cell, Antibiotic resistance, Real-time
monitoring, Rapid diagnosis
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_11,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
109
110 Mingzhou Chen et al.
2 Materials
2.1 Bacteria, Culture 1. Middlebrook 7H9 broth and agar powder (Sigma-Aldrich UK)
Conditions, and was used to make sterile broth and plates.
Contamination Control 2. Middlebrook 7H11 (Sigma Aldrich, UK).
3. The medium should be supplemented with Tween 80 (0.05%
(v/v), Fischer Scientific) and with 2 mL of glycerol for 450 mL
of 7H9 broth.
4. M. smegmatis (NCTC 8159) National Type Culture Collection
(UK Health Security Agency) (see Note 1).
5. Benchtop centrifuge that accepts 20 mL containers and spins at
approximately (find 4500 g).
6. Phosphate Buffered saline.
7. Brain Heart Infusion Agar (Oxoid).
8. Isoniazid (Sigma-Aldrich).
9. Virkon solution.
10. Sterile 96-well plates (Nunc, DK).
11. Incubator (add details).
12. Spectrophotometer (Biochrom UK).
13. 96-well plate reader with the ability to detect at 600 nm absor-
bance (BMG Labtech).
Fig. 1 Outline of the overall arrangement of the acoustic trap and application of the Raman laser
2.3 Raman 1. We use a Raman system similar to that described in our previ-
Spectroscopy System ous work (refs) with the exception that the fluorescence imag-
ing function was not required. The various components of the
system are noted below.
2. A Nikon microscope (Nikon TE2000-E).
3. A tunable Ti:Sa laser (M2 SolsTis lasers, 1 W@785 nm).
4. A spectrometer system (Andor Shamrock SR303i + Newton
camera).
2.5 Precision Syringe 1. Medium is pumped into the chamber at a very slow but con-
Pump trolled speed into the microfluidic chamber by a precision
syringe pump (AL2000, World Precision Instruments), due
to the small volume of the chamber.
114 Mingzhou Chen et al.
3 Methods
3.3 Raman 1. This is the critical part of the experimental process when Raman
Spectroscopy of spectra are taken from the test organism under differing con-
Trapped Bacteria ditions over time. Careful consideration should be taken as to
how the conditions are to be changed. In the instance
described here, we can envisage following the changes in the
trapped culture over time when there is no stress. This is
important to ensure that the experimental setup is not
116 Mingzhou Chen et al.
Table 1
Raman peaks with their assignments and overall trends using three experiments for each condition
over time (up to 8 h of measurement)
Raman peak (cm-1) Associated chemical bonds Potential Raman peak assignments
635 C—C, C—S Tyrosine
783 Nucleic acids
1007 C—C Phenylalanine
1040 Carbohydrates
1080 C—C, C—O—H, C—N, C—O Lipids, carbohydrates, proteins
1130 C—N, C—C, C—O Carbohydrates, proteins, lipids
1150 C—C Carotenoids
1303 CH2 Lipids
1443 CH2, CH3 Lipids, proteins
1523 C〓C Carotenoids
1606 C〓C Tyrosine, phenylalanine
1658 C〓C, C〓O Lipids, proteins, amide I
1750 C〓O Lipids
Monitoring Live Mycobacteria in Real-Time Using a Microfluidic Acoustic. . . 117
4 Notes
References
1. Tortoli E, Cichero P, Piersimoni C et al (1999) microparticles in suspension and reaction mon-
Use of BACTEC MGIT 960 for recovery of itoring using Raman microspectroscopy. Anal
mycobacteria from clinical specimens: multi- Chem 79:7853–7857
center study. J Clin Microbiol 37:3578–3582 11. Tuckermann R, Puskar L, Zavabeti M et al
2. Zhang X, Jiang X, Yang Q et al (2018) Online (2009) Chemical analysis of acoustically levi-
monitoring of bacterial growth with an electri- tated drops by Raman spectroscopy. Anal Bioa-
cal sensor. Anal Chem 90:6006–6011 nal Chem 394:1433–1441
3. Godin M, Delgado FF, Son S et al (2010) 12. Hammarström B, Laurell T, Nilsson J (2012)
Using buoyant mass to measure the growth of Seed particle-enabled acoustic trapping of bac-
single cells. Nat Methods 7:387–390 teria and nanoparticles in continuous flow sys-
4. Gumbo T, Louie A, Deziel MR et al (2004) tems. Lab Chip 12:4296–4304
Selection of a moxifloxacin dose that sup- 13. Lee KS, Palatinszky M, Pereira FC et al (2019)
presses drug resistance in Mycobacterium An automated Raman-based platform for the
tuberculosis, by use of an in vitro pharmacody- sorting of live cells by functional properties.
namic infection model and mathematical mod- Nat Microbiol 4:1035–1048
eling. J Infect Dis 190(1642):1651 14. Gillespie SH, Crook AM, McHugh TD et al
5. Gumbo T, Dona CSWS, Meek C et al (2009) (2014) Four-month moxifloxacin-based regi-
Pharmacokinetics-pharmacodynamics of pyra- mens for drug-sensitive tuberculosis. New
zinamide in a novel in vitro model of tubercu- Engl J Med 371:1577–1587
losis for sterilizing effect: a paradigm for faster 15. Merle CS, Fielding K, Sow OB et al (2014) A
assessment of new antituberculosis drugs. Anti- four-month gatifloxacin-containing regimen
microb Agents Chemother 53(3197):3204 for treating tuberculosis. New Engl J Med
6. Gumbo T, Angulo-Barturen I, Ferrer-Bazaga S 371:1588–1598
(2015) Pharmacokinetic-pharmacodynamic 16. Jindani A, Harrison TS, Nunn AJ et al (2014)
and dose-response relationships of Antituber- High-dose rifapentine with moxifloxacin for
culosis drugs: recommendations and standards pulmonary tuberculosis. New Engl J Med
for industry and academia. J Infect Dis 211: 371:1599–1608
S96–S106 17. Mukamolova GV, Mukamolova GV, Turapov
7. Gutiérrez-Ramos S, Hoyos M, Ruiz-Suárez JC O et al (2010) Resuscitation-promoting factors
(2018) Induced clustering of Escherichia coli reveal an occult population of tubercle bacilli in
by acoustic fields. Sci Rep 8:4668 sputum. Am J Resp Crit Care 181(174):180
8. Bazou D, Kearney R, Mansergh F et al (2011) 18. Sloan DJ, Mwandumba HC, Garton NJ et al
Gene expression analysis of mouse embryonic (2015) Pharmacodynamic modeling of bacil-
stem cells following levitation in an ultrasound lary elimination rates and detection of bacterial
standing wave trap. Ultrasound Med Biol 37: lipid bodies in sputum to predict and under-
321–330 stand outcomes in treatment of pulmonary
9. Baron VO, Chen M, Clark SO et al (2017) tuberculosis. Clin Infect Dis 61:1–8
Label-free optical vibrational spectroscopy to 19. Hammond RJH, Baron VO, Oravcova K et al
detect the metabolic state of M. tuberculosis (2015) Phenotypic resistance in mycobacteria:
cells at the site of disease. Sci Rep 7:9844 is it because I am old or fat that I resist you? J
10. Ruedas-Rama MJ, Domı́nguez-Vidal A, Radel Antimicrob Chemother 70:2823–2827
S et al (2007) Ultrasonic trapping of
Monitoring Live Mycobacteria in Real-Time Using a Microfluidic Acoustic. . . 119
20. Phillips PPJ, Mendel CM, Burger DA et al mycobacteria with a microfluidic acoustic-
(2016) Limited role of culture conversion for Raman platform. Commun Biol 3:236
decision-making in individual patient care and 22. Billington OJ, McHugh TD, Gillespie SH
for advancing novel regimens to confirmatory (1999) Physiological cost of rifampin resistance
clinical trials. BMC Med 14:19 induced in vitro in mycobacterium tuberculo-
21. Baron VO, Chen M, Hammarstrom B et al sis. Antimicrob Agents Chemother 43:1866–
(2020) Real-time monitoring of live 1869
Chapter 12
Abstract
Going back in time through a phylogenetic tree makes it possible to evaluate ancestral genomes and assess
their potential to acquire key polymorphisms of interest over evolutionary time. Knowledge of this kind
may allow for the emergence of key traits to be predicted and pre-empted from currently circulating strains
in the future. Here, we present a novel genome-wide survival analysis and use the emergence of drug
resistance in Mycobacterium tuberculosis as an example to demonstrate the potential and utility of the
technique.
Key words Genome-wide survival analysis, Phylogenetics, Ancestral States, Hazard, Acquisition of
drug resistance, Evolutionary time, Traits, Phenotypes, Prediction
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_12,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
121
122 Arturo Torres Ortiz and Louis Grandjean
Fig. 1 Phylogenetic survival analysis workflow. The proposed approach starts with a genome alignment as
input, from which a phylogenetic tree is inferred. Different phenotypes are then determined for each isolate at
the tips of the tree. For representation, the phenotype of interest is characterized by red color. The ancestral
states for each phenotype are then inferred. To identify genetic signatures of genotypes leading to the
phenotype of interest, we select the last genotype without the phenotype (the pre-resistance node in the
example) and perform a time-to-event analysis where the time is represented by the branch length between
nodes and the event by the presence or absence of the phenotype of interest. Once the phenotype has evolved
within a branch, the rest of nodes are not counted for that clade
1. Genome alignment.
2. Phylogenetic inference.
3. Ancestral state and sequence reconstruction.
4. Survival analysis.
2 Materials
2.2 Software The following R packages are recommended to run the analysis:
(a) Ape [10].
(b) Phangorn [11].
(c) BactDating [12].
(d) Survival [13].
In addition, multiple programs are required to prepare the
data for analysis and for phylogenetic inference. Recom-
mended software is mentioned in later sections.
3 Methods
The proposed workflow takes raw reads as input, and it outputs the
association between each polymorphism to the phenotype of inter-
est. All the scripts are located in our GitHub page (github.com/
arturotorreso/mtb_pre-resistance). The recommended software
and the choice of programming language are optional.
3.1 Pre-processing Before preparing the alignment, the raw reads need to be
pre-processed as outlined in the pseudoseq_pipeline.sh using Trim-
momatic in order to remove illumina adaptors and low-quality
bases [14].
3.2 Genome As genomes are compared at each genomic position, all sequences
Alignment need to be in the same genomic coordinates. Therefore, all
pre-processed reads need to be aligned to the same reference
genome. In order to increase the accuracy of the final alignment,
we recommend preparing the genome alignment using the
pre-processed reads as shown in pseudoseq_pipeline.sh by using
the following software:
(a) Spades (https://github.com/ablab/spades) [15].
First, short-reads from each isolate are assembled into
longer contigs in order to increase the probability that a contig
is properly mapped to the reference genome.
(b) Minimap2 (https://github.com/lh3/minimap2) [16].
Contigs are then mapped against the reference genome
using Minimap2. The output is in the Sequence Alignment/
Map format (SAM/BAM)
(c) SAMtools (https://github.com/samtools/samtools) [17].
Phylogenetic Survival Analysis 125
3.3 Phylogenetic A phylogenetic tree can be inferred using a wide variety of methods.
Inference For this workflow, we propose using a time-calibrated phylogeny,
where the branch length is measured in calendar units. This can be
performed in one step using methods such as BEAST [23],
although this method is limited by the number of isolates in the
analysis and it can become unfeasible with a large sample size.
Alternatively, time calibration can be performed in two steps by
inferring a phylogenetic tree and then time-calibrating it using
methods such as those implemented in BactDating [12], TreeDater
[24], or LSD [25], which allow to obtain a time-calibrated phylog-
eny even in large studies. Instructions about how to run BactDat-
ing can be found at runBactDat.R. When a time-calibrated
phylogenetic tree cannot be accurately inferred, the branch length
of the tree in units of genetic distance can be used as the time
variable within the phylogenetic survival model.
126 Arturo Torres Ortiz and Louis Grandjean
3.5 Ancestral The ancestral sequences of internal nodes of the phylogenetic tree
Sequence can be inferred using similar methods as with the ancestral state
Reconstruction reconstruction previously described (see Note 1).
To perform ancestral sequence reconstruction using Phangorn,
see Note 2.
3.6 Phylogenetic Once the ancestral reconstructions have been finalized, the survival
Survival Analysis analysis can be performed as described in the R scripts within the
GitHub repository survival_analysis.R and survTree_functions.R.
In short, they are summarized as follows:
1. The phylogenetic tree is scanned from the root to the tips.
2. For each branch, the branch length and the phenotypes of the
internal and external nodes are recorded (e.g., 0 for absence
and 1 for presence, or continuous phenotype).
3. Once a phenotype has evolved in a specific branch, subsequent
branches with no phenotype change are not considered.
The time-to-event analysis is performed within R using the
Survival package. A time-to-event model can be fit using a combi-
nation of the Surv function, which creates a survival object, and the
survfit function. Moreover, a Cox Proportional-Hazards Model can
be fitted using the coxph function to compute the instantaneous
hazard of acquiring the phenotype of interest.
Using these functions, the model can be fitted to calculate
differences in survival rates between groups. In our analysis, we
computed differences in drug resistance acquisition between clades,
or between mono-resistant clades in comparison to fully susceptible
groups.
Additionally, a genome-wide association analysis is performed
by fitting a Cox model by selecting the internal node for each
branch where a phenotype change has occurred and the remaining
nodes where the phenotype is not present. This association can be
repeated for each polymorphic site of interest, as indicated in the
scripts gwas1_alignment.R and gwas2_analysis.R.
The Survival and Cox models can be substituted by a lin-
ear model, a logistic regression model or any other suitable
model to compute the association without using the branch length
or time within the model. This can be useful to calculate odds ratios
or to run the association when the branch lengths of the tree are
unreliable.
Phylogenetic Survival Analysis 127
3.7 Output The final output of the analysis consists of the fitted parameters of
the survival and the Cox model, including a p-value reflecting
statistical significance of the differences between groups, both in
the case of clade differences or genotypes.
4 Notes
Acknowledgments
References
1. Gygli SM, Borrell S, Trauner A, Gagneux S methods using simulated genomes and pheno-
(2017) Antimicrobial resistance in Mycobacte- types. Microb Genom 6:e000337. https://doi.
rium tuberculosis: mechanistic and evolution- org/10.1099/mgen.0.000337
ary perspectives. FEMS Microbiol Rev 41:354– 4. Collins C, Didelot X (2018) A phylogenetic
373. https://doi.org/10.1093/femsre/ method to perform genome-wide association
fux011 studies in microbes that accounts for popula-
2. Walker TM, Kohl TA, Omar SV et al (2015) tion structure and recombination. PLoS Com-
Whole-genome sequencing for prediction of put Biol 14:e1005958. https://doi.org/10.
Mycobacterium tuberculosis drug susceptibil- 1371/journal.pcbi.1005958
ity and resistance: a retrospective cohort study. 5. Gagneux S (2018) Ecology and evolution of
Lancet Infect Dis 15:1193–1202. https://doi. Mycobacterium tuberculosis. Nat Rev Micro-
org/10.1016/S1473-3099(15)00062-6 biol 16:202. https://doi.org/10.1038/
3. Saber MM, Shapiro BJ (2020) Benchmarking nrmicro.2018.8
bacterial genome-wide association study
128 Arturo Torres Ortiz and Louis Grandjean
Abstract
Antibiotic resistance is a global challenge likely to cost trillions of dollars in excess costs in the health system
and more importantly, millions of lives every year. A major driver of resistance is the absence of susceptibility
testing at the time a healthcare worker needs to prescribe an antimicrobial. The effect is that many
prescriptions are unintentionally wasted and expose mutable organisms to antibiotics increasing the risk
of resistance emerging. Often simplistic solutions are applied to this growing issue, such as a naı̈ve drive to
increase the speed of drug susceptibility testing. This puts a spotlight on a technological solution and there
is a multiplicity of such candidate DST tests in development. Yet, if we do not define the necessary
information and the speed at which it needs to be available in the clinical decision-making progress as
well as the necessary integration into clinical pathways, then little progress will be made. In this chapter, we
place the technological challenge in a clinical and systems context. Further, we will review the landscape of
some promising technologies that are emerging and attempt to place them in the clinic where they will have
to succeed.
Key words Drug susceptibility testing, Rapid diagnostics, Microbiology automation, Microbiology
systems, Limit of detection (LoD), Photonics, Genome sequencing, Polymerase chain reaction
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_13,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
129
130 Stephen H. Gillespie and Robert J. H. Hammond
when the specimen is produced and sent for culture, loss of viability
in transit, especially true for anaerobes, results in false negative
result and no material for DST.
Thus, in many circumstances, culture is not necessarily the
optimal way of identifying the susceptibility pattern of a causative
pathogen. Although bacteria grow rapidly, dividing in as little as
once every 20 min, the high limit of detection (approx. 108 cfu/mL
Vs. circulating pathogens at ~<100 cfu/mL) means that a positive
result takes many hours to signal positive. This limitation also
applies when these systems are used to perform drug susceptibility
tests. In most circumstances defining the species of the infecting
organism will provide a strong treatment indication to the clinician
based on prior knowledge. For example, detecting the presence of
Staphylococcus aureus suggests a limited number of antibiotics that
will be useful. Recent developments using MALDI-TOF on pri-
mary cultures have proved very useful in generating a report in a
single working shift [1–3]. This is illustrated in Table 1, which
summarizes the time to result for common rapid susceptibility
testing platforms that depend on culture. But, perhaps, the greatest
drawback for many of these systems is the need to start with a pure
culture of the causative organism in most circumstances. This
means that there is usually an initial culture stage imposing a
significant delay of between 8 and 24 h. Despite the recent focus
on antibiotic stewardship, the importance of rapid diagnostics is
often overlooked in the discussions [12].
Table 1
A summary of selected phenotypic drug susceptibility test methods
Time to
positive
Technique Methodology result References
Microscan Lyophilized fluorogenic substrates and pH indicators 16 h [4]
Phoenix Detects change in redox using a colorimetric oxidation- 20 h [5]
reduction indicator
Vitek 2 Identification of biochemical reactions and nutrient usage ~10.4 h [6]
(3.3–17.5 h)
BacT- A colorimetric sensor and reflected light monitor production 17 h [7]
ALERT of carbon dioxide (CO2) that is dissolved in the culture
medium changing the color of the medium
BACTEC Increases in CO2 production due to microbial growth are 12.5 h [7]
detected by non-invasive fluorescence
VersaTREK Detects pressure changes due to gas consumption (O2) ~23.7 h [8–11]
and/or production (CO2, N2, H2, etc.) by microorganisms (11–36.4 h)
in a BC bottle.
Rapid Drug Susceptibility Testing to Preserve Antibiotics 131
Question 1
Question 2
Question 3
failure
Deliver
effective
Identify treatment
susceptibilities
Fig. 1 An overview of the interaction between clinical questions and the microbiological diagnostic process
There are different broad methods to address this issue: detect the
growth of bacteria or detect a drug susceptibility characteristic
directly from the specimen. Despite the difficulties of defining
what is meant by “fast” there is almost universal agreement of the
need to speed up the current laboratory methodologies [15]. These
can be generally divided into three broad categories: detecting the
metabolic effects of growth, detecting growth directly, or detecting
Rapid Drug Susceptibility Testing to Preserve Antibiotics 133
4.1 The BD Phoenix The Phoenix platform determines growth by detecting both the
Automated turbidity of growth and a fluorescent oxidation-reduction indica-
Microbiology System tor. The MIC is determined by using doubling dilutions of
pre-prepared antimicrobials in manufacturer-prepared cards.
There is a degree of flexibility in the content of the cards available
that can be adapted to the antibiotic stewardship program. CLSI-
based breakpoints are employed and interpreted by an expert sys-
tem [18]. The onboard system can be customized to take account
of local reporting algorithms. It provides susceptibility testing
across a wide range of organisms encountered in routine microbio-
logical practice with a simple workflow [19–21].
It could be argued that the reason DST takes so long is the low
sensitivity of the methods to detect and quantify the number of
bacteria present in a sample. Here we summarize some of the
innovative approaches to improve the detection and quantification
of bacteria that are being applied to testing susceptibility. Light-
based solutions (or photonics) have been seen as a way of increasing
the sensitivity of detection of bacterial growth and there is a multi-
plicity of such methods (Table 2). These can detect metabolic
changes, detect reduced or inhibited growth, detect morphological
changes, detect reduced motility, or identify drug-induced death.
5.1 Detecting Detecting metabolic changes has been used to determine suscepti-
Metabolic Components bility by monitoring metabolic changes as bacteria are exposed to
antibiotics. Raman spectroscopy in its many forms has the capacity
to interrogate the chemical components of bacteria [34]. A MIC
methodology based on this approach has been described [35]
although time to result was at least 4 h. Trapping cultures in a
microfluidic acoustic chamber and examining the culture with
wavelength-modified Raman spectroscopy could be used to moni-
tor metabolic changes in real time detecting changes at a molecular
level. This is a promising technique reviewed in this volume but
could be cumbersome to deploy in its current configuration. Its
ability to interrogate the chemistry of an organism means that it
would be able to address problematic organisms such as mycobac-
teria and could be an important methodology for further research
at this time [31].
Rapid Drug Susceptibility Testing to Preserve Antibiotics 135
Table 2
Summary of selected photonics drug susceptibility test methodologies
Time to
Technique positive
name Principle result Target References
Bacterioscan Forward-angle laser light scattering 3h Urinary tract [28, 29]
infection
screening
SLIC Light-scattering and an integrating 10 min UTI, rapid [30]
chamber DST
Raman Monitoring changes in Raman spectrum tbd Research [31]
spectroscopy in bacteria exposed to antibiotic technique
oCelloscope Digital time-lapse angled field <3 h DST [32, 33]
microscopy to detect changes in
number or morphology
Accelerate- Morphokinetic analysis of bacterial 2–7 ha [8–10]
pheno growth and integrated interpretation
a
The lower figure represents time identification and the latter time to AST.
5.2 Detecting Growth Forward angle light scattering is used to detect change in urine
Through Bacterial samples as a means of detecting bacteriuria. The application pro-
Scattering vides a screening system that prevents unnecessary cultures and,
potentially unneeded antibiotic prescription in children and adults
[28, 29]. This being brought to market by Astrego focused on
detection of urinary tract infection.
Another approach the methodology, “Scattered light integrat-
ing collector” (SLIC) has been used to significantly decrease the
limit of bacterial detection to as few as 25 cfu mL-1 [30]. This is
achieved by using an integrating chamber, which captures all the
light of an incident laser beam. As the number of organisms
increases through growth, the amount of scattering increases
making a highly sensitive readout of the number of particles present
in a fluid. Using this technique identifies another important
finding: it is easier to determine whether a bacterium is resistant
than it is to confirm that it is sensitive. When the strain is resistant,
the growth curve follows that of the control, whereas in the case of
susceptibility it is necessary to run the culture until it is certain that
no growth is going to take place.
5.4 Camera Systems Using a complex microfluidic and time-lapse camera system mea-
suring bacterial length in real time, Baltekin et al. detected the
action of antibiotics on E. coli in approximately 15 min [36] .
Digital time-lapse angled field microscopy has been applied to
bacteria growing multi-well plates in the presence and absence of
antibiotic. Changes in size or morphology are plotted in real time
with 95% of results available within 3 h [32, 33]. Choi et al. used
bacterial imaging and microfluidic systems to measure the area
occupied by immobilized microbes. Differing morphological pat-
terns could be used to provide susceptibility results in 3–4 h [37]. A
commercial integrated system directed toward speeding up blood
culture diagnosis detects bacteria on a cartridge and records cell
division, growth patterns, and other changes are analyzed in real
time. This is given the name morphokinetic analysis and delivers
results for identification in 2 h and AST in 7 h [8, 9, 38, 39].
5.5 Atomic Force Atomic force microscopy (AFM) is a powerful technique to image
Microscopy many types of surface including bacteria and fungi. It is capable of
measuring adhesion strength, magnetic forces, and mechanical
properties using a tip 10–20 nm in diameter attached to a cantile-
ver. When the tip moves in response to interaction with the surface
it is measured by a laser. Atomic force microscopy generates a 3d
profile. Samples do not require specialized treatment that would
change the structure of the sample making drug susceptibility
testing possible. Atomic force microscopy has been used in a
research application to monitor the change in the cell wall of
bacteria and fungi while being treated with antibiotics opening
the possibility of using this technique to define susceptibility
[40, 41]. Its application is likely to remain most used in researching
the underlying mechanisms of resistance and response to antibiotics
rather than a diagnostic test [42].
References
negative bacteria isolated from blood cultures. 21. Junkins AD, Arbefeville SS, Howard WJ et al
Microbiol Spectr 9:e01836–e01821 (2010) Comparison of BD phoenix AP work-
10. Chapot V, Effenberg L, Dohmen-Ruetten J flow with Vitek 2. J Clin Microbiol 48:1929–
et al (2021) Evaluation of the accelerate 1931
pheno system for rapid identification and anti- 22. McGregor A, Schio F, Beaton S et al (1995)
microbial susceptibility testing of positive The MicroScan WalkAway diagnostic microbi-
blood culture bottles inoculated with primary ology system – an evaluation. Pathology 27:
sterile specimens from patients with suspected 172–176
severe infections. J Clin Microbiol 59(5): 23. York MK, Brooks GF, Fiss EH (1992) Evalua-
e02637-20 tion of the autoSCAN-W/A rapid system for
11. Samuel LP, Pimentel JD, Tibbetts RJ et al identification and susceptibility testing of
(2011) Comparison of time to positivity of gram-negative fermentative bacilli. J Clin
the VersaTREK® REDOX 80-mL and the Microbiol 30:2903–2910
REDOX EZ draw 40-mL blood culture bottles 24. Bondi A, Curtoni A, Peradotto M et al (2023)
for common bacterial bloodstream pathogens. Performance evaluation of BD Phoenix and
Diagn Microbiol Infect Dis 71:101–105 MicroScan WalkAway plus for determination
12. Hays JP, Mitsakakis K, Luz S et al (2019) The of fosfomycin susceptibility in enterobacterales.
successful uptake and sustainability of rapid Antibiotics 12:1106
infectious disease and antimicrobial resistance 25. Stefaniuk E, Baraniak A, Gniadkowski M et al
point-of-care testing requires a complex ‘mix- (2003) Evaluation of the BD Phoenix auto-
and-match’ implementation package. Eur J mated identification and susceptibility testing
Clin Microbiol Infect Dis 38:1015–1022 system in clinical microbiology laboratory prac-
13. van Belkum A, Bachmann TT, Lüdke G et al tice. Eur J Clin Microbiol Infect Dis 22:479–
(2019) Developmental roadmap for antimicro- 485
bial susceptibility testing systems. Nat Rev 26. Barry J, Brown A, Ensor V et al (2003) Com-
Microbiol 17:51–62 parative evaluation of the VITEK 2 advanced
14. Holmes AH, Moore LSP, Sundsfjord A et al expert system (AES) in five UK hospitals. J
(2016) Understanding the mechanisms and Antimicrob Chemother 51:1191–1202
drivers of antimicrobial resistance. Lancet 27. Khan A, Arias CA, Abbott A et al (2021) Eval-
387:176–187 uation of the Vitek 2, Phoenix, and MicroScan
15. Smith KP, Kirby JE (2019) Rapid susceptibility for antimicrobial susceptibility testing of steno-
testing methods. Clin Lab Med 39:333–344 trophomonas maltophilia. J Clin Microbiol 59:
16. Testing EC on AS European guidelines for e00654–e00621
breakpoint tables for interpretation of MICs 28. Davaro E, Tomaras AP, Chamberland RR et al
and zone diameters Version 11.0, valid from (2020) Evaluation of a novel light scattering
2021-01-01, e: https://www.eucast.org/ methodology for the detection of pathogenic
fileadmin/src/media/PDFs/EUCAST_files/ bacteria in urine. J Appl Lab Med 5:370–376
Breakpoint_tables/v_11.0_Breakpoint_ 29. Hassan F, Bushnell H, Taggart C et al (2019)
Tables.pdf Evaluation of BacterioScan 216Dx in compari-
17. Clinical and Laboratory Standards Institute son to urinalysis as a screening tool for diagno-
(2022) Performance standards for antimicro- sis of urinary tract infections in children. J Clin
bial susceptibility testing, Clinical and Labora- Microbiol 57:e00571-19
tory Standards Institute Wayne 30. Hammond RJH, Falconer K, Powell T et al
18. Winstanley T, Courvalin P (2011) Expert sys- (2022) A simple label-free method reveals bac-
tems in clinical microbiology. Clin Microbiol terial growth dynamics and antibiotic action in
Rev 24:515–556 real-time. Sci Rep-UK 12:19393
19. Carroll KC, Glanz BD, Borek AP et al (2006) 31. Baron VO, Chen M, Hammarstrom B et al
Evaluation of the BD phoenix automated (2020) Real-time monitoring of live mycobac-
microbiology system for identification and teria with a microfluidic acoustic-Raman plat-
antimicrobial susceptibility testing of entero- form. Commun Biol 3:236
bacteriaceae. J Clin Microbiol 44:3506–3509 32. Fredborg M, Rosenvinge FS, Spillum E et al
20. Carroll KC, Borek AP, Burger C et al (2006) (2015) Rapid antimicrobial susceptibility test-
Evaluation of the BD phoenix automated ing of clinical isolates by digital time-lapse
microbiology system for identification and microscopy. Eur J Clin Microbiol Infect Dis
antimicrobial susceptibility testing of staphylo- 34:2385–2394
cocci and enterococci. J Clin Microbiol 44: 33. Fredborg M, Rosenvinge FS, Spillum E et al
2072–2077 (2015) Automated image analysis for
142 Stephen H. Gillespie and Robert J. H. Hammond
susceptibility testing? Lancet Infect Dis 22: mycobacterium tuberculosis from genome
e121–e127 sequence data with Mykrobe. Wellcome Open
58. Cao Y, Parmar H, Gaur RL et al (2021) Xpert Res 4:191
MTB/XDR: a 10-color reflex assay suitable for 65. Sabiiti W, Azam K, Farmer ECW et al (2020)
point-of-care settings to detect isoniazid, fluo- Tuberculosis bacillary load, an early marker of
roquinolone, and second-line-injectable-drug disease severity: the utility of tuberculosis
resistance directly from mycobacterium molecular bacterial load assay. Thorax 75:
tuberculosis-positive sputum. J Clin Microbiol 606–608
59(3):e02314–20 66. Gillespie SH (2002) Evolution of drug resis-
59. Conradie F, Diacon AH, Ngubane N et al tance in mycobacterium tuberculosis: clinical
(2020) Treatment of highly drug-resistant pul- and molecular perspective. Antimicrob Agents
monary tuberculosis. N Engl J Med 382:893– Chemother 46:267–274
902 67. Gygli SM, Keller PM, Ballif M et al (2019)
60. Parcell BJ, Gillespie SH, Pettigrew KA et al Whole genome sequencing in mycobacterium
(2021) Clinical perspectives in integrating tuberculosis. Ann Transl Med 0:e02175–
whole genome sequencing into the investiga- e02118
tion of healthcare and public health 68. Consortium TCr (2022) A data compendium
outbreaks—hype or help? J Hosp Infect J associating the genomes of 12,289 mycobacte-
109:1–9. https://doi.org/10.1016/j.jhin. rium tuberculosis isolates with quantitative
2020.1 resistance phenotypes to 13 antibiotics. PLoS
61. Parcell BJ, Oravcova K, Pinheiro M et al Biol 20:e3001721
(2017) Pseudomonas aeruginosa intensive 69. Project TCrC and the 100000 G (2018) Pre-
care unit outbreak: winnowing of transmissions diction of susceptibility to first-line tuberculo-
with molecular and genomic typing. J Hosp sis drugs by DNA sequencing. N Engl J Med
Infect 98:282–288 379(1403):1415
62. Köser CU, Ellington MJ, Cartwright EJP et al 70. Jouet A, Gaudin C, Badalato N et al (2021)
(2012) Routine use of microbial whole Deep amplicon sequencing for culture-free
genome sequencing in diagnostic and public prediction of susceptibility or resistance to
health microbiology. PLoS Pathog 8: 13 anti-tuberculous drugs. Eur Respir J 57:
e1002824 2002338
63. Bradley P, Gordon NC, Walker TM et al 71. Cabibbe AM, Spitaleri A, Battaglia S et al
(2015) Rapid antibiotic-resistance predictions (2020) Application of targeted next-
from genome sequence data for Staphylococ- generation sequencing assay on a portable
cus aureus and mycobacterium tuberculosis. sequencing platform for culture-free detection
Nat Commun 6:10063 of drug-resistant tuberculosis from clinical
64. Hunt M, Bradley P, Lapierre SG et al (2019) samples. J Clin Microbiol 58(10):e00632–20
Antibiotic resistance prediction for
Chapter 14
Abstract
Mycobacterium tuberculosis is an infectious pathogen that requires biosafety level-3 laboratory for handling.
The risk of transmission is high to laboratory staff, and to manage the organism safely, it is necessary to
construct high containment laboratory facilities at great expense. This limits the application of tuberculosis
diagnostics to areas where there is insufficient capital to invest in laboratory infrastructure. In this method,
we describe a process of inactivating sputum samples by either heat or guanidine thiocyanate (GTC) that
renders them safe without affecting the quantification of viable bacteria. This method eliminates the need
for level 3 containment laboratory for the tuberculosis molecular bacterial load assay (TB-MBLA) and is
applicable in low- and middle-income countries.
Key words Heat inactivation, GTC, Containment level III, TB-MBLA, Tuberculosis, Diagnosis,
Accessibility
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_14,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
145
146 Bariki Mtafya et al.
2 Materials
3 Methods
3.2 RNA Extraction 1. Vortex the samples and transfer the whole liquid into lysing
matrix.
2. Transfer the tubes to a homogenizer (Precellys 24—PeQlab or
FastPrep—MP Biomedicals) (see Note 5).
3. Set the program 6.0 for 40 s on FastPrep or use 6000 rpm for
40 s, if using the Precellys 24 machine for homogenization of
the sample.
4. After homogenization, centrifuge the tubes for 5 min at
12,000 g at room temperature.
5. Leave the tube to stand for 5 min following centrifugation.
6. Label fresh 1.5 mL screw cap tubes (RNase free) with the
corresponding numbers and add 300 μL chloroform into
each clean tube.
7. From the homogenization tubes, carefully transfer the whole
liquid part to the chloroform-containing tube using a fine tip
pipette (be careful not to transfer bits of sample debris or lysing
matrix).
8. Vortex the tubes for 10 s each and incubate at room tempera-
ture for 5 min.
9. Centrifuge at 12000 g at room temperature for 5 min.
10. Label new 1.5 mL micro-centrifuge tubes with the
corresponding numbers.
11. Carefully transfer the upper aqueous phase to the fresh tubes
using 200 μL filter tips, being careful not to transfer any of the
interphase or lower layer.
12. Add 500 μL of 100% ice-cold ethanol to each tube.
13. Mix the contents by inverting the tubes 5 times.
14. Transfer the tubes containing samples to the 20 C freezer
overnight (or to the 80 C freezer for 15 min, if you are
doing the whole extraction on the same day).
15. Chill the microfuge to at least 12 C prior to centrifugation (set
it at 4 C).
16. Centrifuge the samples at 13,000 g for 20 min.
17. Discard the supernatant using a fine pipette tip.
18. Add 500 μL 70% ice-cold ethanol to each tube.
19. Centrifuge as above for a further 10 min.
20. Discard the supernatant using a Pasteur pipette with fine tip.
Use a new pipette for every sample.
21. Dry the RNA at 50 C in a heat hot block (approximately
20–30 min).
Quantifying Viable M. tuberculosis Safely Obviating the Need. . . 149
3.3 DNA Removal 1. This is done using Turbo DNase (Turbo DNA-free kit,
Ambion AM1907) or RNase-Free DNase Set (Qiagen).
2. Make a master mix of the Turbo DNase I; 10 buffer (blue
top) and DNase I enzyme (white top) for the number of
samples plus 2 extra. Per sample, add 10 μL buffer and 2 μL
DNase enzyme to each sample. For example, a 10-sample
master mix would contain 120 μL of 10 buffer and 2 μL of
DNase enzymes.
3. Mix by vortexing and then pipette 12 μL into each tube con-
taining your RNA extract.
4. Spin briefly (5–10 s at 13000 g) to bring everything to the
bottom of the tube. This step eliminates droplets hanging on
the tube wall and ensures that all RNA extract is in contact with
DNase enzyme.
5. Incubate at 37 C for 30 min in the hot block or incubator.
6. After 30 min, add an additional 1 μL of DNase enzyme to each
tube. Ensure all enzyme goes into the sample by pipetting up
and down (3) and then use the tip to stir around and mix the
enzyme with the rest of the mixture.
7. Incubate at 37 C for a further 30 min in the hot-block or
incubator.
8. Thaw the DNase inactivation reagent (white milky substance)
at least 10 min prior the finish of DNase incubation and keep in
the fridge. Resuspend by vortexing.
9. Add 10 μL of DNase inactivation reagent into each RNA
extract and incubate for 5 min at room temperature.
10. Vortex three times during the 5-min incubation step at room
temperature.
11. Centrifuge at 13,000 g for 2 min.
12. Label 1.5 mL screw cap tubes with sample ID and
extraction date.
13. Set your pipette at 105 μL and carefully pipette off the RNA to
a fresh 1.5 mL RNase-free tube without touching any of the
inactivation matrix.
14. Store the RNA at 4 C until use if running qPCR on the same
day or next day.
15. For long-term use, store the RNA at 80 C.
150 Bariki Mtafya et al.
3.4 Preparation of 1. A ready-to-go master mix containing the probes, primers, and
Master Mix Quantitec master mix is supplied in TUBE 1 and RT enzyme is
supplied in TUBE 2. Mix the content from TUBE 1 and
2 following the instructions included in the pack insert to
obtain the RT+ mix or follow the manufacturer’s instructions
if using a different kit (e.g., Life Arc kits).
2. The RT- master mix is supplied as “ready to use”, no additional
manipulation is required.
3. Add 16 μL of the master mix solution to each PCR tube
(if using tubes for the RotorGene Q) or in the 96/384 well
plate other PCR platforms, e.g., Quantistudio, CFX96, or ABI
7500 (see Note 6).
4. Add 4 μL of the RNA sample to each well containing 16 μL of
the master mix (total reaction volume will be 20 μL) as indi-
cated in Table 1.
5. Each sample and controls should be run in duplicate for RT+
reactions (containing the RT enzyme) and in a single reaction
for the master mix without reverse transcriptase (RT-).
6. Seal the tube with PCR tube capes and the plate with the
protective film, respectively.
7. Centrifuge the plate at 1500 rpm for 1 min to ensure all liquids
are at bottom of wells and transfer the plates into the CFX96,
ABI 750, or Quantistudio instrument (see Note 7).
8. No need to centrifuge the tubes if using the Rotorgene. The
machine will spin all the tubes when running.
Table 1
Quantities of the master mixes required and samples
4 Notes
References
1. Lange C, Aarnoutse R, Chesov D, van 5. Joshi R, Reingold AL, Menzies D, Pai M
Crevel R, Gillespie SH, Grobbel HP et al (2006) Tuberculosis among health-care work-
(2020) Perspective for precision medicine for ers in low- and middle-income countries: a
tuberculosis. Front Immunol 11:566608 systematic review. PLoS Med 3(12):
2. Siddiqi SH, Rüsch-Gerdes S (2006) MGIT 2376–2391
Procedure Manual: For BACTEC™ MGIT 6. World Health Organization (2012) Tubercu-
960™ TB System. 2006;(July):1–52 losis laboratory biosafety manual. World
3. World Health Organization (2021) Technical Health Organization Publication. ISBN
Report on critical concentrations for drug sus- 978 92 41504638, Geneva
ceptibility testing of isoniazid and the rifamy- 7. Allen V, Nicol MP, Tow LA (2016) Sputum
cins (rifampicin, rifabutin and rifapentine), processing prior to mycobacterium tuberculo-
1–86 p sis detection by culture or nucleic acid amplifi-
4. World Health Organisation (2022) Global cation testing: a narrative review. Res Rev J
tuberculosis report. Geneva Microbiol Biotechnol 5(1):96–108
152 Bariki Mtafya et al.
8. Mtafya B, Qwaray P, John J, Sichone E, bacillary load, an early marker of disease sever-
Shoo A, Gillespie SH et al (2023) A practical ity: the utility of tuberculosis molecular bacte-
approach to render tuberculosis samples safe rial load assay. Thorax 75:606–608
for application of tuberculosis molecular bacte- 10. Ntinginya NE, Bakuli A, Mapamba D,
rial load assay in clinical settings without a bio- Sabiiti W, Kibiki G, Minja LT et al (2022)
safety level 3 laboratory. Tuberculosis 138 Tuberculosis molecular bacterial load assay
(November 2022):102275 reveals early delayed bacterial killing in patients
9. Sabiiti W, Azam K, Farmer ECW, Kuchaka D, with relapse. Clin Infect Dis 76(3):990–994
Mtafya B, Bowness R et al (2020) Tuberculosis
Chapter 15
Abstract
The diagnosis and monitoring of tuberculosis treatment is difficult as many patients are unable to produce
sputum. This means that many patients are treated on the basis of clinical findings and consequently some
will be exposed to anti-tuberculosis drugs unnecessarily. Moreover, for those appropriately on treatment
and unable to produce a sputum sample, it will be impossible to monitor the response to treatment. We
have shown that stool is a potential alternative sample type for diagnosis of tuberculosis. Currently, available
protocols like the Xpert MTB/RIF use DNA as a target to detect Mycobacterium tuberculosis in stool but
DNA survives long after the organism is dead so it is not certain whether a positive test is from an old or a
partially treated infection. The TB MBLA only detects live organisms and thus, can be used to follow the
response to treatment. In this chapter, we describe a protocol for TB-MBLA, an RNA-based assay, and
apply it to quantify TB bacteria in stool.
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_15,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
153
154 Emmanuel Musisi et al.
2 Materials
2.1 Materials for Stool samples can be collected using a flushable Hystool® stool
Stool Collection collection device (https:/www.hystool.co.uk) (see Note 1).
3 Methods
3.1 Collecting the 1. Attaching the stool collection bag on the toilet.
Stool Sample Peel the backing paper off the adhesive pads on the corners
of the stool collection bag and attach it to the toilet seat using
the adhesive pads.
2. Utilizing the toilet.
Fold the piece of the toilet paper and place it at the bottom
of the stool bag to prevent direct contact between the bag and
stool. Loose stool could burst the bag. Instruct the patient to
sit on the toilet seat and use it normally. Care should be taken
for stool not to pass through the bag. Try as much as possible
not to pee in the stool because this could contaminate the
sample or dissolve it.
156 Emmanuel Musisi et al.
3. Pick a portion of the stool (about two scoops) into the sample
container using a scoop which is attached to the cap of the stool
container.
4. Unpeel the sticky pads from the toilet seat allowing the gab,
sticker, and contents to fall into the toilet and flush normally.
The stool bag used will dissolve in water and the sticker paper
was biodegradable.
3.2 Processing Stool 1. Place about ~6.0 g sterile stool into a pre-labeled stool con-
Samples tainer and homogenize in 10 mL phosphate buffered saline
(PBS). Incubate at room temperature for 15 min.
2. Centrifuge twice the resulting suspension at 3000 g for
20 and 10 min, respectively. Re-suspend the pellet in 6 mL of
phosphate-buffered saline and aliquot portions of 1 mL.
3.3 Preparation of 1. Primers and probes are lyophilized and stable at room temper-
Probes and Primers ature as set out in Table 1.
Notes for the TB-MBLA 2. Make stock primer per the manufacturer’s instructions.
Study
3. Make stock probe per the manufacturer’s instructions.
3.4 RNA Extraction 1. Thaw samples and extraction control (EC) and place on ice
until used.
2. Label same number of homogenization tubes with glass beads.
Label both lids and tubes.
3. Add 100 μL of EC directly into the sample, mix by inverting
20 times. Change tips between samples.
4. Spin at 3000 g for 30 min, carefully tip off GTC into 15 mL
centrifuge tubes.
5. Use a fine pipette tip to remove all the GTC. Avoid the pellet as
much as possible.
Table 1
Preparation of primer mix for the internal control (total volume = 200 μL)
3.5 DNAse Treatment 1. Prepare the master mix by adding 10 μL (10X) DNase buffer,
and 2 μL DNAse enzyme to a sterile microcentrifuge tube (see
Note 6).
2. Add 12 μL of DNAse master mix to each sample.
3. Vortex for 5 s.
4. Spin at 13,000 g for 10 s.
5. Incubate at 37 C for 30 min in the heat hot block.
6. Add 1 μL of DNAse enzyme to each tube. Pipette up and down
three times.
7. Incubate at 37 C for 30 min in the heat hot block.
8. Thaw DNase inactivation buffer (10 min before).
9. Add 10 μL of DNAse inactivation reagent.
158 Emmanuel Musisi et al.
Table 2
Quantities required for the master mix
Volume per RT+ Volume per RT Total vol. of RT+ Total vol. of RT
Master mix reaction μL reaction μL reactions μLa reactions μLa
Quantitec mix 10 10
(PCR mix)
Mtb 16S primer 0.4 0.4
mix (F + R)
IC primer mix 0.4 0.4
(F + R)
Mtb 16S-FAM 0.2 0.2
probe
IC probe 0.2 0.2
RT enzyme 0.2 0
RNAs-free water 4.6 4.8
Total volume 16 16
a
Columns left blank to allow calculations to be inserted
3.7 Reverse 1. For unknown samples, dilute all RNA extracts to be used in a 1:
Transcriptase qPCR 10 ratio in RNase-free water.
2. Mix well by vortexing for 5 s and briefly spin down to remove
any droplets or air bubbles.
3. For standard samples for a standard curve, take the Mtb and EC
RNA standards from the 80 C freezer and thaw at room
temperature. Make seven and six 10-fold dilutions of Mtb and
EC standard samples respectively (see Note 7).
4. Master mix preparation NOTE: Master mix (MM) is a solution
of PCR reagents sufficient to amplify all samples, standards, and
water for a no template control (NTC). The water used as NTC
should be the same water used in the extraction and for
Stool TB MBLA Treatment Monitoring 159
preparing the MM. Ensure that the standards, each RNA sam-
ple, and its decimal dilution are amplified 2 for the reverse
transcriptase positive (RT+) reaction and 1 for the reverse
transcriptase negative (RT) reaction. The RT reaction is a
control to determine the efficiency of DNA removal (see
Note 8).
5. Transfer 16 μL of MM into each PCR reaction tube.
6. Add 4 μL of RNA extract into each RT+ and RT reaction tube
and water into the NTC reaction tubes.
7. Load the reaction tubes into a real-time PCR machine and set
the PCR conditions as follows: 50 C for 30 min, 95 C for
15 min, 40 cycles at 94 C for 45 s, and 60 C for 1 min with
acquisition with fluorophores that absorb in green and yellow
channels.
4 Notes
References
1. Schumacher SG, Denkinger CM (2019) The 6. Musisi E, Ssesolo A, Sloan DJ et al (2022)
impact of Xpert MTB/RIF—do we have a final Detection and quantification of viable mycobac-
answer? The Lancet Global Health 7:e161–e162 terium tuberculosis bacilli in saline-processed
2. Kennedy N, Gillespie SH, Saruni AO et al stool samples by tuberculosis molecular bacterial
(1994) Polymerase chain reaction for assessing load assay: a potential alternative for processing
treatment response in patients with pulmonary stool. Microbiol Spectr 10:e00274–e00222
tuberculosis. J Infect Dis 170(713):716 7. Ruddy M, McHugh TD, Dale JW et al (2002)
3. Friedrich SO, Rachow A, Saathoff E et al (2013) Estimation of the rate of unrecognized cross-
Assessment of the sensitivity and specificity of contamination with mycobacterium tuberculosis
Xpert MTB/RIF assay as an early sputum bio- in London microbiology laboratories. J Clin
marker of response to tuberculosis treatment. Microbiol 40:4100–4104
Lancet Respir Med 1(462):470 8. Mekota A-M, Gillespie SH, Hoelscher M et al
4. Honeyborne I, McHugh TD, Phillips PPJ et al (2022) Building sustainable clinical trial sites in
(2011) Molecular bacterial load assay, a culture- Sub-Saharan Africa through networking, infra-
free biomarker for rapid and accurate quantifica- structure improvement, training and conducting
tion of sputum mycobacterium tuberculosis clinical studies: the PanACEA approach. Acta
bacillary load during treatment. J Clin Microbiol Trop 238:106776
49(3905):3911 9. Mtafya B, Sabi I, John J et al (2022) Systematic
5. Sabiiti W, Azam K, Farmer ECW et al (2020) assessment of clinical and bacteriological mar-
Tuberculosis bacillary load, an early marker of kers for tuberculosis reveals discordance and
disease severity: the utility of tuberculosis molec- inaccuracy of symptom-based diagnosis for
ular bacterial load assay. Thorax 75:606–608 treatment response monitoring. Front Med 9:
992451
Chapter 16
Abstract
Outbreaks are a risk to public health particularly when pathogenic, hypervirulent, and/or multidrug-
resistant organisms (MDROs) are involved. In a hospital setting, vulnerable populations such as the
immunosuppressed, intensive care patients, and neonates are most at risk. Rapid and accurate outbreak
detection is essential to implement effective interventions in clinical areas to control and stop further
transmission. Advances in the field of whole genome sequencing (WGS) have resulted in lowered costs,
increased capacity, and improved reproducibility of results. WGS now has the potential to revolutionize the
investigation and management of outbreaks replacing conventional genotyping and other discrimination
systems. Here, we outline specific procedures and protocols to implement WGS into investigation of
outbreaks in healthcare settings.
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_16,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
161
162 Benjamin J. Parcell et al.
Fig. 1 Recommendations for performance and quality assessment of data and processing in WGS service for
outbreak detection. Quality assessments of data are displayed in dark grey-colored boxes
166 Benjamin J. Parcell et al.
2 Materials
2.1 DNA Extraction • Cell pellet(s) of the bacterial isolate(s) of interest, based on
single colony pick, grown overnight or to late exponential
growth phase if longer in liquid growth medium
(107–109 CFU/mL), obtained by centrifugation (see Note 1).
• DNA extraction kit and associated materials. In this protocol for
the purposes of DNA extraction from Staphylococcus spp. we
suggest the MasterPure Gram Positive DNA Purification Kit
(Cambio), which requires:
– MasterPure Gram Positive DNA Purification Kit (Cambio,
UK).
– TE Buffer.
– Lysostaphin.
– Isopropanol.
– 70% ethanol (freshly prepared).
Fig. 2 An example of a WGS report. The report gives a score on the likelihood of an outbreak and an
explanation of the phylogeny and its clinical implications. Note that for case 1 the patient’s bacteriology
sample isolate was not stored by the microbiology laboratory as at that point in time an outbreak was not
suspected
170 Benjamin J. Parcell et al.
– Vortex.
– Incubator at 37 °C.
– Ice bucket.
2.2 DNA Quality • Qubit Broad Range Assay Kit (ThermoScientific, UK).
Assessment and • Qubit High Sensitivity Assay Kit (ThermoScientific).
Quantification
• Qubit Assay Tubes (ThermoScientific).
• NanoDrop spectrophotometer (ThermoScientific).
• Qubit fluorometer (ThermoScientific).
• Optional—agarose and gel electrophoresis equipment, TAE
Buffer (Tris-acetate-EDTA), high molecular weight DNA lad-
der, DNA loading buffer.
2.3 Library • Nextera XT DNA library preparation kit (Illumina Inc., San
Preparation Diego, CA, USA).
• Nextera XT index kit (Illumina).
• PhiX (Illumina).
• Magnetic DNA-binding beads (such as AMPure XP Beads,
Beckman Coulter, UK).
• Absolute ethanol for molecular biology.
• NaOH.
• PCR grade water.
• 96-well thermal cycler with heated lid.
• Heat block for 1.5 ml centrifuge tubes.
• High-speed microplate shaker.
• Magnetic stand, 96-well format.
• Microplate centrifuge.
• Vortex.
• 96-well PCR plates.
• Multichannel pipettes, single channel pipettes, and sterile
filtered tips.
• 96-well storage plates, round well, 0.8 ml (deep well plate).
• Microseal “B” adhesive seals (BioRad).
• RNase/DNase-free multichannel reagent reservoirs, disposable.
• TruSeq Index Plate Fixture.
• Optional—2100 Bioanalyzer or Tape Station (Agilent, UK).
• Optional—Repeater pipette, 125 ml and 1250 ml syringes.
Application of Pathogen Genomics to Outbreak Investigation 171
2.5 Illumina • MiSeq platform (Illumina Inc., San Diego, CA, USA).
Sequencing • MiSeq Reagent Kit (v2 or v3, Illumina).
3 Methods
3.1 DNA Extraction Following the MasterPure Gram Positive DNA Purification Kit
(Cambio) protocol (see Note 2):
1. Centrifuge 1.0 mL of an overnight Gram-positive bacterial
culture to create a pellet. Discard the supernatant.
2. Add 150 μL of TE Buffer and vortex to resuspend the cell
pellet.
3. Add 1 μL of Ready-Lyse Lysozyme to each resuspended pellet.
4. Incubate at 37 °C for 30 min.
5. Make a mastermix of 1 μL of Proteinase K (50 μg/μL) and
149 μL of Gram-Positive Lysis Solution for each pellet. Add
150 μL of the Proteinase K/Gram Positive Lysis Solution to
each sample and mix thoroughly.
6. Incubate at 65–70 °C for 15 min, vortexing briefly every 5 min.
7. Cool samples to 37 °C and place on ice for 3–5 min.
8. Add 175 μL of MPC Protein Precipitation Reagent to 300 μL
of lysed sample. Vortex for 10 s to mix vigorously. Centrifuge at
4 °C for 10 min at >10,000 × g in a microcentrifuge to pellet
the precipitate. Transfer the supernatant to a clean microcen-
trifuge tube and discard the pellet.
9. Add 1 μL of RNase A (5 μg/μL) to each sample and mix
thoroughly.
10. Incubate at 37 °C for 30 min.
11. Add 500 μL of isopropanol, invert 30–40 times and then
centrifuge at 4 °C for 10 min at >10,000× g in a microcen-
trifuge. Remove the isopropanol without dislodging the DNA
pellet.
12. Rinse the pellet with 70% ethanol. Briefly re-centrifuge if the
pellet becomes dislodged.
13. Resuspend the DNA pellet in 35 μL of TE Buffer.
172 Benjamin J. Parcell et al.
3.2.3 DNA Quantity 1. Prepare Qubit™ solution for samples and 2 standards (1 μL
Qubit reagent and 199 μL Qubit buffer per sample and stan-
dard) and label the tops of the Qubit 0.5 mL tubes accordingly.
2. Distribute the Qubit™ solution into each tube: 190 μL for
tubes containing each standard and 198 μL for sample tubes.
3. Add 10 μL of the standard to its tube and 2 μL of each sample
to its tube.
4. Mix each tube by vortexing for 3–5 s.
5. Incubate at room temperature for 2 min.
6. Read the standards by inserting the tubes into the fluorometer
and after selecting the appropriate procedure (dsDNA BR for
broad range or dsDNA HS for High Sensitivity Qubit kit) press
Read standard.
7. Specify sample volume (2 μL) and units (ng/μL) and read the
fluorescence for each sample.
8. If negative controls are more than 1 ng/μL, repeat the previous
steps to address contamination.
9. We recommend first quantifying each DNA sample using the
Qubit Broad Range Assay Kit, diluting each DNA appropri-
ately if the concentration is too high, then re-quantifying using
the Qubit High Sensitivity Assay Kit. Each DNA sample should
be carefully diluted, in PCR grade water, to 0.4 ng/μL prior to
library preparation (see Note 5).
3.2.4 Sample and For WGS on a MiSeq instrument, it is important to determine how
Sequencing Kit Selection many samples of a given species should be included on a sequencing
run to obtain the maximum amount of quality data. A prepared
library will be sequenced using a MiSeq Reagent Kit (most up-to-
Application of Pathogen Genomics to Outbreak Investigation 173
3.3 Library The Nextera XT library preparation kit uses a single transposase
Preparation enzymatic reaction where sample DNA is simultaneously fragmen-
ted and tagged with sequencing adapters. The library is then ampli-
fied with sample-specific indices, cleaned up, and normalized before
pooling the indexed samples and adding the PhiX control for
sequencing.
3.3.1 Tagmentation 1. Thaw TD and ATM tubes on ice, gently invert 3–5 times and
(see Note 6) spin briefly in microcentrifuge.
2. Pipette 10 μL TD buffer to each well of a labeled 96-well PCR
plate.
3. Using a multichannel pipette, add 5 μL DNA (at 0.4 ng/
μL = 2 ng total) to each sample well of the plate.
4. Add 5 μL ATM mastermix to each well using a multichannel
pipette, and mix gently 5 times by pipette, changing tips
between samples.
5. Cover plate with a PCR seal and centrifuge at 280× g at 20 °C
for 1 min.
6. Incubate in a thermocycler (with heated lid) at 55 °C for 4 min,
followed by a hold step at 10 °C. Once the plate reaches 10 °C,
proceed immediately to the neutralization step.
3.3.2 Neutralization A. Briefly centrifuge the plate, carefully remove the PCR seal, and
add 5 μL chilled NT Buffer to each well using a multichannel
pipette. Pipette gently 5 times to mix, avoiding cross-
contamination.
B. Apply a new PCR seal, briefly centrifuge (280× g at 20 °C for
1 min) and incubate at room temperature for 5 min.
3.3.3 PCR Amplification 1. Thaw NPM and index primers on the bench at room tempera-
(see Note 7) ture, then mix gently by inverting the thawed tubes 3–5 times
and centrifuging them using 1.7 mL microcentrifuge tubes as
adaptors.
174 Benjamin J. Parcell et al.
3.3.4 Post-PCR Clean-Up 1. Remove the magnetic beads from the fridge 30 min prior to
use, to allow them to warm to room temperature (see Note 8).
2. Centrifuge the PCR plate at 280× g at 20 °C for 1 min, remove
the seal and transfer 50 μL of PCR product to a labelled deep-
well plate.
3. Vortex the magnetic beads thoroughly for 30 s until fully
resuspended. Using a repeater pipette, immediately pipette
25 μL of beads into each well of the deep well plate (see Note
9). Aim for the sides of the well to avoid splashing.
4. Shake the deep-well plate on a microplate shaker at 1800 rpm
for 2 min, then incubate at room temperature without shaking
for 5 min.
5. Place the deep well plate on a magnetic stand for 2 min or until
the supernatant has cleared. Using a multichannel pipette,
carefully remove and discard the supernatant (changing tips
between each sample).
6. Keeping the plate on the magnet, pipette 200 μL freshly
prepared 80% ethanol (see Note 10) to each well using a
multichannel pipette. Leave the plate on the magnetic stand
for 30 s, then carefully pipette off and discard the supernatant.
Application of Pathogen Genomics to Outbreak Investigation 175
7. Repeat step 6.
8. To remove any visible residual ethanol, use a P20 multichannel
pipette with fine tips.
9. Keeping the plate on the magnet, allow the pellet to air-dry for
15 min.
10. Remove the plate from magnet and pipette 45 μL RSB into
each well using a repeater pipette. Place the plate on a micro-
plate shaker at 1800 rpm for 2 min.
11. After shaking, incubate the plate at room temperature for
2 min, before returning it on the magnet for a further 2 min
or until the supernatant is clear.
12. Transfer 50 μL of supernatant (this is the cleaned-up library)
from the deep-well plate to a new labeled PCR plate using a
multichannel pipette. Seal the plate and store at -15 to -25 °C
for up to 1 week.
3.3.5 Manual Library 1. Quantify the cleaned libraries using the Qubit High Sensitivity
Normalization and Pooling Assay as previously described in Subheading 3.2.3, to provide
(see Note 11) the concentration of each library in ng/μL.
2. Run 1 μL of the cleaned library on an Agilent High Sensitivity
DNA chip on Bioanalyzer, following standard protocols (see
Note 12), to provide the average fragment size in bp.
3. Calculate the molarity of each library using the following
formula:
Molarity ðnMÞ = ðC=M LÞ 1e6
where
C = concentration (ng/μL, from Qubit)
M = average molar mass (660 g/bp)
L = average fragment length (bp, from Bioanalyzer)
4. Calculate the average concentration of the libraries and pool
them so that each library has a pooled concentration equal to
the average (see Note 13). The following formula is used for
this calculation:
V ind = ðMav 5 μLÞ=Mind
where
Vind = volume of individual library to add to pool (μL)
Mav = average library molarity (nM)
Mind = molarity of individual library (nM, from step C)
5. Keep pooled library in fridge until ready for sequencing. Do
not denature the pooled library until you are ready to begin the
sequencing run (see Note 14).
176 Benjamin J. Parcell et al.
3.5 Preparing the 1. Thaw the reagent cartridge in a basin of lukewarm water. When
MiSeq Instrument immersing the cartridge, do not exceed the “max water level”
mark on the cartridge.
2. In MiSeq Control Software, select “Manage Instrument”,
select “Reboot”.
3. Perform Maintenance Wash: three washes taking approxi-
mately 20 min each.
(a) Prepare fresh wash solution by diluting 25 mL Tween
20 (10%) in 475 mL MilliQ or Elga water.
(b) Add 6 mL wash solution to each reservoir of the wash tray
and transfer the remaining wash solution to the wash
bottle.
(c) In MiSeq Control Software, select “Perform Wash”, select
“Maintenance Wash” and follow the instructions.
4. Open Illumina Experiment Manager, select “Create Sample
Plate” and enter sample IDs and indices. Save and select “Cre-
ate Sample Sheet”. Select “Small Genome Resequencing”,
enter reagent cartridge barcode, and run ID. Set cycle number.
5. Make sure that there are at least 100Gb of free space in the
instrument memory.
6. Prepare fresh 0.2 N NaOH by diluting 200 μL 1 N NaOH in
800 μL sterile MilliQ water.
7. Remove a new flow cell from fridge/cold room. Remove the
flow cell from plastic container, retaining the storage solution
in the container. Flush the flow cell with MilliQ water to wash
off the storage solution, dry carefully with lint-free tissue. Add
several drops of 70% ethanol to lint-free tissue and wipe the
glass part of the flow cell to remove any residue or smears.
Application of Pathogen Genomics to Outbreak Investigation 177
3.6 Denaturing the 1. In a 1.5 μL microtube, add 5 μL pooled library and 5 μL fresh
Pooled Library 0.2 N NaOH.
2. Briefly vortex and centrifuge.
3. Incubate at room temperature for 5 min.
4. Add 990 μL chilled HT1 buffer.
5. Invert several times to mix, then pulse centrifuge.
6. If PhiX control is to be used, either spike in at 1% or at 5% for
low diversity libraries:
Most Libraries (1%):
Denatured and diluted PhiX control 6 μL
Denatured and diluted sample library 594 μL
Low Diversity Libraries (≥ 5%).
Denatured and diluted PhiX control 30 μL
Denatured and diluted sample library 570 μL
7. Pierce the foil seal of the indicated reservoir of the reagent
cartridge with a sterile 1 mL pipette tip and load 600 μL of
denatured library combined with the PhiX library.
3.7 Beginning the 1. Open MiSeq Control Software, select “Sequence” and follow
Sequencing Run instructions to load the flow cell, reagent cartridge, and incor-
poration buffer.
2. Before the run can start, the MiSeq performs several checks:
wait until these have all been performed successfully (indicated
by a green tick), then select “Start Run”.
3. Vibration can negatively affect a MiSeq run. Avoid using any
centrifuges, plate shakers, or vortexes on the same bench as the
MiSeq instrument during the run.
4 Notes
was the least expensive of the kits tested and yielded high
molecular weight DNA with good levels of purity, the Cambio
MasterPure Gram Positive DNA Purification kit gave the high-
est yields for all routine bacterial species. This was particularly
beneficial for time-dependent sequencing of mixed batches of
clinical isolates, containing samples of different species.
3. A quick assessment of the purity of DNA in solution can be
conducted with minimal DNA wastage using a microvolume
spectrophotometer, such as the NanoDrop. The ratio of absor-
bance at 260 and 280 nm provides an insight into the relative
quantity of nucleic acids and proteins, which have absorbance
maxima at 260 nm and 280 nm, respectively. A 260/280 ratio
of approximately 1.8–2.0 is considered an acceptable indicator
of DNA purity. While protein absorbs maximally at 280 nm,
other carryover contaminants such as phenol or guanidine
absorb at 230 nm. Thus, the 260/230 ratio is also checked;
in this case, a ratio of 2.0–2.2 is optimal.
4. Assessing the molecular weight of extracted DNA is of less
importance in short-read sequencing than long-read methods,
where it is a possible determinant of read length. However,
when selecting a DNA extraction method, it is worthwhile
checking the approximate length of extracted DNA fragments.
This can be done cheaply and easily by electrophoresis using a
low percentage agarose gel and identifying excessive degrada-
tion visually under UV light. Alternatively, specialist instru-
ments are available for measuring fragment length of both
genomic DNA and sequencing libraries, such as TapeStation
(Agilent).
5. Quantification of DNA is particularly important for library
preparation. For a sequencing run to yield data that is repre-
sentative of the whole batch of input DNA, all samples must be
of a uniform concentration. Furthermore, the enzymatic pro-
cesses of library preparation require an optimized amount of
input DNA to be present. For these reasons, DNA quantifica-
tion is a major part of the sequencing preparation pipeline.
While spectrophotometry provides a good assessment of
DNA purity, it is strongly advisable to use a DNA-binding
fluorescence-based method for DNA quantification. These
methods use a fluorescent dye which binds to the minor groove
of double-stranded DNA, giving a more specific, sensitive, and
reproducible quantification of DNA concentration. DNA can
then be diluted in sterile PCR-grade water to a uniform and
optimal concentration.
6. As the name suggests, this is a process that combines fragmen-
tation and tagging in a single step. Transposons randomly
fragment the DNA to approximately 200-500 bp lengths,
Application of Pathogen Genomics to Outbreak Investigation 179
Acknowledgments
This work was supported by the Wellcome Trust ISSF award [grant
number 097831/Z/11/Z]; and the Chief Scientist Office through
the Scottish Infection Research Network [SIRN10].
Thanks to Voòng Vı̃nh Phát (OUCRU, Ho Chi Minh City,
Vietnam) for several library preparation optimization suggestions.
References
1. U.S. Department of Health and Human Ser- disease outbreak analysis using the R software.
vices and Centers for Disease Control and Pre- Epidemics 7:28–34
vention (CDC) (2012) Principles of 8. Wellcome Sanger Institute MicroReact
Epidemiology in Public Health Practice, 9. Gladstone RA, Lo SW, Goater R et al (2020)
Third Edition: An Introduction to Applied Visualizing variation within Global Pneumo-
Epidemiology and Biostatistics, coccal Sequence Clusters (GPSCs) and country
2. National Services Scotland (NSS) Chapter 3 – population snapshots to contextualize pneu-
Healthcare Infection Incidents, outbreaks and mococcal isolates. Microb Genom 6:e000357
Data Exceedance. National Infection Preven- 10. Woolhouse M (2017) Quantifying
tion and Control Manual. https://www.nipcm. Transmission 5
hps.scot.nhs.uk/chapter-3-healthcare-infec
tion-incidents-outbreaks-and-data- 11. Johnson E, Raza S University of Cambridge –
exceedance/ What is long read sequencing? https://www.
phgfoundation.org/briefing/what-is-long-
3. Mayhall CG (2011) Hospital Epidemiology read-sequencing
and Infection Control
12. Quick J, Loman NJ, Duraffour S et al (2016)
4. Baxter ICNET Clinical Surveillance Software. Real-time, portable genome sequencing for
https://www.icnetsoftware.com/ Ebola surveillance. Nature 530:228–232
5. Baker MA, Yokoe DS, Stelling J et al (2020) 13. de Jesus JG, Giovanetti M, Rodrigues Faria N,
Automated outbreak detection of hospital- Alcantara LCJ (2019) Acute vector-borne viral
associated pathogens: value to infection pre- infection: Zika and MinION surveillance.
vention programs. Infect Control Hosp Epide- Microbiol Spectr 7(4). https://doi.org/10.
miol 41:1016–1021 1128/microbiolspec.AME-0008-2019.
6. Smith CM, Allen DJ, Nawaz S et al (2019) An PMID: 31400093.
interactive data visualisation application to 14. Cheng VCC, Ip JD, Chu AWH et al (2022)
investigate nosocomial transmission of infec- Rapid spread of severe acute respiratory syn-
tions. Wellcome Open Res 4:100 drome coronavirus 2 (SARS-CoV-2) Omicron
7. Jombart T, Aanensen DM, Baguelin M et al subvariant ba.2 in a single-source community
(2014) OutbreakTools: a new platform for outbreak. Clin Infect Dis 75:e44–e49
182 Benjamin J. Parcell et al.
15. Bai N University of California San Francisco- 28. Koboldt DC, Chen K, Wylie T et al (2009)
NASA’s DNA Sequencing in Space is a Success, VarScan: variant detection in massively parallel
UCSF Researchers Confirm. https://www. sequencing of individual and pooled samples.
ucsf.edu/news/2016/08/404006/nasas- Bioinformatics 25:2283–2285
dna-sequencing-space-success-ucsf- 29. Deatherage DE, Barrick JE (2014) Identifica-
researchers-confirm#:~:text=The%20pocket% tion of mutations in laboratory-evolved
2Dsized%20MinION%20employs,to%20 microbes from next-generation sequencing
decode%20the%20DNA%20sequence data using breseq. Methods Mol Biol 1151:
16. Fan H, Ives AR, Surget-Groba Y et al (2015) 165–188
An assembly and alignment-free method of 30. Seemann T (2014) Prokka: rapid prokaryotic
phylogeny reconstruction from next- genome annotation. Bioinformatics 30:2068–
generation sequencing data. BMC Genomics 2069
16:522 31. Center for Genomic Epidemiology CSI Phy-
17. MacCannell D (2019) Platforms and analytical logeny. https://cge.food.dtu.dk/services/
tools used in nucleic acid sequence-based CSIPhylogeny/
microbial genotyping procedures. Microbiol 32. Leekitcharoenphon P, Kaas RS, Thomsen MCF
Spectr 7:AME-AME-0005-2018 et al (2012) snpTree--a web-server to identify
18. Luheshi L, Raza S, Moorthie S, et al (2015) and construct SNP trees from whole genome
PHG foundation- pathogen genomics into sequence data. BMC Genomics 13(Suppl 7):S6
practice 33. Gustavsson P, Syberfeldt A (2018) A new algo-
19. Deurenberg RH, Bathoorn E, Chlebowicz MA rithm using the non-dominated tree to
et al (2017) Application of next generation improve non-dominated sorting. Evol Comput
sequencing in clinical microbiology and infec- 26:89–116
tion prevention. J Biotechnol 243:16–24 34. Bertels F, Silander OK, Pachkov M et al (2014)
20. Bankevich A, Nurk S, Antipov D et al (2012) Automated reconstruction of whole-genome
SPAdes: a new genome assembly algorithm and phylogenies from short-sequence reads. Mol
its applications to single-cell sequencing. J Biol Evol 31:1077–1088
Comput Biol 19:455–477 35. Naccache SN, Federman S, Veeraraghavan N
21. Zeden MS, Gründling A (2023) Bacterial et al (2014) A cloud-compatible bioinformatics
whole-genome-resequencing analysis: basic pipeline for ultrarapid pathogen identification
steps using the clc genomics workbench soft- from next-generation sequencing of clinical
ware, vol 2023. Cold Spring Harb Protoc, samples. Genome Res 24:1180–1192
p pdb.prot107901 36. Galaxy. https://usegalaxy.org/
22. Zerbino DR, Birney E (2008) Velvet: algo- 37. Caboche S, Even G, Loywick A et al (2017)
rithms for de novo short read assembly using MICRA: an automatic pipeline for fast charac-
de Bruijn graphs. Genome Res 18:821–829 terization of microbial genomes from high-
23. Gardner SN, Slezak T, Hall BG (2015) throughput sequencing data. Genome Biol
kSNP3.0: SNP detection and phylogenetic 18:233
analysis of genomes without genome align- 38. Quijada NM, Rodrı́guez-Lázaro D, Eiros JM
ment or reference genome. Bioinformatics 31: et al (2019) TORMES: an automated pipeline
2877–2878 for whole bacterial genome analysis. Bioinfor-
24. Inouye M, Dashnow H, Raven L-A et al matics 35:4207–4212
(2014) SRST2: rapid genomic surveillance for 39. Chen L, Yang J, Yu J et al (2005) VFDB: a
public health and hospital microbiology labs. reference database for bacterial virulence fac-
Genome Med 6:90 tors. Nucleic Acids Res 33:D325–D328
25. Li H, Durbin R (2009) Fast and accurate short 40. Altschul SF, Madden TL, Sch€affer AA et al
read alignment with Burrows-Wheeler trans- (1997) Gapped BLAST and PSI-BLAST: a
form. Bioinformatics 25:1754–1760 new generation of protein database search pro-
26. Langmead B, Trapnell C, Pop M et al (2009) grams. Nucleic Acids Res 25:3389–3402
Ultrafast and memory-efficient alignment of 41. Florensa AF, Kaas RS, Clausen PTLC et al
short DNA sequences to the human genome. (2022) ResFinder – an open online resource
Genome Biol 10:R25 for identification of antimicrobial resistance
27. McKenna A, Hanna M, Banks E et al (2010) genes in next-generation sequencing data and
The genome analysis toolkit: a mapreduce prediction of phenotypes from genotypes.
framework for analyzing next-generation Microb Genom 8
DNA sequencing data. Genome Res 20: 42. Gupta SK, Padmanabhan BR, Diene SM et al
1297–1303 (2014) ARG-ANNOT, a new bioinformatic
Application of Pathogen Genomics to Outbreak Investigation 183
tool to discover antibiotic resistance genes in 48. Wellcome Connecting Science Courses and
bacterial genomes. Antimicrob Agents Che- Conferences -Genomics and Clinical Microbi-
mother 58:212–220 ology. https://coursesandconferences.
43. Alcock BP, Huynh W, Chalil R et al (2023) wellcomeconnectingscience.org/event/geno
CARD 2023: expanded curation, support for mics-and-clinical-microbiology-2024-
machine learning, and resistome prediction at 20240114/
the comprehensive antibiotic resistance data- 49. Crisan A, McKee G, Munzner T et al (2018)
base. Nucleic Acids Res 51:D690–D699 Evidence-based design and evaluation of a
44. Parcell BJ, Gillespie SH, Pettigrew KA et al whole genome sequencing clinical report for
(2021) Clinical perspectives in integrating the reference microbiology laboratory. PeerJ
whole-genome sequencing into the investiga- 6:e4218
tion of healthcare and public health outbreaks – 50. Ferdinand AS, Kelaher M, Lane CR et al
hype or help? J Hosp Infect 109:1–9 (2021) An implementation science approach
45. Parcell BJ, Oravcova K, Pinheiro M et al to evaluating pathogen whole genome
(2018) Pseudomonas aeruginosa intensive sequencing in public health. Genome Med
care unit outbreak: winnowing of transmissions 13:121
with molecular and genomic typing. J Hosp 51. MasterPure Gram Positive DNA Purification
Infect 98:282–288 Kit Manual. https://cambio.co.uk/library/
46. (2022) WHO guiding principles for pathogen images/html_images/2022%20protocols/
genome data sharing. Presented at the manual_MGP04100_MasterPure-Gram-
47. Johnson S, Parker M (2020) Ethical challenges Positive-DNA.pdf
in pathogen sequencing: a systematic scoping
review. Wellcome Open Res 5:119
Chapter 17
Abstract
Whole genome sequencing (WGS) is becoming an important diagnostic tool for antimicrobial susceptibility
testing of Mycobacterium tuberculosis complex (MTBC) isolates in many countries. WGS protocols usually
start with the preparation of a DNA-library: the critical first step in the process. A DNA-library represents
the genomic content of a DNA sample and consists of unique short DNA fragments. Although available
DNA-library protocols come with manufacturer instructions, details of the entire process, including quality
controls, instrument parameters, and run evaluations, often need to be developed and customized by each
laboratory to implement WGS technology effectively. Here, we provide a detailed workflow for a
DNA-library preparation based on an adapted Illumina protocol optimized for the reduction of reagent
costs.
Key words WGS, Whole genome sequencing, DNA-library, Mycobacterium tuberculosis complex
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_17,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
185
186 Vanessa Mohr et al.
[2]. In 2021, the WHO presented the first catalog of mutations and
their association with antimicrobial resistance which can be inter-
rogated with next generation sequencing (NGS) technology [3].
The first step in the NGS workflow for WGS is the preparation
of a DNA-library. For short-read NGS applications, DNA quality
such as fragmentation is of less concern, but DNA molecular con-
centration needs to be determined accurately, e.g., by fluorometric
assays (see Materials below). A DNA-library is a collection of short
(~500–1200 bp) random DNA fragments representing the geno-
mic content of a biological DNA sample. Although the DNA-
library preparation is a crucial process that can influence the quality
of the sequencing data, there is usually little detailed information in
scientific publications on quality controls, instrument parameters,
and run evaluations.
In fact, there are several library-preparation protocols available,
with state-of-the-art kits requiring as little as 1 ng of input DNA.
For most protocols, the initial DNA undergoes a random fragmen-
tation which can occur via mechanical or enzymatic shearing.
Enzymatic shearing is often the method of choice as it usually
requires less input material, no additional laboratory equipment,
has a significantly shorter turnaround time, and low-cost protocols
can be applied [4]. After or during the shearing procedure, short
oligonucleotides are ligated/attached to the DNA fragments.
These short flanking sequences function as primer binding sites
for the subsequent PCR amplification during which indices (for
sample tracking) and adapters (for attachment to the flow cell, and
sequencing reaction) are integrated into the final DNA-library.
As a showcase method for this chapter, we will present a previ-
ously published library preparation protocol based on Illumina
reagents [4], offering excellent performance at relatively low
costs. Likewise, the NGS sequencing run and run evaluation will
be exemplarily described using the Illumina NextSeq 500 instru-
ment. Most of the workflow principles can be easily transferred to
other common library preparation solutions. For a list of abbrevia-
tions and definitions, see Table 1.
This chapter should be seen in combination with Chap. 18
which presents the workflow for bioinformatic analysis of raw
sequencing data (fastq files) and the detection of mutations impli-
cated in antimicrobial resistance.
2 Materials
Table 1
Abbreviations and definitions
Abbreviation Comment
Fastq Raw sequencing data in text format with base quality values
DNA-library Random collection of DNA fragments
Index Unique sequence used to label a NGS library for multiplexing
Run Multiplexing approach on a DNA sequencing machine
NGS Next generation sequencing
WGS Whole genome sequencing (an NGS application)
MDR-TB Multidrug-resistant tuberculosis
MTBC Mycobacterium tuberculosis complex
Tagmentation The process of tagging and enzymatic fragmentation of genomic DNA
TDE1 Tagment DNA enzyme
TD buffer Tagment DNA buffer
UD index Unique dual index
EB Elution buffer
KAPA A ready-to-use PCR master mix
HT1 Hybridization buffer
PhiX A ready-to-use control library of phage DNA
4. Vortex.
5. Centrifuge for 96-well plates and 1,5 mL microcentrifuge
tubes.
6. Tag DNA Enzyme & Buffer Large Kit (Illumina) (includes
TDE1 and TD buffer). This should be stored at 20 C.
7. IDT for Illumina—Nextera DNA UD Index Set A, B, C, D,
stored at 20 (each kit includes 96 unique index combina-
tions for altogether 192 samples).
8. KAPA HiFi Library Amplification Kit (Roche), stored in
1.5 mL aliquots at 20 C. Serves as PCR mastermix.
9. Magnetic rack like DynaMag™-PCR Magnet (Thermo Fisher
Scientific) and magnetic beads, e.g., MagSi-NGS Prep Plus
(Steinbrenner).
10. EB (5 mM Tris/HCL, pH 8.5) or molecular biology grade
water for elutions and dilutions. Avoid buffers with EDTA, this
will negatively interact with the sequencing chemistry (see
Note 1).
188 Vanessa Mohr et al.
2.2 DNA Sequencing 1. NGS instrument, e.g., NextSeq 500 (see Note 2).
2. Sequencing reagents suitable for the sequencer and sample
number, e.g., NextSeq v2.5 Mid Output 300 Cycles or Next-
Seq v2.5 High Output 300 Cycles.
3. 200 mM Tris-HCL pH 7.
3 Methods
Table 2
Tagmentation program on a PCR machine
55 C Hold
55 C 10 min
10 C Cool down
Table 3
PCR program adding indices and adapters to finalize the DNA-library
72 C 3 min
98 C 5 min
98 C 10 s 14 cycles
63 C 30 s
72 C 30 s
72 C 5 min
10 C Cool down
3.2 DNA Sequencing 1. Plan the sequencing run by selecting the libraries to be multi-
on the NextSeq 500 plexed (see Note 12).
Platform 2. Thaw reagent cartridge in a water bath at room temperature
(~2–3 h).
3. Prepare the sample sheet in the NextSeq 500 local run manager
(see Note 13).
4. Dilute all libraries to 2 nM with EB (see Note 14).
5. 30 Minutes prior to starting, bring the flow cell to room
temperature.
6. Combine 2.5 μL of each library with the same molar concen-
tration in a new tube.
7. Combine in a new tube 10 μL of 2 nM Pool and 10 μL of 0.2 M
NaOH and incubate for 5 min at room temperature to achieve
denaturation.
8. Add 10 μL of 200 mM Tris-pH 7.0 and 970 μL HT1.
9. Dilute the pool to its final loading concentration and add PhiX
(see Note 15).
Use of Whole Genome Sequencing for Mycobacterium tuberculosis. . . 191
3.3 Run Evaluation of In this section one ensures that the products from the run have
a NextSeq 500 Run good quality and meet the pre-defined expected specifications.
Evaluating run performance and sequence data quality is a complex
endeavor. In the following, we list some crucial steps especially for
runs on Illumina instruments.
1. Load the run into the Illumina Sequencing Analysis Viewer
(SAV).
2. On the first tab (Analysis), review the Q30 value. This value
should be within specifications.
3. Open the third tab (Summary). The total yield should be
within specifications for the run chemistry, i.e., for the NextSeq
500 100–120 Gb for a High-Output kit and 32.5–39 Gb for a
Mid-Output kit.
4. Also in the third tab, check whether the density on the flow cell
is in the recommended range, i.e., between 170 and 220 K/
mm2 for the NextSeq 500.
5. Still in the third tab, the value of PF % indicates the number of
flow cell clusters passing quality filters.
6. Open the fourth tab (Indexing). The bars indicate the amount
of sequencing data attributed to each input library. When
DNA-libraries were loaded with equal molarities, number of
reads/bp per sample should be evenly distributed. Obvious
outliers indicate errors during the dilution step to 2 nM (see
Subheading 3.2, step 4).
192 Vanessa Mohr et al.
4 Notes
References
Abstract
Whole genome sequencing of Mycobacterium tuberculosis complex (MTBC) isolates has been shown to
provide accurate predictions for resistance and susceptibility for many first- and second-line anti-tubercu-
losis drugs. However, bioinformatic pipelines and mutation catalogs to predict antimicrobial resistances in
MTBC isolates are often customized and detailed protocols are difficult to access. Here, we provide a step-
by-step workflow for the processing and interpretation of short-read sequencing data and give an overview
of available analysis pipelines.
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_18,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
195
196 Viola Dreyer et al.
has become the de facto standard for NGS data analysis of MTBC
strains [13]. Reaching a final resistance profile from read data
requires a dedicated bioinformatics analysis covering several key
steps:
• Sequence data preprocessing and quality control
• Reference mapping and variant detection
• Annotating detected variants with a resistance marker database
and reporting
Initial sequence data can optionally be subjected to rigorous
quality controls before proceeding with the main analysis pipeline.
These quality controls can include detection of contaminating
non-MTBC sequence information, overall sequence data accuracy,
and overrepresentation of specific oligomers in the tails of the reads.
As modern reference mapping tools can deal with low-sequence
quality regions in the start and end tails of the reads, quality
trimming to cut those regions is usually not necessary.
In the reference mapping step, a suitable program such as BWA
[14] aligns raw read data (fastq files) to a given reference genome.
While there are some potential problems such as reads from highly
repetitive genomic regions, reads of low sequencing accuracy and
hybrid or contaminated reads, this step is relatively simple and
modern algorithms are suitably sophisticated. Routine quality
metrics should also be determined for the resulting reference
mapping. Important quality metrics are average coverage depth,
i.e., the number of reads mapped per position of the genome, and
coverage breadth, i.e., the amount of the reference genome covered
by sequence data (see Note 1). Comparing the sequence of mapped
reads against the reference sequence, e.g., with GATK [15], allows
for the detection of mutations including short indels. It is impor-
tant to note that these are called relative to the reference genome,
and only analysis results using the same reference genome are
therefore directly comparable. For inferring resistance profiles espe-
cially for MDR-MTBC strains, variant detection needs to be as
comprehensive as possible. Pipelines should report a confident
estimate of detected variants and indicate problems with variant
calling caused by insufficient data quality, at least for resistance-
associated regions of the genome. Special consideration must be
given to the detection of heterozygous/mixed variants present in a
minority fraction of the reads at a given position. These can arise
from co-infections with two distinct strains, or from clonal subpo-
pulations evolving in patients during the course of the treatment.
To allow a valid interpretation of resistance-associated variants,
considerable efforts are ongoing to build up databases containing
validated resistance mutations that can then be used to generate
resistance reports from WGS data.
198 Viola Dreyer et al.
2 Methods
2.1 Identifying Drug The following protocol explains how to install MTBseq [12] via
Resistance Conda/Miniconda and provides details for individual data proces-
Determinants with sing steps on a computer running the Linux operating system
MTBseq (tested for Ubuntu 16.04 LTS).
1. Install Conda or Miniconda hereafter install MTBseq with:
Fig. 1 Schematic representation of the MTBseq workflow. Modules encapsulating specific functionality are
shown in blue boxes
end), this module will align the reads of one FastQ-file with
single-end reads or two fastq-files with paired-end reads to a
reference genome. (see Note 3)
Input:
[SampleID]_[LibID]_[*]_[Direction].fastq.gz
Output:
Bam/[SampleID]_[LibID]_[*].bam
Bam/[SampleID]_[LibID]_[*].bai
Bam/[SampleID]_[LibID]_[*].bamlog
3. Refinement of the initial mapping:
To improve the initial mapping, MTBseq employs GATK
for a realignment around small insertions and deletions and a
base call recalibration with default parameters (see Note 4)
[15]:
For base call recalibration, MTBseq employs a set of known
variants to be excluded. The calibration list is stored in the
MTBseq github directory “var/res/MTB_Base_Calibration_-
List.vcf” of the package.
Input:
Bam/[SampleID]_[LibID]_[*].bam
200 Viola Dreyer et al.
Output:
GATK_Bam/[SampleID]_[LibID]_[*].gatk.bam
GATK_Bam/[SampleID]_[LibID]_[*].gatk.bai
GATK_Bam/[SampleID]_[LibID]_[*].gatk.bamlog
GATK_Bam/[SampleID]_[LibID]_[*].gatk.grp
GATK_Bam/[SampleID]_[LibID]_[*].gatk.intervals
4. Summarize base calls for mapped reads
Table 1
Column names with explanations for the file: “Mapping_and_Variant_Statistics.tab”
Column Explanation
Date The date of MTBseq execution
SampleID The sample ID
LibraryID The library ID
FullID Complete dataset name
Total reads The total amount of sequenced reads
Mapped reads (%) The percentage of reads mapped to the reference genome
Genome size The size of the reference genome
Genome GC The GC content of the reference genome
(Any) Total bases (%) Percentage of the reference genome covered by reads
(Any) GC-content GC content of the reference genome covered by reads
(Any) coverage mean Mean coverage depth
(Any) coverage median Median coverage depth
(Unambiguous) Total bases (%) Percentage of the reference genome covered unambiguously
(Unambiguous) GC-content GC content of the reference genome covered unambiguously
(Unambiguous) coverage mean Mean coverage depth of unambiguously covered positions
(Unambiguous) coverage median Median coverage depth of unambiguously covered positions
SNPs Number of detected SNPs
Deletions Number of detected deletions
Insertions Number of detected insertions
Uncovered Positions of the reference genome not covered by a read
Substitutions (including stop codons) Number of substitutions within genes
Use of Whole Genome Sequencing for Mycobacterium tuberculosis. . . 203
Input:
Bam/[SampleID]_[LibID]_[*].bam
Position_Tables/[SampleID]_[LibID]_[*].gatk_position_ta-
ble.tab
Output:
Statistics/Mapping_and_Variant_Statistics.tab
9. Cleaning the data set:
After all samples have been processed the user should check
the Mapping_and_Variant_Statistics.tab file. As it is a tabular
file, it is possible to open it, e.g., with Microsoft Excel or any
other office tools. The most important values are (see also
Note 1):
• Mapped Reads (%), i.e., low percentage can indicate
contaminations
• (Any) Coverage mean, i.e., the overall sequencing depth
• (Unambiguous) Total Bases (%), i.e., the overall sequencing
breadth
10. Extract mutation implicated in drug resistance
Open the variant file in the “Called” folder from the
low-frequency variant detection workflow (point 6, [lowfreq_-
vars].tab) for the sample of interest and apply the following
filters:
(a) Filter the column “ResistanceSNP” for entries, including
positions marked as phylogenetic or resistance mediating
based on the provided input catalog (see Note 8).
(b) Filter the column “InterestingRegion” for additional (yet
unknown) mutations putatively implicated in antimicro-
bial resistance but with an unclear phenotype.
Fig. 2 Example for an antibiotic resistance report form for Mycobacterium tuberculosis complex isolates
2.3 Comparative The MTBseq tool can be also employed for comparative analysis of
Genomic Analysis multiple datasets. However, we have focused here on the steps to
generate a resistance report. The entire pipeline can be run with the
following command:
2.4 Alternative Tools In recent years, several software solutions have been developed for a
for MTBC mostly automated data analysis (with pre-specified variant calling
Antimicrobial parameters) of MTBC NGS data for resistance inference. They
Resistance Prediction differ in their implementation, with some tools representing local
installed software suites and others built as cloud-based solutions.
In the following, we represent a short overview listing several freely
available tools. A crucial point to keep in mind here is which
mutation catalog is used by the software for determining resistance
patterns as this is the major factor for differences in the resulting
drug resistance report [21]. Below, we summarize some of the
advantages and disadvantages.
1. Mykrobe (Predictor TB; https://www.mykrobe.com/):
• No registration necessary.
• Runs offline on a standard Desktop PC or Laptop.
• Desktop Apps available for Mac and Windows 64-bit.
• No Batch mode available.
• Supports Illumina as well as Oxford nanopore data.
• Drugs included: isoniazid, rifampicin, ethambutol, pyrazi-
namide, ofloxacin, moxifloxacin, ciprofloxacin, streptomy-
cin, amikacin, capreomycin, kanamycin.
• Species identification of four MTBC species (M. tuberculosis,
M. bovis, M. africanum, and M. caprae) and 40 nontubercu-
lous mycobacteria (NTM) species (e.g., M. abscessus,
M. avium, and M. intracellulare).
• Sensitivity 82.6%; Specificity 98.5% (independent validation
set n = 1609).
2. TB Profiler (https://tbdr.lshtm.ac.uk/):
• No registration needed.
• Runs as an online service (command-line version available).
• Batch mode available.
• Supports Illumina and Oxford nanopore data.
• Tbdb database (https://github.com/jodyphelan/tbdb)
used as mutation library.
206 Viola Dreyer et al.
2.5 Tool for the Besides M. tuberculosis, there are also several other mycobacterial
Prediction of Antibiotic species that can cause serious disease in humans [22]. These
Resistances in Non- non-tuberculous mycobacteria (NTM) have gained increased
tuberculous attention in the past decades due to reports of emerging global
Mycobacteria (NTMs) infections, increasing number of deaths, and the occurrence of
several outbreaks [22, 23]. NTMs are notoriously difficult to treat
due to intrinsic resistance to a variety of antibiotics [24]. Further-
more, strains can also acquire additional resistances due to muta-
tions in resistance-related genes [22, 24]. Currently, only two
bioinformatic tools are available for whole genome sequencing-
based resistance prediction of NTMs: Mab_ariba (https://github.
com/samlipworth/Mab_ariba) and NTMprofiler (https://github.
com/jodyphelan/ntmdb/tree/main/db/Mycobacteroides_
abscessus). At present, these tools only include mutation catalogs
for resistance prediction of M. abscessus strains to macrolides and
aminoglycosides. Future efforts should therefore focus on the
establishment of additional mutation catalogues for other NTM
species as well as validation of these catalogues using both in vitro
drug susceptibility data, genotypic data and in vivo treatment
outcomes.
3 Notes
References
1. Merker M, Kohl TA, Barilar I et al (2020) 2:e604–e616. https://doi.org/10.1016/
Phylogenetically informative mutations in S2666-5247(21)00175-0
genes implicated in antibiotic resistance in 10. Sonnenkalb L, Carter JJ, Spitaleri A et al
mycobacterium tuberculosis complex. (2023) Bedaquiline and clofazimine resistance
Genome Med 12:27. https://doi.org/10. in mycobacterium tuberculosis: an in-vitro and
1186/s13073-020-00726-5 in-silico data analysis. Lancet Microbe S2666-
2. Allix-Béguec C, Arandjelovic I, Bi L, Beckert P, 5247(23):00002–00002. https://doi.org/10.
Bonnet M, Bradley P et al (2018) Prediction of 1016/S2666-5247(23)00002-2
susceptibility to first-line tuberculosis drugs by 11. Meehan CJ, Goig GA, Kohl TA et al (2019)
DNA sequencing. N Engl J Med 379:1403– Whole genome sequencing of mycobacterium
1 4 1 5 . h t t p s : // d o i . o r g / 1 0 . 1 0 5 6 / tuberculosis: current standards and open
NEJMoa1800474 issues. Nat Rev Microbiol 17:533–545.
3. Walker TM, Kohl TA, Omar SV et al (2015) https://doi.org/10.1038/s41579-019-
Whole-genome sequencing for prediction of 0214-5
mycobacterium tuberculosis drug susceptibil- 12. Kohl TA, Utpatel C, Schleusener V et al (2018)
ity and resistance: a retrospective cohort MTBseq: a comprehensive pipeline for whole
study. Lancet Infect Dis 15:1193–1202. genome sequence analysis of mycobacterium
https://doi.org/10.1016/S1473-3099(15) tuberculosis complex isolates. PeerJ 6:e5895.
00062-6 https://doi.org/10.7717/peerj.5895
4. Miotto P, Tessema B, Tagliani E et al (2017) A 13. Consortium TCr (2022) A data compendium
standardised method for interpreting the asso- associating the genomes of 12,289 mycobacte-
ciation between mutations and phenotypic rium tuberculosis isolates with quantitative
drug resistance in mycobacterium tuberculosis. resistance phenotypes to 13 antibiotics. PLoS
Eur Respir J 50:1701354. https://doi.org/10. Biol 20:e3001721. https://doi.org/10.1371/
1183/13993003.01354-2017 journal.pbio.3001721
5. Grobbel H-P, Merker M, Köhler N et al (2021) 14. Li H, Durbin R (2009) Fast and accurate short
Design of multidrug-resistant tuberculosis read alignment with burrows-wheeler trans-
treatment regimens based on DNA sequenc- form. Bioinforma Oxf Engl 25:1754–1760.
ing. Clin Infect Dis Off Publ Infect Dis Soc https://doi.org/10.1093/bioinformatics/
Am 73:1194–1202. https://doi.org/10. btp324
1093/cid/ciab359 15. McKenna A, Hanna M, Banks E et al (2010)
6. Walker TM, Miotto P, Köser CU et al (2022) The genome analysis toolkit: a MapReduce
The 2021 WHO catalogue of mycobacterium framework for analyzing next-generation
tuberculosis complex mutations associated DNA sequencing data. Genome Res 20:
with drug resistance: a genotypic analysis. Lan- 1297–1303. https://doi.org/10.1101/gr.
cet Microbe 3:e265–e273. https://doi.org/ 107524.110
10.1016/S2666-5247(21)00301-3 16. Finci I, Albertini A, Merker M et al (2022)
7. Kadura S, King N, Nakhoul M et al (2020) Investigating resistance in clinical mycobacte-
Systematic review of mutations associated rium tuberculosis complex isolates with geno-
with resistance to the new and repurposed mic and phenotypic antimicrobial susceptibility
mycobacterium tuberculosis drugs bedaqui- testing: a multicentre observational study. Lan-
line, clofazimine, linezolid, delamanid and pre- cet Microbe 3:e672–e682. https://doi.org/
tomanid. J Antimicrob Chemother 75:2031– 10.1016/S2666-5247(22)00116-1
2043. https://doi.org/10.1093/jac/dkaa136 17. MTBseq (2022) MTBseq. In: MTBseq
8. Battaglia S, Spitaleri A, Cabibbe AM et al Github. https://github.com/ngs-fzb/
(2020) Characterization of genomic variants MTBseq_source. Accessed 20 Mar 2023
associated with resistance to bedaquiline and 18. Homolka S, Projahn M, Feuerriegel S et al
delamanid in naive mycobacterium tuberculo- (2012) High resolution discrimination of clini-
sis clinical strains. J Clin Microbiol 58:e01304– cal mycobacterium tuberculosis complex
e01320. https://doi.org/10.1128/JCM. strains based on single nucleotide polymorph-
01304-20 isms. PLoS ONE 7(7):e39855. https://doi.
9. Ismail N, Rivière E, Limberis J et al (2021) org/10.1371/journal.pone.0039855
Genetic variants and their association with phe- 19. Coll F, McNerney R, Guerra-Assunc JA et al
notypic resistance to bedaquiline in mycobac- (2014) A robust SNP barcode for typing myco-
terium tuberculosis: a systematic review and bacterium tuberculosis complex strains, Nat
individual isolate data analysis. Lancet Microbe
210 Viola Dreyer et al.
Abstract
Genomic sequencing has revolutionized microbial typing methods and transformed high-throughput
methods in reference, clinical, and research laboratories. The detection of antimicrobial-resistant (AMR)
determinants using genomic methods can provide valuable information on the emergence of resistance.
Here we describe an approach to detecting AMR determinants using an open access and freely available
platform which does not require bioinformatic expertise.
Key words Antimicrobial resistance (AMR) determinants, Center for genomic epidemiology (CGE),
Sequencing, KMER, MLST, ResFinder
1 Introduction
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8_19,
© The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature 2024
211
212 Marie Anne Chattaway
14]. Despite multiple research studies and the potential for predict-
ing phenotypic resistance [15–18], further validation is required
before this concept is implemented in clinical laboratories.
The routine use of sequencing technologies and open sharing
of generated data in reference laboratories [4, 19, 20], and research
and clinical laboratories [21] has resulted in a vast amount for
analysis. High-throughput data were initially processed using
bespoke bioinformatic pipelines or specialized bioinformaticians
to process the data. However, open access platforms have now
been developed enabling non-specialized bioinformaticians to pro-
cess raw sequence data files for identification, typing, and detection
of AMR determinants. There are many platforms available for
detection of AMR determinants [1], which will have their own
manuals and publications, and some will be bespoke to a pathogen
of interest, so it is worth researching the options available to suit
your needs. In this chapter, we will describe an AMR platform that
is not pathogen-specific. The aim of the Center for Genomic Epi-
demiology (CGE) platform (http://www.genomicepidemiology.
org/) is to provide access to bioinformatics resources also for those
with limited experience and thereby allow all countries, institutions,
and individuals to take advantage of sequencing technologies. CGE
is entirely non-commercial and operates a number of free online
bioinformatics services. Ideally, you would have an identification
and basic typing of your organism of interest before assessing AMR
determinants. If not, the advantage of using sequence data is that
once you have the data, you can extrapolate different types of
information such as identification to a species level and the
Sequence Type (see Note 1) which are also available on the CGE
platform and will also be covered in this chapter. Bioinformatic
pipelines available on CGE include: the KmerFinder, k-mers are
substrings of length k contained within a biological sequence and
k-mers can be used in the context of computational genomics and
sequence analysis. There is an R package called kmer that enables
fast alignment-free clustering of biological sequences [22]. For the
context of identification of a bacterial pathogen to the species or
sub-species level, a selection of kmers within the large genome have
been selected to develop a publically available database for identifi-
cation which is called KmerFinder. This approach is computation-
ally quicker and more practicable than querying the entire genome.
Further subtyping of a bacterial pathogen can then be performed
using multilocus sequence typing (MLST) and each scheme is
pathogen specific. MLST is an unambiguous procedure for char-
acterizing isolates of bacterial species using the sequences of inter-
nal fragments of (usually) seven house-keeping genes. For each
house-keeping gene, the different sequences present within a bac-
terial species are assigned as distinct alleles and, for each isolate, the
alleles at each of the seven loci define the allelic profile or sequence
type (ST) [23–26]. Using MLST will group strains into complexes
or eBURST groups (Salmonella) and with some bacterial species
Analysis of Whole Genome Sequencing Data for Detection of Antimicrobial. . . 213
2 Materials
3 Methods
3.1 Getting Started 1. CGE runs entirely online and be accessed via the link https://
www.genomicepidemiology.org/.
2. The main page contains an overview of the platform, if you
select the service tab, it will bring you to the full list of services
(see Figs. 1 and 2).
3.2 Identification Unlike some genomic platforms, you do not need to know the
and Typing of Data presumptive identification of your strain to assess the data. The use
of the KMER ID programme so a useful tool to identify your
organism of interest (see Note 2). In order to identify your strain,
you will need to (1) select the program of interest and then
(2) upload your reads, the options given include FASTA files,
FASTQ files, or assembled files.
3.2.1 KmerFinder— 1. Select KmerFinder under the Typing section on the services
Prediction of Bacterial page (see Notes 3 and 4) (Fig. 3).
Species Using a Fast K-Mer 2. Select the host database—bacterial organism.
Algorithm
3. Other options such as scoring method, gene database, identify
threshold, and depth are set at default; these settings are
acceptable for queries but have the option to be changed if
needed.
4. Browse for your files by clicking on the Isolate File Icon.
5. Press upload once you have added your files.
214 Marie Anne Chattaway
6. Enter your email address, the link to your results will be sent
once processed. If there is a delay in receiving your results,
remember to check your junk mail.
7. Once the results are ready, click on the link sent to you via
email.
8. View your data and look at quality metrics to assign identifica-
tion, the first hit will be the best match with the data submitted
and will be your identification. In this example, the result is
Salmonella Typhi.
3.2.2 MLST—Multi 1. Select MLST under the Typing section on the services page (see
Locus Sequence Typing Note 5) (Fig. 4).
(MLST) from an Assembled 2. Select the MLST configuration database—your pathogen of
Genome or from a Set of interest (From the KMER identification).
Reads
3. Select min. Depth for an allele—default option is suggested.
Analysis of Whole Genome Sequencing Data for Detection of Antimicrobial. . . 215
Fig. 3 (a) The front page of the KmerFinder program when the sequence data has been selected for upload. In
this example, pair-end fastq files have been selected. The progress bar will show the progress of the upload.
(b) The pop-up box once your data has been successfully uploaded. The email is entered, and a link is sent to
that email once the results are ready. (c) The output of the KmerFinder result which shows the top hits of the
submitted sequence data against the KmerFinder database. The identification of this isolate is Salmonella
enterica subspecies enterica. The KmerFinder should not be used for serovar Identification
Analysis of Whole Genome Sequencing Data for Detection of Antimicrobial. . . 217
Fig. 4 (a) The front page of the MLST program when the sequence data has been selected for upload. In this
example, pair-end fastq files have been selected. The progress bar will show the progress of the upload. (b)
The output of the MLST result which shows the MLST allelic profile of the submitted sequence data against the
Salmonella MLST database. In this example, the allelic profile is Sequence Type (ST) 36. This Salmonella ST is
associated as being the serovar Typhimurium
3.3 Identification of 1. Select ResFinder under the Phenotyping section on the ser-
Antimicrobial vices page (see Note 6) (Fig. 5).
Resistance 2. Select for mutations of interest such as “Chromosomal muta-
Determinants tions” or “acquired resistance”.
3.3.1 ResFinder— 3. Select the antimicrobial panel of interest for acquired resis-
Detection of Acquired and tance, default is all.
Point Mutations Associated 4. Select the threshold for %ID and minimum length, default is
with Antimicrobial 90% and 60% respectively.
Resistance
5. Select the appropriate species database—based on the KMER
identification.
6. Select the type of your reads—assemble or sequencing specific
output files.
7. Browse for your files by clicking on the Isolate File Icon.
8. Press upload once you have added your files (see Note 7).
218 Marie Anne Chattaway
Fig. 4 (continued)
9. Enter your email address, the link to your results will be sent
once processed. Remember to check your junk mail.
10. Once the results are ready, click on the link sent to your email.
11. View your data to analyze resistance, you will get a line list of
the resistance markers found, what antibiotic they predict resis-
tance to and where the mutations occur. Unknown mutations
should be validated phenotypically (see Note 8).
12. The algorithm and appropriate reference may vary depending
on your study and the pipeline you have selected (see Note 9).
4 Notes
Fig. 5 (a) The front page of the ResFinder program and the selection of the parameters of the antimicrobial
resistance determinants to be included. (b) The front page of the ResFinder program when the sequence data
has been selected for upload. In this example, pair-end fastq files have been selected. The progress bar will
show the progress of the upload. (c) The output of the ResFinder program: in this example, this isolate
contains aac(6′)-Iaa aminoglycoside enzyme encoding genes which predicts resistance to amikacin and
tobramycin, the qnrS1 plasmid which predicts ciprofloxacin resistance, and the CTX-M-15 gene which
predicts resistance to beta-lactams including third-generation cephalosporins
Fig. 5 (continued)
Analysis of Whole Genome Sequencing Data for Detection of Antimicrobial. . . 221
References
1. Hendriksen RS, Bortolaia V, Tate H, Tyson 5. Davies N, Jorgensen F, Willis C, McLauchlin J,
GH, Aarestrup FM et al (2019) Using geno- Chattaway MA (2022) Whole genome
mics to track global antimicrobial resistance, sequencing reveals antimicrobial resistance
vol 7. Front Public Health Rev, p 17 determinants (AMR genes) of Salmonella
2. AMR NGHRUoGSo (2020) Whole-genome enterica recovered from raw chicken and
sequencing as part of national and international ready-to-eat leaves imported into England
surveillance programmes for antimicrobial between 2014 and 2019. J Appl Microbiol
resistance: a roadmap. BMJ Global Health 133(4):2569–2582
5(11):e002244 6. Baker KS, Dallman TJ, Field N, Childs T,
3. Neuert S, Nair S, Day MR, Doumith M, Ash- Mitchell H et al (2018) Genomic epidemiology
ton PM et al (2018) Prediction of phenotypic of Shigella in the United Kingdom shows
antimicrobial resistance profiles from whole transmission of pathogen sublineages and
genome sequences of non-typhoidal salmo- determinants of antimicrobial resistance. Sci
nella enterica. Front Microbiol 9:592 Rep 8(1):7389
4. Chattaway MA, Dallman TJ, Larkin L, Nair S, 7. Chattaway MA, Gentle A, Nair S, Tingley L,
McCormick J et al (2019) The transformation Day M et al (2021) Phylogenomics and anti-
of reference microbiology methods and surveil- microbial resistance of Salmonella Typhi and
lance for salmonella with the use of whole Paratyphi A, B and C in England,
genome sequencing in England and Wales. 2016–2019. Microbial Genomics 7(8):633
Front Public Health 7:317
Analysis of Whole Genome Sequencing Data for Detection of Antimicrobial. . . 223
Stephen H. Gillespie (ed.), Antibiotic Resistance Protocols, Methods in Molecular Biology, vol. 2833,
https://doi.org/10.1007/978-1-0716-3981-8,
© The Editor(s) (if applicable) and The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer
Nature 2024
225
ANTIBIOTIC RESISTANCE PROTOCOLS
226 Index
G M
Gene essentiality...........................................24, 27, 29–31 Macrophages .................................. 2, 67, 95, 96, 99, 100
Gene expressions .................................................... 29, 110 Malaria ............................................................................. 96
Gene silencing ...........................................................24–33 MALDI-TOF ................................................................ 130
Genome ........................................ 24, 138, 163, 195, 211 Mathematical models ......................................... 79–88, 94
Genome alignment .............................................. 123, 124 Maximum likelihood analysis .............................. 122, 126
Genome assembly ......................................................... 164 Microfluidics........................................................ 112, 113,
Genome-wide association studies (GWAS) ........ 121, 122 115, 134–136
GenTB ........................................................................... 206 Microgravity .................................................................... 12
GenXpert ....................................................................... 138 Microscan ............................................................. 130, 134
Granuloma.................................................................67, 96 Microtiter plate .........................................................65–76
Middlebrook medium.......................................15, 38, 39,
H 67, 68, 74
Heat inactivation ........................................................... 147 Minimum inhibitory concentration
(MIC)........................................................... 4, 5, 8,
Helicobacter pylori .................................................. 96, 206
High containment laboratory (BSL3) ......................... 139 24, 37, 44, 45, 81, 101, 114, 133, 134, 198
HIV-1 .............................................................................. 96 Mixed infections.............................................................. 83
Mobile genetic elements
Hollow fiber .......................................................... v, 57–64
Horizontal gene transfer (HGT) ....................80, 81, 122 (MGEs) ......................................80, 163, 166, 211
Hybrid models ................................................................ 88 Model parameterization.................................................. 88
Model variables ...................................................... 97, 125
Hypervirulent ................................................................ 163
Moxifloxacin ........................................................... 66, 205
I MTBseq ...............................................196–205, 207, 208
Multi-drug resistant organisms (MDROs).................. 161
Immune cells .............................................. 80, 82–86, 94, Multi-drug resistant tuberculosis
95, 99, 100, 102, 103 (MDR-TB)....................................... 185, 187, 196
Individual-based modeling ............................................. 83 Multi-locus sequence typing (MLST) ............... 163, 167,
Infection prevention and control (IPC) ............ 161, 162, 212–215, 217, 222
167, 168 Mutant prevention concentration
Insertion deletions (indels)..........................125, 195–197 (MPC)......................................................... 81, 171
Intensive care unit (ICU) ............................................. 163 Mutation rates ........................................................ 80, 122
Isoniazid (INH) ..............................................13, 66, 112, Mutations ................................................... 44, 66, 80, 81,
114, 116, 117, 138, 185, 205, 206 87, 88, 94, 122, 138, 139, 186, 195–198, 203,
205–208, 211, 215, 218, 222, 223
K
Mycobacterium abscessus ............................. 24, 26, 29–31,
Klebsiella pneumoniae ..................................................... 45 36, 38, 205, 207
Kmer ....................................................212–214, 217, 222 Mycobacterium avium complex ...................................... 96
Mycobacterium bovis .................................. 13–18, 20, 205
L Mycobacterium smegmatis ..................... 24, 112, 114–117
Mycobacterium tuberculosis ..................11, 18, 23, 24, 27,
Laboratory information management systems
38, 65–76, 81, 83, 99, 111, 121, 122, 138, 139,
(LIMS) ............................................................... 162
145–151, 153–155, 159, 185–193, 195, 196,
Larvae .......................................................................... 2–10
204–207
Lesion specific acquired resistance ................................. 86
Mykrobe ........................................................................ 205
Library preparation ............................168, 170, 172–175,
178, 179, 186, 208
N
Ligase chain reaction (LCR) ........................................ 137
Lineage classifications .......................................... 201, 206 Non-destructive method .............................................. 110
Lipid-rich mycobacteria ................................................ 111 Non-tuberculosis mycobacteria
Liposomes........................................................................ 36 (NTM) ............................................. 196, 205, 207
Long read sequencing (LRS) ..................... 163, 164, 166 Norfloxacin (NOR)...............................67, 70, 71, 73, 74
Luria-Bertani (LB) agar/broth ................................26, 38 Nucleic acid amplification test (NAAT)....................... 137
ANTIBIOTIC RESISTANCE PROTOCOLS
Index 227
O Sequence Alignment/Map format (SAM) ......... 124, 207
Sequence type (ST) .............................212, 213, 215, 217
Ordinary differential equations Shigella ........................................................................... 213
(ODEs) .........................84, 85, 87, 88, 93, 94, 96 Short-read data (Illumina)......................... 124, 164, 167,
Organic light emitting diodes (OLEDs) .................52, 54 168, 170, 171, 176, 180, 181, 186, 187, 191,
Outbreaks ............................................161–181, 207, 208 193, 196, 205, 206
Single nucleotide polymorphisms (SNPs) ......... 163, 166,
P
195, 196, 201, 202
Paired end reads ............................................................ 199 Spatial models..................................................... 87, 95, 97
Partial differential equations (PDEs) ..................... 84, 87, SPOTi ..............................................24, 27, 30–31, 35–41
88, 93, 94 Staphylococcus aureus .............................. 1–10, 36, 39, 53,
Personal protective equipment (PPE) ......................... 161 54, 130, 137, 163, 173, 206
Phenotypes ........................................ 12, 13, 18, 99, 103, Staphylococcus epidermidis .........................................53, 54
121–124, 126, 127, 195, 203, 211, 223 Stochastic ............................................................ 87, 88, 94
Phenotypic characterization ........................................... 30 Streptococcus Group A ................................................... 163
Phenotypic knockdown .................................................. 24 Survival analysis .................................................... 121–127
PhiX ............................................................ 170, 173, 176, Survival curves...........................................................4, 6, 7
177, 181, 187, 189, 191 Synergy assessment ......................................................... 73
Photodynamic therapy (PDT)........................... 51, 52, 54
Phylogenetic inference......................................... 123–125 T
Phylogenetic survival analysis ......................123, 126–127 Tagmentation .............................173, 179, 187, 189, 192
PhyResSE....................................................................... 206 Target validation .......................................................23–33
Plasmids ............................................................. 24, 26–30, TB molecular bacterial load assay
80, 163, 164, 211, 219 (TB-MBLA).............................................. 153–159
Polymerase chain reaction (PCR) .................... 24, 26–29, TB Profiler ..................................................................... 205
125, 137, 150, 151, 158, 159, 170, 172–175, T cells ................................................ 95, 96, 99, 100, 103
186, 187, 189, 192, 207 Temporal scales ............................................................... 97
Protospacer adjacent motifs (PAM) .................. 24, 27, 32 Thermocycler ...................................................... 173, 174,
Pyrazinamdidase/nicotinamidase (PncA)...................... 66 186, 187, 189, 192
Pyrazinamide (PZA), potentiation..........................65–76, Transformation................................................... 26, 29, 80
205, 206 Transposons...........................................80, 178, 179, 211
Pyrazinoic acid (POA) .................................................... 66 Typing......................................... 163, 166, 212–215, 217
R V
Raman spectroscopy................................... 110, 111, 113, Valinomycin ..................................................................... 66
115, 134, 135 Viability testing .........................................................35–41
Random/pseudo-random migration ........................... 100 Vibrational fingerprints................................................. 111
Reference mapping .............................164, 166, 196, 197 Vitek 2 ......................................................... 130, 133, 136
Resazurin ..................................................... 68, 72, 74, 75
ResFinder.................. 167, 206, 213, 215, 219, 222, 223 W
Resistance determinants....................................... 137, 196
Resistance mutations............................................ 197, 203 Wavelength-modified Raman spectroscopy................. 134
Resistance protocols...................................................... 137 Whole genome sequencing (WGS)........... 121, 137–140,
Reverse transcriptase ................................... 150, 158, 159 163–169, 172, 173, 177, 181, 185–193,
Rifampicin (RIF) ........................................ 13, 66, 67, 70, 195–208, 211–223
71, 73, 74, 138, 154, 185, 205, 206 Within-host modeling....................................... 79, 81–88,
RNA .................... 24, 26, 29, 33, 62, 147–150, 155–159 95, 103
World Health Organisation (WHO)........... 36, 139, 145,
S 154, 167, 185, 186, 196