A comparison of four epidemic waves of COVID-19 in Malawi; an observational cohort study.

BACKGROUND: Compared to the abundance of clinical and genomic information available on patients hospitalised with COVID-19 disease from high-income countries, there is a paucity of data from low-income countries. Our aim was to explore the relationship between viral lineage and patient outcome. METHODS: We enrolled a prospective observational cohort of adult patients hospitalised with PCR-confirmed COVID-19 disease between July 2020 and March 2022 from Blantyre, Malawi, covering four waves of SARS-CoV-2 infections. Clinical and diagnostic data were collected using an adapted ISARIC clinical characterization protocol for COVID-19. SARS-CoV-2 isolates were sequenced using the MinION™ in Blantyre. RESULTS: We enrolled 314 patients, good quality sequencing data was available for 55 patients. The sequencing data showed that 8 of 11 participants recruited in wave one had B.1 infections, 6/6 in wave two had Beta, 25/26 in wave three had Delta and 11/12 in wave four had Omicron. Patients infected during the Delta and Omicron waves reported fewer underlying chronic conditions and a shorter time to presentation. Significantly fewer patients required oxygen (22.7% [17/75] vs. 58.6% [140/239], p < 0.001) and steroids (38.7% [29/75] vs. 70.3% [167/239], p < 0.001) in the Omicron wave compared with the other waves. Multivariable logistic-regression demonstrated a trend toward increased mortality in the Delta wave (OR 4.99 [95% CI 1.0-25.0 p = 0.05) compared to the first wave of infection. CONCLUSIONS: Our data show that each wave of patients hospitalised with SARS-CoV-2 was infected with a distinct viral variant. The clinical data suggests that patients with severe COVID-19 disease were more likely to die during the Delta wave.

Whole genomic comparative analysis of Streptococcus pneumoniae serotype 1 isolates causing invasive and non-invasive infections among children under 5 years in Casablanca, Morocco.

BACKGROUND: Streptococcus pneumoniae serotype 1 remains a leading cause of invasive pneumococcal diseases, even in countries with PCV-10/PCV-13 vaccine implementation. The main objective of this study, which is part of the Pneumococcal African Genome project (PAGe), was to determine the phylogenetic relationships of serotype 1 isolates recovered from children patients in Casablanca (Morocco), compared to these from other African countries; and to investigate the contribution of accessory genes and recombination events to the genetic diversity of this serotype. RESULTS: The genome average size of the six-pneumococcus serotype 1 from Casablanca was 2,227,119 bp, and the average content of coding sequences was 2113, ranging from 2041 to 2161. Pangenome analysis of the 80 genomes used in this study revealed 1685 core genes and 1805 accessory genes. The phylogenetic tree based on core genes and the hierarchical bayesian clustering analysis revealed five sublineages with a phylogeographic structure by country. The Moroccan strains cluster in two different lineages, the five invasive strains clusters altogether in a divergent clade distantly related to the non-invasive strain, that cluster with all the serotype 1 genomes from Africa. CONCLUSIONS: The whole genome sequencing provides increased resolution analysis of the highly virulent serotype 1 in Casablanca, Morocco. Our results are concordant with previous works, showing that the phylogeography of S. pneumoniae serotype 1 is structured by country, and despite the small size (six isolates) of the Moroccan sample, our analysis shows the genetic cohesion of the Moroccan invasive isolates.

Early Signals of Vaccine-driven Perturbation Seen in Pneumococcal Carriage Population Genomic Data.

BACKGROUND: Pneumococcal conjugate vaccines (PCVs) have reduced pneumococcal diseases globally. Pneumococcal genomic surveys elucidate PCV effects on population structure but are rarely conducted in low-income settings despite the high disease burden. METHODS: We undertook whole-genome sequencing (WGS) of 660 pneumococcal isolates collected through surveys from healthy carriers 2 years from 13-valent PCV (PCV13) introduction and 1 year after rollout in northern Malawi. We investigated changes in population structure, within-lineage serotype dynamics, serotype diversity, and frequency of antibiotic resistance (ABR) and accessory genes. RESULTS: In children <5 years of age, frequency and diversity of vaccine serotypes (VTs) decreased significantly post-PCV, but no significant changes occurred in persons ≥5 years of age. Clearance of VT serotypes was consistent across different genetic backgrounds (lineages). There was an increase of nonvaccine serotypes (NVTs)-namely 7C, 15B/C, and 23A-in children <5 years of age, but 28F increased in both age groups. While carriage rates have been recently shown to remain stable post-PCV due to replacement serotypes, there was no change in diversity of NVTs. Additionally, frequency of intermediate-penicillin-resistant lineages decreased post-PCV. Although frequency of ABR genes remained stable, other accessory genes, especially those associated with mobile genetic element and bacteriocins, showed changes in frequency post-PCV. CONCLUSIONS: We demonstrate evidence of significant population restructuring post-PCV driven by decreasing frequency of vaccine serotypes and increasing frequency of few NVTs mainly in children under 5. Continued surveillance with WGS remains crucial to fully understand dynamics of the residual VTs and replacement NVT serotypes post-PCV.

A mosaic tetracycline resistance gene tet(S/M) detected in an MDR pneumococcal CC230 lineage that underwent capsular switching in South Africa.

OBJECTIVES: We reported tet(S/M) in Streptococcus pneumoniae and investigated its temporal spread in relation to nationwide clinical interventions. METHODS: We whole-genome sequenced 12 254 pneumococcal isolates from 29 countries on an Illumina HiSeq sequencer. Serotype, multilocus ST and antibiotic resistance were inferred from genomes. An SNP tree was built using Gubbins. Temporal spread was reconstructed using a birth-death model. RESULTS: We identified tet(S/M) in 131 pneumococcal isolates and none carried other known tet genes. Tetracycline susceptibility testing results were available for 121 tet(S/M)-positive isolates and all were resistant. A majority (74%) of tet(S/M)-positive isolates were from South Africa and caused invasive diseases among young children (59% HIV positive, where HIV status was available). All but two tet(S/M)-positive isolates belonged to clonal complex (CC) 230. A global phylogeny of CC230 (n=389) revealed that tet(S/M)-positive isolates formed a sublineage predicted to exhibit resistance to penicillin, co-trimoxazole, erythromycin and tetracycline. The birth-death model detected an unrecognized outbreak of this sublineage in South Africa between 2000 and 2004 with expected secondary infections (effective reproductive number, R) of ∼2.5. R declined to ∼1.0 in 2005 and <1.0 in 2012. The declining epidemic could be related to improved access to ART in 2004 and introduction of pneumococcal conjugate vaccine (PCV) in 2009. Capsular switching from vaccine serotype 14 to non-vaccine serotype 23A was observed within the sublineage. CONCLUSIONS: The prevalence of tet(S/M) in pneumococci was low and its dissemination was due to an unrecognized outbreak of CC230 in South Africa. Capsular switching in this MDR sublineage highlighted its potential to continue to cause disease in the post-PCV13 era.

Global emergence and population dynamics of divergent serotype 3 CC180 pneumococci.

Streptococcus pneumoniae serotype 3 remains a significant cause of morbidity and mortality worldwide, despite inclusion in the 13-valent pneumococcal conjugate vaccine (PCV13). Serotype 3 increased in carriage since the implementation of PCV13 in the USA, while invasive disease rates remain unchanged. We investigated the persistence of serotype 3 in carriage and disease, through genomic analyses of a global sample of 301 serotype 3 isolates of the Netherlands3-31 (PMEN31) clone CC180, combined with associated patient data and PCV utilization among countries of isolate collection. We assessed phenotypic variation between dominant clades in capsule charge (zeta potential), capsular polysaccharide shedding, and susceptibility to opsonophagocytic killing, which have previously been associated with carriage duration, invasiveness, and vaccine escape. We identified a recent shift in the CC180 population attributed to a lineage termed Clade II, which was estimated by Bayesian coalescent analysis to have first appeared in 1968 [95% HPD: 1939-1989] and increased in prevalence and effective population size thereafter. Clade II isolates are divergent from the pre-PCV13 serotype 3 population in non-capsular antigenic composition, competence, and antibiotic susceptibility, the last of which resulting from the acquisition of a Tn916-like conjugative transposon. Differences in recombination rates among clades correlated with variations in the ATP-binding subunit of Clp protease, as well as amino acid substitutions in the comCDE operon. Opsonophagocytic killing assays elucidated the low observed efficacy of PCV13 against serotype 3. Variation in PCV13 use among sampled countries was not independently correlated with the CC180 population shift; therefore, genotypic and phenotypic differences in protein antigens and, in particular, antibiotic resistance may have contributed to the increase of Clade II. Our analysis emphasizes the need for routine, representative sampling of isolates from disperse geographic regions, including historically under-sampled areas. We also highlight the value of genomics in resolving antigenic and epidemiological variations within a serotype, which may have implications for future vaccine development.

Global Distribution of Invasive Serotype 35D Streptococcus pneumoniae Isolates following Introduction of 13-Valent Pneumococcal Conjugate Vaccine.

A newly recognized pneumococcal serotype, 35D, which differs from the 35B polysaccharide in structure and serology by not binding to factor serum 35a, was recently reported. The genetic basis for this distinctive serology is due to the presence of an inactivating mutation in wciG, which encodes an O-acetyltransferase responsible for O-acetylation of a galactofuranose. Here, we assessed the genomic data of a worldwide pneumococcal collection to identify serotype 35D isolates and understand their geographical distribution, genetic background, and invasiveness potential. Of 21,980 pneumococcal isolates, 444 were originally typed as serotype 35B by PneumoCaT. Analysis of the wciG gene revealed 23 isolates from carriage (n = 4) and disease (n = 19) with partial or complete loss-of-function mutations, including mutations resulting in premature stop codons (n = 22) and an in-frame mutation (n = 1). These were selected for further analysis. The putative 35D isolates were geographically widespread, and 65.2% (15/23) of them was recovered after the introduction of pneumococcal conjugate vaccine 13 (PCV13). Compared with serotype 35B isolates, putative serotype 35D isolates have higher invasive disease potentials based on odds ratios (OR) (11.58; 95% confidence interval[CI], 1.42 to 94.19 versus 0.61; 95% CI, 0.40 to 0.92) and a higher prevalence of macrolide resistance mediated by mefA (26.1% versus 7.6%; P = 0.009). Using the Quellung reaction, 50% (10/20) of viable isolates were identified as serotype 35D, 25% (5/20) as serotype 35B, and 25% (5/20) as a mixture of 35B/35D. The discrepancy between phenotype and genotype requires further investigation. These findings illustrated a global distribution of an invasive serotype, 35D, among young children post-PCV13 introduction and underlined the invasive potential conferred by the loss of O-acetylation in the pneumococcal capsule.

Comparative Genomic Analysis and In Vivo Modeling of Streptococcus pneumoniae ST3081 and ST618 Isolates Reveal Key Genetic and Phenotypic Differences Contributing to Clonal Replacement of Serotype 1 in The Gambia.

Streptococcus pneumoniae serotype 1 is one of the leading causes of invasive pneumococcal disease (IPD) in West Africa, with ST618 being the dominant cause of IPD in The Gambia. Recently however, a rare example of clonal replacement was observed, where the ST3081 clone of serotype 1 replaced the predominant ST618 clone as the main cause of IPD. In the current study, we sought to find the reasons for this unusual replacement event. Using whole-genome sequence analysis and clinically relevant models of in vivo infection, we identified distinct genetic and phenotypic characteristics of the emerging ST3081 clone. We show that ST3081 is significantly more virulent than ST618 in models of invasive pneumonia, and is carried at higher densities than ST618 during nasopharyngeal carriage. We also observe sequence type-specific accessory genes and a unique sequence type-specific fixed mutation in the pneumococcal toxin pneumolysin, which is associated with increased hemolytic activity in ST3081 and may contribute to increased virulence in this clone. Our study provides evidence that, within the same serotype 1 clonal complex, biological properties differ significantly from one clone to another in terms of virulence and host invasiveness, and that these differences may be the result of key genetic differences within the genome.

Epidemiological and Molecular Characterization of an Invasive Group A Streptococcus emm32.2 Outbreak.

An emm32.2 invasive group A streptococcus (iGAS) outbreak occurred in Liverpool from January 2010 to September 2012. This genotype had not previously been identified in Liverpool, but was responsible for 32% (14/44) of all iGAS cases reported during this time period. We performed a case-case comparison of emm32.2 iGAS cases with non-emm32.2 control iGAS cases identified in the Liverpool population over the same time period to assess patient risk factors for emm32.2 iGAS infection. The emm32.2 iGAS cases were confined to the adult population. We show that homelessness, intravenous drug use, and alcohol abuse predisposed patients to emm32.2 iGAS disease; however, no obvious epidemiological linkage between the patients with emm32.2 iGAS could be identified. Comparative whole-genome sequencing analysis of emm32.2 iGAS and non-emm32.2 control isolates was also performed to identify pathogen factors which might have driven the outbreak. We identified 19 genes, five of which had previously been implicated in virulence, which were present in all of the emm32.2 iGAS isolates but not present in any of the non-emm32.2 control isolates. We report that a novel emm32.2 genotype emerged in Liverpool in 2010 and identified a specific subset of genes, which could have allowed this novel emm32.2 genotype to persist in a disadvantaged population in the region over a 3-year period.

Phylogenetic Analysis of Invasive Serotype 1 Pneumococcus in South Africa, 1989 to 2013.

Serotype 1 is an important cause of invasive pneumococcal disease in South Africa and has declined following the introduction of the 13-valent pneumococcal conjugate vaccine in 2011. We genetically characterized 912 invasive serotype 1 isolates from 1989 to 2013. Simpson's diversity index (D) and recombination ratios were calculated. Factors associated with sequence types (STs) were assessed. Clonal complex 217 represented 96% (872/912) of the sampled isolates. Following the introduction of the 13-valent pneumococcal conjugate vaccine (PCV13), ST diversity increased in children <5 years (D, 0.39 to 0.63, P = 0.002) and individuals >14 years (D, 0.35 to 0.54, P < 0.001): ST-217 declined proportionately in children <5 years (153/203 [75%] versus 21/37 [57%], P = 0.027) and individuals >14 years (242/305 [79%] versus 96/148 [65%], P = 0.001), whereas ST-9067 increased (4/684 [0.6%] versus 24/228 [11%], P < 0.001). Three subclades were identified within ST-217: ST-217C1 (353/382 [92%]), ST-217C2 (15/382 [4%]), and ST-217C3 (14/382 [4%]). ST-217C2, ST-217C3, and single-locus variant (SLV) ST-8314 (20/912 [2%]) were associated with nonsusceptibility to chloramphenicol, tetracycline, and co-trimoxazole. ST-8314 (20/912 [2%]) was also associated with increased nonsusceptibility to penicillin (P < 0.001). ST-217C3 and newly reported ST-9067 had higher recombination ratios than those of ST-217C1 (4.344 versus 0.091, P < 0.001; and 0.086 versus 0.013, P < 0.001, respectively). Increases in genetic diversity were noted post-PCV13, and lineages associated with antimicrobial nonsusceptibility were identified.

Understanding pneumococcal serotype 1 biology through population genomic analysis.

BACKGROUND: Pneumococcus kills over one million children annually and over 90 % of these deaths occur in low-income countries especially in Sub-Saharan Africa (SSA) where HIV exacerbates the disease burden. In SSA, serotype 1 pneumococci particularly the endemic ST217 clone, causes majority of the pneumococcal disease burden. To understand the evolution of the virulent ST217 clone, we analysed ST217 whole genomes from isolates sampled from African and Asian countries. METHODS: We analysed 226 whole genome sequences from the ST217 lineage sampled from 9 African and 4 Asian countries. We constructed a whole genome alignment and used it for phylogenetic and coalescent analyses. We also screened the genomes to determine presence of antibiotic resistance conferring genes. RESULTS: Population structure analysis grouped the ST217 isolates into five sequence clusters (SCs), which were highly associated with different geographical regions and showed limited intracontinental and intercontinental spread. The SCs showed lower than expected genomic sequence, which suggested strong purifying selection and small population sizes caused by bottlenecks. Recombination rates varied between the SCs but were lower than in other successful clones such as PMEN1. African isolates showed higher prevalence of antibiotic resistance genes than Asian isolates. Interestingly, certain West African isolates harbored a defective chloramphenicol and tetracycline resistance-conferring element (Tn5253) with a deletion in the loci encoding the chloramphenicol resistance gene (cat pC194), which caused lower chloramphenicol than tetracycline resistance. Furthermore, certain genes that promote colonisation were absent in the isolates, which may contribute to serotype 1's rarity in carriage and consequently its lower recombination rates. CONCLUSIONS: The high phylogeographic diversity of the ST217 clone shows that this clone has been in circulation globally for a long time, which allowed its diversification and adaptation in different geographical regions. Such geographic adaptation reflects local variations in selection pressures in different locales. Further studies will be required to fully understand the biological mechanisms which makes the ST217 clone highly invasive but unable to successfully colonise the human nasopharynx for long durations which results in lower recombination rates.

Comparative Genomic Analysis of Meningitis- and Bacteremia-Causing Pneumococci Identifies a Common Core Genome.

Streptococcus pneumoniae is a nasopharyngeal commensal that occasionally invades normally sterile sites to cause bloodstream infection and meningitis. Although the pneumococcal population structure and evolutionary genetics are well defined, it is not clear whether pneumococci that cause meningitis are genetically distinct from those that do not. Here, we used whole-genome sequencing of 140 isolates of S. pneumoniae recovered from bloodstream infection (n = 70) and meningitis (n = 70) to compare their genetic contents. By fitting a double-exponential decaying-function model, we show that these isolates share a core of 1,427 genes (95% confidence interval [CI], 1,425 to 1,435 genes) and that there is no difference in the core genome or accessory gene content from these disease manifestations. Gene presence/absence alone therefore does not explain the virulence behavior of pneumococci that reach the meninges. Our analysis, however, supports the requirement of a range of previously described virulence factors and vaccine candidates for both meningitis- and bacteremia-causing pneumococci. This high-resolution view suggests that, despite considerable competency for genetic exchange, all pneumococci are under considerable pressure to retain key components advantageous for colonization and transmission and that these components are essential for access to and survival in sterile sites.

High multiple carriage and emergence of Streptococcus pneumoniae vaccine serotype variants in Malawian children.

BACKGROUND: Carriage of either single or multiple pneumococcal serotypes (multiple carriage) is a prerequisite for developing invasive pneumococcal disease. However, despite the reported high rates of pneumococcal carriage in Malawi, no data on carriage of multiple serotypes has been reported previously. Our study provides the first description of the prevalence of multiple pneumococcal carriage in Malawi. METHODS: The study was conducted in Blantyre and Karonga districts in Malawi, from 2008 to 2012. We recruited 116 children aged 0-13 years. These children were either HIV-infected (N = 44) or uninfected (N = 72). Nasopharyngeal samples were collected using sterile swabs. Pneumococcal serotypes in the samples were identified by microarray. Strains that could not be typed by microarray were sequenced to characterise possible genetic alterations within the capsular polysaccharide (CPS) locus. RESULTS: The microarray identified 179 pneumococcal strains (from 116 subjects), encompassing 43 distinct serotypes and non-typeable (NT) strains. Forty per cent (46/116) of children carried multiple serotypes. Carriage of vaccine type (VT) strains was higher (p = 0.028) in younger (0-2 years) children (71 %, 40/56) compared to older (3-13 years) children (50 %, 30/60). Genetic variations within the CPS locus of known serotypes were observed in 19 % (34/179) of the strains identified. The variants included 13-valent pneumococcal conjugate vaccine (PCV13) serotypes 6B and 19A, and the polysaccharide vaccine serotype 20. Serotype 6B variants were the most frequently isolated (47 %, 16/34). Unlike the wild type, the CPS locus of the 6B variants contained an insertion of the licD-family phosphotransferase gene. The CPS locus of 19A- and 20-variants contained an inversion in the sugar-biosynthesis (rmlD) gene and a 717 bp deletion within the transferase (whaF) gene, respectively. CONCLUSIONS: The high multiple carriage in Malawian children provides opportunities for genetic exchange through horizontal gene transfer. This may potentially lead to CPS locus variants and vaccine escape. Variants reported here occurred naturally, however, PCV13 introduction could exacerbate the CPS genetic variations. Further studies are therefore recommended to assess the invasive potential of these variants and establish whether PCV13 would offer cross-protection. We have shown that younger children (0-2 years) are a reservoir of VT serotypes, which makes them an ideal target for vaccination.

Genomic identification of a novel co-trimoxazole resistance genotype and its prevalence amongst Streptococcus pneumoniae in Malawi.

OBJECTIVES: This study aimed to define the molecular basis of co-trimoxazole resistance in Malawian pneumococci under the dual selective pressure of widespread co-trimoxazole and sulfadoxine/pyrimethamine use. METHODS: We measured the trimethoprim and sulfamethoxazole MICs and analysed folA and folP nucleotide and translated amino acid sequences for 143 pneumococci isolated from carriage and invasive disease in Malawi (2002-08). RESULTS: Pneumococci were highly resistant to both trimethoprim and sulfamethoxazole (96%, 137/143). Sulfamethoxazole-resistant isolates showed a 3 or 6 bp insertion in the sulphonamide-binding site of folP. The trimethoprim-resistant isolates fell into three genotypic groups based on dihydrofolate reductase (encoded by folA) mutations: Ile-100-Leu (10%), the Ile-100-Leu substitution together with a residue 92 substitution (56%) and those with a novel uncharacterized resistance genotype (34%). The nucleotide sequence divergence and dN/dS of folA and folP remained stable from 2004 onwards. CONCLUSIONS: S. pneumoniae exhibit almost universal co-trimoxazole resistance in vitro and in silico that we believe is driven by extensive co-trimoxazole and sulfadoxine/pyrimethamine use. More than one-third of pneumococci employ a novel mechanism of co-trimoxazole resistance. Resistance has now reached a point of stabilizing evolution. The use of co-trimoxazole to prevent pneumococcal infection in HIV/AIDS patients in sub-Saharan Africa should be re-evaluated.

Longitudinal analysis within one hospital in sub-Saharan Africa over 20 years reveals repeated replacements of dominant clones of Klebsiella pneumoniae and stresses the importance to include temporal patterns for vaccine design considerations.

BACKGROUND: Infections caused by multidrug-resistant gram-negative bacteria present a severe threat to global public health. The WHO defines drug-resistant Klebsiella pneumoniae as a priority pathogen for which alternative treatments are needed given the limited treatment options and the rapid acquisition of novel resistance mechanisms by this species. Longitudinal descriptions of genomic epidemiology of Klebsiella pneumoniae can inform management strategies but data from sub-Saharan Africa are lacking. METHODS: We present a longitudinal analysis of all invasive K. pneumoniae isolates from a single hospital in Blantyre, Malawi, southern Africa, from 1998 to 2020, combining clinical data with genome sequence analysis of the isolates. RESULTS: We show that after a dramatic increase in the number of infections from 2016 K. pneumoniae becomes hyperendemic, driven by an increase in neonatal infections. Genomic data show repeated waves of clonal expansion of different, often ward-restricted, lineages, suggestive of hospital-associated transmission. We describe temporal trends in resistance and surface antigens, of relevance for vaccine development. CONCLUSIONS: Our data highlight a clear need for new interventions to prevent rather than treat K. pneumoniae infections in our setting. Whilst one option may be a vaccine, the majority of cases could be avoided by an increased focus on and investment in infection prevention and control measures, which would reduce all healthcare-associated infections and not just one.

The Enterics for Global Health (EFGH) Shigella Surveillance Study in Malawi.

Background: Malawi is among 7 countries participating in the Enterics for Global Health (EFGH) Shigella surveillance study, which aims to determine the incidence of medically attended diarrhea attributed to Shigella, a leading bacterial cause of diarrhea in children in low-resource settings. Methods: We describe the EFGH study site in the densely populated informal settlement of Ndirande Township, Blantyre, Malawi. We explore the site's geographical location, demographic characteristics, and the healthcare-seeking behavior of its population, particularly for childhood diarrhea. We also describe the management of childhood diarrhea at the health facility, and the associated challenges to attaining optimum adherence to local and national guidelines at the site. Conclusions: Our overarching aim is to improve global health through understanding and mitigating the impact of diarrhea attributed to Shigella.

Data Management in Multicountry Consortium Studies: The Enterics For Global Health (EFGH) Shigella Surveillance Study Example.

Background: Rigorous data management systems and planning are essential to successful research projects, especially for large, multicountry consortium studies involving partnerships across multiple institutions. Here we describe the development and implementation of data management systems and procedures for the Enterics For Global Health (EFGH) Shigella surveillance study-a 7-country diarrhea surveillance study that will conduct facility-based surveillance concurrent with population-based enumeration and a health care utilization survey to estimate the incidence of Shigella--associated diarrhea in children 6 to 35 months old. Methods: The goals of EFGH data management are to utilize the knowledge and experience of consortium members to collect high-quality data and ensure equity in access and decision-making. During the planning phase before study initiation, a working group of representatives from each EFGH country site, the coordination team, and other partners met regularly to develop the data management systems for the study. Results: This resulted in the Data Management Plan, which included selecting REDCap and SurveyCTO as the primary database systems. Consequently, we laid out procedures for data processing and storage, study monitoring and reporting, data quality control and assurance activities, and data access. The data management system and associated real-time visualizations allow for rapid data cleaning activities and progress monitoring and will enable quicker time to analysis. Conclusions: Experiences from this study will contribute toward enriching the sparse landscape of data management methods publications and serve as a case study for future studies seeking to collect and manage data consistently and rigorously while maintaining equitable access to and control of data.

Diarrhea Case Surveillance in the Enterics for Global Health Shigella Surveillance Study: Epidemiologic Methods.

Background: Shigella is a leading cause of acute watery diarrhea, dysentery, and diarrhea-attributed linear growth faltering, a precursor to stunting and lifelong morbidity. Several promising Shigella vaccines are in development and field efficacy trials will require a consortium of potential vaccine trial sites with up-to-date Shigella diarrhea incidence data. Methods: The Enterics for Global Health (EFGH) Shigella surveillance study will employ facility-based enrollment of diarrhea cases aged 6-35 months with 3 months of follow-up to establish incidence rates and document clinical, anthropometric, and financial consequences of Shigella diarrhea at 7 country sites (Mali, Kenya, The Gambia, Malawi, Bangladesh, Pakistan, and Peru). Over a 24-month period between 2022 and 2024, the EFGH study aims to enroll 9800 children (1400 per country site) between 6 and 35 months of age who present to local health facilities with diarrhea. Shigella species (spp.) will be identified and serotyped from rectal swabs by conventional microbiologic methods and quantitative polymerase chain reaction. Shigella spp. isolates will undergo serotyping and antimicrobial susceptibility testing. Incorporating population and healthcare utilization estimates from contemporaneous household sampling in the catchment areas of enrollment facilities, we will estimate Shigella diarrhea incidence rates. Conclusions: This multicountry surveillance network will provide key incidence data needed to design Shigella vaccine trials and strengthen readiness for potential trial implementation. Data collected in EFGH will inform policy makers about the relative importance of this vaccine-preventable disease, accelerating the time to vaccine availability and uptake among children in high-burden settings.

Exploring Natural Immune Responses to Shigella Exposure Using Multiplex Bead Assays on Dried Blood Spots in High-Burden Countries: Protocol From a Multisite Diarrhea Surveillance Study.

Background: Molecular diagnostics on human fecal samples have identified a larger burden of shigellosis than previously appreciated by culture. Evidence of fold changes in immunoglobulin G (IgG) to conserved and type-specific Shigella antigens could be used to validate the molecular assignment of type-specific Shigella as the etiology of acute diarrhea and support polymerase chain reaction (PCR)-based microbiologic end points for vaccine trials. Methods: We will test dried blood spots collected at enrollment and 4 weeks later using bead-based immunoassays for IgG to invasion plasmid antigen B and type-specific lipopolysaccharide O-antigen for Shigella flexneri 1b, 2a, 3a, and 6 and Shigella sonnei in Shigella-positive cases and age-, site-, and season-matched test-negative controls from all sites in the Enterics for Global Health (EFGH) Shigella surveillance study. Fold antibody responses will be compared between culture-positive, culture-negative but PCR-attributable, and PCR-positive but not attributable cases and test-negative controls. Age- and site-specific seroprevalence distributions will be identified, and the association between baseline antibodies and Shigella attribution will be estimated. Conclusions: The integration of these assays into the EFGH study will help support PCR-based attribution of acute diarrhea to type-specific Shigella, describe the baseline seroprevalence of conserved and type-specific Shigella antibodies, and support correlates of protection for immunity to Shigella diarrhea. These insights can help support the development and evaluation of Shigella vaccine candidates.

Shigella Detection and Molecular Serotyping With a Customized TaqMan Array Card in the Enterics for Global Health (EFGH): Shigella Surveillance Study.

Background: Quantitative polymerase chain reaction (qPCR) targeting ipaH has been proven to be highly efficient in detecting Shigella in clinical samples compared to culture-based methods, which underestimate Shigella burden by 2- to 3-fold. qPCR assays have also been developed for Shigella speciation and serotyping, which is critical for both vaccine development and evaluation. Methods: The Enterics for Global Health (EFGH) Shigella surveillance study will utilize a customized real-time PCR-based TaqMan Array Card (TAC) interrogating 82 targets, for the detection and differentiation of Shigella spp, Shigella sonnei, Shigella flexneri serotypes, other diarrhea-associated enteropathogens, and antimicrobial resistance (AMR) genes. Total nucleic acid will be extracted from rectal swabs or stool samples, and assayed on TAC. Quantitative analysis will be performed to determine the likely attribution of Shigella and other particular etiologies of diarrhea using the quantification cycle cutoffs derived from previous studies. The qPCR results will be compared to conventional culture, serotyping, and phenotypic susceptibility approaches in EFGH. Conclusions: TAC enables simultaneous detection of diarrheal etiologies, the principal pathogen subtypes, and AMR genes. The high sensitivity of the assay enables more accurate estimation of Shigella-attributed disease burden, which is critical to informing policy and in the design of future clinical trials.

Microbiological Methods Used in the Enterics for Global Health Shigella Surveillance Study.

Background: Shigella is a major cause of diarrhea in young children worldwide. Multiple vaccines targeting Shigella are in development, and phase 3 clinical trials are imminent to determine efficacy against shigellosis. Methods: The Enterics for Global Health (EFGH) Shigella surveillance study is designed to determine the incidence of medically attended shigellosis in 6- to 35-month-old children in 7 resource-limited settings. Here, we describe the microbiological methods used to isolate and identify Shigella. We developed a standardized laboratory protocol for isolation and identification of Shigella by culture. This protocol was implemented across all 7 sites, ensuring consistency and comparability of results. Secondary objectives of the study are to determine the antibiotic resistance profiles of Shigella, compare isolation of Shigella from rectal swabs versus whole stool, and compare isolation of Shigella following transport of rectal swabs in Cary-Blair versus a modified buffered glycerol saline transport medium. Conclusions: Data generated from EFGH using culture methods described herein can potentially be used for microbiological endpoints in future phase 3 clinical trials to evaluate vaccines against shigellosis and for other clinical and public health studies focused on these organisms.