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

Deciphering transcript architectural complexity in bacteria and archaea

RNA transcripts are potential therapeutic targets, yet bacterial transcripts have uncharacterized biodiversity. We developed an algorithm for transcript prediction called tp.py using it to predict transcripts (mRNA and other RNAs) in Escherichia coli K12 and E2348/69 strains (Bacteria:gamma-Proteobacteria), Listeria monocytogenes strains Scott A and RO15 (Bacteria:Firmicute), Pseudomonas aeruginosa strains SG17M and NN2 strains (Bacteria:gamma-Proteobacteria), and Haloferax volcanii (Archaea:Halobacteria). From >5 million E. coli K12 and >3 million E. coli E2348/69 newly generated Oxford Nanopore Technologies direct RNA sequencing reads, 2,487 K12 mRNAs and 1,844 E2348/69 mRNAs were predicted, with the K12 mRNAs containing more than half of the predicted E. coli K12 proteins. While the number of predicted transcripts varied by strain based on the amount of sequence data used, across all strains examined, the predicted average size of the mRNAs was 1.6-1.7 kbp, while the median size of the 5'- and 3'-untranslated regions (UTRs) were 30-90 bp. Given the lack of bacterial and archaeal transcript annotation, most predictions were of novel transcripts, but we also predicted many previously characterized mRNAs and ncRNAs, including post-transcriptionally generated transcripts and small RNAs associated with pathogenesis in the E. coli E2348/69 LEE pathogenicity islands. We predicted small transcripts in the 100-200 bp range as well as >10 kbp transcripts for all strains, with the longest transcript for two of the seven strains being the nuo operon transcript, and for another two strains it was a phage/prophage transcript. This quick, easy, and reproducible method will facilitate the presentation of transcripts, and UTR predictions alongside coding sequences and protein predictions in bacterial genome annotation as important resources for the research community.IMPORTANCEOur understanding of bacterial and archaeal genes and genomes is largely focused on proteins since there have only been limited efforts to describe bacterial/archaeal RNA diversity. This contrasts with studies on the human genome, where transcripts were sequenced prior to the release of the human genome over two decades ago. We developed software for the quick, easy, and reproducible prediction of bacterial and archaeal transcripts from Oxford Nanopore Technologies direct RNA sequencing data. These predictions are urgently needed for more accurate studies examining bacterial/archaeal gene regulation, including regulation of virulence factors, and for the development of novel RNA-based therapeutics and diagnostics to combat bacterial pathogens, like those with extreme antimicrobial resistance.

Phenotypic diversity of type III secretion system activity in enteropathogenic Escherichia coli clinical isolates

Introduction. Enteropathogenic Escherichia coli (EPEC) strains pose a significant threat as a leading cause of severe childhood diarrhoea in developing nations. EPEC pathogenicity relies on the type III secretion system (T3SS) encoded by the locus of enterocyte effacement (LEE), facilitating the secretion and translocation of bacterial effector proteins.Gap Statement. While the regulatory roles of PerC (plasmid-encoded regulator) and GrlA (global regulator of LEE-activator) in ler expression and LEE gene activation are well-documented in the EPEC prototype strain E2348/69, understanding the variability in LEE gene expression control mechanisms among clinical EPEC isolates remains an area requiring further investigation.Aim. This study aims to explore the diversity in LEE gene expression control mechanisms among clinical EPEC isolates through a comparative analysis of secretion profiles under defined growth conditions favouring either PerC- or GrlA-mediated activation of LEE expression.Methodology. We compared T3SS-dependent secretion patterns and promoter expression in both typical EPEC (tEPEC) and atypical EPEC (aEPEC) clinical isolates under growth conditions favouring either PerC- or GrlA-mediated activation of LEE expression. Additionally, we conducted promoter reporter activity assays, quantitative real-time PCR and Western blot experiments to assess gene expression activity.Results. Significant differences in T3SS-dependent secretion were observed among tEPEC and aEPEC strains, independent of LEE sequence variations or T3SS gene functionality. Notably, a clinical tEPEC isolate exhibited increased secretion levels under repressive growth conditions and in the absence of both PerC and GrlA, implicating an alternative mechanism in the activation of Ler (LEE-encoded regulator) expression.Conclusion. Our findings indicate that uncharacterized LEE regulatory mechanisms contribute to phenotypic diversity among clinical EPEC isolates, though their impact on clinical outcomes remains unknown. This challenges the conventional understanding based on reference strains and highlights the need to investigate beyond established models to comprehensively elucidate EPEC pathogenesis.

Genomic, transcriptomic, and phenotypic differences among archetype Shigella flexneri strains of serotypes 2a, 3a, and 6

IMPORTANCE

