Bacterial Factors Associated with Lethal Outcome of Enteropathogenic Escherichia coli Infection: Genomic Case-Control Studies

A phylogeny-aware GWAS framework to correct for heritable pathogen effects on infectious disease traits

Infectious diseases are particularly challenging for genome-wide association studies (GWAS) because genetic effects from two organisms (pathogen and host) can influence a trait. Traditional GWAS assume individual samples are independent observations. However, pathogen effects on a trait can be heritable from donor to recipient in transmission chains. Thus, residuals in GWAS association tests for host genetic effects may not be independent due to shared pathogen ancestry. We propose a new method to estimate and remove heritable pathogen effects on a trait based on the pathogen phylogeny prior to host GWAS, thus restoring independence of samples. In simulations, we show this additional step can increase GWAS power to detect truly associated host variants when pathogen effects are highly heritable, with strong phylogenetic correlations. We applied our framework to data from two different host-pathogen systems, HIV in humans and X. arboricola in A. thaliana. In both systems, the heritability and thus phylogenetic correlations turn out to be low enough such that qualitative results of GWAS do not change when accounting for the pathogen shared ancestry through a correction step. This means that previous GWAS results applied to these two systems should not be biased due to shared pathogen ancestry. In summary, our framework provides additional information on the evolutionary dynamics of traits in pathogen populations and may improve GWAS if pathogen effects are highly phylogenetically correlated amongst individuals in a cohort.

Association between Pathogenic Variants of Diarrheagenic Escherichia coli and Growth in Children under 5 Years of Age in the Global Enteric Multicenter Study

There is a lack of information highlighting associations between different pathogenic variants of diarrheagenic Escherichia coli and childhood growth. Pathogenic variants of E. coli from stool samples, collected from 22,567 children enrolled in the Global Enteric Multicenter Study from December 2007 to March 2011, were detected by real-time polymerase chain reaction. We estimated the associations of different pathogenic variants of diarrheagenic E. coli with child growth. The association between an explanatory variable and the outcome variable was assessed using multiple linear regression, where the dependent variables were height-for-age, weight-for-age, and weight-for-height z-scores, and the independent variable was the presence of different pathogenic variants of diarrheagenic E. coli. After adjusting for potential covariates, such as age, gender, diarrhea, breastfeeding status, mother’s education, number of under-5 children, handwashing practice, handwashing material, source of drinking water, wealth index, available toilet facility, copathogens, comorbidity, time, and study site, the multivariable model identified a negative association between different pathogenic variants of diarrheagenic E. coli and child growth. Our analyses may provide the cornerstone for prospective epidemiologic investigation for the development of preventive measures for diarrheagenic E. coli and combat childhood undernutrition.

Site specific incidence rate of genomic subtypes of enteropathogenic Escherichia coli and association with enteric inflammation and child growth

There is a lack of information highlighting the possible association between the genomic subtypes of enteropathogenic Escherichia coli (EPEC) on environmental enteric dysfunction (EED) and on linear growth during childhood. Genomic subtypes of EPEC from stool samples collected from 1705 children enrolled in the MAL-ED birth cohort were detected by TaqMan Array Cards. We measured site-specific incidence rate by using Poisson regression models, identified the risk factors and estimated the association of genomic subtypes of EPEC with the composite EED score and linear growth at 24 months of age. In general, the highest incidence rate (39%) was found among children having aEPEC infection, which was the greatest in Tanzania (54%). Exclusive breastfeeding and having an improved sanitation facility were found to be protective factors against EPEC infection. In the multivariate models, in overall effect after adjusting for the potential covariates aEPEC showed strong positive associations with the EED scores and tEPEC showed a positive association with poor linear growth at 24 months of age. Our analyses may lay the cornerstone for a prospective epidemiologic investigation for a potential vaccine development aimed at reducing the burden of EPEC infections and combat childhood malnutrition.

Enteropathogenic Escherichia coli Infection in Cancer and Immunosuppressed Patients

BACKGROUND

The role of enteropathogenic Escherichia coli (EPEC) as a cause of diarrhea in cancer and immunocompromised patients is controversial. Quantitation of fecal bacterial loads has been proposed as a method to differentiate colonized from truly infected patients.

METHODS

