이 두 표본 멘델 무작위 분석(MR) 연구는 유전적으로 예측된 NETs 관련 염증 특성과 제2형 당뇨병(T2DM) 합병증 간의 잠재적 연관성을 확인하였습니다. 유전적으로 예측된 IL-6는 당뇨병성 관상동맥질환(CAD)과 역상관 관계를 보인 반면, NETs는 T2DM의 신장 및 말초 순환 합병증과 연관이 있는 것으로 나타났습니다. 민감도 분석 결과, 이러한 유의미한 연관성들의 안정성이 대체로 입증되었습니다.
연구 논문
이 두 표본 멘델 무작위 분석(MR) 연구는 유전적으로 예측된 NETs 관련 염증 특성과 제2형 당뇨병(T2DM) 합병증 간의 잠재적 연관성을 확인하였습니다. 유전적으로 예측된 IL-6는 당뇨병성 관상동맥질환(CAD)과 역상관 관계를 보인 반면, NETs는 T2DM의 신장 및 말초 순환 합병증과 연관이 있는 것으로 나타났습니다. 민감도 분석 결과, 이러한 유의미한 연관성들의 안정성이 대체로 입증되었습니다.
호중구 세포외 트랩(Neutrophil extracellular traps, NETs) 관련 염증 특성은 제2형 당뇨병(T2DM) 및 그 합병증에서 중요한 역할을 합니다. 특히, 핵심 염증성 사이토카인인 IL-6는 NETs의 형성과 T2DM 및 합병증의 발병 기전과 밀접하게 연관되어 있습니다. 본 연구는 멘델 무작위 분석(Mendelian Randomization, MR) 접근법을 사용하여 NETs 관련 염증 특성과 T2DM 및 그 합병증 간의 인과적 연관성을 탐구하는 것을 목표로 하였습니다. 본 연구는 T2DM 및 합병증에 대한 대규모 유럽 인구 기반 메타 분석을 포함하는 전전체 연관성 분석(Genome-Wide Association Studies, GWAS) 데이터를 활용한 이표본 MR 설계를 사용하였습니다. 주요 분석 방법으로는 역분산 가중(inverse variance weighted, IVW) 접근법을 사용하였으며, MR-Egger 회귀, 가중 중앙값(weighted median) 및 가중 최빈값(weighted mode) 방법으로 이를 보완하였습니다. 민감도 분석에는 결과의 강건성을 평가하기 위해 MR-Egger, MR-PRESSO, Cochran’s Q 및 leave-one-out 방법이 포함되었습니다. 연구 결과, 유전적으로 예측된 interleukin-6 (IL-6) 수치는 당뇨성 관상동맥 질환(CAD)과 역상관 관계를 보였습니다 (OR = 0.8997, 95% CI: 0.8257–0.9803, P = 0.0158). 또한, NETs는 신장 합병증을 동반한 T2DM (OR=0.97, 95% CI 0.9428–0.998, P = 0.0358) 및 말초 혈관계 합병증을 동반한 T2DM (OR = 1.0342, 95% CI 1.002–1.0673, P = 0.037)과 유의미한 연관성을 나타냈습니다. 유의미한 IVW 연관성에서 이질성이나 수평적 다표현성(horizontal pleiotropy)의 증거는 발견되지 않았습니다. 본 연구는 유전적으로 예측된 NETs 관련 염증 특성이 특정 T2DM 합병증과 연관되어 있음을 시사합니다. 유전적으로 예측된 IL-6는 당뇨성 CAD와 역상관 관계를 보인 반면, NETs는 T2DM의 신장 및 말초 혈관계 합병증과 연관이 있었습니다.
인슐린 저항성을 특징으로 하며 흔히 진행성 베타 세포 기능 부전을 동반하는 제2형 당뇨병(T2DM)은 전 세계적으로 유행병 수준에 도달했습니다1,2. 심혈관 질환, 신부전, 신경병증 및 망막병증을 포함하는 이러한 T2DM 합병증은 개인의 건강과 삶의 질뿐만 아니라, 전 세계 보건 의료 시스템의 자원을 압박하고 비용을 증가시키며 경제 전반에 상당한 부담을 줍니다3,4. T2DM의 다면적인 특성을 고려할 때, 그 병인론적 요인에 대한 세밀한 탐구가 필요합니다. 축적된 증거들은 호중구 세포외 트랩(NETs)이 다양한 대사 장애, 특히 T2DM과 관련된 병태생리에 관여함을 시사하고 있습니다.
호중구 외세포 트랩(Neutrophil extracellular traps, NETs)은 강력한 살균 특성을 가진 탈응축 DNA와 과립 단백질의 외세포 네트워크입니다5,6. cell-free DNA 및 시트룰린화된 histone H3와 같은 마커는 NETosis의 핵심 지표로 작용하며, 다양한 질환에서의 염증 활성을 반영합니다7,8. 축적된 증거들은 당뇨병이 NETs의 형성을 촉진하여 염증과 조직 손상을 확산시킨다는 것을 시사합니다9,10. 임상적으로, 상승된 NETs 마커는 여러 당뇨병 합병증 전반에 걸쳐 고혈당 및 부정적인 결과와 상관관계가 있습니다. 예를 들어, dsDNA 수치의 증가는 STEMI 환자의 포도당 조절 장애와 연관되어 있으며11, NETs 구성 성분은 당뇨병성 망막병증의 진행과 관련이 있고12 당뇨병성 족부 궤양의 절단 위험을 예측합니다13. 이러한 연관성에도 불구하고, NETs와 T2DM 및 그 합병증 발생 사이의 정확한 인과 관계는 아직 완전히 밝혀지지 않았습니다. 축적된 증거들은 대사 염증에서 IL-6와 NETs 사이의 양방향 상호작용을 강조합니다. 고혈당은 호중구를 NETosis가 일어나기 쉬운 상태로 준비시키며, IL-6는 이 에너지 의존적 과정의 강력한 유도제 역할을 합니다14. 임상적으로 IL-6는 제2형 당뇨병에서 상승된 NETs 마커(예: H3Cit 및 cell-free DNA)를 독립적으로 예측하며, 이는 혈전 형성 상태 및 당뇨병성 망막병증과 같은 미세혈관 합병증과 강한 상관관계가 있습니다15,16. 또한, 감염으로 유도된 NETosis는 IL-6의 trans-signaling을 강화하여 상호 조절 루프를 형성할 수 있습니다17. 종합적으로, 이러한 발견들은 특히 당뇨병 환경 내에서 IL-6와 NETs의 밀접하게 결합된 병태생리학적 특성을 강조합니다.
멘델 무작위 분석(Mendelian Randomization, MR)은 자연적으로 발생하는 유전적 변이를 도구 변수로 활용하여 수정 가능한 노출 요인이 건강 결과에 미치는 인과적 효과를 추정하는 강력한 통계적 방법으로, 복잡한 질병 경로에 대한 통찰을 제공합니다18. 멘델의 유전 법칙에 근거한 이 접근 방식은 감수 분열 중 대립유전자의 무작위 배분을 이용하여 관찰 연구에서 인과 관계를 추론할 수 있게 하며, 이를 통해 인과 관계를 설정하는 전통적인 역학적 방법의 일부 한계점을 극복할 수 있게 해줍니다19.
현재까지 유전적으로 예측된 NETs 관련 염증 특성과 제2형 당뇨병(T2DM) 및 그 합병증 사이의 인과 관계를 조사한 MR 연구 보고는 제한적이었습니다. 따라서 본 연구는 이표본 MR 분석 방법을 사용하여 NETs 관련 염증 특성과 T2DM 및 그 합병증 사이의 잠재적인 인과 관계를 탐구하는 것을 목표로 합니다.
연구 설계
본 연구에서는 NETs 관련 염증 특성과 T2DM 및 그 합병증 사이의 인과 관계를 탐구하기 위해 2-표본 MR 접근 방식을 활용하였습니다. MR 프레임워크는 세 가지 핵심 가정을 기반으로 합니다20: 1) 유전적 변이는 노출과 강력하게 연관되어야 하며, 2) 유전적 변이는 오직 노출을 통해서만 결과에 영향을 미쳐야 하고(다면성 배제), 3) 유전적 변이는 어떠한 교란 요인과도 연관되지 않아야 합니다. 본 연구 과정에 대한 자세한 개요는 그림 1을 참조하십시오. 본 원고는 MR 연구 보고를 위한 MR-STROBE 가이드라인을 준수하여 투명성과 재현성의 엄격한 표준을 보장하였습니다21.
데이터 소스
본 MR 분석에 활용된 전장 유전체 연관 분석(GWAS) 데이터는 GWAS 공공 데이터베이스에서 수집되었습니다. T2DM과 관련된 요약 수준의 GWAS의 경우, 74,124명의 환자와 824,006명의 대조군으로 구성된 유럽 인구 대상 최대 규모의 메타 분석과 9,978명의 환자와 12,348명의 대조군이 포함된 유럽 내 전향적 중첩 환자-코호트 연구를 통합하였습니다22. 신장 합병증을 동반한 T2DM 및 말초 혈관 합병증을 동반한 T2DM을 포함한 당뇨병 합병증의 경우, 특정 합병증이 있는 T2DM 환자를 환자군으로, T2DM이 없는 개인을 대조군으로 설정한 연구들의 GWAS 요약 통계량을 추출하였습니다. 데이터 소스에 대한 자세한 정보는 부록 표 1에서 확인할 수 있습니다.
NETs 관련 염증 특성에 대한 GWAS 데이터는 GWAS catalog에서 확보하였다( Supplementary Table 2 참조). 생물학적 특이성의 이질성을 해결하기 위해, 본 연구에서는 호중구 생물학에서의 기능적 역할에 따라 포함된 노출 요인들을 두 가지 뚜렷한 그룹으로 분류하였다: (1) 핵심 NETosis 관련 인자: 이 범주에는 NETs의 구조적 형성에 직접 관여하거나 염색질 응축 해제(chromatin decondensation)의 효소 과정에 관여하는 마커들이 포함된다. 구체적으로 본 연구에서는 NETs 자체, Myeloperoxidase (MPO), Neutrophil Elastase (NE) 및 MPO-DNA 복합체를 포함하였다. MPO와 NE는 히스톤 분해 및 염색질 응축 해제에 필수적인 효소이며, MPO-DNA 복합체는 NETs의 특이적 대리 마커이다23. (2) NETs 관련 염증 매개체: 이 범주에는 NETosis와 밀접하게 연결되어 있으면서 더 광범위한 염증 경로에도 관여하는 상위 조절자 또는 하위 작동체로서의 사이토카인 및 매개체들이 포함된다. 이 그룹에는 Interleukin-6 (IL-6)24, Tumor Necrosis Factor-alpha (TNF-α)25, Neutrophil gelatinase-associated lipocalin levels (NGAL)26 및 Cellular Communication Network Factor 1 (CCN1)27이 포함된다.
본 연구에 사용된 데이터는 오픈 액세스 데이터베이스 또는 이전 연구에서 발표된 GWAS에서 수집되었으므로, 본 연구를 위한 윤리적 승인은 필요하지 않았습니다.
NETs 관련 유전적 도구의 선택
NETs 형성은 기존의 GWAS에서 직접적으로 측정되지 않는 동적인 생물학적 과정이므로, 본 연구에서는 NETs 관련 염증 특성에 대한 유전적 도구를 식별하기 위해 유전자 기반 접근 방식을 채택하였다. MPO와 같이 NETs 형성 및 조절에 결정적인 역할을 하는 것으로 알려진 257개의 포괄적인 유전자 세트를 히스톤 및 호중구 과립 단백질을 인코딩하는 유전자와 함께 기존의 기전 연구 및 발표된 문헌28으로부터 큐레이션하였다.
도구 변수 선택
NETs 관련 염증 특성에 대한 강건한 유전적 도구를 식별하기 위해, 본 연구에서는 GWAS 요약 통계량을 이용한 순차적 필터링 워크플로를 구현하였습니다. 구체적인 운영 단계는 다음과 같습니다:
초기 스크리닝: P < 5 × 10⁻6의 유의성 임계값을 기준으로 노출 형질과 유의하게 연관된 유전 변이들을 추출하였습니다29.
소수 대립유전자 빈도(MAF) 필터링: 통계적 검정력을 확보하기 위해 MAF ≤ 0.01인 단일 염기 다형성(SNPs)은 제외하였다30.
연쇄 불평형(LD) 클럼핑: 연쇄 불평형(LD)으로 인한 교란을 완화하기 위해, 10,000 kb 윈도우 내에서 r2 < 0.001로 설정된 파라미터의 클럼핑 함수를 사용하여 독립적인 SNP를 선정하였다31.
도구 강도 평가: 남은 각 독립변수의 강도는 F-통계량을 사용하여 정량화하였으며, 이는 F = R 값으로 계산되었습니다.2 × (N-2) / (1-R2). F-통계량(F-statistic)을 갖는 SNP만이 > 약한 도구 변수 편향(weak instrument bias)을 최소화하기 위해 10개를 유지하였다.32.
프록시 SNP 식별 및 데이터 조화
결과 데이터셋에서 누락된 SNP를 처리하기 위해, 본 연구에서는 프록시 검색 및 데이터 조화 과정을 수행하였습니다:
프록시 치환: 노출 GWAS의 표적 SNP를 결과 GWAS에서 사용할 수 없는 경우, 본 연구에서는 1,000 Genomes Project 유럽 참조 패널을 기반으로 한 LDproxy() 함수를 활용하였다. 후보 프록시 SNP는 원래의 SNP와 높은 LD(r2> 0.8)를 나타내는 경우에만 선택되었다. 적합한 프록시를 찾지 못한 경우, 해당 SNP는 제외되었다.
데이터 조화: TwoSampleMR 패키지(R version 4.0.5)의 harmonise_data() 함수를 사용하여 노출 및 결과 데이터 세트를 정렬했습니다. 본 연구에서는 parameter action = 2로 설정하여 모든 SNP를 순방향 가닥(forward strand)으로 자동 정렬하고 가닥 방향이 모호한 회문 SNP(palindromic SNPs)를 제거했습니다.
검증: 조화 작업 후, 본 연구에서는 조화된 데이터세트를 수동으로 검토하여 노출 데이터와 결과 데이터 사이의 효과 대립유전자 빈도가 일치하는지 확인하였습니다.
MR 분석 및 민감도 테스트
TwoSampleMR 패키지(R 버전 4.0.5)를 사용하여 인과관계 추론을 수행하였습니다.
1차 및 2차 분석: 본 연구에서는 Inverse Variance Weighted (IVW) 방법을 주요 접근 방식으로 적용하였습니다33. 강건성을 확보하기 위해 MR-Egger34, weighted median 및 weighted mode 방법을 사용하여 보완 분석을 수행하였습니다35.
민감도 평가: IV 간의 이질성은 mr_heterogeneity() 함수를 통한 Cochran’s Q test36를 사용하여 평가하였습니다. 수평적 다표현형성(Horizontal pleiotropy)은 MR-Egger34 절편 테스트(mr_pleiotropy_test())를 통해 평가하였습니다.
이상치 검출: 본 연구에서는 잠재적 이상치를 검출하기 위해 MR-PRESSO 패키지를 사용하였다37. mr_presso() 함수를 1,000회의 시뮬레이션으로 실행하였다. 검출된 이상치 SNP(P < 0.05)를 제거하고, 결과의 안정성을 검증하기 위해 인과 추정치를 재계산하였다. 또한, 인과 관계가 단일 SNP에 의해 주도되지 않았음을 보장하기 위해 leave-one-out 분석을 수행하였다38.
통계적 보정
다중 검정을 고려하기 위해, MR 분석에서 얻은 P-값은 허위 발견율(False Discovery Rate, FDR) 방법을 사용하여 보정되었습니다. 이는 다음을 통해 구현되었습니다. PR의 .adjust() 함수와 method = 파라미터 "허위 발견율(False Discovery Rate, FDR)"보정된 P-값(P허위 발견율) < 0.05 미만인 경우를 통계적으로 유의한 것으로 간주하였다.
선택된 도구 변수
본 연구에서는 NETs 관련 염증 특성을 노출 변수로 하여 MR 분석을 수행하였습니다. 이를 위해 총 106개의 IV가 식별되었으며, 세부 사항은 부록 표 3에 제시되어 있습니다. 이들 IV의 F-통계량 평균값은 26.46이었으며, 범위는 최소 20.6에서 최대 250.4까지였습니다.
T2DM 및 그 다양한 합병증(당뇨성 케토산증, 신장 합병증, 안과적 합병증, 신경학적 합병증 및 말초 순환 합병증)을 결과 변수로 채택했을 때, 3개 SNP에 대한 정보를 요약 데이터에서 일치시킬 수 없었습니다. 본 연구에서는 프록시 SNP(Proxy SNPs)를 활용하였습니다.부록 표 4).
NETs 관련 염증 특성이 제2형 당뇨병(T2DM) 및 그 합병증에 미치는 인과적 연관성
IVW 방법의 주요 결과는 Table 1에 제시되어 있으며, Supplementary Table 5에는 MR-Egger, weighted median 및 weighted mode 접근법을 통해 도출된 보완적 MR 추정치가 상세히 나와 있습니다. 특히, IVW 결과에서는 interleukin-6 (IL-6) 수치와 당뇨병성 CAD 사이에 통계적으로 유의미한 연관성이 나타났으며(OR = 0.8997, 95% CI 0.8257–0.9803, P = 0.0158), 이는 weighted median 방법의 결과와 일치합니다. 이러한 발견은 IL-6가 당뇨병의 심혈관 합병증과 유전적으로 연관되어 있음을 시사합니다. 또한, IVW 분석 결과 NETs는 신장 합병증을 동반한 T2DM (OR = 0.97, 95% CI 0.9428–0.998, P=0.0358) 및 말초 순환 합병증을 동반한 T2DM (OR = 1.0342, 95% CI 1.002–1.0673, P = 0.037)과 유의미한 연관성이 있는 것으로 나타났습니다. 위에서 언급한 양성 결과에 대한 산점도와 포레스트 플롯은 Figure 2에 제시되어 있습니다.
MR Egger 회귀 분석 결과, 분석 내에서 수평적 다면발현의 증거가 나타나지 않았으며, 이에 대한 자세한 내용은 다음에서 확인할 수 있습니다. 표 2Cochran의 Q 검정에서 제시된 이질성 결과에 따라 (표 2) 및 깔때기 도표 (그림 3A–C), 그 중 어느 것도 ~하지 않음이 관찰되었습니다. P Q 통계량의 값이 일반적인 임계값인 0.05 미만으로 나타났으며, 이는 본 연구의 주요 관심 결과에 대한 독립변수들 사이에 유의미한 이질성이 없음을 나타냅니다.
또한, MR-PRESSO 분석(표 3) 결과, CCN1을 노출 변수로, T2DM을 결과 변수로 고려했을 때 이상치(rs139237904)가 식별되었으며, 해당 이상치를 제외한 후에도 결과는 통계적으로 유의하지 않았습니다. 마찬가지로, MPO가 노출 변수이고 T2DM이 결과 변수인 경우 두 개의 이상치(rs10916508 및 rs6544660)가 검출되었으며, 제거 후에도 유의하지 않은 결과가 지속되었습니다. NGAL을 노출 변수로, T2DM을 결과 변수로 했을 때는 하나의 이상치(rs10497219)가 발견되었으며, 제외 후에도 음성 결과가 유지되었습니다. TNF-α를 노출 변수로, 케토산증을 동반한 T2DM을 결과 변수로 한 경우 단일 이상치(rs10103048)가 존재했으며, 이를 제외해도 유의하지 않은 결과에는 변화가 없었습니다. 그러나 NGAL과 T2DM 및 신경학적 합병증의 연관성을 분석했을 때는 두 개의 이상치(rs10497219 및 rs1287048)가 식별되었으며, 이를 제거한 후 결과가 양성으로 전환되어 유의미한 연관성이 있음을 나타냈습니다. 'leave-one-out' 분석(그림 3D–F)은 이러한 발견의 견고함을 뒷받침하는 추가적인 근거를 제공합니다. 요약하자면, 본 MR 분석은 당뇨병 합병증의 병인에 있어 IL-6와 NETs의 잠재적 역할을 시사하며, 이에 대한 추가적인 연구가 필요함을 보여줍니다.
데이터 가용성:
본 연구 중에 생성되거나 분석된 모든 데이터는 본 논문과 부록 표 1, 부록 표 2, 부록 표 3, 부록 표 4 및 부록 표 5에 포함되어 있습니다.

