Journal Information
Share
Share
Download PDF
More article options
Visits
499
Original Article
Full text access
Available online 4 November 2024
Effect of the Antibody-mediated Immune Responses on COPD, Asthma, and Lung Function: A Mendelian Randomization Study
Visits
499
Xu Guixinga, Liu Yilinb, Fan Huayinga, Liang Fanrongb, Li Dehuaa,
Corresponding author
si_ff@163.com

Corresponding author.
a Hospital of Chengdu University of Traditional Chinese Medicine, Chengdu, Sichuan, China
b Acupuncture and Tuina School, Chengdu University of Traditional Chinese Medicine Chinese Medicine, Chengdu, Sichuan, China
This item has received
Received 18 July 2024. Accepted 14 October 2024
Article information
Abstract
Full Text
Bibliography
Download PDF
Statistics
Figures (4)
Show moreShow less
Tables (1)
Table 1. MR Analysis of Antibody-mediated Immune Responses on COPD, Asthma, and Lung Function.
Additional material (1)
Abstract
Introduction

The precise cause of antibody-mediated immune responses on chronic obstructive pulmonary disease (COPD), asthma, and lung function remains unclear. We characterized the relationship between antibody-mediated immune responses to COPD, asthma, and lung function, ultimately achieve the prevention or treatment.

Methods

We obtained summary data from published genome-wide association studies, including antibody-mediated immune responses, COPD, asthma, forced expiratory volume in the first second (FEV1), forced expiratory volume (FVC), and FEV1/FVC. Bidirectional two-sample mendelian randomization (MR) analysis was used to assess causal relationships of antibody-mediated immune responses, COPD, asthma, FEV1, FVC, and FEV1/FVC.

Results

A total of 20 antibody-mediated immune responses were identified have a significant causal effect on COPD, asthma, FEV1, and FVC, with six exhibiting reverse causality. Importantly, the results of the five MR analyses were almost identical with respect to the causal effect of anti-polyomavirus 2 IgG seropositivity and varicella zoster virus glycoprotein E and I antibody levels on the risk of COPD, asthma, FEV1, and FVC.

Conclusions

This study contributes to existing knowledge by investigating the causal relationship between antibody-mediated immune responses and respiratory conditions, including COPD, asthma, and lung function, using a two-sample MR design. The key findings can aid in identifying individuals at risk of these conditions and facilitate early prevention and diagnosis.

Keywords:
Antibody-mediated immune responses
COPD
Asthma
Lung function
Mendelian randomization
Full Text
Introduction

Chronic obstructive pulmonary disease (COPD) and asthma are two common progressive lung diseases, asthma affects more than 300 million people worldwide and its prevalence is increasing,1 and COPD is the third leading cause of death worldwide.2,3 Asthma and COPD are complex inflammatory diseases of the airways, which possibly related to impaired host defenses lead to increased susceptibility to pathogens, pollutants and allergens.4 Decreased lung function is the main clinical feature of asthma and COPD. Immune system may be a key driver in the pathogenesis of COPD, and the immune response is significantly associated with acute exacerbations of COPD.5 Evidence from observational studies suggests that more immune cells are present in the lung tissue of COPD and asthma sufferers compared to healthy individuals, along with a significant immune response.6 Mendelian randomization (MR) studies found that a variety of immune cells and inflammatory factors are risk factors for the development of COPD and asthma.7,8 However, changes in immune cells are the process by which an immune response occurs in the body, and immune cells recognise foreign substances such as antigens and activate, proliferate and differentiate to form effector cells and produce effector molecules to remove the foreign substances. In this process, foreign substances such as antigens are the initiators of the body's immune response, while pathogenic microorganisms are the most common cause of the body's immune response. Antibody-mediated immune responses persist in all stages of asthma and COPD from initial to end-stage,9 but whether there is a clear causal relationship between the development of COPD and asthma remains unclear. It is suggested that antibody-mediated immune responses may be causally related to the development of COPD and asthma. Guillaume et al. measured serological measurements of 20 different microorganisms in 9724 participant serum samples, selected the ones with a seroprevalence of >15% for genome-wide association studies (GWAS).10 Whether these detected antibody-mediated immune responses are causally associated with the development of COPD, asthma, and reduced lung function remains unclear. MR is a potential method of causal inference aimed at estimating the effect of exposure factors on outcomes and being able to control for potential confounders, thus avoiding reverse causation bias.11

The objective of this study was to investigate the causal relationship between antibody-mediated immune responses and COPD, asthma and lung function by bidirectional two-sample MR analysis.

MethodsStudy Design

The study flow chart is shown in Fig. 1. First, we obtained summary data from published GWAS, including features such as antibody-mediated immune responses, COPD, asthma, forced expiratory volume in the first second (FEV1), forced expiratory volume (FVC), and FEV1/FVC ratio (Supplementary Table S1). Second, two-sample MR analysis was used to assess causal relationships of antibody-mediated immune responses to COPD, asthma, FEV1, FVC, and FEV1/FVC ratio. Finally, reversed two-sample MR analysis was used to explore the COPD, asthma, FEV1, FVC, and FEV1/FVC ratio on antibody-mediated immune responses. Our MR study was performed in accordance with STROBE-MR guidelines (Supplementary Table S2).12

Fig. 1.

Study flow chart and the three assumptions that IVs need to fulfil.

