Triagem de pacientes e características da coorte
Um total de 230 pacientes idosos hospitalizados com exacerbação aguda da doença pulmonar obstrutiva crônica (EADPOC) preencheram os critérios de elegibilidade e foram incluídos na análise final. De acordo com a avaliação por ecocardiografia transtorácica, 120 pacientes foram alocados no grupo sem hipertensão pulmonar (sem HP), e 110 pacientes foram alocados no grupo com hipertensão pulmonar (HP) (Figura 1). As características gerais basais foram amplamente comparáveis entre os grupos. Não foram observadas diferenças estatisticamente significativas na distribuição por sexo, idade, histórico de tabagismo, histórico de consumo de álcool, alergia a medicamentos, hipertensão ou diabetes mellitus (todos P > 0,05; Tabela 1). A doença coronariana foi mais frequente no grupo com HP do que no grupo sem HP (13,64% vs. 4,17%, P = 0,017), assim como asma (29,09% vs. 15,00%, P = 0,011), infecção pulmonar (78,18% vs. 1,67%, P < 0,001) e insuficiência cardíaca (65,45% vs. 15,00%, P < 0,001; Tabela 1).
A qualidade ecocardiográfica foi revisada antes da atribuição ao grupo. Entre os pacientes incluídos no final, o sinal Doppler da regurgitação tricúspide foi mensurável em 198 pacientes (86,1%). Nos 32 pacientes restantes (13,9%), a classificação de HP foi baseada em achados estruturais ou funcionais interpretables do coração direito. Vinte e um registros foram revisados por um segundo médico ecocardiografista devido a janelas acústicas limitadas ou sinais de regurgitação tricúspide incompletos, e seis registros permaneceram ininterpretables e foram excluídos antes da construção final da coorte.
| Variável | Grupo sem HP | Grupo com HP (n = 110) | Z/t/χ² | Valor de p |
| (n = 120) |
| Masculino, n (%) | 83 (69,17) | 64 (58,18) | 3,003 | 0,099 |
| Idade, anos | 74,54 ± 7,92 | 74,68 ± 7,17 | -0,144 | 0,886 |
| Histórico de tabagismo, n (%) | 55 (45,83) | 47 (42,73) | 0,224 | 0,691 |
| Histórico de consumo de álcool, n (%) | 46 (38,33) | 37 (33,64) | 0,549 | 0,494 |
| Histórico de alergia a medicamentos, n (%) | 6 (5,00) | 1 (0,91) | 3,255 | 0,122 |
| Hipertensão, n (%) | 42 (35,00) | 37 (33,64) | 0,047 | 0,89 |
| Diabetes mellitus, n (%) | 12 (10,00) | 21 (19,09) | 3,86 | 0,06 |
| Doença cardíaca coronária, n (%) | 5 (4,17) | 15 (13,64) | 6,482 | 0,017 |
| Asma, n (%) | 18 (15,00) | 32 (29,09) | 6,698 | 0,011 |
| Bronquiectasia, n (%) | 14 (11,67) | 22 (20,00) | 3,019 | 0,102 |
| Infecção pulmonar, n (%) | 2 (1,67) | 86 (78,18) | 142,241 | <0,001 |
| Enfisema, n (%) | 3 (2,50) | 0 (0,00) | 2,786 | 0,248 |
| Falência respiratória, n (%) | 33 (27,50) | 41 (37,27) | 2,512 | 0,122 |
| Falência cardíaca, n (%) | 18 (15,00) | 72 (65,45) | 61,338 | <0,001 |
| Albumina, g/L | 36,38 ± 4,41 | 35,77 ± 4,09 | 1,084 | 0,279 |
| Creatinina, µmol/L | 68,75 (58,30, 83,05) | 66,75 (52,28, 82,03) | -0,951 | 0,341 |
| hs-CRP, mg/L | 8,15 (2,15, 48,58) | 12,60 (4,57, 90,03) | -2,11 | 0,035 |
| Dímero D, µg/mL | 0,415 (0,27, 0,68) | 0,61 (0,37, 1,41) | -3,212 | 0,001 |
| Fibrinogênio, g/L | 3,97 (2,90, 4,88) | 4,04 (3,09, 5,04) | -0,94 | 0,347 |
| Contagem de neutrófilos, ×109/L | 5,51 (3,67, 7,58) | 7,31 (4,78, 12,61) | -3,533 | <0,001 |
| Contagem de linfócitos, ×109/L | 1,12 (0,79, 1,59) | 0,85 (0,59, 1,25) | -3,252 | 0,001 |
| RDW-CV, % | 13,60 (13,00, 14,30) | 14,10 (13,30, 14,92) | -2,637 | 0,008 |
| Contagem de plaquetas, ×109/L | 200,00 (163,00, 252,00) | 180,50 (143,00, 209,00) | -2,783 | 0,005 |
| Volume médio de plaquetas, fL | 10,00 (9,10, 11,47) | 10,75 (9,70, 11,72) | -2,793 | 0,005 |
| RNL | 4,49 (2,76, 9,30) | 8,54 (4,73, 15,57) | -4,357 | <0,001 |
| BNP, pg/mL | 59,75 (41,45, 89,40) | 135,50 (64,65, 632,00) | -6,5 | <0,001 |
Tabela 1: Comparação das características clínicas entre os grupos sem HP e com HP. Esta tabela resume as características basais, comorbidades, marcadores inflamatórios, indicadores de coagulação, índices hematológicos e marcadores de estresse cardíaco em pacientes idosos com EPOC agudizada. As variáveis contínuas são apresentadas como média ± desvio padrão ou mediana com intervalo interquartil; as variáveis categóricas são apresentadas como número e porcentagem. Abreviações: AECOPD, exacerbação aguda da doença pulmonar obstrutiva crônica; PH, hipertensão pulmonar; hs-CRP, proteína C-reativa de alta sensibilidade; RDW-CV, coeficiente de variação da amplitude de distribuição eritrocitária; NLR, razão neutrófilo-linfócito; BNP, peptídeo natriurético tipo B.
Diferenças laboratoriais entre os grupos sem PH e com PH
O grupo com HP apresentou maior carga inflamatória, hematológica e de estresse cardíaco relacionada à coagulação (Tabela 1). A mediana da proteína C-reativa de alta sensibilidade foi maior no grupo com HP do que no grupo sem HP (12,60 mg/L vs. 8,15 mg/L, P = 0,035). A RNL também foi maior no grupo com HP (8,54 vs. 4,49, P < 0,001), com uma contagem de neutrófilos mais elevada (7,31 × 109/L vs. 5,51 × 109/L, P < 0,001) e uma contagem de linfócitos mais baixa (0,85 × 109/L vs. 1,12 × 109/L, P = 0,001). O D-dímero foi maior no grupo com HP (0,61 µg/mL vs. 0,415 µg/mL, P = 0,001). A contagem de plaquetas foi menor no grupo com HP (180,50 × 109/L vs. 200,00 × 109/L, P = 0,005), enquanto o volume médio de plaquetas (10,75 fL vs. 10,00 fL, P = 0,005) e o coeficiente de variação da largura de distribuição de hemácias (14,10% vs. 13,60%, P = 0,008) foram maiores. O BNP também diferiu claramente entre os grupos, com um valor mediano de 135,50 pg/mL no grupo com HP e de 59,75 pg/mL no grupo sem HP (P < 0,001). A albumina, a creatinina e o fibrinogênio não diferiram significativamente entre os grupos (todos P > 0,05; Tabela 1).
Seleção do preditor e regressão logística multivariável
Treze variáveis apresentaram diferenças significativas entre os grupos na análise univariada e foram analisadas conjuntamente com preditores clinicamente relevantes. A contagem de neutrófilos e a contagem de linfócitos não foram incluídas junto com a RNL, pois são componentes da razão. A infecção pulmonar e a insuficiência cardíaca foram avaliadas durante a avaliação de redundância devido à sua sobreposição com marcadores inflamatórios e BNP, respectivamente. Nove variáveis foram incluídas no modelo primário de regressão logística multivariada: proteína C-reativa de alta sensibilidade, D-dímero, coeficiente de variação da amplitude de distribuição de hemácias, contagem de plaquetas, volume médio de plaquetas, RNL, BNP, doença coronariana e histórico de asma. Todos os fatores de inflação da variância foram inferiores a 5, indicando ausência de multicolinearidade grave; o valor mais alto foi 2,16 para BNP (Tabela 2).
No modelo primário, quatro variáveis permaneceram independentemente associadas à HP em pacientes idosos com DPOC agudizada (Tabela 2). O dímero D esteve associado a maiores chances de HP (OR = 1,251, IC 95%: 1,001–1,563, P = 0,049), assim como a RNL (OR = 1,043, IC 95%: 1,003–1,083, P = 0,033), o BNP (OR = 1,003, IC 95%: 1,001–1,005, P < 0,001) e o histórico de asma (OR = 3,054, IC 95%: 1,489–6,267, P = 0,002). A proteína C-reativa de alta sensibilidade, o coeficiente de variação da amplitude de distribuição eritrocitária, a contagem de plaquetas, o volume médio de plaquetas e a doença arterial coronariana não se mostraram independentemente associados à HP após ajuste (todos P > 0,05; Tabela 2).
| Variável | B | EP | Wald χ2 | OR (IC 95%) | Valor-p | FIV |
| hs-CRP | -0.002 | 0.002 | 1.099 | 0.998 (0.995, 1.002) | 0.295 | 1.22 |
| D-dímero | 0.224 | 0.114 | 3.877 | 1.251 (1.001, 1.563) | 0.049 | 1.16 |
| RDW-CV | 0.088 | 0.124 | 0.504 | 1.092 (0.856, 1.393) | 0.478 | 1.31 |
| Contagem de plaquetas | -0.003 | 0.002 | 1.281 | 0.997 (0.993, 1.002) | 0.258 | 1.44 |
| Volume médio de plaquetas | 0.096 | 0.088 | 1.197 | 1.101 (0.927, 1.309) | 0.274 | 1.37 |
| RNL | 0.042 | 0.02 | 4.545 | 1.043 (1.003, 1.083) | 0.033 | 1.53 |
| BNP | 0.003 | 0.001 | 13.58 | 1.003 (1.001, 1.005) | <0.001 | 2.16 |
| Doença cardíaca coronariana | 0.662 | 0.637 | 1.078 | 1.938 (0.556, 6.759) | 0.299 | 1.18 |
| Histórico de asma | 1.117 | 0.367 | 9.271 | 3.054 (1.489, 6.267) | 0.002 | 1.11 |
| Constante | -3.285 | 2.13 | 2.378 | — | 0.123 | — |
Tabela 2: Modelo primário de regressão logística multivariável para hipertensão pulmonar. Esta tabela apresenta os coeficientes de regressão, erros padrão, qui-quadrado de Wald χ2 valores, razões de chances, intervalos de confiança de 95%, valores de P e fatores de inflação de variância. D-dímero, RNL, BNP e histórico de asma permaneceram independentemente associados à HP. Abreviações: PH, hipertensão pulmonar; SE, erro padrão; OR, razão de chances; CI, intervalo de confiança; VIF, fator de inflação da variância; NLR, razão neutrófilo-linfócito; BNP, peptídeo natriurético tipo B.
Análise de sensibilidade incluindo infecção pulmonar e insuficiência cardíaca
Um modelo de sensibilidade foi construído adicionando infecção pulmonar e insuficiência cardíaca ao conjunto de preditores primários (Tabela 3). A infecção pulmonar manteve-se fortemente associada à HP (OR = 24,80, IC 95%: 5,62–109,44, P < 0,001), e a insuficiência cardíaca também esteve associada à HP (OR = 3,36, IC 95%: 1,52–7,43, P = 0,003). Após a inclusão dessas duas variáveis, as associações de BNP e NLR foram atenuadas, mas permaneceram com direção consistente. O dímero D e o histórico de asma também mantiveram associações positivas. O modelo de sensibilidade apresentou maior discriminação aparente do que o modelo primário (AUC = 0,91, IC 95%: 0,87–0,95), mas introduziu sobreposição substancial com estados diagnósticos inflamatórios e cardíacos. Portanto, o modelo com quatro variáveis foi mantido como o nomograma primário, enquanto o modelo de sensibilidade foi relatado para ilustrar a influência da infecção pulmonar e da insuficiência cardíaca na interpretação do modelo.
Construção do nomograma
Um nomograma foi construído utilizando os quatro preditores independentes mantidos no modelo primário de regressão logística multivariável: D-dímero, LNR, BNP e histórico de asma (Figura 2A). Cada preditor contribuiu com um valor em pontos de acordo com seu peso de regressão, e a pontuação total foi associada à probabilidade predita individualizada de HP. O caso ilustrativo foi revisado para evitar valores muito distantes da distribuição clínica observada. No exemplo revisado, um paciente idoso hospitalizado com DPAECO apresentava nível de D-dímero de 0,90 µg/mL, LNR de 10,0, nível de BNP de 230 pg/mL e histórico de asma. A pontuação total estimada foi de 173 pontos, correspondendo a uma probabilidade de HP de 78,4% no nomograma (Figura 2B).

