Occupational endotoxin exposure, DNA methylation changes, and immune-cell remodeling: an epigenome-wide association study in two cohorts

기술 요약: 직업적 엔도톡신 노출, DNA 메틸화 변화 및 면역 세포 리모델링

문제 제기
수백만 명의 섬유 및 농업 종사자들이 그람 음성균의 성분인 엔도톡신(리포다당체, LPS)을 포함한 유기 분진에 노출되어 있습니다. 이러한 노출은 기류 제한 및 폐 기능 저하를 포함한 급성 및 만성 호흡기 장애를 유발하는 것으로 잘 알려져 있습니다. 그러나 누적된 엔도톡신 노출이 면역 리모델링과 그에 따른 호흡기 기능 장애로 이어지는 구체적인 분자적 경로는 여전히 불분명합니다. 실험 연구들은 LPS가 후성유전적 조절을 변화시킨다는 점을 시사하지만, 장기적인 직업적 엔도톡신 노출에 따른 혈액 메틸로옴(methylome) 반응은 체계적으로 규명된 바 없습니다.

연구 방법론
저자들은 직업적 엔도톡신 노출과 연관된 DNA 메틸화(DNAm) 변화를 식별하기 위해 두 개의 코호트를 이용한 후성유전체 전장 연관성 연구(EWAS)를 수행하였습니다.

• 발견 코호트 (STWS): 상하이 섬유 종사자 연구(STWS)가 발견 단계로 사용되었습니다. 이 종단적 직업 코호트에는 281명의 엔도톡신 노출 면직물 노동자와 253명의 비노출 실크 노동자가 포함되었습니다. 1996년에 수집된 혈액 샘플은 Illumina MethylationEPIC v2 어레이(930K CpGs)를 사용하여 프로파일링되었습니다. 엔도톡신 노출은 상세한 작업 이력과 현장 구역 측정값(Limulus Amebocyte Lysate 분석법)을 결합하여 누적 엔도톡신 농도(EU/m³-year)를 산출하여 평가되었습니다.
• 검증 코호트 (SWHS): 상하이 여성 건강 연구(SWHS)는 인구 기반의 전향적 코호트로, 표적 병렬 검증을 제공했습니다. 686명의 신규 폐암 사례와 683명의 매칭된 대조군(모두 비흡연자)으로 구성된 중첩 환자-대조군 샘플이 분석되었습니다. 개별 노출 측정값은 사용할 수 없었기에, 노출은 발표된 문헌과 중국 국가 표준 직업 코드를 기반으로 한 직업 유래 서열 노출 강도 순위(저, 중, 고)를 사용하여 평가되었습니다.
• 분석 접근법:
  • EWAS: 선형 혼합 모델을 사용하여 STWS 내 면직물 노동자와 실크 노동자 간의 DNAm 베타 값을 비교하였으며, 연령, 성별, BMI, 흡연, 노출 중단 후 경과 연수 및 EpiDISH로 추정된 면역 세포 분율을 보정하였습니다.
  • 검증: STWS에서 식별된 신호들은 SWHS에서 서열 노출 순위를 사용하여 테스트되었습니다.
  • 다운스트림 분석: 저자들은 경로 풍부성 분석(GSEA 및 ORA), 면역 세포 조성(EpiDISH, CellDMC, eFORGE), DNAm 유래 염증 위험 점수(IRSs), 그리고 후성유전적 연령 가속화(EAA)를 평가하였습니다.
  • 표현형 연관성: STWS에서 엔도톡신 관련 DNAm 바이오마커와 호흡기 표현형(FEV1, FVC, 기류 제한) 사이의 연관성을 탐색하기 위해 종단적 혼합 모델을 사용하였습니다.

주요 기여 및 결과

1. 엔도톡신 관련 DMP 식별:

   • 발견 단계에서는 직업적 엔도톡신 노출과 연관된 36개의 후성유전체 전장 유의한 차별 메틸화 위치(DMPs)를 식별하였습니다 ($P < 5 \times 10^{-8}$).
   • 36개의 DMP 중 35개가 실크 노동자에 비해 면직물 노동자에서 저메틸화(hypomethylation)되어 나타났습니다.
   • 주요 신호로는 cg22740441 (PPP3CA), cg11248361 (AQP4-AS1), cg08009902 (MTA3)가 포함되었습니다.
   • 15개의 DMP의 경우, 노출 중단 후 경과 연수가 노출 효과와 반대 방향으로 메틸화 수준과 연관됨에 따라 메틸화 변화의 부분적 가역성이 관찰되었습니다.

