Genetic determinants of cytokine production in activated human monocytes

본 연구는 366 명의 기증자로부터 유전체, 전사체 및 사이토카인 분비 데이터를 통합하여 단핵구 사이토카인 생성을 조절하고 지질 대사와 질병 위험 간의 연관성을 규명하는 CCR5-Δ32 와 OAS1 을 포함한 주요 유전적 좌위를 규명하였다.

핵심

당신의 몸의 면역 체계를 고도로 훈련된 보안 팀으로 상상해 보세요. 단핵구는 최전방에 서 있는 경비원들입니다. 그들이 침입자 (세균이나 바이러스 등) 를 발견하면 그냥 서 있지 않습니다. 대신 사이토카인이라는 화학적 "사이렌"을 방출하며 지원을 요청합니다. 이 사이렌들은 면역 체계의 다른 부분들이 깨어나 싸우도록 알립니다.

이 논문은 과학자들이 다음과 같은 질문을 던지는 탐정 이야기와 같습니다: "경보가 울렸을 때 왜 어떤 사람의 경비원은 다른 사람보다 더 크게 또는 더 작게 외치는 것일까? 그것은 그들의 훈련 (환경) 때문인가, 아니면 그들의 DNA(유전) 때문인가?"

연구자들이 발견한 바를 간단히 설명해 드리겠습니다:

1. 실험: 경비원들을 깨우기

과학자들은 366 명의 건강한 사람들로부터 혈액 샘플을 채취했습니다. 그들은 단핵구 (경비원) 를 분리하여 실험실 접시에 넣은 후, 두 가지 방식으로 공격을 시뮬레이션했습니다:

• LPS: 세균성 경보 (세균 감염 시뮬레이션).
• IFN-gamma: 바이러스성 경보 (바이러스 감염 시뮬레이션).

그들은 두 가지 것을 측정했습니다:

1. 외침: 접시에 방출된 사이토카인 단백질의 양.
2. 청사진: 그 외침을 만들기 위해 읽히는 RNA(유전 지시서) 의 양.

2. 큰 놀라움: 청사진이 항상 외침과 일치하지는 않음

보통 세포가 더 많은 지시서 (RNA) 를 읽으면 더 크게 외칠 것 (단백질 증가) 으로 예상할 수 있습니다. 하지만 과학자들은 이것이 항상 사실이 아니라는 것을 발견했습니다.

• 비유: 공장을 상상해 보세요. 때로는 관리자가 거대한 주문 목록 (RNA) 을 작성하지만, 공장 바닥이 느리거나 근로자들이 제품을 나중에 비축하기로 결정할 수 있습니다. 다른 때는 관리자가 작은 목록을 작성하지만, 근로자들이 매우 효율적으로 작동하여 막대한 양을 생산하기도 합니다.
• 발견: 많은 경우, 경보가 울리기 전 세포가 가지고 있던 "지시서" (RNA) 의 양은 그것이 얼마나 크게 외칠지 예측하지 못했습니다. 세포의 외침 능력은 종종 지시서가 작성된 후의 단계들, 즉 메시지가 외침으로 번역되는 속도나 외침이 전달되는 방식 등에 의해 통제되었습니다.

3. 유전적 "볼륨 조절기"

연구자들은 366 명 모두의 DNA 를 분석하여 특정 유전적 차이가 이러한 외침을 위한 "볼륨 조절기" 역할을 하는지 확인했습니다. 그들은 사이토카인 방출량을 크게 변화시키는 네 가지 특정 유전 위치를 발견했습니다:

• 조절기 #1 ("PDGF" 스위치): PDGFB 유전자 근처의 유전적 변이는 조광기 스위치처럼 작용했습니다. 이 유전자의 특정 버전을 가진 사람들은 자극을 받았을 때 덜 외쳤습니다. 흥미롭게도, 이 동일한 유전적 스위치는 궤양성 대장염(염증성 장 질환) 과 원발성 담즙성 경화증(간 질환) 의 위험 증가와 연관되어 있습니다. 더 "조용한" 경비원을 갖는 것이 실제로 이러한 특정 질환에 더 취약하게 만들 수 있는 것 같습니다.
• 조절기 #2 ("IL-1RA" 스위치): 또 다른 유전적 위치는 염증을 멈추게 하는 특정 "브레이크" 화학물질 (IL-1RA) 의 방출을 통제했습니다. 이는 세포가 자신의 지시서를 어떻게 읽는지와 연관되어 있었습니다.
• 조절기 #3 ("CCR5" 미스터리): 이는 HIV 에 대한 보호로 알려진 유명한 유전적 돌연변이 (CCR5-Δ32) 입니다. 이 연구는 이 돌연변이를 가진 사람들이 경보가 울렸을 때 두 가지 특정 화학물질 (MIP-1b 및 RANTES) 에 대해 더 크게 외쳤다는 것을 발견했습니다. 그들의 "수용체"(화학물질이 보통 들어가는 문) 가 고장 나 있기 때문에, 화학물질이 흡수되는 대신 외부에 쌓이는 것 같습니다. 이 돌연변이는 일부 질환에 대한 위험은 낮추지만 다른 질환에 대한 위험은 높여, 생물학이 트레이드오프임을 보여줍니다.
• 조절기 #4 ("IP-10" 스위치): IFNB1 유전자와 연관된 유전적 위치는 IP-10 의 방출을 통제했습니다. 이는 초기에 제공된 유전적 지시서 (경보 후 2 시간) 가 훨씬 나중에 (경보 후 24 시간) 외침을 결정한다는 것을 보여주었습니다.

