Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Med., 02 January 2026

Sec. Rheumatology

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1724495

Exploration of factors influencing pulmonary injury in patients with anti-neutrophil cytoplasmic antibody-associated vasculitis and construction of a predictive model

Linlin ZhengLinlin ZhengChongyang Liu
Chongyang Liu*
  • Department of Rheumatology and Immunology, The Third Affiliated Hospital of Chongqing Medical University, Chongqing, China

Background: A bad prognosis results from pulmonary involvement of anti-neutrophil cytoplasmic antibody-associated vasculitis (AAV). A method for early detection of AAV patients at risk of developing pulmonary involvement is still required, despite the existence of predictive models for death. The purpose of this study was to determine risk factors and develop a predictive model for pulmonary involvement in AAV.

Methods: Seventy-one treatment-naïve AAV patients (38 with pulmonary involvement) participated in a retrospective cross-sectional investigation. To choose predictors and create a nomogram, multivariate logistic regression and the least absolute shrinkage and selection operator (LASSO) were employed. Receiver operating characteristic (ROC) curves, calibration plots, and decision curve analysis (DCA) with a 7:3 training-to-validation split were used to assess the model.

Results: Four independent predictors were identified: intermittent symptom pattern (OR = 8.52, 95% CI: 2.17–33.46), limb numbness/paresthesia (OR = 12.21, 95% CI: 1.98–75.50), older age at onset (OR = 1.05 per year, 95% CI: 1.01–1.10), and lower hemoglobin level (OR = 0.97 per g/L, 95% CI: 0.95–0.99). A nomogram incorporating these factors demonstrated good discrimination, with an area under the curve (AUC) of 0.88 (95% CI: 0.77–0.98) in the training set and 0.83 (95% CI: 0.64–1.00) in the validation set. Calibration and decision curve analysis confirmed the model’s clinical utility.

Conclusion: Using standard clinical markers, we created and internally verified a pragmatic prediction model for pulmonary involvement in AAV. Personalized monitoring and early risk categorization may benefit from this approach. To verify its generalizability, external validation in prospective, multi-center cohorts is advised.

Introduction

Anti-neutrophil cytoplasmic antibody-associated vasculitis (AAV) is a systemic autoimmune disease that affects multiple organ systems. Its primary pathogenic mechanism involves the excessive activation of neutrophils by ANCA, leading to abnormal interactions between activated neutrophils and endothelial cells, ultimately causing vasculitic lesions. AAV encompasses three main subtypes: microscopic polyangiitis, granulomatosis with polyangiitis, and eosinophilic granulomatosis with polyangiitis (1). The hallmark features of AAV include systemic inflammatory reactions and fibrinoid necrosis in the walls of small blood vessels, accompanied by positive serum ANCA. However, unlike other vasculitides, immune complex deposition is typically absent in the affected vascular walls (2).

Over the past 30 years, studies have demonstrated a significant increase in the incidence of AAV (2). In 2008, the annual overall incidence of ANCA-associated vasculitis in Europe and North America was approximately 20 cases per million (3). By 2015, the overall incidence had risen to 420/125,000 (2). In China, a sample survey of hospitalized patients reported an incidence rate of approximately 0.025% for AAV (4). It has been reported that untreated granulomatosis with polyangiitis (GPA) has a one-year mortality rate as high as 80%, with an average survival time of about 5 months (5).

AAV predominantly affects middle-aged and elderly individuals and typically has an insidious onset. The disease often alternates between periods of remission and exacerbation, lacks specific clinical manifestations, and primarily involves the lungs, kidneys, and nervous system. In these patients, 85–90% of those with GPA will develop pulmonary involvement, which manifests as pulmonary hemorrhage or interstitial pulmonary fibrosis. In contrast, the incidence of pulmonary damage in microscopic polyangiitis (MPA) is slightly lower, ranging from 25 to 55% (6). Moreover, in patients with MPA, approximately 30–50% of those diagnosed develop pulmonary involvement within 1 year. The corresponding proportions are up to 60% in EGPA (95% including asthma) and 67 to 85% in GPA (7). Pulmonary injury can exhibit significant variability and progress rapidly, ranging from asymptomatic cases to severe pulmonary hemorrhage, dyspnea, and respiratory distress, potentially resulting in patient mortality. In patients with severe pulmonary injury, the 3-month mortality rate is 23%, while the 12-month mortality rate is 26% (8).

