Bacterial Community and Diversity from the Watermelon Cultivated Soils through Next Generation Sequencing Approach

Article information

Plant Pathol J. 2021;37(6):521-532
Publication date (electronic) : 2021 December 1
doi : https://doi.org/10.5423/PPJ.OA.07.2021.0106
Division of Biological Resource Sciences, Interdisciplinary Program in Smart Agriculture, Kangwon National University, Chuncheon 24341, Korea
*Corresponding author: Phone) +82-33-250-6417, FAX) +82-33-259-5558, E-mail) younslee@kangwon.ac.kr
Handling Editor: Kihyuck Choi
Received 2021 July 4; Revised 2021 September 13; Accepted 2021 September 17.

Abstract

Knowledge and better understanding of functions of the microbial community are pivotal for crop management. This study was conducted to study bacterial structures including Acidovorax species community structures and diversity from the watermelon cultivated soils in different regions of South Korea. In this study, soil samples were collected from watermelon cultivation areas from various places of South Korea and microbiome analysis was performed to analyze bacterial communities including Acidovorax species community. Next generation sequencing (NGS) was performed by extracting genomic DNA from 92 soil samples from 8 different provinces using a fast genomic DNA extraction kit. NGS data analysis results revealed that, total, 39,367 operational taxonomic unit (OTU), were obtained. NGS data results revealed that, most dominant phylum in all the soil samples was Proteobacteria (37.3%). In addition, most abundant genus was Acidobacterium (1.8%) in all the samples. In order to analyze species diversity among the collected soil samples, OTUs, community diversity, and Shannon index were measured. Shannon (9.297) and inverse Simpson (0.996) were found to have the highest diversity scores in the greenhouse soil sample of Gyeonggi-do province (GG4). Results from NGS sequencing suggest that, most of the soil samples consists of similar trend of bacterial community and diversity. Environmental factors play a key role in shaping the bacterial community and diversity. In order to address this statement, further correlation analysis between soil physical and chemical parameters with dominant bacterial community will be carried out to observe their interactions.

Microbial communities are group of microorganisms that are inhabited in soil sharing a common space. Soil microbes may directly or indirectly affect the soil properties. Bacteria, actinomycetes, fungi, algae, and protozoa are the microbes found diverse in soil (Subba Rao, 1999). Soil bacteria could fix nitrogen as a result released oxygen into the atmosphere (Canfield, 2014; Subba Rao, 1999) which affects soil structure and fertility. Soil microbes has vital role in the soil environment. Change of soil environment is often linked with the variation of microbial populations (Kuang et al., 2018). Microbe’s diversity in soil plays a major role in maintaining healthy soil because these microbes are engaged in many vital functions like crucial cycles of C, N, P, formation of soil, removal of toxin, and so on (Fatima et al., 2014). Microbial composition and diversity are hugely known as prime factors in processing ecological activities (Bauman et al., 2013; Hernandez-Raquet et al., 2013). Numerous studies have found significant effects of biotic and abiotic factors on underground microbial composition and land management (Jangid et al., 2008).

Recently several studies have been focused in soil microbial community to modify and better use to promote plant health (Li et al., 2017). Soil microbiome analysis can predetermine whether plants are resistant or susceptible to diseases. Variation in soil microbiome composition and functioning can determine the reaction of plant-pathogen interaction under natural field conditions (Wei et al., 2019). Watermelon (Citrullus lanatus) is a major cucurbitaceous crop widely grown in tropical and subtropical regions across the world (Chomicki and Renner, 2015). Recently, plants have given account a holobiont, a unit representing both the host and its associated microbiome (Vandenkoornhuyse et al., 2015). Particularly, most roots and stem colonizing microorganisms seem to have originated from the rhizosphere while some are able to move systematically within the plant (Rosenblueth and Martinez-Romero, 2006). However, some microorganisms are restricted to belowground parts of plants (Compant et al., 2011). Such, plant-associated microbiome is highly diverse and holds a wide range of taxa. In order to study those microorganisms associated with plant parts and soil, composition and diversity analysis of those microorganisms is utmost.

Previous studies on the diversity of microbial communities rely on the isolation from soil sample using media selective technique for their identification morphologically. These media-based techniques for microbial community’s diversity study have been lead to biased method due to the unculturability of many microbes. Some soil microbes need specific media to grow in peculiar conditions. It has been reported that only 1% of soil microbes can be cultured in the laboratory using traditional culture techniques (Marchesi et al., 1998). In order to address these issues, microbiologists have adopted culture-independent techniques which include molecular biology based methods which consists of extracted DNA from soil samples is subjected to PCR amplification (Roh et al., 2006). These methods provide the detailed insights of species richness, evenness, and composition of bacterial community.

In order to provide insights in soil microbial activity and diversity, soil microbiome is the recently used technique which provides soil microbial community structure by using next generation sequencing (NGS) technology. This sequencing technology has shown a step way change in facilitating the characterization of soil prokaryotic and eukaryotic microbes (Bates et al., 2013; Fierer and Jackson, 2006; Ramirez et al., 2018; Tedersoo et al., 2014; Wu et al., 2011). Such sequencing technology for microbial community analysis uses ribosomal gene fragments from shotgun metagenomics (Guo, 2016). Microbiome research has crucial applications for understanding and manipulating ecosystem processes such as organic matter turnover, nutrient cycling, and the development or inhibition of soil pathogens (Philippot et al., 2013). Research in microbiome area has great ability to know the impacts of soil microbes on the productivity of natural plant communities and agroecosystems (McPherson et al., 2018).

