전 세계적으로 고령화가 진행되면서 단순히 오래 사는 것보다 기능과 삶의 질을 유지하는 ‘건강한 노화’가 중요한 사회적·정책적 목표가 되고 있다. 그런데 노년기의 우울과 불안은 인지 저하, 허약, 혈관질환, 사망 등 여러 부정적 건강 결과와 관련되는 것으로 알려져 있다. 문제는 이러한 정서적 증상 부담이 신체의 생물학적 노화와도 연결되는지, 특히 서구권이 아닌 아시아 노년층에서도 그런 관련이 나타나는지가 충분히 연구되지 않았다는 점이다. 이 논문은 바로 이 공백을 겨냥한다.

이 연구의 핵심 결과변수는 후성유전학적 나이 가속, 즉 EAA이다. 배경지식으로 설명하면, 후성유전학은 DNA 염기서열 자체를 바꾸지 않으면서 유전자 발현과 관련된 조절이 일어나는 방식을 가리킨다. 대표적인 메커니즘이 DNA 메틸화이다. 나이가 들수록 DNA 메틸화 패턴도 변화하기 때문에, 여러 CpG 부위의 메틸화 정보를 이용해 개인의 생물학적 나이를 추정하는 ‘후성유전학적 시계’가 개발되었다. 이 논문은 Horvath, Hannum, PhenoAge, GrimAge, DunedinPACE 등 여러 세대의 시계를 활용했다. Horvath, Hannum, PhenoAge, GrimAge 계열 지표는 주성분 기반으로 계산된 뒤 생활연령의 영향을 제거한 EAA로 분석되었고, DunedinPACE는 노화 ‘속도’에 가까운 지표로 다루어졌다. EAA는 후성유전학적 나이를 실제 생활연령에 회귀시킨 뒤 남는 잔차로 계산된다. 쉽게 말해, 같은 나이의 다른 사람들과 비교했을 때 후성유전학적으로 더 늙어 보이면 EAA가 양수이고, 더 젊어 보이면 음수이다.

연구진의 주요 질문은 두 가지로 정리할 수 있다. 첫째, 노년기의 우울 증상과 불안 증상은 여러 후성유전학적 노화 지표에서 EAA와 관련되는가? 둘째, 시간이 지나면서 한 개인의 우울 또는 불안 증상이 변할 때, 그 사람의 EAA도 함께 변하는가? 연구진은 특히 우울 증상과 2세대 후성유전학적 시계인 PhenoAge 기반 지표, 즉 PCPhenoEAA의 관련을 주된 분석으로 설정했다. 다른 후성유전학적 시계나 불안 관련 분석은 탐색적 분석으로 다루었다.

분석 대상은 싱가포르의 전향적 코호트 연구인 DaHA에서 DNA 메틸화 자료가 있는 지역사회 거주 노년층이다. 원래 DaHA에는 2011년부터 2017년 사이 60세 이상 싱가포르 시민 또는 영주권자 1010명이 등록되었고, 이 중 620명이 2017년부터 2020년 사이 재조사를 받았다. 본 연구의 최종 분석 표본은 DNA 메틸화 자료가 있는 672명이다. 평균 연령은 70세였고 표준편차는 5.8세였다. 전체 표본 중 484명은 여성, 188명은 남성이었다. 결혼 상태인 참여자는 68.2%, 4~5룸 HDB에 해당하는 중간 주거층 참여자는 72.9%, 정규교육을 받지 않은 참여자는 35.6%였다. 평균 교육연수는 6.5년이었고, BMI 평균은 24, 만성질환 수 평균은 2.7개였다. 흡연과 음주는 대체로 낮은 수준이어서 정기 흡연자는 3.9%, 정기 음주자는 10.4%였다. 신체활동은 매일 한다고 응답한 사람이 37.6%로 가장 많았고, 사회활동은 주 1~3회가 43.2%로 가장 많았다. 전반적으로 우울과 불안 증상 수준은 낮았다. GDS-15 평균은 1.30, GAI-20 평균도 1.30이었다. 우울 선별 기준인 GDS 5점 이상에 해당한 사람은 56명, 8.3%였고, 불안 선별 기준인 GAI 8점 이상에 해당한 사람은 28명, 4.2%였다. 남녀 구성에는 차이가 있었다. 남성은 평균 교육연수가 더 길고 결혼 상태 비율이 높았으며, 여성은 무학력과 사별 비율이 더 높았다. 흡연과 음주도 남성에서 더 흔했다. 이런 표본 특성은 성별, 교육, 결혼상태, 생활습관을 보정해야 하는 이유를 잘 보여준다.

