Diversity, Nesting Ecology and Worker-level Morphometric Structure of Stingless Bees (Hymenoptera: Apidae: Meliponini) in the Bimodal Forest Zone of Central Cameroon

기술 요약: 중앙 카메룬의 무가슴벌(Stingless Bees)의 다양성, 둥지 생태학 및 일벌 수준의 형태측정 구조

문제 제기
아프리카 열대 지역의 무가슴벌(Meliponini 부족)은 생태적, 경제적으로 중요한 수분 매개자임에도 불구하고, 네오트로피카(Neotropical)나 아시아 지역에 비해 연구 기록이 매우 미비하다. 카메룬의 경우, 기존 연구가 고산 지대에 집중되어 있어, 센터(Centre) 지역의 이봉형 산림 지대에 대한 조사는 거의 이루어지지 않았다. 또한, Hypotrigona 및 Liotrigona와 같은 주요 속(genus)의 분류 체계는 미해결 상태이며, 잠재적인 은폐종(cryptic diversity)의 존재 가능성이 있다. 전통적인 형태측정법은 벌의 크기와 형태를 기술하는 실용적인 도구이지만, 아프리카 종을 위한 지역별 형태측정 기준선(baseline)이 부족하며, 환경이나 알로메트리(allometry, 상대 성장)에 의한 표현형 가소성과 실제 유전적 분화를 구분할 필요가 있다. 본 연구는 이 특정 농업 생태 구역 내 무가슴벌의 생물 다양성 데이터, 둥지 생태학 및 형태적 특성 규명에 대한 공백을 메우고자 한다.

