📄 농업

Genetic architecture and genomic prediction of stem forking across two generations of loblolly pine (Pinus taeda L.)

이 연구는 로블로이 소나무의 줄기 갈라짐(stem forking)이 가계 수준에서는 상당한 가산적 유전 통제 하에 있어 결함을 줄이기 위한 가계 기반 게놈 선발이 효과적인 반면, 개별 나무에서의 발현은 환경 요인의 강한 영향을 받아 개체 유전율이 낮고 예측 정확도를 확보하기 어렵다는 점을 밝히고 있다.

원저자: Pinar Guner, M. Nasir Shalizi, Trevor D. Walker, Fikret Isik

게시일 2026-10-23

📖 3 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Pinar Guner, M. Nasir Shalizi, Trevor D. Walker, Fikret Isik

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

미국 남부의 숲에서 로블로이 소나무(loblolly pine)는 목재 산업의 일꾼입니다. 이 나무들은 빠르게 자라고 다양한 토양에 적응하며, 건설용 목재부터 종이 제품에 이르기까지 모든 것을 위한 목재를 제공합니다. 수십 년 동안 임업인들은 이 나무들을 더 높고 곧게 자라게 하기 위해 개량해 왔으며, 그 결과 야생림보다 훨씬 더 생산적인 거대한 조림지를 만들어냈습니다. 하지만 자연은 좀처럼 완벽하지 않으며, 가장 잘 개량된 나무조차 때때로 특정한 결함인 '줄기 갈라짐(stem fork)' 현상을 보이기도 합니다. 이는 주 줄기가 두 개 이상의 경쟁하는 주 줄기로 갈라지거나, 큰 가지가 날카롭고 어색한 각도로 자랄 때 발생합니다. 이러한 변형은 나무에 커다란 옹이와 불규칙한 결을 만들어내어, 고부가가치 제품을 만들기에는 부적합하게 만들고 폭풍 시 부러질 가능성을 높입니다. 갑작스러운 추위나 물리적 손상이 나무를 갈라지게 할 수 있다는 사실은 잘 알려져 있지만, 과학자들은 이 결함이 유전자에 얼마나 기록되어 있는지, 그리고 얼마나 단순히 날씨에 의한 불운인지에 대해 오랫동안 의문을 품어왔습니다. 이러한 균형을 이해하는 것은 다음 세대의 숲을 위한 최적의 부모 나무를 선택하고자 하는 육종가들에게 매우 중요합니다.

노스캐롤라이나 주립대학교의 연구팀은 두 세대에 걸친 대규모 로블로이 소나무 육종 프로그램을 조사하여 이 수수께 수를 풀기 위해 나섰습니다. 그들은 약 19,000그루의 나무 데이터를 수집하여 각 나무가 갈라진 줄기를 가졌는지 아니면 곧은 줄기를 가졌는지를 기록했으며, 이를 상세한 DNA 마커로부터 얻은 유전 정보와 결합했습니다. 연구진은 이 형질이 이진적(binary)이라는 통계적 과제에 직면했습니다. 즉, 나무는 줄기가 갈라져 있거나 그렇지 않거나 둘 중 하나라는 것입니다. 이를 처리하기 위해 그들은 환경의 영향과 유전의 영향을 분리할 수 있는 고급 컴퓨터 모델을 사용했습니다. 그들은 결과의 견고함을 보장하기 위해 수년간 사용되어 온 방식과 이 유형의 데이터에 더 정확하도록 설계된 새로운 접근 방식을 비교하며 다양한 계산 방식을 테스트했습니다.

연구 결과, 개별 나무와 가족 집단 사이의 극명한 대조가 드러났습니다. 단일 나무를 관찰했을 때, 연구진은 유전적 구성이 줄기 갈라짐 여부에 미치는 영향이 매우 작다는 것을 발견했습니다. 결과는 주로 서리나 토양 조건과 같은 환경적 요인에 의해 결정되었습니다. 개별 나무에 대한 유전적 영향은 너무 낮아서 특정 묘목이 갈라질지 예측하는 것이 어려웠습니다. 그러나 전체 나무 가족의 평균적인 성능을 살펴보았을 때, 양상은 완전히 바뀌었습니다. 가족 단위에서는 매우 강력한 유전적 패턴이 나타났는데, 이는 어떤 가족은 유전적으로 갈라진 나무를 생산하기 쉬운 반면, 다른 가족은 자연적으로 저항력이 있다는 것을 의미합니다. 이는 단일 나무가 갈라질지 여부를 쉽게 예측할 수는 없지만, 어떤 가족이 평균적으로 갈라진 나무를 적게 생산할지는 신뢰할 수 있게 예측할 수 있음을 시사합니다.

