Comparing genome-wide significant and chemosensory variants as instruments for dietary patterns in Mendelian randomisation

1. 문제 제기 (Problem Statement)

• 인과성 추론의 어려움: 식이는 심대사 질환의 주요 위험 요인이지만, 관찰 역학 연구는 교란 변수 (confounding) 와 역인과관계 (reverse causation) 의 영향을 받아 인과성을 입증하기 어렵습니다. 무작위 대조 시험 (RCT) 은 윤리적, 실용적 한계가 있습니다.
• 식이 패턴의 복잡성: 단일 영양소나 식품이 아닌 '식이 패턴' (예: 건강식, 불건강식, 페스카타리안 등) 은 구성 성분 간의 상호작용이 복잡하여, 기존 GWAS 에서 도출된 유전적 변이가 특정 식이뿐만 아니라 다른 대사 경로 (예: FTO, APOE 유전자 영역) 에도 영향을 미치는 수평적 다면성 (horizontal pleiotropy) 문제를 야기합니다.
• 도구변수 (IV) 의 한계: 기존 GWAS 기반의 유전적 변이는 통계적 유의성 (P < 5×10⁻⁸) 만을 기준으로 선정되어, 건강 상태에 따른 식이 변경 (역인과관계) 이나 교란 요인과 연관될 가능성이 있습니다. 반면, 화학수용체 (맛과 냄새) 관련 유전 변이는 식이 선택에 대한 생물학적 메커니즘이 명확하여 다면성 위험이 낮을 것으로 기대되지만, 복합 식이 패턴을 설명하는 데 충분한 통계적 검정력 (power) 을 가지는지 검증이 필요했습니다.

2. 방법론 (Methodology)

이 연구는 2-샘플 멘델 무작위화 (Two-sample MR) 프레임워크를 사용했으며, 주요 단계는 다음과 같습니다.

• 데이터 소스:
  • 노출 (Exposure): UK Biobank 기반의 4 가지 주성분 분석 (PCA) 도출 식이 패턴 (Unhealthy, Healthy, Meat-based, Pescatarian) 에 대한 GWAS 요약 통계.
  • 결과 (Outcome): 12 가지 심대사 지표 (BMI, 관상동맥질환, 지질 프로필, 혈압, 제 2 형 당뇨병, 당화혈색소, 공복 인슐린 등).
  • 인구: 유럽 계통 (European ancestry) 으로 제한하여 인구 구조 편향을 최소화.
• 도구변수 (IV) 선정 전략 (두 가지 접근법 비교):
  1. 전통적 IV (Conventional IVs): 식이 패턴 GWAS 에서 통계적으로 유의한 (P < 5×10⁻⁸) 독립 변이 선정.
  2. 생물학적 IV (Receptor IVs): 미각 및 후각 수용체 유전자의 비동일성 (nonsynonymous) SNP 1,214 개 중 식이 패턴과 연관된 변이 선정.
• 품질 관리 (Quality Control) 및 필터링:
  • PheWAS (Phenome-wide association study): 후보 변이가 초기 생활 특성 (출생, 아동기 체형), 사회경제적 지위 (SES), 생활 습관 (알코올, 흡연 등) 과 연관되는지 확인. 연관된 변이는 제거하거나 다변량 MR(MVMR) 로 보정.
  • Steiger 필터링: 변이가 결과보다 노출 (식이 패턴) 을 더 잘 설명하는지 확인하여 역인과관계 배제.
  • 다변량 MR (MVMR): 교육 수준, 소득, 신체 활동, 흡연, 아동기 체형 등을 공변량으로 포함하여 교란 효과 보정.
• 통계 분석:
  • 주 분석: 역가중치 (IVW) 방법 사용.
  • 민감도 분석: MR-Egger, 가중 중앙값 (Weighted Median), MR-PRESSO (이상치 제거) 를 통해 다면성과 이질성 평가.
  • Bonferroni 보정 적용 (α = 0.00152).

3. 주요 기여 (Key Contributions)