The information regarding the bacterial community composition and diversity of watermelon cultivated soils in different geographical locations in Korea is limited and scarce. In order to address this, we have conducted this study with the following objectives: (1) to provide the insight information to researchers to describe the bacterial diversity present in soils in different geographical locations where the watermelon is cultivated, (2) to understand the soil surface ecosystem response to bacterial communities in varying environmental conditions.

Materials and Methods

Field site description and soil sampling

The study was conducted on August 2020 in South Korea. Detailed geographical coordinates and sampling locations have been presented in Table 1 and Supplementary Fig. 1. A total of 92 soil samples were taken from the use of soil auger (0–15 cm depth) from watermelon cultivated soils. Soil samples were collected in zipper bag with proper labeling. Each sample was placed in a plastic zipper bag and thoroughly homogenized before analysis. In the laboratory, the samples were sieved through a 2-mm mesh to remove roots and other debris and subsequently divided into two parts. Soil samples were kept at −20C for the subsequent DNA extraction and molecular analysis of community composition.

Soil samples collection sites with their GPS coordinates

Soil DNA extraction

Soil DNA extraction was carried out from the watermelon cultivated soils to sterilized buffer (NaCl, 4.25 g; KH2PO4, 0.15 g; Na2HPO4, 0.3 g; MgSO4 0.1 g; gelatin, 0.05 g), mixed, and allowed to stand for 10 minutes. Soil solution was extracted using Fast DNA Spin Kit for soil (MP, Biomedicals, Seoul, Korea): Fast DNA SPIN Kit for soil according to the extraction method. A quantitative polymerase chain reaction (qPCR) amplification was conducted after the total DNA extraction from the soil samples. A qPCR amplification of the prokaryotic microbes (16S rRNA genes V4 region) was performed using the forward primer 515F (5′-GTGYCAGCMGCCGCGGTAA-3′) and the reverse primer 806R (5′-GGACTACHVGGGTWTCTAAT-3′) (Bates et al., 2011). Amplicons were sequenced on Illumina-MiSeq platform (Illumina Inc., San Diego, CA, USA) at SHBIO Technology (Shanghai, China).

Analysis of sequencing data for bacteria community composition

