# Prevalence of Enteropathogens and Virulence Traits in Brazilian Children With and Without Diarrhea

Victor R Merino, Viviane Nakano, Sabine Delannoy, Patrick Fach, Gabriela G F Alberca, Mauricio J Farfan, Roxane M F Piazza, Mario J Avila-Campos

► To cite this version:

Victor R Merino, Viviane Nakano, Sabine Delannoy, Patrick Fach, Gabriela G F Alberca, et al. Prevalence of Enteropathogens and Virulence Traits in Brazilian Children With and Without Diarrhea. Frontiers in Cellular and Infection Microbiology, 2020, 10, 10.3389/fcimb.2020.549919. anses-03715607

HAL Id: anses-03715607

https://anses.hal.science/anses-03715607

Submitted on 6 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Prevalence of Enteropathogens and Virulence Traits in Brazilian Children With and Without Diarrhea

OPEN ACCESS

**Edited by:**
Jorge Eugenio Vidal,
University of Mississippi Medical
Center, United States

**Reviewed by:**
Javier Sanchez-Villamil,
The University of Texas Medical
Branch at Galveston, United States
Miguel Angel Garcia Bereguiain,
University of the Americas, Ecuador

*Correspondence:
Roxane M. F. Piazza
roxane.piazza@butantan.gov.br
Mauricio J. Farfan
mfarfan@med.uchile.cl

In Memoriam:
This work is dedicated to Prof. Dr.
Mario J. Avila-Campos who conceived
and supervised this project

†These authors have contributed
equally to this work

‡Present address:
Viviane Nakano,
GeneOne, DASA, São Paulo, Brazil

**Specialty section:**
This article was submitted to
Clinical Microbiology,
a section of the journal
Frontiers in Cellular and Infection
Microbiology

**Received:** 07 April 2020
**Accepted:** 25 August 2020
**Published:** 25 September 2020

**Citation:**
Merino VR, Nakano V, Delannoy S,
Fach P, Alberca GGF, Farfan MJ,
Piazza RMF and Avila-Campos MJ
(2020) Prevalence of Enteropathogens
and Virulence Traits in Brazilian
Children With and Without Diarrhea.
Front. Cell. Infect. Microbiol. 10:549919.
doi: 10.3389/fcimb.2020.549919

Victor R. Merino¹†, Viviane Nakano ¹†‡, Sabine Delannoy², Patrick Fach², Gabriela G. F. Alberca¹, Mauricio J. Farfan³*, Roxane M. F. Piazza⁴*and Mario J. Avila-Campos¹

¹Laboratório de Anaeróbios, Departamento de Microbiologia, Instituto de Ciências Biomédicas, Universidade de São Paulo, São Paulo, Brazil, ²Agence Nationale de Sécurité Sanitaire de l'Alimentation, de l'Environnement et du Travail, Food Safety Laboratory, Université Paris-Est, Champs-sur-Marne, France, ³Departamento de Pediatria y Cirugia Infantil, Facultad de Medicina, Campus Oriente-Hospital Dr. Luis Calvo Mackenna, Universidad de Chile, Santiago, Chile, ⁴Laboratório de Bacteriologia, Instituto Butantan, São Paulo, Brazil

The use of molecular diagnostics for pathogen detection in epidemiological studies have allowed us to get a wider view of the pathogens associated with diarrhea, but the presence of enteropathogens in asymptomatic individuals has raised several challenges in understanding the etiology of diarrhea, and the use of these platforms in clinical diagnosis as well. To characterize the presence of the most relevant bacterial enteropathogens in diarrheal episodes, we evaluated here the prevalence of diarrheagenic *E. coli* pathotypes, *Salmonella* spp., and *Yersinia enterocolitica* in stool samples of children with and without diarrhea using real-time quantitative PCR (qPCR). We found that the presence of genetic markers associated with bacterial pathogens was significantly higher in stool samples from the diarrhea group compared to the control (*P* < 0.001). Bacterial loads in samples positive for eae and *aggR* markers were also determined. Compared to samples from asymptomatic children, a significantly higher number of copies of the eae gene were found in diarrhea samples. Also, the presence of genetic markers associated with STEC strains with clinical significance was evaluated in eae-positive samples by high-throughput real-time PCR. The data presented herein demonstrated that asymptomatic children of an urban area in Brazil might be enteropathogen reservoirs, especially for STEC.

**Keywords:** enteric pathogens, *Escherichia coli* pathotypes, diarrheal children, healthy children, molecular diagnostics

## INTRODUCTION

Diarrheal disease is the second most important cause of child morbidity and mortality in the world, and is frequently due to contaminated food and water. Worldwide, 780 million individuals lack access to safe drinking water and 2.5 billion lack adequate sanitation. Diarrhea is widespread throughout developing countries. In low-income countries, children under 3 years old experience on average three episodes of diarrhea every year (World Health Organization, 2017).
Rotavirus and diarrheagenic *Escherichia coli* (DEC) are the two most common etiological agents of moderate-to-severe diarrhea in low-income countries (World Health Organization, 2017). Other pathogens such as Cryptosporidium, Shigella, *Campylobacter*, *Salmonella*, and *Yersinia* species also cause human bacterial gastroenteritis, but location-specific etiological patterns also need to be considered (Bottone, 1999; Freitas Neto et al., 2010; Marinou et al., 2012; Nunes et al., 2012; World Health Organization, 2017; Lima F. M. et al., 2019; Leli et al., 2020).

Enterotoxigenic *E. coli* (ETEC), enteroinvasive *E. coli* (EIEC), typical and atypical enteropathogenic *E. coli* (tEPEC and aEPEC), typical and atypical enteroaggregative *E. coli* (tEAEC and aEAEC), and Shiga toxin-producing *E. coli* (STEC) and its subgroup enterohemorrhagic *E. coli* (EHEC), enteroinvasive *E. coli* (EIEC), and diffusely adherent *E. coli* (DAEC) comprise DEC pathotypes (Kaper et al., 2004; Gomes et al., 2016). These pathogens use many common features to colonize the intestinal mucosa and trigger disease, but its initiation, progression, and complications vary considerably (Croxen et al., 2013). Thus, ETEC isolates carry the enterotoxins LT and/or ST on plasmids, as well as colonization factors (CFs). tEPEC, aEPEC, and LEE-positive STEC strains (EHEC) carry the locus of enterocyte effacement (LEE) which contains *eae*, the gene that encodes intimin, but only tEPEC has the bundle-forming pilus gene (bfp). STEC strains, LEE-positive and -negative, carry the Shiga toxin genes (stx1 and/or stx2). Some virulence factors for EAEC for a certain group of isolates are found on the pAA plasmid, and the presence of the *aggR* gene differentiates tEAEC from aEAEC. EIEC and also *Shigella spp.* have the ability to invade cells mainly through the pINV plasmid and acquired additional virulence attributes mainly from chromosomal pathogenicity islands (PAIs) (Croxen et al., 2013).