We studied 77 adult cancer and immunosuppressed patients with diarrhea and EPEC identified in stools by FilmArray, 25 patients with pathogen-negative diarrhea, and 21 healthy adults without diarrhea. Stools were studied by quantitative polymerase chain reaction (qRT-PCR) for EPEC genes eaeA and lifA/efa-1 and strains characterized for virulence factors and adherence to human intestinal enteroids (HIEs).

RESULTS

Patients with EPEC were more likely to have community-acquired diarrhea (odds ratio, 3.82 [95% confidence interval, 1.5-10.0]; P = .008) compared with pathogen-negative cases. Although EPEC was identified in 3 of 21 (14%) healthy subjects by qPCR, the bacterial burden was low compared to patients with diarrhea (≤55 vs median, 6 × 104 bacteria/mg stool; P < .001). Among EPEC patients, the bacterial burden was higher in those who were immunosuppressed (median, 6.7 × 103 vs 55 bacteria/mg; P < .001) and those with fecal lifA/ifa-1 (median, 5 × 104 vs 120 bacteria/mg; P = .015). Response to antimicrobial therapy was seen in 44 of 48 (92%) patients with EPEC as the sole pathogen. Antimicrobial resistance was common and strains exhibited distinct patterns of adherence with variable cytotoxicity when studied in HIEs. Cancer care was delayed in 13% of patients.

CONCLUSIONS

Immunosuppressed cancer patients with EPEC-associated diarrhea carry high burden of EPEC with strains that are resistant to antibiotics, exhibit novel patterns of adherence when studied in HIEs, and interfere with cancer care.

Comparative Genomics of Atypical Enteropathogenic Escherichia coli from Kittens and Children Identifies Bacterial Factors Associated with Virulence in Kittens

Typical enteropathogenic Escherichia coli (tEPEC) is a leading cause of diarrhea and associated death in children worldwide. Atypical EPEC (aEPEC) lacks the plasmid encoding bundle-forming pili and is considered less virulent, but the molecular mechanism of virulence is poorly understood. We recently identified kittens as a host for aEPEC where intestinal epithelial colonization was associated with diarrheal disease and death. The purposes of this study were to (i) determine the genomic similarity between kitten aEPEC and human aEPEC isolates and (ii) identify genotypic or phenotypic traits associated with virulence in kitten aEPEC. We observed no differences between kitten and human aEPEC in core genome content or gene cluster sequence identities, and no distinguishing genomic content was observed between aEPEC isolates from kittens with nonclinical colonization (NC) versus those with lethal infection (LI). Variation in adherence patterns and ability to aggregate actin in cultured cells mirrored descriptions of human aEPEC. The aEPEC isolated from kittens with LI were significantly more motile than isolates from kittens with NC. Kittens may serve as a reservoir for aEPEC that is indistinguishable from human aEPEC isolates and may provide a needed comparative animal model for the study of aEPEC pathogenesis. Motility seems to be an important factor in pathogenesis of LI associated with aEPEC in kittens.

Factors associated with typical enteropathogenic Escherichia coli infection among children <5 years old with moderate-to-severe diarrhoea in rural western Kenya, 2008-2012

Typical enteropathogenic Escherichia coli (tEPEC) infection is a major cause of diarrhoea and contributor to mortality in children <5 years old in developing countries. Data were analysed from the Global Enteric Multicenter Study examining children <5 years old seeking care for moderate-to-severe diarrhoea (MSD) in Kenya. Stool specimens were tested for enteric pathogens, including by multiplex polymerase chain reaction for gene targets of tEPEC. Demographic, clinical and anthropometric data were collected at enrolment and ~60-days later; multivariable logistic regressions were constructed. Of 1778 MSD cases enrolled from 2008 to 2012, 135 (7.6%) children tested positive for tEPEC. In a case-to-case comparison among MSD cases, tEPEC was independently associated with presentation at enrolment with a loss of skin turgor (adjusted odds ratio (aOR) 2.08, 95% confidence interval (CI) 1.37–3.17), and convulsions (aOR 2.83, 95% CI 1.12–7.14). At follow-up, infants with tEPEC compared to those without were associated with being underweight (OR 2.2, 95% CI 1.3–3.6) and wasted (OR 2.5, 95% CI 1.3–4.6). Among MSD cases, tEPEC was associated with mortality (aOR 2.85, 95% CI 1.47–5.55). This study suggests that tEPEC contributes to morbidity and mortality in children. Interventions aimed at defining and reducing the burden of tEPEC and its sequelae should be urgently investigated, prioritised and implemented.

