The impact of nationwide folic acid fortification on genetic variants associated with conotruncal heart defects

기술 요약: 전국적인 엽산 강화가 원추부 심장 결함(Conotruncal Heart Defects)과 관련된 유전적 변이에 미치는 영향

문제 제 ステ이트먼트 (Problem Statement)
선천성 심장 질환(CHD)은 선천적 기형 관련 이환율 및 사망률의 주요 원인이지만, 대부분의 사례에 대한 병인을 명확히 알 수 없으며, 이는 다유전자적 위험(polygenic risk)과 환경적 요인 사이의 복잡한 상호작용을 포함하는 것으로 보인다. 엽산 보충이 신경관 결손 및 원추부 심장 결함(CTHD)의 위험을 줄이는 데 효과적이라는 사실은 확립되어 있으나, 구체적인 메커니즘은 여전히 부분적으로만 이해되고 있다. 엽산 관련 유전자 변이와 CTHD를 조사한 기존의 유전 역학 연구들은 모체의 엽산 결핍을 설명하는 데 어려움이 있고 표본 크기가 작다는 이유로 결정적인 결과를 내놓지 못했다. 특히, 1998년 미국의 곡물 및 시리얼 엽산 강화 의무화 전후에 태어난 CTHD 아동들 사이에서 유전적 변이 분포를 직접 비교한 연구는 지금까지 없었다. 본 연구는 인구 수준의 영양 중재(엽산 강화)가 CTHD에 대한 유전적 위험 변이의 침투도(penetrance) 또는 표현도(expressivity)를 어떻게 변화시켰는지에 대한 이해의 공백을 다룬다.

방법론 (Methodology)
본 연구는 dbGaP 저장소를 통해 접근 가능한 소아 심장 유전체 컨소시엄(PCGC)의 유전체 시퀀싱 데이터를 활용하였다.

• 코호트 정의: 분석에는 유럽계 혈통으로 추정되는 1,285명이 포함되었다. 코호트는 세 그룹으로 계층화되었다:
  1. 전체 분석 그룹 (모든 참가자, n=1,285).
  2. CTHD 하위 집단 (n=534): 팔로사징(tetralogy of Fallot), 동맥간(truncus arteriosus), 대동맥간 전위(D-transposition of the great arteries), 우심실 유출로 이중 돌출(double outlet right ventricle), 그리고 B형 대동맥궁 단절(interrupted aortic arch type B)을 포함함.
  3. 비-CTHD 하위 집단 (n=751): 중격 결손(septal defects), 좌측 폐쇄성 병변(left-sided obstruction lesions), 단심실 결손(single ventricle defects)을 포함함.