• 비교 연구 설계: 통계적 유의성 기반의 전통적 IV 와 생물학적 메커니즘 (화학수용체) 기반의 IV 를 동일한 식이 패턴과 결과에 대해 직접 비교하여 MR 연구에서의 도구변수 선정 전략을 평가했습니다.
• 엄격한 필터링 프로토콜 적용: 다면성과 역인과관계를 제거하기 위해 PheWAS, Steiger 필터링, MVMR 을 종합적으로 적용하여 MR 연구의 엄격성을 높였습니다.
• 방법론적 통찰: PCA 기반의 복합 식이 패턴을 MR 로 분석할 때, 화학수용체 변이가 가진 생물학적 특이성 (specificity) 이 통계적 검정력 부족으로 인해 실제 인과 효과 검출에 한계가 있음을 실증했습니다.

4. 연구 결과 (Results)

• 전통적 IV (Conventional IVs) 분석 결과:
  • Pescatarian (생선 위주) 식이 패턴: 공복 인슐린 (Fasting Insulin) 감소를 유의하게 예측 (βIVW = -0.10, P = 1.19×10⁻³). 이 결과는 민감도 분석 (MR-Egger, Weighted Median 등) 에서도 일관되게 나타났습니다.
  • Unhealthy (불건강) 식이 패턴: 수축기/이완기 혈압 및 당화혈색소 (HbA1c) 증가와 연관되었으나, 사용 가능한 IV 수가 매우 적어 (2 개) 다면성 보정이 불가능했고, 일부 변이가 카페인 섭취와 강하게 연관되어 있어 결과 해석에 주의가 필요했습니다.
  • 기타: 다른 식이 패턴과 심대사 지표 간에는 유의한 인과 관계가 확인되지 않았습니다.
• 화학수용체 IV (Receptor IVs) 분석 결과:
  • 결과: 모든 식이 패턴 - 결과 쌍에서 유의한 연관성이 발견되지 않았습니다 (Null findings).
  • 원인: 화학수용체 변이는 개별 식품 선택에는 강력하지만, 복합적인 '식이 패턴' 전체를 설명하는 분산 (variance explained) 이 매우 낮아 통계적 검정력이 부족했기 때문입니다.
• 다변량 MR (MVMR): Pescatarian 식이와 공복 인슐린의 연관성은 약한 도구변수 편향 (weak instrument bias) 이 존재했음에도 불구하고 방향성이 일관되었으나, Unhealthy 식이의 경우 신뢰구간이 0 을 포함하여 불확실성이 컸습니다.

5. 의의 및 결론 (Significance & Conclusion)

• MR 의 유효성 확인: 엄격하게 필터링된 전통적 IV 를 사용하여 잘 알려진 식이 - 질병 관계 (페스카타리안 식이와 인슐린 민감도) 를 재확인함으로써, MR 이 영양 역학에서 인과 추론을 강화하는 도구로 유용함을 입증했습니다.
• 생물학적 IV 의 한계: 화학수용체 변이는 다면성 위험이 낮다는 이론적 장점이 있었으나, 복합 식이 패턴을 설명하는 데 필요한 설명력이 부족하여 MR 분석에 활용하기에는 통계적 검정력이 미흡함이 드러났습니다. 이는 개별 식품이 아닌 '패턴'을 연구할 때 생물학적 변이 선정 전략의 한계를 보여줍니다.
• 향후 방향성:
  • PCA 기반 패턴 대신 지침 기반 (Guideline-based) 식이 점수 (예: Healthy Eating Index) 를 GWAS 에 활용하여 임상적 해석과 적용 가능성을 높여야 함.
  • 개별 식품에 대한 생물학적 변이 (예: 유당 분해 효소, 특정 맛 수용체) 를 집계하여 전체 식이 패턴을 모델링하는 새로운 통계적 방법 개발 필요.
  • 식이 요인 단독이 아닌, 다른 위험 요인들과의 복합적 상호작용을 고려한 다변량 접근의 중요성 강조.

요약: 본 연구는 엄격한 MR 분석을 통해 페스카타리안 식이가 공복 인슐린을 낮춘다는 인과적 증거를 제시했으나, 화학수용체 기반의 생물학적 도구변수는 통계적 검정력 부족으로 인해 복합 식이 패턴 연구에는 적합하지 않음을 보여주었습니다. 이는 영양 역학 연구에서 도구변수 선정과 식이 노출의 정의 (개별 식품 vs. 패턴) 가 인과 추론의 성패를 좌우함을 시사합니다.