Nevertheless, the designation of certain isolates into a pathotype is complicated due to the plasticity of the *E. coli* genome resulting in the emergence of hybrid strains; besides, DEC epidemiology varies in different regions of the world and may depend on the host and environmental factors (Croxen et al., 2013). It is important to note that DEC pathotypes have been largely disregarded because of the lack of laboratory capabilities (Piazza et al., 2010; Miliwebsky et al., 2016). Moreover, in several case-control studies, it shows almost equal distribution, especially in EPEC and EAEC pathotypes (Bueris et al., 2007; Hernandes et al., 2009; Gomes et al., 2016; Imdad et al., 2018; Dias et al., 2020).

Besides LEE-positive pathogens and *Shigella*, *Salmonella* also uses the type III secretion system (T3SS) apparatus to deliver effector proteins into the host cell, and *invA*, an inner-membrane component, is critical to the functioning of *Salmonella* T3SS (Galán et al., 1992; Ginocchio and Galán, 1995). In addition, the most common virulence-associated gene in pathogenic Yersinia enterocolitica proved to be *ystA*, which can therefore be considered the best target gene to be amplified to evaluate the presence of pathogenic biotypes (Peruzy et al., 2017).

The use of multiple pathogen detection platforms in epidemiological studies have allowed us to get a broader view of the pathogens associated with diarrhea, but the presence of enteropathogens in asymptomatic individuals has raised

several challenges in understanding the etiology of diarrhea, as well as the use of these platforms in clinical diagnosis (Liu et al., 2012; Lima A. A. M et al., 2019). The use of real-time quantitative PCR (qPCR) to determine bacterial loads in stool samples appears as an alternative to associate a pathogen with diarrheal cases.

To determine the most prevalent bacterial pathogens of diarrhea in children, we used qPCR to detect and quantify the presence of EPEC, ETEC, EAEC, STEC, EIEC, *Shigella*, *Salmonella*, and *Y. enterocolitica*. Furthermore, in a subset of *eae*-positive samples, we analyzed by high-throughput real-time PCR the presence of genetic markers related to EHEC strains with clinical significance.

# MATERIALS AND METHODS

## Sample Collection, Ethical, and Biosafety Procedures

Stool samples were obtained from 110 children with diarrhea (54 boys and 56 girls) and 150 children without diarrhea (71 boys and 79 girls), aged from 1 month to 7 years old. Stool samples were collected from March 2008 through November 2010 at Hospital Municipal do Tatuapé, Hospital Infantil Cândido Fontoura, Hospital Infantil Menino Jesus and Centro Educacional Unificado (CEU-Butantã) in Sao Paulo city, SP, Brazil. Children did not display illnesses or comorbidities other than diarrhea, and they were not under antibiotic treatment for at least 3 months prior to sample collection. Diarrhea was defined as three or more unformed stools in the 24 h prior to enrollment. The Ethics Committee of the Biomedical Science Institute at the University of Sao Paulo approved this study (2006/743). All procedures followed the management of biosafety and biosecurity, since Biomedical Science Institute at São Paulo University guarantees excellent conditions for workers and population, and environmental protection from biohazards exposure.

## DNA Extraction

Total DNA from 200 mg of stool was obtained by using a QIAamp DNA Stool Mini Kit (QIAGEN) according to the manufacturer's instructions, resuspended in 100 µL of ultrapure water and stored at -80°C until use. DNA concentrations were determined by spectrophotometer (NanoDrop 2000, Thermo Scientific, USA), and 5 µL of each DNA sample were checked for integrity on 1% agarose gel.

## Identification of the Main Enteric Pathogens by qPCR

qPCR assays for genes *eae* (attaching and effacing lesions), *ipaH* (enteroinvasive mechanism), *lt* (heat-labile toxin), *st* (heat-stable toxin), *aggR* (aggregative fimbriae regulator), *invA* [inner membrane component of the *Salmonella* T3SS apparatus], and *ystA* (enterotoxin YstA *Yersinia* virulence plasmid) were carried out in duplicate and performed in a total volume of 25 µL, containing 2X TaqMan universal master mix (Applied Biosystems, USA), 10 µM of each primer, 10 µM of TaqMan probe, and 2 ng of DNA. Amplifications were performed in a
thermal cycler programmed as follows: denaturation at 95°C for 10 min, followed by 45 cycles of two steps: denaturation at 95°C for 15 s and an annealing temperature at 60°C for 1 min. The primer/probe sets used are shown in **Supplementary Table 1**. A standard curve was also derived using 10-fold DNA dilutions from the reference strains with their respective primer pairs. Amplifications were adjusted to $R^2 > 0.900$. A *P* < 0.05 was considered statistically significant. A sample was considered positive for a target gene when the detected fluorescence generated a curve above the background fluorescence, which was established by the Rotor-Gene Q6000 analytical software (Qiagen, Brazil).

## Differentiation of Enteric Pathogens by PCR

Samples positive for *eae* (*n* = 57) and *ipaH* (*n* = 15) by real-time PCR were further screened for the *eae*, *bfpA* (bundle-forming pilus for typical EPEC), *ipaH*, *lacY* (lactose permease), and *stx1* and *stx2* (Shiga toxins) genes by multiplex PCR according to the protocol of Aranda et al. (2004) and Pavlovic et al. (2011) for differentiation between EPEC/STEC and EIEC/Shigella samples.

## Characterization of the *eae*-positive Samples by High-Throughput Real-Time PCR