While the variability in pulmonary involvement and its associated high mortality rates highlight the clinical significance of early intervention, the lack of consensus on the risk factors contributing to this condition presents a major challenge. Understanding these factors is crucial for the development of predictive tools that could facilitate timely diagnosis and targeted treatment, ultimately improving patient outcomes. Therefore, we designed a retrospective study to explore the correlations between various clinical and pathological indicators and the occurrence of pulmonary involvement in AAV patients, and to establish a predictive model that can aid clinicians in making timely and informed medical decisions.

Materials and methods

Study design

This retrospective cross-sectional study was conducted using clinical data from patients diagnosed with AAV who were hospitalized at the Third Affiliated Hospital of Chongqing Medical University between August 2018 and October 2024. The study protocol was approved by the hospital’s Ethics Committee (Approval No.: 2025-Kelun-2). Informed consent was waived due to the retrospective nature of the study and the use of fully anonymized clinical data. Our final cohort included 71 eligible treatment-naïve AAV patients (38 with pulmonary involvement). Model development and reporting followed the TRIPOD guidelines (9).

Participant selection

A total of 71 AAV patients were included, comprising 38 with pulmonary involvement and 33 without. Diagnosis of AAV was based on the 2012 Chapel Hill Consensus Conference criteria (1, 10, 11), supported by serological ANCA positivity and/or histopathological evidence of necrotizing vasculitis. Both ANCA-positive and ANCA-negative patients with histologically confirmed renal vasculitis were eligible.

Inclusion criteria: (1) Age ≥ 18 years; (2) New diagnosis or diagnosis within 1 year of AAV, with no prior treatment with glucocorticoids, immunosuppressants, or biologics; (3) Complete clinical records available.

Exclusion criteria: (1) Age < 18 years; (2) Secondary vasculitis due to drugs, malignancies, or infections (e.g., tuberculosis); (3) Coexisting autoimmune diseases (e.g., systemic lupus erythematosus, rheumatoid arthritis, anti-GBM nephritis); (4) Pre-existing pulmonary diseases (e.g., tumors) or history of pulmonary injury due to infection, drugs, or environmental factors prior to AAV onset.

Diagnostic criteria for pulmonary involvement

Pulmonary involvement was defined according to established composite criteria (7, 12, 13): (a) Clinical symptoms (e.g., hemoptysis, cough, dyspnea); (b) Imaging findings consistent with alveolar hemorrhage, interstitial changes, nodules, cavitation, or tracheobronchial inflammation; (c) Pulmonary function impairment, defined as forced expiratory volume in 1 s (FEV₁) < lower limit of normal (LLN) or < 70% predicted, FEV₁/forced vital capacity (FVC) < LLN, or diffusing capacity of the lung for carbon monoxide (DLCO) < 80% predicted or below LLN (14). All imaging studies were reviewed by a multidisciplinary team including radiologists, pulmonologists, and rheumatologists.

Data collection

Data were systematically extracted from electronic medical records using a standardized data collection form. All variables were predefined and categorized as follows:

General information: Sex (recorded as male or female), age at onset (in years), smoking history (categorized as never smoker or ever smoker, where ever smoker was defined as having smoked at least 100 cigarettes in their lifetime), hypertension (no/yes), diabetes (no/yes), renal function markers, and urinary protein levels (classified as negative, 1+, 2+, or 3+).

Biochemical parameters were assessed using standardized laboratory procedures: Hemoglobin (HGB, reported in g/L) was measured by colorimetry using automated hematology analyzers (Sysmex XN-series) as an indicator of anemia. White blood cell count (WBC) and platelet count (PLT) were quantified via flow cytometry and impedance-based methods, respectively, on the same platform to evaluate inflammatory status and coagulation function; High-sensitivity C-reactive protein (hs-CRP) was analyzed by immunoturbidimetry (Roche Cobas c501) as a marker of systemic inflammation; Alanine aminotransferase (ALT), and aspartate aminotransferase (AST) were determined enzymatically (Roche Cobas c501) to assess liver function; Serum creatinine was measured using an enzymatic assay for renal function evaluation; Complement components C3 and C4 were detected via immunoturbidimetry (Siemens BN II) to monitor complement activity; Rheumatoid factor (RF) was quantified by immunonephelometry (Siemens BN II) as an indicator of autoimmune activity; Neutrophil percentage (%) was derived from automated differential counts of peripheral blood smears to estimate acute inflammation.

Clinicopathological factors: The following clinicopathological factors were documented. (1) Lung infection within 1 year prior to AAV diagnosis was recorded as present or absent based on corresponding ICD-10 diagnostic codes (15). (2) Symptom frequency during treatment was classified as either intermittent or persistent. Intermittent was defined by symptoms occurring in discrete episodes with intervening remission periods of at least 2 weeks, whereas persistent was defined by symptoms that were continuously present. This assessment covered the six-month period preceding diagnosis. (3) Limb numbness or paresthesia was recorded as present or absent based on documentation in the medical records.

