PGAP3 regulates human bronchial epithelial cell mRNAs present in asthma and respiratory virus reference data sets

Eric Leslie, Marina Miller, Allison Lafuze, Sofya Svyatskaya, Gil-Soon Choi, David H. Broide
doi: https://doi.org/10.1101/2024.07.03.24309917

Abstract

PGAP3 is a glycosylphosphatidylinositol (GPI) phospholipase gene localized within chromosome 17q12-21, a region highly linked to asthma. Although much is known about the function of other chromosome 17q12-21 genes expressed at increased levels in bronchial epithelium such as ORMDL3 and GSDMB, little is known about the function of increased PGAP3 expression in bronchial epithelium in the context of asthma. The aim of this study was therefore to determine whether increased PGAP3 expression in human bronchial epithelial cells regulated expression of mRNA pathways important to the pathogenesis of asthma by utilizing RNA-sequencing and bioinformatic analysis. We performed RNA-sequencing on normal human bronchial epithelial cells transfected with PGAP3 for 24 and 48 hours. PGAP3 regulated genes were compared to asthma and respiratory virus (influenza A, rhinovirus, respiratory syncytial virus) reference data sets to identify PGAP3 target genes and pathways. Approximately 9% of the upregulated PGAP3-induced genes were found in an asthma reference data set, 41% in a rhinovirus reference data set, 33% in an influenza A reference data set, and 3% in a respiratory syncytial virus reference data set. PGAP3 significantly upregulated the expression of several genes associated with the innate immune response and viral signatures of respiratory viruses associated with asthma exacerbations. Two of the highest expressed genes induced by PGAP3 are RSAD2, OASL, and IFN-λ, which are anti-viral genes associated with asthma. PGAP3 also upregulated the antiviral gene BST2, which like PGAP3 is a GPI-anchored protein. We conclude that PGAP3 expression in human bronchial epithelial cells regulates expression of genes known to be linked to asthma, and also regulates the bronchial epithelial expression of genes pertinent to the pathogenesis of respiratory viral triggered asthma exacerbations.

Introduction

Chromosome 17q12-21 is highly linked to asthma as has been observed in multiple genetic epidemiologic studies1–3. The biologic function of two of the nine genes located in chromosome 17q12-21 (ORMDL3 and GSDMB)2–5 have been well studied and both genes have been found to promote airway remodeling, which could contribute to the development of asthma. The chr:17q12-21 region contains proximal (PGAP3-ERBB2), core (IKZF3-ZPBP2-GSDMB-ORMDL3), and distal (GSDMA) regions, which have been suggested as independent regions associated with asthma2. Studies of the relationship of PGAP3 SNPs and asthma have predominantly focused on three PGAP3 SNPs, namely rs2517954, rs2517955, and rs29415046–11, which are subsequently abbreviated 954, 955, and 504. Several studies have demonstrated that these three PGAP3 SNP risk alleles are associated with increased asthma susceptibility, increased asthma severity, and increased longitudinal asthma exacerbations6–8. SNP 955 in PGAP3 is significantly associated with longitudinal asthma exacerbations after multiple-test adjustment10. In addition, two SNPs in PGAP3 (955, 504) have risk alleles associated with increased PGAP3 expression10. eQTL studies of PGAP3 demonstrated that PGAP3 SNP 954 has an independent association with increased PGAP3 lung expression11, while PGAP3 SNP 955 is a strong eQTL for PGAP3 in lung tissue11. Although there is thus clear evidence from eQTL studies11 as well as genetic association studies6–8 that the risk alleles of PGAP3 associate with PGAP3 expression, asthma susceptibility, asthma severity, and asthma exacerbations, there are currently no published studies of the biology of PGAP3 in asthma to explain why there is this association, which is the focus of the current study.