Virulence-Related Genes and Coenteropathogens Associated with Clinical Outcomes of Enteropathogenic Escherichia coli Infections in Children from the Brazilian Semiarid Region: a Case-Control Study of Diarrhea

Enteropathogenic Escherichia coli (EPEC) is a major cause of diarrhea in children from developing countries and presents high genetic variability. We aimed to characterize the EPEC virulence-related gene (VRG) distribution and copathogens associated with diarrhea and nutrition-related outcomes in children from the low-income Brazilian semiarid region. A cross-sectional case-control study of diarrhea was conducted in 1,191 children aged 2 to 36 months from the northeast region of Brazil. Stool samples were collected and clinical, epidemiological, and anthropometric data were identified from each child. A broad molecular evaluation of enteropathogens was performed, and EPEC-positive samples were further investigated for 18 VRGs using five multiplex PCRs. EPEC was detected in 28.2% of the study population, with similar proportions among cases and controls. Typical EPEC (tEPEC) infections were more often associated with diarrhea than atypical EPEC (aEPEC) infections, while aEPEC infections presented a higher prevalence. The VRG ler, a negative regulator of the locus of enterocyte effacement, was associated with the absence of diarrhea in aEPEC-positive children; espB, a major component of the type 3 secretion system, was associated with diarrhea in tEPEC-positive children; the presence of procolonization VRGs-the combination of cesT positivity, espP negativity, and the presence of the map gene-was associated with undernutrition; and Campylobacter spp., norovirus, and enteroaggregative E. coli (EAEC) coinfections were associated with increased clinical severity in EPEC-infected children. These data identified tEPEC strains associated with diarrhea and specific VRGs of EPEC (ler, espB, cesT, and map genes) and Campylobacter spp., norovirus, and EAEC to be major contributors to diarrhea and undernutrition in children from a low-income Brazilian region.

Temporal Variability of Escherichia coli Diversity in the Gastrointestinal Tracts of Tanzanian Children with and without Exposure to Antibiotics

This study increases the number of resident Escherichia coli genome sequences, and explores E. coli diversity through longitudinal sampling. We investigate the genomes of E. coli isolated from human gastrointestinal tracts as part of an antibiotic treatment program among rural Tanzanian children. Phylogenomics demonstrates that resident E. coli are diverse, even within a single host. Though the E. coli isolates of the gastrointestinal community tend to be phylogenomically similar at a given time, they differed across the interrogated time points, demonstrating the variability of the members of the E. coli community in these subjects. Exposure to antibiotic treatment did not have an apparent impact on the E. coli community or the presence of resistance and virulence genes within E. coli genomes. The findings of this study highlight the variable nature of specific bacterial members of the human gastrointestinal tract.

The stability of the Escherichia coli populations in the human gastrointestinal tract is not fully appreciated, and represents a significant knowledge gap regarding gastrointestinal community structure, as well as resistance to incoming pathogenic bacterial species and antibiotic treatment. The current study examines the genomic content of 240 Escherichia coli isolates from 30 children, aged 2 to 35 months old, in Tanzania. The E. coli strains were isolated from three time points spanning a six-month time period, with and without antibiotic treatment. The resulting isolates were sequenced, and the genomes compared. The findings in this study highlight the transient nature of E. coli strains in the gastrointestinal tract of these children, as during a six-month interval, no one individual contained phylogenomically related isolates at all three time points. While the majority of the isolates at any one time point were phylogenomically similar, most individuals did not contain phylogenomically similar isolates at more than two time points. Examination of global genome content, canonical E. coli virulence factors, multilocus sequence type, serotype, and antimicrobial resistance genes identified diversity even among phylogenomically similar strains. There was no apparent increase in the antimicrobial resistance gene content after antibiotic treatment. The examination of the E. coli from longitudinal samples from multiple children in Tanzania provides insight into the genomic diversity and population variability of resident E. coli within the rapidly changing environment of the gastrointestinal tract of these children.

