Abstract An unusual decline and collapse of young established trees known as “rapid apple decline” (RAD) has become a major concern for apple growers, particularly in the northeastern United States. This decline is characterized by stunted growth, pale yellow to reddish leaves, and tree collapse within weeks after onset of symptoms. We studied declining apple trees to identify potential involvement of abiotic and biotic stresses. We used 16S and ITS to profile bacterial and fungal communities in the soil, rhizosphere, roots, and shoots and tested for the presence of six viruses in scions and rootstocks of symptomatic and asymptomatic trees. The viruses detected were not associated with RAD symptoms. Bacterial and fungal populations were highly variable in plant tissue, soil and rhizosphere samples, with bacteroidetes, firmicutes, proteobacteria, acidobacteria, and actinobacteria the predominant bacterial classes in various samples. ‘Alphaproteobacteria-rickettsiales’, a bacterial class usually reduced in water-limiting soils, had significantly low abundance in root samples of symptomatic trees. Basidiomycota and Ascomycota fungal classes were the most common fungal classes observed, but neither showed differential enrichment between symptomatic and asymptomatic trees. Analyzing weather data showed an extremely cold winter followed by drought in 2015–2016, which likely weakened the trees to make them more susceptible to varied stresses. In addition, similar physical and nutritional soil composition from symptomatic and asymptomatic trees rules out the role of nutritional stress in RAD. Necrotic lesions and wood decay symptoms dispersing from bark or vascular cambium towards the heartwood were observed primarily below the graft union of declining apple trees, suggesting that the rootstock is the originating point of RAD. We speculate that differences in abiotic factors such as moisture levels in declining roots in combination with extreme weather profiles might cause RAD but cannot clearly rule out the involvement of other factors.

Citation: Singh J, Silva KJP, Fuchs M, Khan A (2019) Potential role of weather, soil and plant microbial communities in rapid decline of apple trees. PLoS ONE 14(3): e0213293. https://doi.org/10.1371/journal.pone.0213293 Editor: Ricardo Aroca, Estacion Experimental del Zaidin, SPAIN Received: October 14, 2018; Accepted: February 18, 2019; Published: March 6, 2019 Copyright: © 2019 Singh et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: The read sequences from 16S bacterial data and ITS fungal data have been deposited in the National Center for Biotechnology Information (NCBI) sequence read archive (SRA) database under Bioproject PRJNA495655 and PRJNA495767, respectively. Funding: This research was supported by the National Institute of Food and Agriculture, U.S. Department of Agriculture, Hatch project NYC-625410 under accession number 1011771. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.