A LightCycler® 1536 (Roche, Meylan, France) was used to perform high-throughput real-time PCR amplifications as described previously (Delannoy et al., 2012b), except that 1 µL of sample DNA was used in each reaction for a final reaction volume of 2 µL. The thermal profile was modified as follows: 95°C for 1 min, followed by 45 cycles of 95°C for 0 s, and 58°C for 30 s. All ramp rates were set to 2°C/s. The primers and probes used targeted genes encoding intimin (*eae*, *eae-alpha*, *eae-beta*, *eae-gamma*, *eae-epsilon*, and *eae-theta*), O group-associated genes for the top 7 serogroups (*wzxO26*, *rfbE*_{0157}, *wzy*_{0145}, *wzx*_{0103}, *wbdI*_{011}, *wzx*_{0121}, *wzx*_{045}), flagellar antigens H11, H19, H2, H28, H7, and H8 (*fliCH*11, *fliCH*19, *fliCH*2, *fliCH*28, *fliCH*7, and *fliCH*8), urease (*ureD*), effector proteins translocated by T3SS (EspL2 [*ent* {or *espL2*}], NleB [*nleB*], NleE [*nleE*], NleA [*nleA*], NleF [*nleF*], EfA1 [*efa1*], EspN [*espN*], EspK [*espK*], and EspM1 [*espM1*]), and other genetic markers including Z6065, Z2098, Z2099, *terE*, *hlyA*, *ehxA*, and *pagC* and markers related to the CRISPR loci *SP_O157_A*, *SP_O157_B*, *SP_O157_C*, *SP_O26_C*, *SP_O26_D*, *SP_O26_E* (Perelle et al., 2004; Fratamico et al., 2009, 2010; Bugarel et al., 2010; Delannoy et al., 2012a, 2013, 2015, 2016; Piazza et al., 2013). The *wecA* gene was used as a reference genetic marker for *E. coli* (Bugarel et al., 2011). An inhibition control (IC) was introduced in each sample to check for potential inhibition of the PCR reaction due to intrinsic characteristics of the sample. IC is a recombinant pBluescriptISK+ plasmid containing the *dsb* gene from *Ehrlichia canis* (Michelet et al., 2014). The plasmid was added to each sample at a concentration of ~0.3 pg/µL. Primers and probe specific for the *Ehrlichia canis dsb* gene were used to detect IC (Michelet et al., 2014).

## Statistical Analyses

We compared the presence of pathogen between groups using the Fisher exact test. A comparison of bacterial loads in *eae-* and *aggR*-positive samples between samples from children with and without diarrhea was performed with a Mann-Whitney *U*-test. A *P* < 0.05 was considered statistically significant. All statistical analyses were performed in Stata v.12.1.

## RESULTS

### Prevalence of Enteropathogens in Stool Samples Obtained From Children With and Without Diarrhea

A total of 260 stool samples (110 from children with diarrhea and 150 children without diarrhea) were collected from public hospitals and schools, respectively, in Sao Paulo city (the largest city in South America) from March 2008 through November 2010. The groups studied consisted of children between 1 month and 7 years old, representing much of childhood. Each group was further divided into 3 age groups: 0–11 months, 12–23 months, and older than 24 months. The pathogens detected through their representative genetic markers in the stool samples from children with and without diarrhea of different age are reported in Table 1.

No genetic markers were detected in 39 of 110 (35%) children with diarrhea and 92 out of 150 (61%) without diarrhea. The presence of genetic markers associated with bacterial pathogens was significantly higher in stool samples from the diarrhea group (71/110) compared to the control (58/150) (*P* < 0.001). A significant increase in the number of samples positive for one pathogen was found in the diarrhea samples (46/110) compared to asymptomatic group (39/150) (*P* = 0.0036). Although we observed statistical differences between groups in samples with 3 or 4 pathogens, in both groups, the number of samples is low, not allowing further conclusions (Table 2).

### Bacterial Loads in *eae-* and *aggR*-Positive Samples

Quantification of bacterial loads was shown to be a promising tool to identify the etiological agent of diarrheal diseases. Of all the genes evaluated, *aggR* and *eae* were the markers most found in diarrheal and asymptomatic children (71 and 56 samples, respectively). We compared the copy number of *eae* and *aggR* genes in all the samples positive for these markers in both diarrhea and asymptomatic groups (Figure 1). A significantly higher copy number of *eae* was found in diarrhea samples (*P* < 0.039); no differences in *aggR* copy number were found between groups (*P* = 0.43).

### Prevalence of EHEC-Associated Virulence Markers in *eae*-positive Stool Samples From Children With and Without Diarrhea

In the control group, we detected 13% *eae*-positive samples and 11% *stx*-positive samples, while in the diarrhea group, 34% *eae*-positive samples, and 27%-*stx* positive samples were present,
<table><caption>**TABLE 1** | Distribution of virulence genes from enteric pathogens in fecal samples from children with and without diarrhea.</caption><thead><tr><th rowspan="2">Pathogens detected</th><th rowspan="2">Genetic markers</th><th colspan="4">Asymptomatic group (n = 150)</th><th colspan="4">Diarrhea group (n = 110)</th><th rowspan="2">Total</th><th rowspan="2">P-value</th></tr><tr><th>0-11 months (n = 54)</th><th>12-23 months (n = 59)</th><th>&gt; 24 months (n = 37)</th><th>Total</th><th>0-11 months (n = 55)</th><th>12-23 months (n = 50)</th><th>&gt; 24 months (n = 5)</th></tr></thead><tbody><tr><td>STEC</td><td>eae, stx1/stx2</td><td>2</td><td>4</td><td>0</td><td>6</td><td>2</td><td>9</td><td>0</td><td>11</td><td>0.19</td></tr><tr><td>EAEC</td><td>aagR</td><td>15</td><td>9</td><td>1</td><td>25</td><td>10</td><td>8</td><td>0</td><td>18</td><td>0.016</td></tr><tr><td>EPEC Typ</td><td>eae, bfp</td><td>0</td><td>1</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0.44</td></tr><tr><td>EPEC atyp</td><td>eae</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>1</td><td>0</td><td>4</td><td>0.033</td></tr><tr><td>ETEC</td><td>lt/st</td><td>3</td><td>1</td><td>0</td><td>4</td><td>1</td><td>3</td><td>0</td><td>4</td><td>0.38</td></tr><tr><td>EIEC</td><td>ipaH, LacY</td><td>0</td><td>1</td><td>2</td><td>3</td><td>4</td><td>3</td><td>1</td><td>8</td><td>0.088</td></tr><tr><td>STEC EAEC</td><td>eae, aggR, stx1/stx2</td><td>3</td><td>1</td><td>0</td><td>4</td><td>1</td><td>2</td><td>0</td><td>3</td><td>0.25</td></tr><tr><td>STEC EPECtyp</td><td>eae, bfp, stx1/stx2</td><td>1</td><td>1</td><td>0</td><td>2</td><td>3</td><td>5</td><td>0</td><td>8</td><td>0.049</td></tr><tr><td>STEC ETEC</td><td>eae, stx1/stx2, lt/st</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.13</td></tr><tr><td>EAEC ETEC</td><td>aggR, lt/st</td><td>3</td><td>3</td><td>0</td><td>6</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0.013</td></tr><tr><td>EAEC EPECtyp</td><td>aggR, eae, bfp</td><td>2</td><td>1</td><td>0</td><td>3</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.02</td></tr><tr><td>EAEC EIEC</td><td>aggR, ipaH, LacY</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0.182</td></tr><tr><td>EAEC EPECtyp</td><td>aggR, eae</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td><td>2</td><td>0.098</td></tr><tr><td>EAEC Salmonella</td><td>aagR, invA</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>1</td><td>0.18</td></tr><tr><td>STEC EAEC ETEC</td><td>eae, stx1/stx2, aggR, lt/st</td><td>1</td><td>1</td><td>0</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0.057</td></tr><tr><td>STEC EAEC EPECtyp</td><td>eae, aggR, stx1/stx2, bfp</td><td>0</td><td>1</td><td>0</td><td>1</td><td>2</td><td>0</td><td>1</td><td>3</td><td>0.2</td></tr><tr><td>STEC EPECtyp EIEC</td><td>eae, stx1/stx2, bfp, ipaH, LacY</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0.18</td></tr><tr><td>STEC EPECtyp ETEC</td><td>eae, stx1/stx2, bfp, lt/st</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0.18</td></tr><tr><td>STEC EPECtyp Yersinia</td><td>eae, stx1/stx2, bfp, invA</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0.18</td></tr><tr><td>STEC ETEC EIEC</td><td>eae, stx1/stx2, lt/st, ipaH, LacY</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>0.18</td></tr><tr><td>STEC EPECtyp ETEC EIEC</td><td>eae, stx1/stx2, bfp, lt/st, ipaH, Lacy</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>0</td><td>2</td><td>0.098</td></tr></tbody></table>

