- Open Access
Family with sequence similarity 13 member A mediates TGF-β1-induced EMT in small airway epithelium of patients with chronic obstructive pulmonary disease
Respiratory Research volume 22, Article number: 192 (2021)
To explore the role of family with sequence similarity 13 member A (FAM13A) in TGF-β1-induced EMT in the small airway epithelium of patients with chronic obstructive pulmonary disease (COPD).
Small airway wall thickness and protein levels of airway remodeling markers, EMT markers, TGF-β1, and FAM13A were measured in lung tissue samples from COPD and non-COPD patients. The correlations of FAM13A expression with COPD severity and EMT marker expression were evaluated. Gain- and loss-of-function assays were performed to explore the functions of FAM13A in cell proliferation, motility, and TGF-β1-induced EMT marker alterations in human bronchial epithelial cell line BEAS-2B.
Independent of smoking status, lung tissue samples from COPD patients exhibited significantly increased small airway thickness and collagen fiber deposition, along with enhanced protein levels of remodeling markers (collagen I, fibronectin, and MMP-9), mesenchymal markers (α-SMA, vimentin, and N-cadherin), TGF-β1, and FAM13A, compared with those from non-COPD patients. FAM13A expression negatively correlated with FEV1% and PO2 in COPD patients. In small airway epithelium, FAM13A expression negatively correlated with E-cadherin protein levels and positively correlated with vimentin protein levels. In BEAS-2B cells, TGF-β1 dose-dependently upregulated FAM13A protein levels. FAM13A overexpression significantly promoted cell proliferation and motility in BEAS-2B cells, whereas FAM13A silencing showed contrasting results. Furthermore, FAM13A knockdown partially reversed TGF-β1-induced EMT marker protein alterations in BEAS-2B cells.
FAM13A upregulation is associated with TGF-β1-induced EMT in the small airway epithelium of COPD patients independent of smoking status, serving as a potential therapeutic target for anti-EMT therapy in COPD.
Chronic obstructive pulmonary disease (COPD) is a progressive and partially reversible disease characterized by airflow obstruction resulting from a chronic inflammatory response of the airways and lungs to noxious particles or gases . The prevalence of COPD ranges from 8.80 to 14.53% among different countries , and COPD accounts for over 2 million deaths worldwide every year . It has been predicted that COPD will rise from the fourth to the third leading cause of global death from 2013 to 2020 . Therefore, it is important to understand the underlying mechanisms of COPD development and identify novel therapeutic targets for COPD treatment.
Small airways (< 2 mm internal diameter) are the major sites of airflow obstruction in patients with COPD . Small airway remodeling occurs in the early stage of COPD and largely contributes to airway obstruction in the disease . Peribronchiolar fibrosis, airway thickening resulting from the excess amount of extracellular matrix (ECM), loss of alveolar attachments, and inflammatory cell infiltration are considered markers of small airway remodeling in COPD . Recent studies have suggested epithelial–mesenchymal transition (EMT), the conversion of adherent epithelial cells into migratory mesenchymal cells, as a core pathological factor in airway remodeling during COPD development [8, 9]. The airway epithelium from COPD patient exhibits persistent activation of EMT, with decreased protein levels of epithelial markers such as E-cadherin and zona occludence-1 (ZO-1) and increased protein levels of mesenchymal markers, including N-cadherin, vimentin, and α-smooth muscle actin (α-SMA) [10,11,12]. Investigators have demonstrated that transforming growth factor-beta 1 (TGF-β1) expression is upregulated in the airway epithelium of COPD patients and that TGF-β1 induces EMT in the bronchial epithelial cells [13, 14]. However, the underlying mechanism underlying the regulation of TGF-β1 in EMT during COPD small airway remodeling is not fully understood.
Although smoking is the primary risk factor for COPD, genetic variation also contributes to the COPD genesis . Studies have reported that variants in the family with sequence similarity 13 member A (FAM13A) are associated with chronic lung diseases, including COPD, asthma, lung cancer, and pulmonary fibrosis [16,17,18,19]. We previously demonstrated that a genetic variant in FAM13A correlates with the increased risk of COPD in never-smokers . This suggests that FAM13A might be involved in the genesis and development of COPD independent of smoking status. FAM13A protein contains a rho GTPase activating protein domain that inactivates EMT-associated GTPases with converting GTP to GDP . FAM13A has been shown to regulate EMT in human bronchial epithelial cells and pulmonary artery endothelial cells [19, 22]. Based on these findings, we hypothesized that FAM13A might participate in TGF-β1-induced EMT during small airway remodeling in COPD.
To test our hypothesis, we measured the FAM13A expression in the small airway epithelium of lung tissue samples from COPD and non-COPD patients. We evaluated the correlations of FAM13A expression with COPD severity and EMT marker protein levels. We also investigated the roles of FAM13A in cell proliferation and motility and TGFβ1-induced EMT in human bronchial epithelial cell line BEAS-2B. Our findings identify FAM13A as a novel mediator of TGFβ1-induced EMT in the small airways, providing FAM13A as a potential therapeutic target for COPD treatment.
Patients and specimens
A total of 69 patients were recruited in the Department of Thoracic Surgery, General Hospital of Ningxia Medical University between January 2018 and December 2019. These patients were 35 patients with lung adenocarcinoma, 17 patients with lung squamous cell carcinoma, 12 patients with benign lung nodules, and 5 patients with bullous emphysema. Before the study, all the patients with lung adenocarcinoma, lung squamous cell carcinoma, and benign lung nodules have undergone lobectomy. The patients with bullous emphysema have undergone lung volume reduction surgery.
The inclusion criteria were: (1) lung cancer patients presenting a single nodule without lung metastasis in preoperative CT scan; (2) COPD patients with a post-bronchodilator forced expiratory volume in 1 s (FEV1)/forced volume vital capacity (FVC) ratio < 70%; (3) smokers (smoking history > 5 years and quitting smoking < 2 years) with a post-bronchodilator FEV1/FVC ratio > 70%; (4) never-smokers with a post-bronchodilator FEV1/FVC ratio > 70%; (5) pulmonary function testing was performed within 3 months before surgery.
The exclusion criteria were: (1) chronic lung diseases, such as asthma, tuberculosis, pneumonia, and interstitial lung disease; (2) acute COPD or discontinuing treatment < 2 months; (3) declining research participation. The patients were divided into non-smokers without COPD (n = 22), smokers without COPD (n = 23), non-smokers with COPD (n = 10), and smokers with COPD (n = 14) groups. The clinical characteristics of the patients were summarized in Table 1.
Lung tissue samples (5 cm × 5 cm; > 5 cm from the tumor margin) were obtained from the patients during lobectomy or lung volume reduction surgery. All tissue samples were immediately processed upon surgical removal. After fixing with 4% paraformaldehyde, the tissue samples were dehydrated in a graded ethanol series, paraffin-embedded, and cut into 4-µm-thick sections. This study was approved by the Ethics Committee of General Hospital of Ningxia Medical University. Written consent was obtained from all patients.
Cell culture and TGF-β1 stimulation
BEAS-2B cells were provided by Professor Fuqiang Wen at West China Hospital, Sichuan University (Chengdu, Sichuan, China). The cells were maintained in DMEM/F12 (1:1; Gibco, Thermo Fisher Scientific, Waltham, MA, USA) supplemented with 10% fetal bovine serum (FBS; Gibco) and 1% penicillin/streptomycin (Solarbio, Beijing, China) in a humidified atmosphere of 5% CO2 at 37 °C. TGF-β1 in various concentrations (PepproTech, Rocky Hill, NJ, USA) was used to treat BEAS-2B cells (see concentrations in later methods).
Construction of lentiviral expression vectors and transfection
For the overexpression of FAM13A, we collected RNA samples from clinical COPD patients and examined the mRNA expression of the five transcripts of FAM13A according to entries in NCBI and identified the mRNA level of the FAM13A transcript (NM_001015045.2) as higher than those of the other four transcripts. Therefore, the coding sequence of human FAM13A (NM_001015045.2) was subsequently cloned into a lentiviral vector GV492 (Genechem, Shanghai, China). Double-stranded DNA-coding small hairpin RNA (shRNA) for FAM13A (5′-ACACCAAACAGCAGAGAAA-3′, synthesized by Genechem) was cloned into a lentiviral vector GV248 (Genechem). Lentivirus was packaged in HEK293T cells. For lentiviral transduction, BEAS-2B cells were transduced with 15–multiplicity of infection of lentiviral vectors expressing FAM13A (LV-FAM13A) or shRNA against FAM13A (shFAM13A) for 48–72 h (Additional file 1: Fig. S1). BEAS-2B cells transfected with the empty GV492 vector and the empty GV248 vector were used as the negative control for overexpression and knockdown, respectively.
Measurement of small airway wall thickness
The deparaffined sections were stained with hematoxylin and eosin using standard methods for histological examination. The absence of tumor tissue was confirmed by two different pathologists who determined more accurately by observing the morphology of the cells in the HE-stained sections. The airway wall thickness was measured in 3–5 small airways with a basement membrane perimeter < 1000 μm in each section, as previously described . Image-pro Plus (IPP) 6.0 software (Media Cybernetics, Bethesda, MA, USA) conferred quantification of stained samples, which can quantify a length directly, quantify an area by outlining the staining and calculation, or quantify the staining by weighing the optical density and calculation. Using IPP 6.0, the length of the basement membrane (Pbm), the area defined by the basement membrane (Abm), and the area defined by the total adventitial perimeter (Ao) were all quantified. The calculation formula used was Thickness = total bronchial wall area (WAt)/Pbm, which was calculated as Ao − Abm.
Measurement of collagen area
The collagen fibers in the lung tissue samples from patients were stained using a Masson’s trichrome staining kit (Solarbio, Beijing, China) following the manufacturer’s instruction. The Blue-stained collagen fibers were observed under an Olympus microscope (Olympus, Japan). The percentage of the collagen area was measured as previously described . Briefly, the total collagen area (stained by Masson's trichrome) was quantified using the IPP6.0 software. The calculation formula for the percentage of the collagen area was: Collagen area (%) = Area/Pbm.
Immunohistochemical (IHC) staining
The protein levels of small airway remodeling markers and EMT markers were determined by immunohistochemical staining. Briefly, the paraffin-embedded sections were dewaxed in xylene and dehydrated in ethanol, followed by incubation with 3% H2O2 for 10 min. After antigen retrieval, the sections were blocked with 10% normal goat serum, then incubated with anti-FAM13A (1:800; HPA038109, Sigma, Dorset, UK), anti-fibronectin (1:500; ab2413, Abcam, Cambridge, UK), anti-E-cadherin (1:400; ab76055, Abcam), anti-α-SMA (1:500; ab7817, Abcam), anti-ZO-1 (1:200; ab216880, Abcam), anti-MMP9 (1:500; ab76003, Abcam), anti-N-cadherin (1:500; ab18203, Abcam), anti-collagen I (1:500; ab34710, Abcam), or anti-TGF-β1 (1:1000; ab27969, Abcam) antibody. Detection of the antigen–antibody complex was performed using a secondary antibody and a DAB detection kit. All the secondary antibodies and the DAB detection kit were purchased from Beijing Zhong Shan-Golden Bridge Biological Technology (Beijing, China). The results were examined and photographed under an Olympus microscope (Olympus, Japan) at a magnification of × 400. Except for vimentin, the expression of target proteins was expressed as the average optical density (AOD) of the positive staining in 5 randomly selected fields as previously described . Briefly, the integrated optical density (IOD) and the total area of staining (Area) were quantified by IPP6.0. The calculation formula of AOD was: AOD = IOD/Area. Vimentin expression was calculated as vimentin-positive cell numbers/Pbm, as previously described . The former was manually counted while Pbm was quantified using IPP 6.0 as already described.
Immunocytochemical and immunofluorescent staining
FMA13A expression in response to TGF-β1 stimulation was detected by immunocytochemical and immunofluorescent staining. Briefly, BEAS-2B cells were seeded on a coverslip in a 24-well plate at a density of 1 × 104 cells/well and incubate overnight. Cells were treated with TGF-β1 (10 ng/mL; PepproTech) for 48 h. For immunocytochemical staining, cells were fixed with 4% paraformaldehyde for 20 min at room temperature. After blocking with 10% normal goat serum for 20 min, cells were incubated with anti-FAM13A (1:100, Sigma-Aldrich, St. Louis, MO, USA) at 4 °C overnight, followed by incubation with HRP-conjugated secondary antibody (Beijing Zhong Shan-Golden Bridge Biological Technology) for 20 min at room temperature and for 10 min counter-staining by using haematin and hydrochloric acid 1% in 70% alcohol. After rinsing with tap water, the slides were dried and sealed with neutral resin. Images were acquired under an Olympus microscope.
For immunofluorescent staining, cells were fixed with 4% paraformaldehyde for 20 min at room temperature. For FAM13A detection, after blocking with 5% bovine serum albumin in phosphate-buffered saline (PBS), cells were incubated with anti-FAM13A at 4 °C overnight with a dilution of 1:100, followed by 3 washes with PBS. The cells were then incubated with a FITC-conjugated secondary antibody (Goat Anti-Rabbit IgG H&L, ab150081, Abcam) in the dark for 1 h at room temperature. For filamentous actin (F-actin) detection, after blocking with 1% bovine serum albumin, cells were incubated with phalloidin-iFluor 647 Conjugate (40734ES75, Yeasen biotech, Shanghai, China) in the dark for 30 min at room temperature. The fluorescence was visualized using an Olympus fluorescence microscope.
Bromodeoxyuridine (BrdU) staining
A total of 10 mg BrdU powder (ab142567, Beijing Zhong Shan-Golden Bridge Biological Technology) was dissolved in 10 mL pure water (injection-grade) to make a BrdU solution at a stock concentration of 1.0 mg/mL were stored at 4 ℃ in the dark. BEAS-2B cells were prepared in a 12-well plate and cultured, as previously described, until the cell density of 50%. The BrdU stock solution was further diluted to a final concentration of 0.03 μg/mL and filled into the wells were incubated at 37 ℃ for 4 h before BrdU was removed, followed by washing with PBS for three times. The cells were fixed with 4% paraformaldehyde for 30 min and washed using PBS for three times (5 min each time). Next, HCL (2 mol/L) was applied to the cells for 30 min at room temperature to dissolve DNA. After washing with PBS for three times (5 min each time), the cells were incubated with the blocked membranes in 5% BSA for 1 h at room temperature. The cells were subsequently incubated overnight with the primary antibody (ab6326, 1:200, Abcam) at 4 ℃ and washed with PBS for three times. This was followed by incubation with the secondary antibody (1:200, A-21209, Invitrogen, USA) for 2 h and washing with PBS for three times. The 4′, 6-diamidino-2-phenylindole (DAPI, 1:500, Sigma, USA) was then added to stand for 2 min before PBS washing for three times. Lastly, a fluorescent quenching agent (Solarbio, China) was added for photo-taken.
Cell migration and invasion assays
Cell migration and invasion were measured with the wound healing assay and Transwell chamber assay, respectively. For the wound healing assay, BEAS-2B cells transfected with LV-FAM13A, shFAM13A, or corresponding negative control were seeded into a 6-well plate at a density of 1 × 105 cells per well and cultured overnight. A 200 µL micropipette tip was used to generate a 2 mm-wide scratch line in the cell monolayer. Cells were allowed to migrate for 24 h. Images were captured under an inverted light microscope (Olympus, Japan).
For the invasion assay, BEAS-2B cells transfected with LV-FAM13A, shFAM13A, or corresponding negative control were loaded into a Matrigel-coated upper chamber (Corning, Corning, NY, USA) filled with serum-free DMEM. The lower chamber was filled with 500 μL DMEM containing 20% FBS. After incubation at 37 °C for 48 h, nonmigrating or non-invading cells remaining in the upper chamber were removed with a cotton swab. The migrating or invading cells adhering to the lower surface were fixed and stained with crystal violet (0.1%). The stained cells were counted in 3–5 randomly selected fields under an inverted light microscope (Olympus) at × 200 magnification. Images were acquired using an Olympus camera (Olympus, Japan).
Cell proliferation and apoptosis analyses
Cell proliferation and apoptosis were measured with flow cytometry analysis. For cell proliferation analysis, BEAS-2B cells (1–5 × 105/mL) were harvested at 48 h after transfection and fixed with 80% methanol for 5 min. Cells were washed with PBS and fixed with 4% paraformaldehyde for 20 min at room temperature, followed by penetration with 0.1% Triton X-100. After 30 min of blocking with 10% normal goat serum (Solarbio), cells were incubated with anti-Ki67 (1:100; ab16667; Abcam) for 30 min at 22 ℃. Cells were then incubated with secondary antibody (1:500, Goat Anti-Rabbit IgG H&L ab150083; Abcam) for 30 min at 22 ℃, followed by resuspension in PBS and flow cytometry analysis on a BD Accuri flow cytometer (Becton Dickinson, San Jose, CA, USA).
Cell apoptosis was measured using an Annexin V-PE/7-AAD apoptosis kit (Liankebio, Hangzhou, Zhejiang, China) following the manufacturer’s instruction. Briefly, BEAS-2B cells (1–3 × 106/mL) were harvested at 48 h after transfection. Cells were digested with enzymes free of ethylenediaminetetraacetic acid (EDTA) to prepare into single-cell suspensions. Following counting, the cells were re-suspended with Annexin binding buffer (containing HEPES buffer, calcium chloride, and Ca-free), supplemented with Annexin V reagents and PE in the kit. After washing with ice-cold PBS, cells were reconstituted in the binding buffer and stained with Annexin V-PE/7-AAD, followed by flow cytometry analysis on a BD Accuri flow cytometer.
Quantitative real-time PCR (qRT-PCR)
FAM13A mRNA expression in BEAS-2B cells after transfection was determined using qRT-PCR. Cells were seeded in a 6-well plate at a density of 2 × 105 cells/well and grown overnight. Cells were transfected with LV-FAM13A, shFAM13A, or corresponding negative control for 48 h. Total RNA was isolated using Trizol reagent (Abcom, USA) following the manufacturer’s instruction. cDNA was synthesized using a reverse transcription kit (Takara, Japan). qPCR was performed using SYBR Green (Takara) and gene-specific primers (FAM13A, forward: 5′-AAGTCTGCCTCATTGGCTGTGG-3′, reverse: 5′-CAGCCCATAACAGCAGATAGGC-3′; GAPDH, F: 5′-CACCATTGGCAATGAGCGGTTC-3′, reverse: 5′-AGGTCTTTGCGGATGTCCACGT-3′) on a quantitative PCR instrument (Roche, Germany). Gene expression was quantified using the 2△△Ct method.
Western blot analysis
BEAS-2B cells or human lung tissue samples were lysed using the protein isolation kit (#KGP250, KeyGEN BioTECH, Nanjing, China) or homogenized with lysis buffer containing phenylmethylsulphonyl fluoride on ice, followed by centrifugation at 12,000 rpm for 5 min. Protein concentrations were measured using a BCA quantification kit (Keygentec, Nanjing, Jiangsu, China). Proteins were separated on 10% SDS gel and transferred to a polyvinylidene fluoride membrane, followed by 1 h of blocking with 5% skim milk. The membrane was then incubated with anti-N-cadherin (1:100; ab18203; Abcam), anti-vimentin (1:5000; ab92547; Abcam), anti-α-SMA (1:500; ab7817; Abcam), anti-E-cadherin (1:1000; ab76055; Abcam), or anti-GAPDH (1:10,000; PAB36264; Bio swamp, Wuhan, Hubei, China) antibody overnight at 4 °C, followed by 3 washes with Tris-buffered saline containing 0.1% Tween 20 (TBST). The membrane was then incubated with horseradish peroxidase-conjugated secondary antibody (1:10,000; SAB43710 or SAB43711; Bio swamp) for 2 h at room temperature. After an additional three washes with TBST and the addition of exposure solution (Millipore, #WBKLS0500, USA), the protein bands were visualized using an enhanced chemiluminescence reagent (Life iLab, Shanghai, China) and analyzed using the Image J software (NIH, Bethesda, MD, USA).
Statistical analyses were carried out using Prism 6.0 software (Graphpad, San Diego, CA, USA). Data are presented as mean ± standard deviation. Differences between categorical variables were compared using the χ2 test. Comparison among groups were conducted using analysis of variance (ANOVA), followed by the Student–Newman–Keuls test. The association between FAM13A expression and variables of interest was analyzed by multivariate Cox proportional hazards analysis. The correlation analysis between FAM13A expression level and FEV1/FVC%, FEV1%pre, PO2, PCO2, E-cadherin, and vimentin were analyzed by partial correlation analysis the patient's age was the control variable. A value of P < 0.05 was considered statistically significant.
EMT and EMT inducer TGF-β1 are associated with small airway remolding in COPD patients independent of smoking status
To evaluate the correlation of small airway remodeling with COPD, we measured the remodeling markers in lung tissue surrounding the small airways of the patients. As shown in Additional file 1: Fig. S1A–C, the lung tissue samples of COPD patients exhibited remarkably increased small airway wall thickness, collagen fiber deposition, and the protein levels of collagen I, fibronectin, and MMP9, as compared with those of non-COPD patients. No significant differences were observed in these parameters between non-smokers and smokers with COPD. In non-COPD patients, we only observed significantly increased small airway wall thickness in the smokers than those in the non-smokers. These results suggest that COPD patients undergo small airway remodeling independent of smoking status.
Considering the crucial role of EMT in airway remodeling, we detected the EMT markers in the lung tissue samples from the patients. As shown in Additional file 1: Fig. S2A and S2B, compared with those of non-COPD patients, the lung tissue samples of COPD patients had significantly attenuated E-cadherin and ZO-1 protein levels and markedly enhanced α-SMA, vimentin, and N-cadherin protein levels. No significant differences were observed between the non-smokers and smokers with or without COPD. Similar results were observed in Western blot analysis, except for significant differences in α-SMA and N-cadherin levels among non-smokers and smokers without COPD (Additional file 1: Fig. S2C). These results suggest that EMT is activated in the small airways of COPD patients independent of smoking status.
Because TGF-β1 is a potent inducer of EMT and plays an important role in the pathogenesis of airway remodeling , we further detected TGF-β1 protein levels in the lung tissue samples from the patients. Additional file 1: Figure S2D shows that in non-smokers, the lung tissue samples of COPD patients had significantly upregulated TGF-β1 protein levels.
FAM13A is upregulated in the small airway epithelium in COPD patients and correlates with COPD severity and EMT marker expression
To explore the involvement of FAM13A in small airway remodeling in COPD, we measured FAM13A protein levels in the small airway epithelium of the patients. IHC staining showed that FAM13A protein levels were considerably upregulated in the small airway epithelium of COPD patients independent of smoking status (Fig. 1A).
Then, we assessed the correlations of FAM13A protein levels with COPD severity. As shown in Fig. 1B and C, FAM13A protein levels in the small airway epithelium were significantly and negatively correlated with FEV1/FVC ratio, pre-bronchodilator FEV1, and PO2 of the patients. An opposite trend was observed in the correlation of FAM13A protein levels with PCO2, despite a statistical non-significance. These data suggest that FAM13A protein levels in the small airway epithelium are positively correlated with COPD severity. Furthermore, FAM13A protein expression was significantly and negatively correlated with E-cadherin protein expression, but positively correlated with vimentin expression (Fig. 1D), suggesting that FAM13A is associated with EMT in small airways.
TGF-β1 stimulation upregulates FAM13A expression in BESA-2B cells
To understand the mechanism underlying the involvement of FAM13A in lung EMT, we determined FAM13A protein expression in BEAS-2B cells in response to TGF-β1 stimulation. As a result, TGF-β1 stimulation significantly enhanced FAM13A protein expression in a dose-dependent manner (Fig. 2A) in BEAS-2B cells. Immunocytochemical (Fig. 2B) and immunofluorescent (Fig. 2C) staining demonstrated consistent results. This finding suggests that FAM13A might be a downstream effector of TGF-β1 in the regulation of EMT.
FAM13A promotes cell proliferation, migration, and invasion and is required for suppressing apoptosis in BEAS-2B cells
Because epithelial cells acquire migratory and invasive properties during EMT , we sought to explore the functions of FAM13A in the proliferation and motility of BEAS-2B cells. Figure 2D and E show the transfection efficiency of lentiviral vectors overexpressing FAM13A or shFAM13A. Ki67 staining revealed that FAM13A overexpression significantly promoted cell proliferation, whereas FAM13A silencing significantly inhibited cell proliferation in BEAS-2B cells (Fig. 3A). Consistent results were observed by BrdU staining (Fig. 3B). In terms of cell apoptosis, although FAM13A overexpression failed to inhibit BEAS-2B cell apoptosis, FAM13A knockdown noticeably promoted BEAS-2B cell apoptosis (Fig. 3C). This finding suggests that FAM13A is required for the suppression of BEAS-2B cell apoptosis. In addition, FAM13A overexpression promoted cell migration and invasion, whereas FAM13A silencing inhibited cell motility in BEAS-2B cells (Fig. 4A and B). Taken together, these data suggest that FAM13A promotes the phenotype switch from epithelial to mesenchymal cells in airway epithelial cells.
Knockdown of FAM13A abrogates TGF-β1-induced alterations in EMT marker protein expression in BEAS-2B cells
Next, we explored whether FAM13A is required for TGF-β1-induced EMT in airway epithelial cells. Western blotting revealed that FAM13A overexpression significantly attenuated E-cadherin and ZO-1 protein levels while enhancing vimentin, α-SMA, and N-cadherin protein levels in BEAS-2B cells (Fig. 4C). Although FAM13A knockdown did not affect EMT marker protein levels in BEAS-2B cells in the absence of TGF-β1, FAM13A knockdown partially reversed TGF-β1-induced alterations in the EMT marker protein levels in BEAS-2B cells (Fig. 4D). This finding suggests that TGF-β1 induces EMT in BEAS-2B cells at least partially via FAM13A.
Past clinical research has suggested that smoking is the most crucial risk factor for COPD. However, genetic factors also play an essential role in the susceptibility to COPD. In this study, we demonstrated that regardless of smoking status, the small airways of COPD patients exhibited active remodeling and EMT compared with non-COPD patients. EMT has been reported as an active process in the airways of COPD patients . TGF-β1 induces EMT in the bronchial epithelial cells [29,30,31]. However, it remains unclear how TGF-β1 promotes EMT in small airway remodeling in COPD. This study shows that FAM13A can promote the proliferation, migration, and invasion of BEAS-2B cells, and knockout of FAM13A can partially reverse the changes in the expression of EMT markers induced by TGF-β1. These are consistent with previous reports concluding that FAM13A may be a potential EMT-correlated airway obstruction susceptibility gene [19, 32]. In their studies, FAM13A was demonstrated to increase the susceptibility to airway obstruction independent of smoking status  and to closely relate to the EMT of cystic fibrosis airway . In summary, this study implies that FAM13A mediates TGF-β1-induced EMT in small airways in COPD, providing FAM13A as a potential therapeutic target for anti-EMT treatment of COPD.
Small airway remodeling in COPD involves fibroblast proliferation and migration, excessive collagen deposition, and an imbalance between ECM synthesis and degradation, resulting in airway wall thickening and luminal narrowing . Increased deposition of ECM components, such as collagen fibers, fibronectin, α-SMA, and MMP-9, have been extensively used as airway remodeling markers by many clinical and experimental studies [33,34,35]. Our results showed that these remodeling markers were remarkably increased in the small airway epitheliums of COPD patients than non-COPD patients. Since EMT has been suggested as a core pathological factor in airway remodeling during COPD development [8, 9], we measured EMT marker expression in the small airway epithelium. As expected, the increases in the remodeling markers were accompanied by the upregulation of mesenchymal marker expression and downregulation of epithelial marker expression. These findings collectively indicate active remodeling and EMT processes in the small airways of COPD patients.
Interestingly, although studies have reported that smoking activates EMT and airway remodeling in COPD and other lung diseases [10, 36], our results showed that the alterations in these markers were independent of smoking status in COPD patients. This suggests that other factors contribute to EMT and airway remodeling, at least in COPD. Some investigators have demonstrated that in addition to environmental determinants, genetic risk factors, including genetic variants, epigenetic changes, and somatic mutations, may cause COPD [37,38,39]. Based on previous studies, including ours, demonstrating the associations of FAM13A with the increased risk of COPD in never-smokers and with lung EMT [19, 20, 22, 40], we sought to investigate the expression pattern and the roles of FAM13A in COPD. We found that the FAM13A protein level was elevated in the small airway epithelium in COPD patients compared with non-COPD patients, independent of smoking status, consistent with a previous report that FAM13A protein levels are elevated in very severe human COPD lungs . In addition, we noticed that FAM13A protein expression negatively correlated with FEV1% and PO2 in COPD patients. This was consistent with a previous study that a FAM13A variant is associated with lower FEV1/FVC levels in never-smokers with COPD . This finding suggests that high FAM13A expression is associated with severe COPD.
During COPD development, repetitive injury resulting from chronic exposure to environmental or genetic factors may cause persistent activation of repair processes in airway epithelium, such as EMT and altered proliferation and migration of bronchial epithelial cells . In this study, FAM13A overexpression significantly promoted cell proliferation, migration, and invasion, whereas FAM13A silencing significantly inhibited these cell behaviors in BEAS-2B cells. Interestingly, it appears the regulation of FAM13A of cell proliferation is complicated and associated with environmental factors. Knockout of FAM13A in mice without exposure to smoke resulted in accelerated lung cell proliferation compared with FAM13A overexpression. Under the smoke condition, overexpression of FAM13A augmented the proliferation of alveolar septal cells . In this study, we also attributed the proliferation of BEAS-2B cells to FAM13A and altering the expression of FAM13A induced the TGF-β1/ Smad signaling. Although in several in vitro studies, TGF-β1 has been demonstrated to be inhibitory to the growth of epithelial and tumor cells. Other reports support that TGF-β1 can confer bidirectional regulation of cell proliferation under different circumstances [42,43,44]. For instance, overexpression of TGF-β1 can cause the proliferation of epithelia in the head and neck, thereby promoting squamous cell carcinoma . From our viewpoint, the involvement of TGF-β1/Smad in multiple signaling pathways cannot be ruled out, and how FAM13A functions and interact with TGF-β1/Smad signaling remains to be further elucidated. This will inspire us to perform future experiments with enlarged sample numbers, multivariable analyses, and in vivo approaches to increase knowledge of the molecular pathways.
Furthermore, FAM13A upregulation was implicated in small airway EMT, suggested by its significant correlations with attenuated E-cadherin expression and enhanced vimentin expression in the small airway epithelium. Hence, we speculated that FAM13A might facilitate COPD progression via EMT in small airways. These indicate that FAM13A is involved in the phenotype switch from epithelial to mesenchymal cells in airway epithelial cells. Despite its positive function in embryonic development and tissue repair, EMT has been extensively shown to play a key role in tumor formation and metastasis . Recent studies also indicated that EMT might be vital to the pathogenesis of airway remodeling in COPD and other fibrotic diseases . It can increase the number of fibroblasts and myofibroblasts in the airway wall and thus mediate airway remodeling. More than 50% of patients with COPD are accompanied by lung cancer [47, 48].
Furthermore, the expression of FAM13A is significantly increased in both COPD and lung cancer tissues [40, 49]. Some common mechanisms are likely shared by these two diseases [8, 12, 50, 51]. FAM13A has been reported as a downstream effector of EMT inducer TGF-β in A549 lung cancer cells . Meanwhile, FAM13A is a susceptibility gene for lung cancer . Its rs9224 site can change the susceptibility of squamous cell carcinoma by affecting the binding of miRNA-22-5p . Our in vitro results showed that TGF-β1 stimulation upregulated FAM13A expression in BEAS-2B cells in a dose-dependent manner. In the meantime, TGF-β1 expression was significantly upregulated in lung tissue surrounding the small airways of COPD patients compared with non-COPD patients, independent of smoking status. These findings suggest a presence of TGF-β1/FAM13A signaling in the regulation of small airway EMT in COPD. Indeed, loss of FAM13A resulted in a partial but significant reversal of TGF-β1-induced alterations in EMT marker expression in BEAS-2B cells, suggesting that TGF-β1 regulates EMT in COPD airway remodeling at least partially via FAM13A. Previous studies have shown that FAM13A regulates RhoA activity through the rhoGAP domain, which causes changes in EMT markers of bronchial epithelial cells in lung cystic fibrosis . In addition, in the pulmonary hypertension model, FAM13A can alleviate EMT of endothelial cells by inhibiting the β-catenin signaling of pulmonary artery endothelial cells . As such, it will be of interest in a future study to find if the participation of FAM13A in the occurrence of EMT is a common mechanism linked with both COPD and lung cancer.
The limitation of this study is that even though two different pathologists determined the absence of lung cancer by observing the morphology of each section. However, due to unavailability, no specific marker was used to identify tumor cells.
In conclusion, in this study, we demonstrated that FAM13A expression was upregulated in the small airway epithelium of COPD patients independent of smoking status. FAM13A promoted cell proliferation and motility and mediated TGF-β1-induced EMT marker alterations in BEAS-2B cells. Our results identify FAM13A as a novel mediator of TGF-β1 signaling in small airway remodeling, providing FAM13A as a potential therapeutic target for anti-EMT therapy against COPD.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Family with sequence similarity 13 member A
Chronic obstructive pulmonary disease
α-Smooth muscle actin
Transforming growth factor-beta 1
Small hairpin RNA
Lentiviral vectors expressing FAM13A
Vestbo J, Hurd SS, Agustí AG, Jones PW, Vogelmeier C, Anzueto A, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am J Respir Crit Care Med. 2013;187:347–65.
Varmaghani M, Dehghani M, Heidari E, Sharifi F, Moghaddam SS, Farzadfar F. Global prevalence of chronic obstructive pulmonary disease: systematic review and meta-analysis. East Mediterr Health J. 2019;25:47–57.
Global, regional, and national age-sex specific mortality for 264 causes of death, 1980–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet. 2017;390:1151–210.
Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380:2095–128.
Wang Q, Wang Y, Zhang Y, Zhang Y, Xiao W. The role of uPAR in epithelial-mesenchymal transition in small airway epithelium of patients with chronic obstructive pulmonary disease. Respir Res. 2013;14:67.
Hogg JC, Paré PD, Hackett TL. The contribution of small airway obstruction to the pathogenesis of chronic obstructive pulmonary disease. Physiol Rev. 2017;97:529–52.
Sköld CM. Remodeling in asthma and COPD–differences and similarities. Clin Respir J. 2010;4(Suppl 1):20–7.
Mahmood MQ, Sohal SS, Shukla SD, Ward C, Hardikar A, Noor WD, et al. Epithelial mesenchymal transition in smokers: large versus small airways and relation to airflow obstruction. Int J Chron Obstruct Pulmon Dis. 2015;10:1515–24.
Sohal SS. chronic obstructive pulmonary disease (COPD) and lung cancer: epithelial mesenchymal transition (EMT), the missing link? EBioMedicine. 2015;2:1578–9.
Milara J, Peiró T, Serrano A, Cortijo J. Epithelial to mesenchymal transition is increased in patients with COPD and induced by cigarette smoke. Thorax. 2013;68:410–20.
Sohal SS, Reid D, Soltani A, Ward C, Weston S, Muller HK, et al. Reticular basement membrane fragmentation and potential epithelial mesenchymal transition is exaggerated in the airways of smokers with chronic obstructive pulmonary disease. Respirology. 2010;15:930–8.
Sohal SS, Reid D, Soltani A, Ward C, Weston S, Muller HK, et al. Evaluation of epithelial mesenchymal transition in patients with chronic obstructive pulmonary disease. Respir Res. 2011;12:130.
Doerner AM, Zuraw BL. TGF-beta1 induced epithelial to mesenchymal transition (EMT) in human bronchial epithelial cells is enhanced by IL-1beta but not abrogated by corticosteroids. Respir Res. 2009;10:100.
Papaporfyriou A, Loukides S, Kostikas K, Simoes DCM, Papatheodorou G, Konstantellou E, et al. Increased levels of osteopontin in sputum supernatant in patients with COPD. Chest. 2014;146:951–8.
Schwartz AG. Genetic epidemiology of cigarette smoke-induced lung disease. Proc Am Thorac Soc. 2012;9:22–6.
Cho MH, Boutaoui N, Klanderman BJ, Sylvia JS, Ziniti JP, Hersh CP, et al. Variants in FAM13A are associated with chronic obstructive pulmonary disease. Nat Genet. 2010;42:200–2.
Li X, Howard TD, Moore WC, Ampleford EJ, Li H, Busse WW, et al. Importance of hedgehog interacting protein and other lung function genes in asthma. J Allergy Clin Immunol. 2011;127:1457–65.
Young RP, Hopkins RJ, Hay BA, Whittington CF, Epton MJ, Gamble GD. FAM13A locus in COPD is independently associated with lung cancer - evidence of a molecular genetic link between COPD and lung cancer. Appl Clin Genet. 2011;4:1–10.
Corvol H, Rousselet N, Thompson KE, Berdah L, Cottin G, Foussigniere T, et al. FAM13A is a modifier gene of cystic fibrosis lung phenotype regulating rhoa activity, actin cytoskeleton dynamics and epithelial-mesenchymal transition. J Cyst Fibros. 2018;17:190–203.
Zhang Y, Qiu J, Zhang P, Zhang J, Jiang M, Ma Z. Genetic variants in FAM13A and IREB2 are associated with the susceptibility to COPD in a Chinese rural population: a case-control study. Int J Chron Obstruct Pulmon Dis. 2018;13:1735–45.
Amin E, Jaiswal M, Derewenda U, Reis K, Nouri K, Koessmeier KT, et al. Deciphering the molecular and functional basis of RHOGAP family proteins: a systematic approach toward selective inactivation of rho family proteins. J Biol Chem. 2016;291:20353–71.
Rinastiti P, Ikeda K, Rahardini EP, Miyagawa K, Tamada N, Kuribayashi Y, et al. Loss of family with sequence similarity 13, member A exacerbates pulmonary hypertension through accelerating endothelial-to-mesenchymal transition. PLoS ONE. 2020;15:e0226049.
Bracke KR, D’Hulst AI, Maes T, Demedts IK, Moerloose KB, Kuziel WA, et al. Cigarette smoke-induced pulmonary inflammation, but not airway remodelling, is attenuated in chemokine receptor 5-deficient mice. Clin Exp Allergy. 2007;37:1467–79.
Cho JY, Miller M, Baek KJ, Han JW, Nayar J, Rodriguez M, et al. Immunostimulatory DNA inhibits transforming growth factor-beta expression and airway remodeling. Am J Respir Cell Mol Biol. 2004;30:651–61.
Guo X, Fan Y, Cui J, Hao B, Zhu L, Sun X, et al. NOX4 expression and distal arteriolar remodeling correlate with pulmonary hypertension in COPD. BMC Pulm Med. 2018;18:111.
Kamitani S, Yamauchi Y, Kawasaki S, Takami K, Takizawa H, Nagase T, et al. Simultaneous stimulation with TGF-β1 and TNF-α induces epithelial mesenchymal transition in bronchial epithelial cells. Int Arch Allergy Immunol. 2011;155:119–28.
Lamouille S, Xu J, Derynck R. Molecular mechanisms of epithelial–mesenchymal transition. Nat Rev Mol Cell Biol. 2014;15:178–96.
Gohy ST, Hupin C, Fregimilicka C, Detry BR, Bouzin C, Gaide Chevronay H, et al. Imprinting of the COPD airway epithelium for dedifferentiation and mesenchymal transition. Eur Respir J. 2015;45:1258–72.
Ladak SS, Roebuck E, Powell J, Fisher AJ, Ward C, Ali S. The role of miR-200b-3p in modulating TGF-β1-induced injury in human bronchial epithelial cells. Transplantation. 2019;103:2275–86.
Yang ZC, Yi MJ, Ran N, Wang C, Fu P, Feng XY, et al. Transforming growth factor-β1 induces bronchial epithelial cells to mesenchymal transition by activating the Snail pathway and promotes airway remodeling in asthma. Mol Med Rep. 2013;8:1663–8.
Pattarayan D, Sivanantham A, Krishnaswami V, Loganathan L, Palanichamy R, Natesan S, et al. Tannic acid attenuates TGF-β1-induced epithelial-to-mesenchymal transition by effectively intervening TGF-β signaling in lung epithelial cells. J Cell Physiol. 2018;233:2513–25.
van der Plaat DA, de Jong K, Lahousse L, Faiz A, Vonk JM, van Diemen CC, et al. Genome-wide association study on the FEV(1)/FVC ratio in never-smokers identifies HHIP and FAM13A. J Allergy Clin Immunol. 2017;139:533–40.
Grzela K, Litwiniuk M, Zagorska W, Grzela T. Airway remodeling in chronic obstructive pulmonary disease and asthma: the role of matrix metalloproteinase-9. Arch Immunol Ther Exp (Warsz). 2016;64:47–55.
Sun C, Zhu M, Yang Z, Pan X, Zhang Y, Wang Q, et al. LL-37 secreted by epithelium promotes fibroblast collagen production: a potential mechanism of small airway remodeling in chronic obstructive pulmonary disease. Lab Invest. 2014;94:991–1002.
Annoni R, Lanças T, Yukimatsu Tanigawa R, de Medeiros MM, de Morais FS, Bruno A, et al. Extracellular matrix composition in COPD. Eur Respir J. 2012;40:1362–73.
Hou W, Hu S, Li C, Ma H, Wang Q, Meng G, et al. Cigarette smoke induced lung barrier dysfunction, EMT, and tissue remodeling: a possible link between COPD and lung cancer. Biomed Res Int. 2019;2019:2025636.
Decramer M, Janssens W, Miravitlles M. Chronic obstructive pulmonary disease. Lancet. 2012;379:1341–51.
Yao H, Rahman I. Role of histone deacetylase 2 in epigenetics and cellular senescence: implications in lung inflammaging and COPD. Am J Physiol Lung Cell Mol Physiol. 2012;303:L557–66.
Tzortzaki EG, Dimakou K, Neofytou E, Tsikritsaki K, Samara K, Avgousti M, et al. Oxidative DNA damage and somatic mutations: a link to the molecular pathogenesis of chronic inflammatory airway diseases. Chest. 2012;141:1243–50.
Jiang Z, Lao T, Qiu W, Polverino F, Gupta K, Guo F, et al. A chronic obstructive pulmonary disease susceptibility gene, FAM13A, regulates protein stability of β-catenin. Am J Respir Crit Care Med. 2016;194:185–97.
Ganesan S, Sajjan US. Repair and remodeling of airway epithelium after injury in chronic obstructive pulmonary disease. Curr Respir Care Rep. 2013;2.
Grande JP. Role of transforming growth factor-beta in tissue injury and repair. Proc Soc Exp Biol Med. 1997;214:27–40.
Lu SL, Reh D, Li AG, Woods J, Corless CL, Kulesz-Martin M, et al. Overexpression of transforming growth factor beta1 in head and neck epithelia results in inflammation, angiogenesis, and epithelial hyperproliferation. Cancer Res. 2004;64:4405–10.
Roberts AB, Anzano MA, Wakefield LM, Roche NS, Stern DF, Sporn MB. Type beta transforming growth factor: a bifunctional regulator of cellular growth. Proc Natl Acad Sci USA. 1985;82:119–23.
Kalluri R, Weinberg RA. The basics of epithelial-mesenchymal transition. J Clin Invest. 2009;119:1420–8.
Jolly MK, Ward C, Eapen MS, Myers S, Hallgren O, Levine H, et al. Epithelial–mesenchymal transition, a spectrum of states: role in lung development, homeostasis, and disease. Dev Dyn. 2018;247:346–58.
Young RP, Hopkins RJ, Christmas T, Black PN, Metcalf P, Gamble GD. COPD prevalence is increased in lung cancer, independent of age, sex and smoking history. Eur Respir J. 2009;34:380–6.
Spyratos D, Papadaki E, Lampaki S, Kontakiotis T. Chronic obstructive pulmonary disease in patients with lung cancer: prevalence, impact and management challenges. Lung Cancer (Auckl). 2017;8:101–7.
Eisenhut F, Heim L, Trump S, Mittler S, Sopel N, Andreev K, et al. FAM13A is associated with non-small cell lung cancer (NSCLC) progression and controls tumor cell proliferation and survival. Oncoimmunology. 2017;6:1256526.
Sohal SS, Mahmood MQ, Walters EH. Clinical significance of epithelial mesenchymal transition (EMT) in chronic obstructive pulmonary disease (COPD): potential target for prevention of airway fibrosis and lung cancer. Clin Transl Med. 2014;3:33.
Thiery JP, Sleeman JP. Complex networks orchestrate epithelial–mesenchymal transitions. Nat Rev Mol Cell Biol. 2006;7:131–42.
Ziółkowska-Suchanek I, Mosor M, Gabryel P, Grabicki M, Żurawek M, Fichna M, et al. Susceptibility loci in lung cancer and COPD: association of IREB2 and FAM13A with pulmonary diseases. Sci Rep. 2015;5:13502.
Yu Y, Mao L, Lu X, Yuan W, Chen Y, Jiang L, et al. Functional Variant in 3’UTR of FAM13A Is Potentially Associated with Susceptibility and Survival of Lung Squamous Carcinoma. DNA Cell Biol. 2019;38:1269–77.
This work was supported by the National Natural Science Foundation of China (Grant No. 81760011), Ningxia Natural Science Foundation of China (Grant No. 2020AAC03378), the National Key Research and Development Program (Grant No. 2018YFC1313704). The funders had no role in study design, data collection and analysis, decision to publish, or manuscript preparation.
Ethics approval and consent to participate
This study was approved by the Ethics Committee of General Hospital of Ningxia Medical University. Written consent was obtained from all patients.
Consent for publication
Written consent was obtained from all patients.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhu, J., Wang, F., Feng, X. et al. Family with sequence similarity 13 member A mediates TGF-β1-induced EMT in small airway epithelium of patients with chronic obstructive pulmonary disease. Respir Res 22, 192 (2021). https://doi.org/10.1186/s12931-021-01783-z
- Epithelial–mesenchymal transition
- Chronic obstructive pulmonary disease
- Small airway epithelium