Introduction “Rapid apple decline” (RAD) describes a decline and collapse of young apple trees. This decline has become a major concern for growers in apple producing areas in the central, northeastern, and northwestern United States and in Ontario, Canada in recent years [1–4]. A similar situation was first described in north-central Washington State orchards in 1983 [5], and fifteen years later in southern British Columbia, Canada [6]. RAD is usually characterized by stunted tree growth, chlorotic canopy, and tree collapse within weeks after the onset of symptom development [2–4]. Symptoms generally first appear on one limb, with small, rolled leaves and reduced terminal growth, followed by the full tree canopy manifesting pale yellow to reddish leaves. Cankers and shedding are visible at the graft union, and wood necrosis progresses upstream to the trunk of the tree. The root system generally appears healthy, except that small feeder roots are absent [2, 5]. No common root rot pathogens or nutrient deficiency has been associated with RAD. After originating in one point of an orchard, the disorder seems to spread to adjacent trees. Symptomatic trees are usually removed from orchards due to poor productivity; however, if kept in the orchard, the symptoms may spread throughout the tree within a single year and ultimately lead to tree death. Despite extensive speculation on the potential causes of RAD, the causative agent is still unknown. For example, cultivar-rootstock incompatibility, extreme weather conditions, wood-boring insects, and pathogen infection have all been proposed as possible causes. Trees presenting incompatible grafts are frequently described as displaying breaks or malformations at the graft union, leaf chlorosis, early defoliation, plant wilt, and premature death. Unfavorable weather conditions have also been suggested to be involved in RAD [7–9]. For instance, freezing temperatures can cause direct injury to plant tissue, making them vulnerable to secondary abiotic or biotic stresses [8–9]. Likewise, drought or flooding may potentially cause retarded shoot and leaf growth, leaf chlorosis and defoliation, root necrosis, wilting, and eventually plant senescence [10–16]. Wood-boring insects can also cause serious damage or death of apple trees and have been identified in trees with RAD [17]. Insect infested trees usually have a sickly appearance, a sparse and pale-colored foliage, and can die with a heavy fruit crop during the fruit maturation stage. Although insect borer injury to the graft union or trunk has been reported in RAD-symptomatic areas, it is believed that insects take advantage of already declining trees, and compound the injury by providing an entryway for destructive fungi [2, 4, 17]. The involvement of plant pathogens in RAD is still a matter of speculation. Many important diseases of apple trees are caused by pathogens that initiate infections at wounds caused by insects, humans, machinery, fire, lightning, wind, hail, animals, or nutritional and physiological disorders [18–19]. Symptoms from other microbes such as wood-rotting and Phytophthora crown rot pathogens and mycoplasma-like organisms showed partial matches to RAD symptoms [5, 20–22]. In contrast, the necrotic symptoms of the inner bark on RAD-affected trees were more likely to be from latent apple viruses such as apple stem pitting virus and apple stem grooving virus [2, 4]. Although apple luteovirus 1 (ALV1) was recently characterized from trees with RAD [23], the association between the presence of ALV1 and RAD symptoms is weak. Therefore, no virus has been confirmed as the causative agent for RAD [2]. It might be possible that RAD is caused by an as yet unknown pathogen. The identification of microorganisms responsible for plant disease has relied mainly on culture-dependent techniques and PCR amplification of genomic DNA [24]. However, these techniques are specific and can miss certain infectious microorganisms or groups of organisms that cause diseases. Microbial communities have been shown to have synergistic effects by improving agronomical features, such as limiting or preventing attacks by phytopathogens [25–27] or causing the establishment and development of plant diseases [20, 24, 27–31]. High throughput sequencing approaches are powerful tools to determine the involvement of known and new organisms in disease etiology and permit the investigation of complex microbiomes [20, 27, 32, 33]. Routine 16S and ITS sequencing allow rapid identification of soil, rhizosphere, or endophyte bacterial and fungal communities, respectively [34–39]. For example, a bacterial community potentially interacting with a colonizing fungus was isolated from the grapevine trunk fungal disease “esca” [33]. Similarly, three fungi typically not considered root pathogens have been reported as causal agents of apple replant disease [20]. A comparative analysis between healthy and declining trees can help identify differentially abundant classes of bacteria and fungus and therefore assist in the identification of their putative role in RAD. We carried out a study in an orchard with a block of ‘HoneyCrisp’ trees showing symptoms of RAD to identify the potential role of different abiotic and biotic stresses. We characterized the morphological status of wood from RAD symptomatic trees to localize the origin of necrosis, compared microbial community profiles of soil, rhizosphere, roots, and shoots of symptomatic and asymptomatic trees to identify biotic agents involved in RAD, and analyzed weather data and physical characteristics of soil to investigate the involvement of abiotic stresses.

