Introduction

미생물은 토양 입단화를 유도하고 (Tisdall and Oades, 1982) 유기물 분해를 통해 탄소 순환에 기여하며 (Miltner et al., 2012; Hemkemeyer et al., 2021) 인산 가용화 (Khan et al., 2016) 및 토양 비옥도에 관여하는 등 토양 생태계에서 다양하고 중요한 역할을 하는 것으로 알려져 있다 (Lynch and Bragg, 1985; Marschner et al., 2004; Hemkemeyer et al., 2021). 토양 뿐만 아니라 식물에도 큰 영향을 미치는데, 식물이 사용 가능한 형태로 질소를 변환시키는 등 토양 질소 순환에 중요한 역할을 한다 (Bohlool et al., 1992). 또한 식물 생장에 간접적으로 도움을 주기도 하며 식물의 병 저항성을 유도하거나 미생물 간 상호작용으로 병을 억제하기도 한다 (Compant et al., 2010; Cha et al., 2016; Dewen et al., 2017; Vlot et al., 2021). 따라서 토양 건전성을 증가시키고 지속가능한 농업을 이루기 위해 최근 토양 미생물에 대한 연구가 활발히 이루어지고 있다. 그러나 토양 미생물은 흙 1 g에 수 천 개 이상의 종이 존재하며 복잡한 군집을 이루고 있어 아직까지 알려지지 않은 부분이 많다 (Whitman et al., 1998; Torsvik et al., 2002; Fierer, 2017).

일반적으로 미생물이 풍부하고 군집의 구성이 다양할수록 토양과 작물에 긍정적인 영향을 미치는 것으로 알려져 있다 (Matos et al., 2005; Lamb et al., 2011). 여러 요인들이 토양미생물의 군집에 영향을 미치는데, 그 중 하나로 식물 뿌리 주변에 형성되는 근권 (rhizosphere)이 있으며 식물 뿌리는 다양한 생리활성 화합물을 근권으로 배출하고 이 삼출물에 이끌려 미생물들이 모여들게 된다 (Sørensen, 1997; Brimecombe et al., 2000; Bais et al., 2006). 따라서 식물의 종류에 따라 미생물 군집이 변화하며 근권은 특정 미생물의 비율이 높은 반면 비근권 토양은 균일하고 풍부한 미생물 군집 형태를 나타낸다 (Schreiter et al., 2014; Bakker et al., 2015). 그러나 근권이 아닌 서식지 유형이나 토지 이용 형태 등에 따라 미생물 분포 또한 크게 변화된다. 최근 연구에 따르면 프랑스나 뉴질랜드에서 산림, 초원 등 토양 유형에 따라 토양 특성과 미생물 간 상관관계가 있었고 미생물 풍부도와 네트워크에 차이가 있었다 (Hermans et al., 2017; Karimi et al., 2019). 담수 환경인 우리나라 벼논은 혐기 조건에서 많이 발견되는 Chloroflexota의 상대풍부도가 높았고 밭, 과수원, 시설재배지에 비해 미생물 네트워크가 더 많이 형성되는 것으로 밝혀졌다 (Petriglieri et al., 2018; Lee et al., 2020). 토양 유형 뿐만 아니라 토양 관리 방법에 따라 달라지는 토양 화학성은 미생물 군집 변화에 상당한 영향을 미친다고 잘 알려져 있다 (Aziz et al., 2013; Bartram et al., 2014; Yoon et al., 2022; Kim et al., 2023). 이전 연구들은 pH, 전기 전도도, 유기물 함량, 가용성 인산 등 변화에 따라 미생물 군집 구성이 달라지는 것을 확인하였고 (Chodak et al., 2015; Kim et al., 2016; Xu et al., 2020), 특히 pH는 다른 이화학성에 비해 미생물 군집 구조와 다양성에 큰 영향을 미친다 (Lauber et al., 2009; Rousk et al., 2010; Kim et al., 2016; Lee et al., 2020). 토양 pH는 산성비와 같은 기후 요인에 의해 감소하거나 농경지의 경우 암모늄 비료 사용에 의해 산성화되기도 하는데 (Tian and Niu, 2015; Goulding, 2016), Tripathi et al. (2018)의 연구에 따르면 토양 pH가 산성 또는 염기성일수록 하나의 환경 필터로 작용하여 특정 계통으로 분화를 유도하지만 pH가 중성일수록 미생물 군집의 분포 비율은 균등해진다 (Lauber et al., 2009; Tripathi et al., 2018). 또한 공간적 거리보다 pH 차이가 미생물 군집에 더 큰 영향을 주었다. 현재까지 토양 pH가 미생물 군집 변화에 영향을 미치며 문 (phylum) 수준에서 미생물 변화에 대한 많은 연구가 이루어졌지만 pH에 따른 속 (genus) 수준에서의 핵심 미생물에 관한 연구는 미미하다. 따라서 본 연구는 속 수준에서 전국 과수원의 중성 토양과 산성 토양을 비교하여 미생물 군집 분석을 수행하였고, 이를 통해 우리는 토양 pH에 따라 민감하게 반응하는 미생물 군집 특성을 활용하여 토양 건전성 평가 지표 개발을 위한 기반을 마련하고자 한다.