우울 증상은 15문항 노인우울척도인 GDS-15로 측정되었다. 총점은 0점에서 15점이며, 점수가 높을수록 우울 증상이 크다. 연구진은 GDS 5점 이상을 임상적으로 의미 있는 우울 선별 기준으로도 사용했다. 불안 증상은 20문항 노인불안척도인 GAI-20으로 측정되었고, 총점은 0점에서 20점이다. GAI 8점 이상은 불안장애 가능성을 시사하는 기준으로 사용되었다. DNA 메틸화는 금식 후 정맥혈에서 채취한 혈액을 이용해 측정되었고, Illumina 850K array를 사용했다. 연구진은 품질관리 절차를 거쳐 여러 후성유전학적 시계를 계산했으며, 기술적 잡음을 줄이기 위해 HorvathAge, Horvath2Age, HannumAge, PhenoAge, GrimAge에는 주성분 기반 나이 지표를 적용했다. 기술통계에서 EAA 평균이 전체 표본에서 0에 가깝게 나타나는 것은 생활연령을 제거한 잔차 지표라는 계산 방식과 관련된다. 다만 회귀분석의 β는 주로 표준화된 계수로 제시되므로, 본문에서 말하는 0.087이나 0.244 같은 값은 ‘년 단위’가 아니라 표준편차 단위의 차이로 이해해야 한다.

주요 분석에는 선형혼합효과모형이 사용되었다. 일부 참여자는 두 번의 자료를 제공했기 때문에, 같은 사람에게서 나온 반복 관측치들이 서로 독립이라고 가정하면 안 된다. 이를 고려하기 위해 연구진은 무작위 절편을 포함해 개인 내 군집성을 조정했다. 분석의 참여자 기준 표본 수는 672명이었다. 모형은 단계적으로 구성되었다. 보정하지 않은 모형, 인구사회학적·건강 관련 요인을 보정한 모형, 그리고 생활습관 변수까지 포함한 완전 보정 모형을 비교했다. 보정 변수에는 성별, 교육연수, 민족, 결혼상태, 주거 유형, BMI, 만성질환 수, 신체활동과 사회활동 빈도, 흡연, 음주 등이 포함되었다. 생활연령은 EAA 계산 과정에서 이미 분리되었기 때문에 추가 공변량으로 넣지 않았다. 기초 상관구조를 보면 생활연령은 여러 후성유전학적 시계 산출값 및 관련 생체표지와 뚜렷하게 함께 움직였고, 교육연수와 몇몇 주성분 변수도 여러 지표와 상관을 보였다. 반면 우울·불안 점수의 단순상관은 지표별로 방향과 강도가 일정하지 않았다. 따라서 단순상관만으로 해석하기보다, 생활연령을 제거한 EAA와 여러 공변량을 고려한 혼합모형을 사용한 점이 방법론적으로 중요하다.

연구진은 DNA 메틸화 자료가 두 시점 모두 있는 116명, 총 232개 관측치를 대상으로 개인 내 변화 분석도 수행했다. 이 분석에서는 추적조사 값에서 기초조사 값을 뺀 변화량을 계산했다. 예를 들어 우울 증상의 변화량과 PCPhenoEAA 변화량이 같은 개인 안에서 함께 움직이는지를 살펴본 것이다. 이 방식은 사람들 사이의 고정된 차이, 예를 들어 원래부터 건강 상태나 사회경제적 조건이 다른 차이만으로 결과가 설명되는지를 어느 정도 줄여 볼 수 있다는 장점이 있다. 그러나 여전히 관찰자료이므로 인과관계를 확정하는 분석은 아니다.

가장 중요한 결과는 우울 증상이 높을수록 PCPhenoEAA가 높다는 점이다. 세 단계의 보정수준을 나란히 보면, 우울 증상을 연속변수로 넣었을 때 PCPhenoEAA와의 양의 관련은 보정하지 않은 모형, 부분 보정 모형, 완전 보정 모형에서 모두 유지되었다. 완전 보정 모형에서 우울 증상이 1표준편차 증가할 때 PCPhenoEAA는 0.087표준편차 증가했다. 추정치는 β = 0.087, 95% 신뢰구간은 [0.023, 0.151], p = 0.008이었다. 이는 효과 크기가 매우 큰 결과라기보다는, 통계적으로 유의하고 일관된 작은 규모의 관련이라고 해석하는 것이 적절하다. 또한 GDS 5점 이상으로 우울 선별 기준을 충족한 참여자는 그렇지 않은 참여자에 비해 평균적으로 PCPhenoEAA가 0.244표준편차 높았다. 이 결과의 추정치는 β = 0.244, 95% 신뢰구간 [0.027, 0.461], p = 0.030이었다. 즉 우울 점수를 연속적인 증상 부담으로 보아도, 임상적으로 의미 있는 선별 기준으로 나누어 보아도 핵심 패턴은 같은 방향이었다.