2. 검증:

   • 어레이 차이로 인해 SWHS에서 테스트 가능한 30개의 DMP 중 18개가 일치하는 효과 방향을 보였습니다.
   • 4개의 DMP가 검증 코호트에서 명목적 유의성($P < 0.05$)에 도달하였습니다: cg11869499 (POLG), cg16290931 (TSKS), cg05019530 (EPIC1), cg26033526 (PSMB9/TAP1).

3. 면역 리모델링 및 경로 풍부성:

   • 면역 세포 조성: 면직물 노동자는 실크 노동자에 비해 유의하게 높은 호중구 비율(63.5% 대 58.9%, $P = 7.52 \times 10^{-10}$)과 낮은 NK 세포 및 기타 림프구 아형 비율을 보였습니다.
   • 세포 유형 특이성: CellDMC 분석은 주로 호중구($n=203$)와 조절 T 세포($n=1$)에서 나타나는 204개의 세포 유형 특이적 신호를 식별하였습니다.
   • 경로: GSEA는 면역 및 신호 전달 경로를 강조하였습니다. 가장 상위의 KEKK 경로는 IgSF 세포 부착 분자 신호 전달(FDR = $8.2 \times 10^{-6}$)이었습니다. 주요 경계 유전자는 ITGB1 및 PIK3CA였습니다. DMP에 매핑된 유전자들의 과잉 표현 분석(ORA) 또한 면역 및 항원 처리 경로를 가리켰습니다.

4. 염증 위험 점수(IRS) 및 호흡기 표현형:

   • 누적 엔도톡신 노출은 CRP 기반 DNAm IRS(예: IRS(Ligthart), IRS(Hillary))와 양의 상관관계를 보였으나, IL-6 기반 IRS(Stevenson)와는 상관관계가 없었습니다.
   • 높은 기저 CRP 기반 IRS는 더 낮은 FEV1 및 FVC와 연관되었습니다. 종단적 분석 결과, 기저 IRS는 낮은 FEV1과 연관되어 있었으나, FVC와의 연관성은 추적 관찰 기간 동안 점진적으로 더 부정적으로 변했습니다 (IRS가 1-SD 높을 때마다 연간 약 -1.8 ~ -2.1 mL 감소).
   • 여러 엔도톡신 관련 DMP 및 CRP 기반 IRS는 동시적 및 전향적 기류 제한 및 폐 기능 저하와 탐색적 연관성을 보였습니다.

의의 및 주장
본 논문은 직업적 엔도톡신 노출이 면역 리모델링(특히 호중구 비율 증가와 면역 신호 전달 및 세포 부착에 관여하는 부위의 저메틸화로 특징지어지는)을 반영하는 조정된 DNAm 변화와 연관되어 있다고 주장합니다. 이러한 결과는 엔도톡신이 CRP 기반 염증 위험 점수와의 연관성에서 알 수 있듯이, 부분적으로 염증 경로를 통해 작용할 수 있음을 시사합니다.

저자들은 혈액 DNAm이 과거 노출에 대한 상대적으로 안정적인 후성유전적 기록 역할을 하여 더 근접한 단백질체 마커를 보완한다는 "만남의 중간 지점(meet-in-the-middle)" 프레임워크 내에서 이 발견들을 배치합니다. 그들은 일부 DMP의 노출 중단 후 부분적 가역성이 이러한 후성유전적 표지의 역동적인 특성을 뒷받침한다고 언급합니다.

연구는 누적 엔도톡신 노출이 면역 세포 리모델링 및 염증 경로와 연결되어 있으며, 이것이 다시 호흡기 손상과 연관된다고 결론짓습니다. 그러나 저자들은 독립적인 재현과 기전 연구가 특정 기능과 임상적 관련성을 확립하기 위해 필요하다는 점을 명시하며 신중한 태도를 유지합니다. 또한, 본 연구는 아직 동료 검토를 거치지 않은 프리프린트이며 임상 지침으로 사용되어서는 안 된다는 점을 분명히 하며, 이 바이오마커들이 임상용으로 준비되었다는 주장을 명시적으로 피하고 있습니다.