연구 방법
본 연구는 2021년 12월부터 2022년 6월까지 카메룬 센터 지역의 5개 행정 구역(Méfou-et-Akono, Lékié, Mfoundi, Méfou-et-Afamba, Nyong-et-So'o)에서 수행되었다.

• 샘플링 및 인벤토리: 연구진은 현지 양봉가 및 주민들의 도움을 받아 107개의 군집을 찾아냈다. 둥지 기질을 기록하였으며, 형태측정 분석을 위해 일벌들을 채집하였다.
• 형태측정: 24가지 일벌 형질(몸, 머리, 날개, 다리 치수 및 앞날개 감각모(hamuli) 개수 포함)을 스테레오 현미경과 ZEN 소프트웨어를 사용하여 3회 반복 측정하였다. 분석은 가장 흔한 두 종인 Hypotrigona gribodoi(n=229)와 Dactylurina staudingeri(n=201)에 집중되었다.
• 통계 분석:
  • 군집 구조: 발견된 군집 및 일벌 수를 기준으로 샤논(Shannon), 심슨(Simpson), 피엘루(Pielou) 지수를 계산하였다.
  • 단변량/다변량 검정: ANOVA, Kruskal–Wallis, MANOVA를 통해 지역 효과를 평가하였다.
  • 형태형(Morphotype) 구분: 차원 축소를 위해 주성분 분석(PCA)을 사용하였다. 응집 계층적 클러스터링(Ward's linkage)을 통해 후보 그룹을 제안하였고, 이후 leave-one-out 교차 검증을 포함한 판별 함수 분석(DFA)을 통해 이를 검증하였다.
  • 알로메트리 및 지리학: 크기와 형태를 구분하기 위해 로그-형태 변환(log-shape transformation)을 사용하여 크기를 보정하였다. 맨텔 검정(Mantel tests)을 통해 형태측정적 거리(Mahalanobis)와 지리적 거리 간의 상관관을 평가하였다.
  • 검증: 식별된 클러스터가 연속적인 경사(gradient)나 무작위 분할의 결과물이 아님을 보장하기 위해 순열 검정(permutation tests), 부트스트랩 재표본 추출, 갭 통계량(gap statistics)을 사용하였다.

주요 결과

• 다양성 및 둥지: 4개의 분류군이 기록되었다: H. gribodoi, D. staudingeri, Meliponula ferruginea, 그리고 미확인 Hypotrigona-유사 형태종. H. gribodoi가 압도적으로 우세하였으며(발견된 군집의 93.5%), 주로 건물 틈새(벽, 창문)에 둥지를 틀었다. D. staudingeri(5개 군집)는 노출된 수목 둥지를 구축하였다. M. ferruginea와 미확인 형태종은 드물게 나타났으며, 단일 지역에서만 발견되었다.
• 군집 지수: H. gribodoi의 우점으로 인해 군집 기반 다양성은 낮았다(Shannon H′ = 0.28). 일벌 기반 다양성은 더 높게 나타났다(H′ = 1.20), 이는 은폐형 둥지 조성자에 대한 탐지 편향을 강조한다.
• 형태측정 프로필: D. stestingeri는 H. gribodoi보다 선형 크기가 약 두 배 컸다. 앞날개 감각모(hamuli)의 수는 H. gribodoi의 경우 5.00으로 불변이었고, D. staudingeri의 경우 6.61 ± 0.49로 거의 일정하여 진단적 형질 역할을 하였다.
• 형태형 구조:
  • H. gribodoi: PCA와 클러스터링을 통해 두 가지 형태형이 식별되었다. 교차 검증된 DFA 정확도는 99.1%였다. 그러나 형태형 내 변이가 지배적이었으며(78.9%), 이 그룹들은 별개의 개체가 아닌 크기 구조적 연속체(continuum)를 형성하였다.
  • D. staudingeri: 네 가지 형태형이 식별되었으며, 교차 검증 정확도는 97.5%였다. 이들은 크기가 점진적으로 증가하는 경사(gradient)를 형성하였다. 형태형 내 변이가 여전히 지배적이었으나(64.2%), 형태형 간 변이 성분이 H. gribodoi보다는 컸다.
  • 알로메트리: 제1주성분(PC1)은 전체 몸 크기와 강하게 상관되었다(r ≈ 0.95). 결정적으로, 등척성(isometric) 크기 효과를 제거한 후에도 분류 정확도가 높게 유지되어( H. gribodoi 99.1%, D. staudingeri 90.0%), 크기 변화와 함께 실제적인 형태 차이가 존재함을 나타냈다.
• 지리학적 구조: 맨텔 검정 결과, 두 종 모두에서 형태측정적 분화와 지리적 거리 사이의 유의미한 상관관계가 나타나지 않았다(p > 0.28). 이는 관찰된 형태형들이 이 데이터셋을 기준으로 지리적으로 구조화된 계통이 아님을 시사한다.

주요 기여

1. 최초의 지역별 기준선 제공: 본 연구는 중앙 카메론의 이봉형 산림 지대 내 무가슴벌에 대한 최초의 포괄적인 형태측정 기준선을 제공하며, 두 주요 종에 대해 24가지 형질을 기록하였다.
2. 진단 형질 검증: 앞날개 감각모(hamulus)의 수가 H. gribodoi(5)와 D. staudingeri(~6.6)를 구분하는 견고하고 잠정적인 진단 형질임을 확인하였다.
3. 표현형 대 유전적 구분: 본 연구는 통계적으로 분리 가능한 형태형이 존재하지만, 이들이 확립된 유전적 계통이라기보다는 크기와 국지적 조건에 의해 구조화된 정적인 표현형 다양성을 나타낸다는 것을 엄격하게 입증하였다. 연구는 이 그룹들을 분자적 확인이 필요한 "운영 형태형(operational morphotypes)"으로 명시하였다.
4. 생태적 통찰: H. gribodoi의 인간 거주 구조물 내 둥지 조성 특성(친인간적 성향)과 D. staudingeri의 수목 의존성을 기록함으로써, 구체적인 생태적 맥락을 제공하였다.

의의 및 주장
저자들은 본 연구가 카메론 무가슴벌의 보전과 활용을 위한 기초적인 참조 자료를 구축한다고 주장한다. "지역별 형태측정 기준선"을 제공함으로써, 본 연구는 종의 신뢰할 수 있는 식별을 가능하게 하며, 이는 멜리포니컬처(meliponiculture, 무가슴벌 양봉)를 위한 모든 육종, 종 선택 또는 품질 보증 프로그램의 전제 조건이다.

본 논문은 식별된 형태형이 재현 가능한 통계적 그룹이지만, 분자 데이터(예: COI 바코딩) 없이 이를 별개의 생물학적 종으로 해석하는 데에는 주의가 필요하다고 겸허히 주장한다. 또한, H. gribodoi는 그 둥지 습성을 통해 저투입 가축화 가능성을, D. staudingeri는 그 방어적 특성을 통해 가치를 제시하는 동시에, 미확인 Hypotrigona-유사 형태종이 분자적 확인이 필요한 후보 분류군임을 언급하였다. 궁극적으로, 본 연구는 아프리카 무가슴벌의 보전 및 수분 전략이 각 종의 독특한 둥지 생태와 인간과의 상호작용을 고려한 종 특이적이어야 한다는 논거를 뒷받침한다.