All data extraction was conducted independently by two researchers. Any discrepancies were resolved through consensus or by consultation with a third reviewer. The primary outcome was the presence of pulmonary involvement at the time of AAV diagnosis, defined according to composite criteria incorporating clinical symptoms, imaging results, and pulmonary function tests, as detailed in Section 2.3 (Diagnostic criteria for pulmonary involvement).

Statistical analysis

Categorical variables are presented as frequencies and percentages, continuous variables as mean ± standard deviation (SD) or median with interquartile range (IQR), based on distribution. Group comparisons used chi-square or Fisher’s exact test for categorical variables, and Student’s t-test or Mann–Whitney U test for continuous variables. A two-tailed p < 0.05 was considered statistically significant.

Variable selection was performed using Least Absolute Shrinkage and Selection Operator (LASSO) regression with 10-fold cross-validation (R version 4.2.3, glmnet package 4.1.8). Univariate logistic regression was used to further investigate the predictors found by LASSO. In order to evaluate their independent associations with pulmonary involvement, variables with p < 0.05 in univariate analysis were added to a multivariate logistic regression model (SPSS 25.0). The results were presented as odds ratios (ORs) with 95% confidence intervals (CIs).

Restricted cubic spline (RCS) analysis with 3 knots was used to examine the non-linear relationship between hemoglobin (continuous) and pulmonary involvement, adjusted for age (R package rms). A nomogram was constructed based on the final multivariate model. Model performance was evaluated using receiver operating characteristic (ROC) curves, calibration plots, and decision curve analysis (DCA). The dataset was randomly split into training (70%) and validation (30%) sets. To further mitigate overfitting and obtain robust performance estimates, internal validation was performed using 500 bootstrap samples on the training set.

Results

General characteristics and univariate analysis of AAV patients with or without pulmonary involvement

A total of 71 patients diagnosed with AAV were included in the analysis, of whom 38 had developed pulmonary involvement and 33 had not (see Table 1). The age of onset was significantly higher in the pulmonary involvement group (66.61 ± 8.95 years) than in the non-pulmonary involvement group (55.09 ± 19.69 years; t = −3.09, p = 0.003). Hemoglobin levels also differed significantly between the groups, with the pulmonary involvement group displaying lower levels (89.97 ± 22.12 g/L) compared to the non-pulmonary involvement group (105.61 ± 27.19 g/L, t = 2.67, p = 0.009). In contrast, no significant differences were identified for other variables, including complement C4 levels, quantitative immunoglobulin (IgG), and erythrocyte sedimentation rate (ESR) (all p > 0.05).

Table 1
www.frontiersin.org

Table 1. General characteristics and clinicopathological factors of patients with and without pulmonary involvement in AAV.

LASSO regression analysis

LASSO regression analysis was conducted to identify the optimal risk factors for pulmonary involvement in AAV patients (see Figure 1). The analysis determined a lambda.min of 0.032, representing the most regularized model within one standard error of the minimum cross-validation error. The coefficient profiles and the selection of lambda are illustrated in Figure 1, providing a clear visualization of the process used to identify the most impactful predictors.

Figure 1
Panel A shows a coefficient plot against log lambda with multiple curves labeled from one to ten. Panel B displays a binomial deviance plot against log lambda, with a red dotted curve representing deviance and vertical lines indicating variability.

Figure 1. Optimal risk factors identified through LASSO regression analysis. (A) LASSO coefficient profile plot for the 10 categorical variables. The coefficient profile is plotted against the logarithmic sequence of λ [log(λ)]. (B) Selection of the optimal λ parameter using 10-fold cross-validation in the LASSO model. The left dashed line represents the lambda value with the minimum error (lambda.min), while the right dashed line represents the lambda value at one standard error from the minimum (lambda.1-SE).

Multivariable logistic regression analysis

Multivariable logistic regression analysis revealed that intermittent symptom frequency during treatment and the presence of limb numbness or paresthesia were significantly associated with an increased risk of pulmonary involvement (see Table 2). Patients experiencing intermittent symptoms had an over eight-fold increased risk (OR = 8.52, 95% CI: 2.17–33.46, p = 0.002), while those with limb numbness or paresthesia exhibited a twelve-fold increase in risk (OR = 12.21, 95% CI: 1.98–75.50, p = 0.007). Additionally, age at onset and hemoglobin levels were significant predictors. Each additional year of age increased the risk by 5% (OR = 1.05, 95% CI: 1.01–1.10, p = 0.043), and for each 1 g/L decrease in hemoglobin, the risk increased by 3% (OR = 0.97, 95% CI: 0.95–0.99, p = 0.038).