• 노출 분류: 참가자들은 엽산 강화 정책에 따른 출생 연도를 기준으로 분류되었다. 1997년 또는 그 이전에 태어난 경우를 "강화 전(pre-fortification)"으로, 2000년 또는 그 이후에 태어난 경우를 "강화 후(post-fortification)"로 분류하였다. 정책 전환기를 고려하여 1998년과 1999년 출생자는 제외하였다.
• 데이터 처리: 시퀀싱 데이터는 GRCh38/hg38에 정렬되었으며, 배치 효과(batch effects)를 최소화하기 위해 공동으로 처리되었다. 품질 관리(QC) 필터로는 유전자형 품질 >20, 대립유전자 읽기 깊이 >10, 미세 대립유전자 빈도 >0.01, Hardy-Weinberg 평형 p > 1x10⁻⁶를 적용하였다.
• 통-계적 분석: PLINK v2.0을 사용하여 가법 유전 모델(additive genetic model) 하에서 로지스틱-Firth 하이브리드 회귀를 이용한 전전체 유전체 연관 분석(GWAS)을 수행하였다. 인구 층화(population stratification)를 통제하기 위해 상위 10개의 주성분(principal components)을 공변량으로 포함하였다. 암시적 유의성(suggestive significance)은 p < 1 x 10⁻⁵로 정의하였다.
• 기능적 농축 (Functional Enrichment): 각 모델의 상위 100개 순위 유전자를 대상으로 DAVID (v6.8)를 사용하여 유전자 세트 농축 분석을 수행하였다. 유의성은 Benjamini-Hochberg 허위 발견율(FDR) 조정을 거친 수정된 Fisher의 정확 검정(modified Fisher's exact tests)을 통해 평가되었다.

주요 결과 (Key Results)

• GWAS 결과:
  • 전체 그룹 (모델 1): 94개의 변이가 명목상으로 농축되었다 (p < 10⁻⁵). 대부분의 변이는 강화 전 시대에 농축되었으며 (OR < 1), 이는 특정 유전적 변이가 강화 이후 CCHD 사례에서 덜 빈번해지는 인구 전체의 변화를 시사한다.
  • CTHD 하위 집단 (모델 2): 26개의 변이가 식별되었다. 특히, DHRS3 유전자의 두 변이(rs755_1703, OR 2.10; rs6541043, OR 2.21)와 PPARGC1β 로커스(특히 인트론 내 miRNA인 MIR378A)의 네 변이는 강화 후 그룹에서 명목상으로 농축되었다 (OR 4.47–4.51).
  • 비-CTHD 하위 집단 (모델 3): 59개의 변이가 식별되었다. MIR548A1HG 및 NKX1-1의 변이는 모델 1과 유사하게 강화 전 그룹에서 농축되었으며, 이는 이들이 CTHD 특이적 위험보다는 더 넓은 범위의 심장 표현형에 대한 공유 위험 요인임을 시사한다.
  • 엽산 경로 유전자: 클래식하게 엽산 운반 또는 1-탄소 대사와 관련된 변이(예: MTHFR, MTRR, MTR)는 어떤 모델에서도 검출되지 않았다.
• 기능적 농축:
  • CTHD 그룹에서는 Rho 단백질 신호 전달, 세포 골격 역학, 세포 간 접합, 그리고 전사 조절을 포함한 농축된 생물학적 경로가 나타났으며, 이는 심장 발달 과정과 일치한다.
  • 비-CTHD 그룹에서는 염색질 결합(chromatin binding) 및 산화적 인산화(oxidative phosphorylation)가 농축된 용어로 나타났다.
  • 저자들은 최종 해석 시 p-값이 다중 검정(multiple testing)을 위해 교정되지 않았음을 언급하며, 이러한 농축 결과가 가설 생성적(hypothesis-generating)임을 명시하였다.

의의 및 주장 (Significance and Claims)
본 논문은 CTHD 사례가 1998년 미국 엽산 강화 의무화 전후로 유전적 변이를 직접 비교한 첫 번째 연구라고 주장한다. 본 연구는 엽산 강화가 엽산 반응성 로커스의 게놈 지형을 수정했을 수 있는 "자연 실험(natural experiment)" 역할을 한다고 상정한다.

• 유전자-영양 상호작용: 강화 후 CTHD 그룹에서 특이적으로 농축된 DHRS3 및 PPARGC1β 변이의 발견은 잠재적인 유전자-영양 상호작용을 시사한다. DHRS3는 레티노이드 항상성에 관여하며, 저자들은 엽산 결핍과 레티노산 불균형 모두 CTHD와 연관되어 있다는 점을 들어 엽산이 레티노이드 경로를 조절할 수 있다고 언급한다. PPARGC1β는 심장 대사 성숙과 관련이 있으며, 이의 연관성은 시간이 흐름에 따라 개선된 임상적 결과에 따른 생존 편향(survival bias) 또는 변화된 영양 환경과의 특정한 상호작용을 반영할 수 있다.
• 인구 변화: 광범위한 코호트 전반에서 많은 변이가 강화 후 그룹에서 덜 빈번하게 나타난 것은, 강화가 인구 수준에서 유전적 위험 인자의 침투도를 변화시켰을 수 있음을 시사한다.
• 신중한 해석: 저자들은 사용된 암시적 유의성 임계값(p < 10⁻⁵)과 개별 수준의 모체 엽산 데이터 부재를 근거로, 결과가 확증적이라기보다 "가설 생성적"임을 명시적으로 밝히고 있다. 또한, 관찰된 변화가 연구 기간 동안의 생존율 향상이나 모체의 특성 변화(예: 비만, 흡연)와 같은 세속적 경향(secular trends)의 영향을 받을 수 있음을 인정한다.

결론적으로, 본 연구는 엽산 강화가 CTHD에 대한 유전적 위험 변이의 분포에 영향을 미쳤을 수 있다는 증거를 제공하며, 이는 모체의 엽산 상태, 배아 레티노이드 수준, 그리고 심장 발달 경로 사이의 복잡한 상호작작용을 강조한다.