우울 선별 기준을 적용했을 때 DunedinPACE에서도 노화 속도가 더 빠른 방향의 관련이 관찰되었다. 다만 이 관련은 다중검정 보정 후에는 유의하지 않은 결과로 해석되었다. 나머지 후성유전학적 시계들에서는 우울 증상과의 관련이 일관되게 유의하지 않았다. 일부 1세대 시계나 다른 지표에서 추정치가 양의 방향을 보이거나 명목상 관련이 관찰되더라도, 신뢰구간이 0을 포함하거나 다중검정 보정을 통과하지 못해 확정적 결과로 보기는 어렵다. 따라서 이 연구의 핵심 결과는 여러 지표 전반에 걸친 일반적 노화 가속이라기보다, PCPhenoEAA에서 가장 뚜렷하게 나타난 우울 관련 신호라고 보는 편이 정확하다. 반면 불안 증상은 분석된 후성유전학적 노화 지표 전반에서 일관된 유의한 관련을 보이지 않았다. 불안 점수를 연속변수로 보아도, GAI 8점 이상 선별 기준으로 보아도 대부분의 추정치는 0에 가까웠고 신뢰구간도 넓게 0을 포함했다. 연구진은 표본 내 불안 증상의 변동성이 작았기 때문에 관련을 탐지할 통계적 여지가 제한되었을 가능성도 언급한다.

반복 DNA 메틸화 자료가 있는 116명에 대한 종단 분석에서도 흥미로운 결과가 나왔다. 전체적으로는 추적 시점에서 우울과 불안 증상이 감소하는 경향이 있었고, EAA는 개인 내에서 비교적 안정적이었다. 그럼에도 개인 내 우울 증상 변화는 PCPhenoEAA 변화와 함께 움직였다. 시간에 따라 우울 증상이 1표준편차 증가할 때 PCPhenoEAA는 0.19표준편차 증가했다. 추정치는 β = 0.19, 95% 신뢰구간 [0.01, 0.36], p = 0.043이었다. 혈액 세포 구성 변화까지 추가로 보정한 민감도 분석에서도 이 관련은 유지되었다. 이는 우울-EAA 관련이 단순히 우울한 사람과 그렇지 않은 사람 사이의 차이만은 아닐 가능성을 시사한다.

연구진은 결과의 강건성을 확인하기 위해 여러 민감도 분석을 수행했다. 혈액 세포 구성 보정, 인지 상태 추가 보정, 이상치 영향을 줄이는 강건회귀, 비선형 관계 점검 등이 포함되었다. 우울 증상과 PCPhenoEAA의 관련은 이러한 분석에서도 대체로 유지되었다. 모형 진단에서도 선형성·정규성 가정에 큰 문제가 없고, 뚜렷한 이분산성 증거가 없다는 점이 확인되었다. 이상치의 영향을 줄이는 강건회귀에서도 대부분의 관측치가 원래 분석과 거의 같은 가중치를 받았고, 추정 패턴도 원래 결과와 대체로 비슷했다. 혈액 세포 구성 관련 주성분이나 인지 상태를 추가로 넣어도 우울-PCPhenoEAA 관련은 방향이 유지되고 통계적으로도 유의했다. 성별과의 상호작용도 탐색적으로 검토되었는데, 일부 분석은 우울과 후성유전학적 노화의 관련이 성별에 따라 다를 수 있음을 시사했다. 논문은 이러한 가능성을 해석할 때 CRF 수용체 결합 방식이나 DNA 메틸화 효소 발현 차이처럼 성별에 따라 다르게 나타날 수 있는 생물학적 경로를 함께 언급한다. 다만 이러한 상호작용 결과는 모든 지표와 모든 우울 변수화 방식에서 일관되지는 않았고, 해당 생물학적 경로들이 DaHA 자료에서 직접 측정되어 매개분석으로 검증된 것도 아니다. 따라서 성별 차이는 결론이라기보다 후속 연구에서 검증할 가설로 보는 것이 적절하다.