Table 2
www.frontiersin.org

Table 2. Univariate and multivariate logistic regression analysis of significant factors associated with pulmonary involvement in AAV.

Nomogram establishment and evaluation of risk prediction model

The performance of the predictive model was evaluated using ROC curves, calibration plots, and decision curve analysis (see Figure 2). The nomogram (Figure 2A) visually represents the contributions of key predictors, such as age at onset, hemoglobin levels, symptom frequency, and limb numbness or paresthesia, allowing for individualized risk estimation of pulmonary involvement in AAV patients. The model demonstrated strong discriminative power, with an AUC of 0.88 (95% CI: 0.77–0.98) in the training dataset (Figure 2B) and 0.83 (95% CI: 0.64–1.00) in the validation dataset (Figure 2E). Calibration curves for both datasets (Figures 2C,F) showed excellent agreement between predicted and observed outcomes, with Hosmer-Lemeshow p-values of 0.814 (training) and 0.403 (validation), indicating good model calibration. Decision curve analysis (Figures 2D,G) demonstrated a high net benefit across a range of risk thresholds, supporting the model’ s clinical utility in identifying high-risk patients for tailored interventions.

Figure 2
Panel A shows a clinical scoring system for predicting risk based on age of onset, hemoglobin levels, symptom frequency, and limb numbness. B and E are ROC curves showing model sensitivity and specificity with AUC values of 0.88 and 0.83, respectively. C and F display calibration plots with Hosmer-Lemeshow tests, indicating model fit. D and G show decision curves comparing net benefit for different risk thresholds for models, all patients, and no patients.

Figure 2. Nomogram and performance evaluation parameters for the risk prediction model of pulmonary involvement in AAV. (A) Nomogram for predicting the risk of pulmonary involvement in AAV. To use the nomogram: for each predictor (e.g., Age at onset), locate the patient’s value on the corresponding axis, draw a vertical line upward to the “Points” line to read the assigned points. Sum the points for all predictors to obtain the “Total Points.” Finally, locate the total points on the bottom axis and draw a vertical line downward to read the predicted probability of pulmonary involvement on the “Risk” axis. (B) ROC curve for the training set in the risk prediction model. (C) Calibration curve for the training set in the risk prediction model. (D) Decision curve analysis (DCA) for the training set in the risk prediction model. (E) ROC curve for the validation set in testing the risk prediction model. (F) Calibration curve for the validation set in testing the risk prediction model. (G) DCA curve for the validation set in testing the risk prediction model.

Restricted cubic spline analysis for age and hemoglobin levels

The relationship between hemoglobin levels and pulmonary involvement outcomes in AAV patients was further investigated using restricted cubic spline (RCS) analysis. The RCS plot, presented in Figure 3, demonstrates a non-linear association, with lower hemoglobin levels corresponding to an increased risk of pulmonary involvement. Notably, hemoglobin levels below 100 g/L were significantly associated with higher risk, and a steep escalation in risk was observed at levels under 90 g/L. This analysis, adjusted for age as a covariate, highlights the critical threshold of hemoglobin levels in influencing pulmonary outcomes in AAV patients.

Figure 3
A line graph shows the relationship between hemoglobin levels (g/L) and odds ratios (95% confidence interval) from 40 to 160. The shaded area represents the confidence interval. The odds ratio peaks above 2.0 around 60 g/L and declines sharply, stabilizing below 1.0 beyond 100 g/L. P values for overall and nonlinear effects are 0.206 and 0.725, respectively.

Figure 3. RCS analysis of the relationship between hemoglobin levels and pulmonary involvement outcomes in AAV.

Discussion

Pulmonary involvement in patients with AAV considerably increases mortality risk, making early identification of high-risk individuals critical for improved prognosis (16, 17). In this study, regression analysis identified age at onset, hemoglobin level, intermittent symptom pattern, and limb numbness or paresthesia as independent predictors of pulmonary involvement. A clinical risk nomogram was created based on these variables. In both training and validation cohorts, the model showed strong calibration and outstanding discriminatory performance (AUC 0.88 vs. 0.83), and decision curve analysis verified its clinical utility over a broad threshold range. By including commonly available clinical signs, this useful tool may help doctors identify high-risk patients early on, allowing for more proactive intervention or closer monitoring and ultimately improving outcomes in this severe subgroup.