(0.13MB).
Data Sources

The GWAS summary statistics of antibody-mediated immune responses were retrieved from the large single homogenous population cohort, more than 500,000 British adults were recruited by CDC between 2006 and 2010, of whom 9724 provided samples of geranium for the measurement of 20 different micro groups.13

The GWAS summary statistics of COPD were retrieved from UK biobank with 2115 European ancestry cases, and 454,233 European ancestry controls.14 Considering that COPD subtypes emphysema and bronchiectasis may have different pathogenesis, we also analyzed the causal relationship between emphysema, bronchiectasis and exposure.14 Also, asthma has shown association with COPD and may have the same pathogenesis and pathological basis, so we also analyzed the causal relationship between asthma and exposure.14

The GWAS summary statistics of lung function, including FEV1, FVC, and FEV1/FVC, were retrieved from 475,645 European ancestry individuals, which is the largest global genomics study of lung function to date, with genome-wide association analyses of FEV1, FVC, and FEV1/FVC in 49 cohorts, adjusting for age, age squared, sex, and height, and applying genomic controls using the linkage disequilibrium (LD) score regression intercept.15

Instrumental Variable Selection

The causal effects of MR analysis are based on three basic assumptions of genetic variation, including: (1) instrumental variables (IVs) are related to exposure factors; (2) IVs are not related to confounding factors; (3) IVs are not related to outcome variables and only affects the outcome variable through exposure factors.

The locus-wide significance threshold of SNP was P<1×10−5, and set P<1×10−5 if the number of single nucleotide polymorphism (SNP) screened after de-linked imbalance is less than 5 for antibody-mediated immune responses, and COPD, asthma, emphysema, bronchiectasis, FEV1, FVC, and FEV1/FVC ratio.16 The SNPs with impact of linkage disequilibrium were removed by clumping process (r2<0.001, window size=10,000kb). Even if each IV satisfies the three assumptions of the instrumental variables, chance may cause the estimates of the IVs to differ significantly from the expected results. Therefore, we used the MR pleiotropy residual sum and outlier (MR-PRESSO) test, and the outliers removed, respectively, to detect horizontal pleiotropy (MR-PRESSO outlier test P-values <0.05), and eliminated the pleiotropy effect (the SNPs were eliminated one by one according to horizontal pleiotropy until the MR-PRESSO global test P-value >0.05).17 SNPs related to outcome variables were excluded (P<0.05). The strength of the selected SNPs was evaluated using F-statistic, where SNPs with F-statistic <10 were excluded to avoid weak instrument bias in the MR analysis.18

Statistical AnalysisTwo-sample Mendelian Randomization

Five MR method was used to evaluate the causal relationship among the antibody-mediated immune responses, and COPD, asthma, emphysema, bronchiectasis, FEV1, FVC, and FEV1/FVC ratio. The inverse variance-weighted (IVW) was primary method to identify the causal relationship due to IVW usually provides the highest statistical power. The IVW method, inspired by meta-analytic methods, statistically combines or weighted linear regressions of the Wald ratio estimates across IVs, with less variance through weighted averaging in favor of the most statistically powerful results.19 And, in response to the possible heterogeneity and genetic pleiotropy of IV in MR analyses, the IVW approach can also assess the differences between causality estimates specific to each locus. In the absence of multinomially, IVW provides virtually unbiased estimates of causality.20 To obtain more conservative and robust estimates, we all used the random effects IVW test. The maximum likelihood,21 MR-Egger,22 weighted median,23,24 and weighted mode were as second methods to infer causal relationship.