Materials and Methods

토양 시료 채취

본 연구를 위해 제주를 포함한 전국 과수원 199지점을 선정하였다. 토양 시료는 2022년 4 - 6월에 과수원 내 비근권 토양 3곳의 표토 (0 - 15 cm)를 채취하여 하나로 혼합하였으며, 온도변화에 따른 토양 시료의 변화를 최소화하기 위해 스티로폼박스에 보관하여 실험실로 운반하였다. 수집된 토양 시료는 2 mm 체로 거른 뒤 미생물 활성의 균일화를 위해 22°C에서 6 - 7일간 배양하였다. 미생물 군집분석을 위한 토양은 -80°C 냉동고에 보존하면서 사용하였다.

토양 화학성 분석

토양 화학성은 풍건 후 2 mm 체로 거른 후 농촌진흥청 분석법에 따라 분석하였다 (NAAS, 2010). 토양 pH와 EC는 1:5법 (w/w)을 이용하여 pH meter, EC meter로 측정하였고, EC값은 측정값에 5를 곱하여 구하였다. 토양 유기물 함량은 Tyurin법에 따라 0.4 N K₂Cr₂O₇ 용액으로 분해 후 0.2 N Fe(NH₄)₂(SO₄)₂ ‧ 6H₂O 용액으로 역적정하여 토양 탄소함량을 구하고 보정계수 (1.724)를 곱하였다. 유효인산은 Lancaster법으로 추출 및 발색 후 비색계로 흡광도를 측정하여 분석하였다. 치환성 양이온 함량은 1 N NH₄CH₃CO₂ (pH 7.0) 1:10 (w/v) 추출액을 inductively coupled plasma optical emission spectroscopy (ICP-OES)로 측정하였다.

DNA 추출 및 시퀀싱

토양시료로부터 DNeasyPowerSoil kit (Qiagen, Germany)를 사용하여 제조사 프로토콜에 따라 토양 미생물 DNA를 추출하였다. 추출된 DNA는 Quant-IT PicoGreen (Invitrogen, USA)을 이용하여 정량하였다.