In our group, limb numbness or paresthesia was found to be a significant independent predictor of pulmonary involvement (OR = 12.21, 95% CI: 1.98–75.50). The interpretation of this connection necessitates careful evaluation of our study context, even if it has received less attention in previous literature that specifically focused on AAV-related pulmonary involvement. Due to the very small sample size and the modest number of episodes for this particular symptom, the large confidence range indicates significant uncertainty. However, the association’s strength calls for a mechanistic analysis. It could be explained by similar pathways of systemic microvascular injury, where endothelial damage in the pulmonary vasculature occurs parallel to brain ischemia from peripheral or central nervous system vasculitic involvement (18). These neurological symptoms frequently indicate a more severe and widespread kind of vasculitis, which may exacerbate pulmonary damage by increasing immunological activation and systemic inflammation (19). Therefore, this data implies that a straightforward bedside evaluation of neurological symptoms could be a therapeutically effective trigger to raise suspicion for concurrent pulmonary risk, even with the noted statistical uncertainty. This discovery merits validation in larger, prospective cohorts.

Additionally, there was a significant correlation between an intermittent symptom pattern and a higher risk of pulmonary involvement (OR = 8.52, 95% CI: 2.17–33.46). This pattern, which is marked by periods of remission mixed with recurrent episodes of symptoms, may indicate a more erratic or relapsing form of disease activity. Such instability may expose the pulmonary vasculature to recurrent inflammatory assaults, resulting in cumulative endothelial damage and compromised repair mechanisms, which may eventually increase the risk of overt pulmonary involvement. This idea is corroborated by imaging studies, where greater relapse rates in AAV have been independently linked to specific active CT patterns (such as ground-glass opacities) that indicate persistent inflammation (20). A prolonged pattern of symptoms, on the other hand, may be indicative of a more indolent, chronic inflammatory state that develops gradually and may be linked to various pathophysiological adaptations, which could account for the reduced observed risk of acute pulmonary problems. The idea of disease fluctuation as a marker of severity is consistent with observations in other autoimmune conditions, where relapsing–remitting courses are frequently associated with worse long-term outcomes and greater organ damage accrual, even though it has not been thoroughly investigated in the context of AAV-related pulmonary involvement. Recent findings that acute exacerbation is a substantial risk factor for mortality in patients with AAV-associated interstitial lung disease further support this connection (21). Because they can be a subset with increased systemic vasculitic activity and a higher risk of pulmonary involvement, patients with changing symptoms must be constantly monitored.

A lower hemoglobin level was found to be independently linked to a higher probability of pulmonary involvement in our multivariable analysis. There are probably two mechanisms at work in this connection. Lower hemoglobin level often reflects greater disease burden and systemic inflammation in AAV. More importantly, oxygen supply is directly hampered by its lowering (22). Through mechanisms such hypoxia-inducible factor (HIF) activation, the ensuing tissue hypoxia can increase pro-inflammatory cytokines and worsen oxidative stress, creating a vicious cycle that encourages vascular remodeling and parenchymal damage (23, 24). A crucial clinical threshold was highlighted by our data, which showed a sharp increase in risk at hemoglobin levels below 100 g/L, especially under 90 g/L. This implies that aggressively treating anemia might go beyond supportive treatment and possibly break this harmful pathological cycle.

Our investigation indicated that older age at onset was a major risk factor, which is consistent with earlier results (25). There are a number of interconnected reasons why older patients are more vulnerable. First, a significant influence is probably played by immunosenescence, which is defined by a dysregulated inflammatory response and compromised immunological homeostasis. Second, systemic inflammation may be exacerbated by a greater load of comorbidities like diabetes and hypertension. Lastly, the threshold for inflammatory damage may be lowered by age-related structural alterations in the lung (26, 27).

The nomogram developed in this study integrates these readily available clinical indicators to estimate an individual patient’s risk of pulmonary involvement at diagnosis. It is important to emphasize that the probability generated by the model serves as a decision-support tool to inform initial assessment strategies, not as a definitive diagnostic cutoff. The continuous quantitative risk estimate it provides aims to complement, not replace, comprehensive clinical judgment. Based on the risk distribution in our cohort, a predicted probability exceeding 70% may help identify a high-risk subgroup in whom more extensive pulmonary evaluation—such as prompt high-resolution CT imaging and intensified monitoring—could be prioritized. A lower predicted probability may align with a standard assessment pathway. It is critical to note that this suggested threshold is exploratory and should be validated and calibrated in external populations; clinical decision-making must always incorporate the full patient context. Notably, by incorporating indicators such as neurological symptoms, the model enhances its explanatory power and relevance for bedside evaluation. This tool is designed to aid early risk stratification and support individualized management. If further validated in multicenter studies, its integration into electronic health record systems could improve healthcare workflow efficiency; however, its real-world impact must be evaluated through implementation research in actual clinical practice.