PGAP3 (Post-GPI Attachment to Proteins 3) (also known as PERLD19 and PER112) is a chr:17q12-21 localized gene that encodes a glycosylphosphatidylinositol (GPI)-specific phospholipase that primarily localizes to the Golgi apparatus12, but can also be detected in the plasma membrane and cytosol (data from Human Protein Atlas v23.proteinatlas.org: https://www.proteinatlas.org/ENSG00000161395-PGAP3)13. PGAP3 encodes a seven-transmembrane protein in the Golgi that removes unsaturated fatty acids from the sn-2 position of GlycosylPhosphatidyliInositol (GPI)12. GPI anchoring is a posttranslational modification that tethers proteins to plasma membranes, and thus plays a key role in protein sorting and trafficking12. There are over 150 human GPI-anchored proteins (GPI-APs), including receptors, adhesion molecules, and enzymes12, which underscores the importance of the genes including PGAP3 known to be involved in the biosynthesis and remodeling of the GPI anchor. The GPI anchor is formed in the ER, transported to the Golgi for fatty acid remodeling and cellular export to lipid rafts in the cell membrane12,13. The PGAP3 SNP 504 linked to asthma9 has also been associated with peripheral blood mononuclear cell surface and secreted GPI-AP levels9, underscoring a potential mechanistic link between PGAP3, GPI-AP, and lipid rafts. Although there are no studies of the function of PGAP3 related to asthma, PGAP3 is known to regulate lipid rafts14 that contain GPI-AP, receptors, adhesion molecules, signaling molecules, and enzymes of potential importance to asthma.

As increased expression of other chr 17q12-21 genes (i.e. ORMDL3, GSDMB)2–5 in bronchial epithelial cells has identified genes and pathways important to the pathogenesis of asthma, the aim of this study was to determine whether increased PGAP3 expression in normal human bronchial epithelial cells (NHBE) regulated expression of mRNA pathways important to the pathogenesis of asthma by utilizing RNA-sequencing and bioinformatic analysis.

Methods

Normal Human Bronchial Epithelial Cell Cultures (NHBE)

NHBE (ATCC) were cultured in 75 mL flasks with 25 mL epithelial-specific media comprised of 98% bronchial epithelial cell media, 1% epithelial cell growth serum, and 1% 10,000 U/mL penicillin (ScienCell). For experiments, NHBE were passaged and subdivided for continued stock culturing (ζ 1x106 cells) or transferred to standard 6-well plates (1x106 per well in 2 mL media) for 24 hours at 37°C prior to experiments. Experiments were completed within 3 – 7 passages.

PGAP3 Transfection of NHBE

To study the effects on biological pathways of increased PGAP3 expression in NHBE (which do not express significant levels of PGAP3), plasmid PGAP3 (PGAP3) was transfected into NHBE using similar methods as has previously been used to study the biology of increased ORMDL3 or increased GSDMB expression in NHBE15,16. Transfection solutions were prepared with an expression vector containing PGAP3 plasmid or with only the expression vector (empty vector) obtained from Origene (PGAP3 Human Untagged Clone, Accession No. NM_033419, vector: pCMV6-XL4, https://www.origene.com), lipofectamine (Invitrogen), and serum-free media (Opti-MEM, Cat No. 11058021, Invitrogen). Solutions were prepared with the following final volumes per well: 1 μL plasmid, 15 μL lipofectamine, and 234 μL DMEM media. After solutions were incubated for 5 minutes, epithelial-specific media was aspirated, cells were washed with 1x PBS (Gibco), and 250 μL transfection solution was added to each well. An additional 750 μL DMEM media was added to each well for a final volume of 1 mL per well. NHBE were transfected for 24 hours and 48 hours prior to RNA-sequencing.

NHBE RNA Extraction and cDNA Synthesis

After transfection with PGAP3, RNA was extracted from NHBE using a RNeasy Mini Kit (50, Cat No. 74104, Qiagen). Following the manufacturer’s instructions, cells were washed with PBS prior to adding 350 mL lysate buffer (Buffer RLT) to each well. After obtaining RNA diluted in RNase-free water, 2 μL of sample were used to measure RNA quantity and quality using a NanoDrop 2000/2000c spectrophotometer (ND2000, Thermo Fisher Scientific). Samples were then diluted in RNase-free water to standardize concentrations of 1 μg RNA in a final volume of 20 μL, added to cDNA EcoDry Premix (Cat No. 639543, TaKaRa) and reverse transcribed using an Eppendorf Mastercycler (MilliporeSigma) set to 42°C for 60 min and 70°C for 10 min.

RNA-Sequencing and Statistical Analysis of PGAP3-transfected NHBE

Confirmation of PGAP3 transfection was assessed via quantitative polymerase chain reaction (qPCR, PGAP3 primer Assay ID: Hs.PT.58.19820740, IDT). Three PGAP3-transfected and three control samples were sent to GENEWIZ (Azenta Life Sciences) for RNA-sequencing where differential gene expression was obtained from their automated workflow. In brief, sequence reads were assessed for quality, trimmed using Trimmomatic v.0.36, and mapped to the human reference genome GRCh38 using the STAR aligner v.2.5.2b. Unique gene hit counts were calculated using the Subread package v.1.5.2 and filtered for reads that fell within exon regions. Differential expression analysis was performed using DESeq2. Significance and log2 fold change (log2FC) were determined via the Wald test.

Bioinformatic Analysis of PGAP3-Regulated Genes and Comparison to Reference Data Sets

We performed bioinformatic and statistical analysis on the following data sets:

1. PGAP3 transfected NHBE vs Empty Vector transfected NHBE (described above).

2. PGAP3 transfected NHBE vs Asthma reference data set

Differentially expressed genes (p < 0.05) were analyzed by overrepresentation, gene set enrichment, and network topology-based analysis using WebGestalt (https://www.webgestalt.org) (Liao 2019). For overrepresentation analysis, an asthma reference data set was used that included genes associated with 61 asthma-specific loci discovered by Pividori et al. (2019)17 in a genome-wide association study of the UK BioBank data from 37,846 self-reported asthma patients. Specifically, we identified protein-coding genes associated with the 61 asthma loci from a search in National Center for Biotechnology Information (NCBI) Gene Database (https://www.ncbi.nlm.nih.gov/gene/?term=) on March 28, 2023 searching with the 61 asthma loci and then with the advanced search function ‘Default Map Location’ and results filtered for “current” and “protein-coding.” This search identified 2326 protein-coding genes, which we used as our asthma reference data set in our analysis. For network topology analysis, the protein-protein interaction network from BIOGRID was used. Significant Gene Ontology (GO) biological processes of non-redundant genes and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were obtained from overrepresentation and gene set enrichment analyses while seed genes were obtained from network topology analysis. Inclusion parameters used for analysis (as applicable) included 5-2000 genes per category, false discovery rate significance level defined as < 0.05, and Benjamin-Hochberg multiple test adjustment.

3. PGAP3 transfected NHBE vs Respiratory Virus reference data sets

Based on bioinformatics analyses, differentially expressed genes were also compared to other RNA-seq data sets regarding NHBE infected with influenza A (Gao et al. 2021)18, rhinovirus (Helling et al. 2020)19, or respiratory syncytial virus (Mayer et al. 2007)20 to compare PGAP3-induced genes with respiratory viruses known to cause exacerbations of asthma. The influenza A virus reference data set acquired from Gao et al. (2021)18 contains significantly differentially expressed genes (FDR < 0.05), determined by single-cell RNA-sequencing, of human nasal upper airway epithelial cells collected from adult patients with active influenza A virus infection. While nasal epithelium may not directly reflect the expression of specific genes within bronchial epithelium, nasal samples have previously been studied to determine the effects of 17q21 risk alleles21 and in another study the functions of differentially expressed genes in asthmatic bronchial and nasal epithelium were shown to be similar22. Therefore, the reference data set from Gao et al. (2021) was included for further analysis regarding any mechanistic changes due to increased PGAP3 expression. The rhinovirus reference data set acquired from Helling et al. (2020)19 contains differentially expressed genes (FDR < 0.05), determined by RNA-sequencing, of human asthmatic vs. vehicle-treated bronchial epithelial cells infected with rhinovirus for 24 hours. The respiratory syncytial virus reference data set acquired from Mayer et al. (2007)17 contains significantly regulated genes, determined by microarray analysis, of human BEAS-2B cells (an immortalized human bronchial epithelial cell line) infected with respiratory syncytial virus for 24 hours.

GPI-Anchored Proteins among Significantly Differentially Expressed Genes

Given PGAP3 being a GPI protein and its role in GPI-AP remodeling, GPI-APs were searched for among differentially expressed genes (p < 0.05). For this, a reference list of GPI-APs was obtained from UniProt (https://www.uniprot.org)23 on August 23, 2023. This list was obtained by searching for “GPI-anchored protein” and filtering for Human, resulting in a reference list of 168 GPI-APs.

Data Availability Statement

The RNA-sequencing data set generated and analyzed in the current study are available on Figshare: https://doi.org/10.6084/m9.figshare.24666630.v1.

Results

Increased PGAP3 upregulated multiple genes in NHBE

In total, 654 genes at 24 hrs and 308 genes at 48 hrs were upregulated (p < 0.05) after NHBE were transfected with PGAP3 compared to empty vector control (Table 1). Unidentified genes were denoted with “NA” and removed from further analysis.

View this table:
Table 1. PGAP3-Regulated Genes in Normal Human Bronchial Epithelial Cells

Comparison of PGAP3-upregulated genes in NHBE to reference asthma data set

To better understand the relationship between PGAP3-induced differential gene expression to asthma and to respiratory viruses that promote asthma exacerbations and/or the development of asthma (rhinovirus, influenza A, respiratory syncytial virus), PGAP3 upregulated genes were compared to reference data sets. In the UK BioBank asthma reference data set17, 62 PGAP3 upregulated genes expressed at 24 hrs were detected, and 36 PGAP3 upregulated genes expressed at 48 hrs were detected (Table 2).

View this table:
Table 2. PGAP3-Regulated Genes Found in Reference Data Sets

Comparison of PGAP3 Upregulated genes in NHBE to reference respiratory virus data sets

In comparing the PGAP3 upregulated genes in NHBE to the respiratory virus reference data sets, a large number of PGAP3 upregulated genes were found in the influenza A reference data set (217 genes at 24 hrs and 147 genes, as well as the rhinovirus reference data set (267 genes at 24 hrs and 163 genes at 48 hrs) (Table 2). Relatively fewer PGAP3 regulated genes were also found in the respiratory syncytial virus reference data set (22 genes at 24 hrs and 9 genes at 48 hrs) (Table 2).

Top 10 PGAP3 Upregulated and Top 10 PGAP3 Downregulated genes in NHBE

Tables 3 and 4 lists the top 10 upregulated and top 10 downregulated differentially expressed (p < 0.05) genes in PGAP3 transfected NHBE and provides information regarding in which reference data sets these PGAP3 induced genes were also found. The most notable findings are found in the top 10 PGAP3 upregulated genes where several interferon-stimulating genes (IFIT1, IFI6, IFI27, and IFI44L), RSAD2 (which encodes the antiviral IFN-induced protein viperin24, and the chemokine CXCL11 are also found in both the influenza A and rhinovirus reference data sets. This suggests that some of the most upregulated genes induced by PGAP3 are involved with the innate immune response to respiratory viruses with predominately antiviral activity. Other notable genes induced by PGAP3 include OASL, which was detected in three reference data sets (asthma, influenza A, rhinovirus), IFNL1 (also called IFN-λ), a prominent antiviral gene expressed after rhinovirus infection19, and BST2, which is the only GPI-AP found in the top 10 upregulated genes at 24 or 48 hours and is also found in the influenza A and rhinovirus reference data sets. Among the most downregulated genes, there are several indications of dysregulated calcium signaling through lower expression of genes involved with calcium binding (PAMR1) and calcium-dependent cell adhesion (CDH15), along with other adhesion molecules (CD36, THY1), and HLA-DQA1, a well-known Class II antigen presenting molecule associated with asthma25.

View this table:
Table 3. Top 10 PGAP3-Upregulated Genes in NHBE
View this table:
Table 4. Top 10 PGAP3-Downregulated Genes in NHBE

Top PGAP3 Upregulated genes found in the asthma and respiratory viral reference data sets

Given a greater number of genes differentially expressed by PGAP3 vs empty vector control at 24 hours and an upregulation of the most significant pathways (discussed in the bioinformatics analysis results below), a closer analysis of PGAP3-induced upregulated genes at 24 hours was used to find connections with the asthma and respiratory viral reference data sets. Tables 5-7 provide the list of top PGAP3-induced upregulated genes (p < 0.05) found in the asthma and respiratory viral reference data sets. These tables highlight several additional genes commonly associated with asthma. PGAP3 regulated genes found in the UK BioBank asthma reference data set include the chemokine CCL5 (also known as RANTES), the STAT inhibitor SOCS1, and the IL-15 receptor IL15RA (Table 5). Several additional PGAP3 regulated antiviral genes are shared between the influenza A (Table 6) and rhinovirus (Table 7) reference data sets, including MX1, MX2, OAS1, and OAS2. The respiratory syncytial virus reference data set shares the following four PGAP3 regulated genes with the other data sets of potential importance to asthma: DDX58 (also known as RIG-I, found in the rhinovirus and influenza A reference data sets), CXCL3 (found in the influenza A reference data set), SOCS3 (found in the influenza A reference data set), and IL6 (found in the asthma and rhinovirus reference data sets) (Table 8).

View this table:
Table 5. PGAP3-Upregulated Genes Found in the UK BioBank Asthma Reference Data Set
View this table:
Table 6. Top 25 PGAP3-Upregulated Genes Found in Influenza A Reference Data Set
View this table:
Table 7. Top 25 PGAP3-Upregulated Genes Found in Rhinovirus Reference Data Set
View this table:
Table 8. PGAP3-Upregulated Genes Found in Respiratory Syncytial Virus Reference Data Set

PGAP3 upregulated genes found in more than one respiratory viral reference data set

Table 9 provides the top 25 PGAP3 upregulated genes also found in the influenza A, rhinovirus, and respiratory syncytial virus reference data sets. In addition, Figure 1 provides a diagram to illustrate the overlap of PGAP3-differentially expressed genes found in the three individual respiratory viral reference data sets (rhinovirus, influenza A, respiratory syncytial virus). When analyzing which genes are unique or shared between the data sets, the largest proportion consists of shared genes between the rhinovirus and influenza A reference data sets (43%) while < 2% of the genes were shared between rhinovirus and respiratory syncytial virus reference data sets and between influenza A and respiratory syncytial virus reference data sets. The rhinovirus reference data set had the largest proportion of unique genes (32%), followed by influenza A (18%) and respiratory syncytial virus (< 2%). Seven genes (< 2%) were shared between all three respiratory virus reference data sets: IFIT1, IFI44, DDX58, PMAIP1, IRF1, FAM46A, and PLAUR.

Figure 1.
Figure 1.

Venn diagram of unique and shared PGAP3 regulated genes found in the three reference data sets of human bronchial epithelial cells infected with respiratory viruses known to cause exacerbations of asthma: rhinovirus, influenza A, and respiratory syncytial virus.

View this table:
Table 9. Top 25 PGAP3-Upregulated Genes Found in Both Influenza A and Rhinovirus Reference Data Sets

PGAP3 regulation of GPI-Anchored Proteins

Table 10 focuses only on the 13 GPI-anchored proteins within the RNA-Seq data set (Table 10A; 6 upregulated: Table 10B; 7 downregulated). There were 6 PGAP3 upregulated GPI-anchored proteins (p < 0.05) at 24 hours and 4 PGAP3 upregulated GPI-anchored proteins (p < 0.05) at 48 hours. The upregulated GPI-anchored proteins BST2 (7.14 log2 fold change) and PGAP3 had noticeably higher fold changes compared to the other GPI-anchored proteins. All 6 GPI-anchored protein genes upregulated by PGAP3 were found in at least one reference data set (1 in the asthma reference data set, 4 in the influenza A virus reference data set, 5 in the rhinovirus reference data set, and 1 in the respiratory syncytial virus reference data set). The upregulation of the GPI-anchored protein LY6E at both 24 and 48 hours is also notable given its role in enhancing influenza A infection and other viruses, via entry into the cell, in an evolutionarily conserved manner26. Overall, these results demonstrate that PGAP3 influenced expression of a subset GPI-APs, in particular BST2 and LY6E.

View this table:
Table 10. PGAP3-Upregulated GPI-Anchored Proteins

PGAP3 regulated pathway analysis by GSEA

Results of the gene set enrichment analysis (GSEA) include 27 positively regulated and 0 negatively regulated GO biological processes at 24 hours while there were 14 positively regulated and 5 negatively regulated GO biological processes at 48 hours. Regarding KEGG pathways, there were multiple positively regulated pathways (9 at 24 hrs; 8 at 48 hours), with no negatively regulated pathways. The top 10 GO biological processes and KEGG pathways based on enrichment score are provided in Figure 2. Significant KEGG pathways (Figure 2A) generally demonstrated an upregulation of innate immune responses in defense against viral RNA and peptide fragments of infected cells (Cytosolic DNA-sensing as well as RIG-I, NOD-like, and Toll-like receptor signaling) as well as collective changes in gene expression common to several viruses. A key viral pathway detected is influenza A, which is known to promote wheezing in asthma patients or those predisposed to develop asthma27. The five downregulated pathways are all involved in the cell cycle (cell cycle G2/M phase transition, regulation of cell cycle phase transition, cytokinesis, mitotic cell cycle phase transition, and chromosome segregation). The GO biological process results (Figure 2B) are consistent with the significant KEGG pathways by their associations with innate immunity that includes, but is not limited to, type 1 interferon and cellular defense responses. When the UK BioBank asthma reference data set was used as a background in the overrepresentation analysis (Figure 3), the top KEGG pathways (Figure 3A) and GO biological processes (Figure 3B) are overall consistent with the GSEA since many are involved with viral signatures and the innate immune system.

Figure 2.
Figure 2. Gene Set Enrichment Analysis.

Illustrated are the top 10 KEGG pathways (2A) and Gene Ontology (GO) biological processes (2B) by enrichment score in human bronchial epithelial cells transfected with PGAP3 for 24 and 48 hours.

Figure 3.
Figure 3. Overrepresentation Analysis.

Illustrated are the top 10 KEGG pathways (3A) and Gene Ontology (GO) biological processes (3B) by enrichment score in human bronchial epithelial cells transfected with PGAP3 for 24 and 48 hours. The UK BioBank asthma reference data set was used, which included 2,326 genes associated with 61 asthma-specific loci identified by Pividori et al. (2019)17.

PGAP3 regulated Pathway Analysis by Highlighted Seeds and Network Topology Analysis

Highlighted seeds with a log2FC ≥ 2.0 were included to identify the most connected or potentially influential seeds (genes) from the network topology analysis (Table 11). TRIM25, ISG15, and HLA-C were three of the four highlighted seeds that met these criteria. They are potentially involved in novel pathways that we identified that PGAP3 regulates including RIG-1 signaling and antigen presentation (TRIM25 and ISG15 are involved with RIG-I signaling and HLA-C is involved with antigen presentation).

View this table:
Table 11. Network Topology Analysis Highlighted Seeds

Discussion

As asthma is a disease of the airways the bronchial epithelium plays an important role in the pathogenesis of asthma as inhaled allergens and viruses can trigger the airway epithelium to express pathways important to the pathogenesis of asthma28. Increased expression of chromosome 17q12-21 genes such as ORMDL3 and GSDMB have been detected in the bronchial epithelium in asthma and are associated with expression of pathways pertinent to the pathogenesis of asthma16,29. In this study, we have identified that increased expression of another chromosome 17q12-21 gene, namely PGAP3, upregulates expression of 654 genes in bronchial epithelium. We compared the genes induced by PGAP3 to an asthma reference data set based on asthma-specific loci discovered from the UK BioBank to determine whether any of the genes induced by PGAP3 were also present in this asthma data set. Approximately 9% of the upregulated genes induced by PGAP3 in NHBE (62 out of 654 genes upregulated by PGAP3) were identified in the asthma reference data set, which further supports PGAP3’s link to asthma. Interestingly, as respiratory viruses are also linked to asthma exacerbations and/or the development of asthma we identified that approximately 41% of the upregulated PGAP3-induced genes in epithelium were found in a rhinovirus reference data set, 33% in an influenza A reference data set, and 3% in a respiratory syncytial virus reference data set. PGAP3 significantly upregulated the expression of several genes associated with the innate immune response and viral signatures of respiratory viruses associated with asthma exacerbations. Two of the highest expressed genes induced by PGAP3 are RSAD230 and OASL31, which are anti-viral genes associated with asthma. PGAP3 also upregulated the antiviral gene BST232, which like PGAP3 is a GPI-anchored protein. These studies suggest that PGAP3 expression in human bronchial epithelial cells regulates expression of genes known to be linked to asthma, and that PGAP3 also regulates the bronchial epithelial responses to viruses pertinent to the pathogenesis of respiratory viral triggered asthma exacerbations. Interestingly, the pathways induced by PGAP3 in bronchial epithelium have differences from those induced by other chromosome 17q12-21 genes linked to asthma such as ORMDL3 or GSDMB, but share with ORMDL3 and GSDMB the ability to induce a few selected genes such as OAS4,16.

The 62 genes upregulated by PGAP3 we identified as also being present in the UK BioBank asthma reference data set (Table 5), include the chemokine CCL5 (also known as RANTES), the STAT inhibitor SOCS1, the IL-15 receptor IL15RA. CCL5 is localized to the same chromosome loci as PGAP3 (17q12-21) and the CCL5 -403G/A polymorphism has been linked to asthma33. Recently, high mRNA expression of CCL5 has been correlated with the Th1 cytokines CXCL9 and CXCL10 (also known as IP-10) in sputum from patients with asthma34. In mouse models of asthma there is a link between CCL5 and tissue-resident memory T-cell activation as inhibition of CCR5 (the CCL5 receptor) blunts the effect34. Several chemokines in addition to CCL5 were upregulated by PGAP3 including CXCL9, and CXCL10, which may indicate expression of PGAP3 promotes recruitment of inflammatory cells into the airway. Another PGAP3 upregulated gene of note is IL15RA, which is a receptor that binds IL-15 and is expressed in the sputum of patients with asthma35. IL-15 has also been shown to promote Th2 differentiation and allergic sensitization in a mouse model of asthma36. Finally, induction of SOCS1 by PGAP3 may also be important to asthma as RNA-sequencing analysis of human nasal epithelial cells from children with asthma found significant expression of SOCS137 that is known to regulate cytokine signaling in both human bronchial epithelial cells from patients with asthma38 and in the lungs of a mouse model of asthma39.

Respiratory viruses are important triggers of asthma exacerbations40. In addition, RV has been shown to increase the risk of the development of childhood asthma in subjects with a SNP linked to chromosome 17q12-2141. Our current study demonstrates that increased expression of PGAP3 in NHBE through induction of antiviral pathways could potentially influence the duration and/or severity of viral infection in epithelium through several mechanisms. For example, PGAP3 expression in NHBE induced robust expression of antiviral genes commonly associated with asthma (RSAD2, OASL, OAS1, OAS2, OAS3, MX1, and MX2), and anti-viral pathways (Cytosolic DNA sensing, RIG-I signaling, NOD-like receptor signaling, and Toll-like receptor signaling). Several other genes including RSAD2, MX1, and BST2 along with ISG15 as highlighted by the network topology analysis have been previously linked to anti-viral responses in epithelial cells 42. Induction of these genes and innate pathways by PGAP3 could all potentially limit the extent of a viral infection in epithelium. The connections from the RNA-Seq analysis to RIG-I signaling is particularly notable. RIG-I signaling is known to induce the antiviral genes OASL and the GPI-AP gene BST243,44, which are two of the highest expressed upregulated genes in the PGAP3 data set. OASL is also known to further increase RIG-I signaling, eliciting a positive feedback loop, to improve the antiviral response31. While not found in the pathway analysis, IFN-λ was one of the top 10 upregulated genes (FC = 7.43) and has been shown to be positively associated with peak levels of RV infection, despite asthmatic epithelial cells showing a delayed IFN-λ response to rhinovirus45. However, since a delayed IFN-λ response to rhinovirus is not always shown46,47, whether PGAP3 induces IFN-λ in response, or in this case an anticipated response, to viral infection requires further investigation. The ability of PGAP3 to induce BST2 expression, along with other anti-viral genes, suggests PGAP3 regulates antiviral responses in NHBE. Interestingly, TRIM25, a key RIG-I activator48, had the most gene interactions in the network topology analysis, while RIG-I signaling is the second most enriched KEGG pathway. In addition to these PGAP3 regulated anti-viral responses, a smaller number of PGAP3 regulated pro-viral genes were also detected, including LY6E. Given that LY6E enhances cellular infection of influenza A virus among several others (Chikungunya virus, dengue virus, human immunodeficiency virus, West Nile virus, and yellow fever virus)26, PGAP3 may contribute to increased susceptibility to viral infection at the same time it coordinates cellular antiviral responses to protect against viruses.

The enriched GO biological processes and KEGG pathways obtained from the PGAP3 overrepresentation analysis, which incorporated the UK BioBank asthma reference data set, largely reflect those found in the gene set enrichment analysis. This suggests that many of the cellular effects of PGAP3 in NHBE include genes from asthma-specific loci. The presence of TRIM25 and ISG15 among the highlighted seeds (FC = 2.0+) further demonstrates PGAP3’s potential connection to antiviral responses via RIG-I signaling. The bioinformatic analyses provided supporting evidence linking PGAP3 overexpression in NHBE to asthma and respiratory viral signatures.

Although we have identified that PGAP3 upregulated 654 genes of which 9% have been identified in an asthma reference data set, and that PGAP3 upregulated genes and pathways linked to innate immunity and antiviral response, there are limitations to the study in terms of the need for future investigations to examine the biological consequences of PGAP3-induced upregulated genes and pathways, which are currently only inferred from the RNA-Seq and bioinformatic analysis. Studies of bronchial epithelial cells from asthmatics with SNPs linked to increased PGAP3 expression would also be of interest as compared to studies of NHBE. However, recruitment of asthmatics with SNPs linked to increased PGAP3 for an invasive bronchoscopy would be needed.

In summary, we have identified that increased expression of the Chr17q12-21 gene PGAP3 in bronchial epithelium upregulates expression of 654 genes in bronchial epithelium approximately 9% of which genes were identified in the UK Biobank asthma reference data set, supporting PGAP3’s link to asthma. The PGAP3 regulated genes detected in the UK Biobank asthma reference data set include RANTES (CCL5), and OASL. Interestingly, we also identified that PGAP3 regulated genes were identified in data sets of respiratory viruses linked to asthma (rhinovirus and influenza A). Two of the highest expressed genes induced by PGAP3 are RSAD2 and OASL, which are antiviral genes associated with asthma30,31. PGAP3 also upregulated the antiviral gene BST2, which like PGAP3 is a GPI-anchored protein. These studies suggest that PGAP3 expression in human bronchial epithelial cells regulates expression of genes known to be linked to asthma, and that PGAP3 also regulates the bronchial epithelial responses to viruses pertinent to the pathogenesis of respiratory viral triggered asthma exacerbations. These studies underscore the potential importance of PGAP3 to the biology of asthma.

Author Contributions

Conception or design of the work (DHB, EL, MM), acquisition (EL, MM, AL, SS, GSC), analysis or interpretation of the work (EL, MM, DB), drafting or revising manuscript for submission (EL, DB).

Ethical Approval

No ethical approval was obtained because this study does not use human subjects or animals.

Informed Consent

No informed consent was obtained because this study does not use human subjects or animals.

Funding

DHB is supported by NIH grants AI 070535, AI 107779, and AI 242236. EL is supported by T32 AI 007469 to DHB.

Statement of Competing Interests

Authors did not receive any third-party payments or services that influenced the submitted work.

Footnotes

  • Author affiliations were updated. All information was added for University of California San Diego and Dr. Choi's second affiliation was included.

References

  1. 1.↵
    Moffatt MF, Kabesch M, Liang L, et al. Genetic variants regulating ORMDL3 expression contribute to the risk of childhood asthma. Nature. 2007;48(7152):470-473. doi:10.1038/nature06014
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    Stein MM, Thompson EE, Schoettler N, et al. A decade of research on the 17q12-21 asthma locus: Piecing together the puzzle. J Allergy Clin Immunol. 2018;142(3):749–764.e3. doi:10.1016/j.jaci.2017.12.974
    OpenUrlCrossRefPubMed
  3. 3.↵
    Das S, Miller M, Broide DH. Chromosome 17q21 genes ORMDL3 and GSDMB in asthma and immune diseases. Adv Immunol. 2017;135:1–52. doi:10.1016/BS.AI.2017.06.001
    OpenUrlCrossRef
  4. 4.↵
    Miller M, Tam AB, Cho JY, et al. ORMDL3 is an inducible lung epithelial gene regulating metalloproteases, chemokines, OAS, and ATF6. Proc Natl Acad Sci U S A. 2012;109(41):16648–16653. doi:10.1073/PNAS.1204151109
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    Pham AK, Miller M, Rosenthal P, et al. ORMDL3 expression in ASM regulates hypertrophy, hyperplasia via TPM1 and TPM4, and contractility. JCI Insight. 2021;6(7):e136911. doi:10.1172/jci.insight.136911
    OpenUrlCrossRef
  6. 6.↵
    Anantharaman R, Andiappan AK, Nilkanth PP, Suri BK, Wang DY, Chew FT. Genome-wide association study identifies PERLD1 as asthma candidate gene. BMC Med Genet. 2011;12:170. doi:10.1186/1471-2350-12-170
    OpenUrlCrossRef
  7. 7.
    Halapi E, Gudbjartsson DF, Jonsdottir GM, et al. A sequence variant on 17q21 is associated with age at onset and severity of asthma. E J Hum Genet. 2010;18(8):902–908. doi:10.1038/ejhg.2010.38
    OpenUrlCrossRefPubMed
  8. 8.↵
    Nicodemus-Johnson J, Myers RA, Sakabe NJ, et al. DNA methylation in lung cells is associated with asthma endotypes and genetic risk. JCI Insight. 2016;1(20):e90151. doi:10.1172/JCI.INSIGHT.90151
    OpenUrlCrossRef
  9. 9.↵
    Sio YY, Anantharaman R, Lee SQE, Matta SA, Ng YT, Chew FT. The Asthma-associated PER1-like domain-containing protein 1 (PERLD1) Haplotype Influences Soluble Glycosylphosphatidylinositol Anchor Protein (sGPI-AP) Levels in Serum and Immune Cell Proliferation. Sci Rep. 2020;10(1):715. doi:10.1038/s41598-020-57592-9
    OpenUrlCrossRef
  10. 10.↵
    Li X, Christenson SA, Modena B, et al. Genetic analyses identify GSDMB associated with asthma severity, exacerbations, and antiviral pathways. J Allergy Clin Immunol. 2021;147(3):894–909. doi:10.1016/j.jaci.2020.07.030
    OpenUrlCrossRef
  11. 11.↵
    Ober C, McKennan CG, Magnaye KM, et al. Expression quantitative trait locus fine mapping of the 17q12–21 asthma locus in African American children: a genetic association and gene expression study. Lancet Respir Med. 2020;8(5):482–492. doi:10.1016/S2213-2600(20)30011-4
    OpenUrlCrossRef
  12. 12.↵
    Fujita M, Umemura M, Yoko-o T, Jigami Y. PER1 is required for GPI-phospholipase A2 activity and involved in lipid remodeling of GPI-anchored proteins. Mol Biol Cell. 2006;17(12):5253–5264. doi:10.1091/mbc.e06-08-0715
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    Maeda Y, Tashima Y, Houjou T, et al. Fatty acid remodeling of GPI-anchored proteins is required for their raft association. Mol Biol Cell. 2007;18(4):1497–1506. doi:10.1091/mbc.e06-10-0885
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    Murakami H, Wang Y, Hasuwa H, Maeda Y, Kinoshita T, Murakami Y. Enhanced response of T lymphocytes from Pgap3 knockout mouse: Insight into roles of fatty acid remodeling of GPI anchored proteins. Biochem Biophys Res Commun. 2012;417(4):1235–1241. doi:10.1016/j.bbrc.2011.12.116
    OpenUrlCrossRefPubMed
  15. 15.↵
    Chen J, Miller M, Unno H, Rosenthal P, Sanderson M, Broide DH. Orosomucoid-like 3 (ORMDL3) upregulates airway smooth muscle proliferation, contraction, and Ca2+ oscillations in asthma. J Allergy Clin Immunol. 2018;142(1):207–218.e6. doi:10.1016/j.jaci.2017.08.015
    OpenUrlCrossRef
  16. 16.↵
    Das S, Miller M, Beppu AK, et al. GSDMB induces an asthma phenotype characterized by increased airway responsiveness and remodeling without lung inflammation. Proc Natl Acad Sci U S A. 2016;113(46):13132–13137. doi:10.1073/PNAS.1610433113
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    Pividori M, Schoettler N, Nicolae DL, Ober C, Im HK. Shared and distinct genetic risk factors for childhood-onset and adult-onset asthma: genome-wide and transcriptome-wide studies. Lancet Respir Med. 2019;7(6):509–522. doi:10.1016/S2213-2600(19)30055-4
    OpenUrlCrossRef
  18. 18.↵
    Gao KM, Derr AG, Guo Z, et al. Human nasal wash RNA-Seq reveals distinct cell-specific innate immune responses in influenza versus SARS-CoV-2. JCI Insight. 2021;6(22):e152288. doi:10.1172/JCI.INSIGHT.152288
    OpenUrlCrossRef
  19. 19.↵
    Helling BA, Sobreira DR, Hansen GT, et al. Altered transcriptional and chromatin responses to rhinovirus in bronchial epithelial cells from adults with asthma. Commun Biol. 2020;3(1):678. doi:10.1038/s42003-020-01411-4
    OpenUrlCrossRef
  20. 20.↵
    Mayer AK, Muehmer M, Mages J, et al. Differential recognition of TLR-dependent microbial ligands in human bronchial epithelial cells. J Immunol. 2007;178(5):3134–3142. doi:10.4049/JIMMUNOL.178.5.3134
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    Jakwerth CA, Weckmann M, Illi S, et al. 17q21 variants disturb mucosal host defense in childhood asthma. Germany Am J Respir Crit Care Med. 2024;209(12):1519. doi:10.1164/rccm.202305-0934OC
    OpenUrlCrossRef
  22. 22.↵
    Wagener AH, Zwinderman AH, Luiten S, et al. The impact of allergic rhinitis and asthma on human nasal and bronchial epithelial gene expression. PLoS One. 2013;8(11):e80257. doi:10.1371/journal.pone.0080257
    OpenUrlCrossRef
  23. 23.↵
    UniProt Consortium. UniProt: the Universal Protein Knowledgebase in 2023. Nucleic Acids Res. 2023;51(D1):D523–D531. doi:10.1093/nar/gkac1052
    OpenUrlCrossRefPubMed
  24. 24.↵
    Wiedemann GM, Geary CD, Lau CM, Sun JC. Cutting Edge: STAT1-Mediated Epigenetic Control of Rsad2 Promotes Clonal Expansion of Antiviral NK Cells. J Immunol. 2020;205(1):21–25. doi:10.4049/JIMMUNOL.2000086
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    Schoettler N, Rodríguez E, Weidinger S, Ober C. Advances in asthma and allergic disease genetics: Is bigger always better? J Allergy Clin Immunol. 2019;144(6):1495–1506. doi:10.1016/j.jaci.2019.10.023
    OpenUrlCrossRef
  26. 26.↵
    Mar KB, Rinkenberger NR, Boys IN, et al. LY6E mediates an evolutionarily conserved enhancement of virus infection by targeting a late entry step. Nat Commun. 2018;9(1):3603. doi:10.1038/s41467-018-06000-y
    OpenUrlCrossRefPubMed
  27. 27.↵
    Chang Y, Kang JS, Jung K, et al. OASL1-mediated inhibition of type I IFN reduces influenza A infection-induced airway inflammation by regulating ILC2s. Allergy Asthma Immunol Res. 2022;14(1):99–116. doi:10.4168/aair.2022.14.1.99
    OpenUrlCrossRef
  28. 28.↵
    Akdis CA. Does the epithelial barrier hypothesis explain the increase in allergy, autoimmunity and other chronic conditions? Nat Rev Immunol. 2021;21(11):739–751. doi:10.1038/s41577-021-00538-7
    OpenUrlCrossRef
  29. 29.↵
    Miller M, Tam AB, Mueller JL, et al. Cutting Edge: Targeting epithelial ORMDL3 increases, rather than reduces, airway responsiveness and is associated with increased sphingosine-1-phosphate. J Immunol. 2017;198(8):3017–3022. doi:10.4049/JIMMUNOL.1601848
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    Wang F, Shang Y. Radical S-adenosyl methionine domain-containing 2, a potential target of D-tryptophan in asthma treatment, regulates T helper cell type 1/2 balance. Int Immunopharmacol. 2024;129:111581. doi:10.1016/J.INTIMP.2024.111581
    OpenUrlCrossRef
  31. 31.↵
    Zhu J, Ghosh A, Sarkar SN. OASL-a new player in controlling antiviral innate immunity. Curr Opin Vriol. 2015;12:15–19. doi:10.1016/j.coviro.2015.01.010
    OpenUrlCrossRefPubMed
  32. 32.↵
    Berry KN, Kober DL, Su A, Brett TJ. Limiting respiratory viral infection by targeting antiviral and immunological functions of BST-2/Tetherin: Knowledge and gaps. Bioessays. 2018;40(10):e1800086. doi:10.1002/bies.201800086
    OpenUrlCrossRef
  33. 33.↵
    Zhang YG, Huang J, Zhang J, et al. RANTES gene polymorphisms and asthma risk: A meta-analysis. Arch Med Res. 2010;41(1):50–58. doi:10.1016/J.ARCMED.2010.01.002
    OpenUrlCrossRefPubMed
  34. 34.↵
    Gauthier M, Kale SL, Oriss TB, et al. CCL5 is a potential bridge between type 1 and type 2 inflammation in asthma. J Allergy Clin Immunol. 2023;152(1):94–106.e12. doi:10.1016/J.JACI.2023.02.028
    OpenUrlCrossRef
  35. 35.↵
    Hastie AT, Steele C, Dunaway CW, et al. Complex association patterns for inflammatory mediators in induced sputum from subjects with asthma. Clin Exp Allergy. 2018;48(7):787–797. doi:10.1111/CEA.13129
    OpenUrlCrossRefPubMed
  36. 36.↵
    Rückert R, Brandt K, Braun A, et al. Blocking IL-15 prevents the induction of allergen-specific T cells and allergic inflammation in vivo. J Immunol. 2005;174(9):5507–5515. doi:10.4049/JIMMUNOL.174.9.5507
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    Xu Z, Forno E, Acosta-Pérez E, et al. Differential gene expression in nasal airway epithelium from overweight or obese youth with asthma. Pediatr Allergy Immunol. 2022;33(4):e13776. doi:10.1111/PAI.13776
    OpenUrlCrossRef
  38. 38.↵
    Magnaye KM, Naughton KA, Huffman J, et al. A-to-I editing of miR-200b-3p in airway cells is associated with moderate-to-severe asthma. Eur Respir J. 2021;58(1):2003862. doi:10.1183/13993003.03862-2020
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    Lee C, Kolesnik TB, Caminschi I, et al. Suppressor of cytokine signalling 1 (SOCS1) is a physiological regulator of the asthma response. Clin Exp Allergy. 2009;39(6):897–907. doi:10.1111/J.1365-2222.2009.03217.X
    OpenUrlCrossRefPubMed
  40. 40.↵
    Jartti T, Bønnelykke K, Elenius V, Feleszko W. Role of viruses in asthma. Semin Immunopathol. 2020;42(1):61–74. doi:10.1007/s00281-020-00781-5
    OpenUrlCrossRef
  41. 41.↵
    Çalışkan M, Bochkov YA, Kreiner-Møller E, et al. Rhinovirus wheezing illness and genetic risk of childhood-onset asthma. N Engl J Med. 2013;368(15):1398–1407. doi:10.1056/NEJMOA1211592
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    Brice DC, Figgins E, Yu F, Diamond G. Type I interferon and interferon-stimulated gene expression in oral epithelial cells. Mol Oral Microbiol. 2019;34(6):245–253. doi:10.1111/OMI.12270
    OpenUrlCrossRef
  43. 43.↵
    Zhu J, Zhang Y, Ghosh A, et al. Antiviral activity of human OASL protein is mediated by enhancing signaling of the RIG-I RNA sensor. Immunity. 2014;40(6):936–948. doi:10.1016/j.immuni.2014.05.007
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    Wang SM, Huang KJ, Wang CT. BST2/CD317 counteracts human coronavirus 229E productive infection by tethering virions at the cell surface. Virology. 2013;449:287–296. doi:10.1016/j.virol.2013.11.030
    OpenUrlCrossRef
  45. 45.↵
    Krammer S, Sicorschi Gutu C, Grund JC, Chiriac MT, Zirlik S, Finotto S. Regulation and function of interferon-lambda (IFNλ) and its receptor in asthma. Front Immunol. 2021;12:731807. doi:10.3389/fimmu.2021.731807
    OpenUrlCrossRef
  46. 46.↵
    Bosco A. Novel role of the GSDMB/IFNG axis in childhood asthma. Am J Respir Crit Care Med. 2024;209(8):899–900. doi:10.1164/rccm.202312-2259ED
    OpenUrlCrossRef
  47. 47.↵
    Ritchie AI, Jackson DJ, Edwards MR, Johnston SL. Airway epithelial orchestration of innate immune function in response to virus infection. A focus on asthma. Ann Am Thorac Soc. 2016;13 Suppl 1:S55–S63. doi:10.1513/AnnalsATS.201507-421MG
    OpenUrlCrossRef
  48. 48.↵
    Rehwinkel J, Gack MU. RIG-I-like receptors: their regulation and roles in RNA sensing. Nat Rev Immunol. 2020;20(9):537–551. doi:10.1038/s41577-020-0288-3
    OpenUrlCrossRefPubMed
Back to top
Posted July 06, 2024.
Download PDF
Email
PGAP3 regulates human bronchial epithelial cell mRNAs present in asthma and respiratory virus reference data sets
Eric Leslie, Marina Miller, Allison Lafuze, Sofya Svyatskaya, Gil-Soon Choi, David H. Broide
medRxiv 2024.07.03.24309917; doi: https://doi.org/10.1101/2024.07.03.24309917
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area