P < 0.05 was considered statistically significant.

<table><caption>**TABLE 2** | Association of enteric pathogens in fecal samples from children with and without diarrhea determined by virulence genes presence.</caption><thead><tr><th rowspan="2">Number of pathogens detected</th><th colspan="2">Asymptomatic group (n = 150)</th><th colspan="2">Diarrhea group (n = 110)</th></tr><tr><th>Total</th><th>Total</th><th>Total</th><th>P-value</th></tr></thead><tbody><tr><td>None</td><td>92</td><td>39</td><td></td><td>&gt;0.001</td></tr><tr><td>Positive for 1 or more pathogens</td><td>58</td><td>71</td><td></td><td>&gt;0.001</td></tr><tr><td>1 pathogen</td><td>39</td><td>46</td><td></td><td>0.0036</td></tr><tr><td>2 pathogens</td><td>16</td><td>16</td><td></td><td>0.1734</td></tr><tr><td>3 pathogens</td><td>3</td><td>7</td><td></td><td>0.0354</td></tr><tr><td>4 pathogens</td><td>0</td><td>2</td><td></td><td>0.0486</td></tr></tbody></table>

and thus, the presence of these genes were significantly higher in the diarrhea group. Surprisingly, our results indicated higher percentages of the majority of the additional genetic markers that are reported to be mostly predominant in EHEC in asymptomatic children (Table 3). Comparing all genetic markers, no differences were observed between the groups (P = 0.312), but the presence of the genes stx2, ehxA, nleB, nleE, nleF, espM1, pagC, ent, efa1, espV, Z2098, Z2099, Z6065, terE, ureD, and SP_O157 A, B, and C were statistically significant in the asymptomatic children compared to the diarrhea groups (P-values varying from 0.0001 to 0.0035) (Table 3).

DISCUSSION

A better knowledge of the distribution of pathogens linked to stool contamination is important for evaluating the pathogen loads in the environment as well as the associated health risks (Zhang et al., 2016). Over the last few decades, new molecular methods have become essential in allowing the rapid detection and characterization of pathogens (Miliwebsky et al., 2016).

The present study was performed to determine the molecular epidemiology of enteropathogens among children under seven years old with and without acute gastroenteritis in Sao Paulo, SP, Brazil. Virulence markers for DEC pathotypes were largely detected confirming some other surveys in Latin America and China (Ochoa et al., 2009; Moreno et al., 2010; Zhang et al., 2016; Lima A. A. M et al., 2019; Lima F. M. et al., 2019). Different results were reported in Western Kenya among bacterial pathogens, where two species of *Shigella* were the most prevalent (Swierczewski et al., 2013).

Our results showed that the genes encoding virulence factors for *Salmonella* and *Yersinia* were detected in very low frequencies, in agreement with another survey in counties near Sao Paulo city (Lima F. M. et al., 2019).

Herein, DEC-associated genetic markers were detected in 71 out of 110 (65%) children with diarrhea and 58 out of 150 (39%) without diarrhea and significant differences were observed between the two groups. These results are in agreement with other studies in South American countries. In Northeast Brazil, the presence of EPEC and EAEC were consistently associated
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Quantification of *eae*- and *aggR*- markers in stool. Positive stool samples for *aggR* (left) or *eae* (right) markers from children with (diarrhea group) or without (control) diarrhea were analyzed by quantitative qPCR. Data are presented as median log copy number/mg of stool sample ± IQR. \*P < 0.05 (Mann-Whitney *U*-test).</figcaption></figure>

with diarrhea (Maranhão et al., 2008; Moreno et al., 2010). Furthermore, EAEC and STEC showed by multivariate logistic regression analysis significant odds of being associated with a risk for diarrheal diseases as well as norovirus, adenovirus, rotavirus, and *Giardia* infection. Furthermore, EAEC enteric infections are associated with a high clinical severity score plus moderate-to-severe dehydration (Lima A. A. M et al., 2019).

The shifting in DEC epidemiology has been described in different studies and geographical regions (Pérez et al., 2010; Torres, 2017; Imdad et al., 2018; Pérez-Corrales and Leandro-Sandí, 2019) including Brazil, where until the 1990s, EPEC, mainly the typical EPEC serotypes producing EPEC adherence factor (EAF), was the main cause of infantile diarrhea (Gomes et al., 1991; Rosa et al., 1998; Regua-Mangia et al., 2004), and also *Salmonella* species (8%), ETEC (7%), and *Shigella* species (5%) were associated with diarrhea (Gomes et al., 1991), but it seems that they are becoming more and more rare, as presented herein and described elsewhere (Franzolin et al., 2005; Bueris et al., 2007; Moreno et al., 2010; Lima A. A. M et al., 2019; Lima F. M. et al., 2019; Ori et al., 2018).