This study has several limitations. First, the single-center retrospective design may introduce selection bias, and the modest sample size constrains statistical power and estimate precision, particularly for predictors with wide confidence intervals. However, this challenge is not uncommon in focused studies of ANCA-associated vasculitis, a rare disease with a recently estimated global pooled incidence of 17.2 per million person-years (28). Comparable exploratory prediction model studies in this field have operated with similar cohort sizes (29). To mitigate overfitting within these constraints, we employed rigorous methods including LASSO regularization and bootstrap internal validation. Second, due to the inherent constraints of retrospective data collection, several important variables—including specific ANCA subtypes (MPO- vs. PR3-ANCA) and genetic markers—were not routinely available in our clinical records and therefore could not be analyzed. This may affect the model’s performance across different patient subtypes. Third, the cross-sectional design precludes assessment of causality or long-term predictive ability. Finally, while internally validated, the model requires external validation in prospective, multi-center cohorts to confirm its generalizability and clinical utility.

Conclusion

In this study, we identified key risk factors for pulmonary involvement in AAV, including intermittent symptoms, limb numbness or paresthesia, age at onset, and hemoglobin levels, and constructed a predictive model based on these factors. The model demonstrated good discriminatory ability, offering a potential tool for integration into clinical workflows to aid in early risk identification and guide personalized treatment. While future multi-center studies are required to validate its applicability across diverse populations, these findings offer useful insights for optimizing management strategies in patients with AAV.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Ethics Committee of the Third Affiliated Hospital of Chongqing Medical University. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants’ legal guardians/next of kin because the waiver of written informed consent was due to this being a retrospective cross-sectional study using fully anonymized existing clinical data of antineutrophil cytoplasmic antibody-associated vasculitis (AAV) patients, practical unfeasibility in contacting all participants, and approval from the Ethics Committee of the Third Affiliated Hospital of Chongqing Medical University.

Author contributions

LZ: Data curation, Conceptualization, Supervision, Writing – original draft, Resources. CL: Software, Investigation, Project administration, Supervision, Writing – review & editing, Formal analysis.

Funding

The author(s) declared that financial support was not received for this work and/or its publication.

Acknowledgments

The authors would like to thank all patients who participated in this study and the members of the research team for their dedicated efforts in clinical data collection, collation, and statistical analysis throughout the research process.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Jennette, JC, Falk, RJ, Bacon, PA, Cid, MC, Ferrario, F, Flores-Suarez, LF, et al. 2012 revised international Chapel Hill consensus conference nomenclature of vasculitides. Arthritis Rheum. (2013) 65:1–11. doi: 10.1002/art.37715,

PubMed Abstract | Crossref Full Text | Google Scholar

2. Mohammad, AJ. An update on the epidemiology of ANCA-associated vasculitis. Rheumatology (Oxford). (2020) 59:iii42–50. doi: 10.1093/rheumatology/keaa089,

PubMed Abstract | Crossref Full Text | Google Scholar

3. Kronbichler, A, Bajema, I, Bruchfeld, A, Mastroianni Kirsztajn, G, and Stone, JH. Diagnosis and management of ANCA-associated vasculitis. Lancet. (2024) 403:683–98. doi: 10.1016/S0140-6736(23)01736-1,

PubMed Abstract | Crossref Full Text | Google Scholar

4. Li, J, Cui, Z, Long, JY, Huang, W, Wang, JW, Wang, H, et al. The frequency of ANCA-associated vasculitis in a national database of hospitalized patients in China. Arthritis Res Ther. (2018) 20:226. doi: 10.1186/s13075-018-1708-7,

PubMed Abstract | Crossref Full Text | Google Scholar

5. Tan, JA, Dehghan, N, Chen, W, Xie, H, Esdaile, JM, and Avina-Zubieta, JA. Mortality in ANCA-associated vasculitis: ameta-analysis of observational studies. Ann Rheum Dis. (2017) 76:1566–74. doi: 10.1136/annrheumdis-2016-210942

Crossref Full Text | Google Scholar