IMPORTANCE This study increases the number of resident Escherichia coli genome sequences, and explores E. coli diversity through longitudinal sampling. We investigate the genomes of E. coli isolated from human gastrointestinal tracts as part of an antibiotic treatment program among rural Tanzanian children. Phylogenomics demonstrates that resident E. coli are diverse, even within a single host. Though the E. coli isolates of the gastrointestinal community tend to be phylogenomically similar at a given time, they differed across the interrogated time points, demonstrating the variability of the members of the E. coli community in these subjects. Exposure to antibiotic treatment did not have an apparent impact on the E. coli community or the presence of resistance and virulence genes within E. coli genomes. The findings of this study highlight the variable nature of specific bacterial members of the human gastrointestinal tract.

Responses of the Human Gut Escherichia coli Population to Pathogen and Antibiotic Disturbances

Research on human-associated E. coli tends to focus on pathogens, such as enterotoxigenic E. coli (ETEC) strains, which are a leading cause of diarrhea in developing countries. However, the severity of disease caused by these pathogens is thought to be influenced by the microbiome. The nonpathogenic E. coli community that resides in the human gastrointestinal tract may play a role in pathogen colonization and disease severity and may become a reservoir for virulence and antibiotic resistance genes. Our study used whole-genome sequencing of E. coli before, during, and after challenge with an archetype ETEC isolate, H10407, and antibiotic treatment to explore the diversity and resiliency of the resident E. coli population in response to the ecological disturbances caused by pathogen invasion and antibiotic treatment.

Studies of Escherichia coli in the human gastrointestinal tract have focused on pathogens, such as diarrhea-causing enterotoxigenic E. coli (ETEC), while overlooking the resident, nonpathogenic E. coli community. Relatively few genomes of nonpathogenic E. coli strains are available for comparative genomic analysis, and the ecology of these strains is poorly understood. This study examined the diversity and dynamics of resident human gastrointestinal E. coli communities in the face of the ecological challenges presented by pathogen (ETEC) challenge, as well as of antibiotic treatment. Whole-genome sequences obtained from E. coli isolates from before, during, and after ETEC challenge were used in phylogenomic and comparative genomic analyses to examine the diversity of the resident E. coli communities, as well as the dynamics of the challenge strain, H10407, a well-studied ETEC strain (serotype O78:H11) that produces both heat-labile and heat-stable enterotoxins. ETEC failed to become the dominant E. coli clone in two of the six challenge subjects, each of whom exhibited limited or no clinical presentation of diarrhea. The E. coli communities of the remaining four subjects became ETEC dominant during the challenge but reverted to their original, subject-specific populations following antibiotic treatment, suggesting resiliency of the resident E. coli population following major ecological disruptions. This resiliency is likely due in part to the abundance of antibiotic-resistant ST131 E. coli strains in the resident populations. This report provides valuable insights into the potential interactions of members of the gastrointestinal microbiome and its responses to challenge by an external pathogen and by antibiotic exposure.

IMPORTANCE Research on human-associated E. coli tends to focus on pathogens, such as enterotoxigenic E. coli (ETEC) strains, which are a leading cause of diarrhea in developing countries. However, the severity of disease caused by these pathogens is thought to be influenced by the microbiome. The nonpathogenic E. coli community that resides in the human gastrointestinal tract may play a role in pathogen colonization and disease severity and may become a reservoir for virulence and antibiotic resistance genes. Our study used whole-genome sequencing of E. coli before, during, and after challenge with an archetype ETEC isolate, H10407, and antibiotic treatment to explore the diversity and resiliency of the resident E. coli population in response to the ecological disturbances caused by pathogen invasion and antibiotic treatment.

Transcriptional Variation of Diverse Enteropathogenic Escherichia coli Isolates under Virulence-Inducing Conditions

