Exploring genetic polymorphism in exon 1 region of PGLYRP-1 gene and its association with age at first calving and production traits in Sahiwal cattle

기술 요약: Sahiwal 소의 PGLYRP-1 엑손 1 내 유전적 다형성 및 생산 형질과의 연관성

문제 제기
인도 아대륙의 저명한 토착 낙농 품종인 Sahiwal 소는 열대 적응력과 질병 내성으로 가치를 인정받고 있다. 그러나 표현형 평가에 의존하는 전통적인 육종 프로그램은 경제적으로 중요한 형질의 낮은~중간 정도의 유전력, 형질 발현의 지연, 그리고 후손 검정(progeny testing)의 높은 비용으로 인해 한계에 직면해 있다. 유전적 개량을 가속화하기 위해 생산 및 번식 효율과 연관된 분자 마커를 식별하는 후보 유전자 접근법이 제안되고 있다. Peptidoglycan Recognition Protein 1 (PGLYRP-1) 유전자는 보형 염색체 18번에 위치하며, 선천성 면역, 유방 건강에서의 역할, 그리고 다른 소 집단(예: 홀스타인 및 무라 검소)에서 확인된 체세포 점수 및 유량과의 연관성 때문에 관심 대상인 후보 유전자이다. 그럼에도 불구하고, Sahiwal 소에서 PGLYRP-1 다형성에 관한 데이터는 제한적이다. 본 연구는 Sahiwal 소 집단에서 PGLYRP-1 유전자의 엑손 1 내 다형성을 조사하고, 생산 형질(305일 유량, 총 유량, 피크 유량, 락테이션 기간, 건유 기간, 효율 지수) 및 번식 형질(초산 연령)과의 연관성을 평가하는 것을 목적으로 하였다.

연구 방법론

• 연구 대상: 본 연구는 인도 라라 라지파트 라이 수의학 및 동물과학 대학교(LUVAS, Hisar)의 조직화된 농장의 Sahiwal 소 100두를 대상으로 하였다.
• 데이터 수집: 첫 번째 락테이션 형질에 대한 표현형 데이터를 농장 이력서에서 추출하였다. 분만 연도는 두 시기(2013–2020년 및 2021–2025년)로 그룹화하였고, 분만 계절은 네 가지 범주(여름, 몬순, 가을, 겨울)로 분류하였다.
• 분자 분석:
  • DNA 추출: Maxwell RSC Genomic DNA Extraction Kit를 사용하여 경정맥 혈액 샘플로부터 게놈 DNA를 분리하였다.
  • PCR 증폭: Wang 등(2013)이 이전에 보고한 프라이머를 사용하여 PGLYRP-1 유전자의 엑손 1을 포함하는 796 bp 단편을 증폭하였다.
  • 유전자형 결정: G+102C SNP를 검출하기 위해 ApaI 제한 효소를 이용한 PCR-제한 효소 절단 단편 길이 다형성(PCR-RFLP)을 수행하였다. 절단 패턴은 2% 아가로스 겔 전기영동을 통해 시각화되었다.
  • 염기서열 분석: PCR-RFLP 결과를 확인하고 Bos indicus 참조 서열(JN085440.1)에 대한 뉴클레오타이드 치환을 검증하기 위해 10개의 대표적인 앰플리콘에 대해 생거 시퀀싱(Sanger sequencing)을 수행하였다.
• 통계 분석:
  • 유전자형 및 대립유전자 빈도를 계산하였으며, 카이제곱 검정을 사용하여 하디-바인베르크 평형을 평가하였다.
  • 비유전적 요인의 영향은 최소제곱 혼합 모델(least-squares mixed model)을 사용하여 분석하였다.
  • 유전자형-형질 연관성은 일반 선형 모델(GLM)을 사용하여 평가하였다.

주요 결과