6. Alba, MA, Flores-Suárez, L, Henderson, AG, Xiao, H, Hu, P, Nachman, PH, et al. Interstital lung disease in ANCA vasculitis. Autoimmun Rev. (2017) 16:722–9. doi: 10.1016/j.autrev.2017.05.008,

PubMed Abstract | Crossref Full Text | Google Scholar

7. Sacoto, G, Boukhlal, S, Specks, U, Flores-Suárez, LF, and Cornec, D. Lung involvement in ANCA-associated vasculitis. Presse Med. (2020) 49:104039. doi: 10.1016/j.lpm.2020.104039,

PubMed Abstract | Crossref Full Text | Google Scholar

8. Chalkia, A, Flossmann, O, Jones, R, Nair, JR, Simpson, T, Smith, R, et al. Avacopan for ANCA-associated vasculitis with hypoxic pulmonary haemorrhage. Nephrol Dial Transplant. (2024) 39:1473–82. doi: 10.1093/ndt/gfae020

Crossref Full Text | Google Scholar

9. Collins, GS, Reitsma, JB, Altman, DG, and Moons, KGM. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. BMJ. (2015) 350:g7594–4. doi: 10.1136/bmj.g7594,

PubMed Abstract | Crossref Full Text | Google Scholar

10. Ntatsaki, E, Carruthers, D, Chakravarty, K, D’Cruz, D, Harper, L, Jayne, D, et al. BSR and BHPR guideline for the management of adults with ANCA-associated vasculitis. Rheumatology (Oxford). (2014) 53:2306–9. doi: 10.1093/rheumatology/ket445

Crossref Full Text | Google Scholar

11. Hunder, A. W., and Bloch, D A L Calabrese, H, Fauci, A S, Fries, J F, et al. The American College of Rheumatology 1990 criteria for the classification of vasculitis. Introduction. Arthritis Rheum, 33, 1065–1067 (1990) doi: 10.1002/art.1780330802

Crossref Full Text | Google Scholar

12. Hellmich, S-AB, Schirmer, JH, Hellmich, B, Sanchez-Alamo, B, Berti, A, Blockmans, D, et al. EULAR recommendations for the management of ANCA-associated vasculitis: 2022 update. Ann Rheum Dis. (2024) 83:30–47. doi: 10.1136/ard-2022-223764,

PubMed Abstract | Crossref Full Text | Google Scholar

13. Dieterich, J, Hellmich, B, Mahrhold, J, Feng, YS, El Rai, A, Nessyt, F, et al. Pulmonary function in patients with ANCA-associated vasculitis. Sarcoidosis Vasc Diffuse Lung Dis. (2024) 41:e2024025. doi: 10.36141/svdld.v41i2.15577

Crossref Full Text | Google Scholar

14. Stanojevic, S, Kaminsky, DA, Miller, MR, Thompson, B, Aliverti, A, Barjaktarevic, I, et al. ERS/ATS technical standard on interpretive strategies for routine lung function tests. Eur Respir J. (2022) 60:2101499. doi: 10.1183/13993003.01499-2021

Crossref Full Text | Google Scholar

15. World Health Organization. International statistical classification of diseases and related health problems, 10th revision (ICD-10). Geneva: World Health Organization (2016).

Google Scholar

16. Zhou, P, Li, Z, Gao, L, Que, C, Li, H, Ma, J, et al. Pulmonary involvement of ANCA-associated vasculitis in adult Chinese patients. BMC Pulm Med. (2022) 22:35. doi: 10.1186/s12890-022-01829-y,

PubMed Abstract | Crossref Full Text | Google Scholar

17. Matsuda, KT, Kotani, T, Suzuka, T, Kiboshi, T, Fukui, K, Wakama, M, et al. Evaluation of poor prognostic factors of respiratory related death in microscopic polyangiitis complicated by interstitial lung disease. Sci Rep. (2021) 11:1490. doi: 10.1038/s41598-021-81311-7

Crossref Full Text | Google Scholar

18. Suppiah, HR, Batra, R, Suppiah, R, Hadden, RDM, Arden, NK, Collins, MP, et al. Peripheral neuropathy in ANCA-associated vasculitis: outcomes from the European Vasculitis study group trials. Rheumatology (Oxford). (2011) 50:2214–22. doi: 10.1093/rheumatology/ker266,

PubMed Abstract | Crossref Full Text | Google Scholar

19. Fouka, DF, and Steiropoulos, P. Pathogenesis of pulmonary manifestations in ANCA-associated Vasculitis and Goodpasture syndrome. Int J Mol Sci. (2024) 25:5278. doi: 10.3390/ijms25105278,