The distribution of genes encoding either for one single pathogen in the diarrhea group (66%) and in control (67%) or in multiple coinfections, 34 and 33%, respectively, was similar in the present study. However, except for *Salmonella* plus *E. coli* and *Yersinia* plus *E. coli*, we cannot assure that we detected different strains, since they were not isolated from stools. Nevertheless, coinfections have been described, suggesting associations among potential enteropathogens in the etiology of diarrhea (Orlandi et al., 2006; Lima F. M. et al., 2019). Indeed, some hybrid strains possessing typical genes from different pathotypes (heteropathogens) have recently emerged (Santos et al., 2020). For example, the O104:H4 strain associated with the severe German outbreak in 2011 was a hybrid of both EAEC and EHEC (Brzuszkiewicz et al., 2011). More recently, enterohemorrhagic *E. coli* hybrid pathotype O80:H2 has emerged in Europe as an

important serotype responsible for hemolytic uremic syndrome (HUS) associated with extraintestinal infections in children (Cointe et al., 2018). Also, coinfections have been reported, detecting in stool samples the presence of EPEC and EAEC and *S. enterica houtenae* and EAEC (Lima F. M. et al., 2019).

The incidence of intestinal infection caused by STEC is low in the Brazilian population (de Souza et al., 2011; Ori et al., 2018; Lima F. M. et al., 2019). Thus, it was surprising to find more EHEC-associated genetic markers in the asymptomatic group than in the diarrhea group, since these genetic markers are usually more associated with *E. coli* strains of greater clinical significance (Delannoy et al., 2012a,b, 2016; Piazza et al., 2013). This suggests that these potential EHEC strains are not the etiological agent of diarrhea in this sampled group, which remains to be investigated.

Current methods for either screening or detecting of enteropathogens by PCR, multiplex PCR, and qPCR have been used worldwide, given the high sensitivity and specificity, which are satisfactory parameters to ensure an adequate diagnosis (Panchalingam et al., 2012). However, the presence of several pathogens in diarrhea samples, as we observed, and the high frequency of enteropathogens in asymptomatic children has made it difficult to use these methods as a tool to identify the etiological agent of diarrhea.

Here, we found a significantly higher copy number of *eae* in diarrhea samples compared to samples from asymptomatic children, suggesting that LEE-containing pathogens could be the etiological agents responsible for diarrhea in children. Even though we did not find differences between groups in the number of copies for *aggR*, our data must be complemented with a larger number of samples and/or selection of other EAEC markers to support the use of qPCR for molecular diagnosis for this pathogen.

Our study had the following limitations. The number of samples was limited and we did not search for other
<table><caption>**TABLE 3** | Virulence profile determination of *eae*+ *Escherichia coli* from fecal samples of children with and without diarrhea*.</caption><thead><tr><th>Genetic markers</th><th>32 fecal samples with diarrhea (%)</th><th>19 fecal samples without diarrhea (%)</th><th>P-value</th></tr></thead><tbody><tr><td>bfp</td><td>1 (3.1)</td><td>0 (0)</td><td>0.2192</td></tr><tr><td>stx1</td><td>1 (3.1)</td><td>3 (15.8)</td><td>0.0513</td></tr><tr><td>stx2</td><td>8 (25.0)</td><td>12 (63.2)</td><td>0.0035</td></tr><tr><td>ehxA</td><td>6 (18.8)</td><td>13 (68.4)</td><td>0.0002</td></tr><tr><td>hlyA</td><td>12 (37.5)</td><td>7 (37.0)</td><td>0.4858</td></tr><tr><td>nleA</td><td>5 (15.6)</td><td>2 (10.5)</td><td>0.3043</td></tr><tr><td>nleB</td><td>12 (37.5)</td><td>14 (73.7)</td><td>0.0062</td></tr><tr><td>nleE</td><td>12 (37.5)</td><td>15 (79.0)</td><td>0.0020</td></tr><tr><td>nleF</td><td>13 (40.6)</td><td>15 (79.0)</td><td>0.0039</td></tr><tr><td>espM1</td><td>12 (37.5)</td><td>14 (73.7)</td><td>0.0062</td></tr><tr><td>pagC</td><td>6 (18.8)</td><td>15 (79.0)</td><td>0.0001</td></tr><tr><td>ent</td><td>7 (22.0)</td><td>15 (79.0)</td><td>0.0001</td></tr><tr><td>efa1</td><td>3 (9.4)</td><td>9 (47.4)</td><td>0.0010</td></tr><tr><td>espK</td><td>6 (18.8)</td><td>4 (21.0)</td><td>0.4242</td></tr><tr><td>espN</td><td>1 (3.1)</td><td>1 (5.3)</td><td>0.3477</td></tr><tr><td>espV</td><td>8 (25.0)</td><td>13 (68.4)</td><td>0.0012</td></tr><tr><td>Z2098</td><td>7 (22.0)</td><td>12 (63.2)</td><td>0.0016</td></tr><tr><td>Z2099</td><td>12 (37.5)</td><td>16 (84.2)</td><td>0.0006</td></tr><tr><td>Z6065</td><td>12 (37.5)</td><td>13 (68.4)</td><td>0.0164</td></tr><tr><td>terE</td><td>29 (91.0)</td><td>18 (94.7)</td><td>0.3151</td></tr><tr><td>ureD</td><td>14 (44.0)</td><td>15 (79.0)</td><td>0.0073</td></tr><tr><td>SP_O157-A</td><td>6 (18.8)</td><td>11 (58.0)</td><td>0.0021</td></tr><tr><td>SP_O157-B</td><td>1 (3.1)</td><td>7 (37.0)</td><td>0.0007</td></tr><tr><td>SP_O157-C</td><td>6 (18.8)</td><td>11 (58.0)</td><td>0.0021</td></tr><tr><td>O26_C</td><td>4 (12.5)</td><td>1 (5.3)</td><td>0.2017</td></tr><tr><td>O26_D</td><td>2 (6.2)</td><td>1 (5.3)</td><td>0.4474</td></tr></tbody></table>

P < 0.05 was considered statistically significant. *General P value between groups 0.312.

relevant pathogens, such as rotavirus, norovirus, protozoa, and helminthes. For bacterial loads, we only evaluated two markers, which might be insufficient, considering that DEC harbors several virulence determinants that can be used as a target for diagnosis. Another limitation of the study was that we did not isolate the pathogenic strains. This would be helpful to identify if asymptomatic children might shed the bacteria and possibly be contagious.

In conclusion, we found that the presence of genetic markers associated with bacterial pathogens was significantly higher in stool samples from the diarrhea group compared to the control. Also, compared to samples from asymptomatic children, there was a significantly higher copy number of the *eae* gene in

diarrhea samples. The data presented herein demonstrated that asymptomatic children of an urban area in Brazil might be enteropathogen reservoirs, especially for STEC, thus favoring their widespread.

## DATA AVAILABILITY STATEMENT

All datasets generated for this study are included in the article/Supplementary Material.

## ETHICS STATEMENT