• 유전적 변이: 본 연구에서는 GG와 GC 두 가지 유전자형을 식별하였다. CC 유전자형은 관찰되지 않았다.
  • 유전자형 빈도: GG (0.97), GC (0.03).
  • 대립유전자 빈도: G (0.985), C (0.015).
  • 하디-바인베르크 평형: 해당 집단은 하디-바인베르크 평형 상태에 있는 것으로 나타났다 (카이제곱 값은 유의하지 않음).
  • 서열 확인: 생거 시퀀싱을 통해 참조 서열 대비 g.345 위치에서의 C > G 뉴클레오타이드 전위(transversion)를 확인하였으며, 이를 통해 ApaI 제한 효소 부위 다형성을 검증하였다.
• 비유전적 요인: 분만 시기는 305일 유량(305DMY), 총 유량(TMY), 락테이션 기간(LL), 효율 지수 및 초산 연령(AFC)에 유의미한 영향을 미쳤다. 초기 시기(2013–2020년)에 분만한 개체들은 후기 시기에 비해 높은 생산성과 낮은 AFC를 보였다. 분만 계절은 피크 유량과 효율 지수에는 유의미한 영향을 미쳤으나, AFC나 건유 기간에는 유의미한 영향을 미치지 않았다.
• 유전자형-형질 연관성:
  • 통계적 유의성: PGLYRP-1 유전자형과 평가된 생산 또는 번식 형질 사이에서 통계적으로 유의미한 연관성은 발견되지 않았다.
  • 수치적 경향: 통계적 유의성에도 불구하고 수치적 경향은 관찰되었다. GG 개체(n=97)는 GC 개체(n=3)에 비해 락테이션 기간과 건유 기간이 짧고, 피크 유량이 높으며, 우유 생산 효율이 좋고, 초산 연령이 빨랐다. 반대로, 세 마리의 GC 개체는 수치상으로 더 높은 305일 유량과 총 유량을 나타냈다.
  • 통계적 검정력의 한계: 저자들은 GC 유전자형의 극도로 낮은 빈도(단 3마리)와 CC 유전자형의 완전한 부재가 유전자 효과를 탐지하기 위한 통계적 검정력을 심각하게 제한했음을 명시적으로 언급하였다.

의의 및 주장
본 논문은 PGLYRP-1 엑손 1 영역이 Sahiwal 소에서 다형성을 나타내지만, 이 특정 집단에서는 현재의 데이터가 생산 또는 번식 형질과의 통계적으로 유의미한 연관성을 뒷받침하지 못한다고 결론지었다. 저자들은 이러한 유의미한 결과의 부재를 주로 불균형한 유전자형 분포(GG의 압도적 우세)로 인한 것이라고 설명하였다.

본 연구의 의의에 대한 주장은 다음과 같다:

1. 기초 데이터: 중국 홀스타인과 같은 다른 품종에 비해 G 대립유전자가 우세함을 밝힘으로써, Sahiwal 소에서 PGLYRP-1 엑손 1의 다형성 빈도를 최초로 보고하였다.
2. 방법론적 검증: 이 특정 품종에서 이 로커스(locus)를 유전자형 결정하기 위한 PCR-RFLP 및 생거 시퀀싱의 유용성을 확인하였다.
3. 결과의 예비적 성격: 저자들은 관찰된 수치적 경향이 예비적임을 명시적으로 밝히고 있다. 그들은 이 패턴이 이형접합체 집단의 작은 표본 크기로 인해 기능적 유전자 효과의 증거로 해석될 수 없음을 경고한다.
4. 향-방향: 본 논문은 PGLYRP-1이 이 품종에서 마커 보조 선발(marker-assisted selection)을 위한 유효한 마커로 고려되기 위해서는 더 균형 잡힌 유전자형 분포를 가진 더 큰 규모의 다중 농장 Sahiwal 집단에서의 검증이 필요하다고 주장한다. 본 연구는 즉각적인 육종 프로그램 적용을 주장하는 것이 아니라, 이 로커스의 유용성을 결정하기 위한 추가 조사의 필요성을 강조하는 데 목적이 있다.