시퀀싱 라이브러리 구축은 Illumina 16S Metagenomic Sequencing Library 프로토콜에 따라 진행하였다. 5X reaction buffer, 1 mM의 dNTP, 500 nM의 Forward/Reverse PCR 프라이머와 Heruclase II fusion DNA polymerase (Agilent Technologies, USA)를 이용하여 2 ng의 genomic DNA를 증폭하였다. 세균 분석을 위해 V3 - V4 영역의 16S rRNA을 대상으로 하는 프라이머 쌍 (V337F: 5’-T TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3’, V805R: 5’- GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVGGGTATCTAATCC-3)을, 진균 분석을 위해 ITS3 - ITS4 영역을 대상으로 하는 프라이머 쌍 (ITS3-F: 5’-CCTACGGGNGGCWGCAG-3’, ITS4-R: 5’-GACTACHVGGGTATCTAATCC-3’)을 사용하였으며 PCR 조건은 다음과 같다. 첫 변성은 95°C에서 3분 수행하였고 95°C에서 30초, 55°C에서 30초, 72°C에서 30초의 25사이클 반복 후 72°C에서 5분간 수행하였다. 이렇게 얻어진 PCR 산물은 AMPure 비드 (Agencourt Bioscience, USA)를 사용하여 정제하였고, 정제한 PCR 산물 2 µL를 인덱스가 포함된 최종 라이브러리를 구축하기 위해 NexteraXT Indexed Primer를 사용하여 앞선 PCR 조건의 25사이클 대신 10사이클 조건으로 변경하여 2차 PCR 증폭하였다. 2차 PCR 산물도 동일한 AMPure 비드로 정제하였으며, 최종 정제 산물은 Illumina Sequencing 플랫폼용 KAPA 라이브러리 정량 키트의 프로토콜에 따라 정량 하였다. 2 × 300 bp Paired-end 시퀀싱은 Macrogen Inc. (Republic of Korea)에서 MiSeq™ 플랫폼 (Illumina, USA)을 이용하였다.

MiSeq 시퀀싱 데이터 분석