Given the genomic diversity between S. flexneri serotypes and the paucity of data to support serotype-specific phenotypic differences, we applied in silico and in vitro functional analyses of archetype strains of 2457T (Sf2a), J17B (Sf3a), and CH060 (Sf6). These archetype strains represent the three leading S. flexneri serotypes recommended for inclusion in multivalent vaccines. Characterizing the genomic and phenotypic variation among these clinically prevalent serotypes is an important step toward understanding serotype-specific host-pathogen interactions to optimize the efficacy of multivalent vaccines and therapeutics. This study underpins the importance for further large-scale serotype-targeted analyses.

Genomics and pathotypes of the many faces of Escherichia coli

Escherichia coli is the most researched microbial organism in the world. Its varied impact on human health, consisting of commensalism, gastrointestinal disease, or extraintestinal pathologies, has generated a separation of the species into at least eleven pathotypes (also known as pathovars). These are broadly split into two groups, intestinal pathogenic E. coli (InPEC) and extraintestinal pathogenic E. coli (ExPEC). However, components of E. coli's infinite open accessory genome are horizontally transferred with substantial frequency, creating pathogenic hybrid strains that defy a clear pathotype designation. Here, we take a birds-eye view of the E. coli species, characterizing it from historical, clinical, and genetic perspectives. We examine the wide spectrum of human disease caused by E. coli, the genome content of the bacterium, and its propensity to acquire, exchange, and maintain antibiotic resistance genes and virulence traits. Our portrayal of the species also discusses elements that have shaped its overall population structure and summarizes the current state of vaccine development targeted at the most frequent E. coli pathovars. In our conclusions, we advocate streamlining efforts for clinical reporting of ExPEC, and emphasize the pathogenic potential that exists throughout the entire species.

Type III secretion system effector subnetworks elicit distinct host immune responses to infection

Citrobacter rodentium, a natural mouse pathogen which colonises the colon of immuno-competent mice, provides a robust model for interrogating host-pathogen-microbiota interactions in vivo. This model has been key to providing new insights into local host responses to enteric infection, including changes in intestinal epithelial cell immunometabolism and mucosal immunity. C. rodentium injects 31 bacterial effectors into epithelial cells via a type III secretion system (T3SS). Recently, these effectors were shown to be able to form multiple intracellular subnetworks which can withstand significant contractions whilst maintaining virulence. Here we highlight recent advances in understanding gut mucosal responses to infection and effector biology, as well as potential uses for artificial intelligence (AI) in understanding infectious disease and speculate on the role of T3SS effector networks in host adaption.

Intraspecies Transcriptional Profiling Reveals Key Regulators of Candida albicans Pathogenic Traits

The human commensal and opportunistic fungal pathogen Candida albicans displays extensive genetic and phenotypic variation across clinical isolates. Here, we performed RNA sequencing on 21 well-characterized isolates to examine how genetic variation contributes to gene expression differences and to link these differences to phenotypic traits. C. albicans adapts primarily through clonal evolution, and yet hierarchical clustering of gene expression profiles in this set of isolates did not reproduce their phylogenetic relationship. Strikingly, strain-specific gene expression was prevalent in some strain backgrounds. Association of gene expression with phenotypic data by differential analysis, linear correlation, and assembly of gene networks connected both previously characterized and novel genes with 23 C. albicans traits. Construction of de novo gene modules produced a gene atlas incorporating 67% of C. albicans genes and revealed correlations between expression modules and important phenotypes such as systemic virulence. Furthermore, targeted investigation of two modules that have novel roles in growth and filamentation supported our bioinformatic predictions. Together, these studies reveal widespread transcriptional variation across C. albicans isolates and identify genetic and epigenetic links to phenotypic variation based on coexpression network analysis.IMPORTANCE Infectious fungal species are often treated uniformly despite clear evidence of genotypic and phenotypic heterogeneity being widespread across strains. Identifying the genetic basis for this phenotypic diversity is extremely challenging because of the tens or hundreds of thousands of variants that may distinguish two strains. Here, we use transcriptional profiling to determine differences in gene expression that can be linked to phenotypic variation among a set of 21 Candida albicans isolates. Analysis of this transcriptional data set uncovered clear trends in gene expression characteristics for this species and new genes and pathways that were associated with variation in pathogenic processes. Direct investigation confirmed functional predictions for a number of new regulators associated with growth and filamentation, demonstrating the utility of these approaches in linking genes to important phenotypes.

Hybrid Atypical Enteropathogenic and Extraintestinal Escherichia coli (aEPEC/ExPEC) BA1250 Strain: A Draft Genome