그림 1: 본 MR 연구의 연구 설계. (A) MR 분석은 세 가지 핵심 가정에 기반합니다. (B) 연구 설계 개요. NETs: neutrophil extracellular traps; T2DM: 제2형 당뇨병; SNP: 단일 염기 다형성; MR: 멘델 무작위 분석; IVW: 역분산 가중법; MR-PRESSO: Mendelian Randomization Pleiotropy RESidual Sum and Outlier 이 그림의 확대 버전을 보시려면 여기를 클릭하십시오.

그림 2: T2DM 및 그 합병증에 대한 NETs 관련 염증 특성의 인과적 연관성. (A) IL-6 수치와 당뇨성 CAD 간의 연관성을 보여주는 산점도. (B) NETs와 신장 합병증을 동반한 T2D 간의 연관성을 보여주는 산점도. (C) NETs와 말초 혈관 합병증을 동반한 T2D 간의 연관성을 보여주는 산점도. (D) 당뇨성 CAD에서 IL-6 수치의 포레스트 플롯. (E) 신장 합병증을 동반한 T2D에서 NETs의 포레스트 플롯. F: 말초 혈관 합병증을 동반한 T2D에서 NETs의 포레스트 플롯. 이 그림의 더 큰 버전을 보시려면 여기를 클릭하십시오.