The studies involving human participants were reviewed and approved by Ethics Committee of the Biomedical Science Institute at the University of Sao Paulo. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

## AUTHOR CONTRIBUTIONS

VM, VN, SD, PF, RP, and MA-C participated in the design of the study. VM, VN, SD, PF, MF, RP, and MA-C participated in data analysis. MF, RP, SD, and PF participated in writing the manuscript. VM and VN carried out the qPCR experiments. GA carried out the multiplex PCR experiments. SD and PF carried out the high-throughput real-time PCR experiments. All authors read and approved the final manuscript.

## FUNDING

This study was supported by grants from CNPq 158799/2012-7, FAPESP 2013/17739-9, 2017/14821-7, and FONDECYT 1200994.

## ACKNOWLEDGMENTS

The authors thank Mrs. Marcia H. Fukugaiti for her technical support. Standard strains were kindly donated by Dr. Juliana Pfrimer Falcão, Faculdade de Ciências Farmacêuticas de Ribeirão Preto (FCFRP-USP), SP, Brazil and Dr. Marcia Regina Franzolin from Instituto Butantan, SP, Brazil. Dr. A. Leyva (USA) helped with English editing of the manuscript.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2020.549919/full#supplementary-material

## REFERENCES

Aranda, K. R. S., Fagundes-Neto, U., and Scaletsky, I. C. A. (2004). Evaluation of multiplex PCRs for diagnosis of infection with diarrheagenic *Escherichia coli* and *Shigella* spp. *J. Clin. Microbiol.* 42, 5849-5853. doi: 10.1128/JCM.42.12.5849-5853.2004

Bottone, E. J. (1999). *Yersinia enterocolitica*: overview and epidemiologic correlates. *Microbes Infect.* 1, 323-333. doi: 10.1016/S1286-4579(99) 80028-8

Brzuszkiewicz, E., Thürmer, A., Schuldes, J., Leimbach, A., Liesegang, H., Meyer, F. D., et al. (2011). Genome sequence analyses of two isolates from the recent *Escherichia coli* outbreak in Germany reveal the emergence of a new pathotype:
entero-aggregative-haemorrhagic *Escherichia coli* (EAHEC). *Arch. Microbiol.* 193, 883-891. doi: 10.1007/s00203-011-0725-6

Bueris, V., Sircili, M. P., Taddei, C. R., dos Santos, M. F., Franzolin, M. R., Martinez, M. B., et al. (2007). Detection of diarrheagenic *Escherichia coli* from children with and without diarrhea in Salvador, Bahia, Brazil. *Mem. Inst. Oswaldo Cruz.* 102, 839-844. doi: 10.1590/S0074-02762007005000116

Bugarel, M., Beutin, L., and Fach, P. (2010). Low-density macroarray targeting non-locus of enterocyte effacement effectors (nle genes) and major virulence factors of Shiga toxin-producing *Escherichia coli* (STEC): a new approach for molecular risk assessment of STEC isolates. *Appl. Environ. Microbiol.* 76, 203-211. doi: 10.1128/AEM.01921-09

Bugarel, M., Beutin, L., Scheutz, F., Loukiadis, E., and Fach, P. (2011). Identification of genetic markers for differentiation of Shiga toxin-producing, enteropathogenic, and avirulent strains of *Escherichia coli* O26. *Appl. Environ. Microbiol.* 77, 2275-2281. doi: 10.1128/AEM.02832-10

Cointe, A., Birgy, A., Mariani-Kurkdjian, P., Liguori, S., Courroux, C., Blanco, J., et al. (2018). Emerging multidrug-resistant hybrid pathotype Shiga toxin-producing *Escherichia coli* O80 and related strains of clonal complex 165, Europe. *Emerg. Infect Dis.* 24, 2262-2269. doi: 10.3201/eid2412.180272

Croxen, M. A., Law, R. J., Scholz, R., Keeney, K. M., Wlodarska, M., and Finlay, B. B. (2013). Recent advances in understanding enteric pathogenic *Escherichia coli*. *Clin. Microbiol. Rev.* 26, 822-880. doi: 10.1128/CMR.00022-13

de Souza, R. L., Carvalhaes, J. T. A., Nishimura, L. S., de Andrade, M. C., and Guth, B. E. C. (2011). Hemolytic uremic syndrome in pediatric intensive care units in São Paulo, Brazil. *Open Microbiol J.* 5, 76-82. doi: 10.2174/1874285801105010076

Delannoy, S., Beutin, L., Burgos, Y., and Fach, P. (2012a). Specific detection of enteroaggregative hemorrhagic *Escherichia coli* O104:H4 strains by use of the CRISPR locus as a target for a diagnostic real-time PCR. *J. Clin. Microbiol.* 50, 3485-3492. doi: 10.1128/JCM.01656-12

Delannoy, S., Beutin, L., and Fach, P. (2012b). Use of clustered regularly interspaced short palindromic repeat sequence polymorphisms for specific detection of enterohemorrhagic *Escherichia coli* strains of serotypes O26:H11, O45:H2, O103:H2, O111:H8, O121:H19, O145:H8, and O157:H7 by real time PCR. *J. Clin. Microbiol.* 50, 4035-4040. doi: 10.1128/JCM.02097-12

Delannoy, S., Beutin, L., and Fach, P. (2013). Towards a molecular definition of enterohemorrhagic *Escherichia coli* (EHEC): detection of genes located on O island 57 as markers to distinguish EHEC from closely related enteropathogenic *Escherichia coli* strains. *J. Clin. Microbiol.* 51, 1083-1088. doi: 10.1128/JCM.02864-12

Delannoy, S., Chaves, B. D., Ison, S. A., Webb, H. E., Beutin, L., Delaval, J., et al. (2016). Revisiting the STEC testing approach: using *espK* and *espV* to make enterohemorrhagic *Escherichia coli* (EHEC) detection more reliable in beef. *Front. Microbiol.* 7:1. doi: 10.3389/fmicb.2016.00001

Delannoy, S., Mariani-Kurkdjian, P., Bonacorsi, S., Liguori, S., Ison, S. A., and Fach, P. (2015). Draft genome sequences of human-pathogenic *Escherichia coli* O26:H11 strains carrying the *stx2* gene only and circulating in France. *Genome Announc.* 3:e00852-e00815. doi: 10.1128/genomeA.00852-15

Dias, R. C. B., Tanabe, R. H. S., Vieira, M. A., Cergole-Novella, M. C., Dos Santos, L. F., Gomes, T. A. T., et al. (2020). Analysis of the virulence profile and phenotypic features of typical and atypical enteroaggregative *Escherichia coli* (EAEC) isolated from diarrheal patients in Brazil. *Front. Cell. Infect. Microbiol.* 10:144. doi: 10.3389/fcimb.2020.00144