4. "혈소판" 단서

과학자들은 두 가지 특정 화학물질 (BDNF 및 PDGF-BB) 에 대한 외침 양이 혈소판(응고를 돕는 작은 혈액 세포) 이 단핵구에 붙어 있는 수와 연관되어 있음을 발견했습니다.

• 비유: 경비원 (단핵구) 이 건설 노동자 (혈소판) 무리 옆에 서 있는 것과 같습니다. 근처에 건설 노동자가 많을수록 경비원은 더 크게 외쳤습니다. 이는 "소음"이 경비원 자신의 DNA 에서만 나오는 것이 아니라, 그들이 함께한 회사에서 비롯된 것임을 시사합니다.

5. "지질" 연결

연구자들은 RNA 와 단백질 외침을 어떻게 연결하는지 통제하는 유전자들을 살펴봤을 때 놀라운 패턴을 발견했습니다: 지방 대사였습니다.

• 비유: 외침의 볼륨을 통제하는 "관리자"들은 종종 세포가 지방 (지질) 을 처리하는 방식과 관련된 유전자들이었습니다. 마치 세포의 외침 능력이 "연료 탱크"가 얼마나 잘 관리되는지에 달려 있는 것처럼요. 이는 면역 체계의 외침을 신체가 지방을 처리하는 방식과 직접적으로 연결시켜, 심장 질환과 염증을 이해하는 데 큰 의미가 있습니다.

6. "코로나" 연결

마지막으로, 그들은 OAS1이라는 유전자를 살펴보았는데, 이는 중증 COVID-19의 위험 인자로 알려져 있습니다.

• 발견: 이 유전자의 특정 버전은 RNA 생성량을 바꾸지는 않았지만, RNA 와 단백질 외침 사이의 연결을 끊었습니다. 이 버전을 가진 사람들에서 세포는 지시서를 정상적으로 읽었지만, 실제 외침과의 연결이 10 가지 다른 사이토카인에 대해 끊어졌습니다. 이는 이 유전자가 COVID-19 를 더 악화시키는 이유가 면역 체계의 "외침" 타이밍이나 조율을 혼란시키기 때문일 수 있음을 시사합니다.

요약

이 논문은 면역 체계의 반응이 DNA 기반의 단순한 "켜기/끄기" 스위치가 아님을 보여줍니다. 그것은 다음과 같은 복잡한 오케스트라입니다:

1. 유전학이 볼륨 조절기를 설정합니다.
2. 타이밍이 중요합니다 (초기에 일어나는 일이 나중에 일어나는 일에 영향을 줍니다).
3. 맥락이 중요합니다 (주변에 있는 다른 세포들이 소음을 바꿉니다).
4. 지방 대사가 볼륨을 통제하는 숨겨진 역할을 합니다.

가장 중요한 점은, 이 연구가 질병과 싸우는 신체의 방식을 이해하려면 단순히 "지시서" (RNA) 만 보는 것만으로는 부족하며, 전체 그림을 보려면 실제 "외침" (단백질) 을 들어야 함을 증명한다는 것입니다.

기술 요약

기술 요약: 활성화된 인간 단핵구에서의 사이토카인 생산에 대한 유전적 결정인자

문제 및 동기
단핵구는 선천성 면역 반응의 중추 매개체로서, 사이토카인 분비를 통해 염증을 조절합니다. 전장 유전체 연관 분석 (GWAS) 과 발현 양적 형질 유전자좌 (eQTL) 매핑은 유전적 변이가 유전자 발현 및 질병 위험과 어떻게 연결되는지를 성공적으로 규명해 왔으나, 유전적 변이와 기능적 사이토카인 분비 간의 관계는 여전히 잘 규명되지 않았습니다. 이전 연구들은 조절 결정인자가 세포 유형과 활성화 상태에 따라 크게 달라지는 맥락 의존적 특성을 가지며, RNA 발현이 종종 단백질 양을 정확하게 예측하지 못한다고 지적했습니다. 또한, 전혈을 대상으로 한 대규모 단백질 QTL(pQTL) 연구는 세포 유형 특이적 또는 활성화 특이적 조절 기전을 놓칠 수 있습니다. 본 연구는 선천성 면역 활성화 후 주요 인간 단핵구에서의 사이토카인 분비에 대한 유전적 및 전사적 결정인자를 매핑하고, 유전체, 전사체, 단백질체 데이터를 통합하여 면역 반응의 조절 구조를 정의하는 것을 목표로 합니다.