그림 3: 민감도 분석 결과. (A) 당뇨성 CAD에서 IL-6 수치의 깔때기 도표. (B) 신장 합병증을 동반한 T2D에서 NETs의 깔때기 도표. (C) 말초 순환 합병증을 동반한 T2D에서 NETs의 깔때기 도표. (D) 당뇨성 CAD에서 IL-6 수치의 하나-제외(Leave-one-out) 민감도 분석. (E) 신장 합병증을 동반한 T2D에서 NETs의 하나-제외 민감도 분석. (F) 말초 순환 합병증을 동반한 T2D에서 NETs의 하나-제외 민감도 분석. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
표 1: T2DM 및 그 합병증에 대한 NETs 관련 염증 특성의 일차적 인과 관계 이 표를 다운로드하려면 여기를 클릭하십시오.
표 2: 도구 변수의 이질성 및 다표현성 결과 이 표를 다운로드하려면 여기를 클릭하십시오.
표 3: MR-PRESSO 결과 이 표를 다운로드하려면 여기를 클릭하십시오.
보충 표 1: T2DM 및 그 합병증의 데이터 소스 이 파일을 다운로드하려면 여기를 클릭하십시오.
보충 표 2: NETs 관련 염증성 특성의 데이터 출처 이 파일을 다운로드하려면 여기를 클릭하십시오.
보충 표 3: 도구 변수(IVs)의 선택 이 파일을 다운로드하려면 여기를 클릭하십시오.
보조 표 4: MR 분석에 사용된 프록시 SNP 목록 이 파일을 다운로드하려면 여기를 클릭하십시오.
보충 표 5: NETs 관련 염증 특성과 T2DM 합병증 간의 인과적 연관성에 대한 상보적 MR 추정치 이 파일을 다운로드하려면 여기를 클릭하십시오.
본 연구는 두 표본 멘델 무작위 분석(MR) 접근 방식을 사용하여 NETs 관련 염증 특성과 T2DM 및 그 합병증 사이의 연관성을 조사합니다. 주요 결과에 따르면, 유전적으로 예측된 IL-6 수치와 당뇨병성 CAD 사이에 통계적으로 유의미한 역연관성이 있음이 시사되었습니다. 이러한 결과는 혈중 IL-6가 당뇨병의 심혈관 합병증에 직접적인 보호 효과를 나타내기보다는, IL-6 신호 전달 경로의 복잡한 유전적 구조를 반영하는 것일 수 있습니다. 또한, 본 연구는 T2DM 환자의 신장 및 말초 순환 합병증을 포함한 다양한 당뇨병 합병증과 NETs 사이의 연관성을 밝혔습니다. 이러한 연관성은 당뇨병 관련 합병증의 병태생리에서 NETs의 중요성을 강조하며, 이러한 염증 경로를 표적으로 하는 향후 치료 전략에 정보를 제공할 수 있습니다.
IL-6와 당뇨병성 CAD 사이의 연관성은 상당한 논쟁의 대상이 되어 왔으며, 여러 연구에서 IL-6의 역할에 대해 다양한 관점을 제시하고 있습니다. 본 연구는 IL-6가 당뇨병 관련 심혈관 합병증과 복잡하고 부정적인 연관성이 있을 수 있음을 시사합니다. Cupido et al.39는 IL-6 신호 전달의 교란이 심혈관 대사 질환 위험의 감소와 연결된다는 점을 들어 이러한 결과를 뒷받침합니다. 반대로, Indumathi et al.40는 당뇨병 환자에게서 IL-6의 저메틸화와 과발현을 관찰하였으며, 이는 CAD에서 병원성 역할을 한다는 점을 시사합니다. 이는 IL-6 수치의 상승을 당뇨병성 CAD의 전염증성 유도 인자 및 예후 마커로 식별한 Shrivastav et al.41 및 Chittaragi et al.42의 결과와 일치합니다. 그러나 Liu et al.43는 죽상경화증에서 IL-6와 같은 아디포카인의 이중적인 성격을 강조하며, 이러한 상황 의존적 이중성이 연구마다 상충하는 결과가 나타나는 이유일 가능성이 높다고 설명합니다. 이러한 결과의 불일치는 여러 요인에 기인할 수 있습니다. 첫째, 연구 설계, 환자군 및 IL-6 수치 측정 방법의 이질성이 서로 다른 결론으로 이어졌을 수 있습니다44. 둘째, IL-6와 CAD 발생 사이의 시간적 관계는 복잡하며 인과적 성분과 반응적 성분을 모두 포함할 수 있는데, 이는 단면 연구나 일부 종단 연구에서는 완전히 포착되지 않았을 수 있습니다45.
본 연구와 여러 다른 연구들이 IL-6와 당뇨병성 CAD 사이의 유의미한 연관성을 시사하고 있지만, 이러한 관계의 정확한 성격은 여전히 복잡하며 다양한 생물학적 및 방법론적 요인에 의해 영향을 받을 수 있습니다. 본 MR 분석에서 유전적으로 예측된 IL-6 수치와 당뇨 관련 CAD 사이에 통계적으로 유의미한 음의 연관성(OR=0.8997)이 나타났다는 점은 주목할 만하며, 이는 IL-6가 심혈관 병리를 촉진한다는 일반적인 견해와 대조됩니다. 이러한 불일치는 염증 및 조직 복구 과정에서 IL-6가 갖는 복잡한 이중 역할46 때문이거나, MR 분석에 사용된 유전적 도구가 IL-6의 전신적 효과와는 다른 특정 생물학적 신호 전달 경로를 반영하고 있을 가능성에 기인할 수 있습니다. 본 연구에서 관찰된 음의 연관성을 단순히 IL-6가 심장 보호 효과를 가진다고 해석해서는 안 되며, 오히려 이는 IL-6 수용체 신호 전달의 기저에 있는 복잡한 유전적 메커니즘을 반영하는 것일 가능성이 높습니다. 또는, 관찰된 음의 연관성이 기존의 관찰 연구에서는 완전히 포착되지 않은 염증 네트워크 내의 보상 기전이나 음성 피드백 루프의 영향을 받았을 수도 있습니다. 향후 연구에서는 잘 설계된 종단적 연구와 CAD 발병 기전의 다요인적 성격을 고려함으로써 이러한 복잡성을 규명하는 것을 목표로 해야 합니다.
당뇨병 합병증, 특히 신장 질환의 맥락에서 NETs의 역할은 염증과 조직 손상의 중요한 매개체로 점점 더 인식되고 있습니다. Magaña-Guerrero et al.12은 당뇨망막병증에서 NETs를 고혈당 및 신부전과 연결 지으며, 이를 질병 진행의 표지자로 정의했습니다. 당뇨병성 신장 질환(DKD)의 경우, Liu et al.47은 midkine의 상향 조절이 NETs 형성을 촉진하여 신장 손상을 유발함을 입증하였고, Zheng et al.48은 NETs가 사구체 내피 기능 장애와 파이롭토시스를 유도하며, 이러한 효과가 DNase I 처리에 의해 완화됨을 보여주었습니다. 이러한 실험적 결과 외에도, 최근의 MR 및 전사체 분석을 통해 특정 NETs 관련 유전자(예: CLIC3, GBP2)가 증식성 당뇨망막병증의 인과적 위험 요인임이 확인되었습니다28. 흥미롭게도, 이러한 유전적 근거는 신장 합병증과의 부정적인 연관성을 보인 앞선 결과들과 뚜렷한 대조를 이루며, NETs의 조직 특이적 메커니즘을 부각시킵니다. 종합적으로, 본 연구는 T2DM에서 NETs와 말초 순환 합병증 사이의 유의미한 인과 관계를 입증함으로써 이러한 전반적인 연구 지형을 보완합니다. 당뇨병성 신병증에서 NETs의 역할에 영향을 미칠 수 있는 요인들의 복잡한 상호작용을 밝히기 위해서는 추가적인 연구가 필수적입니다.
본 연구에서는 IL-6와 당뇨병성 CAD 사이의 잠재적 연관성뿐만 아니라, NETs와 다른 당뇨병 합병증 사이의 유의미한 상관관계를 확인하였으나, 이 두 염증 매개체 간의 상호작용에 대해서는 추가적인 조사가 필요합니다. 본 연구에서는 IL-6에서 NETs로 이어지는 인과 경로를 검증하기 위한 전용 MR 분석을 수행하지 않았으나, 상당한 양의 실험적 근거는 IL-6가 NETs 관련 염증 특성의 상위 구동 인자로 작용할 수 있음을 시사합니다.
생물학적으로 IL-6는 급성기 반응의 중추적인 사이토카인입니다.49 그리고 호중구가 NETs를 형성하도록 프라이밍하는 것으로 나타났습니다.50당뇨병의 고혈당 미세환경에서 흔히 관찰되는 IL-6 수치의 상승은51호중구 내의 JAK/STAT3 신호 전달 경로를 활성화할 수 있으며52이러한 활성화는 peptidylarginine deiminase 4 (PAD4)의 전사 상향 조절에 있어 결정적인 단계입니다.53히스톤 시트룰린화 및 염색질 응축 해제를 담당하는 효소인54, NETosis의 특징54따라서 전신적 IL-6 상향 조절이 NETs의 과도한 형성을 촉진하고, 이것이 다시 혈관 내피 세포 손상과 혈전증을 악화시켜 당뇨병성 관상동맥질환(CAD) 및 신장 병변을 포함한 다양한 당뇨병 합병증으로 이어진다는 순차적 병리 축을 가설로 세우는 것은 생물학적으로 타당합니다. IL-6의 유해한 영향이 NETs를 통해 매개되는 정도를 공식적으로 정량화하기 위해서는 다변량 멘델 무작위 분석(MR) 또는 매개 분석을 이용한 향후 연구가 필요합니다.
본 연구에서 제시한 분석 프레임워크는 향후 응용 가능성이 매우 높습니다. 방법론적으로, 이 워크플로우를 다변량 멘델 무작위 분석(Multivariable MR, MVMR)으로 확장하여 상관관계가 있는 노출 요인들의 독립적인 인과적 효과를 분석함으로써, 복잡한 생물학적 경로를 더욱 포괄적으로 다룰 수 있습니다. 또한, 향후 연구 단계에서 공국소화(colocalization) 분석을 통합한다면, 노출 요인과 결과 변수가 동일한 인과적 변이를 공유하는지 확인하여 연쇄 불평형 편향(linkage disequilibrium bias)의 위험을 더욱 줄일 수 있을 것입니다. 임상적으로, 이 접근법은 약물 타겟 검증을 위한 강력한 청사진을 제공합니다. 수평적 다표현성(horizontal pleiotropy)을 체계적으로 필터링하고 도구 변수의 타당성을 확보함으로써, 본 프레임워크는 임상적 성공 확률이 높은 치료 타겟의 우선순위를 정할 수 있습니다. 궁극적으로, 이러한 유전적 통찰력을 위험 예측 모델로 전환함으로써, 조기 개입의 혜택을 가장 많이 받을 수 있는 고위험 하위 그룹을 식별하는 정밀 의료 전략을 촉진할 수 있을 것입니다.
본 연구는 유전적 변이를 도구 변수(IV)로 활용하여 인과관계를 규명하는 2-표본 멘델 무작위 분석(MR) 접근 방식을 사용하며, 이를 통해 NETs 관련 염증 특성과 T2DM 및 그 합병증 사이의 관계를 탐구하기 위한 강력한 프레임워크를 제공한다. 이러한 접근 방식은 MR을 사용하여 교란 요인과 역인과관계를 완화함으로써 관찰 연구의 한계를 개선한다. IVW 추정치에만 의존하는 표준 MR 구현과 달리, 본 연구는 강건한 민감도 분석(가중 중앙값, MR-Egger)과 다면발현 확인 과정을 통합하였다. 이러한 다중 방법론적 프레임워크는 인과 추정치가 유효하지 않은 도구 변수에 의해 좌우되지 않도록 보장하여, 단일 방법 분석보다 더 높은 타당성을 제공한다. 또한, 데이터 조화 및 품질 관리를 위한 표준화되고 투명한 프로토콜을 통해 이 접근 방식의 재현성을 보장하며, 이는 복잡한 유전자 분석에서 자주 발생하는 "재현성 위기" 문제를 해결한다. 효율성과 적용 가능성 측면에서, 이 간소화된 파이프라인은 수동 데이터 큐레이션과 관련된 계산 부담을 줄여 여러 표현형을 신속하게 스크리닝할 수 있게 한다. 이러한 확장성 덕분에 본 프레임워크는 일반적인 분석뿐만 아니라 고처리량 표현체 전반 연관성 스캔(PheWAS)에도 광범위하게 적용 가능하며, 경직된 단일 단계 관찰 접근 방식보다 복잡한 생물학적 네트워크를 탐구하는 데 더 다재다능한 도구를 제공한다.
하지만 본 연구에는 몇 가지 제한점이 있습니다. 첫째, 유전적 도구를 사용하는 것은 선택된 변이들이 관심 있는 노출 및 결과와 단일 경로를 통해 연관되어 있다는 가정에 의존하므로, 생물학적 시스템의 복잡성을 완전히 포착하지 못할 수 있습니다. 둘째, 본 연구의 MR 분석은 주로 유럽인 집단의 데이터에 의존했기 때문에, 연구 결과를 다른 인종이나 지리적 집단으로 일반화하는 데 한계가 있을 수 있습니다. 셋째, 본 연구에서 NETs 관련 염증 특성과 T2DM 및 그 합병증 사이의 잠재적 연관성을 확인했지만, 구체적인 분자 조절 네트워크와 기저 생물학적 메커니즘은 향후 실험적 연구를 통해 추가적인 규명이 필요합니다. 유의미한 인과적 연관성은 당뇨병성 CAD, 신장 및 말초 순환 합병증과 같은 특정 합병증에 국한되었으며, 테스트한 모든 당뇨병 결과에서 일관되게 관찰된 것은 아니라는 점에 유의해야 합니다. 또한, MR은 교란 요인을 줄일 수 있지만, 측정되지 않은 환경적 요인이나 집단 층화와 같은 모든 잠재적 편향 원인을 완전히 설명할 수는 없습니다.
또한, 연구 결과의 일반화 가능성은 유전 데이터가 도출된 특정 인구 집단에 의해 제한될 수 있으며, 이러한 결과를 확정하기 위해서는 다양한 인구 집단에서의 재현 연구가 필요합니다. 본 연구는 대규모 GWAS 요약 통계량을 활용하였으나, 임상 코호트나 실험 모델로부터 얻은 보충적인 멀티오믹스 시퀀싱 데이터가 부재하여, 관찰된 인과적 연관성의 기저에 있는 분자적 기전 및 전사 변화를 검증하는 능력에 한계가 있습니다. 본 MR 워크플로우는 분석상의 어려움, 특히 대립유전자 조화(allele harmonization) 및 다면발현(pleiotropy) 문제를 해결하기 위해 특정 조정이 필요했습니다. 회문 구조 SNP(palindromic SNPs)의 가닥 모호성(strand ambiguity)을 완화하기 위해, 본 연구에서는 중간 대립유전자 빈도(>0.42)를 가진 변이들을 엄격히 제외하도록 파이프라인을 수정하였습니다. 나아가, 잠재적인 수평적 다면발현 문제를 해결하기 위해 MR-PRESSO 이상치 검사를 워크플로우에 통합하였습니다. 이를 통해 유효하지 않은 도구 변수를 반복적으로 제거함으로써, 최종 인과 추정치가 강건하게 도출되고 다면발현 이상치에 의해 좌우되지 않도록 보장하였습니다.
MR 결과의 타당성은 몇 가지 중요한 분석 단계에 크게 의존합니다. 첫째, IV의 엄격한 선택이 무엇보다 중요합니다. 엄격한 전유전체 유의성 임계값(P<5×10-6)을 준수하고 연쇄 불평형 클럼핑(r2<0.001)을 수행함으로써, 본 연구는 약한 도구 변수 편향의 위험을 최소화하고 선택된 SNP가 노출과 강력하게 연관되어 있음을 보장했습니다. 둘째, 인과 추정치의 강건성은 수평적 다형성 해결 여부에 달려 있습니다. 본 연구는 특히 MR-Egger, 가중 중앙값(weighted median), MR-PRESSO와 같이 서로 보완적인 여러 방법을 사용하여 필요한 안전장치를 마련했습니다. 이러한 방법들 간의 일관성과 유의하지 않은 MR-Egger 절편 값은 IV가 오직 노출을 통해서만 결과에 영향을 미친다는 가정을 뒷받침합니다. 유전적 변이는 NET 관련 염증 특성의 대리 지표 역할을 할 수 있지만, 생식세포 유전자형에서 전사 활성화 및 단백질 방출로 이어지는 조절 경로는 엄격하게 선형적이지 않은 복잡한 생물학적 과정을 포함하므로, 생체 내 NETosis의 역동적이고 실시간적인 변동을 완전히 반영하지 못할 수 있습니다. 재현성과 관련하여, 대규모의 공개 GWAS 요약 통계량을 사용함으로써 본 분석은 독립적으로 검증될 수 있습니다. 조화(harmonization) 및 이상치 제거를 위한 특정 파라미터를 포함하여 분석 파이프라인의 투명성을 확보함으로써, 이러한 결과의 정확한 복제가 가능합니다. 이러한 단계들은 종합적으로 본 연구의 내부 타당성을 강화하며, 향후 연구를 위한 재현 가능한 프레임워크를 제공합니다.
결론적으로, 본 연구는 2-표본 MR 접근법을 사용하여 NET 관련 염증 특성과 T2DM 및 그 합병증 사이의 연관성을 탐구하였습니다. 연구 결과는 IL-6 수치와 당뇨병성 CAD 사이에 유의미한 연관성이 있음을 시사하며, 이는 당뇨병의 심혈관 합병증에서 IL-6의 잠재적인 역할을 강조합니다. 또한, 본 연구는 NETs와 신장 및 말초 순환 합병증을 포함한 다양한 당뇨병 합병증 사이의 연관성을 밝혀냈으며, 이는 당뇨병 관련 질환의 병인에서 NET 관련 염증 특성의 중요성을 뒷받침합니다. 이러한 통찰은 T2DM에서 염증, 포도당 대사, 그리고 혈관 건강 사이의 복잡한 상호작용을 이해하는 데 기여합니다. 본 연구는 당뇨병 합병증에서 NETs의 역할에 대한 가치 있는 증거를 제공하는 한편, 관련된 생물학적 기전의 복잡성과 정확한 경로 및 잠재적 치료 표적을 규명하기 위한 추가 연구의 필요성을 인정합니다.
저자들은 상충하는 이해관계가 없음을 선언합니다.
본 연구는 Zhongshan 사회과학기술 연구 중점 프로젝트(Grant number: 2021B3013)와 Zhongshan 사회과학기술 연구 프로젝트(Grant number: 2022B1080)의 지원을 받았습니다. 연구 지원 기관은 연구 설계, 데이터 수집 및 분석, 출판 결정 또는 원고 작성에 관여하지 않았습니다.
| 이름 | 회사 | 카탈로그 번호 | 댓글 |
|---|---|---|---|
| CCN1 | GWAS 카탈로그 | GCST90277805 | http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90277001-GCST90278000/GCST90277805/harmonised/GCST90277805.h.tsv.gz |
| 당뇨병성 관상동맥질환 | GWAS 카탈로그 | GCST006405 | http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST006001-GCST007000/GCST006405/Fall_gwas_results_any_chd_tab.txt.gz |
| IL-6 | GWAS 카탈로그 | GCST90010146 | http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90010001-GCST90011000/GCST90010146/harmonised/33303764-GCST90010146-EFO_0007937.h.tsv.gz |
| MPO | GWAS 카탈로그 | GCST90087967 | http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90087001-GCST90088000/GCST90087967/harmonised/35078996-GCST90087967-EFO_0011039.h.tsv.gz |
| MPO-DNA | GWAS 카탈로그 | GCST90013658 | http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90013001-GCST90014000/GCST90013658/harmonised/33717105-GCST90013658-EFO_0005243.h.tsv.gz |
| MR-PRESSO 패키지 | GitHub (rondolab) | 버전 1.0 | 다면발현 보정 인과 추정치 검출에 사용되는 R 패키지. |
| NE | GWAS 카탈로그 | GCST90162366 | http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90162001-GCST90163000/GCST90162366/harmonised/GCST90162366.h.tsv.gz |
| 호중구 세포외 트랩 (NETs) | GWAS 카탈로그 | GCST90137414 | http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90137001-GCST90138000/GCST90137414/harmonised/GCST90137414.h.tsv.gz |
| NGAL (중성 호중구 지질글리코칼릭스 유사단백질) | GWAS 카탈로그 | GCST90161504 | http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90161001-GCST90162000/GCST90161504/harmonised/GCST90161504.h.tsv.gz |
| R 소프트웨어 | R Foundation for Statistical Computing | 버전 4.0.5 | 모든 분석에 사용된 주요 통계 컴퓨팅 환경. |
| 제2형 당뇨병 | FINN R9 | 제2형 당뇨병 | https://storage.googleapis.com/finngen-public-data-r9/summary_stats/finngen_R9_T2D.gz |
| 케토산증을 동반한 제2형 당뇨병 | FINN R9 | E4_DM2KETO | https://storage.googleapis.com/finngen-public-data-r9/summary_stats/finngen_R9_E4_DM2KETO.gz |
| 신경계 합병증을 동반한 제2형 당뇨병 | FINN R9 | E4_DM2NEU | https://storage.googleapis.com/finngen-public-data-r9/summary_stats/finngen_R9_E4_DM2NEU.gz |
| 안과적 합병증을 동반한 제2형 당뇨병 | FINN R9 | E4_DM2OPTH | https://storage.googleapis.com/finngen-public-data-r9/summary_stats/finngen_R9_E4_DM2OPTH.gz |
| 말초 순환 장애를 동반한 제2형 당뇨병 | FINN R9 | E4_DM2PERIPH | https://storage.googleapis.com/finngen-public-data-r9/summary_stats/finngen_R9_E4_DM2PERIPH.gz |
| 신장 합병증을 동반한 제2형 당뇨병 | FINN R9 | E4_DM2REN | https://storage.googleapis.com/finngen-public-data-r9/summary_stats/finngen_R9_E4_DM2REN.gz |
| TNF-a | GWAS 카탈로그 | GCST004426 | http://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST004001-GCST005000/GCST004426/harmonised/27989323-GCST004426-EFO_0004684.h.tsv.gz |
| TwoSampleMR 패키지 | GitHub (MRCIEU) | 버전 0.6.22 | 멘델 무작위 분석(Mendelian Randomization analysis)을 수행하는 데 사용되는 R 패키지. |
이 JoVE 논문의 텍스트 또는 그림 재사용 허가 요청
허가 요청