이 논문은 생물학적 메커니즘을 직접 검증하지는 않았다. 다만 연구진은 기존 문헌을 바탕으로 우울 증상이 스트레스 반응, 면역·염증, 대사 조절, 보상회로, 누적 생리적 부담과 관련될 수 있다고 설명한다. 가능한 경로는 크게 몇 가지로 정리할 수 있다. 첫째, 우울과 만성 스트레스는 HPA 축 조절 이상, 코르티솔 반응, 글루코코르티코이드 관련 조절, 산화스트레스 및 염증성 사이토카인 변화와 연결될 수 있다. 둘째, 도파민성 기능 이상은 보상 민감도 저하, 운동 기능이나 보행 속도 저하 같은 노화 관련 기능 변화와 함께 논의될 수 있다. 셋째, 반복되는 스트레스와 신체 부담이 누적되면 allostatic load가 높아지고, IL-6, TNF, CRP 같은 염증 표지, 허약, 인지 저하와 연결될 수 있다는 설명도 제시된다. 넷째, 성별에 따른 CRF 관련 신호전달이나 DNA 메틸화 효소 발현 차이도 가능한 조절 요인으로 논의된다. 중요한 점은 이 경로들이 이번 분석에서 직접 입증된 인과 메커니즘이 아니라, 우울 증상과 PCPhenoEAA의 통계적 관련을 해석하기 위한 생물학적으로 그럴듯한 틀이라는 점이다. PhenoAge 같은 2세대 후성유전학적 시계는 단순히 생활연령을 예측하는 데 그치지 않고 건강수명이나 사망 위험과 관련된 생의학적 표지를 반영하도록 훈련된 지표이다. 따라서 앞서 확인된 것처럼 우울 관련 신호가 모든 시계에서 똑같이 나타나기보다 PCPhenoEAA에서 가장 뚜렷했던 것은, 우울이 단순한 시간 경과보다 염증, 생리적 부담, 기능 저하와 가까운 노화 측면에 더 민감하게 연결되었기 때문일 가능성이 있다.

이 연구의 의의는 노년기 정신건강과 생물학적 노화 사이의 연결을 아시아 지역사회 표본에서 검토했다는 데 있다. 특히 우울 증상이 PCPhenoEAA와 관련된다는 결과는 건강한 노화 정책에서 우울 선별과 개입을 더 중요하게 다룰 필요가 있음을 시사한다. 노년기 우울은 비교적 평가와 개입이 가능한 요인이다. 따라서 지역사회 노인 보건서비스, 1차 의료, 정신건강 프로그램에서 우울 증상을 조기에 확인하고 적절한 지원을 제공하는 것은 생물학적 노화와 관련된 건강위험을 낮추는 데 도움이 될 가능성이 있다. 또한 우울이 스트레스 조절, 염증, 허약, 인지 기능 같은 여러 건강 영역과 맞물릴 수 있다는 해석은 정신건강 개입을 신체건강 관리, 사회적 지지, 활동 유지와 분리해서 보지 말아야 함을 시사한다. 다만 여기서 말하는 것은 정책적 가능성이지, 이 연구 하나만으로 우울 치료가 EAA를 낮춘다고 입증되었다는 뜻은 아니다.

가장 중요한 한계는 관찰연구라는 점이다. 연구진은 다양한 공변량을 보정하고 개인 내 변화 분석까지 수행했지만, 측정되지 않은 교란요인이 남아 있을 수 있다. 예를 들어 염증 표지자, 수면 습관, 심리사회적 스트레스, 회복탄력성 같은 요인은 충분히 조정되지 않았을 수 있다. 또한 두 시점 모두 DNA 메틸화 자료가 있는 종단 분석 표본은 116명으로 제한되었고, 반복 혈액 채취에 참여하지 않은 사람들의 선택 과정이 결과에 영향을 주었을 가능성도 있다. 표본 자체도 비교적 건강하고 우울·불안 증상 수준이 낮았기 때문에, 임상적 우울증이나 불안장애 환자 집단에 그대로 일반화하기는 어렵다. 특히 불안 선별 기준을 충족한 사람이 4.2%에 그쳤다는 점은 불안과 EAA의 관련을 탐지하기 어려웠을 가능성을 함께 고려하게 한다. 더불어 논문에서 제시된 HPA 축, 도파민성 기능, allostatic load, 성별 특이적 신호전달이나 DNA 메틸화 조절 경로는 결과 해석을 돕는 개념적 설명이지, 본 연구가 모두 직접 측정해 검증한 매개경로는 아니다. 따라서 이 연구의 결론은 ‘우울 증상과 PCPhenoEAA가 관련된다’는 수준에 머물러야 하며, 그 생물학적 경로와 임상적 개입 가능성은 별도의 연구 설계로 확인해야 한다.