Franzolin, M. R., Alves, R. C., Keller, R., Gomes, T. A., Beutin, L., Barreto, M. L., et al. (2005). Prevalence of diarrheagenic *Escherichia coli* in children with diarrhea in Salvador, Bahia, Brazil. *Mem. Inst. Oswaldo Cruz.* 100, 359-363. doi: 10.1590/S0074-0276200500040004

Fratamico, P. M., DebRoy, C., Miyamoto, T., and Liu, Y. (2009). PCR detection of enterohemorrhagic *Escherichia coli* O145 in food by targeting genes in the *E. coli* O145 O-antigen gene cluster and the Shiga toxin 1 and Shiga toxin 2 genes. *Foodborne Pathog. Dis.* 6, 605-611. doi: 10.1089/fpd.2008.0254

Fratamico, P. M., Yan, X. H., Liu, Y. H., DebRoy, C., Byrne, B., Monaghan, A., et al. (2010). *Escherichia coli* serogroup O2 and O28ac O-antigen gene cluster sequences and detection of pathogenic *E. coli* O2 and O28ac by PCR. *Can. J. Microbiol.* 56, 308-316. doi: 10.1139/W10-010

Freitas Neto, O. C., Penha Filho, R. A. C., Barrow, P., and Berchieri Junior, A. (2010). Sources of human non-typhoid Salmonellosis: a review. *Rev. Bras. Cienc. Avic.* 12, 1-11. doi: 10.1590/S1516-635X2010000100001

Galán, J. E., Ginocchio, C., and Costeas, P. (1992). Molecular and functional characterization of the *Salmonella* invasion gene InvA homology of InvA to members of a new protein family. *J. Bacteriol.* 174, 4338-4349. doi: 10.1128/JB.174.13.4338-4349.1992

Ginocchio, C. C., and Galán, J. E. (1995). Functional conservation among members of the *Salmonella typhimurium* InvA family of proteins. *Infect. Immun.* 63, 729-732. doi: 10.1128/IAI.63.2.729-732.1995

Gomes, T. A., Elias, W. P., Scaletsky, I. C., Guth, B. E., Rodrigues, J. F., Piazza, R. M. F., et al. (2016). Diarrheagenic *Escherichia coli*. *Braz. J. Microbiol.* 47(Suppl. 1), 3-30. doi: 10.1016/j.bjm.2016.10.015

Gomes, T. A., Rassi, V., MacDonald, K. L., Ramos, S. R., Trabulsi, L. R., Vieira, M. A., et al. (1991). Enteropathogens associated with acute diarrheal disease in urban infants in São Paulo, Brazil. *J. Infect. Dis.* 164, 331-337. doi: 10.1093/infdis/164.2.331

Hernandes, R. T., Elias, W. P., Vieira, M. A., and Gomes, T. A. (2009). An overview of atypical enteropathogenic *Escherichia coli*. *FEMS. Microbiol. Lett.* 297, 137-149. doi: 10.1111/j.1574-6968.2009.01664.x

Imdad, A., Foster, M. A., Iqbal, J., Fonnesbeck, C., Payne, D. C., Zhang, C., et al. (2018). Diarrheagenic *Escherichia coli* and acute gastroenteritis in children in Davidson County, Tennessee, United States: a case-control study. *Pediatr. Infect. Dis. J.* 37, 543-548. doi: 10.1097/INF.0000000000001908

Kaper, J. B., Nataro, J. P., and Mobley, H. L. (2004). Pathogenic *Escherichia coli*. *Nat. Rev. Microbiol.* 2, 123-140. doi: 10.1038/nrmicro818

Leli, C., Di Matteo, L., Gotta, F., Vay, D., Cavallo, V., Mazzeo, R., et al. (2020). Evaluation of a multiplex gastrointestinal PCR panel for the aetiological diagnosis of infectious diarrhoea. *Infect. Dis.* 52, 114-120. doi: 10.1080/23744235.2019.1688861

Lima, A. A. M., Oliveira, D. B., Quetz, J. S., Havt, A., Prata, M. M. G., Lima, I. F. N., et al. (2019). Etiology and severity of diarrheal diseases in infants at the semiarid region of Brazil: a case-control study. *PLoS Negl. Trop. Dis.* 13:e0007154. doi: 10.1371/journal.pntd.0007154

Lima, F. M., de Monteiro, F., Mucci, E. R., Ahagon, C. M., Dos Santos Carmo, A. M., Eterovic, A., et al. (2019). Epidemiology and genetic screening of diarrheagenic *Escherichia coli* among symptomatic and asymptomatic children. *J. Med. Microbiol.* 68, 1033-1041. doi: 10.1099/jmm.0.001020

Liu, L., Johnson, H. L., Cousens, S., Perin, J., Scott, S., Lawn, J. E., et al. (2012). Global, regional, and national causes of child mortality: an updated systematic analysis for 2010 with time trends since 2000. *Lancet* 379, 2151-2161. doi: 10.1016/S0140-6736(12)60560-1

Maranhão, H. S., Medeiros, M. C., Scaletsky, I. C., Fagundes-Neto, U., and Morais, M. B. (2008). The epidemiological and clinical characteristics and nutritional development of infants with acute diarrhoea, in north-eastern Brazil. *Ann. Trop. Med. Parasitol.* 102, 357-365. doi: 10.1179/136485908X 278865

Marinou, I., Bersimis, S., Ioannidis, A., Nicolaou, C., Mitroussia-Ziouva, A., Legakis, N. J., et al. (2012). Identification and antimicrobial resistance of Campylobacter species isolated from animal sources. *Front. Microbiol.* 3:58. doi: 10.3389/fmicb.2012.00058

Michelet, L., Delannoy, S., Devillers, E., Umhang, G., Aspan, A., Juremalm, M., et al. (2014). High-throughput screening of tick-borne pathogens in Europe. *Front. Cell. Infect. Microbiol.* 4:103. doi: 10.3389/fcimb.2014.00103

Miliwebsky, E., Schelotto, F., Varela, G., Luz, D., Chinen, I., and Piazza, R. M. F. (2016). "Human diarrheal infections: diagnosis of diarrheagenic *Escherichia coli* pathotypes," in *Escherichia coli in the Americas*, ed A. G. Torres (Cham: Springer), 343-369. doi: 10.1007/978-3-319-45092-6_15

