<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2023.1173023</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Plant Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Assessment and comparison of rhizosphere communities in cultivated <italic>Vaccinium</italic> spp. provide a baseline for study of causative agents in decline</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Kawash</surname><given-names>Joseph</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2274595"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Oudemans</surname><given-names>Peter V.</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1809614"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Erndwein</surname><given-names>Lindsay</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Polashock</surname><given-names>James J.</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>*</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/89632"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>USDA-ARS, Genetic Improvement of Fruits and Vegetables Laboratory</institution>, <addr-line>Chatsworth, NJ</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Plant Biology, Rutgers University</institution>, <addr-line>New Brunswick, NJ</addr-line>, <country>United States</country></aff>
<aff id="aff3"><sup>3</sup><institution>Oak Ridge Institute for Science and Education (ORISE) Postdoctoral Scholar, USDA-ARS, Genetic Improvement of Fruits and Vegetables Laboratory</institution>, <addr-line>Chatsworth, NJ</addr-line>, <country>United States</country></aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Irena Mavri&#x10d; Plesko, Agricultural Institute of Slovenia, Slovenia</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Rishi Ram Burlakoti, Agriculture and Agri-Food Canada (AAFC), Canada; Mura Jyostna Devi, USDA, United States</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: James J. Polashock, <email xlink:href="mailto:james.polashock@usda.gov">james.polashock@usda.gov</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>27</day>
<month>06</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1173023</elocation-id>
<history>
<date date-type="received">
<day>24</day>
<month>02</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>29</day>
<month>05</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Kawash, Oudemans, Erndwein and Polashock</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Kawash, Oudemans, Erndwein and Polashock</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>It has long been recognized that the community of organisms associated with plant roots is a critical component of the phytobiome and can directly or indirectly contribute to the overall health of the plant. The rhizosphere microbial community is influenced by a number of factors including the soil type, the species of plants growing in those soils, and in the case of cultivated plants, the management practices associated with crop production. <italic>Vaccinium</italic> species, such as highbush blueberry and American cranberry, are woody perennials that grow in sandy, acidic soils with low to moderate levels of organic matter and a paucity of nutrients. When properly maintained, fields planted with these crops remain productive for many years. In some cases, however, yields and fruit quality decline over time, and it is suspected that degenerating soil health and/or changes in the rhizosphere microbiome are contributing factors. Determining the assemblage of bacterial and fungal microorganisms typically associated with the rhizosphere of these crops is a critical first step toward addressing the complex issue of soil health. We hypothesized that since blueberry and cranberry are in the same genus and grow in similar soils, that their associated rhizosphere microbial communities would be similar to each other. We analyzed the eukaryotic (primarily fungal) and bacterial communities from the rhizosphere of representative blueberry and cranberry plants growing in commercial fields in New Jersey. The data presented herein show that while the bacterial communities between the crops is very similar, the fungal communities associated with each crop are quite different. These results provide a framework for examining microbial components that might contribute to the health of <italic>Vaccinium</italic> spp. crops in New Jersey and other parts of the northeastern U.S.</p>
</abstract>
<kwd-group>
<kwd><italic>Vaccinium corymbosum</italic>
</kwd>
<kwd><italic>Vaccinium macrocarpon</italic>
</kwd>
<kwd>phytobiome</kwd>
<kwd>acidic soil</kwd>
<kwd>soil health</kwd>
</kwd-group>
<counts>
<fig-count count="8"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="51"/>
<page-count count="11"/>
<word-count count="4919"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Plant Pathogen Interactions</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<title>Introduction</title>
<p>Highbush Blueberry (<italic>Vaccinium corymbosum</italic> L.) and large-fruited or American cranberry (<italic>V. macrocarpon</italic> Ait.) are long-lived woody perennials that are native to North America. Cultivars of both species are clonally propagated. Blueberry is propagated from hardwood or softwood cuttings or <italic>in vitro</italic> culture and rooted prior to field planting (<xref ref-type="bibr" rid="B34">Mainland, 1966</xref>). Cranberry cuttings can be pressed directly into field soil or rooted prior to transplanting into the field (<xref ref-type="bibr" rid="B19">Eck, 1990</xref>). Establishment time can vary depending on planting location, but in the northeastern U.S., it typically takes blueberry fields 4-7 years to reach maturity. New plantings of cranberry similarly take about 4-5 years to reach full production. Once established, plantings may remain productive for many years. However, the yield of mature plantings of blueberry and cranberry sometimes declines over time. In areas where decline is prevalent, and replanting does not solve the problem, the primary cause is likely to be soil-related.</p>
<p>It is well established that long-term monoculture can reduce soil health (<xref ref-type="bibr" rid="B32">Ketcheson, 1980</xref>; <xref ref-type="bibr" rid="B41">Porter et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B24">Garside et&#xa0;al., 2001</xref>; <xref ref-type="bibr" rid="B5">Bennett et&#xa0;al., 2012</xref>). This decline of soil health can be complex and includes, but is not limited to, a buildup of salts and/or other toxic materials, reduction in organic matter, and a buildup of detrimental soil organisms (<xref ref-type="bibr" rid="B45">Tewoldemedhin et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B37">Mazzola and Manici, 2012</xref>). Poor performance of replanted blueberry fields, termed Blueberry replant disease (BRD) in the southeastern U.S., was found to be primarily due to pathogenic nematodes (<xref ref-type="bibr" rid="B29">Jagdale et&#xa0;al., 2013</xref>), while cranberry decline in Wisconsin was associated with <italic>Phytophthora</italic> spp. (<xref ref-type="bibr" rid="B30">Jeffers, 1988</xref>). In contrast, blueberry and cranberry decline in New Jersey, does not seem to be associated with either <italic>Phytophthora</italic> spp. or nematodes, although in some replanted blueberry fields where experimental nematicides were tested, plant growth improved (Oudemans, unpublished). Thus, it is likely that the cause(s) of decline in either crop varies by locality and may be due to complex community interactions, rather than a single pathogenic organism.</p>
<p>Blueberry and cranberry naturally grow in well-drained, sandy acidic (pH 4.0-5.5) soils that tend to be low in nutrients and with low to moderate (2-7 percent) organic matter (<xref ref-type="bibr" rid="B20">Eck and Childers, 1966</xref>; <xref ref-type="bibr" rid="B18">Eck, 1988</xref>; <xref ref-type="bibr" rid="B19">Eck, 1990</xref>, Pavlis et&#xa0;al., unpublished). Since the soils in which these crops grow is not typical, as compared to the soils in which most crops are grown, the assemblage of microorganisms in these soils is likely to be unique. While we expect to find a range of acidophilic organisms, <xref ref-type="bibr" rid="B23">Fierer (2017)</xref> suggested that abiotic stressors, including acidic conditions, can limit microbial survival and growth. He further notes that soil bacterial communities and abundances of specific taxa can be predicted from site characteristics which include soil pH (<xref ref-type="bibr" rid="B23">Fierer, 2017</xref>).</p>
<p>In addition to soil type, the plants growing in the soil have a profound effect on the rhizosphere microbiome (<xref ref-type="bibr" rid="B13">Costa et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B6">Berendsen et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B36">Marschner, 2012</xref>; <xref ref-type="bibr" rid="B11">Chaparro et&#xa0;al., 2014</xref>). We expected that because highbush blueberry and cranberry are in the same genus (<italic>Vaccinium</italic>) and grow naturally in similar soils and environments, that there would be significant commonalities in the rhizosphere communities associated with these crops. However, plant habit (1-2&#xa0;m high shrub for blueberry; low-growing 0.2 -0.6&#xa0;m high vines for cranberry) and management practices differ between these crops and these differences might drive some distinctions in the rhizosphere communities of these crops.</p>
<p>Management practices differ greatly between blueberry and cranberry. Highbush blueberry in the northeastern U.S. is still grown primarily in clean culture, i.e. all competing vegetation is controlled by tilling and herbicide use (<xref ref-type="bibr" rid="B3">Ballinger, 1966</xref>). The plants are pruned to remove old wood and to encourage the growth of 2-3 year old productive shoots. The plants are deciduous and go dormant for the winter. In contrast, cranberry is grown in solid stands of vines in beds of varying size with little or no exposed soil. Weeds are pulled by hand or controlled with herbicides to form a near-monoculture of cranberry vines. Cranberry plants retain most of their leaves at the end of the growing season, but still go dormant for the winter.</p>
<p>In this study, we sought to determine and compare the rhizosphere microbiome (bacteria and fungi) of cultivated cranberry and blueberry fields in New Jersey as a prerequisite to exploring the cause(s) of decline. While we predicted commonalities among the microbial communities based on similarity of soil type and plant genus, management practices and plant growth habit were expected to cause a divergence in at least some aspects of the microbial communities.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<title>Materials and methods</title>
<sec id="s2_1">
<title>Plant material and collection</title>
<p>Blueberry rhizosphere soil samples were collected from 16 mature (over 10 years old) fields (2 samples/field) across 4 different commercial farms (3 in Atlantic County, NJ, USA and 1 in Burlington County, NJ, USA), for a total of 32 samples. An additional paired blueberry soil sample was collected from a 10-year old research field at the P.E. Marucci Blueberry and Cranberry Research Station (Burlington County, NJ, USA) (<xref ref-type="supplementary-material" rid="ST1"><bold>Table S1</bold></xref>). Rhizosphere soil was collected, with a small shovel, from the root zone (between 30 and 40&#xa0;cm from the crown) of 3-4 mature blueberry plants per sample location to a depth of 16&#xa0;cm and homogenized in buckets. The homogenized soil in each bucket was considered to be one sample. Cranberry soils from mature beds (over 10 years old) were collected from 8 fields (4 samples/field) from a commercial farm (E, <xref ref-type="supplementary-material" rid="ST1"><bold>Table S1</bold></xref>) in Burlington County, NJ, USA (32 samples). One extra paired sample was collected from commercial farm F (<xref ref-type="supplementary-material" rid="ST1"><bold>Table S1</bold></xref>). The 34 cranberry soil samples were collected with a soil probe to a depth of 10&#xa0;cm. The core, containing the roots of 3-4 samples per field were homogenized and the roots removed before DNA isolation. All blueberry and cranberry samples were collected on 22 or 23 July, 2015.</p>
</sec>
<sec id="s2_2">
<title>DNA extraction and amplification</title>
<p>Total DNA was extracted from the soil samples using the PowerLyzer PowerSoil DNA Isolation kit (MoBio Laboratories, Carlsbad, CA, USA), but with protocol modifications suggested by the manufacturer for soils with low biomass and based on our own experience. Briefly, 0.5&#xa0;g of each soil sample was added to the dry glass bead tubes, followed by 500 &#xb5;l Bead Solution, 200 &#xb5;l of phenol:chloroform (1:1) pH 8.0, and 60 &#xb5;l of C1 solution. Samples were processed in a TissueLyser II (Qiagen, Germantown, MD, USA) at a frequency of 20&#xa0;Hz for 2 x 5&#xa0;min and the slurry was transferred to Phase Lock Gel Light 2&#xa0;ml tubes (5 Prime, Hilden, Germany). The tubes were centrifuged at 10,000 x g for 1&#xa0;min and the supernatants were transferred to clean 2&#xa0;ml tubes. The remainder of the protocol was as per the manufacturer&#x2019;s standard protocol.</p>
<p>The DNA samples were processed for sequencing by Molecular Research LP (Shallowater, TX, USA). Diagnostic portions of the DNA were amplified using PCR. Specifically, the 16S rRNA gene V4 variable region primers 515/806 (<xref ref-type="bibr" rid="B10">Caporaso et&#xa0;al., 2011</xref>) and the rDNA ITS region primers ITS1/4 (<xref ref-type="bibr" rid="B49">White et&#xa0;al., 1990</xref>) were used to assay bacteria and fungi respectively.</p>
</sec>
<sec id="s2_3">
<title>DNA sequencing and computational processing</title>
<p>All amplified DNA fragments were sequenced using the Illumina MiSeq (Illumina, San Diego, CA, USA) platform. Operational taxonomic unit (OTU) picking was performed by Molecular Research LP (Shallowater, TX, USA) using an in-house pipeline developed by the provider (Molecular Research) from protocols described in; (<xref ref-type="bibr" rid="B17">Dowd et&#xa0;al., 2008a</xref>; <xref ref-type="bibr" rid="B16">Dowd et&#xa0;al., 2008b</xref>; <xref ref-type="bibr" rid="B21">Edgar, 2010</xref>; <xref ref-type="bibr" rid="B8">Capone et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B22">Eren et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B44">Swanson et&#xa0;al., 2011</xref>). Reads were separated by sample barcodes, followed by removal of sequences not meeting quality criteria, including short sequences &lt; 200bp, sequences with ambiguous base calls, and sequences with homopolymer runs exceeding 6bp. Singletons and chimeras were also removed. The reads were then separated into bins and the binned sequences were clustered into OTUs at 97% similarity. Final OTUs were taxonomically classified from a database derived from GreenGenes (<xref ref-type="bibr" rid="B15">DeSantis et&#xa0;al., 2006</xref>), RDPII, (<ext-link ext-link-type="uri" xlink:href="http://rdp.cme.msu.edu">http://rdp.cme.msu.edu</ext-link>), and NCBI (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov">www.ncbi.nlm.nih.gov</ext-link>) using BLASTn. 16S and 18S OTU tables for lowbush blueberry in Nova Scotia, Canada were kindly provided by S. Yurgel (<xref ref-type="bibr" rid="B50">Yurgel et&#xa0;al., 2017</xref>). 18S OTU tables for cultivated grape were provided by M. Chou (<xref ref-type="bibr" rid="B12">Chou et&#xa0;al., 2018</xref>).</p>
</sec>
<sec id="s2_4">
<title>Statistical analysis and graph generation</title>
<p>Alpha diversity was calculated using the QIIME analysis toolkit (<xref ref-type="bibr" rid="B9">Caporaso et&#xa0;al., 2010</xref>), specifically, the alpha_diversity.py script using the Shannon index. OTU tables for bacteria and fungi were processed individually. Analysis of variance (ANOVA) was used to compare alpha diversity between blueberry and cranberry for each sample group.</p>
<p>Differential abundance and statistical significance (assigned at 0.05) of OTUs between blueberry and cranberry were calculated using the Analysis of Compositions of Microbiomes (ANCOM) R package utilizing the provided stringent correction for multiple testing (<xref ref-type="bibr" rid="B35">Mandal et&#xa0;al., 2015</xref>). OTU tables for the bacteria and fungi were analyzed independently. Subsequent boxplots of differentially abundant OTUs were generated using a python script with NumPy and Matplotlib modules (<xref ref-type="bibr" rid="B46">Van der Walt et&#xa0;al., 2011</xref>).</p>
<p>Assignment of guilds for fungi was performed using a python script to query the FUNGuild database (<xref ref-type="bibr" rid="B39">Nguyen et&#xa0;al., 2016</xref>). OTU tables were OTU tables were parsed using python, specifying a fungal database search. All guild assignments were filtered as being plant-relevant unless otherwise noted.</p>
<p>Diversity graphs were generated with a python script utilizing the NetworkX package (<xref ref-type="bibr" rid="B26">Hagberg et&#xa0;al., 2008</xref>) and Matplotlib (<xref ref-type="bibr" rid="B28">Hunter, 2007</xref>) modules for creating the display. Each entry on the taxonomic tree starting with an individual phylum present in the collective OTU set is represented by a node, with nodes linked through descent of the taxonomic hierarchy to the family rank. Hence, following a given phylum are all of the classes that are associated with that phylum present in the OTU table. Unique OTUs found within a given taxon were counted and attributed to node area. Therefore, nodes with increased diversity will be larger, having several connecting nodes with a variety of lower taxa associated, or with few connecting edges, indicating a large number of less diverse taxa. Shading of the node ranged from blue (r,g,b:0,0,1) for blueberry and red (r,g,b:1,0,0) for cranberry. The corresponding node color was therefore a ratio of the number of OTUs found in cranberry versus blueberry (cranberry_OTUs/total_OTUs,0, blueberry_OTUs/total_OTUs). Taxa that are more diverse in cranberry samples are shaded increasingly red, while taxa more diverse in blueberry samples are shaded blue. The arrangement of nodes for display was established by the NetworkX fruchterman_reingold _layout method.</p>
<p>Circular graphs of the percent abundance of genera were generated using GraPhiAn (<xref ref-type="bibr" rid="B1">Asnicar et&#xa0;al., 2015</xref>). The imputed taxonomic trees were generated using the hierarchy of the OTU assignments from each OTU table (bacterial or fungal). The datasets were filtered to display only OTUs from the respective kingdom that were classifiable to the genus level. Corresponding graph values for each OTU represent the percent (as a number of total reads assigned to that particular OTU, out of all assigned reads) abundance of a given taxonomic level. Abundance equal to or exceeding 0.5% was annotated for graphical representation of contributing taxa.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<title>Results</title>
<sec id="s3_1">
<title>Fungi</title>
<p>The fungal analysis resulted in 1.26M reads matched to OTUs for blueberry and 1.35M reads for cranberry, averaging 38.2K reads and 39.7K reads per sample respectively. Both groups were dominated by the phyla Ascomycota, 63.0% and 88.6%, and Basidiomycota, 20.0% and 6.9%, in blueberry and cranberry respectively (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1A</bold></xref>). The percent abundance of major fungal phyla varied between blueberry and cranberry, and included; Chytridiomycota, Mucoromycota, Cryptomycota, Etomophthoromycota, and Glomeromycota (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1A</bold></xref>). However, no significant differential abundance between blueberry and cranberry was identified by ANCOM for any of these phyla.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Percent abundance of fungal <bold>(A)</bold> and bacterial <bold>(B)</bold> phyla in the rhizosphere of cultivated highbush blueberry and cranberry in New Jersey (N.J.). Data for lowbush blueberry from Nova Scotia (<xref ref-type="bibr" rid="B50">Yurgel et&#xa0;al., 2017</xref>) shown for comparison.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1173023-g001.tif"/>
</fig>
<p>There were 753 fungal genera identified using these methods, but the relative abundance of the majority (722) was below 1% (<xref ref-type="supplementary-material" rid="SF1"><bold>Figures S1</bold></xref>, <xref ref-type="supplementary-material" rid="SF2"><bold>S2</bold></xref>). The alpha diversity, based on the Shannon Index, was found to be significantly different between blueberry and cranberry (<xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2</bold></xref>). Most of the fungi that were classified using the FUNGuild database, at the genus level, with an abundance of 0.5% or greater, were saprotrophs (22) followed by symbionts (16), plant pathogens (8), endophytes (7), and epiphytes (3) (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3</bold></xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Alpha diversity (Shannon index) in the rhizosphere of cultivated highbush blueberry and cranberry in New Jersey (N.J.) for fungi and bacteria sampling. Horizontal line within each box represents the mean. Whiskers extend to the last data points that fall within 1.5 times the interquartile range from either the first or third quartile. Asterisk indicates significant difference between samples (<italic>P</italic> &lt; 0.001).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1173023-g002.tif"/>
</fig>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Boxplots of relative abundance for classified fungi in highbush blueberry (blue) and cranberry (red) in New Jersey. The vertical line within each box represents the median value. The shaded area extends from the first through the third quartiles. Whiskers extend to the last data points that fall within 1.5 times the interquartile range from either the first or third quartile. All data beyond these bounds are deemed outliers and are represented by individual circles.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1173023-g003.tif"/>
</fig>
<p>Within the shared (blueberry and cranberry) rhizosphere fungi, significant differential abundance was found in the pathogenic genera <italic>Fusarium</italic> (7.0% and 0.0%), <italic>Phoma</italic> (1.0% and 0.0%), <italic>Rhizophydium</italic> (0.5% and 0.7%), <italic>Coniochaeta (3.6% and 0.0)</italic>, and <italic>Mollisia</italic> (3.2% and 0.0%) (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4</bold></xref>). Blueberry also contained several fungal pathogens nearly uniquely (extremely low evidence, below 0.1% in cranberry) in its rhizosphere microbiome, though each represented less than 1% of the genus content. These pathogens included <italic>Athelia</italic> (0.7%), <italic>Alternaria</italic> (0.4%) and <italic>Cladosporium</italic> (0.3%).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Differentially abundant fungal genera with mean read counts greater than 200 in the rhizosphere of cultivated highbush blueberry (blue) and cranberry (red) in New Jersey, according to ANCOM. The horizontal line within each box represents the median value, the shaded area extends from the first through the third quartiles. Whiskers extend to the last points that fall within 1.5 times the interquartile range from either the first or third quartile. All data beyond these bounds are deemed outliers and are represented by individual circles.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1173023-g004.tif"/>
</fig>
<p>Of the symbionts in the blueberry and cranberry soils (16 genera were 0.5% or greater), 4 (<italic>Phialophora</italic> 0.1% and 0.5%, <italic>Rhizoscyphus</italic> 0.2% and 32.4%<italic>, Oidiodendron</italic> 0.0% and 0.1%, and <italic>Hymenoscyphus</italic> 0.0% and 0.5%<italic>)</italic> are known ericoid mycorrhizae. Nine other genera of mycorrhizal (ecto and arbuscular) fungi were also detected at abundances of 0.5% or greater, that are not generally thought to be associated with ericoid plants. Differentially abundant non-ericoid mycorrhizal fungi included (blueberry and cranberry), <italic>Ambispora</italic> (1.2% and 0.0%), <italic>Peziza</italic> (2.1% and 0.0%), <italic>Russula</italic> (5.3% and 0.1%), <italic>Sphaerosporella</italic> (0.1% and 24.8%), (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4</bold></xref>). <italic>Sarcosphaera</italic> (0.4% and 0.1%) was found most abundant in blueberry soil, but not significantly so. Endophyte and epiphyte classified genera were also detected with significant differential abundance including <italic>Cryptococcus</italic> (13.8% and 5.6%), and <italic>Mollisia</italic> (3.2% and 0.0%) (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4</bold></xref>).</p>
<p>We also compared differential abundance between both blueberry and cranberry soils with grape, taking note in particular of significant pathogenic and symbiotic fungi from cranberry and blueberry soils. Profiles of grape symbiotic fungal abundance were also relatively low and did not trend towards either blueberry or cranberry. The grape profile was similar to cranberry, for <italic>Russula</italic> and similar to blueberry, for <italic>Oidiodendron</italic> and <italic>Hymenoscyphus</italic> (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>). Of the pathogenic fungi, grape shared a similar <italic>Fusarium</italic> profile with blueberry, having a significantly high abundance compared to cranberry. However, the grape profile for the pathogenic fungus <italic>Phoma</italic>, had little similarity with blueberry, and greater resemblance to cranberry, although all three groups were significantly different from one another (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>). <italic>Rusulla</italic>, a genus containing ectomycorrhizal mushrooms was found to significantly differ across all three crops (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5C</bold></xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Differentially abundant fungal genera as a percentage of reads identified in the rhizosphere of cultivated highbush blueberry (left) and cranberry (center) in New Jersey, and cultivated grape (right) in New York (<xref ref-type="bibr" rid="B12">Chou et&#xa0;al., 2018</xref>). Fungal genera are categorized by: ericoid mycorrhizae <bold>(A)</bold>, pathogens <bold>(B)</bold>, and ectomycorrhizae <bold>(C)</bold>. The horizontal line within each box represents the median value, the boxed area extends from the first through the third quartiles. Whiskers extend to the last points that fall within 1.5 times the interquartile range from either the first or third quartile. All data beyond these bounds are deemed outliers and are represented by individual circles. Asterisk indicates significant difference between samples (<italic>P</italic> &lt; 0.05).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1173023-g005.tif"/>
</fig>
</sec>
<sec id="s3_2">
<title>Bacteria</title>
<p>In the bacterial analysis, 2.38M reads were matched with OTUs in blueberry, averaging 72.1K reads per sample, and 2.31M reads in cranberry, averaging 67.9K reads per sample. There was little difference in the makeup of the bacterial microbiome at the phylum level for blueberry vs. cranberry (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1B</bold></xref>). The order of dominance, however, differed between the two (blueberry vs. cranberry), though Proteobacteria (42.4% and 50.0%) and Acidobacteria (19.3% and 23.1%) were the most prevalent for both N.J. crops. High numbers of bacteria, from 5 other phyla [Planctomyces (6.2% and 5.4%), Verrucomicobia (1.6% and 1.9%), Bacteriodetes (4.4% and 3.5%), Firmicutes (7.6% and 5.8%), and Chloroflexi (5.2% and 1.9%)] were also detected. Any unique phyla found in either group represented less than 1% abundance in the rhizosphere microbiome. Even at the genus level, unique taxa were largely identified at an abundance below 1% (<xref ref-type="supplementary-material" rid="SF3"><bold>Figures S3</bold></xref>, <xref ref-type="supplementary-material" rid="SF4"><bold>S4</bold></xref>). Alpha diversity was not significantly different between blueberry and cranberry for bacteria OTUs (<xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2</bold></xref>).</p>
<p>Though no known bacterial pathogens were identified at the genus level as being significantly differentially present between blueberry and cranberry at any reasonable abundance, <italic>Pseudomonas</italic> (0.2% and 0.1%), <italic>Burkholderia</italic> (0.4% and 0.3%), <italic>Candidatus</italic> (genus 2.9% and 1.9%), <italic>Saccharimonas</italic> (0.2% and 0.3%), and <italic>Rhizobium</italic> (2.1% and 1.3%) were all found in the rhizosphere microbiome.</p>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>Discussion</title>
<p>The fungi were remarkably similar at the phylum level, as there were no notable omissions between the rhizosphere microbiome of blueberry or cranberry. They only differed in the hierarchical abundance within phyla groups (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1A</bold></xref>), and not significantly so. Although blueberry had a significantly higher alpha diversity compared to cranberry, potentially due to distinguishing agricultural practices and possibly plant growth habit. Cranberry is typically cultivated in wetlands and this may explain why the chytridiomycota, a phylum common in aquatic environments (<xref ref-type="bibr" rid="B4">Barr, 1990</xref>; <xref ref-type="bibr" rid="B25">Gleason et&#xa0;al., 2008</xref>), is prevalent in the cranberry soils. In fact, the chytrid <italic>Synchytrium vaccinii</italic> is a known pathogen of cranberry and requires free water to spread (<xref ref-type="bibr" rid="B43">Stretch and Carris, 2017</xref>). However, fungi from this phylum were not exclusive to cranberry as they were also found in both the New Jersey (N.J.) and the Nova Scotia (N.S.) blueberry soils (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1A</bold></xref>). These findings parallel other reports of <italic>Vaccinium</italic> spp. soils, especially the identification of ascomycota and basidiomycota as major composition groups, though not with absolute ranking similarities (<xref ref-type="bibr" rid="B33">Li et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B38">Morvan et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B51">Zhou et&#xa0;al., 2022</xref>).</p>
<p>Many fungal genera were classified that have not been previously recognized as inhabiting the soils associated with <italic>Vaccinium</italic> spp. in the northeastern U.S. (<xref ref-type="supplementary-material" rid="SF1"><bold>Figures S1</bold></xref>, <xref ref-type="supplementary-material" rid="SF2"><bold>S2</bold></xref>). Detection of unique genera is not surprising considering; 1) that these methods are more sensitive than traditional isolation and culturing methods and 2) attention is typically directed toward organisms that are pathogenic or beneficial to the crop in question, thus, neutral organisms and those found in low abundance are often not characterized and/or reported. We classified the fungal genera into guilds where possible, including pathotrophs (plant, fungal, and animal), saprotrophs, symbionts, epiphytes, and endophytes. While most of the detected organisms were saprotrophs, as would be expected to colonize organic soil matter, several plant pathogens were also identified at abundances of 0.5% or greater. The presence of pathogens in the soil, even at low levels, may signal a potential disease problem since these pathogens may serve as a reservoir for an outbreak that could express under disease-conducive conditions.</p>
<p>Generally, of the classifiable fungal genera, blueberry contained a higher pathogen load. For example, <italic>Coniochaeta</italic> spp. are known to be pathogens of woody perennials and thus is a potential pathogen of these crops (<xref ref-type="bibr" rid="B14">Damm et&#xa0;al., 2010</xref>). <italic>Athelia rolfsii</italic> causes Southern blight of solanaceous crops, beans and many other plants, but has not been reported as a pathogen of <italic>Vaccinium</italic> spp. (<xref ref-type="bibr" rid="B2">Aycock, 1961</xref>). <italic>Alternaria</italic> spp. and <italic>Cladosporium</italic> spp. are known to be pathogens of blueberry and cranberry (<xref ref-type="bibr" rid="B40">Polashock et&#xa0;al., 2017</xref>). Overall, N.J. blueberry contained a higher percent abundance of the identifiable potential pathogens <italic>Coniochaeta</italic>, <italic>Fusarium</italic>, <italic>Mollisia</italic>, and <italic>Phoma</italic> compared to cranberry, with several samples harboring statistical outliers of pathogen abundance far above average (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4</bold></xref>, represented as circles). In addition, <italic>Blumeria</italic>, a species of which causes powdery mildew on grasses, was uniquely identified in N.S. blueberries (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6</bold></xref>). N.J. cranberry however, had a comparatively low pathogen load, and low pathogen diversity (<xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3</bold></xref>, <xref ref-type="fig" rid="f4"><bold>4</bold></xref>, <xref ref-type="fig" rid="f6"><bold>6</bold></xref>).</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>The average percent abundance of the identified plant pathogenic fungal genera in the rhizosphere of cultivated highbush blueberry (blue) and cranberry (red) in New Jersey (N.J.). Data for lowbush blueberry (green) from Nova Scotia (<xref ref-type="bibr" rid="B50">Yurgel et&#xa0;al., 2017</xref>) shown for comparison.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1173023-g006.tif"/>
</fig>
<p>Several symbionts were also identified in the rhizosphere microbiome of N.J. blueberry and N.J. cranberry, including ericoid, ecto, and arbuscular mycorrhizae. These may represent important, but previously unknown, <italic>Vaccinium</italic> spp.- associated mycorrhizae. N.J. cranberry had a particularly high average abundance of <italic>Rhizoscyphus</italic> (note that <italic>Rhizoscyphus ericae</italic> = <italic>Hymenoscyphus ericae</italic>, but our data generally are not able to reach the species level and both genera are still listed in the databases used- hence both are used in this paper), compared to N.J. blueberry samples (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3</bold></xref>). Interestingly, it has been reported that mycorrhizae that proliferate in response to crop monoculture, may in fact be detrimental to that crop (<xref ref-type="bibr" rid="B31">Johnson et&#xa0;al., 1992</xref>). While endophytes and epiphytes were also detected, their importance in blueberry and cranberry production is not currently known. Similar to southern highbush blueberry, neither N.J. blueberry or N.J. cranberry carried a particularly high load of <italic>Oidiodendron</italic>, though cranberry carried significantly more than blueberry (<xref ref-type="bibr" rid="B33">Li et&#xa0;al., 2020</xref> and <xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5</bold></xref>).</p>
<p>We also compared known symbionts and pathogens with another small fruit crop grown in the northeast, grape. Similar to blueberry, grape exhibited very low levels of both <italic>Oidiodendron</italic> sp. and <italic>Hymenoscyphus</italic> sp., both of which are known to be ericoid mycorhizzae (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>), though no evidence of <italic>Russula</italic> sp. was found (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5C</bold></xref>). The analysis of pathogens showed that grape and blueberry show a similarly high abundance of <italic>Fusarium</italic> sp. (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>), which is shown to be a potential cause of decline in grape (<xref ref-type="bibr" rid="B27">Highet and Nair, 1995</xref>; <xref ref-type="bibr" rid="B7">Bustamante et&#xa0;al., 2022</xref>).</p>
<p>When compared to findings in southern highbush, as described in <xref ref-type="bibr" rid="B33">Li et&#xa0;al. (2020)</xref>, Hyaloscyphaceae, Leotiaceae, Pezoloma, and Hyaloscypha were not found in meaningful abundance in N.J. blueberry and N.J. cranberry.</p>
<p>Bacterial composition was similar between N.J. blueberry, N.J. cranberry, and N.S. blueberry. One of the top phyla however, chloroflexi, was found to differ significantly in abundance between blueberry and cranberry soils. Similarity in composition continued at the genus level, where little difference was found between N.J. blueberry and N.J. cranberry at abundances greater than 1.0%. Both N.J. blueberry and N.J. cranberry were also found to have similar bacterial alpha diversity indexes. As compared to southern highbush blueberry (<xref ref-type="bibr" rid="B33">Li et&#xa0;al., 2020</xref>), both N.J. blueberry and N.J. cranberry contained high levels of proteobacteria, acidobacteria, and actinobacteria making up the majority of bacterial phyla. This trend followed other studies that identified similarly high abundance of proteobacteria, acidobacteria, and actinobacteria in <italic>Vaccinium</italic> soils (<xref ref-type="bibr" rid="B33">Li et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B38">Morvan et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B51">Zhou et&#xa0;al., 2022</xref>). However, both N.J. blueberry and N.J. cranberry contained a heavier load of Firmicutes than southern highbush blueberry.</p>
<p>Only a small number of diseases of blueberry and cranberry are known to be caused by bacteria, with some of the agents being; <italic>Pseudomonas syringae</italic>, <italic>Xylella fastidiosa</italic>, <italic>Burkholderia andropogonis</italic>, <italic>Agrobacterium</italic> spp., <italic>Rhizobium</italic> spp. and <italic>Candidatus</italic> spp. (which cause phytoplasma diseases) (<xref ref-type="bibr" rid="B40">Polashock et&#xa0;al., 2017</xref>). At the genus level, several of these pathogens were identified at similar abundance in N.J. blueberry and N.J. cranberry rhizospheres including <italic>Pseudomonas, Burkholderia, Rhizobium, and Candidatus</italic>, although generally at low average abundances (&lt;0.5%) (<xref ref-type="supplementary-material" rid="SF3"><bold>Figures S3</bold></xref>, <xref ref-type="supplementary-material" rid="SF4"><bold>S4</bold></xref>).</p>
<p>Rhizophere microbiome diversity is thought to be an important indicator of soil health (<xref ref-type="bibr" rid="B6">Berendsen et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B48">Wallenstein, 2017</xref>; <xref ref-type="bibr" rid="B42">Schloter et&#xa0;al., 2018</xref>), and comparison in diversity can help to elucidate potential pathogenic influences, nutritional weaknesses, or microbial niches that have developed in differing soil samples. To help in visualizing these complex communities, we generated representations of microbial diversity for the largest and most diverse bacterial and fungal phyla (<xref ref-type="fig" rid="f7"><bold>Figures&#xa0;7</bold></xref>, <xref ref-type="fig" rid="f8"><bold>8</bold></xref>). The Proteobacteria exhibited a high level of diversity, as indicated by the relatively high number of large hubs, of various intermediate colors, and a high level of connectivity between nodes. There do remain, however, a few nodes that appear to be dominated in either the cranberry or blueberry groups, suggesting that niches have developed for particular taxa, or that the ecosystem does not support the other associated organisms. Conversely in Firmicutes, there are few noteworthy classes that are attached to the parent phylum. Though the phylum Firmicutes may not be as large as that of Proteobacteria, of the two largest classes represented here, Firmicutes appears to be more dominated by the cranberry group. This potentially indicates that there is a lack of diversity of Firmicutes in blueberry.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Diversity graphs of unique OTUs found in the bacteria, Proteobacteria <bold>(A)</bold> and Firmicutes <bold>(B)</bold>, from the rhizosphere soil of highbush blueberry and cranberry from New Jersey. Head nodes are labeled by phylum, lower nodes follow taxonomic hierarchy identified in available OTUs to the family level. Node area is proportional to the number of unique OTUs identified. Node color is a proportion of the number of unique highbush blueberry (blue) OTUs as compared to unique cranberry (red) OTUs.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1173023-g007.tif"/>
</fig>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Diversity graphs of unique OTUs found in fungi, Ascomycota <bold>(A)</bold> and Basidiomycota <bold>(B)</bold>, from the rhizosphere soil of highbush blueberry and cranberry from New Jersey. Head nodes are labeled by phylum, lower nodes follow taxonomic hierarchy identified in available OTUs to the family level. Node area is proportional to the number of unique OTUs identified. Node color is a proportion of the number of unique highbush blueberry (blue) OTUs as compared to unique cranberry (red) OTUs. Notable pathogen and symbiont families mentioned in the text are also labeled.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1173023-g008.tif"/>
</fig>
<p>Sequencing methods and bioinformatics tools are evolving rapidly and application of this technology allows us to gain insight into the rhizosphere microbiome, a key component of the phytobiome that impacts plant health. The OTUs generated from these samples were used to identify primarily fungi and bacteria that comprise the microbial communities that are associated with the rhizosphere of blueberries and cranberries in New Jersey. While we found similar communities between blueberry and cranberry in New Jersey, we noted differences between our data and that of the lowbush blueberry rhizosphere in Nova Scotia. Although lowbush blueberry (<italic>Vaccinium angustifolium</italic>) is in the same genus, the environmental conditions, management practices and soil characteristics are very different between the two locations (N.J. and Nova Scotia, Canada). Where similarity may still occur, is in characterizing the factors associated with poor plant performance as many of the diseases that are known to affect these crops overlap (<xref ref-type="bibr" rid="B40">Polashock et&#xa0;al., 2017</xref>). Using these methods, allows the detection and subsequent targeting of specific organisms using chemical control and/or modification of cultural practices. At the same time, enhancement of beneficial organisms, such as ericoid mycorrhizae, can also be monitored.</p>
<p>Many years ago, Waksman (<xref ref-type="bibr" rid="B47">Waksman, 1918</xref>) stated &#x201c;The cranberry soils are so distinctly different from ordinary soils that it was thought for a long time that no very large number of bacteria can exist in them and that the microbial population consists predominantly of molds&#x201d;. Though the uniqueness of the soils in which cranberries (and blueberries) holds true, we have shown that there are many phyla of bacteria and fungi that are adapted to these crops and the environments in which they grow. As we further characterize this myriad of organisms, we ultimately hope to characterize the microorganism(s) specifically associated with crop decline.</p>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The data presented in the study are deposited in the NCBI Sequence Read Archive under BioProject ID: PRJNA956268</p>
</sec>
<sec id="s6" sec-type="author-contributions">
<title>Author contributions</title>
<p>The project was conceived and directed by JP and PO. PO directed field collections and JP directed all laboratory work. JK and LE performed various data analyses and created the figures. JP and JK wrote the manuscript draft with contributions and editing by PO and LE. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s7" sec-type="funding-information">
<title>Funding</title>
<p>This work was, in part, supported by funding from the New Jersey Blueberry and Cranberry Research Council and USDA-NRCS-CIG grant 69-2B29-15-246 to P. V. Oudemans.</p>
</sec>
<sec id="s8" sec-type="COI-statement">
<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>
<p>The reviewer MD declared a shared affiliation with the authors JP, JK to the handling editor at the time of review.</p>
</sec>
<sec id="s9" sec-type="disclaimer">
<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 id="s10" sec-type="supplementary-material">
<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/fpls.2023.1173023/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2023.1173023/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Image_1.tiff" id="SF1" mimetype="image/tiff">
<label>Supplementary Figure&#xa0;1</label>
<caption>
<p>Diversity of fungal genera from the rhizosphere soil of highbush blueberry from New Jersey. Genera with an average relative abundance greater than 0.5% across samples are labeled. The height of the bars of the outer ring correspond to the average relative abundance. Each color represents a different phylum, and phyla with an average abundance greater than 0.5% are labeled. Tree nodes follow taxonomic hierarchy from kingdom (central) to genus (exterior).</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Image_2.tiff" id="SF2" mimetype="image/tiff">
<label>Supplementary Figure&#xa0;2</label>
<caption>
<p>Diversity of fungal genera from the rhizosphere soil of cranberry from New Jersey. Genera with an average relative abundance greater than 0.5% across samples are labeled. The heights of the bars of the outer ring correspond to the average relative abundance. Each color represents a different phylum and phyla with an average abundance greater than 0.5% are labeled. Tree nodes follow taxonomic hierarchy from kingdom (central) to genus (exterior).</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Image_3.tiff" id="SF3" mimetype="image/tiff">
<label>Supplementary Figure&#xa0;3</label>
<caption>
<p>Diversity of bacterial genera from the rhizosphere soil of highbush blueberry from New Jersey. Genera with an average relative abundance greater than 0.5% across samples are labeled. The height of the bars of the outer ring correspond to the average relative abundance. Each color represents a different phylum and phyla with an average abundance greater than 0.5% are labeled. Tree nodes follow taxonomic hierarchy from kingdom (central) to genus (exterior).</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Image_4.tiff" id="SF4" mimetype="image/tiff">
<label>Supplementary Figure&#xa0;4</label>
<caption>
<p>Diversity of bacterial genera from the rhizosphere soil of cranberry from New Jersey. Genera with an average relative abundance greater than 0.5% across samples are labeled. The heights of the bars of the outer ring correspond to the average relative abundance. Each color represents a different phylum and phyla with an average abundance greater than 0.5% are labeled. Tree nodes follow taxonomic hierarchy from kingdom (central) to genus (exterior).</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Table_1.docx" id="ST1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document">
<label>Supplementary Table&#xa0;1</label>
<caption>
<p>Locations, coordinates, and soil types for blueberry and cranberry rhizosphere soil collections.</p>
</caption>
</supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Asnicar</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Weingart</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Tickle</surname> <given-names>T. L.</given-names>
</name>
<name>
<surname>Huttenhower</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Segata</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Compact graphical representation of phylogenetic data and metadata with GraPhlAn</article-title>. <source>PeerJ</source> <volume>3</volume>, <elocation-id>e1029</elocation-id>. doi: <pub-id pub-id-type="doi">10.7717/peerj.1029</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aycock</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>1961</year>). <article-title>Summation-symposium on <italic>Sclerotium rolfsii</italic>
</article-title>. <source>Phytopathology</source> <volume>51</volume>, <fpage>107</fpage>&#x2013;<lpage>128</lpage>.</citation>
</ref>
<ref id="B3">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Ballinger</surname> <given-names>W. E.</given-names>
</name>
</person-group> (<year>1966</year>). &#x201c;<article-title>Soil management, nutrition, and fertilizer practices</article-title>,&#x201d; in <source>Blueberry culture</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Eck</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Childers</surname> <given-names>N. F.</given-names>
</name>
</person-group> (<publisher-loc>New Brunswick, NJ</publisher-loc>: <publisher-name>Rutgers University Press</publisher-name>), <fpage>132</fpage>.</citation>
</ref>
<ref id="B4">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Barr</surname> <given-names>D. J. S.</given-names>
</name>
</person-group> (<year>1990</year>). &#x201c;<article-title>Phylum chytridiomycota</article-title>,&#x201d; in <source>Handbook of Protoctista</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Margulis</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Corliss</surname> <given-names>J. O.</given-names>
</name>
<name>
<surname>Melkonian</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Chapman</surname> <given-names>D. J.</given-names>
</name>
</person-group> (<publisher-loc>Sudbury, MA, USA</publisher-loc>: <publisher-name>Jones and Bartlett</publisher-name>), <fpage>454</fpage>&#x2013;<lpage>466</lpage>.</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bennett</surname> <given-names>A. J.</given-names>
</name>
<name>
<surname>Bending</surname> <given-names>G. D.</given-names>
</name>
<name>
<surname>Chandler</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Hilton</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Mills</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Meeting the demand for crop production: the challenge of yield decline in crops grown in short rotations</article-title>. <source>Biol. Rev.</source> <volume>87</volume>, <fpage>52</fpage>&#x2013;<lpage>71</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1469-185X.2011.00184.x</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Berendsen</surname> <given-names>R. L.</given-names>
</name>
<name>
<surname>Pieterse</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>Bakker</surname> <given-names>P. A.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>The rhizosphere microbiome and plant health</article-title>. <source>Trends Plant Sci.</source> <volume>17</volume>, <fpage>478</fpage>&#x2013;<lpage>486</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.tplants.2012.04.001</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bustamante</surname> <given-names>M. I.</given-names>
</name>
<name>
<surname>Elfar</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Smith</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Bettiga</surname> <given-names>L. J.</given-names>
</name>
<name>
<surname>Tian</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Torres</surname> <given-names>G. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>First report of fusarium annulatum associated with young vine decline in California</article-title>. <source>Plant Dis.</source> <volume>106</volume> (<issue>10</issue>), <fpage>2752</fpage>. doi: <pub-id pub-id-type="doi">10.1094/PDIS-12-21-2790-PDN</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Capone</surname> <given-names>K. A.</given-names>
</name>
<name>
<surname>Dowd</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Stamatas</surname> <given-names>G. N.</given-names>
</name>
<name>
<surname>Nikolovski</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Diversity of the human skin microbiome early in life</article-title>. <source>J. Investig. Dermatol.</source> <volume>131</volume>, <fpage>2026</fpage>&#x2013;<lpage>2032</lpage>. doi: <pub-id pub-id-type="doi">10.1038/jid.2011.168</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Caporaso</surname> <given-names>J. G.</given-names>
</name>
<name>
<surname>Kuczynski</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Stombaugh</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Bittinger</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Bushman</surname> <given-names>F. D.</given-names>
</name>
<name>
<surname>Costello</surname> <given-names>E. K.</given-names>
</name>
<etal/>
</person-group>. (<year>2010</year>). <article-title>QIIME allows analysis of high-throughput community sequencing data</article-title>. <source>Nat. Methods</source> <volume>7</volume>, <fpage>335</fpage>&#x2013;<lpage>336</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nmeth.f.303</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Caporaso</surname> <given-names>J. G.</given-names>
</name>
<name>
<surname>Lauber</surname> <given-names>C. L.</given-names>
</name>
<name>
<surname>Walters</surname> <given-names>W. A.</given-names>
</name>
<name>
<surname>Berg-Lyons</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Lozupone</surname> <given-names>C. A.</given-names>
</name>
<name>
<surname>Turnbaugh</surname> <given-names>P. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>108</volume> (<supplement>Supplement 1</supplement>), <fpage>4516</fpage>&#x2013;<lpage>4522</lpage>. doi: <pub-id pub-id-type="doi">10.1073/pnas.1000080107</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chaparro</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Badri</surname> <given-names>D. V.</given-names>
</name>
<name>
<surname>Vivanco</surname> <given-names>J. M.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Rhizosphere microbiome assemblage is affected by plant development</article-title>. <source>ISME J.</source> <volume>8</volume>, <fpage>790</fpage>&#x2013;<lpage>803</lpage>. doi: <pub-id pub-id-type="doi">10.1038/ismej.2013.196</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chou</surname> <given-names>M. Y.</given-names>
</name>
<name>
<surname>Vanden Heuvel</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Bell</surname> <given-names>T. H.</given-names>
</name>
<name>
<surname>Panke-Buisse</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kao-Kniffin</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Vineyard under-vine floor management alters soil microbial composition, while the fruit microbiome shows no corresponding shifts</article-title>. <source>Sci. Rep.</source> <volume>8</volume> (<issue>1</issue>), <fpage>1</fpage>&#x2013;<lpage>9</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-018-29346-1</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Costa</surname> <given-names>R.</given-names>
</name>
<name>
<surname>G&#xf6;tz</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Mrotzek</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Lottmann</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Berg</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Smalla</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Effects of site and plant species on rhizosphere community structure as revealed by molecular analysis of microbial guilds</article-title>. <source>FEMS Microbiol. Ecol.</source> <volume>56</volume>, <fpage>236</fpage>&#x2013;<lpage>249</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1574-6941.2005.00026.x</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Damm</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Fourie</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Crous</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title><italic>Coniochaeta</italic> (<italic>Lecythophora</italic>), collophora gen. nov. and <italic>Phaeomoniella</italic> species associated with wood necroses of <italic>Prunus</italic> trees</article-title>. <source>Persoonia</source> <volume>24</volume>, <fpage>60</fpage>. doi: <pub-id pub-id-type="doi">10.3767/003158510X500705</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>DeSantis</surname> <given-names>T. Z.</given-names>
</name>
<name>
<surname>Hugenholtz</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Larsen</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Rojas</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Brodie</surname> <given-names>E. L.</given-names>
</name>
<name>
<surname>Keller</surname> <given-names>K.</given-names>
</name>
<etal/>
</person-group>. (<year>2006</year>). <article-title>Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB</article-title>. <source>Appl. Environ. Microbiol.</source> <volume>72</volume>, <fpage>5069</fpage>&#x2013;<lpage>5072</lpage>. doi: <pub-id pub-id-type="doi">10.1128/AEM.03006-05</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dowd</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Callaway</surname> <given-names>T. R.</given-names>
</name>
<name>
<surname>Wolcott</surname> <given-names>R. D.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>McKeehan</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Hagevoort</surname> <given-names>R. G.</given-names>
</name>
<etal/>
</person-group>. (<year>2008</year>b). <article-title>Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP)</article-title>. <source>BMC Microbiol.</source> <volume>8</volume>, <fpage>125</fpage>. doi: <pub-id pub-id-type="doi">10.1186/1471-2180-8-125</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dowd</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wolcott</surname> <given-names>R. D.</given-names>
</name>
<name>
<surname>Domingo</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Carroll</surname> <given-names>J. A.</given-names>
</name>
</person-group> (<year>2008</year>a). <article-title>Bacterial tag&#x2013;encoded FLX amplicon pyrosequencing (bTEFAP) for microbiome studies: bacterial diversity in the ileum of newly weaned salmonella-infected pigs</article-title>. <source>Foodborne Pathog. Dis.</source> <volume>5</volume>, <fpage>459</fpage>&#x2013;<lpage>472</lpage>. doi: <pub-id pub-id-type="doi">10.1089/fpd.2008.0107</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Eck</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>1988</year>). <source>Blueberry science</source> (<publisher-loc>New Brunswick</publisher-loc>: <publisher-name>Rutgers University Press</publisher-name>).</citation>
</ref>
<ref id="B19">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Eck</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>1990</year>). <source>The American cranberry</source> (<publisher-loc>New Brunswick, NJ</publisher-loc>: <publisher-name>Rutgers University Press</publisher-name>).</citation>
</ref>
<ref id="B20">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Eck</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Childers</surname> <given-names>N. F.</given-names>
</name>
</person-group> (<year>1966</year>). <source>Blueberry culture</source> (<publisher-loc>New Brunswick</publisher-loc>: <publisher-name>Rutgers University Press</publisher-name>).</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Edgar</surname> <given-names>R. C.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Search and clustering orders of magnitude faster than BLAST</article-title>. <source>Bioinformatics</source> <volume>26</volume>, <fpage>2460</fpage>&#x2013;<lpage>2461</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btq461</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Eren</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Zozaya</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Taylor</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>Dowd</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Martin</surname> <given-names>D. H.</given-names>
</name>
<name>
<surname>Ferris</surname> <given-names>M. J.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Exploring the diversity of gardnerella vaginalis in the genitourinary tract microbiota of monogamous couples through subtle nucleotide variation</article-title>. <source>PloS One</source> <volume>6</volume>, <elocation-id>e26732</elocation-id>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0026732</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fierer</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Embracing the unknown: disentangling the complexities of the soil microbiome</article-title>. <source>Nat. Rev. Microbiol.</source> <volume>15</volume>, <fpage>579</fpage>&#x2013;<lpage>590</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nrmicro.2017.87</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garside</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Bell</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Magarey</surname> <given-names>R. C.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Monoculture yield decline-fact not fiction</article-title>. <source>Proc. Int. Soc Sugar Cane Technol.</source> <volume>24</volume>, <fpage>16</fpage>&#x2013;<lpage>21</lpage>.</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gleason</surname> <given-names>F. H.</given-names>
</name>
<name>
<surname>Kagami</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Lefevre</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Sime-Ngando</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>The ecology of chytrids in aquatic ecosystems: roles in food web dynamics</article-title>. <source>Fungal Biol. Rev.</source> <volume>22</volume>, <fpage>17</fpage>&#x2013;<lpage>25</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.fbr.2008.02.001</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hagberg</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Schult</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Swart</surname> <given-names>P. J.</given-names>
</name>
</person-group> (<year>2008</year>). &#x201c;<article-title>Exploring network structure, dynamics, and function using NetworkX</article-title>,&#x201d; in <source>Proceedings of the 7th Python in science conference (SciPy2008)</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Varoquaux</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Vaught</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Millman</surname> <given-names>J.</given-names>
</name>
</person-group> (<publisher-name>Pasadena, CA, USA</publisher-name>), <fpage>11</fpage>&#x2013;<lpage>15</lpage>.</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Highet</surname> <given-names>A. S.</given-names>
</name>
<name>
<surname>Nair</surname> <given-names>N. G.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title><italic>Fusarium oxysporum</italic> associated with grapevine decline in the hunter valley, NSW, Australia</article-title>. <source>Aust. J. Grape Wine Res.</source> <volume>1</volume> (<issue>1</issue>), <fpage>48</fpage>&#x2013;<lpage>50</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1755-0238.1995.tb00077.x</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hunter</surname> <given-names>J. D.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Matplotlib: a 2D graphics environment</article-title>. <source>Comput. Sci. Eng.</source> <volume>9</volume>, <fpage>90</fpage>&#x2013;<lpage>95</lpage>. doi: <pub-id pub-id-type="doi">10.1109/MCSE.2007.55</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jagdale</surname> <given-names>G. B.</given-names>
</name>
<name>
<surname>Holladay</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Brannen</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Cline</surname> <given-names>W. O.</given-names>
</name>
<name>
<surname>Agudelo</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Nyczepir</surname> <given-names>A. P.</given-names>
</name>
<etal/>
</person-group>. (<year>2013</year>). <article-title>Incidence and pathogenicity of plant-parasitic nematodes associated with blueberry (<italic>Vaccinium</italic> spp.) replant disease in Georgia and north Carolina</article-title>. <source>J. Nematol.</source> <volume>45</volume>, <fpage>92</fpage>&#x2013;<lpage>98</lpage>.</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jeffers</surname> <given-names>S. N.</given-names>
</name>
</person-group> (<year>1988</year>). <article-title><italic>Phytophthora</italic> species associated with a cranberry decline syndrome in Wisconsin</article-title>. <source>Phytopathology</source> <volume>78</volume>, <fpage>1572</fpage>.</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Johnson</surname> <given-names>N. C.</given-names>
</name>
<name>
<surname>Copeland</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Crookston</surname> <given-names>R. K.</given-names>
</name>
<name>
<surname>Pfleger</surname> <given-names>F. L.</given-names>
</name>
</person-group> (<year>1992</year>). <article-title>Mycorrhizae: possible explanation for yield decline with continuous corn and soybean</article-title>. <source>Agron. J.</source> <volume>84</volume>, <fpage>387</fpage>&#x2013;<lpage>390</lpage>. doi: <pub-id pub-id-type="doi">10.2134/agronj1992.00021962008400030007x</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ketcheson</surname> <given-names>J. W.</given-names>
</name>
</person-group> (<year>1980</year>). <article-title>Long-range effects of intensive cultivation and monoculture on the quality of southern Ontario soils</article-title>. <source>Can. J. Soil Sci.</source> <volume>60</volume>, <fpage>403</fpage>&#x2013;<lpage>410</lpage>. doi: <pub-id pub-id-type="doi">10.4141/cjss80-045</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Mavrodi</surname> <given-names>O. V.</given-names>
</name>
<name>
<surname>Hou</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Blackmon</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Babiker</surname> <given-names>E. M.</given-names>
</name>
<name>
<surname>Mavrodi</surname> <given-names>D. V.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Comparative analysis of rhizosphere microbiomes of southern highbush blueberry (<italic>Vaccinium corymbosum</italic> l.), darrow&#x2019;s blueberry (<italic>V. darrowii</italic> camp), and rabbiteye blueberry (<italic>V. virgatum</italic> aiton)</article-title>. <source>Front. Microbiol.</source> <volume>11</volume>, <elocation-id>370</elocation-id>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2020.00370</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Mainland</surname> <given-names>C. M.</given-names>
</name>
</person-group> (<year>1966</year>). &#x201c;<article-title>Propagation and planting</article-title>,&#x201d; in <source>Blueberry culture</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Eck</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Childers</surname> <given-names>N. F.</given-names>
</name>
</person-group> (<publisher-loc>New Brunswick</publisher-loc>: <publisher-name>Rutgers University Press</publisher-name>), <fpage>111</fpage>&#x2013;<lpage>131</lpage>.</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mandal</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Van Treuren</surname> <given-names>W.</given-names>
</name>
<name>
<surname>White</surname> <given-names>R. A.</given-names>
</name>
<name>
<surname>Eggesb&#xf8;</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Knight</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Peddada</surname> <given-names>S. D.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Analysis of composition of microbiomes: a novel method for studying microbial composition</article-title>. <source>Microb. Ecol. Health Dis.</source> <volume>26</volume>, <fpage>27663</fpage>. doi: <pub-id pub-id-type="doi">10.3402/mehd.v26.27663</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Marschner</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2012</year>). &#x201c;<article-title>Rhizosphere biology</article-title>,&#x201d; in <source>Marschner's mineral nutrition of higher plants</source>, <edition>3rd ed.</edition> Ed. <person-group person-group-type="editor">
<name>
<surname>Marschner</surname> <given-names>P.</given-names>
</name>
</person-group> (<publisher-loc>Oxford, UK</publisher-loc>: <publisher-name>Academic Press</publisher-name>), <fpage>369</fpage>&#x2013;<lpage>388</lpage>.</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mazzola</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Manici</surname> <given-names>L. M.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Apple replant disease: role of microbial ecology in cause and control</article-title>. <source>Annu. Rev. Phytopathol.</source> <volume>50</volume>, <fpage>45</fpage>&#x2013;<lpage>65</lpage>. doi: <pub-id pub-id-type="doi">10.1146/annurev-phyto-081211-173005</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Morvan</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Meglouli</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Sahraoui</surname> <given-names>Loun&#xe8;s-HadjA.</given-names>
</name>
<name>
<surname>Hijri</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Into the wild blueberry (<italic>Vaccinium angustifolium</italic>) rhizosphere microbiota</article-title>. <source>Environ. Microbiol.</source> <volume>22</volume> (<issue>9</issue>), <fpage>3803</fpage>&#x2013;<lpage>3822</lpage>. doi: <pub-id pub-id-type="doi">10.1111/1462-2920.15151</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nguyen</surname> <given-names>N. H.</given-names>
</name>
<name>
<surname>Song</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Bates</surname> <given-names>S. T.</given-names>
</name>
<name>
<surname>Branco</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Tedersoo</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Menke</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild</article-title>. <source>Fungal Ecol.</source> <volume>20</volume>, <fpage>241</fpage>&#x2013;<lpage>248</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.funeco.2015.06.006</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Polashock</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>Caruso</surname> <given-names>F. L.</given-names>
</name>
<name>
<surname>Averill</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Schilder</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2017</year>). <source>Compendium of blueberry, cranberry, and lingonberry diseases and pests</source> (<publisher-loc>St. Paul, MN</publisher-loc>: <publisher-name>APS Press, American Phytopathological Society</publisher-name>).</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Porter</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Crookston</surname> <given-names>R. K.</given-names>
</name>
<name>
<surname>Ford</surname> <given-names>J. H.</given-names>
</name>
<name>
<surname>Huggins</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Lueschen</surname> <given-names>W. E.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>Interrupting yield depression in monoculture corn: comparative effectiveness of grasses and dicots</article-title>. <source>Agron. J.</source> <volume>89</volume>, <fpage>247</fpage>&#x2013;<lpage>250</lpage>. doi: <pub-id pub-id-type="doi">10.2134/agronj1997.00021962008900020015x</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schloter</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Nannipieri</surname> <given-names>P.</given-names>
</name>
<name>
<surname>S&#xf8;rensen</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>van Elsas</surname> <given-names>J. D.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Microbial indicators for soil quality</article-title>. <source>Biol. Fertil. Soils</source> <volume>54</volume> (<issue>1</issue>), <fpage>1</fpage>&#x2013;<lpage>10</lpage>. doi: <pub-id pub-id-type="doi">10.1007/s00374-017-1248-3</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Stretch</surname> <given-names>A. W.</given-names>
</name>
<name>
<surname>Carris</surname> <given-names>L. M.</given-names>
</name>
</person-group> (<year>2017</year>). &#x201c;<article-title>Red gall</article-title>,&#x201d; in <source>Compendium of blueberry, cranberry, and lingonberry diseases and pests</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Polashock</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>Caruso</surname> <given-names>F. L.</given-names>
</name>
<name>
<surname>Averill</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Schilder</surname> <given-names>A.</given-names>
</name>
</person-group> (<publisher-loc>Saint Paul, MN</publisher-loc>: <publisher-name>American Phytopathological Society Press</publisher-name>), <fpage>108</fpage>&#x2013;<lpage>109</lpage>.</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Swanson</surname> <given-names>K. S.</given-names>
</name>
<name>
<surname>Dowd</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Suchodolski</surname> <given-names>J. S.</given-names>
</name>
<name>
<surname>Middelbos</surname> <given-names>I. S.</given-names>
</name>
<name>
<surname>Vester</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Barry</surname> <given-names>K. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Phylogenetic and gene-centric metagenomics of the canine intestinal microbiome reveals similarities with humans and mice</article-title>. <source>ISME J.</source> <volume>5</volume>, <fpage>639</fpage>&#x2013;<lpage>649</lpage>. doi: <pub-id pub-id-type="doi">10.1038/ismej.2010.162</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tewoldemedhin</surname> <given-names>Y. T.</given-names>
</name>
<name>
<surname>Mazzola</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Labuschagne</surname> <given-names>I.</given-names>
</name>
<name>
<surname>McLeod</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>A multi-phasic approach reveals that apple replant disease is caused by multiple biological agents, with some agents acting synergistically</article-title>. <source>Soil Biol. Biochem.</source> <volume>43</volume>, <fpage>1917</fpage>&#x2013;<lpage>1927</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.soilbio.2011.05.014</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Van der Walt</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Colbert</surname> <given-names>S. C.</given-names>
</name>
<name>
<surname>Varoquaux</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>The NumPy array: a structure for efficient numerical computation</article-title>. <source>Comput. Sci. En.</source> <volume>13</volume>, <fpage>22</fpage>&#x2013;<lpage>30</lpage>. doi: <pub-id pub-id-type="doi">10.1109/MCSE.2011.37</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Waksman</surname> <given-names>S. A.</given-names>
</name>
</person-group> (<year>1918</year>). <article-title>The occurrence of azotobacter in cranberry soils</article-title>. <source>Science</source> <volume>48</volume>, <fpage>653</fpage>&#x2013;<lpage>654</lpage>. doi: <pub-id pub-id-type="doi">10.1126/science.48.1252.653</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wallenstein</surname> <given-names>M. D.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Managing and manipulating the rhizosphere microbiome for plant health: a systems approach</article-title>. <source>Rhizosphere</source> <volume>3</volume>, <fpage>230</fpage>&#x2013;<lpage>232</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.rhisph.2017.04.004</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>White</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Burns</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Taylor</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>1990</year>). &#x201c;<article-title>Amplification and direct sequencing of fungal ribsomal RNA genes for phylogenetics</article-title>,&#x201d; in <source>PCR protocols</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Innis</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Gelfand</surname> <given-names>D. H.</given-names>
</name>
<name>
<surname>Sninsky</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>White</surname> <given-names>T. J.</given-names>
</name>
</person-group> (<publisher-loc>New York</publisher-loc>: <publisher-name>Academic Press Inc.</publisher-name>), <fpage>315</fpage>&#x2013;<lpage>322</lpage>.</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yurgel</surname> <given-names>S. N.</given-names>
</name>
<name>
<surname>Douglas</surname> <given-names>G. M.</given-names>
</name>
<name>
<surname>Comeau</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Mammoliti</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Dusault</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Percival</surname> <given-names>D.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Variation in bacterial and eukaryotic communities associated with natural and managed wild blueberry habitats</article-title>. <source>Phytobiomes</source> <volume>1</volume>, <fpage>102</fpage>&#x2013;<lpage>113</lpage>. doi: <pub-id pub-id-type="doi">10.1094/PBIOMES-03-17-0012-R</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Shao</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Long</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Effects of soil properties and microbiome on highbush blueberry (<italic>Vaccinium corymbosum</italic>) growth</article-title>. <source>Agronomy</source> <volume>12</volume> (<issue>6</issue>), <fpage>1263</fpage>. doi: <pub-id pub-id-type="doi">10.3390/agronomy12061263</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>