Enteropathogenic Escherichia coli (EPEC) bacteria are a diverse group of pathogens that cause moderate to severe diarrhea in young children in developing countries. EPEC isolates can be further subclassified as typical EPEC (tEPEC) isolates that contain the bundle-forming pilus (BFP) or as atypical EPEC (aEPEC) isolates that do not contain BFP. Comparative genomics studies have recently highlighted the considerable genomic diversity among EPEC isolates. In the current study, we used RNA sequencing (RNA-Seq) to characterize the global transcriptomes of eight tEPEC isolates representing the identified genomic diversity, as well as one aEPEC isolate. The global transcriptomes were determined for the EPEC isolates under conditions of laboratory growth that are known to induce expression of virulence-associated genes. The findings demonstrate that unique genes of EPEC isolates from diverse phylogenomic lineages contribute to variation in their global transcriptomes. There were also phylogroup-specific differences in the global transcriptomes, including genes involved in iron acquisition, which had significant differential expression in the EPEC isolates belonging to phylogroup B2. Also, three EPEC isolates from the same phylogenomic lineage (EPEC8) had greater levels of similarity in their genomic content and exhibited greater similarities in their global transcriptomes than EPEC from other lineages; however, even among closely related isolates there were isolate-specific differences among their transcriptomes. These findings highlight the transcriptional variability that correlates with the previously unappreciated genomic diversity of EPEC. IMPORTANCE Recent studies have demonstrated that there is considerable genomic diversity among EPEC isolates; however, it is unknown if this genomic diversity leads to differences in their global transcription. This study used RNA-Seq to compare the global transcriptomes of EPEC isolates from diverse phylogenomic lineages. We demonstrate that there are lineage- and isolate-specific differences in the transcriptomes of genomically diverse EPEC isolates during growth under in vitro virulence-inducing conditions. This study addressed biological variation among isolates of a single pathovar in an effort to demonstrate that while each of these isolates is considered an EPEC isolate, there is significant transcriptional diversity among members of this pathovar. Future studies should consider whether this previously undescribed transcriptional variation may play a significant role in isolate-specific variability of EPEC clinical presentations.

Insights into enterotoxigenic Escherichia coli diversity in Bangladesh utilizing genomic epidemiology

Enterotoxigenic Escherichia coli (ETEC) cause more than 500,000 deaths each year in the developing world and are characterized on a molecular level by the presence of genes that encode the heat-stable (ST) and/or heat-labile (LT) enterotoxins, as well as surface structures, known as colonization factors (CFs). Genome sequencing and comparative genomic analyses of 94 previously uncharacterized ETEC isolates demonstrated remarkable genomic diversity, with 28 distinct sequence types identified in three phylogenomic groups. Interestingly, there is a correlation between the genomic sequence type and virulence factor profiles based on prevalence of the isolate, suggesting that there is an optimal combination of genetic factors required for survival, virulence and transmission in the most successful clones. A large-scale BLAST score ratio (LS-BSR) analysis was further applied to identify ETEC-specific genomic regions when compared to non-ETEC genomes, as well as genes that are more associated with clinical presentations or other genotypic markers. Of the strains examined, 21 of 94 ETEC isolates lacked any previously identified CF. Homology searches with the structural subunits of known CFs identified 6 new putative CF variants. These studies provide a roadmap to exploit genomic analyses by directing investigations of pathogenesis, virulence regulation and vaccine development.

High Prevalence of Virulence Genes in Specific Genotypes of Atypical Enteropathogenic Escherichia coli

Atypical enteropathogenic Escherichia coli (aEPEC) strains are emerging enteropathogens that have been detected worldwide. A collection of 228 aEPEC strains (121 from diarrheal patients, 27 from healthy carriers, 47 from animals and 33 from raw meats) were investigated for serotypes, virulence gene profiles and phylogenetic relationships. Sixty-six O serogroups were identified. Serogroup O51 was the most prevalent, followed by O119, O26 and O76. For the 20 virulence genes detected, statistically significant differences were observed in the overall prevalence of efa1 (lifA), nleB, nleE, set/ent, paa, and ehxA genes among strains from diarrheal patients, healthy carriers, animals and raw meats, respectively. Strains from diarrheal patients had significantly higher levels of efa1 (lifA) (29.8 vs. 0%, P = 0.0002), nleB (41.3 vs. 7.4%, P = 0.0004), nleE (43.8 vs. 7.4%, P = 0.0002) and set/ent (41.3 vs. 7.4%, P = 0.0004) genes than strains obtained from healthy carriers. The paa gene was identified more often in isolates from raw meats (63.6 vs. 14.8%, P < 0.0001), animals (42.6 vs. 14.8%, P < 0.0122), and diarrheal patients (36.4 vs. 14.8%, P < 0.0225) than in strains obtained from healthy carriers. The ehxA gene was detected more frequently in strains from raw meats than in strains from diarrheal patients (27.3 vs. 2.5%, P = 0.0000) and healthy carriers (27.3 vs. 7.4%, P = 0.0474). The phylogenetic marker, yjaA, was more frequently observed in strains among healthy carriers than in diarrheal patient strains. Among the 228 aEPEC strains, 79 sequence types (STs) were identified. The prominent STs, which comprised strains carrying the four OI-122 genes and lpfA, were ST40, ST328, and ST29. Overall, the results indicate that aEPEC strains isolated in China are highly heterogeneous. aEPEC strains that are potentially more pathogenic appear to be related to specific STs or clonal complexes and serotypes. The high prevalence of diarrhea-associated genes in animal or raw meat strains suggests a zoonotic transmission pathway for potentially human pathogenic aEPEC.