Figura 2: Nômograma para estimar o risco de hipertensão pulmonar em pacientes idosos com EPOC agudizada. (A) O nômograma foi construído utilizando D-dímero, LNR, BNP e histórico de asma. Cada preditor recebeu um valor em pontos, e a pontuação total foi associada à probabilidade predita de HP. (B) Um exemplo revisado mostrou D-dímero = 0.90 µg/mL, LNR = 10.0, BNP = 230 pg/mL e histórico de asma = Sim, correspondendo a 173 pontos e uma probabilidade estimada de HP de 78,4%. Abreviaturas: AECOPD, exacerbação aguda da doença pulmonar obstrutiva crônica; PH, hipertensão pulmonar; NLR, razão neutrófilo-linfócito; BNP, peptídeo natriurético tipo B. Clique aqui para visualizar uma versão maior desta figura.
Discriminação do modelo de predição
O nomograma de quatro variáveis alcançou uma área sob a curva característica de operação do receptor de 0,82, com um intervalo de confiança de 95% de 0,760–0,874 (Figura 3A). Em comparação, o D-dímero isolado apresentou uma área sob a curva de 0,623, a LCR isolada apresentou uma área sob a curva de 0,666 e o BNP isolado apresentou uma área sob a curva de 0,748. Um modelo que combina os três marcadores quantitativos resultou em uma área sob a curva de 0,793 (Figura 3B). Assim, o nomograma demonstrou melhor discriminação do que cada marcador individual e discriminação ligeiramente superior à do modelo com três marcadores.