토양 시료 내 미생물 분석을 위해 16S rRNA 유전자와 ITS 영역의 염기서열 판독은 operational taxonomic unit (OTU)를 기반으로 UPARSE 파이프라인을 사용하여 수행하였다 (Edgar, 2013). 두 쌍의 염기서열을 병합하고 300 bp 이하의 짧은 염기서열은 제거하였다. 이후 중복 염기서열을 제거하였으며, UPARSE-OTU 알고리즘 (https://www.drive5.com)을 통해 상동성 97%를 기준으로 클러스터링하고 키메라를 제거하였다. 위 과정은 USERACH9.1.13_i86linux64을 사용하여 수행하였다. 정제한 OTU 파일은 Mothur 1.46.1 프로그램 (Schloss et al., 2009)을 이용하여 가장 적은 염기서열 수를 기준으로 표준화 하였으며, 종 풍부도 (Chao1, ACE)와 종 균등도 그리고 종 다양성 (Shannon, Invsimpson)을 계산하였다. 또한 생물학적 분석을 위해 세균과 진균은 각각 Ribosomal Database Project (RDP, 062020 version)와 UNITE (10.05.2021 version) 데이터베이스를 사용하여 Mothur 프로그램으로 분류하였다.

미생물 군집 분석

과수원 토양의 세균과 진균의 생명정보분석을 위해 R 프로그램 (4.1.3 version)을 사용하였다. 농촌진흥청 흙토람 플랫폼에서 제공하는 과수원 토양의 적정 pH 기준을 참고하여 중성 범위 (pH 6 - 7)와 산성 범위 (pH 3 - 5)로 나누어 각각 33개, 56개를 선발하여 분석하였다. pH의 중성 범위와 산성 범위에 해당하는 과수원 토양의 미생물 군집 차이를 확인하고자 OTU를 이용한 Bray-Curtis 거리를 기반으로 vegan 패키지의 principal coordinate analysis (PCoA)를 사용하여 측정하였다. 시각화를 위해서는 ggplot2 패키지 (3.3.5 version)를 이용하였다. 두 pH 범위 간 서로 다르게 분포하는 미생물을 탐색하기 위해 analysis of composition of microbiomes with bias correction (ANCOM-BC)를 사용하여 계산하였다 (Lin and Peddada, 2020). 그 전에 계산이 용이하도록 Phyloseq 패키지를 이용하여 OTU 데이터, 생물학적 분류 데이터 그리고 시료 변수 정보 데이터를 병합 및 전처리 하였다 (McMurdie and Holmes, 2013). 결과는 ggplot2 패키지의 volcano plot으로 시각화 하였다. 두 pH 범위 간 다르게 분포하는 미생물들의 상대풍부도는 ANCOM-BC 결과를 기반으로 differential abundance (DA)를 통해 분석하였다. 군집 내 미생물 간 상관관계를 기반으로 영향력을 분석하기 위한 네트워크 분석은 로그 치환된 구성요소 간의 선형 Pearson 상관을 추정하는 sparse correlations for compositional data (SparCC)를 사용하였다. 마찬가지로 ANCOM-BC 결과에서 서로 다른 토양 pH 범위에 유의한 차이를 나타내는 미생물에 한하여 상관 관계를 추론하였다.

통계 분석

통계분석은 R 프로그램 (4.1.3 version)을 사용하여 수행하였다. 알파다양성 수치의 정규분포 여부를 확인하기 위해 샤피로 검정 (Shapiro test)을 하였고, 정규 분포일 때 등분산 검정 (F test)을 하였다. 모집단의 분산이 같을 때와 같지 않을 때 각각 Student T test와 Welch’s T test를 수행하였다. 정규 분포이지 않을 때 비모수 검정으로 Wilcoxon test를 수행하였다. PCoA의 유의성 검정을 위해 vegan 패키지의 adonis 기능을 사용하여 permutational multivariate analysis of variance (PERMANOVA)를 수행하였다.

Results and Discussion

토양의 이화학성 중에 토양 미생물 군집에 가장 큰 영향을 미치는 요인은 pH로 알려져 있다 (Fierer and Jackson, 2006; Constancias et al., 2015). 본 연구는 농촌진흥청 흙토람 플랫폼에서 제공하는 과수원 토양의 화학성 기준을 토대로 적정 pH인 중성 범위와 산성 범위의 두 개 pH범위를 선정하였다. 전국 과수원 토양 199지점 중 중성 범위 63지점과 산성 범위 26지점을 선별하여 토양 내 세균과 진균의 알파다양성을 비교하였다 (Fig. 1). 그 결과, 중성 범위의 세균 종 풍부도 (Chao1 지수, Welch T test, P < 0.001***)와 종 다양성 (Shannon 지수,Wilcoxon test, P < 0.001***) 평균치는 각각 4,583과 6.7로 산성 범위의 3,079와 6.1에 비해 유의성 있게 높았다. 진균의 경우는 중성 범위의 종 풍부도 평균치는 879로 산성 범위 701에 비해 유의성 있게 높은 반면 Shannon 지수 (Student T test, P = 0.3)는 유의한 차이가 없었다. 이는 선행 연구에서 토양 pH가 중성에서 세균의 종 풍부도와 종 다양성이 산성일 때보다 높다라는 결과와 유사하다 (Fierer and Jackson, 2006; Lee et al., 2020). 세균의 경우 대부분이 세포 내 pH가 중성에 가깝기 때문에 토양의 산성 pH는 특정 세균을 제외한 대부분의 세균에 상당한 스트레스를 줄 수 있다 (Lauber et al., 2009). 따라서 산성 pH 범위에서 미생물의 성장과 활성이 감소하여 종 풍부도와 종 다양성이 낮았던 것으로 판단된다. 작물 또한 중성 토양에서 생육이 활발하기 때문에 (Läuchli and Grattan, 2012), 중성 pH에서 미생물 특성이 작물과 직접적인 연관성이 있는 지 근권 토양에 대한 추가적인 연구를 통해 확인할 필요가 있다. 그러나 비근권 토양이더라도 토양 내 미생물이 풍부하고 다양할수록 토양 탄소 및 질소의 순환이 원활하고 토양 병이 억제되는 등 긍정적인 간접 효과가 있기 때문에 pH에 따른 미생물 군집의 알파 다양성 변화는 충분히 토양 건전성을 평가할 수 있을 것으로 판단된다.

Fig. 1.

Comparison of microbial alpha diversity in orchard soils according to pH range. Box plots represent richness and diversity of microorganisms in the range of soil pH 3 - 5 (acidic Soil) and 6 - 7 (neutral Soil). Chao1 and Shannon index indicate microbial richness and diversity, respectively. (A) Chao1, (B) Shannon of bacteria and (C) Chao1 of fungi were significantly higher in the neutral soil (P < 0.001***), but (D) Shannon of fungi was not (P = 0.3).

PCoA 분석

토양 내 미생물 분포를 비교하고 차이를 파악하기 위해 PCoA를 사용하여 세균과 진균의 과 (family) 수준에서 두 pH 그룹 간 군집 조성과 거리를 측정하였다 (Fig. 2). 두 그룹은 세균의 경우 첫 번째 요인 (PCo1: 15.2%)에 의해 분리되었고, 진균의 경우 마찬가지로 첫 번째 요인 (PCo1: 4%)에 의해 분리되었다. 이 결과는 adonis 기능을 이용한 PERMANOVA (P = 0.001***)에서 통계적으로 유의한 차이가 있는 것으로 나타났다. 따라서 토양 pH에 따라 세균과 진균 군집이 달라지는 것을 확인하였고 이전 연구 결과와 유사한 경향을 보였다 (Fierer and Jackson, 2006). 본 연구는 수종의 영향이 매우 적을 것으로 판단되나 조사 대상 지점 중 귤을 재배하는 제주 토양의 경우, 제주 화산회 토양이 내륙 토양과 달라 미생물 군집 구조에 영향을 주었을 것으로 추정된다.

Fig. 2.

PCoA for the analysis of differences in microbial communities according to pH range. The red and blue nodes represent microbial communities at the family level in the acidic and neutral soil, respectively. Both (A) bacteria and (B) fungi showed a significant difference between the soil pH range as a result of the PERMANOVA test (P = 0.001).

ANCOM-BC

우리는 PCoA를 통해 과수원 토양 pH에 따른 두 그룹의 미생물 군집이 유의한 차이가 있는 것을 확인하였고, 두 그룹 토양 내 서로 다르게 분포하는 미생물들을 파악하기 위해 속 (genus) 수준에서 ANCOM-BC를 수행하였다 (Fig. 3). 그래프의 x축은 Log Fold Change (LFC) 값으로, 판별 효과를 나타내며 y축은 유의성을 나타내는 q값 (adjust p value)이다. 붉은 선은 결과 해석을 용이하게 하기 위한 한계 값이며, LFC 값은 세균의 경우 ±2, 진균의 경우 ±1로 설정하였고 q값은 세균과 진균 모두 0.01로 설정하였다. LFC 값이 클수록 즉, 양쪽으로 멀리 떨어질수록 두 그룹 간 다르게 분포하는 미생물로 해석할 수 있다. 그 결과, 세균은 Chujaibacter, Streptococcus, Ktedonobacter 등 8개 속이 중성 범위에 비해 산성 범위에서 많이 분포하였고, Acidobacteria에 속하는 Gp4_unclassified와 Gp11_unclassified 그리고 Povalibacter와 Chryseolinea 등 16개 속이 산성 범위에 비해 중성 범위에서 많이 분포하였다. 진균은 Saitozyma, Umbelopsis, Penicillium을 비롯한 10개 속이 중성 범위에 비해 산성 범위에서 많이 분포하였고, Gibellulopsis, Titaea, Mrakia을 비롯한 15개 속이 산성 범위에 비해 중성 범위에서 많이 분포하였다.

Fig. 3.

Volcano plot to identify differences in microorganisms of soil pH range at the genus level. The plots based on ANCOM-BC analysis of (A) bacteria and (B) fungi in orchard soils. The red line represents the difference in expression of microorganisms in the pH range and regulation according to significance (q < 0.001). Red and blue nodes suggest genus of microorganism significantly abundant in the acidic and neutral soil, respectively, by the red line.

ANCOM-BC를 통해 도출된 미생물 속의 해당 pH 범위 토양 내 분포 비율을 확인하고자 DA 분석을 수행하였다 (Fig. 4). 이를 통해 실제 토양에선 분포 비율이 매우 낮더라도 ANCOM-BC 결과에서 높은 LFC와 유의성을 나타내어 큰 영향을 미치는 주요 미생물로 잘못 해석할 수 있는 오류를 방지하였다. 그 결과, 토양 내 세균 및 진균의 상대풍부도가 0.5% 이상인 속은 각각 Povalibacter, Chujaibacter, Saitozyma, Lophotrichus 등이고, 이들은 중성 pH 범위와 산성 범위 간 최소 약 2배에서 최대 10배 이상 차이가 났다. 흥미롭게도 세균 속인 Streptococcus와 진균 속인 Umbelopsis는 중성 pH 범위 대비 산성 범위에서 각각 34배, 15배 이상 많이 존재하였다. 이 결과는 Streptococcus 일부 종의 풍부도가 토양 pH와 음의 상관관계를 나타낸다는 이전 연구와 유사하였다 (Zhao et al., 2015). 또한 Streptococcus는 대장균과 함께 토양 및 수질오염의 지표 역할을 하는데 (Van Donsel et al., 1967; Kibbey et al., 1978; Massa et al., 2001; Zhao et al., 2015; Al-Taai, 2021), 본 연구의 산성 토양 pH 범위에서 많이 분포하는 결과의 원인은 해당 지점의 사용 비료 종류 및 시용량의 영향이 있을 것으로 추정된다. 하지만 pH 외의 많은 환경 변수가 영향을 미치기 때문에 자세한 연구를 통한 원인 분석이 필요하다.

Fig. 4.

Differential abundance (DA) analysis based on ANCOM-BC. It shows the distribution of microbial genera with significant differences in each soil pH range. Red and blue bars mean the distribution ratio (%) of microorganisms, which are (A) bacteria and (B) fungi, in the acidic and neutral soil, respectively.

Network analysis

미생물 keystone 분류군은 시공간적으로 항상 존재하면서 미생물 군집 구조와 기능에 상당한 영향을 미치는 분류군으로 정의된다 (Banerjee et al., 2018; Mamet et al., 2019). 따라서 특정 미생물이 두 pH 범위 간 유의한 차이로 존재하며 상대풍부도가 높고 미생물 간 상관관계가 높으면 토양 pH에 따른 주요 미생물일 것이라는 가정 하에 SparCC 네트워크 분석을 하였다 (Fig. 5). 초록 선은 정의 상관관계를, 두께는 상관관계 값을 나타낸다 (|Correlation| > 0.7, P value < 0.001). 반면 주황 선은 부의 상관관계를 의미하지만 존재하지 않았다. 파란 원형 노드는 ANCOM-BC 기반 중성 pH 범위에 유의적으로 분포하는 분류군이며 빨간 원형 노드는 산성 pH 범위에 분포하는 분류군이다. Eigenvector centrality는 노드의 전이적 영향을 측정하는 알고리즘으로, 네트워크 연결을 점수로 할당하여 해당 노드의 영향을 측정하는 개념이다 (Bonacich, 2007; Lee et al., 2022). 노드의 크기는 Eigenvector centrality의 수치 (0.00 - 1.00)를 반영하는데, 클수록 더 많은 네트워크를 형성하는 분류군 임을 의미하며 빨간 테두리로 표시하였다 (Eigenvector centrality > 0.6). 그 결과, 세균의 경우 146개 속의 네트워크를 확인하였다. 중성 범위와 산성 범위에 분포하는 속은 각각 16개, 7개였으며 주요 속 (keystone taxa)은 15개였다. 진균의 경우 62개 속의 네트워크를 확인하였다. 중성 범위와 산성 범위에 분포하는 속은 각각 12개, 3개였으며 주요 속은 6개였다. 이 중에서 중성 범위 토양의 주요 세균 속은 9개 (Terrimonas, Povalibacter, Gp11_unclassified, Gp22_unclassified, Chloroflexia_unclassified, Sterolibacteriaceae_unclassified, Hyphomicrobiaceae_unclassified, GP17_unclassified, Gp4_unclassified)가 확인되었지만 산성 범위 토양의 주요 세균 속은 확인되지 않았다. 또한 진균은 중성 및 산성 범위 모두 주요 속은 확인되지 않았다. 미생물 간 많은 네트워크를 형성하는 미생물은 핵심 분류군으로서 많은 영향을 미칠 것이라는 이전 연구에 의하면 (Layeghifard et al., 2017; Mamet et al., 2019), Terrimonas와 Povalibacter를 비롯한 9개의 세균 속은 과수원 토양에서 적절한 pH 범위의 주요 미생물로 판단된다. 비록 분류는 되지 않았지만 Acidobacteria 문에 속하는 Gp 분류군 중 Gp1 - 3은 pH 3.5 - 5.0인 산림 토양에서 풍부한 반면 Gp4 - 7 및 Gp10은 pH 4.5 - 9.5에서 우세하다는 Kim et al. (2021)의 연구 결과와 유사한 경향이었다. 특히 Bacteroidota에 속하는 Terrimonas는 산림 토양에 풍부하여 다른 분류군과 많이 상호 작용하는 것으로 알려져 있는그람 음성, 호기성 균이다 (Fierer et al., 2007). Lu et al. (2023)에 따르면 높은 토양 유기 탄소 농도는 Terrimonas를 풍부하게 하는 경향이 있었고, Karimi et al. (2019)은 Terrimonas를 비교적 유기물 함량이 높고 미생물 네트워크 형성이 많은 산림과 초원 토양의 핵심 미생물로 보고한 바 있다. Gammaproteobacteria에 속하는 Povalibacter는 Siles and Margesin (2017)에 의하면 토양 유기물 (soil organic matter, SOM)과 토양 유기 탄소가 통계적으로 유의하게 높은 토양에서 풍부하였다. 이전 연구와 더불어 본 연구 결과에 따르면 Terrimonas와 Povalibacter가 이화학적 측면에서 안정적인 토양의 미생물 지표로 활용될 수 있음을 시사한다. 게다가 산성 토양일수록 칼슘, 아연 등 작물이 이용하는 이온 수가 적어지며 작물의 원소 불균형을 초래하는 점을 고려하였을 때 (Sims, 1986; Zhao et al., 2011), 중성 토양의 주요 미생물 9개는 과수원 토양의 좋은 세균 지표가 될 것으로 판단된다. 앞으로 데이터베이스 최신화에 따라 미분류군의 동정을 통한 자세한 분류가 필요하며 미생물의 역할과 기능에 관한 추가적인 연구가 필요하다.

Fig. 5.

Network analysis based on ANCOM-BC and DA analysis to explore key microorganism genera in acidic and neutral soils. In the community network of (A) bacteria and (B) fungi, red and blue nodes represent of microorganisms distributed in the acidic and neutral soil, respectively. The node size appears larger as the eigenvector centrality is higher, and 0.7 or more is marked with a red border.

Conclusions

본 연구는 토양 미생물과 토양 pH 간의 관계를 분석하기 위하여 전국 과수원의 중성 pH 범위 토양 56개와 산성 pH 범위 토양 23개의 미생물 군집을 비교 분석하였다. 중성 범위의 세균 풍부도와 다양성, 진균 풍부도는 산성 범위에 비해 통계적으로 유의하게 높았다 (P < 0.001). PCoA와 PERMANOVA에 따르면 두 범위 간 미생물 군집 조성과 거리는 크게 달랐으며 (P < 0.001), ANCOM-BC는 중성 토양에서 16개 세균과 15개 곰팡이 속의 높은 풍부도를 나타냈다. 마지막으로 Terrimonas와 Povalibacter를 포함한 9개 속이 중성 범위에서 다른 미생물 분류군에 중요한 영향을 미치는 핵심 미생물로 확인되었다. 앞으로 토양 pH 외 물리화학성 또는 작물 생육 등의 요인을 기준으로 적합 또는 부적합 토양을 나누어 토양 미생물 군집을 분석하고, 단일 요인 또는 다중 요인 별 핵심 미생물을 탐색한다면 보다 정확한 토양질을 평가하는 미생물 지표를 개발할 수 있을 것으로 판단된다.