Age-Dependent Susceptibility to Enteropathogenic Escherichia coli (EPEC) Infection in Mice

Enteropathogenic Escherichia coli (EPEC) represents a major causative agent of infant diarrhea associated with significant morbidity and mortality in developing countries. Although studied extensively in vitro, the investigation of the host-pathogen interaction in vivo has been hampered by the lack of a suitable small animal model. Using RT-PCR and global transcriptome analysis, high throughput 16S rDNA sequencing as well as immunofluorescence and electron microscopy, we characterize the EPEC-host interaction following oral challenge of newborn mice. Spontaneous colonization of the small intestine and colon of neonate mice that lasted until weaning was observed. Intimate attachment to the epithelial plasma membrane and microcolony formation were visualized only in the presence of a functional bundle forming pili (BFP) and type III secretion system (T3SS). Similarly, a T3SS-dependent EPEC-induced innate immune response, mediated via MyD88, TLR5 and TLR9 led to the induction of a distinct set of genes in infected intestinal epithelial cells. Infection-induced alterations of the microbiota composition remained restricted to the postnatal period. Although EPEC colonized the adult intestine in the absence of a competing microbiota, no microcolonies were observed at the small intestinal epithelium. Here, we introduce the first suitable mouse infection model and describe an age-dependent, virulence factor-dependent attachment of EPEC to enterocytes in vivo.

Bacteria-host relationship: ubiquitin ligases as weapons of invasion

Eukaryotic cells utilize the ubiquitin (Ub) system for maintaining a balanced functioning of cellular pathways. Although the Ub system is exclusive to eukaryotes, prokaryotic bacteria have developed an armory of Ub ligase enzymes that are capable of employing the Ub systems of various hosts, ranging from plant to animal cells. These enzymes have been acquired through the evolution and can be classified into three main classes, RING (really interesting new gene), HECT (homologous to the E6-AP carboxyl terminus) and NEL (novel E3 ligases). In this review we describe the roles played by different classes of bacterial Ub ligases in infection and pathogenicity. We also provide an overview of the different mechanisms by which bacteria mimic specific components of the host Ub system and outline the gaps in our current understanding of their functions. Additionally, we discuss approaches and experimental tools for validating this class of enzymes as potential novel antibacterial therapy targets.

Genomic diversity of EPEC associated with clinical presentations of differing severity

Enteropathogenic Escherichia coli (EPEC) are diarrhoeagenic E. coli, and are a significant cause of gastrointestinal illness among young children in developing countries. Typical EPEC are identified by the presence of the bundle-forming pilus encoded by a virulence plasmid, which has been linked to an increased severity of illness, while atypical EPEC lack this feature. Comparative genomics of 70 total EPEC from lethal (LI), non-lethal symptomatic (NSI) or asymptomatic (AI) cases of diarrhoeal illness in children enrolled in the Global Enteric Multicenter Study was used to investigate the genomic differences in EPEC isolates obtained from individuals with various clinical outcomes. A comparison of the genomes of isolates from different clinical outcomes identified genes that were significantly more prevalent in EPEC isolates of symptomatic and lethal outcomes than in EPEC isolates of asymptomatic outcomes. These EPEC isolates exhibited previously unappreciated phylogenomic diversity and combinations of virulence factors. These comparative results highlight the diversity of the pathogen, as well as the complexity of the EPEC virulence factor repertoire.