<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="brief-report" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Environ. Sci.</journal-id>
<journal-title>Frontiers in Environmental Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Environ. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-665X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1638383</article-id>
<article-id pub-id-type="doi">10.3389/fenvs.2025.1638383</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Environmental Science</subject>
<subj-group>
<subject>Brief Research Report</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>
<italic>Escherichia coli</italic> thrives in soil 24 months after grazing exclusion in a rainfed Mediterranean biodiverse pasture</article-title>
<alt-title alt-title-type="left-running-head">Verheijen et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenvs.2025.1638383">10.3389/fenvs.2025.1638383</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Verheijen</surname>
<given-names>Frank</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2036418/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Jelin&#x10d;i&#x107;</surname>
<given-names>Antun</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3097085/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jeffery</surname>
<given-names>Simon</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Domingos</surname>
<given-names>Tiago</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Khodaparast</surname>
<given-names>Zahra</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2548292/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Bastos</surname>
<given-names>Ana Catarina</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Environment and Planning, Centre for Environmental and Marine Studies (CESAM), University of Aveiro</institution>, <addr-line>Aveiro</addr-line>, <country>Portugal</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Crop and Environment Sciences Department</institution>, <institution>Harper Adams University</institution>, <addr-line>Newport</addr-line>, <country>United Kingdom</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Marine</institution>, <institution>Environment and Technology Centre (MARETEC)</institution>, <institution>LARSyS</institution>, <institution>Instituto Superior T&#x00E9;cnico</institution>, <institution>University of Lisbon</institution>, <addr-line>Lisbon</addr-line>, <country>Portugal</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2531310/overview">Lyudmyla Symochko</ext-link>, University of Coimbra, Portugal</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2024143/overview">Xiaomeng Chen</ext-link>, Northeast Agricultural University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2168756/overview">Mitiku Mihiret Seyoum</ext-link>, University of Arkansas, United States</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Antun Jelin&#x10d;i&#x107;, <email>antunjelincic@gmail.com</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>29</day>
<month>08</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="ecorrected">
<day>26</day>
<month>09</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="collection">
<year>2025</year>
</pub-date>
<volume>13</volume>
<elocation-id>1638383</elocation-id>
<history>
<date date-type="received">
<day>30</day>
<month>05</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>24</day>
<month>07</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2025 Verheijen, Jelin&#x10d;i&#x107;, Jeffery, Domingos, Khodaparast and Bastos.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Verheijen, Jelin&#x10d;i&#x107;, Jeffery, Domingos, Khodaparast and Bastos</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<p>
<italic>Escherichia coli</italic> is a facultative anaerobic bacterium mainly found in the lower intestines of mammals. It is transferred to soils via animal feces and can splash onto plant surfaces during rainfall or be taken up through roots into the plant tissue. Subsequently, it can pose a risk to farm workers handling the soil or the vegetation, or groundwater quality through leaching. Many inoculation studies show that <italic>E. coli</italic> can survive in soils for weeks to months. The few field studies that have been published suggest <italic>E. coli</italic> can survive for longer in cold and wet soils. We analyzed 20 topsoil samples from a relatively warm and aerobic Mediterranean pasture soil where grazing, and therefore manure inputs, had been excluded for 24&#xa0;months, using the shotgun DNA sequencing method. We found that <italic>E. coli</italic> had not only survived for 2&#xa0;years, but it had also comparatively thrived, as evidenced by its relative abundance of 17.0% in the control treatment and 20.2% in the biochar-amended treatment (no significant difference). Within-treatment variation was low for the control soil (coefficient of variation of 0.15) and approximately double that for the biochar-amended soil (0.28). There were no correlations between the relative abundance of <italic>E. coli</italic> and soil organic matter content, available P, soil texture, pH, or other soil nutrients. Our results suggest the existence of a naturalized <italic>E. coli</italic> soil population, with potentially important implications for farm management, as well as for environmental and public health under the One Health umbrella.</p>
</abstract>
<kwd-group>
<kwd>soil</kwd>
<kwd>biochar</kwd>
<kwd>
<italic>E. coli</italic>
</kwd>
<kwd>grassland</kwd>
<kwd>shotgun DNA</kwd>
</kwd-group>
<counts>
<page-count count="7"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Toxicology, Pollution and the Environment</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>
<italic>Escherichia coli</italic> is a facultative anaerobic bacterium. It is a mesophilic species with optimum growth from 37&#xb0;C to 42&#xb0;C, with growth greatly impeded at lower temperatures and declining rapidly at higher temperatures. <italic>E. coli</italic> has been intensively studied as it serves as a model organism for understanding fundamental biological processes, including gene expression, replication, and metabolism. It is also widely used in medicine for producing proteins and enzymes, as well as containing several pathogenic strains.</p>
<p>
<italic>E. coli</italic> in soils has mostly been investigated in survival studies under laboratory conditions, with soil samples from the field inoculated with specific strains of the bacterium. In these studies, <italic>E. coli</italic> abundance, commonly measured in CFUs, tends to decrease rapidly after inoculation before levelling off until the bacterium is found to be below detection limits. Reported survival periods range from weeks to months to 336&#xa0;days, with longer survival times found for soils that were wetter and for generic rather than pathogenic strains (<xref ref-type="bibr" rid="B16">Murphy et al., 2024</xref>).</p>
<p>Traditionally, <italic>E. coli</italic> is considered an indicator of recent fecal contamination. Survival rates of <italic>E. coli</italic> in soils have been found to decline with high temperatures and aerobic conditions and low soil clay and dissolved organic carbon (DOC) contents (<xref ref-type="bibr" rid="B24">van Elsas et al., 2011</xref>). Microbial community structure and soil management practices were also found to have complex influences on <italic>E. coli</italic> survival rates. Nonetheless, evidence has suggested that certain strains of <italic>E. coli</italic> can persist in soil environments, particularly in relatively wet and cold pasture soils of Ireland (<xref ref-type="bibr" rid="B26">Brennan et al., 2010a</xref>) and the Alps (<xref ref-type="bibr" rid="B23">Texier et al., 2008</xref>). Studies have shown that <italic>E. coli</italic> can survive for extended periods of time in soils amended with manure, suggesting potential for establishment of naturalized populations. The timeframe of survival is an important consideration, given <italic>E. coli&#x2019;s</italic> pathogenic potential, and that restrictions on manure and sewage sludge application to soil are time restricted when growing fresh produce (i.e., soft fruits and vegetables; <xref ref-type="bibr" rid="B29">Department for Environment, Food &#x0026; Rural Affairs, 2018</xref>). Factors such as soil moisture, clay content, vegetation density, and rooting can contribute to enhanced persistence of <italic>E. coli</italic> in soils. Some hydrochars have been shown to adsorb <italic>E. coli</italic> in column experiments (<xref ref-type="bibr" rid="B7">Chung et al., 2017</xref>). However, the impact of biochar on <italic>E. coli</italic> survival in soil is variable, ranging from no effect with six different biochar types (<xref ref-type="bibr" rid="B2">Alegbeleye and Sant&#x2019;Ana, 2023</xref>) to having biocidal properties against the pathogenic strain <italic>E. coli</italic> O157:H7, with the efficacy of the biocidal properties being dependent on biochar production temperature and biochar age (<xref ref-type="bibr" rid="B9">Gurtler et al., 2020</xref>). In rainfed pastures, the presence of grazing livestock can contribute further to the introduction and maintenance of <italic>E coli</italic> populations in soil (<xref ref-type="bibr" rid="B23">Texier et al., 2008</xref>; <xref ref-type="bibr" rid="B15">Muirhead, 2023</xref>). Naturalization of such communities may affect farm soil productivity, as well as have broader off-farm implications for environmental quality and public health.</p>
<p>In their mini-review, <xref ref-type="bibr" rid="B24">van Elsas et al. (2011)</xref> questioned the assumption that <italic>E. coli</italic> shows natural declines in open environments and stated that any environmental reservoir of the bacterium can cause risks to human health: &#x201c;The factors that determine the rate of survival of particular <italic>E. coli</italic> strains, and thus the risks (to human health), are not (yet) easily predictable, and our capability to understand the effects of the secondary habitat (especially soil and water resources) on <italic>E. coli</italic> behaviour will be paramount to our abilities to manage the organism from both environmental and public health perspectives.&#x201d;</p>
<p>Therefore, to help reduce this knowledge gap, in this study, we aimed to quantify the relative abundance of <italic>E. coli</italic> in topsoil samples of a rainfed biodiverse pasture, collected 24&#xa0;months after exclusion of grazing animals, and to investigate correlations between <italic>E. coli</italic> survival and soil physicochemical variables, including biochar amendment, in a Mediterranean pasture soil.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>2 Methods and materials</title>
<sec id="s2-1">
<title>2.1 Study site and auxiliary data</title>
<p>The study site was located in the rain shadow of the Serra da Estrela mountain range in Portugal, that is, in the foothills directly to the east of Covilh&#xe3; on a 500&#xa0;ha livestock production farm called Quinta da Fran&#xe7;a, which is managed both commercially and as a pilot farm by Terraprima (<ext-link ext-link-type="uri" xlink:href="https://www.terraprima.pt/en/">https://www.terraprima.pt/en/</ext-link>). The study site has a hot-summer Mediterranean climate (Csa in the K&#xf6;ppen classification) with wet and mild winters (annual rainfall in the range 600&#x2013;1,200&#xa0;mm). Soils are formed from granitic parent material, resulting in relatively deep (70&#x2013;80&#xa0;cm) dystric Cambisols (<xref ref-type="bibr" rid="B10">IUSS, 2014</xref>) of a sandy loam texture, with a topsoil organic matter content of 2.1%&#x2013;3.6%, and strongly acidic (pH in 1&#xa0;M KCl: 3.8&#x2013;4.3).</p>
<p>The field site (40&#xb0;16&#x2032;25.51&#x2033;N, 7&#xb0;25&#x2032;22.32&#x2033;W) covered an area of 0.25&#xa0;ha and was located at 440&#xa0;m.a.s.l. on a 2&#xb0; slope facing south. It was part of a vineyard until about 1985, when it was converted into a permanent pasture with grazing predominantly by cows. Since 2002, it has been managed as a sown biodiverse pasture a.k.a. SBP (see <xref ref-type="bibr" rid="B21">Teixeira et al., 2015</xref> for a detailed description). Antibiotics were used sparingly in the cattle; that is, on average, an estimated one cow per year (in a herd of &#x223c;100 animals) is treated with antibiotics, prescribed by the veterinarian, which are administered for 3&#xa0;days while keeping the cow in the stables, including for an extra week after administration.</p>
<p>In April 2022, the field site was completely fenced off, thereby excluding grazing animals. After approximately 6&#xa0;months, we noticed some rabbit droppings in a few places, so we rabbit-proofed the fence with a finer mesh to a height of 0.4&#xa0;m. Since then, we have not observed rabbit droppings again, with only mouse droppings found occasionally inside plots (estimated, but not recorded, to be &#x3c;1% cover). In October 2022, the fenced area was disked to 20&#xa0;cm depth, and 20 plots of 3&#xa0;m &#xd7; 3&#xa0;m were installed in a randomized block design, 10 control plots and 10 biochar-amended plots where biochar was mixed with topsoil (3% v/v) in a cement mixer to ensure a homogeneous distribution before being placed back into each plot (see <xref ref-type="bibr" rid="B13">Jongen et al., 2025</xref> for a full description of the installation procedure). Blocks were 7&#xa0;m apart, and plots within blocks were at least 3&#xa0;m apart. Each plot was instrumented with a microclimate probe (TOMST extreme) in October 2023, which recorded topsoil temperature and volumetric soil moisture content at &#x2212;6&#xa0;cm, both logged at 15-min intervals. During the wet part of the year&#x2014;October to April&#x2014;the topsoil temperature ranged from 2.5&#xb0;C to 24.3&#xb0;C. The median volumetric topsoil moisture content during this period was 37%. We did not record topsoil temperatures during the dry part of the year&#x2014;May to September&#x2014;but we did record temperatures in other plots at the same field site. The median topsoil temperature was 24&#xb0;C, while the maximum was 42&#xb0;C (the topsoil reached temperatures warmer than 35&#xb0;C on 43 days). The median volumetric topsoil moisture content during the summer period was 16%.</p>
</sec>
<sec id="s2-2">
<title>2.2 Soil sampling</title>
<p>During the approximate peak of biomass development of the SBP, on 24 April 2024, that is, 24&#xa0;months after grazing exclusion, we took six topsoil samples (0&#x2013;20&#xa0;cm) from each of the plots using an auger, taking two samples under each plant functional group; that is, grasses (patches with <italic>Lolium multiflorum</italic>, <italic>L. perenne</italic>, <italic>Vulpia bromoides</italic>, <italic>Agrostis</italic> sp.), legumes (patches with <italic>Trifolium subterraneum</italic>), and forbs (patches with <italic>Chamaemelum</italic> sp., <italic>Leontodon taraxacoides</italic>, <italic>Coleostephus myconis</italic>, and <italic>Bellardia trixago</italic>), and mixed them into one composite sample per plot (20 samples in total). In the biochar-amended plots, this corresponded to the biochar incorporation depth. Each composite sample was passed through a 2-mm sieve, homogenized again with a sterile spoon, and &#x223c;50&#xa0;g of the &#x3c;2&#xa0;mm subsample was transferred to a plastic zip lock bag, closed, and frozen at &#x2212;25&#xb0;C in a mobile field freezer. After each plot was sampled, the auger and sieve were cleaned with 70% ethanol before continuing to the next plot. We used a portable electricity-powered freezer to keep the samples at &#x2212;25&#xb0;C until arrival at the laboratory. Upon arrival, samples were frozen at &#x2212;75&#xb0;C. After 30&#xa0;days, samples were sent for DNA extraction and shotgun metagenomic analysis to Macrogen Inc., South Korea, in a box with four icepacks.</p>
</sec>
<sec id="s2-3">
<title>2.3 Shotgun metagenomic sequencing</title>
<p>The <italic>E. coli</italic> relative abundance (RA) data came as part of a metagenomic shotgun sequencing. The sequencing aimed at quantifying the relative abundances of bacteria and Archaea present in the topsoil during the peak biomass development of SBP, and therefore, it does not inform on their potential pathogenicity or virulence. The sequencing library was prepared by DNA fragmentation and ligation using ligate adapters. Adapter-ligated fragments were amplified by PCR. Clusters were generated using bridge amplification, and sequencing was performed with Illumina synthesis (SBS) technology. After sequencing was completed, Illumina HiSeq/NovaSeq raw data were demultiplexed by index sequences, and paired-end FASTQ files were generated for each sample. Adapter sequences and data with an average phred quality score less than 20 were removed using Trimmomatic (v0.39) (<xref ref-type="bibr" rid="B5">Bolger et al., 2014</xref>) of the Kneaddata (v0.10.0) (<xref ref-type="bibr" rid="B3">Beghini et al., 2021</xref>) pipeline (option: SLIDINGWINDOW:4:20). The pre-processed data were analyzed using MetaPhlAn4 (v4.0.0) (<xref ref-type="bibr" rid="B4">Blanco-Miguez et al., 2023</xref>) for &#x223c;1 million microorganisms composed of NCBI reference genomes and species-level genome bins (SGBs). The reads were mapped to the specific marker genes of the microbial species using Bowtie2 (<xref ref-type="bibr" rid="B14">Langmead and Salzberg, 2012</xref>), and the species abundances were calculated based on the average number of reads mapped to the marker genes. At this time, marker genes that mapped to less than 33% were removed as the species were considered absent.</p>
</sec>
<sec id="s2-4">
<title>2.4 Statistical analysis</title>
<p>A linear mixed-effects model (LME) with restricted maximum likelihood was used to analyze the response of <italic>E. coli</italic> relative abundance to treatments with blocks as a random factor, using the R packages lme4 and lmerTest. Estimated marginal means, along with their confidence intervals, were extracted from the LME model using the &#x201c;modelbased&#x201d; package and visualized as a violin plot using the &#x201c;ggplot2&#x201d; and &#x201c;see&#x201d; packages. Correlation graphs between <italic>E. coli</italic> RA, soil moisture, and temperature parameters were made using the &#x201c;ggpubr&#x201d; package.</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<p>
<italic>E. coli</italic> was present in composite topsoil samples from all 20 plots (<xref ref-type="sec" rid="s13">Supplementary Figure S1</xref>). The estimated marginal means were 17.0% and 20.2% for control and biochar-amended treatments, respectively, without any significant difference between treatments (<xref ref-type="fig" rid="F1">Figure 1a</xref>). No significant correlations were found between <italic>E. coli</italic> relative abundance (RA) and soil pH or soil organic matter contents before or after biochar application (<xref ref-type="fig" rid="F1">Figures 1b&#x2013;d</xref>) or metrics of soil temperature during the growing season (<xref ref-type="sec" rid="s13">Supplementary Figure S1</xref>). For biochar-amended plots, <italic>E. coli</italic> RA significantly increased with increasing soil moisture content, whereas there was no significant trend for the control plots. Potential indicators of <italic>E. coli</italic> inputs from animal feces did not show any significant correlations with <italic>E. coli</italic> RA (<xref ref-type="fig" rid="F1">Figures 1e,f</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>
<bold>(a)</bold> Violin plot of <italic>E. coli</italic> RA in topsoil samples. Yellow dots represent the estimated marginal means of its relative abundance in each treatment, along with confidence intervals, extracted from the linear mixed effects model; <bold>(b)</bold> scatterplot between soil pH measured in KCl and <italic>E. coli</italic> RA; <bold>(c)</bold> scatterplot between soil organic matter content and <italic>E. coli</italic> RA; <bold>(d)</bold> scatterplot between soil organic matter content in the plots before biochar application and <italic>E. coli</italic> RA; <bold>(e)</bold> scatterplot between total available N in the soil and <italic>E. coli</italic> RA; <bold>(f)</bold> scatterplot between plant available P (represented as P<sub>2</sub> from P<sub>2</sub>O<sub>5</sub>) and <italic>E. coli</italic> RA.</p>
</caption>
<graphic xlink:href="fenvs-13-1638383-g001.tif">
<alt-text content-type="machine-generated">Graphs analyzing the relative abundance of *Escherichia coli* depict various conditions: a violin plot shows control versus biochar treatment; scatter plots illustrate relationships with pH (b), soil organic matter (c and d), and plant nutrients nitrogen (e) and phosphorus (f). Data points are color-coded by treatment type.</alt-text>
</graphic>
</fig>
<p>One of the samples was a statistical outlier (2.7 standard deviations removed from the mean in the biochar-amended treatment). Within-treatment variation was low for the control soil (coefficient of variation (CV) of 0.15) and approximately double that for the biochar-amended soil (0.28), after omission of the outlier.</p>
</sec>
<sec id="s4">
<title>4 Discussion and recommendations</title>
<p>The prevalence of <italic>E. coli</italic> (100%) in our pasture topsoil samples is higher than in previous work, which did not use a &#x201c;without grazing&#x201d; exclusion, that found forest and pasture land use to be the most prevalent for the presence of <italic>E. coli</italic> with 73% prevalence, compared to 10% prevalence in cropland and 15% in shrubland (Dusek et al., 2018). The analysis used in this study can only report on relative abundance rather than population growth <italic>per se</italic>. Nevertheless, the results demonstrate that the <italic>E. coli</italic> population was still maintained at significant levels, as a proportion of the total microbial community investigated in this study over timeframes previously thought to be untenable without further inoculation from inputs such as manures or sewage sludges, especially in such a warm environment (<xref ref-type="bibr" rid="B22">Tran et al., 2020</xref>). The higher prevalence in our study may be because the previous study focused on the landscape scale, whereas our study was in a relatively homogeneous, small area, that is, 0.25&#xa0;ha. In addition, direct observations of runoff during an intense rainfall event, as well as observed soil surface evidence of erosion/deposition at the &#x3c;1&#xa0;m<sup>2</sup> scale throughout the field site following rainfall, suggest that bacteria may have been distributed from cow pats to other areas of the field site before grazing exclusion, thereby contributing to the observed <italic>E. coli</italic> prevalence. This is also supported by the lack of observed significant correlations between indicators of animal feces inputs (available N and P) and <italic>E. coli</italic> RA. <xref ref-type="bibr" rid="B23">Texier et al. (2008)</xref> found &#x3e;97% <italic>E. coli</italic> prevalence in Alpine grassland soils without grazing exclusion using the 16S rRNA method. Another grazing-exclusion field study, in Ireland, using lysimeters, found the prevalence of <italic>E. coli</italic> in leachate samples to vary from 3% to 48%, with the higher value being for a poor-draining Stagnic Luvisol (<xref ref-type="bibr" rid="B27">Brennan et al., 2010b</xref>). One possible reason may be that their use of the CFU method underestimated <italic>E. coli</italic> presence when it was in a viable but non-culturable (VBNC) state, to which our DNA sequencing method is not sensitive.</p>
<p>One of the biochar-amended treatment plots was a clear statistical extreme outlier for which we could not find any explanation in our auxiliary data. It is also unlikely to be caused by a point source, for example, wild animal excrement, because of the sampling design of six subsamples across each 9-m<sup>2</sup> plot. Potentially, it was caused by an auxiliary variable that we did not measure, or by a DNA analysis error. Apart from the extreme outlier, the remaining 19 plots showed relatively consistent <italic>E. coli</italic> RA, that is, 9.3%&#x2013;24.1%, with no significant difference between the control and biochar-amended treatments. Other biochar amendment studies have found increases, decreases, and no effect in <italic>E. coli</italic> survival depending on combinations of biochar and soil properties (<xref ref-type="bibr" rid="B25">Ye et al., 2018</xref>; <xref ref-type="bibr" rid="B1">Ahmadimoghadam et al., 2023</xref>; <xref ref-type="bibr" rid="B2">Alegbeleye and Sant&#x2019;Ana, 2023</xref>). Although the biochar treatment in our study had similar <italic>E. coli</italic> RA after 24&#xa0;months of grazing/manure exclusion as the control treatment, it did have double the range and almost double the within-treatment variation (after omitting the outlier). However, increases in range and within-treatment variation are not uncommon for studies where the treatment consists of incorporating a foreign material into the soil (<xref ref-type="bibr" rid="B19">Olmo et al., 2016</xref>), and a CV of 0.28 for a field study is still quite low.</p>
<p>Our findings do not conform with the reported environmental causes linked to <italic>E. coli</italic> decline, that is, high temperatures, aerobic conditions, low soil clay contents, and low DOC contents (<xref ref-type="bibr" rid="B24">van Elsas et al., 2011</xref>). If anything, two factors indicate reduced survivability. Specifically, clay content, which was 18% in the soil studied here, can be considered relatively low in this context, because soils with &#x3c;18% clay are classified as &#x201c;coarse&#x201d; (<xref ref-type="bibr" rid="B8">Feeney et al., 2024</xref>). Aerobic conditions can be considered relatively high because the soil is not compacted (exclusion of grazing animals), and the good drainage conditions cause a quick decrease in SMC during dry periods. Nevertheless, we still found evidence indicating that in this relatively aerobic soil, there was a trend for <italic>E. coli</italic> to increase with increasing soil moisture content during the growing season. We only found this trend for biochar-amended plots, which may be in part because of the twice-as-large range in <italic>E. coli</italic> RA compared to the control treatment. More field data, covering a wider range of soil moisture contents, are required to explore whether this trend is specific to biochar-amended soil or also occurs in control soil. Although we do not have data on DOC, the relatively low SOM content (before biochar amendment) would likely mean that DOC was also likely low, with this factor also likely to decrease <italic>E. coli</italic> abundance in the soil over time.</p>
<p>Due to the restricted inputs through exclusion of grazing by livestock and rabbits, the observed 17%&#x2013;20% relative abundance of <italic>E. coli</italic> after 24 months is suggestive of not only survival, but naturalization of <italic>E. coli</italic>. The field site has been under the same management since its conversion from a vineyard around 40&#xa0;years before this study. Antibiotic use in the grazing animals is so sparse that the potential of an antibiotic-resistant strain of <italic>E. coli</italic> entering the soil that way is extremely unlikely. Various bacterial adaptation mechanisms have been suggested recently to potentially favor <italic>E. coli</italic>, for example, stress resistance in nutrient-poor environments (Somorin et al., 2016), such as the soil environment included in this study, and other microbiological DNA or transcriptional mechanisms (<xref ref-type="bibr" rid="B20">Sazykin et al., 2024</xref>; <xref ref-type="bibr" rid="B17">Nakamoto et al., 2025</xref>), which are beyond the scope of this brief research report.</p>
<p>While our results demonstrate that <italic>E. coli</italic> can persist and even thrive relative to other microorganisms in rainfed pasture soils long after the removal of grazing animals, we acknowledge that our study did not include strain-level identification or functional gene profiling. As such, we cannot determine whether the <italic>E. coli</italic> populations observed harbor virulence factors or antibiotic resistance genes, which are critical for assessing their potential risks to public and environmental health. As such, any discussion of environmental and public health implications is, therefore, speculative at this point. Nevertheless, the long-term persistence of <italic>E. coli</italic> in soils raises important questions about its ecological adaptation and the possibility of these populations becoming naturalized components of the soil microbiome.</p>
<p>While the specific impacts are not fully understood, <italic>E coli</italic> naturalization in these pastoral ecosystems can have both on-farm and off-farm implications (<xref ref-type="bibr" rid="B12">Jones, 1999</xref>; <xref ref-type="bibr" rid="B15">Muirhead, 2023</xref>), which should be taken into consideration under a One Health approach. At the farm level, this could reflect in altered microbial and nutrient dynamics, thus impacting soil productivity and other soil-based ecosystem services (<xref ref-type="bibr" rid="B6">Cheong et al., 2024</xref>). Environmental health implications include potential soil/water contamination (<xref ref-type="bibr" rid="B12">Jones, 1999</xref>), splashing during rainfall events, as well as surface runoff into downstream freshwater systems and leaching into groundwater reservoirs (<xref ref-type="bibr" rid="B15">Muirhead, 2023</xref>). In addition to directly affecting human communities that rely on these resources, other direct human exposure routes are linked to occupational health and safety for individuals working in or near contaminated soil or water via ingestion, inhalation, or dermal contact, thus potentially leading to gastrointestinal illness if pathogenic strains are present. Concerns in relation to food safety may also arise upon <italic>E. coli</italic> entering the food chain, particularly in integrated crop-livestock systems (<xref ref-type="bibr" rid="B12">Jones, 1999</xref>; <xref ref-type="bibr" rid="B15">Muirhead, 2023</xref>; <xref ref-type="bibr" rid="B6">Cheong et al., 2024</xref>).</p>
<p>Specific recommendations for research and management involve filling current knowledge gaps in elucidating the potential long-term consequences to soil ecosystems, livestock, and public health, as well as the underlying mechanisms for <italic>E. coli</italic> survival and human/animal exposure. It is also necessary to further understand the ecological role of naturalized <italic>E. coli</italic> populations in soil, their interactions with native microbes, and the factors influencing their persistence and pathogenicity. Specifically, we recommend that future studies incorporate metagenomic, whole-genome sequencing, or functional gene assays to explore the phylogenetic diversity, pathogenic potential, and resistance profiles of soil-resident <italic>E. coli</italic> strains. Such approaches would significantly strengthen risk assessments and guide appropriate monitoring and management strategies.</p>
<p>The findings of this study support the tentative conclusions from two previous field studies that <italic>E. coli</italic> may become naturalized in the soil (<xref ref-type="bibr" rid="B26">Brennan et al. , 2010a</xref>; <xref ref-type="bibr" rid="B23">Texier et al., 2008</xref>). This is contrary to a wealth of laboratory and greenhouse studies that report weeks to months of <italic>E. coli</italic> survival in soil (<xref ref-type="bibr" rid="B24">van Elsas et al., 2011</xref>; <xref ref-type="bibr" rid="B16">Murphy et al., 2024</xref>). Considering the implications for OneHealth, it would seem prudent to further investigate <italic>E. coli</italic> presence and survival in field soils. Specifically, the hypothesis that <italic>E. coli</italic> survives for longer, or has an increased chance of establishing itself in the soil microbial community, in soils that are year-round predominantly wet, is questioned by the results from this study, in which the soils are dry for extended periods of time. Therefore, we recommend that <italic>E. coli</italic>, including generic and pathogenic strains, be monitored in soils of other Mediterranean field sites, under any land use, that receive organic fertilizers where <italic>E. coli</italic> may be present (manures, slurries, and low-grade treated sewage), and pastures with and without the presence of grazers up to 2&#xa0;years after the final amendment with those fertilizers. In addition, monitoring <italic>E. coli</italic> in the same samples with multiple measurement techniques, for example, culturing and various types of differential abundance methods (<xref ref-type="bibr" rid="B18">Nearing et al., 2022</xref>), may provide further insights into how these quantifications relate to each other.</p>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>We conclude that <italic>E. coli</italic> not only survived, but apparently thrived, either growing or maintaining its prevalence relative to the abundance of other microorganisms characterized in this study, in both the control and biochar-amended pasture soil, 24&#xa0;months after grazing/manure input exclusion. This has implications for monitoring and management requirements as well as policy implications, whereby currently the exclusion period for growing fresh produce globally varies between 90&#xa0;days and 365&#xa0;days following manure application. Future research incorporating genomic or functional profiling is key for determining whether soil-resident <italic>E. coli</italic> represent a benign, naturalized community or a latent risk to human and animal health.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The authors acknowledge that the data presented in this study must be deposited and made publicly available in an acceptable repository, prior to publication. Frontiers cannot accept a manuscript that does not adhere to our open data policies.</p>
</sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>FV: Formal Analysis, Writing &#x2013; original draft, Methodology, Writing &#x2013; review and editing, Conceptualization, Validation. AJ: Writing &#x2013; original draft, Methodology, Data curation, Conceptualization, Investigation, Visualization, Validation, Writing &#x2013; review and editing, Formal Analysis. SJ: Writing &#x2013; review and editing, Writing &#x2013; original draft, Supervision. TD: Writing &#x2013; review and editing, Writing &#x2013; original draft. ZK: Data curation, Writing &#x2013; original draft, Writing &#x2013; review and editing. AB: Resources, Writing &#x2013; original draft, Funding acquisition, Conceptualization, Methodology, Writing &#x2013; review and editing, Supervision.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research and/or publication of this article. This work is funded by national funds through FCT&#x2013;Funda&#xe7;&#xe3;o para a Ci&#xea;ncia e a Tecnologia I.P., under the project/grant UID/50006 &#x2b; LA/P/0094/2020 (<ext-link ext-link-type="uri" xlink:href="http://doi.org/10.54499/LA/P/0094/2020">doi.org/10.54499/LA/P/0094/2020</ext-link>). We acknowledge the support of projects POLLINATE (<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.54499/PTDC/EAM-AMB/1509/2021">https://doi.org/10.54499/PTDC/EAM-AMB/1509/2021</ext-link>) and HABITAT project 15341 (COMPETE2030- FEDER-00901600) through COMPETE 2030, FEDER, and FCT. We also acknowledge funding of authors FV (<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.54499/2023.06689.CEECIND/CP2840/CT0011">https://doi.org/10.54499/2023.06689.CEECIND/CP2840/CT0011</ext-link>) and A.C. Bastos through OE (<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.54499/DL57/2016/CP1482/CT0006">https://doi.org/10.54499/DL57/2016/CP1482/CT0006</ext-link>). We also wish to thank the European Commission Joint Program SOIL for the funding of the project TRUESOIL (<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.54499/EJPSoils/0001/2021">https://doi.org/10.54499/EJPSoils/0001/2021</ext-link>). La Caixa Foundation 100010434 is acknowledged for funding A. Jelin&#x10d;i&#x107; through the PhD fellowship Inphinit Incoming 2022 (PASTURE&#x2b;, LCF/BQ/DI22/11940011).</p>
</sec>
<ack>
<p>We wish to acknowledge Oscar Gonzalez-Pelayo and Victoria Arcenegui-Baldo for assistance with sampling.</p>
</ack>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>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.</p>
</sec>
<sec id="s10">
<title>Correction note</title>
<p>A correction has been made to this article. Details can be found at: <ext-link ext-link-type="uri" xlink:href="http://doi.org/10.3389/fenvs.2025.1708183">10.3389/fenvs.2025.1708183</ext-link>.</p>
</sec>
<sec sec-type="ai-statement" id="s11">
<title>Generative AI statement</title>
<p>The author(s) declare that no Generative AI was used in the creation of this manuscript.</p>
</sec>
<sec sec-type="disclaimer" id="s12">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec sec-type="supplementary-material" id="s13">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fenvs.2025.1638383/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fenvs.2025.1638383/full&#x23;supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Supplementaryfile1.docx" id="SM1" mimetype="application/docx" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ahmadimoghadam</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Tabatabaei</surname>
<given-names>S. H.</given-names>
</name>
<name>
<surname>Nourmahnad</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Hydrophobic biochar enhances Escherichia coli retention in a sandy loam soil influenced by irrigation management</article-title>. <source>Int. J. Environ. Sci. Technol.</source> <volume>20</volume> (<issue>8</issue>), <fpage>8399</fpage>&#x2013;<lpage>8410</lpage>. <pub-id pub-id-type="doi">10.1007/s13762-023-05002-8</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alegbeleye</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Sant&#x2019;Ana</surname>
<given-names>A. S.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Survival behavior of six enterotoxigenic <italic>Escherichia coli</italic> strains in soil and biochar-amended soils</article-title>. <source>Environ. Res.</source> <volume>223</volume>, <fpage>115443</fpage>. <pub-id pub-id-type="doi">10.1016/j.envres.2023.115443</pub-id>
<pub-id pub-id-type="pmid">36781014</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Beghini</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>McIver</surname>
<given-names>L. J.</given-names>
</name>
<name>
<surname>Blanco-M&#xed;guez</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Dubois</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Asnicar</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Maharjan</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2021</year>). <article-title>Integrating taxonomic, functional, and strain-level profiling of diverse microbial communities with bioBakery 3</article-title>. <source>Elife</source> <volume>10</volume>, <fpage>e65088</fpage>. <pub-id pub-id-type="doi">10.7554/elife.65088</pub-id>
<pub-id pub-id-type="pmid">33944776</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Blanco-Miguez</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Beghini</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Cumbo</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>McIver</surname>
<given-names>L. J.</given-names>
</name>
<name>
<surname>Thompson</surname>
<given-names>K. N.</given-names>
</name>
<name>
<surname>Zolfo</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Extending and improving metagenomic taxonomic profiling with uncharacterized species using MetaPhlAn 4</article-title>. <source>Nat. Biotechnol.</source> <volume>41</volume>, <fpage>1633</fpage>&#x2013;<lpage>1644</lpage>. <pub-id pub-id-type="doi">10.1038/s41587-023-01688-w</pub-id>
<pub-id pub-id-type="pmid">36823356</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bolger</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Marc</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Bjoern</surname>
<given-names>U.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Trimmomatic: a flexible trimmer for illumina sequence data</article-title>. <source>Bioinformatics</source> <volume>30</volume>, <fpage>2114</fpage>&#x2013;<lpage>2120</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btu170</pub-id>
<pub-id pub-id-type="pmid">24695404</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brennan</surname>
<given-names>F. P.</given-names>
</name>
<name>
<surname>Abram</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Chinalia</surname>
<given-names>F. A.</given-names>
</name>
<name>
<surname>Richards</surname>
<given-names>K. G.</given-names>
</name>
<name>
<surname>O&#x2019;Flaherty</surname>
<given-names>V.</given-names>
</name>
</person-group> (<year>2010a</year>). <article-title>Characterization of environmentally persistent Escherichia coli isolates leached from an Irish soil</article-title>. <source>Applied and environmental microbiology</source>, <volume>76</volume>(<issue>7</issue>), <fpage>2175</fpage>&#x2013;<lpage>2180</lpage>.<pub-id pub-id-type="pmid">20154122</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brennan</surname>
<given-names>F. P.</given-names>
</name>
<name>
<surname>O&#x2019;Flaherty</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Kramers</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Grant</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Richards</surname>
<given-names>K. G.</given-names>
</name>
</person-group> (<year>2010b</year>). <article-title>Long-term persistence and leaching of Escherichia coli in temperate maritime soils</article-title>. <source>Applied and environmental microbiology</source>, <volume>76</volume>(<issue>5</issue>), <fpage>1449</fpage>&#x2013;<lpage>1455</lpage>.<pub-id pub-id-type="pmid">20038692</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cheong</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Jay-Russell</surname>
<given-names>M. T.</given-names>
</name>
<name>
<surname>Chandler-Khayd</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Di</surname>
<given-names>F. J.</given-names>
</name>
<name>
<surname>Haghani</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Aminanadi</surname>
<given-names>P.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Presence of foodborne pathogens and survival of generic Escherichia coli in an organic integrated crop-livestock system</article-title>. <source>Front. Sustain. Food Syst.</source> <volume>8</volume>. <pub-id pub-id-type="doi">10.3389/fsufs.2024.1343101</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chung</surname>
<given-names>J. W.</given-names>
</name>
<name>
<surname>Edewi</surname>
<given-names>O. C.</given-names>
</name>
<name>
<surname>Foppen</surname>
<given-names>J. W.</given-names>
</name>
<name>
<surname>Gerner</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Krebs</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Lens</surname>
<given-names>P. N. L.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Removal of <italic>Escherichia coli</italic> by intermittent operation of saturated sand columns supplemented with hydrochar derived from sewage sludge</article-title>. <source>Appl. Sci.</source> <volume>7</volume> (<issue>8</issue>), <fpage>839</fpage>. <pub-id pub-id-type="doi">10.3390/app7080839</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="book">
<collab>Department for Environment, Food &#x0026; Rural Affairs</collab>. (<year>2018</year>). <source>Annual report and accounts 2017&#x2013;2018 <italic>(HC 1263)</italic>
</source>. <publisher-name>UK Government</publisher-name>. <comment>Available online at: <ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="https://www.gov.uk/government/publications/defras-annual-report-and-accounts-2017-to-2018">https://www.gov.uk/government/publications/defras-annual-report-and-accounts-2017-to-2018</ext-link>.</comment>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dusek</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Hewitt</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Schmidt</surname>
<given-names>K. N.</given-names>
</name>
<name>
<surname>Bergholz</surname>
<given-names>P. W.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Landscape-scale factors affecting the prevalence of Escherichia coli in surface soil include land cover type, edge interactions, and soil pH</article-title>. <source>Applied and environmental microbiology</source>, <volume>84</volume> (<issue>10</issue>), <fpage>e02714</fpage>&#x2013;<lpage>17</lpage>.<pub-id pub-id-type="pmid">29523546</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feeney</surname>
<given-names>C. J.</given-names>
</name>
<name>
<surname>Bentley</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>De Rosa</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Panagos</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Emmett</surname>
<given-names>B. A.</given-names>
</name>
<name>
<surname>Thomas</surname>
<given-names>A.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Benchmarking soil organic carbon (SOC) concentration provides more robust soil health assessment than the SOC/clay ratio at European scale</article-title>. <source>Sci. Total Environ.</source> <volume>951</volume>, <fpage>175642</fpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2024.175642</pub-id>
<pub-id pub-id-type="pmid">39163941</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gurtler</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Mullen</surname>
<given-names>C. A.</given-names>
</name>
<name>
<surname>Boateng</surname>
<given-names>A. A.</given-names>
</name>
<name>
<surname>Ma&#x161;ek</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Camp</surname>
<given-names>M. J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Biocidal activity of fast pyrolysis biochar against <italic>Escherichia coli</italic> O157: H7 in soil varies based on production temperature or age of biochar</article-title>. <source>J. Food Prot.</source> <volume>83</volume> (<issue>6</issue>), <fpage>1020</fpage>&#x2013;<lpage>1029</lpage>. <pub-id pub-id-type="doi">10.4315/0362-028x.jfp-19-331</pub-id>
<pub-id pub-id-type="pmid">32438397</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="book">
<collab>International Union of Soil Sciences (IUSS)</collab> (<year>2014</year>). <source>World reference base for soil resources 2014: international soil classification system for naming soils and creating legends for soil maps. world soil resources reports no. 106</source>. <publisher-loc>Rome</publisher-loc>: <publisher-name>FAO</publisher-name>.</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jones</surname>
<given-names>D. L.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Potential health risks associated with the persistence of Escherichia coli O157 in agricultural environments</article-title>. <source>Soil Use Manag.</source> <volume>15</volume> (<issue>2</issue>), <fpage>76</fpage>&#x2013;<lpage>83</lpage>. <pub-id pub-id-type="doi">10.1111/j.1475-2743.1999.tb00069.x</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Jongen</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bastos</surname>
<given-names>A. C.</given-names>
</name>
<name>
<surname>Jelin&#x10d;i&#x107;</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Sim&#xf5;es</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Gonzalez Pelayo</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Niemeyer</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2025</year>). <article-title>Biochar, plant productivity and community structure: a two-year study in a newly sown legume-rich pasture in Portugal under ambient and reduced precipitation</article-title>. <comment>Available online at: <ext-link ext-link-type="uri" xlink:href="https://ssrn.com/abstract=5215014">https://ssrn.com/abstract&#x3d;5215014</ext-link>. (Accessed August, 2025).</comment>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Langmead</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Salzberg</surname>
<given-names>S. L.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Fast gapped-read alignment with bowtie 2</article-title>. <source>Nat. Methods</source> <volume>9</volume> (<issue>4</issue>), <fpage>357</fpage>&#x2013;<lpage>359</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.1923</pub-id>
<pub-id pub-id-type="pmid">22388286</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Muirhead</surname>
<given-names>R. W.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Escherichia coli runoff from sheep and dairy cow grazed pasture: a plot scale simulation</article-title>. <source>J. Environ. Qual.</source> <volume>52</volume> (<issue>2</issue>), <fpage>287</fpage>&#x2013;<lpage>295</lpage>. <pub-id pub-id-type="doi">10.1002/jeq2.20448</pub-id>
<pub-id pub-id-type="pmid">36645734</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Murphy</surname>
<given-names>C. M.</given-names>
</name>
<name>
<surname>Weller</surname>
<given-names>D. L.</given-names>
</name>
<name>
<surname>Bardsley</surname>
<given-names>C. A.</given-names>
</name>
<name>
<surname>Ingram</surname>
<given-names>D. T.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Oryang</surname>
<given-names>D.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Survival of twelve pathogenic and generic Escherichia coli strains in agricultural soils as influenced by strain, soil type, irrigation regimen, and soil amendment</article-title>. <source>J. food Prot.</source> <volume>87</volume> (<issue>10</issue>), <fpage>100343</fpage>. <pub-id pub-id-type="doi">10.1016/j.jfp.2024.100343</pub-id>
<pub-id pub-id-type="pmid">39147099</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nakamoto</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Kobayashi</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Watanabe</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Kikuta</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Imamura</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Shimada</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Identification of a comprehensive set of transcriptional regulators involved in the long-term survivability of Escherichia coli in soil</article-title>. <source>Sci. Rep.</source> <volume>15</volume> (<issue>1</issue>), <fpage>4279</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-025-85609-8</pub-id>
<pub-id pub-id-type="pmid">39905026</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nearing</surname>
<given-names>J. T.</given-names>
</name>
<name>
<surname>Douglas</surname>
<given-names>G. M.</given-names>
</name>
<name>
<surname>Hayes</surname>
<given-names>M. G.</given-names>
</name>
<name>
<surname>MacDonald</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Desai</surname>
<given-names>D. K.</given-names>
</name>
<name>
<surname>Allward</surname>
<given-names>N.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Microbiome differential abundance methods produce different results across 38 datasets</article-title>. <source>Nat. Commun.</source> <volume>13</volume> (<issue>1</issue>), <fpage>342</fpage>. <pub-id pub-id-type="doi">10.1038/s41467-022-28034-z</pub-id>
<pub-id pub-id-type="pmid">35039521</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Olmo</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Lozano</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Barr&#xf3;n</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Villar</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Spatial heterogeneity of soil biochar content affects soil quality and wheat growth and yield</article-title>. <source>Sci. Total Environ.</source> <volume>562</volume>, <fpage>690</fpage>&#x2013;<lpage>700</lpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2016.04.089</pub-id>
<pub-id pub-id-type="pmid">27110980</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sazykin</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Naumova</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Azhogina</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Klimova</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Karchava</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Khmelevtsova</surname>
<given-names>L.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Glyphosate effect on biofilms formation, mutagenesis and stress response of E. &#x441;oli</article-title>. <source>J. Hazard. Mater.</source> <volume>461</volume>, <fpage>132574</fpage>. <pub-id pub-id-type="doi">10.1016/j.jhazmat.2023.132574</pub-id>
<pub-id pub-id-type="pmid">37748310</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Somorin</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Abram</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Brennan</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>O&#x2019;Byrne</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The general stress response is conserved in long-term soil-persistent strains of Escherichia coli</article-title>. <source>Applied and Environmental Microbiology</source>, <volume>82</volume> (<issue>15</issue>), <fpage>4628</fpage>&#x2013;<lpage>4640</lpage>.<pub-id pub-id-type="pmid">27235429</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Texier</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Prigent-Combaret</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Gourdon</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Poirier</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Faivre</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Dorioz</surname>
<given-names>J. M.</given-names>
</name>
<etal/>
</person-group> (<year>2008</year>). <article-title>Persistence of culturable Escherichia coli fecal contaminants in dairy alpine grassland soils</article-title>. <source>Journal of Environmental Quality</source>, <volume>37</volume>(<issue>6</issue>), <fpage>2299</fpage>&#x2013;<lpage>2310</lpage>.<pub-id pub-id-type="pmid">18948484</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Teixeira</surname>
<given-names>R. F.</given-names>
</name>
<name>
<surname>Proenca</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Crespo</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Valada</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Domingos</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>A conceptual framework for the analysis of engineered biodiverse pastures</article-title>. <source>Ecol. Eng.</source> <volume>77</volume>, <fpage>85</fpage>&#x2013;<lpage>97</lpage>. <pub-id pub-id-type="doi">10.1016/j.ecoleng.2015.01.002</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tran</surname>
<given-names>D. T.</given-names>
</name>
<name>
<surname>Bradbury</surname>
<given-names>M. I.</given-names>
</name>
<name>
<surname>Van Ogtrop</surname>
<given-names>F. F.</given-names>
</name>
<name>
<surname>Bozkurt</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Jones</surname>
<given-names>B. J.</given-names>
</name>
<name>
<surname>McConchie</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Environmental drivers for persistence of Escherichia coli and salmonella in manure-amended soils: a meta-analysis</article-title>. <source>J. food Prot.</source> <volume>83</volume> (<issue>7</issue>), <fpage>1268</fpage>&#x2013;<lpage>1277</lpage>. <pub-id pub-id-type="doi">10.4315/0362-028x.jfp-19-460</pub-id>
<pub-id pub-id-type="pmid">32577760</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Van Elsas</surname>
<given-names>J. D.</given-names>
</name>
<name>
<surname>Semenov</surname>
<given-names>A. V.</given-names>
</name>
<name>
<surname>Costa</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Trevors</surname>
<given-names>J. T.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Survival of Escherichia coli in the environment: fundamental and public health aspects</article-title>. <source>ISME J.</source> <volume>5</volume> (<issue>2</issue>), <fpage>173</fpage>&#x2013;<lpage>183</lpage>. <pub-id pub-id-type="doi">10.1038/ismej.2010.80</pub-id>
<pub-id pub-id-type="pmid">20574458</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ye</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Jiao</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Targeted inactivation of antibiotic-resistant Escherichia coli and Pseudomonas aeruginosa in a soil-lettuce system by combined polyvalent bacteriophage and biochar treatment</article-title>. <source>Environ. Pollut.</source> <volume>241</volume>, <fpage>978</fpage>&#x2013;<lpage>987</lpage>. <pub-id pub-id-type="doi">10.1016/j.envpol.2018.04.070</pub-id>
<pub-id pub-id-type="pmid">30029332</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>