Figura 3: Curvas características de operação do receptor para o nomograma e preditores quantitativos. (A) O nomograma de quatro variáveis alcançou uma AUC de 0,82, com um IC de 95% de 0,760–0,874. (B) Os valores de AUC foram 0,623 para D-dímero, 0,666 para RNL, 0,748 para BNP e 0,793 para o modelo combinado de marcadores quantitativos. Abreviações: ROC, característica de operação do receptor; AUC, área sob a curva; RNL, razão neutrófilo-linfócito; BNP, peptídeo natriurético tipo B. Clique aqui para visualizar uma versão maior desta figura.
Calibração e validação interna
A validação interna foi realizada utilizando reamostragem por bootstrap com 1.000 réplicas. A análise de calibração comparou as probabilidades preditas de HP com os desfechos observados de HP (Figura 4). O intercepto de calibração aparente foi 0,0000 e a inclinação de calibração aparente foi 1,0000. Após correção por bootstrap, o intercepto de calibração foi -0,08 e a inclinação de calibração foi 0,86. O erro médio absoluto de calibração foi 0,054 e o escore de Brier foi 0,1878. O teste de Hosmer-Lemeshow foi estatisticamente significativo (P = 0,0005), indicando desvio residual de calibração. Portanto, apesar da boa discriminação, os resultados de calibração foram interpretados com cautela, e uma recalibração pode ser necessária antes da utilização em outros contextos.