방법론
본 연구는 366 명의 건강한 유럽계 기원 기증자 코호트를 활용했습니다. 주요 단핵구를 분리하여 네 가지 조건에 노출시켰습니다: 무처리 (UT), 2 시간 지질다당류 (LPS) 자극, 24 시간 LPS 자극, 그리고 24 시간 인터페론-감마 (IFN$\gamma$) 자극.

• 다중 오믹스 통합: 각 샘플에 대해 저자는 멀티플렉스 Luminex 분석을 사용하여 28 가지 사이토카인의 분비를 정량화하고, Illumina HumanHT-12 v4 BeadChip 어레이를 통해 유전자 발현을 측정했으며, Illumina HumanOmniExpress-12v1.0 을 이용한 전장 유전체 유전형 분석을 수행한 후 Haplotype Reference Consortium 에 대한 임피테이션을 실시했습니다.
• 통계적 분석:
  • 사이토카인 역학: 조건 간 사이토카인 분비 및 RNA 발현의 유의한 변화를 식별하기 위해 선형 회귀 모델을 사용했습니다.
  • 유전적 매핑: 사이토카인 분비와 관련된 유전적 로커스를 식별하기 위해 다변량 (MANOVA) 및 단변량 GWAS 를 수행했습니다. 유의성 임계값은 MANOVA 의 경우 $P < 1.7 \times 10^{-8}$, 단변량 검정의 경우 $P < 6.0 \times 10^{-10}$로 설정했습니다.
  • 전사적 상관관계: RNA 발현이 단백질 분비를 예측하는 능력을 선형 회귀를 통해 평가했습니다.
  • 공발현 QTL(coExQTL): 저자는 cis-eQTL 을 매핑하고 유전형과 사이토카인 분비 간의 상호작용을 테스트하여 RNA 전사와 단백질 분비 간의 관계를 수정하는 유전적 변이를 식별했습니다.
  • 공위치 분석: 사이토카인 분비 QTL 이 알려진 질병 GWAS 로커스와 공통된 원인 변이를 공유하는지 확인하기 위해 베이지안 공위치 분석을 수행했습니다.

주요 결과

1. 사이토카인 분비의 역학:

   • LPS 자극은 28 가지 사이토카인 중 19 가지의 변화를 유도한 반면, IFN$\gamma$는 5 가지를 유도했습니다.
   • 두 가지 뚜렷한 시간적 패턴이 관찰되었습니다: 초기 지속적 변화 (예: TNF$\alpha$, IL-1$\beta$) 와 후기 고립적 변화 (예: GRO$\alpha$, MCP-1).
   • 사이토카인 분비는 매우 조율되어 있으며, 계층적 클러스터링을 통해 공동 조절되는 사이토카인 그룹이 드러났습니다.
   • RNA-단백질 불일치: 기저 RNA 발현은 유도된 사이토카인 분비와 거의 또는 전혀 상관관계가 없었습니다. 자극된 RNA 발현은 일부 쌍 (예: CCL2:MCP-1) 에 대해 단백질 분비를 예측했으나, 분산의 일부만 설명했습니다 ($R^2$ 최대 0.56). 특히 일부 사이토카인 (예: IL-2, IL-18) 은 유의미한 동시 RNA 유도 없이 분비되었고, 일부 전사체 (예: PDGF-BB, IL-8) 는 해당 단백질 변화 없이 증가하여 번역, 운송, 분해와 같은 전사 후 조절이 주요 결정인자임을 시사했습니다.

2. 전사 상태 예측인자:

   • 기저 전사 상태는 일반적으로 사이토카인 분비를 예측하지 못했으나, IFN$\gamma$ 자극 후 IL-1$\beta$ 분비를 예측하는 JUN 발현을 제외했습니다 (아마도 인플라마솜 프라이밍을 통해).
   • 혈소판 특이적 유전자로 풍부해진 기저 서명 (단핵구 - 혈소판 응집체, MPAs 를 나타냄) 은 모든 조건에서 BDNF 및 PDGF-BB 분비와 유의하게 상관관계가 있었습니다.