연구팀은 또한 줄기 결함의 책임이 있는 특정 유전자를 찾기 위해 나무의 DNA를 탐색하며, 쉽게 추적할 수 있는 '마스터 스위치'를 찾기를 희망했습니다. 그들은 특정 가족 혈통 내에서만 줄기 갈라짐과 강하게 연결된 하나의 특정 유전적 지점을 식별했습니다. 해당 특정 가족의 경우, 특정 버전의 유전 마커가 나무가 갈라짐 현상을 보일 확률을 높였습니다. 그러나 연구진이 이 마커를 사용하여 전체 육종 집단의 예측력을 개선하려고 시도했을 때, 이는 도움이 되지 않았습니다. 그 마커는 전체 집단에 활용되기에는 특정 가족 계보에 너무 국한되어 있었습니다. 이 발견은 줄기 갈라짐이 모든 나무에 영향을 미치는 하나의 주요 유전자에 의해 발생하는 것이 아니라, 환경 조건에 의해 쉽게 압도될 수 있는 수많은 작은 유전적 효과들이 함께 작용하여 발생하는 것임을 시사합니다.

산림 경영에 주는 시사점은 명확합니다. 줄기 갈라짐 형질은 개별 수준에서 환경의 영향을 매우 많이 받기 때문에, 나무 자체의 외형만을 기준으로 최고의 단일 나무를 선택하는 것은 그리 효과적이지 않습니다. 대신, 가장 성공적인 전략은 전체 가족의 평균 성능을 바탕으로 가족을 선택하는 것입니다. 유전체 도구를 사용하여 가족의 유전적 잠재력을 예측함으로써, 육종가들은 다음 세대의 개별 나무들이 여전히 예측 불가능한 날씨에 직면하더라도, 갈라지는 경향을 물려줄 가능성이 낮은 부모를 선택할 수 있습니다. 연구는 또한 성장과 결함 사이의 미묘한 연관성도 찾아냈습니다. 더 빨리 자라는 나무일수록 갈라짐의 유전적 위험이 약간 더 높은 경향이 있었습니다. 이는 육종가들이 단순히 가장 빨리 자라는 나무만을 쫓아서는 안 되며, 미래의 숲이 생산성과 품질을 모두 갖출 수 있도록 성장과 형태 사이의 균형을 맞춰야 함을 의미합니다. 궁극적으로, 이 연구는 우리가 날씨를 통제할 수는 없지만, 유전학을 사용하여 날씨에 더 잘 견딜 수 있는 숲을 구축할 수 있음을 확인시켜 줍니다.

기술 요약: 로블로리 소나무(Loblolly Pine)의 줄기 분지(Stem Forking)에 대한 유전 구조 및 게놈 예측

문제 제기
주줄기가 두 개의 우점 줄기로 나뉘거나 큰 가지가 나타나는 현상인 줄기 분지는 Pinus taeda(로블로리 소나무)에서 목재 품질을 저하시키고 경제적 가치를 감소시키는 중요한 결함이다. 서리나 기계적 손상과 같은 환경적 요인은 잘 알려져 있으나, 이 이진 형질(binary trait)의 유전 구조는 여전히 명확히 밝혀지지 않았다. 기존 연구들은 적은 표본 크기, 정확한 분산 성분 추정을 방해하는 실험 설계, 그리고 분지가 다유전자 효과(polygenic effects)에 의해 조절되는지 아니면 주요 효과 유전자좌(major-effect loci)에 의해 조절되는지에 대한 불명확성으로 인해 한계를 보여왔다. 또한, 대규모 육종 집단 내 이진 형질 분석 시 발생하는 통계적 과제, 특히 동물 모델(animal models)에서 벌칙 유사-우도법(Penalized Quasi-Likelihood, PQL)이 가질 수 있는 편향성 문제도 이 맥락에서 충분히 다뤄지지 않았다.

연구 방법론
본 연구는 미국 남부의 여러 현장에서 평가된 19,165그루(1세대 15,289그루, 2세대 3,876그루)의 2세대 육종 집단을 활용하였다. 표현형 데이터는 6년생(G1) 및 4년생(G2)에서 이진 형질(분지 있음 vs. 없음)로 기록되었다. 유전형 분석은 Pita50K Axiom™ SNP 어레이를 사용하여 수행되었으며, 품질 관리를 거쳐 최종적으로 15,848개의 SNP 데이터셋을 확보하였다. 그러나 하위 게놈 분석은 4,003그루(G1 2,067그루, G2 1,337그루)의 유전형 분석 대상 개체만을 대상으로 진행되었다.

분석 프레데임워크는 다음과 같은 고급 통계 및 게놈 접근법을 통합하였다:

유전 파라미터 추정: 저자들은 조건부 PQL(logit 및 probit 연결 함수 사용)을 이용한 일반화 선형 혼합 모델(GLMM)을 Echidna 소프트웨어로 구현된 한계 PQL(marginal PQL, mPQL) 방식과 비교하였다. 잠재적 책임성 척도(latent liability scale)에서의 유전력을 추정하기 위해 개체 단위(animal) 모델과 가계 단위(일반 결합 능력, GCA) 모델을 모두 적합시켰다.
게놈 예측: 단일 단계 게론믹 BLUP(ssGBLUP)을 사용하여 혈통 및 게놈 관계 행렬( $\mathbf{H}$ 행렬)을 혼합하여 육종가(breeding values)를 예측하였다. 교차 세대 예측 능력은 G1에서 모델을 학습시키고 G2에서 검증하는 방식으로 평가되었다.
QTL 및 연관 매핑: G1 집단에서 MLMM, FarmCPU, BLINK 방법을 사용하여 전전체 연관 분석(GWAS)을 수행하였다. 또한, 특정 부모를 공유하는 6개의 완전 형제 가족(full-sib families)을 대상으로 가족 기반 연관 매핑을 수행하여 양적 형질 유전자좌(QTL)를 식별하였다.
마커 통합: GWAS 및 연관 매핑을 통해 식별된 유의미한 마커를 ssGBLUP 모델의 고정 공변량으로 테스트하여 예측 성능에 미치는 영향을 평가하였다.
다변량 분석: 다변량 ssGBLUP 모델을 사용하여 분지, 수고(tree height), 흉고 직경(DBH) 간의 가산적 유전 상관관계를 추정하였다.

주요 결과

유전 구조: 개체 단위의 협의 유전력( $h^2$ )은 0.09–0.13으로 낮게 나타났는데, 이는 개별 나무의 분지 발현이 환경 및 발달 요인의 영향을 크게 받는다는 것을 의미한다. 반면, 가계 평균 유전력은 0.93으로 높게 나타나, 가계 수준에서는 유의미한 가산적 유전 제어가 존재함을 확인하였다.
모델 성능: probit PQL 모델이 logit 및 marginal 모델보다 약간 높은 유전력 추정치를 보였다. 그러나 모든 모델 구성에서 예측된 육종가는 거의 완벽하게 상관되어 있었으며, 이는 모델 선택이 나무나 가계의 순위 결정에 미치는 영향이 미미함을 시사한다.
QTL 식별: GWAS와 연관 매핑을 통해 연관 그룹 3(linkage group 3)에서 마커 SNP505839와 관련된 유의미한 계통 특이적 연관성을 발견하였다. 이 마커는 부모 N11210의 자손들 사이에서만 분지와 유의미한 연관성을 보였으며, C 대립유전자의 복제수가 하나씩 증가할 때마다 분지 발생 오즈(odds)가 45% 증가하였다(오즈비 = 1.45).
게놈 예측: 세대 간 가계 수준 예측 능력은 중간 수준( $r = 0.40$ )이었다. G2의 개체 단위 예측은 약했다(AUC = 0.55). 결정적으로, 유의미한 마커 SNP505839를 고정 공변량으로 포함시켜도 표준 ssGBLUP 모델과 비교했을 때 가계 수준의 예측 능력이나 개체 단위의 변별력(AUC)이 개선되지 않았다.
유전 상관관계: 분지는 수고( $r_a = 0.36$ ) 및 흉고 직경( $r_a = 0.24$ )과 완만한 양의 가산적 유전 상관관계를 보였으며, 이는 성장만을 위한 선발이 의도치 않게 분지 위험을 높일 수 있음을 시사한다.

의의 및 주장
본 논문은 로블로리 소나무의 줄기 분지 기저 유전은 다유전자적(polygenic)인 반면, 그 표현형 발현은 환경 요인에 의해 강력하게 조절된다고 주장한다. 본 연구는 개별 나무의 표현형이 유전적 가치를 판단하는 지표로는 부적합하지만, 가계 평균 평가는 매우 효과적임을 입증하였다. 따라서 저자들은 게놈 선발을 가계 수준에서 적용하는 것이 육종 집단 내 줄기 분지 발생을 줄이는 데 효과적인 전략이 될 것이라고 결론짓는다.

방법론 측면에서, 본 연구는 이진 형질에 대한 marginal PQL의 비교 평가를 통해, 해당 방식이 대안적인 근사치를 제공하기는 하지만, 발생 빈도가 높고 표본 크기가 큰 본 연구의 맥락에서는 유전력 추정치나 순위에 유의미한 변화를 주지 않았음을 밝혀냈다. 또한, 특정 혈통(N11210) 내에서 탐지된 주요 효과 QTL이 반드시 전체 집단 수준의 예측력 향상으로 이어지지는 않는다는 점을 강조하며, 이는 특정 고정 효과 마커에 의존하지 않고도 다유전자 효과를 포착하는 데 있어 밀집된 마커 커버리지를 가진 ssGBLUP의 유용성을 뒷받침한다. 본 연구의 결과는 성장 및 줄기 형태와 함께 분지를 다형질 선발 체계에 통합함으로써, 상관된 부정적 반응의 위험을 완화할 수 있다는 근거를 제공한다.

유사한 논문