Raw sequences received after PCR amplicon sequencing were analyzed and processed mainly in three steps: data trimming, taxonomy assignment and normalization, and diversity analysis. Briefly, raw sequences were filtered using FastQC run using the software fastqc v0.11.2 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc). Forward and reverse direction of fastQC files were assimilated from software vsearch v2.10.3 (https://www.github.com/torognes/vsearch), followed by filtering of the chimeric sequences. This process leads to the removal of unassembled and low-quality sequences with >97% similarity could be clustered into operational taxonomic units (OTU). OTU and the representative sequences for each OTU were constructed using QIIME (Quantitative Insights into Molecular Ecology) software (version 1.17) package (Jiang et al., 2017) to annotate the taxonomic information for each representative sequence.

Analysis of soil bacterial diversity

Alpha and beta diversity analyses were carried out based on the normalized OTU data. Alpha diversity indices, including Chao1, Shannon, and Simpson indices, were calculated with QIIME (version 1.7.0) (Caporaso et al., 2010) to assess the richness and species diversity for each soil sample of different geographical locations. The beta-diversity patterns of the bacterial communities (Bray-Curtis dissimilarities) (Sokal and Michener, 1958) were characterized using principal coordinate analysis (PCoA) and the unweighted pair-group method with arithmetic mean (UPGMA) clustering based on unweighted UniFrac distances. The unweighted Uni-Frac (UUF, qualitative measure) assesses the presence or absence of taxa and is, thus, more sensitive to rare taxa.

Results and Discussion

Bacterial community composition across the different geographical locations

In order to find the bacterial fruit blotch causing pathogen, Acidovorax citrulli from the soil samples collected from different locations cultivated with watermelon plants NGS of the hypervariable V1 and V4 regions of bacterial 16S ribosomal RNA was performed. Results showed diverse bacterial community composition in all the locations. Abundancy analysis results showed that, 6 (Proteobacteria, Acidobacteria, Actinobacteria, Bacteroidetes, Gemmatimonadetes, and Verrucomicrobia) bacterial phylum dominantly presents in all of the 92 soil samples (Fig. 1). However, Proteobacteria, Bacteroidetes, Acidobacteria, Gemmatimonadetes, Verrucomicrobia, Actinobacteria, Chlorofexi, Firmicutes, Ignavibacteriae, Planctomycetes, and Nitrospirae were the dominant phylum detected in two different land use type soils (open field and greenhouses) from all the provinces. Among these, Proteobacteria was the most abundantly dominant bacterial phylum in all the soil samples. Soil samples from Gangwon-do, comprise about 50%, 17%, and 10% of Proteobacteria, Bacteroidetes, and Acidobacteria, respectively (Supplementary Fig. 2). Soil samples from Gyeongsangbuk-do also mostly dominated by Proteobacteria, Acidobacteria, and Bacteroidetes with around 44%, 27%, and 22%, respectively (Supplementary Fig. 3). In case of soil samples from Gyeongsangnam-do, Proteobacteria, Actinobacteria and Firmicutes were the most dominant in the ratio about 48%, 35%, and 14%, respectively (Supplementary Fig. 4). Dominant phylum coverage of the bacteria from the soil of Jeollabuk-do were Proteobacteria (32–51%), Bacteroidetes (11–22%), and Acidobacteria (8–10%) (Supplementary Fig. 5). Soil samples from Jeollanam-do hold Proteobacteria, Acidobacteria, and Bacteroidetes as the most dominant bacteria phylum with the abundancy ratio of about 37%, 22%, and 14%, respectively (Supplementary Fig. 6). Abundance ratios of around 63%, 54%, and 10% of Bacteroidetes, Proteobacteria, and Acidobacteria was obtained from soil samples of Chungcheongbuk-do (Supplementary Fig. 7). Soil samples from Chungcheongnam-do were highly accounted for Proteobacteria, Bacteroidetes, and Acidobacteria with the dominant ratios of around 43%, 20%, and 18%, respectively (Supplementary Fig. 8). Moreover, soil samples from Gyeonggi-do consisted of about 46%, 20%, and 10% of Proteobacteria, Bacteroidetes, and Actinobacteria, respectively (Supplementary Fig. 9). Monitoring soil bacterial communities in soils from different geographical conditions is a major challenge. It helps implement a reliable sustainable management of soils in agriculture (Bender et al., 2016). Effective knowledge regarding the geographical distribution of microbial communities is critical to better understand the ways soils under different climate conditions (Plassart et al., 2019). In addition, pattern of bacterial diversity is markedly different from that of plant species and most soil bacteria and may alter with environmental and geographical distance (Wang et al., 2015). However, primary factors affecting microbial spatial distribution are still under debate (Lauber et al., 2009). The distribution of soil bacterial community can be altered by various factors like; soil properties, litter quality, root exudates, and also by temperature and precipitation in different geographical conditions (Nielsen et al., 2010). In this study, we observed slight variation in dominant bacterial composition across the different geographical locations (Fig. 1). However, factors affecting such variations in bacterial composition in different geographical locations need to be studied further.

Fig. 1

Top most dominant phylum along with the abundance ratio detected across all of the 92 soil samples. GW, Gangwon-do; GB, Gyeongsangbuk-do; KN, Gyeongsangnam-do; JB, Jeollabuk-do; JN, Jeollanam-do; CB, Chungcheongbuk-do; CN, Chungcheongnamdo; GG, Gyeonggi-do.

Moreover, in this study, we are also trying to assess whether our collected watermelon cultivated soil samples consist A. citrulli or not from various geographical locations. In this regard, we could not observe any A. citrulli isolate in our studied soil samples as this bacterial pathogen is a seed-borne disease affecting the embryo and endosperm of seeds (Rahimi-Midani et al., 2018). Furthermore, previous reports suggest that, A. citrulli actively colonize and move through xylem vessels then reach fruits and seeds (Bahar et al., 2009) and the chance of moving this pathogen into soil is very rare except by the other external influencing factor. Not detection of A. citrulli in our studied soil samples might be due to the above mentioned reasons.

Bacterial diversity across the different geographical locations

Bacterial diversity was analyzed in soils of different provinces. Alpha diversity was assessed using Chao1, Shannon indices, and beta diversity was analyzed using PCoA. The OTU in all the analyzed soil samples was in range of 276 to 2216. OTU is an operational definition used to classify groups of closely related individuals and ‘OTU’ is simply the group of organisms currently being studied. The number of OTUs in this study was found lowest from the greenhouse soil of Jeollabuk-do (JB4) and highest from greenhouse soil of Gyeonggi-do (GG4) (Fig. 2). Bacterial diversity according in land use types also was investigated in this study. Twenty-six soil samples were from open-field (GB1-GB22, JB1 and JB3, JN1-JN3) and 66 were from greenhouse (GW1-GW10, KN1-KN10, CB1–14, CN1-CN17, and GG1-GG12). OTUs from open-field soil was recorded highest in GB12 and lowest in JB1 (Supplementary Fig. 10). Soils from greenhouse of GG4 was the highest OTUs as compared with other greenhouse soils (Supplementary Fig. 11). For several years, OTUs have been the most commonly used units of diversity, especially when analyzing small subunit 16S for prokaryotes marker gene sequence datasets (Sommer et al., 2017). In order to relate this with our study, 16S rRNA metagenomics sequencing revealed the JB4 with the lowest OTUs and GG4 with the highest OTUs. In total soil samples, Shannon and inverse Simpson diversity indices were found lower in JB4 (Shannon, 3.53; Inverse Simpson, 0.73) and highest in GG4 (Shannon, 9.29; Inverse Simpson, 0.995) (Fig. 3). Moreover, diversity indices (Shannon and inverse Simpson) according to land use types were also observed. Results suggest that, among the open field soils, diversity was higher in GB12 and lower in JB1 (Supplementary Fig. 12). Whereas, among the greenhouse soils, diversity was higher in GG4 and lower in JB1 (Supplementary Fig. 13). The Shannon diversity index represents the microbial community richness, indicating the more abundant species in the community (Hill et al., 2003). Microbial diversity has been closely associated to soil health. Microbial community diversity enhances the well-functioning of the soil ecosystem. Rao 1999 reported that high microbial diversity is an indicator of a healthy soil, while low microbial diversity indicates a stressed condition of soil. Moreover, the microbial diversity of soil is the most crucial component that reflects the quality and stability of soil (Cassán and Diaz-Zorita, 2016) and has the ability to set various non-biological indicators. Hence, the exploration of soil microbial diversity is critical factor for the estimation of health and fertility of soil using high-throughput DNA sequencing technology (Roesch et al., 2007). In this study, however, all the soils samples were from watermelon cultivated soils, bacterial community diversity were found more diversified. Bacterial species richness and diversity (OTUs, Chao1, Shannon, Inverse Simpson) was assessed lowest in soil of Jeollabukdo (JB4) and highest in soil of Gyeonggi-do (GG4).

Fig. 2

Operational taxonomic units (OTUs) of bacterial community obtained from total 92 soil samples from watermelon cultivated soil across the different provinces. X-axis indicates soil samples from various geographical locations: GW, Gangwon-do; KB and GB, Gyeongsangbuk-do; KN, Gyeongsangnam-do; JB, Jeollabuk-do; JN, Jeollanam-do; CB, Chungcheongbuk-do; CN, Chungcheongnamdo; GG, Gyeonggi-do.

Fig. 3

Bacterial community richness and diversity obtained from total 92 soil samples from watermelon cultivated soils across the different provinces. X-axis indicates soil samples from various geographical locations: GW, Gangwon-do; KB and GB, Gyeongsangbuk-do; KN, Gyeongsangnam-do; JB, Jeollabuk-do; JN, Jeollanam-do; CB, Chungcheongbuk-do; CN, Chungcheongnam-do; GG, Gyeonggi-do.

Variation in bacterial community structure in different geographical locations

UPGMA and PCoA plot was used to visualize the variation in bacterial community structure in watermelon cultivated soils across the different geographical locations. The unweighted UniFrac (UUF, qualitative measure) provides the status of taxa which is more sensitive to rare taxa. The UPGMA tree in our study revealed that, bacterial communities from soil samples of same province lie in nearby clusters (Fig. 4). Qualitative measurement of bacterial community among different provinces mostly exhibited separate clade in UPGMA tree (Fig. 5). However, some bacterial community taxa of soil samples of same province showed slight variation as seen in the UPGMA tree (Fig. 5). Interestingly, unweighted and weighted unifrac tree showed that, only bacterial community from the soils of Gyeongsangnam-do (KN1-KN10) lies in same clade (Figs. 4 and 5). In contrast, bacterial community from the soils of Chungcheongbuk-do (CB1-CB14) were lying in different clade in both unweighted and weighted unifrac tree (Figs. 4 and 5). Measures of variation in microbial community composition in different soils were broadly performed by β-diversity. Examining the patterns of β-diversity (similarities/dissimilarities) across different land uses soils provides insights into the dynamics of microbial communities (Anderson et al., 2011). In our study, we have tried to examine the patterns of bacterial community composition in watermelon cultivated soils in different locations in South Korea. PCoA plot of weighted and unweighted unifrac tree exhibited significant variation in bacterial community structure among the land use types across different geographical locations (Figs. 6 and 7). Our study distinctly showed that the dominant bacterial community composition was slightly varied among the locations. Index for each bacterial community composition in respective geographical locations has been indicated in the PCoA graph of Figs. 6 and 7.

Fig. 4

Weighted unifrac tree showing hierarchical clustering of average linkage of bacterial community from 92 soil samples.

Fig. 5

Unweighted unifrac tree showing hierarchical clustering of average linkage of bacterial community from 92 soil samples.

Fig. 6

Principal coordinate analysis (PCoA) plot of weighted unifrac showing variation based on the group of the bacterial communities encountered across the soil samples. Different colored plots indicate the variation in bacterial community composition according to land use types and geographical locations.

Fig. 7

Principal coordinate analysis (PCoA) plot of unweighted unifrac showing variation of the bacterial communities encountered based on the group across the soil samples. Different colored plots indicate the variation in bacterial community composition according to land use types and geographical locations.

Conclusion

The highest bacterial diversity was observed in the greenhouse soil from Gyeonggi-do province (GG4) and the lowest bacterial diversity was detected in the greenhouse soil from Jeollabuk-do province (JB1). We have found slight variation in bacterial diversity in open field soils and greenhouse soils of different provinces. However, bacterial community composition and diversity was somewhat similar from the soil samples of same land use types. Significant reason behind the variation in bacterial diversity is might be due to the soils from different locations and land use types. Alpha diversity indices (OTUs, Shannon and Inverse Simpson) and beta diversity indices (UPGMA tree, dissimilarity) are the key factors to determine overall diversity of microbes present in the soil. In this study, we observed watermelon cultivated soil samples from nearby area holds similar bacterial community composition and diversity as compare to different areas through UPGMA tree using distance matrix values. In contrast, we could not find the A. citrulli from our soil samples and further soil samples and watermelon plant samples will be analyzed in order to find the A. citrulli using NGS approach. However, further studies regarding the correlation of chemical attributes shaping the bacterial community in different land use type are further needed to draw the general trend. Thus, to the best of our knowledge, this is the first report regarding bacterial community composition and diversity from the watermelon cultivated soil in various geographical locations in South Korea. Also, we expect that the soil microbiome analysis could provide some biological indicators for the assessment of soil health and ecosystem services, but for this purpose, further studies with a holistic approach would be worthwhile.

Acknowledgments

This work was supported by Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, and Forestry (IPET) through Project No. 120088-05-1-SB010 funded by Ministry of Agriculture, Food and Rural Affairs (MAFRA). This research has been worked partially with the support of a research grant of Kangwon National University in 2019 and 2020.

Notes

Conflict of interest

No potential conflict of interest relevant to this article was reported.

Electronic Supplementary Material

Supplementary materials are available at The Plant Pathology Journal website (http://www.ppjonline.org/).

References

Anderson MJ, Crist TO, Chase JM, Vellend M, Inouye BD, Freestone AL, Sanders NJ, Cornell HV, Comita LS, Davies KF, Harrison SP, Kraft NJB, Stegen JC, Swenson NG. 2011;Navigating the multiple meanings of β diversity: a roadmap for the practicing ecologist. Ecol Lett 14:19–28.
Bahar O, Goffer T, Burdman S. 2009;Type IV Pili are required for virulence, twitching motility, and biofilm formation of Acidovorax avenae subsp. citrulli . Mol Plant-Microbe Interact 22:909–920.
Bates ST, Berg-Lyons D, Caporaso JG, Walters WA, Knight R, Fierer N. 2011;Examining the global distribution of dominant archaeal populations in soil. ISME J 5:908–917.
Bates ST, Clemente JC, Flores GE, Walters WA, Parfrey LW, Knight R, Fierer N. 2013;Global biogeography of highly diverse protistan communities in soil. ISME J 7:652–659.
Baumann K, Dignac M-F, Rumpel C, Bardoux G, Sarr A, Steffens M, Maron PA. 2013;Soil microbial diversity affects soil organic matter decomposition in a silty grassland soil. Biogeochemistry 114:201–212.
Bender SF, Wagg C, van der Heijden MGA. 2016;An underground revolution: biodiversity and soil ecological engineering for agricultural sustainability. Trends Ecol Evol 31:440–452.
Canfield DE. 2014. Oxygen: a four billion year history Princeton University Press. Princeton, NJ, USA: p. 216.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, Mc-Donald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. 2010;QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336.
Cassán F, Diaz-Zorita M. 2016; Azospirillium sp. in current agriculture: from the laboratory to the field. Soil Biol Biochem 103:117–130.
Chomicki G, Renner SS. 2015;Watermelon origin solved with molecular phylogenetics including Linnaean material: another example of museomics. New Phytol 205:526–532.
Compant S, Mitter B, Colli-Mull JG, Gangl H, Sessitsch A. 2011;Endophytes of grapevine flowers, berries, and seeds: identification of cultivable bacteria, comparison with other plant parts, and visualization of niches of colonization. Microb Ecol 62:188–197.
Fatima F, Pathak N, Rastogi Verma S. 2014;An improved method for soil DNA extraction to study the microbial assortment within rhizospheric region. Mol Biol Int 2014:518960.
Fierer N, Jackson RB. 2006;The diversity and biogeography of soil bacterial communities. Proc Natl Acad Sci U S A 103:626–631.
Guo J. 2016. Rhizosphere metagenomics of three biofuel crops. PhD thesis Michigan State University; East Lansing, MI, USA: 108.
Hernandez-Raquet G, Durand E, Braun F, Cravo-Laureau C, Godon J-J. 2013;Impact of microbial diversity depletion on xenobiotic degradation by sewage-activated sludge. Environ Microbiol Rep 5:588–594.
Hill TCJ, Walsh KA, Harris JA, Moffett BF. 2003;Using ecological diversity measures with bacterial communities. FEMS Microbiol Ecol 43:1–11.
Jangid K, Williams MA, Franzluebbers AJ, Sanderling JS, Reeves JH, Jenkins MB, Endale DM, Coleman DC, Whitman WB. 2008;Relative impacts of land-use, management intensity and fertilization upon soil microbial community structure in agricultural systems. Soil Biol Biochem 40:2843–2853.
Jiang J, Song Z, Yang X, Mao Z, Nie X, Guo H, Peng X. 2017;Microbial community analysis of apple rhizosphere around Bohai Gulf. Sci Rep 7:88918.
Kuang S, Su Y, Wang H, Yu W, Lang Q, Matangi R. 2018;Soil microbial community structure and diversity around the aging oil sludge in Yellow River Delta as determined by high-throughput sequencing. Archae 2018:7861805.
Lauber CL, Hamady M, Knight R, Fierer N. 2009;Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl Environ Microbiol 75:5111–5120.
Li S, Peng M, Liu Z, Shah SS. 2017;The role of soil microbes in promoting plant growth. Mol Microbiol Res 7:30–37.
Marchesi JR, Sato T, Weightman AJ, Martin TA, Fry JC, Hiom SJ, Dymock D, Wade WG. 1998;Design and evaluation of useful bacterium-specific PCR primers that amplify genes coding for bacterial 16S rRNA. Appl Environ Microbiol 6:795–799.
McPherson MR, Wang P, Marsh EL, Mitchell RB, Schachtman DP. 2018;Isolation and analysis of microbial communities in soil, rhizosphere, and roots in perennial grass experiments. J Vis Exp 137:57932.
Nielsen UN, Osler GHR, Campbell CD, Burslem DFRP, van der Wal R. 2010;The influence of vegetation type, soil properties and precipitation on the composition of soil mite and microbial communities at the landscape scale. J Biogeogr 37:1317–1328.
Philippot L, Spor A, Hénault C, Bru D, Bizouard F, Jones CM, Sarr A, Maron P-A. 2013;Loss in microbial diversity affects nitrogen cycling in soil. ISME J 7:1609–1619.
Plassart P, Prévost-Bouré NC, Uroz S, Dequiedt S, Stone D, Creamer R, Griffiths RI, Bailey MJ, Ranjard L, Lemanceau P. 2019;Soil parameters, land use, and geographical distance drive soil bacterial communities along a European transect. Sci Rep 9:605.
Rahimi-Midani A, Lee YS, Kang S-W, Kim M-K, Choi T-J. 2018;First isolation and molecular characterization of bacteriophages infecting Acidovorax citrulli, the causal agent of bacterial fruit blotch. Plant Pathol J 34:59–64.
Ramirez KS, Knight CG, de Hollander M, Brearley FQ, Constantinides B, Cotton A, Creer S, Crowther TW, Davison J, Delgado-Baquerizo M, Dorrepaal E, Elliott DR, Fox G, Griffiths RI, Hale C, Hartman K, Houlden A, Jones DL, Krab EJ, Maestre FT, McGuire KL, Monteux S, Orr CH, van der Putten WH, Roberts IS, Robinson DA, Rocca JD, Rowntree J, Schlaeppi K, Shepherd M, Singh BK, Straathof AL, Bhatnagar JM, Thion C, van der Heijden MGA, de Vries FT. 2018;Detecting macroecological patterns in bacterial communities across independent studies of global soils. Nat Microbiol 3:189–196.
Roesch LFW, Fulthorpe RR, Riva A, Casella G, Hadwin AKM, Kent AD, Daroub SH, Camargo FAO, Farmerie WG, Triplett EW. 2007;Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J 1:283–290.
Roh C, Villatte F, Kim B-G, Schmid RD. 2006;Comparative study of methods for extraction and purification of environmental DNA from soil and sludge samples. Appl Biochem Biotechnol 134:97–112.
Rosenblueth M, Martínez-Romero E. 2006;Bacterial endophytes and their interactions with hosts. Mol Plant-Microbe Interact 19:827–837.
Sokal RR, Michener CD. 1958;A statistical method for evaluating systematic relationships. Univ Kansas Sci Bull 38:1409–1438.
Sommer SA, Van Woudenberg L, Lenz PH, Cepeda G, Goetze E. 2017;Vertical gradients in species richness and community composition across the twilight zone in the North Pacific Subtropical Gyre. Mol Ecol 21:6136–6156.
Subba Rao NS. 1999. Soil microbiology 4th edth ed. Science Publishers. Enfield, NH, USA: p. 407.
Tedersoo L, Bahram M, Põlme S, Kõljalg U, Yorou NS, Wijesundera R, Villarreal Ruiz L, Vasco-Palacios AM, Thu PQ, Suija A, Smith ME, Sharp C, Saluveer E, Saitta A, Rosas M, Riit T, Ratkowsky D, Pritsch K, Põldmaa K, Piepenbring M, Phosri C, Peterson M, Parts K, Pärtel K, Otsing E, Nouhra E, Njouonkou AL, Nilsson RH, Morgado LN, Mayor J, May TW, Majuakim L, Lodge DJ, Lee SS, Larsson KH, Kohout P, Hosaka K, Hiiesalu I, Henkel TW, Harend H, Guo LD, Greslebin A, Grelet G, Geml J, Gates G, Dunstan W, Dunk C, Drenkhan R, Dearnaley J, De Kesel A, Dang T, Chen X, Buegger F, Brearley FQ, Bonito G, Anslan S, Abell S, Abarenkov K. 2014;Fungal biogeography. Global diversity and geography of soil fungi. Science 346:1256688.
Vandenkoornhuyse P, Quaiser A, Duhamel M, Le Van A, Dufresne A. 2015;The importance of the microbiome of the plant holobiont. New Phytol 206:1196–1206.
Wang X, Van Nostrand JD, Deng Y, Lü X, Wang C, Zhou J, Han X. 2015;Scale-dependent effects of climate and geographic distance on bacterial diversity patterns across northern China’s grasslands. FEMS Microbiol Ecol 91:fiv133.
Wei Z, Gu Y, Friman V-P, Kowalchuk GA, Xu Y, Shen Q, Jousset A. 2019;Initial soil microbiome composition and functioning predetermine future plant health. Sci Adv 5:eaaw0759.
Wu T, Ayres E, Bardgett RD, Wall DH, Garey JR. 2011;Molecular study of worldwide distribution and diversity of soil animals. Proc Natl Acad Sci U S A 108:17720–17725.

Article information Continued

Fig. 1

Top most dominant phylum along with the abundance ratio detected across all of the 92 soil samples. GW, Gangwon-do; GB, Gyeongsangbuk-do; KN, Gyeongsangnam-do; JB, Jeollabuk-do; JN, Jeollanam-do; CB, Chungcheongbuk-do; CN, Chungcheongnamdo; GG, Gyeonggi-do.

Fig. 2

Operational taxonomic units (OTUs) of bacterial community obtained from total 92 soil samples from watermelon cultivated soil across the different provinces. X-axis indicates soil samples from various geographical locations: GW, Gangwon-do; KB and GB, Gyeongsangbuk-do; KN, Gyeongsangnam-do; JB, Jeollabuk-do; JN, Jeollanam-do; CB, Chungcheongbuk-do; CN, Chungcheongnamdo; GG, Gyeonggi-do.

Fig. 3

Bacterial community richness and diversity obtained from total 92 soil samples from watermelon cultivated soils across the different provinces. X-axis indicates soil samples from various geographical locations: GW, Gangwon-do; KB and GB, Gyeongsangbuk-do; KN, Gyeongsangnam-do; JB, Jeollabuk-do; JN, Jeollanam-do; CB, Chungcheongbuk-do; CN, Chungcheongnam-do; GG, Gyeonggi-do.

Fig. 4

Weighted unifrac tree showing hierarchical clustering of average linkage of bacterial community from 92 soil samples.

Fig. 5

Unweighted unifrac tree showing hierarchical clustering of average linkage of bacterial community from 92 soil samples.

Fig. 6

Principal coordinate analysis (PCoA) plot of weighted unifrac showing variation based on the group of the bacterial communities encountered across the soil samples. Different colored plots indicate the variation in bacterial community composition according to land use types and geographical locations.

Fig. 7

Principal coordinate analysis (PCoA) plot of unweighted unifrac showing variation of the bacterial communities encountered based on the group across the soil samples. Different colored plots indicate the variation in bacterial community composition according to land use types and geographical locations.

Table 1

Soil samples collection sites with their GPS coordinates

S.N Collection site Sample code GPS coordinates Land use type
1 Gangwon-do GW1 38°08′25.6″N, 128°00′28.2″E Greenhouse soil
2 Gangwon-do GW2 38°08′16.5″N, 128°00′29.1″E Greenhouse soil
3 Gangwon-do GW3 38°15′47.2″N, 128°07′13.9″E Greenhouse soil
4 Gangwon-do GW4 38°16′15.5″N, 128°07′23.8″E Greenhouse soil
5 Gangwon-do GW5 38°17′54.3″N, 128°08′04.0″E Greenhouse soil
6 Gangwon-do GW6 38°18′13.9″N, 128°08′02.4″E Greenhouse soil
7 Gangwon-do GW7 38°12′04.9″N, 128°02′40.7″E Greenhouse soil
8 Gangwon-do GW8 38°11′50.1″N, 128°02′40.3″E Greenhouse soil
9 Gangwon-do GW9 38°11′55.3″N, 128°02′20.8″E Greenhouse soil
10 Gangwon-do GW10 38°10′49.6″N, 128°02′21.1″E Greenhouse soil
11 Gyeongsangbuk-do GB1 36°50′26.7″N, 128°54′40.4″E Open field soil
12 Gyeongsangbuk-do GB2 36°51′58.0″N, 128°52′41.6″E Open field soil
13 Gyeongsangbuk-do GB3 36°48′17.4″N, 128°57′49.1″E Greenhouse soil
14 Gyeongsangbuk-do GB4 36°49′04.5″N, 128°57′45.7″E Open field soil
15 Gyeongsangbuk-do GB5 36°54′34.2″N, 128°52′43.6″E Open field soil
16 Gyeongsangbuk-do GB6 36°53′54.0″N, 128°52′18.9″E Open field soil
17 Gyeongsangbuk-do GB7 36°53′17.7″N, 128°49′04.3″E Open field soil
18 Gyeongsangbuk-do GB8 36°52′02.1″N, 128°45′23.6″E Open field soil
19 Gyeongsangbuk-do GB9 36°49′00.4″N, 128°45′13.2″E Open field soil
20 Gyeongsangbuk-do GB10 36°52′22.8″N, 128°50′01.4″E Open field soil
21 Gyeongsangbuk-do GB11 36°50′30.5″N, 128°54′19.5″E Open field soil
22 Gyeongsangbuk-do GB12 36°49′55.9″N, 128°58′46.7″E Open field soil
23 Gyeongsangbuk-do GB13 36°53′29.4″N, 129°00′15.4″E Open field soil
24 Gyeongsangbuk-do GB14 36°53′37.5″N, 129°00′21.6″E Open field soil
25 Gyeongsangbuk-do GB15 36°56′28.1″N, 129°01′17.9″E Open field soil
26 Gyeongsangbuk-do GB16 36°51′01.2″N, 129°01′39.2″E Open field soil
27 Gyeongsangbuk-do GB17 36°49′37.2″N, 128°57′52.7″E Open field soil
28 Gyeongsangbuk-do GB18 36°49′51.3″N, 128°57′13.4″E Open field soil
29 Gyeongsangbuk-do GB19 36°49′56.5″N, 128°53′50.2″E Open field soil
30 Gyeongsangbuk-do GB20 36°50′38.6″N, 128°54′10.3″E Open field soil
31 Gyeongsangbuk-do GB21 36°44′12.9″N, 128°51′50.1″E Open field soil
32 Gyeongsangbuk-do GB22 36°44′17.0″N, 128°53′09.0″E Open field soil
33 Gyeongsangnam-do KN1 35°18′08.3″N, 128°19′06.4″E Greenhouse soil
34 Gyeongsangnam-do KN2 35°18′29.0″N, 128°18′50.2″E Greenhouse soil
35 Gyeongsangnam-do KN3 35°19′01.5″N, 128°18′52.1″E Greenhouse soil
36 Gyeongsangnam-do KN4 35°18′41.8″N, 128°17′59.0″E Greenhouse soil
37 Gyeongsangnam-do KN5 35°18′32.8″N, 128°17′52.9″E Greenhouse soil
38 Gyeongsangnam-do KN6 35°21′05.4″N, 128°20′06.8″E Greenhouse soil
39 Gyeongsangnam-do KN7 35°21′09.5″N, 128°20′05.6″E Greenhouse soil
40 Gyeongsangnam-do KN8 35°21′08.7″N, 128°20′10.9″E Greenhouse soil
41 Gyeongsangnam-do KN9 35°21′10.0″N, 128°25′08.0″E Greenhouse soil
42 Gyeongsangnam-do KN10 35°21′00.1″N, 128°24′50.2″E Greenhouse soil
43 Jeollabuk-do JB1 35°20′33.9″N, 126°35′45.4″E Open field soil
44 Jeollabuk-do JB2 35°21′08.1″N, 126°35′45.3″E Greenhouse soil
45 Jeollabuk-do JB3 35°21′25.1″N, 126°35′44.3″E Open field soil
46 Jeollabuk-do JB4 35°22′49.2″N, 126°35′46.2″E Greenhouse soil
47 Jeollanam-do JN1 34°58′28.2″N, 126°47′24.4″E Open field soil
48 Jeollanam-do JN2 34°58′29.6″N, 126°47′22.3″E Open field soil
49 Jeollanam-do JN3 35°10′52.9″N, 126°59′59.0″E Open field soil
50 Chungcheongbuk-do CB1 36°38′22.3″N, 127°28′11.4″E Greenhouse soil
51 Chungcheongbuk-do CB2 36°38′22.8″N, 127°28′11.4″E Greenhouse soil
52 Chungcheongbuk-do CB3 36°50′53.6″N, 127°32′37.6″E Greenhouse soil
53 Chungcheongbuk-do CB4 36°53′30.3″N, 127°31′58.9″E Greenhouse soil
54 Chungcheongbuk-do CB5 36°56′03.2″N, 127°33′41.3″E Greenhouse soil
55 Chungcheongbuk-do CB6 36°57′00.4″N, 127°31′14.6″E Greenhouse soil
56 Chungcheongbuk-do CB7 36°56′23.7″N, 127°30′19.2″E Greenhouse soil
57 Chungcheongbuk-do CB8 36°57′04.4″N, 127°29′03.1″E Greenhouse soil
58 Chungcheongbuk-do CB9 36°56′42.0″N, 127°26′45.6″E Greenhouse soil
59 Chungcheongbuk-do CB10 36°56′57.7″N, 127°27′57.3″E Greenhouse soil
60 Chungcheongbuk-do CB11 36°59′11.5″N, 127°28′39.9″E Greenhouse soil
61 Chungcheongbuk-do CB12 36°58′43.6″N, 127°29′25.9″E Greenhouse soil
62 Chungcheongbuk-do CB13 36°59′31.9″N, 127°28′47.2″E Greenhouse soil
63 Chungcheongbuk-do CB14 37°11′56.2″N, 127°44′26.1″E Greenhouse soil
64 Chungcheongnam-do CN1 36°19′22.7″N, 126°56′08.9″E Greenhouse soil
65 Chungcheongnam-do CN2 36°19′22.7″N, 126°56′08.9″E Greenhouse soil
66 Chungcheongnam-do CN3 36°20′29.2″N, 126°53′52.2″E Greenhouse soil
67 Chungcheongnam-do CN4 36°19′48.7″N, 126°54′11.1″E Greenhouse soil
68 Chungcheongnam-do CN5 36°19′58.8″N, 126°54′29.8″E Greenhouse soil
69 Chungcheongnam-do CN6 36°20′12.4″N, 126°53′52.3″E Greenhouse soil
70 Chungcheongnam-do CN7 36°20′25.7″N, 126°54′49.0″E Greenhouse soil
71 Chungcheongnam-do CN8 36°20′15.7″N, 126°54′56.3″E Greenhouse soil
72 Chungcheongnam-do CN9 36°19′37.2″N, 126°55′28.6″E Greenhouse soil
73 Chungcheongnam-do CN10 36°19′38.4″N, 126°55′10.6″E Greenhouse soil
74 Chungcheongnam-do CN11 36°19′52.3″N, 126°54′50.1″E Greenhouse soil
75 Chungcheongnam-do CN12 36°16′46.8″N, 127°02′39.6″E Greenhouse soil
76 Chungcheongnam-do CN13 36°16′47.9″N, 127°02′35.1″E Greenhouse soil
77 Chungcheongnam-do CN14 36°16′40.7″N, 127°02′53.0″E Greenhouse soil
78 Chungcheongnam-do CN15 36°16′50.8″N, 127°02′45.6″E Greenhouse soil
79 Chungcheongnam-do CN16 36°15′30.2″N, 126°51′16.7″E Greenhouse soil
80 Chungcheongnam-do CN17 36°14′43.2″N, 126°51′54.8″E Greenhouse soil
81 Gyeonggi-do GG1 37°32′13.0″N, 127°41′21.1″E Greenhouse soil
82 Gyeonggi-do GG2 37°31′55.7″N, 127°41′21.0″E Greenhouse soil
83 Gyeonggi-do GG3 37°31′35.0″N, 127°41′14.7″E Greenhouse soil
84 Gyeonggi-do GG4 37°31′22.3″N, 127°41′20.3″E Greenhouse soil
85 Gyeonggi-do GG5 37°31′36.5″N, 127°41′13.3″E Greenhouse soil
86 Gyeonggi-do GG6 37°32′54.4″N, 127°43′54.9″E Greenhouse soil
87 Gyeonggi-do GG7 37°32′24.8″N, 127°44′31.4″E Greenhouse soil
88 Gyeonggi-do GG8 37°32′32.2″N, 127°44′51.2″E Greenhouse soil
89 Gyeonggi-do GG9 37°32′16.6″N, 127°35′27.7″E Greenhouse soil
90 Gyeonggi-do GG10 37°32′09.2″N, 127°35′36.4″E Greenhouse soil
91 Gyeonggi-do GG11 37°32′21.5″N, 127°41′40.4″E Greenhouse soil
92 Gyeonggi-do GG12 37°32′30.5″N, 127°41′21.0″E Greenhouse soil