PubMed Abstract | Crossref Full Text | Google Scholar

20. Vogt, K, Fink, CB, Schreibing, TM, Krämer, S, Reinartz, S, and Rauen, T. Distinct pulmonary patterns in ANCA-associated vasculitides: insights from a retrospective single center cohort study. Rheumatol Int. (2024) 44:2435–43. doi: 10.1007/s00296-024-05664-8,

PubMed Abstract | Crossref Full Text | Google Scholar

21. He, X, Yuan, W, Yang, Y, Ji, J, and Chen, X. Risk factors for poor prognosis in ANCA-associated vasculitis with interstitial lung disease: a systematic review and meta-analysis. Clin Rheumatol. (2025) 44:1675–89. doi: 10.1007/s10067-025-07378-z,

PubMed Abstract | Crossref Full Text | Google Scholar

22. Buehler, PW, Baek, JH, Lisk, C, Connor, I, Sullivan, T, Kominsky, D, et al. Free hemoglobin induction of pulmonary vascular disease: evidence for an inflammatory mechanism. Am J Physiol Lung Cell Mol Physiol. (2012) 303:L312–26. doi: 10.1152/ajplung.00074.2012

Crossref Full Text | Google Scholar

23. Gresnigt, RA, Rasid, O, Savers, A, Jouvion, G, Dannaoui, E, Parlato, M, et al. Reducing hypoxia and inflammation during invasive pulmonary aspergillosis by targeting the interleukin-1 receptor. Sci Rep. (2016) 6:26490. doi: 10.1038/srep26490

Crossref Full Text | Google Scholar

24. Angelini, DJ, Su, Q, Yamaji-Kegan, K, Fan, C, Skinner, JT, Champion, HC, et al. Hypoxia-induced mitogenic factor (HIMF/FIZZ1/RELMalpha) induces the vascular and hemodynamic changes of pulmonary hypertension. Am J Physiol Lung Cell Mol Physiol. (2009) 296:L582–93. doi: 10.1152/ajplung.90526.2008

Crossref Full Text | Google Scholar

25. Kumar, BE. Age-related decline in immunity: implications for vaccine responsiveness. Expert Rev Vaccines. (2008) 7:467–79. doi: 10.1586/14760584.7.4.467,

PubMed Abstract | Crossref Full Text | Google Scholar

26. Pera, A, Campos, C, López, N, Hassouneh, F, Alonso, C, Tarazona, R, et al. Immunosenescence: implications for response to infection and vaccination in older people. Maturitas. (2015) 82:50–5. doi: 10.1016/j.maturitas.2015.05.004

Crossref Full Text | Google Scholar

27. Cho, WK, Lee, CG, and Kim, LK. COPD as a disease of Immunosenescence. Yonsei Med J. (2019) 60:407–13. doi: 10.3349/ymj.2019.60.5.407,

PubMed Abstract | Crossref Full Text | Google Scholar

28. Redondo-Rodriguez, R, Mena-Vázquez, N, Cabezas-Lucena, AM, Manrique-Arija, S, Mucientes, A, and Fernández-Nebro, A. Systematic review and Metaanalysis of worldwide incidence and prevalence of Antineutrophil cytoplasmic antibody (ANCA) associated Vasculitis. J Clin Med. (2022) 11:573. doi: 10.3390/jcm11092573,

PubMed Abstract | Crossref Full Text | Google Scholar

29. Guo, Q, Shao, Y, Yu, L, Zhang, X, Shang, J, Feng, X, et al. Development of a nomogram model to predict mortality in ANCA-associated Vasculitis patients with pulmonary involvement. Clin Respir J. (2025) 19:e70067. doi: 10.1111/crj.70067,

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: ANCA-associated vasculitis, pulmonary involvement, predictive model, risk factors, LASSO regression, hemoglobin, intermittent symptoms

Citation: Zheng L and Liu C (2026) Exploration of factors influencing pulmonary injury in patients with anti-neutrophil cytoplasmic antibody-associated vasculitis and construction of a predictive model. Front. Med. 12:1724495. doi: 10.3389/fmed.2025.1724495

Received: 14 October 2025; Revised: 09 December 2025; Accepted: 09 December 2025;
Published: 02 January 2026.

Edited by:

Gianluca Sambataro, Kore University of Enna, Italy

Reviewed by:

Carmen Silvia Valente Barbas, University of São Paulo, Brazil
Elena Treppo, University of Udine, Italy

Copyright © 2026 Zheng and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chongyang Liu, NjUxMjAyQGhvc3BpdGFsLmNxbXUuZWR1LmNu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.