Moreno, A. C., Fernandes-Filho, A., Gomes, T. A., Ramos, S. T. S., Montemor, L. P., Tavares, V. C., et al. (2010). Etiology of childhood diarrhea in the northeast of Brazil: significant emergent diarrheal pathogens. *Diagn. Microbiol. Infect. Dis.* 66, 50-57. doi: 10.1016/j.diagmicrobio.2008.03.017

Nunes, M. R., Magalhães, P. P., Penna, F. J., Nunes, J. M., and Mendes, E. N. (2012). Diarrhea associated with *Shigella* in children and susceptibility to antimicrobials. *J. Pediatr.* 88, 125-128. doi: 10.2223/JPED.2131

Ochoa, T. J., Ecker, L., Barletta, F., Mispireta, M. L., Gil, A. I., Contreras, C., et al. (2009). Age-related susceptibility to infection with diarrheagenic *Escherichia coli* among infants from periurban areas in Lima, Peru. *Clin. Infect. Dis.* 49, 1694-1702. doi: 10.1086/648069

Ori, E. L., Takagi, E. H., Andrade, T. S., Miguel, B. T., Cergole-Novella, M. C., Guth, B. E. C., et al. (2018). Diarrhoeagenic *Escherichia coli* and *Escherichia albertii* in
Brazil: pathotypes and serotypes over a 6-year period of surveillance. *Epidemiol. Infect.* 147:e10. doi: 10.1017/S0950268818002595

Orlandi, P. P., Magalhães, G. F., Matos, N. B., Silva, T., Penatti, M., Nogueira, P. A., et al. (2006). Etiology of diarrheal infections in children of Porto Velho (Rondonia, Western Amazon region, Brazil). *Braz. J. Med. Biol. Res.* 39, 507-517. doi: 10.1590/S0100-879X2006000400011

Panchalingam, S., Antonio, M., Hossain, A., Mandomando, I., Ochieng, B., Oundo, J., et al. (2012). Diagnostic microbiologic methods in the GEMS-1 case/control study. *Clin. Infect. Dis.* 55(Suppl. 4), S294–302. doi: 10.1093/cid/cis754

Pavlovic, M., Luze, A., Konrad, R., Buger, A., Sing, A., Busch, U., et al. (2011). Development of a duplex real time PCR for differentiation between *Escherichia coli* and *Shigella* spp. *J. Appl. Microbiol.* 110, 1245–1251. doi: 10.1111/j.1365-2672.2011.04973.x

Perelle, S., Dilasser, F., Grout, J., and Fach, P. (2004). Detection by 5'- nuclease PCR of Shiga-toxin producing *Escherichia coli* O26, O55, O91, O103, O111, O113, O145 and O157:H7, associated with the world's most frequent clinical cases. *Mol. Cell Probes.* 18, 185-192. doi: 10.1016/j.mcp.2003.12.004

Pérez, C., Gómez-Duarte, O. G., and Arias, M. L. (2010). Diarrheagenic *Escherichia coli* in children from Costa Rica. *Am. J. Trop. Med. Hyg.* 83, 292–297. doi: 10.4269/ajtmh.2010.09-0412

Pérez-Corrales, C., and Leandro-Sandi, K. (2019). Diarrheagenic *Escherichia coli* in Costa Rican children: a 9-year retrospective study. *BMC Res. Notes.* 12:297. doi: 10.1186/s13104-019-4313-1

Peruzy, M. F., Murru, N., Perugini, A. G., Capuano, F., Delibato, E., Mercogliano, R., et al. (2017). Evaluation of virulence genes in *Yersinia enterocolitica* strains using SYBR Green real-time PCR. *Food Microbiol.* 65, 231–235. doi: 10.1016/j.fm.2017.03.004

Piazza, R. M. F., Abe, C. M., Horton, D. S. P. Q., Miliwebsky, E., Chinen, I., Vaz, T. M. I., et al. (2010). "Detection and subtyping methods of diarrheagenic *Escherichia coli* strains," in *Pathogenic Escherichia coli* in Latin America, ed A. G. Torres (Sharjah: Bentham Science Publishers), 95–115.

Piazza, R. M. F., Delannoy, S., Fach, P., Saridakis, H. O., Pedroso, M. Z., Rocha, L. B., et al. (2013). Molecular and phenotypic characterization of *Escherichia coli* O26:H8 among diarrheagenic *Escherichia coli* O26 strains isolated in Brazil. *Appl. Environ. Microbiol.* 79, 6847–6854. doi: 10.1128/AEM.01693-13

Regua-Mangia, A. H., Gomes, T. A. T., Vieira, M. A. M., Andrade, J. R., Irino, K., and Teixeira, L. M. (2004). Frequency and characteristics of diarrhoeagenic *Escherichia coli* strains isolated from children with and without diarrhoea in Rio de Janeiro, Brazil. *J. Infect.* 48, 161–167. doi: 10.1016/S0163-4453(03)00138-5

Rosa, A. C., Mariano, A. T., Pereira, A. M., Tibana, A., Gomes, T. A. T., and Andrade, J. R. (1998). Enteropathogenicity markers in *Escherichia coli* isolated from infants with acute diarrhoea and healthy controls in Rio de Janeiro, Brazil. *J. Med. Microbiol.* 47, 781–790. doi: 10.1099/00222615-47-9-781

Santos, A. C. M., Santos, F. F., Silva, R. M., and Gomes, T. A. T. (2020). Diversity of hybrid- and hetero-pathogenic *Escherichia coli* and their potential implication in more severe diseases. *Front. Cell. Infect. Microbiol.* 10:339. doi: 10.3389/fcimb.2020.00339

Swierczewski, B. E., Odundo, E. A., Koech, M. C., Ndonye, J. N., Kirera, R. K., Odhiambo, C. P., et al. (2013). Enteric pathogen surveillance in a case-control study of acute diarrhoea in the town of Kisii, Kenya. *J. Med. Microbiol.* 62, 1774–1776. doi: 10.1099/jmm.0.059139-0

Torres, A. G. (2017). *Escherichia coli* diseases in Latin America-a “One Health” multidisciplinary approach. *Pathog. Dis.* 1:75. doi: 10.1093/femspd/ftx012

World Health Organization (2017). *Diarrhoeal Disease*. Available online at: https://www.who.int/news-room/fact-sheets/detail/diarrhoeal-disease (accessed July 21, 2020).

Zhang, S.-X., Yang, C.-L., Gu, W.-P., Ai, L., Serrano, E., Yang, P., et al. (2016). Case-control study of diarrheal disease etiology in individuals over 5 years in Southwest China. *Gut Pathog.* 8:58. doi: 10.1186/s13099-016-0141-1

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Merino, Nakano, Delannoy, Fach, Alberca, Farfan, Piazza and Avila-Campos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.