Figura 4: Curva de calibração do nomograma de predição de hipertensão pulmonar. A curva de calibração compara as probabilidades preditas e observadas de HP. As curvas aparente e corrigida por bootstrap foram geradas utilizando 1.000 reamostragens por bootstrap. O escore de Brier foi de 0,1878 e o teste de Hosmer-Lemeshow foi significativo (P = 0,0005), indicando desvio residual de calibração. Abreviações: HP, hipertensão pulmonar; HL, Hosmer-Lemeshow. Clique aqui para visualizar uma versão maior desta figura.
Utilidade clínica avaliada pela análise de curva de decisão
A análise da curva de decisão mostrou que a curva do nomograma permaneceu acima das estratégias de referência de tratar todos e não tratar nenhum, em uma faixa de probabilidade de limiar de aproximadamente 1%–70% (Figura 5). Esse achado sugeriu um benefício líquido potencial na identificação de pacientes idosos com DPOC exacerbada que possam necessitar de revisão ecocardiográfica mais rigorosa, monitoramento cardiopulmonar ou avaliação diagnóstica adicional. Como o modelo foi validado apenas internamente, os resultados da curva de decisão foram interpretados como evidência preliminar de utilidade clínica, e não como comprovação de benefício na implementação.

Figura 5: Análise de curva de decisão do nomograma de predição de hipertensão pulmonar. A análise de curva de decisão mostrou que o nomograma proporcionou benefício líquido potencial ao longo de probabilidades de limiar variando aproximadamente de 1% a 70%, em comparação com as estratégias de tratar todos e não tratar nenhum. Abreviações: DCA, análise de curva de decisão; PH, hipertensão pulmonar. Clique aqui para visualizar uma versão maior desta figura.
Resultados gerais do modelo
O risco de HP em pacientes idosos hospitalizados com DPOC exacerbação aguda esteve associado a informações relacionadas à inflamação, coagulação, estresse cardíaco e doenças das vias aéreas. NLR, BNP, dímero D e histórico de asma associaram-se independentemente à HP e foram integrados a um nomograma com boa discriminação. A análise de sensibilidade mostrou que infecção pulmonar e insuficiência cardíaca estiveram fortemente associadas à HP, mas sua inclusão aumentou a sobreposição com marcadores inflamatórios e de estresse cardíaco. A calibração não foi totalmente ideal, especialmente considerando o teste de Hosmer-Lemeshow significativo. Esses achados sustentam o nomograma como uma ferramenta de estratificação de risco validada internamente, que necessita de validação externa e possível recalibração antes de uso clínico mais amplo.
DISPONIBILIDADE DE DADOS:
O conjunto de dados que sustenta as descobertas deste estudo foi depositado no repositório Figshare e está disponível publicamente em: https://doi.org/10.6084/m9.figshare.32765667.v1 .
| Item de Análise | Nomograma Primário com Quatro Variáveis | Modelo de Sensibilidade Incluindo Infecção Pulmonar e Insuficiência Cardíaca |
| Previsores incluídos | D-dímero, NLR, BNP, histórico de asma | D-dímero, NLR, BNP, histórico de asma, infecção pulmonar, insuficiência cardíaca |
| OR de infecção pulmonar (IC 95%) | Não incluído | 24,80 (5,62; 109,44) |
| OR de insuficiência cardíaca (IC 95%) | Não incluído | 3,36 (1,52; 7,43) |
| Direção do D-dímero | Positiva | Positiva |
| Direção do NLR | Positiva | Positiva, mas atenuada |
| Direção do BNP | Positiva | Positiva, mas atenuada |
| Direção do histórico de asma | Positiva | Positiva |
| AUC (IC 95%) | 0,82 (0,760; 0,874) | 0,91 (0,87; 0,95) |
| Pontuação de Brier | 0,1878 | 0,146 |
| Intercepto aparente de calibração | 0 | 0 |
| Inclinação aparente de calibração | 1 | 1 |
| Intercepto de calibração corrigido por bootstrap | -0,08 | -0,05 |
| Inclinação de calibração corrigida por bootstrap | 0,86 | 0,89 |
| Erro médio absoluto de calibração | 0,054 | 0,041 |
| Valor-p de Hosmer-Lemeshow | 0,0005 | 0,021 |
| Intervalo de limiar da curva de decisão com benefício líquido | Aproximadamente 1%–70% | Aproximadamente 2%–75% |
| Função do modelo | Nomograma primário | Apenas análise de sensibilidade |
Tabela 3: Análise de sensibilidade e indicadores de validação interna. Esta tabela compara o nomograma primário de quatro variáveis com o modelo de sensibilidade que incluiu adicionalmente infecção pulmonar e insuficiência cardíaca. São apresentados discriminação, calibração, teste de Hosmer-Lemeshow, amplitude da curva de decisão e indicadores corrigidos por reamostragem bootstrap. Abreviações: AUC, área sob a curva característica de operação do receptor; OR, razão de chances; CI, intervalo de confiança; NLR, razão neutrófilo-linfócito; BNP, peptídeo natriurético tipo B.