Materials and methods Site description and weather data The experimental orchard is located in a commercial apple orchard in Newark, Wayne County, New York in the United States (43°7’30”N, 77°5’31”W, and 170m elevation). The orchard has more than 30 cultivars of apples grafted on different rootstocks. The majority of declining trees were observed in a block of ‘HoneyCrisp’. The block had 1,700 ‘HoneyCrisp’ trees grafted onto the Malling 9 (M.9 NIC 29) rootstock. This block was established in 2010, in a high-density planting system, with 1.25m x 4.25m spacing between trees and rows. The soil type is characterized as silty loam and belongs to the hydrologic soil group B/D with moderate to very slow infiltration rate and moderate potential for frost action according to the United States Department of Agriculture (USDA)—Natural Resources Conservation Service Soil Climate Analysis Network (NRCS) (USDA-NRCS, 2017). The water movement in the most restrictive layer and the shrink-swell potential is low in this area. The area is considered moderately well drained, with low probability of flooding and ponding. Annual precipitation in the area varies from 787 to 1,447mm, and the frost-free period varies from 100 to 190 days. Weather datasets from 2013 to 2017 were accessed from three weather stations located in Phelps, Farmington, and Sodus, New York using the Network for Environment and Weather Applications (NEWA) (http://newa.cornell.edu) [40]. These three locations were within 15 miles of the ‘HoneyCrisp’ apple orchard selected for this study. Weather data included temperature (maximum, average, and minimum), precipitation, relative humidity, leaf wetness, and wind speed. There were no major differences in weather parameters among the three stations. Therefore, datasets from three weather stations were used to calculate monthly mean values for each year. Sampling for soil, microbiome and virus analysis Asymptomatic and symptomatic (declining) apple trees were dispersed across the experimental orchard. Four rows (R1-R4) in the north end of the block had clearly declining as well as healthy-looking trees. The two central rows (R2 and R3) were selected for soil and tree tissue sampling, avoiding outer open rows to minimize experimental error (S1A Fig). Three sets of trees (two asymptomatic trees on either side of a symptomatic tree) were sampled from R2 (S1B Fig). In R3, 20 (10 asymptomatic and 10 symptomatic) trees were randomly selected and sampled (S1C Fig). A total of 87 samples were collected from shoots (13 symptomatic, 16 asymptomatic), roots (13 symptomatic, 16 asymptomatic), and the rhizosphere (13 symptomatic, 16 asymptomatic) of ‘HoneyCrisp’ trees selected in R2 and R3. For sampling, shoots (30cm in length) were cut with pruning shears from three different positions in the tree canopy (top, middle, and base) and pooled together (Fig 1A). Soil surrounding the tree rootstock was dug out to expose lateral roots for root and rhizosphere sampling. Five root pieces (5-20cm) adjacent to the rootstock trunk were collected using pruning shears and a chisel (Fig 1B). Pruning shears, shovels, chisels, and augers used during sampling were cleaned with 50% bleach between uses to avoid cross contamination, and gloves were changed between samples. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 1. Schematic overview of the sampling strategy of soil and tree samples in a ‘HoneyCrisp’ orchard block with rapid apple decline (RAD). Sampling and pooling of shoots from three different positions in the tree canopy (top, middle, and base) (A), Sampling of root tissue and rhizosphere (B). Soil sampling for the analysis of soil physical and chemical properties (C). https://doi.org/10.1371/journal.pone.0213293.g001 A total of nine bulk soil samples were collected from six asymptomatic and three symptomatic trees in R2 using a soil auger (4.5 cm Ø). For each sample, six soil cores (~20 cm depth) were collected around the tree, approximately 20cm from the center of the tree trunk. The first 2 cm of each sample was discarded, and the remaining soil was combined, sieved with a 2-mm mesh, and homogenized (Fig 1C). One portion of all soil samples was kept in -80°C for DNA extraction. Another portion of each sample was used for soil physicochemical analysis at the Soil Lab, Cornell University in Ithaca, New York. Plant material and soil samples were kept on ice and carefully transported to the laboratory in a cold container. In the lab, roots were gently shaken inside sampling bags and the rhizosphere soil was carefully detached and transferred to a sterile container. Shoot, root, and rhizosphere samples were processed within 24 hours. Visual characterization of rootstocks of symptomatic tress At the end of the growing season in 2017, eighteen symptomatic apple trees, a subset of the trees that were sampled for microbiome and virus testing, were pulled out from rows R2 and R3. The trunks were cut 30cm above the graft union and brought to the lab for visual characterization of symptoms. Three parallel cuts (5 cm each) were made above and below the graft union to verify the presence of abnormal wood, as well as signs and/or symptoms of potential pathogen and insect damage. High-resolution images of the cross section of each cut piece were taken to analyze necrotic lesions as the percentage of necrotic wood area to the total cross section area. A macro in ImageJ (version 1.51 h - http://imagej.nih.gov/ij) was used to convert pictures to a binary format, assigning white and black colors to healthy and pigmented tissue, respectively. The outputs were used to estimate healthy and pigmented necrotic wood area. The image data were further evaluated to localize the origin and progress of the necrosis using the following scale: DS = 1, external wood discoloration, DS = 2, internal wood discoloration, DS = 3, co-occurrence of external and internal wood discoloration, DS = 4, edge-shaped discoloration (<50%), DS = 5, edge-shaped discoloration (>50%), and DS = 6, circular discoloration. Virus detection Shoot and root samples of asymptomatic and symptomatic trees were tested for the presence of viruses using double-antibody sandwich enzyme-linked immunosorbent assay (DAS-ELISA) and commercial antibodies (Bioreba, Reinach, Switzerland) according to the manufacturer's instructions. Viruses assayed were apple stem grooving virus (ASGV), apple stem pitting virus (ASPV), apple chlorotic leafspot virus (ACLSV), tomato ringspot virus (ToRSV), apple mosaic virus (APMV), and tobacco ringspot virus (TRSV). Briefly, shoot and root samples (0.5g) were homogenized in Bioreba universal extraction bags containing extraction buffer (5ml) using a HOMEX 6 (Bioreba, Reinach, Switzerland). DAS-ELISA was conducted as follow: Specific polyclonal antibodies were diluted 1,000 times in coating buffer, dispensed into wells of a microtiter plate (100 μl/well), and incubated overnight at 4°C. Plates were washed three times with phosphate-buffered saline (PBS) containing 0.05% Tween 20 and 500 μl of sample extracts were added to two wells and kept overnight at 4°C. After three washes with PBS-Tween, 100 μl of conjugate antibody diluted 1,000 times were added and incubated 4h at 30°C. After adding the substrate, absorbance readings were recorded using a microplate reader (EL800, Biotech. Instruments, USA) at 405 nm. Samples were considered virus-positive when the average absorbance values (A 405nm ) were at least three times higher than those of the negative control. Sample processing and DNA extraction Shoot and root samples were processed prior to DNA extraction. Shoots and roots were first washed with sterile water, blotted dry on sterile absorbent paper, and transferred to a laminar flow hood. The bark of the shoots and the cortex of the roots were removed to trim the sample ends. Wood samples were cut horizontally, and small wood chips were quickly collected using a sterile scalpel and transferred to sterile bags. Samples were immediately frozen in liquid nitrogen. Approximately 100 mg of shoot or root tissue were disrupted and homogenized using a tissue lyser II (Qiagen, Crawley, UK). Three 5 mm stainless steel beads were used per sample. Three rounds of 45 s of grinding were performed at a frequency of 30 Hz. Total genomic DNA was extracted from shoots and roots using DNeasy Plant Mini Kit (Qiagen, Germany), and microbial DNA was extracted from bulk soil and rhizosphere samples using the MoBio Power Soil DNA isolation kit (Qiagen, Germany), according to the manufacturer’s instructions. Quality and quantity of extracted DNA was assessed by NanoDrop (absorbance ratio at both 260/280 and 230/260 nm) and by electrophoresis on 1.5% agarose gels. Concentration of extracted DNA was adjusted to 5ng/μl and samples were stored at -80°C until use. Characterization of microbial communities Bacterial and fungal communities were characterized in soil, shoot, root, and rhizosphere samples. The gene‐specific primers V3-357F (5’-CCTACGGGNGGCWGCAG-3’) and V4-805R (5’-GACTACHVGGGTATCTAATCC-3’) targeting the V3 and V4 regions of 16S rRNA gene were used to study bacterial communities (Klindworth et al., 2013; Yim et al., 2015). The gene‐specific primers ITS1F (5’-TCCGTAGGTGAACCTGCGG-3’) and ITS4R (5’-TCCTCCGCTTATTGATATGC-3’), targeting the two transcribed intergenic spacers (ITS) ITS1 and ITS4 rDNA regions were used to study fungal communities (Manter and Vivanco, 2007). The forward and reverse Illumina overhang adapter sequences added to locus‐specific sequences were: 5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-3’ and 5’-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG-3’, respectively. DNA was amplified by PCR in a reaction mixture (25 μl final volume) consisting of 2.5 μl of microbial DNA (5 ng/μl), 5 μl of each amplicon PCR primer (1 μM), and 12.5 μl of 2x KAPA HiFi HotStart ReadyMix (KAPA Biosystems). The PCR amplification conditions for 16S rRNA region were as follows: 95°C for 3 min; 25 cycles of denaturation at 95°C for 30 s; hybridization at 55°C for 30 s; extension at 72°C for 30 s; and a final extension at 72°C for 5 min. For amplification of ITS rDNA region, PCR conditions were: 95°C for 5 min; 25 cycles of denaturation at 95°C for 30 s; hybridization at 55°C for 1 min; extension at 72°C for 1 min; and a final extension at 72°C for 10 min. The amplicons were sequenced using paired-end sequencing on an Illumina MiSeq instrument at the Institute of Biotechnology at Cornell University in Ithaca New York, United States. Raw sequence reads were de-multiplexed, low quality read ends were trimmed using Trimmomatic [41], and low-quality sequences were removed. QIIME2 (https://qiime2.org/) was used to perform the downstream diversity and taxonomy composition analysis. Corresponding paired end reads were merged and un-joined reads were discarded. Another quality processing was performed using ‘Deblur’ plugin in QIIME2 to remove chimeric sequences. Read quality profiles were visualized to retain high-quality sequences for 16S and ITS datasets. The remaining sequences were used to determine differences in bacterial and fungal communities between asymptomatic and symptomatic samples from root, shoot, soil, and rhizosphere, and to calculate the Shannon diversity index to obtain alpha and beta diversity statistics. Sequences were grouped to obtain operational taxonomic units (OTUs) with 97% similarity. The resulting OTUs were compared against the trained full-length Greengenes 13_8 OTUs database (http://greengenes.secondgenome.com/) for bacterial taxonomic classification. For ITS taxonomy analysis, a database was trained using the fungus sequences in UNITE (Fungal ITS) (https://unite.ut.ee/). Fungus OTUs were compared against the resulting database for taxonomic analysis of ITS sequences. A differential abundance analysis of the OTU was independently performed for bacterial and fungal datasets using a Kruskal-Wallis test with multiple correlation testing (FDR). The test was performed by separately comparing the asymptomatic and symptomatic samples within root, shoot, soil, and rhizosphere. The sample-specific differentially abundant OTUs were determined using a threshold p<0.01. The output files were visualized in QIIME2. Statistical analysis The means and standard deviations of various soil elements from symptomatic and asymptomatic trees were analyzed using student’s T-test assuming unequal variances. A p-value of 0.05 was used as a threshold to declare significant differences. Similarly, the monthly means from different weather variables were compared against five-year averages using a student’s T-test, assuming unequal variances with significant p-value threshold of less than 0.05. The differential abundance of bacterial and fungal classes were analyzed using the Kruskal-Wallis test with multiple correlation testing (FDR). Class enrichment data was retrieved from QIIME2 output and a principal component analysis (PCA) was performed in R (https://www.r-project.org/) to visualize sample-specific variation. The ggplot2 package in R was used to plot PAC biplots.

Discussion Tree decline/death is often caused by diseases such as rootstock fire blight, phytophthora root and crown rots, or apple replant disease, or abiotic factors like soil wetness, extreme cold, spring and fall freezes, or drought, as well as injuries to trunk, graft union, crown and roots by rodents or insect borers. However, rapid or sudden decline of established apple trees (RAD), a recent concern to the apple industry, is diagnosed when symptoms do not match any of the above. We addressed this concern by conducting a comprehensive study to identify the role of soil nutrition, weather conditions, fungal and bacterial communities, as well as viruses in rapid decline of established ‘HoneyCrisp’ trees grafted onto the rootstock M.9 (clone NIC 29). Analysis of weather variables indicated abnormal trends at the experimental orchard site over the five-year growth period with December 2014, January, February and March 2015 being exceptionally cold. According to National Climate Report from the National Oceanic and Atmosphoric Administration (NOAA) National Centers for Environmental Information of the US, February 2015 was the third-coldest February on record in the region since 1934. In 2015, there were heavy rains between June-July followed by a severe drought in 2016 during the apple growing season. A severe temperature drop followed by a moderately warmer winter may have caused direct damage to apple trees or indirectly made them more susceptible to biotic and abiotic stresses [8, 9]. Both significant loss of trees and decrease in yield have previously occurred due to freezing temperatures in the United States and Canada [7, 42–45]. Cold injury symptoms become apparent in the spring following a hard winter with trees exhibiting stunted growth, wilting, and death. Rainfall data also indicated one drought instance in the area in 2016. Although no study has linked RAD with drought, general water stress can cause severe damage by reducing tree growth, root damage, and senescence [12–16, 46]. The occurrence of severe cold followed by drought, or either individually, might not directly cause RAD, but could have weakened the trees and led to the proliferation of insects and infection by opportunistic pathogens. Although the general soil fertility status was not optimal for commercial apple production [47], differences between soil nutrition profiles from symptomatic and asymptomatic apple trees were not significant, ruling out any potential role of nutritional stress in RAD in the experimental orchard. The concentrations of organic matter, total nitrogen, phosphorus, potassium, magnesium, and boron were all low in the experimental orchard area. Meanwhile, sulfur concentration was nearly two-fold over the recommended amount [47]. Nutritional differences might not directly contribute, but can exacerbate the impact of other stress factors such as weather extremes, insects, and pathogens, therefore increasing the chances of decline of apple trees. In some rootstock and scion combinations, a weak graft union can be impacted under extreme abiotic stresses, leading to a slow collapse of the tree. Also, rootstocks with shallow root-systems in high density plantings may have limited access to nutrients and water and thus be unable to support the heavy crop, foliage, and biomass under extreme weather, leading to decline and death of trees. The negative effect can be exacerbated in soils with poor water holding capacity. The presence of necrotic lesions and wood decay were among the potential signs of declining apple trees in the orchard. These symptoms followed a specific dispersion pattern from bark or vascular cambium towards the heartwood. Also, the wood decay and discoloration mainly occurred below the graft union and progressed downwards. These observations indicated that the rootstock seems to be the originating point of the decline. Previous reports have implicated incompatibility between rootstock and scion as a potential trigger for RAD [2]. Poor vascular connections, phloem degeneration, and vascular discontinuity [48] can cause such graft incompatibility. However, the pattern of wood rot and discoloration in the trees examined in this study appears to be different from that seen with scion-rootstock incompatibility, which is generally equally dispersed above and below the graft union. The exact cause of necrotic lesions and discoloration only in rootstock wood is yet unclear, but it could involve multiple pathogens in RAD. Some microbes can also cause root rots and internal clotting of vascular tissues that lead to tree decline and ultimately death [22, 49, 50]. Rootstock fire blight caused by Erwinia amylovora and root or crown rot caused by Phytophthora and Pythium species can lead to tree death by infecting the root system. However, we did not identify symptoms or pathogens of common soil-borne root diseases or apple rootstock blight in any of the declining ‘HoneyCrisp’ trees. Apple pathogenic viruses can potentially result in decline of trees with similar symptoms as observed for RAD. Tree decline from latent viruses, usually transmitted by grafting and top working, is an emerging commercial problem [4]. There were no differences between healthy-looking and declined plants for the presence of the stem pitting, stem grooving, apple chlorotic leaf spot, or apple mosaic, the most common latent viruses of apples. Tomato ringspot virus, which can be involved in apple tree decline and death, was not found in any of the declining ‘HoneyCrisp’ trees. Similarly, tobacco ringspot virus was not identified in this study. Fungal pathogens cause several soil-borne diseases in tree crops, resulting in symptoms similar to those observed in declining apple trees. For example, a vascular fungus, Ceratocystis fimbriate can kill mango trees two months after the initial infection, causing mango sudden decline or mango wilt [51]. In our study, several classes of bacterial and fungal communities were present in root, shoot, soil, and rhizosphere samples from apple trees. Similar bacterial and fungal classes have been detected previously in different apple tissues and soil types [38, 52–53]. For instance, both Basidiomycota and Ascomycota fungus classes were detected in microbiomes of the shoots and soil samples in apple trees [38, 53]. Similarly, bacterial classes related to proteobacteria, actinobacteria, and acidobacteria were also detected in these studies. These probably represent the most common bacterial and fungal classes in apple orchards from New York soils. However, future studies can more clearly define the enrichment of these microbiome populations in the apple cultivation area. Differential abundance analysis of ITS sequences of shoot, root, soil, and rhizosphere samples showed no significant enrichment of any fungal species in RAD affected apple trees. However, 16S sequencing identified a single bacterial class, ‘alphaproteobacteria-rickettsiales’, with differential abundance in the roots of healthy-looking and declined apple trees. These bacteria types were less abundant in the declined samples. The relative abundance of proteobacteria usually declines in water-limiting soils [54–56], suggesting a difference in the moisture level in roots of healthy-looking and declined trees. However, it is difficult to determine at this point whether roots from declining apple trees endured low moisture at the time of sampling or low moisture resulted from the severe drought in 2016. Nonetheless, these observations indicate the role of water-limiting conditions towards the rapid decline of apple trees. This can provide a focus for future RAD research to further explore the roles of soil and plant water status, taking into consideration the drought tolerance of different rootstocks and their susceptibility to rapid decline. Moreover, differences in root system architecture of the rootstocks can also be studied under extreme water stress conditions in different soil types for its role in RAD. For instance, shallow rooted rootstocks can be at a disadvantage in non-irrigated orchards with a low water table, whereas rootstocks with deeper and more vigorous roots may be more tolerant to drought stress conditions. Parasitic nematodes such as root lesion, root knot, and dagger nematodes can also infect the root system of apple trees, resulting in leaf chlorosis, stunted growth of trees and poor yields. However, we did not evaluate the presence of pathogenic nematodes in RAD affected orchards. RAD symptoms are distinct from those of apple replant disease. A complex of multiple species of fungi and oomycetes including Rhizoctonia, Phytophthora, and Pythium species, and root lesion nematodes causes apple replant disease [31], characterized by reduced productivity in orchards repeatedly planted with the same or closely related fruit trees. Symptoms of replant disease, including uneven or stunted growth of trees with short internodes, are visible shortly after planting new trees. When a tree is uprooted, discolored roots, root tip necrosis, and reduced root biomass can be seen. Young trees may die within the first year. Many will survive but overall fruit production and quality are reduced [57]. In conclusion, we did not find any statistically significant differences in soil and weather profiles of healthy-looking and declined trees. Similarly, no particular fungi and viruses were associated with the symptomatic trees. A single class of proteobacteria showed differential abundance between symptomatic and asymptomatic samples, suggesting a possible role of water-limiting conditions. Similarly, the role of different opportunistic or previously unknown pathogens should not be excluded. We speculate that the onset of RAD symptoms is much later than the actual cause of the decline. The findings in this study, however, will require further validation in different declining orchards that have diverse scion-rootstock combinations and different soil types and weather conditions.

Acknowledgments We thank David Strickland for collecting data and Dr. Ali Srour from University of Southern Illinois, Carbondale for discussions about 16S and ITS sampling and analysis.