3. 사이토카인 분비의 유전적 결정인자:
   본 연구는 네 가지 전장 유전체 유의성 로커스를 식별했습니다:

   • PDGFB (cis): PDGFB 상류의 로커스 (rs5757584) 는 IFN$\gamma$ 및 LPS 자극 후 PDGF-BB 분비 감소를 associat 했습니다. 이 변이는 자극된 단핵구에서 PDGFB의 eQTL 로 작용하며, 궤양성 대장염 및 일차성 담관 경화증에 대한 GWAS 위험 로커스와 공위치합니다.
   • IL1RN (cis): IL1RN 근처의 로커스 (rs11123160) 는 LPS 자극 24 시간 시점에 특이적으로 IL-1RA 분비와 연관되었습니다. 이 변이는 자극된 단핵구에서 IL1RN의 eQTL 입니다.
   • CCR5-Δ32 (trans): CCR5-Δ32 결실을 태그하는 rs113010081 변이는 전사 조절자로 작용했습니다. 이는 IFN$\gamma$ 후 MIP-1$\beta$ 분비와 LPS 후 RANTES 분비를 증가시켰습니다. 이 효과는 CCL4/CCL5 전사에 독립적이었으며, 수용체 매개 피드백 또는 제거와 관련된 기전을 시사합니다.
   • IFNB1/CXCL10 (trans): IFNB1 근처의 로커스 (rs13296842) 는 LPS 24 시간 시점의 IP-10(CXCL10) 분비와 연관되었습니다. 이는 변이가 초기 (2 시간) 에 IFNB1의 cis-eQTL 로 작용하여 하류 CXCL10 발현 및 단백질 분비를 후기 (24 시간) 에 유도하는 전사 조절 캐스케이드를 반영합니다.

4. RNA-단백질 관계의 유전적 조절 (coExQTL):

   • 저자는 유전자가 RNA 발현과 사이토카인 분비 간의 상관관계를 수정하는 239 개의 유의한 coExQTL 을 식별했습니다.
   • 지질 대사 풍부화: 사이토카인 분비와 연결된 coExQTL 을 가진 유전자들은 지질 대사 과정에 유의하게 풍부했습니다.
     • PLTP (rs7270354): IFN$\gamma$ 후 6 가지 사이토카인 (IL-6, IL-10 포함) 의 분비와 PLTP RNA 를 연결하는 coExQTL.
     • VLDLR (rs6145): LPS 후 5 가지 사이토카인 (TNF, IL-1$\beta$ 포함) 과 VLDLR RNA 를 연결하는 coExQTL.
   • OAS1 아이소형 특이성: COVID-19 위험 로커스 rs10774671(OAS1의 스플라이스 변이) 은 coExQTL 로 작용하여 LPS 자극 후 OAS1 ENST00000452357 아이소형의 발현을 10 가지 사이토카인의 분비와 분리시켰습니다. 이 효과는 아이소형 특이적이며 RNA-RNA 수준에서는 관찰되지 않아, 이 변이가 단핵구 활성화에서 전사 후 조절 역할을 함을 강조합니다.

의의 및 주장
저자들은 본 연구가 맥락 특이적 방식으로 사이토카인 분비의 유전적 및 전사적 결정인자에 대한 포괄적인 규명을 제공한다고 주장합니다. 주요 주장은 다음과 같습니다:

• 맥락 특이성: 활성화된 단핵구에서의 사이토카인 분비 유전적 결정인자는 안정 상태 전혈 pQTL 연구에서 관찰된 것과 자주 구별되며, 질병 기전을 이해하기 위해 활성화된 상태에서 면역 세포를 연구할 필요성을 강조합니다.
• 전사 후 조절: 데이터는 RNA 발현이 종종 단백질 분비의 poor 예측인자임을 보여주며, 번역, 운송, 분해가 면역 반응에서 중요한 조절 층임을 나타냅니다.
• 질병 관련성: 식별된 로커스 (예: PDGFB, CCR5-Δ32, OAS1) 는 자가면역 질환 (궤양성 대장염, 일차성 담관 경화증) 및 중증 COVID-19 와 공위치하거나 알려진 위험 인자이며, 단핵구 사이토카인 생산의 조절 불균형이 유전적 위험과 질병 병리 사이의 기전적 연결고리임을 시사합니다.
• 지질 대사: 공발현 네트워크에서 지질 대사 유전자의 풍부화는 지단백 대사가 선천성 면역 사이토카인 반응을 조절하는 데 중추적인 역할을 함을 시사합니다.

저자들은 분석된 사이토카인의 수가 적고, 스플라이싱 분석을 제한하는 마이크로어레이 사용, 그리고 건강한 유럽계 기원 성인에 국한된 점 등 한계를 인정하며, 향후 연구는 이러한 발견을 다른 인구 집단 및 질병 상태로 확장해야 한다고 언급했습니다.