The method of reverse MR was used to explore the reverse causality of COPD, asthma, emphysema, bronchiectasis, FEV1, FVC, and FEV1/FVC ratio on antibody-mediated immune responses, and the methods and steps of the analysis were the same as before. All MR analyses were performed in R (version 4.3.1) software, using the “TwoSampleMR” (version 0.5.7) (https://github.com/MRCIEU/TwoSampleMR) and “MR-PRESSO” (version 1.0) (https://github.com/rondolab/MR-PRESSO) packages.17

ResultsCausal Effects of Antibody-mediated Immune Responses on COPD

A total of 13 antibody-mediated immune responses were identified as causally associated with COPD (PIVW<0.05, Fig. 2), and six did not show heterogeneity (PIVW>0.05, PMR-Egger>0.05), and horizontal pleiotropy (P>0.05). Importantly, all five analyses showed that varicella zoster virus glycoproteins E and I antibody levels (odds ratio (OR)=1.91, 95% CI (1.73–2.12), P<0.001), and anti-polyomavirus 2 IgG seropositivity (OR=1.22, 95% CI (1.15–1.29), P<0.001) were a significant risk factor for COPD (Table 1 and Fig. 3; Supplementary Tables S3 and S4). The results of Wald ratio analysis of individual SNPs are shown in Supplementary Table S5.

Fig. 2.

Causal estimates of bidirectional MR between antibody-mediated immune responses, COPD, asthma, emphysema, bronchiectasis, FEV1, FVC, and FEV1/FVC. Top: Estimates from the IVW analysis of antibody-mediated immune responses on COPD, asthma, emphysema, bronchiectasis, FEV1, FVC, and FEV1/FVC. Bottom: Estimates from the IVW analysis of COPD, asthma, emphysema, bronchiectasis, FEV1, FVC, and FEV1/FVC on antibody-mediated immune responses. The antibody-mediated immune responses features underlined in red were related to COPD, asthma, FEV1, and FVC. COPD: chronic obstructive pulmonary disease; FEV1: forced expiratory volume in the first second; FVC: forced vital capacity; MR: Mendelian randomization.

(0.54MB).
Table 1.

MR Analysis of Antibody-mediated Immune Responses on COPD, Asthma, and Lung Function.

Exposure  Outcome  Method  nSNP  Beta  SE  P-value  r2  F-statistic 
Human herpes virus 7 U14 antibody levels  COPD  IVW  0.399  0.141  0.005  0.004  33.645 
Anti-polyomavirus 2 IgG seropositivity  COPD  IVW  19  0.195  0.029  0.000  0.008  70.984 
Polyomavirus 2 JC VP1 antibody levels  COPD  IVW  0.384  0.111  0.001  0.010  50.764 
Anti-Merkel cell polyomavirus IgG seropositivity  COPD  IVW  21  0.190  0.035  0.000  0.007  59.298 
Merkel cell polyomavirus VP1 antibody levels  COPD  IVW  11  0.314  0.071  0.000  0.009  50.017 
Toxoplasma gondii p22 antibody levels  COPD  IVW  38  0.055  0.028  0.046  0.017  21.600 
Anti-varicella zoster virus IgG seropositivity  COPD  IVW  29  0.116  0.048  0.015  0.003  25.067 
Varicella zoster virus glycoproteins E and I antibody levels  COPD  IVW  27  0.649  0.051  0.000  0.007  51.944 
Epstein-Barr virus EA-D antibody levels  COPD  IVW  14  −0.415  0.149  0.005  0.006  43.227 
Epstein-Barr virus EBNA-1 antibody levels  COPD  IVW  44  −0.219  0.055  0.000  0.010  73.901 
Epstein-Barr virus VCA p18 antibody levels  COPD  IVW  20  −0.422  0.082  0.000  0.006  49.217 
Epstein-Barr virus ZEBRA antibody levels  COPD  IVW  32  −0.252  0.079  0.001  0.009  68.712 
Human herpes virus 6 IE1B antibody levels  COPD  IVW  32  −0.190  0.081  0.019  0.003  22.618 
Chlamydia trachomatis PorB antibody levels  Emphysema  IVW  89  0.014  0.009  0.106  0.077  22.003 
Anti-polyomavirus 2 IgG seropositivity  Bronchiectasis  IVW  19  −0.296  0.058  0.000  0.008  70.984 
Polyomavirus 2 JC VP1 antibody levels  Bronchiectasis  IVW  −0.581  0.156  0.000  0.010  50.764 
Anti-varicella zoster virus IgG seropositivity  Bronchiectasis  IVW  29  −0.180  0.063  0.004  0.003  25.067 
Varicella zoster virus glycoproteins E and I antibody levels  Bronchiectasis  IVW  27  −0.259  0.103  0.012  0.007  51.944 
Epstein-Barr virus EBNA-1 antibody levels  Bronchiectasis  IVW  44  0.377  0.082  0.000  0.010  73.901 
Epstein-Barr virus VCA p18 antibody levels  Bronchiectasis  IVW  20  0.446  0.132  0.001  0.006  49.217 
Epstein-Barr virus ZEBRA antibody levels  Bronchiectasis  IVW  32  −0.299  0.099  0.002  0.009  68.712 
Anti-human herpes virus 6 IE1B IgG seropositivity  Bronchiectasis  IVW  22  −0.201  0.072  0.005  0.003  21.847 
Herpes simplex virus 2 mgG-1 antibody levels  Asthma  IVW  23  0.223  0.107  0.037  0.016  22.226 
Anti-polyomavirus 2 IgG seropositivity  Asthma  IVW  19  0.445  0.079  0.000  0.008  70.984 
Polyomavirus 2 JC VP1 antibody levels  Asthma  IVW  1.097  0.250  0.000  0.010  50.764 
Anti-Merkel cell polyomavirus IgG seropositivity  Asthma  IVW  21  0.305  0.081  0.000  0.007  59.298 
Anti-varicella zoster virus IgG seropositivity  Asthma  IVW  29  0.268  0.083  0.001  0.003  25.067 
Varicella zoster virus glycoproteins E and I antibody levels  Asthma  IVW  27  0.568  0.147  0.000  0.007  51.944 
Chlamydia trachomatis PorB antibody levels  Asthma  IVW  89  0.070  0.026  0.007  0.077  22.003 
Epstein-Barr virus EBNA-1 antibody levels  Asthma  IVW  44  −0.982  0.098  0.000  0.010  73.901 
Epstein-Barr virus VCA p18 antibody levels  Asthma  IVW  20  −0.973  0.194  0.000  0.006  49.217 
Epstein-Barr virus ZEBRA antibody levels  Asthma  IVW  32  0.305  0.142  0.032  0.009  68.712 
Human herpes virus 7 U14 antibody levels  FEV1  IVW  −0.073  0.034  0.032  0.004  35.716 
Anti-polyomavirus 2 IgG seropositivity  FEV1  IVW  22  −0.045  0.006  0.000  0.008  66.528 
Polyomavirus 2 JC VP1 antibody levels  FEV1  IVW  −0.094  0.018  0.000  0.011  56.507 
Anti-Merkel cell polyomavirus IgG seropositivity  FEV1  IVW  22  −0.035  0.007  0.000  0.007  59.810 
Anti-varicella zoster virus IgG seropositivity  FEV1  IVW  30  −0.030  0.008  0.000  0.003  24.974 
Varicella zoster virus glycoproteins E and I antibody levels  FEV1  IVW  30  −0.164  0.009  0.000  0.007  52.539 
Chlamydia trachomatis momp A antibody levels  FEV1  IVW  36  −0.004  0.002  0.011  0.023  21.738 
Epstein-Barr virus EA-D antibody levels  FEV1  IVW  15  0.094  0.027  0.000  0.006  42.686 
Epstein-Barr virus EBNA-1 antibody levels  FEV1  IVW  45  0.059  0.009  0.000  0.010  73.699 
Epstein-Barr virus VCA p18 antibody levels  FEV1  IVW  22  0.108  0.015  0.000  0.006  49.075 
Human herpes virus 6 IE1A antibody levels  FEV1  IVW  33  0.030  0.012  0.013  0.003  23.058 
Anti-polyomavirus 2 IgG seropositivity  FVC  IVW  22  −0.046  0.005  0.000  0.008  66.528 
Polyomavirus 2 JC VP1 antibody levels  FVC  IVW  −0.096  0.009  0.000  0.011  56.507 
Anti-Merkel cell polyomavirus IgG seropositivity  FVC  IVW  22  −0.039  0.005  0.000  0.007  59.810 
Merkel cell polyomavirus VP1 antibody levels  FVC  IVW  14  −0.029  0.011  0.011  0.008  48.161 
Anti-varicella zoster virus IgG seropositivity  FVC  IVW  30  −0.013  0.005  0.006  0.003  24.974 
Varicella zoster virus glycoproteins E and I antibody levels  FVC  IVW  30  −0.100  0.010  0.000  0.007  52.539 
Chlamydia trachomatis momp A antibody levels  FVC  IVW  36  −0.004  0.002  0.016  0.023  21.738 
Epstein-Barr virus EBNA-1 antibody levels  FVC  IVW  45  0.055  0.007  0.000  0.010  73.699 
Epstein-Barr virus VCA p18 antibody levels  FVC  IVW  22  0.087  0.012  0.000  0.006  49.075 
Anti-BK polyomavirus IgG seropositivity  FVC  IVW  16  0.004  0.002  0.016  0.002  21.063 
Human herpes virus 7 U14 antibody levels  FEV1/FVC  IVW  −0.068  0.030  0.025  0.004  35.716 
Merkel cell polyomavirus VP1 antibody levels  FEV1/FVC  IVW  14  0.041  0.015  0.005  0.008  48.161 
Anti-varicella zoster virus IgG seropositivity  FEV1/FVC  IVW  30  −0.039  0.010  0.000  0.003  24.974 
Varicella zoster virus glycoproteins E and I antibody levels  FEV1/FVC  IVW  30  −0.166  0.010  0.000  0.007  52.539 
Cytomegalovirus pp150 antibody levels  FEV1/FVC  IVW  25  −0.011  0.005  0.019  0.004  21.167 
Epstein-Barr virus EA-D antibody levels  FEV1/FVC  IVW  15  0.170  0.027  0.000  0.006  42.686 
Epstein-Barr virus VCA p18 antibody levels  FEV1/FVC  IVW  22  0.061  0.025  0.015  0.006  49.075 
Human herpes virus 6 IE1A antibody levels  FEV1/FVC  IVW  33  0.038  0.015  0.014  0.003  23.058 
Anti-human herpes virus 6 IE1B IgG seropositivity  FEV1/FVC  IVW  22  −0.021  0.011  0.043  0.003  21.737 

Note: COPD, asthma, FEV1, FVC. COPD: chronic obstructive pulmonary disease; FEV1: forced expiratory volume in the first second; FVC: forced expiratory volume; MR: Mendelian randomization; IVW: inverse variance-weighted; SE: standard error.

Fig. 3.

Forest plots for causal effects of antibody-mediated immune responses on COPD, asthma, emphysema, bronchiectasis, FEV1, FVC, and FEV1/FVC. The horizontal bars correspond to the estimated OR with 95% CI using the IVW and MR-Egger method. COPD, asthma, FEV1, FVC. COPD: chronic obstructive pulmonary disease; FEV1: forced expiratory volume in the first second; FVC: forced vital capacity; MR: Mendelian randomization; OR: odds ratio; CI, confidence interval.

(1.74MB).

Only one antibody-mediated immune responses were identified as causally associated with emphysema (PIVW<0.05, Fig. 2), did not show heterogeneity (PIVW>0.05, PMR-Egger>0.05), and horizontal pleiotropy (P>0.05) (Table 1; S3 and S4 in Supplementary). The Chlamydia trachomatis PorB antibody levels were risk factor on emphysema, but there are differences in the results of different analyses (Fig. 3). The results of Wald ratio analysis of individual SNPs are shown in Supplementary Table S5.

A total of eight antibody-mediated immune responses were identified as causally associated with bronchiectasis (PIVW<0.05, Fig. 2), and six did not show heterogeneity (PIVW>0.05, PMR-Egger>0.05), and horizontal pleiotropy (P>0.05) (Table 1; S3 and S4 in Supplementary). Importantly, all five analyses showed that anti-polyomavirus 2 IgG seropositivity (OR=0.74, 95% CI (0.66–0.83), P<0.001) was a significant protect factor for bronchiectasis (Fig. 3). The results of Wald ratio analysis of individual SNPs are shown in Supplementary Table S5.

Causal Effects of Antibody-mediated Immune Responses on Asthma

A total of 10 antibody-mediated immune responses were identified as causally associated with asthma (PIVW<0.05, Fig. 2), and all trait did not show heterogeneity (PIVW>0.05, PMR-Egger>0.05), and horizontal pleiotropy (P>0.05) (Table 1; S3 and S4 in Supplementary). Importantly, all five analyses showed that anti-polyomavirus 2 IgG seropositivity (OR=0.74, 95% CI (0.66–0.83), P<0.001) was a significant risk factor for asthma, and Epstein-Barr virus EBNA-1 antibody levels (OR=0.38, 95% CI (0.31–0.45), P<0.001) was a significant protect factor for asthma (Fig. 3). The results of Wald ratio analysis of individual SNPs are shown in Supplementary Table S5.

Causal Effects of Antibody-mediated Immune Responses on FEV1, FVC, and FEV1/FVC Ratio

A total of 11 antibody-mediated immune responses were identified as causally associated with FEV1 (PIVW<0.05, Fig. 2), and only one did not show heterogeneity (PIVW>0.05, PMR-Egger>0.05), and horizontal pleiotropy (P>0.05) (Table 1; S3 and S4 in Supplementary). All five analyses showed that anti-polyomavirus 2 IgG seropositivity (OR=0.96, 95% CI (0.94–0.97), P<0.001), and varicella zoster virus glycoproteins E and I antibody levels (OR=0.85, 95% CI (0.83–0.86), P<0.001) were a significant risk factor for FEV1 (Fig. 3). The results of Wald ratio analysis of individual SNPs are shown in Supplementary Table S5.

A total of 10 antibody-mediated immune responses were identified as causally associated with FVC (PIVW<0.05, Fig. 2), and two did not show heterogeneity (PIVW>0.05, PMR-Egger>0.05), and horizontal pleiotropy (P>0.05) (Table 1; S3 and S4 in Supplementary). All five analyses showed that anti-polyomavirus 2 IgG seropositivity (OR=0.95, 95% CI (0.95–0.96), P<0.001), and varicella zoster virus glycoproteins E and I antibody levels (OR=0.90, 95% CI (0.89–0.92), P<0.001) were a significant risk factor for FVC, and Epstein-Barr virus VCA p18 antibody levels (OR=1.09, 95% CI (1.07–1.12), P<0.001) was a significant protect factor for FVC (Fig. 3). The results of Wald ratio analysis of individual SNPs are shown in Supplementary Table S5.

A total of nine antibody-mediated immune responses were identified as causally associated with FEV1/FVC (PIVW<0.05, Fig. 2), and one did not show heterogeneity (PIVW>0.05, PMR-Egger>0.05), and horizontal pleiotropy (P>0.05) (Table 1; S3 and S4 in Supplementary). All five analyses showed that varicella zoster virus glycoproteins E and I antibody levels (OR=0.85, 95% CI (0.83–0.86), P<0.001) was a significant risk factor for FVC, and Merkel cell polyomavirus VP1 antibody levels (OR=1.04, 95% CI (1.01–1.07), P=0.005) was a significant protect factor for FVC (Fig. 3). The results of Wald ratio analysis of individual SNPs are shown in Supplementary Table S5.

Causal Effects of COPD, Asthma, and Lung Function on Antibody-mediated Immune Responses

The reverse MR analysis found only a small number of causal relationships of COPD, asthma, and lung function on antibody-mediated immune responses, and these varied considerably when analyzed by different analytical methods (Figs. 2 and 4; Supplementary Tables S6 and S7), the results of Wald ratio analysis of individual SNPs are shown in Supplementary Table S8. COPD, emphysema, bronchiectasis, asthma, FEV1, FVC, and FEV1/FVC were identified as causally associated with 1, 2, 3, 7, 2, 3, and 4 antibody-mediated immune responses (PIVW<0.05, Fig. 4), respectively. The results of Wald ratio analysis of individual SNPs are shown in Supplementary Table S8.

Fig. 4.

Forest plots for causal effects of COPD, asthma, emphysema, bronchiectasis, FEV1, FVC, and FEV1/FVC on antibody-mediated immune responses. The horizontal bars correspond to the estimated OR with 95% CI using the IVW and MR-Egger method. COPD, asthma, FEV1, FVC. COPD: chronic obstructive pulmonary disease; FEV1: forced expiratory volume in the first second; FVC: forced vital capacity; MR: Mendelian randomization; OR: odds ratio; CI: confidence interval.

(0.72MB).
Discussion

To our knowledge, this is the first study to analyses the causal relationship between antibody-mediated immune responses and COPD patients, asthmatics, and lung function. The results show that a total of 20 antibody-mediated immune responses that could have a significant causal effect on COPD, asthma, FEV1, and FVC. There was a total of seven antibody-mediated immune responses that could have a significant effect on COPD, asthma, FEV1, FVC at the same time, namely anti-Merkel cell polyomavirus IgG seropositivity, anti-polyomavirus 2 IgG seropositivity, anti-varicella zoster virus IgG seropositivity, Epstein-Barr virus EBNA-1 antibody levels, Epstein-Barr virus VCA p18 antibody levels, polyomavirus 2 JC VP1 antibody levels, and varicella zoster virus glycoproteins E and I antibody levels. More importantly, we find little difference between the results of the different analyses in terms of the risk causal impact of anti-polyomavirus 2 IgG seropositivity, and varicella zoster virus glycoproteins E and I antibody levels on COPD, asthma, FEV1, and FVC.

Infectious diseases are a major contributor to the global burden of disease and can have some impact on systemic diseases, such as fatigue, dyspnoea, pain, anxiety, depression and cognitive dysfunction following infection with COVID-19.25 And as this class of viruses is usually preventable, clarifying the causal relationship between viruses and COPD, asthma, and lung function will contribute to better clinical management and treatment of patients with COPD and asthma. The results of this study showed that polyomavirus and varicella zoster virus-induced immune responses were risk factors for COPD, asthma, FEV1, and FVC, and the results were consistent across the five MR analysis methods.

Polyomaviruses are small DNA viruses that are widely distributed in nature. Polyomaviruses usually infect immunocompromised individuals and can be transmitted through the respiratory system. Their infections are common and they do not show symptoms after infection.26 It has been shown that the CD8+ T cell response is specific for polyomavirus and is associated with clinical improvement in patients with polyomavirus infection.27 CD8+ T cells were found to be increased in mild to moderate COPD lungs and may contribute to inflammation before severe disease.28 Recent findings have identified the presence of polyomavirus infections in COPD patients.29 It is suggested that polyomavirus may have elicited the CD8+ T cell response that led to exacerbation in COPD patients. Polyomavirus has been thought to be closely associated with CD4+ T with B cells,30 and there is an interactive relationship between CD4+ T and CD8+ T,31 as well as CD4+ T driving B cell differentiation. And B cells may be potential transmitters of polyomavirus or involved in the control of polyomavirus.32 CD4+ T can be found to be crucial in the cellular immune response to polyomavirus, and Cao et al. found that CD4+ T is a protective factor against COPD.8 Some findings suggest that polyomavirus infection or persistence may not lead to the development of COPD.33,34 This is inconsistent with the results of the present study, and we believe that the possible reason for this is that all their studies were observational, and there may be errors in the findings.

The varicella zoster virus is left in the body after experiencing chickenpox and reactivates years later in a form called shingles. It can often occur in people whose immune systems are not working properly, such as in respiratory diseases like asthma and COPD.35 For varicella zoster virus, several guidelines have recommended vaccination against varicella zoster virus, and herpes zoster vaccination rates may be low in patients aged 50 years ≥COPD, which may indicate a lack of awareness of herpes zoster virus risk factors among clinicians and patients.36 The results of the present study identified varicella zoster virus as a risk factor for COPD, asthma, and reduced lung function, and therefore varicella zoster virus vaccination may be appropriate for all immunocompromised populations. Importantly, the increased economic burden and morbidity of herpes zoster, its prevention should be done with the use of vaccines, but specific vaccinations are not discussed in this study.

Unlike the present study, risk factors for immune responses to varicella zoster virus-induced immune responses were observed only in asthma. A meta-analysis reports that patients with asthma and COPD have a 24% and 41% increased risk of developing herpes zoster, respectively, compared with healthy controls.37 Asthma and COPD increase the risk of herpes zoster and related complications at any age and may be further elevated in patients receiving inhaled corticosteroids. A Spanish observational, retrospective, non-interventional study found that the most prevalent respiratory diseases (asthma, COPD, and lung cancer) were associated with a higher risk of herpes zoster and postherpetic neuralgia in a standard clinical practice setting.38 The mechanism may be a Th1/Th2 immune imbalance. Similarly, many studies have described an association between asthma and herpes zoster,39 which is consistent with our study. It has been shown that COPD is also a risk factor for herpes zoster virus infection and that the incidence of acute exacerbations of COPD is higher before and after herpes zoster virus attacks.40 However, the results of this study did not find COPD as a risk factor for herpes zoster virus infection.

This study still has some limitations. Firstly, the data for this study were all derived from European populations and further validation with more data from other ethnic populations is still needed to improve the extrapolation of the findings. Second, this study did not find statistical differences in the sensitivity analyses of multiple causal relationships, which may have led to false-positive results in this study. Third, lung diseases usually include both airway and alveolar lesions, and the lack of alveolar-related GWAS data led this study to focus primarily on lung diseases with airway alterations.

Conclusions

This study provides a comprehensive assessment of the causal relationship between antibody-mediated immune responses phenotypes and COPD, asthma, and lung function, and that the relationship is not confounded by other factors. Seven antibody-mediated immune responses phenotypes that may simultaneously have a causal effect on COPD, asthma, FEV1, and FVC, and anti-polyomavirus 2 IgG seropositivity, and varicella zoster virus glycoproteins E and I antibody levels had a clear role. The results of this study can help to screen people at risk of COPD, asthma, and for early prevention and preemptive diagnosis of pre-COPD, asthma condition.

CRediT Authorship Contribution Statement

Drs. LDH and LFR had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. Concept and design: LDH and LFR. Acquisition, analysis, or interpretation of data: XGX, LYL. Drafting of the manuscript: XGX. Statistical analysis: XGX. Obtained funding: LDH and LFR. Supervision: FHY. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Innovation Team and Talents Cultivation Program of the National Administration of Traditional Chinese Medicine (Grant No. ZYYCXTD-D-202003); Natural Science Foundation of China (No: 82374602). The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Declaration of Competing Interest

The authors declare that they have no competing interests.

Acknowledgements

We are grateful for the public sharing of GWAS summary data for COPD, asthma, lung function, and antibody-mediated immune responses, which made this work possible.

Data Availability

The data used in this study can be obtained from GWAS Catalog (https://www.ebi.ac.uk/gwas/; antibody-mediated immune responses GWAS ID: GCST90006909 to GCST90006884; COPD: GCST90044074; emphysema: GCST90044073; bronchiectasis: GCST90044075; asthma: GCST90044072; FEV1: GCST90292609; FVC: GCST90292610; FEV1/FVC: GCST90292611), and the code scripts used in this study are available upon reasonable request from the corresponding author.

Appendix B
Supplementary Data

The following are the supplementary data to this article:

References
[1]
M.I. Asher, L. García-Marcos, N.E. Pearce, D.P. Strachan.
Trends in worldwide asthma prevalence.
Eur Respir J, 56 (2020), pp. 2002094
[2]
P. Venkatesan.
GOLD COPD report: 2023 update.
Lancet Respir Med, 11 (2023), pp. 18
[3]
GBD 2019 Diseases and Injuries Collaborators.
Global burden of 369 diseases, injuries in 204 countries, territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019.
Lancet, 396 (2020), pp. 1204-1222
[4]
F. Schleich, N. Bougard, C. Moermans, M. Sabbe, R. Louis.
Cytokine-targeted therapies for asthma and COPD.
Eur Respir Rev, 32 (2023), pp. 220193
[5]
K.R. Bracke, F. Polverino.
Blunted adaptive immune responses and acute exacerbations of COPD: breaking the code.
Eur Respir J, 62 (2023), pp. 2301030
[6]
P.J. Barnes.
Inflammatory mechanisms in patients with chronic obstructive pulmonary disease.
J Allergy Clin Immunol, 138 (2016), pp. 16-27
[7]
H. Mikkelsen, E.M. Landt, M. Benn, B.G. Nordestgaard, M. Dahl.
Causal risk factors for asthma in Mendelian randomization studies: a systematic review and meta-analysis.
Clin Transl Allergy, 12 (2022), pp. e12207
[8]
Z. Cao, T. Wu, Y. Fang, F. Sun, H. Ding, L. Zhao, et al.
Dissecting causal relationships between immune cells, plasma metabolites, and COPD: a mediating Mendelian randomization study.
Front Immunol, 15 (2024), pp. 1406234
[9]
S.M. O’Connor, C.E. Taylor, J.M. Hughes.
Emerging infectious determinants of chronic diseases.
Emerg Infect Dis, 12 (2006), pp. 1051
[10]
G. Butler-Laporte, D. Kreuzer, T. Nakanishi, A. Harroud, V. Forgetta, J.B. Richards.
Genetic determinants of antibody-mediated immune responses to infectious diseases agents: a genome-wide and HLA association study. Open forum infectious diseases.
Oxford University Press, (2020),
[11]
C.A. Emdin, A.V. Khera, S. Kathiresan.
Mendelian randomization.
JAMA, 318 (2017), pp. 1925-1926
[12]
V.W. Skrivankova, R.C. Richmond, B.A.R. Woolf, N.M. Davies, S.A. Swanson, T.J. VanderWeele, et al.
Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration.
Br Med J, 375 (2021), pp. n2233
[13]
G. Butler-Laporte, D. Kreuzer, T. Nakanishi, A. Harroud, V. Forgetta, J.B. Richards.
Genetic determinants of antibody-mediated immune responses to infectious diseases agents: a genome-wide and HLA Association Study.
Open Forum Infect Dis, 7 (2020),
[14]
L. Jiang, Z. Zheng, H. Fang, J. Yang.
A generalized linear mixed model association tool for biobank-scale data.
Nat Genet, 53 (2021), pp. 1616-1621
[15]
N. Shrine, A.G. Izquierdo, J. Chen, R. Packer, R.J. Hall, A.L. Guyatt, et al.
Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk.
Nat Genet, 55 (2023), pp. 410-422
[16]
S. Sanna, N.R. van Zuydam, A. Mahajan, A. Kurilshikov, A. Vich Vila, U. Võsa, et al.
Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases.
Nat Genet, 51 (2019), pp. 600-605
[17]
M. Verbanck, C.Y. Chen, B. Neale, R. Do.
Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases.
Nat Genet, 50 (2018), pp. 693-698
[18]
B.L. Pierce, H. Ahsan, T.J. Vanderweele.
Power and instrument strength requirements for Mendelian randomization studies using multiple genetic variants.
Int J Epidemiol, 40 (2011), pp. 740-752
[19]
S. Burgess, A. Butterworth, S.G. Thompson.
Mendelian randomization analysis with multiple genetic variants using summarized data.
Genet Epidemiol, 37 (2013), pp. 658-665
[20]
G. Hemani, J. Zheng, B. Elsworth, K.H. Wade, V. Haberland, D. Baird, et al.
The MR-Base platform supports systematic causal inference across the human phenome.
eLife, 7 (2018), pp. e34408
[21]
B.L. Pierce, S. Burgess.
Efficient design for Mendelian randomization studies: subsample and 2-sample instrumental variable estimators.
Am J Epidemiol, 178 (2013), pp. 1177-1184
[22]
J. Bowden, G. Davey Smith, S. Burgess.
Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression.
Int J Epidemiol, 44 (2015), pp. 512-525
[23]
J. Bowden, G. Davey Smith, P.C. Haycock, S. Burgess.
Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator.
Genet Epidemiol, 40 (2016), pp. 304-314
[24]
F.P. Hartwig, G. Davey Smith, J. Bowden.
Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption.
Int J Epidemiol, 46 (2017), pp. 1985-1998
[25]
A. Nalbandian, A.D. Desai, E.Y. Wan.
Post-COVID-19 condition.
[26]
N. Ahsan, K.V. Shah.
Polyomaviruses and human diseases.
Adv Exp Med Biol, 577 (2006), pp. 1-18
[27]
I. Cortese, P. Muranski, Y. Enose-Akahata, S.K. Ha, B. Smith, M. Monaco, et al.
Pembrolizumab treatment for progressive multifocal leukoencephalopathy.
N Engl J Med, 380 (2019), pp. 1597-1605
[28]
A.B. Villaseñor-Altamirano, D. Jain, Y. Jeong, J.A. Menon, M. Kamiya, H. Haider, et al.
Activation of CD8(+) T cells in chronic obstructive pulmonary disease lung.
Am J Respir Crit Care Med, 208 (2023), pp. 1177-1195
[29]
P. Russo, F. Milani, A. De Iure, S. Proietti, D. Limongi, C. Prezioso, et al.
Effect of cigarette smoking on clinical and molecular endpoints in COPD patients.
Int J Mol Sci, 25 (2024), pp. 5834
[30]
J. Bayliss, T. Karasoulos, C.A. McLean.
Immunosuppression increases JC polyomavirus large T antigen DNA load in the brains of patients without progressive multifocal leukoencephalopathy.
J Infect Dis, 207 (2013), pp. 133-136
[31]
S. Gheuens, E. Bord, S. Kesari, D.M. Simpson, R.T. Gandhi, D.B. Clifford, et al.
Role of CD4+ and CD8+ T-cell responses against JC virus in the outcome of patients with progressive multifocal leukoencephalopathy (PML) and PML with immune reconstitution inflammatory syndrome.
J Virol, 85 (2011), pp. 7256-7263
[32]
C.L. Bennett, D. Focosi, M.P. Socal, J.C. Bian, C. Nabhan, W.J. Hrushesky, et al.
Progressive multifocal leukoencephalopathy in patients treated with rituximab: a 20-year review from the Southern Network on Adverse Reactions.
Lancet Haematol, 8 (2021), pp. e593-e604
[33]
F.C. Ringshausen, M. Heckmann, B. Weissbrich, F. Neske, I. Borg, U. Knoop, et al.
No evidence for WU polyomavirus infection in chronic obstructive pulmonary disease.
Infect Agent Cancer, 4 (2009), pp. 12
[34]
L. Ren, R. Gonzalez, Z. Xie, J. Zhang, C. Liu, J. Li, et al.
WU and KI polyomavirus present in the respiratory tract of children, but not in immunocompetent adults.
J Clin Virol, 43 (2008), pp. 330-333
[35]
E. Safonova, B.P. Yawn, T. Welte, C. Wang.
Risk factors for herpes zoster: should people with asthma or COPD be vaccinated?.
Respir Res, 24 (2023), pp. 35
[36]
S. Simon, O. Joean, T. Welte, J. Rademacher.
The role of vaccination in COPD: influenza, SARS-CoV-2, pneumococcus, pertussis, RSV and varicella zoster virus.
Eur Respir Rev, 32 (2023), pp. 230034
[37]
R.P. Insinga, R.F. Itzler, J.M. Pellissier, P. Saddier, A.A. Nikas.
The incidence of herpes zoster in a United States administrative database.
J Gen Intern Med, 20 (2005), pp. 748-753
[38]
D. Morena, S. Lumbreras, J.M. Rodríguez, C. Campos, M. Castillo, M. Benavent, et al.
Chronic respiratory diseases as a risk factor for herpes zoster infection.
Arch Bronconeumol, 59 (2023), pp. 797-804
[39]
F. Marra, K. Parhar, B. Huang, N. Vadlamudi.
Risk Factors for herpes zoster infection: a meta-analysis.
Open Forum Infect Dis, 7 (2020),
[40]
P. Ghaswalla, P. Thompson-Leduc, W.Y. Cheng, C. Kunzweiler, M.J. Wang, M. Bogart, et al.
Increased health care resource utilization and costs associated with herpes zoster among patients aged ≥50 years with chronic obstructive pulmonary disease in the United States.
Chronic Obstr Pulm Dis, 8 (2021), pp. 502-516
Copyright © 2024. The Authors
Archivos de Bronconeumología
Article options
Tools

Are you a health professional able to prescribe or dispense drugs?