Diarrheagenic Escherichia coli is the major bacterial etiological agent of severe diarrhea and a major concern of public health. These pathogens have acquired genetic characteristics from other pathotypes, leading to unusual and singular genetic combinations, known as hybrid strains and may be more virulent due to a set of virulence factors from more than one pathotype. One of the possible combinations is with extraintestinal E. coli (ExPEC), a leading cause of urinary tract infection, often lethal after entering the bloodstream and atypical enteropathogenic E. coli (aEPEC), responsible for death of thousands of people every year, mainly children under five years old. Here we report the draft genome of a strain originally classified as aEPEC (BA1250) isolated from feces of a child with acute diarrhea. Phylogenetic analysis indicates that BA1250 genome content is genetically closer to E. coli strains that cause extraintestinal infections, other than intestinal infections. A deeper analysis showed that in fact this is a hybrid strain, due to the presence of a set of genes typically characteristic of ExPEC. These genomic findings expand our knowledge about aEPEC heterogeneity allowing further studies concerning E. coli pathogenicity and may be a source for future comparative studies, virulence characteristics, and evolutionary biology.

Comparative Transcriptomics of Shiga Toxin-Producing and Commensal Escherichia coli and Cytokine Responses in Colonic Epithelial Cell Culture Infections

Ingestion of Shiga toxin-producing Escherichia coli (STEC) can result in a range of illness severity from asymptomatic to hemorrhagic colitis and death; thus risk assessment of STEC strains for human pathogenicity is important in the area of food safety. Illness severity depends in part on the combination of virulence genes carried in the genome, which can vary between strains even of identical serotype. To better understand how core genes are regulated differently among strains and to identify possible novel STEC virulence gene candidates that could be added to the risk assessment repertoire, we used comparative transcriptomics to investigate global gene expression differences between two STEC strains associated with severe illness and a commensal E. coli strain during in vitro intestinal epithelial cell (IEC) infections. Additionally, we compared a wide array of concomitant cytokine levels produced by the IECs. The cytokine expression levels were examined for a pattern representing STEC pathogenicity; however, while one STEC strain appeared to elicit a proinflammatory response, infection by the other strain produced a pattern comparable to the commensal E. coli. This result may be explained by the significant differences in gene content and expression observed between the STEC strains. RNA-Seq analysis revealed considerable disparity in expression of genes in the arginine and tryptophan biosynthesis/import pathways between the STEC strains and the commensal E. coli strain, highlighting the important role some amino acids play in STEC colonization and survival. Contrasting differential expression patterns were observed for genes involved in respiration among the three strains suggesting that metabolic diversity is a strategy utilized to compete with resident microflora for successful colonization. Similar temporal expression results for known and putative virulence genes were observed in the STEC strains, revealing strategies used for survival prior to and after initial adherence to IECs. Additionally, three genes encoding hypothetical proteins located in mobile genetic elements were, after interrogation of a large set of E. coli genomes, determined to likely represent novel STEC virulence factors.

The population genetics of pathogenic Escherichia coli

Escherichia coli is a commensal of the vertebrate gut that is increasingly involved in various intestinal and extra-intestinal infections as an opportunistic pathogen. Numerous pathotypes that represent groups of strains with specific pathogenic characteristics have been described based on heterogeneous and complex criteria. The democratization of whole-genome sequencing has led to an accumulation of genomic data that render possible a population phylogenomic approach to the emergence of virulence. Few lineages are responsible for the pathologies compared with the diversity of commensal strains. These lineages emerged multiple times during E. coli evolution, mainly by acquiring virulence genes located on mobile elements, but in a specific chromosomal phylogenetic background. This repeated emergence of stable and cosmopolitan lineages argues for an optimization of strain fitness through epistatic interactions between the virulence determinants and the remaining genome.

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.

Pathovar Transcriptomes

Archetypal pathogenic bacterial strains are often used to elucidate regulatory networks of an entire pathovar, which encompasses multiple lineages and phylogroups. With enteropathogenic Escherichia coli (EPEC) as a model system, Hazen and colleagues (mSystems 6:e00024-17, 2017, https://doi.org/10.1128/mSystems.00024-17) used 9 isolates representing 8 lineages and 3 phylogroups to find that isolates with similar genomic sequences exhibit similarities in global transcriptomes under conditions of growth in medium that induces virulence gene expression, and they found variation among individual isolates.

Archetypal pathogenic bacterial strains are often used to elucidate regulatory networks of an entire pathovar, which encompasses multiple lineages and phylogroups. With enteropathogenic Escherichia coli (EPEC) as a model system, Hazen and colleagues (mSystems 6:e00024-17, 2017, https://doi.org/10.1128/mSystems.00024-17) used 9 isolates representing 8 lineages and 3 phylogroups to find that isolates with similar genomic sequences exhibit similarities in global transcriptomes under conditions of growth in medium that induces virulence gene expression. They also found variation among individual isolates. Their work illustrates the importance of moving beyond observing regulatory phenomena of a limited number of regulons in a few archetypal strains, with the possibility of correlating clinical symptoms to key transcriptional pathways across lineages and phylogroups.