Abstract
Aging is one of the key contributing factors for chronic obstructive pulmonary diseases (COPD) and other chronic inflammatory lung diseases. Cigarette smoke is a major etiological risk factor that has been shown to alter cellular processes involving mitochondrial function, cellular senescence and telomeric length. Here we determined how aging contribute to the alteration in the gene expression of above mentioned cellular processes that play an important role in the progression of COPD and IPF. We hypothesized that aging may differentially alter the expression of mitochondrial, cellular senescence and telomere genes in smokers and patients with COPD and IPF compared to non-smokers. Total RNA from human lung tissues from non-smokers, smokers, and patients with COPD and IPF were processed and analyzed based on their ages (younger: <55 yrs and older: >55 yrs). NanoString nCounter panel was used to analyze the gene expression profiles using a custom designed codeset containing 112 genes including 6 housekeeping controls (mitochondrial biogenesis and function, cellular senescence, telomere replication and maintenance). mRNA counts were normalized, log2 transformed for differential expression analysis using linear models in the limma package (R/Bioconductor). Data from non-smokers, smokers and patients with COPD and IPF were analyzed based on the age groups (pairwise comparisons between younger vs. older groups). Several genes were differentially expressed in younger and older smokers, and patients with COPD and IPF compared to non-smokers which were part of the mitochondrial biogenesis/function (HSPD1, FEN1, COX18, COX10, UCP2 & 3), cellular senescence (PCNA, PTEN, KLOTHO, CDKN1C, TNKS2, NFATC1 & 2, GADD45A) and telomere replication/maintenance (PARP1, SIRT6, NBN, TERT, RAD17, SLX4, HAT1) target genes. Interestingly, NOX4 and TNKS2 were increased in the young IPF as compared to the young COPD patients. Genes in the mitochondrial dynamics and other quality control mechanisms like FIS1 and RHOT2 were decreased in young IPF compared to their age matched COPD subjects. ERCC1 (Excision Repair Cross-Complementation Group 1) and GADD45B were higher in young COPD as compared to IPF. Aging plays an important role in various infectious diseases. Elderly patients with chronic lung disease and smokers were found to have high incidence and mortality rates in the current pandemic of SARS-CoV-2 infection. Immunoblot analysis in the lung homogenates of smokers, COPD and IPF subjects revealed increased protein abundance of important proteases and spike proteins like TMPRSS2, furin and DPP4 in association with a slight increase in SARS-CoV-2 receptor ACE2 levels. This may further strengthen the observation that smokers, COPD and IPF subjects are more prone to COVID-19 infection. Overall, these findings suggest that altered transcription of target genes that regulate mitochondrial function, cellular senescence, and telomere attrition add to the pathobiology of lung aging in COPD and IPF and other smoking-related chronic lung disease in associated with alterations in SARS-CoV-2 ACE2-TMPRSS2-Furin-DPP4 axis for COVID-19 infection.
Introduction
Aging is an important factor influencing the overall lung health and function 1,2. Lung function declines with progress in age after lung maturation. Evidences suggest that a significant contribution of various environmental factors influence the ageing lung 3. According to the Behavioral Risk Factor Surveillance System (BRFSS) data from 2017, 6.2% (age-adjusted) US adults were reported to have chronic obstructive pulmonary disease (COPD) 4. Further, there has been an increasing reports of young COPD population visiting for different hospital services 5. Ageing influences many chronic lung diseases like COPD, idiopathic pulmonary fibrosis (IPF) and asthma. Also chronic lung diseases like asthma and IPF share some of the common yet distinct features compared to COPD 6. Environmental stress factors like smoking remains a common influencing factor for the disease progression in all these three cases.
Cigarette smoke (CS) is one of the strongest contributing risk factors in the pathogenesis of COPD along with the decline in lung function 7. Cellular senescence is a process of irreversible cell cycle arrest, having both beneficial and harmful effects depending on the cell state 8. CS plays a role in advancing the lung aging by altering the process of cellular senescence 9. Several factors including oxidative stress influence the process of cellular senescence. Telomeres and mitochondria play a major role in influencing the process of cellular senescence and are often associated with maintaining the lung health 8,10,11. Earlier reports from our laboratory and others have shown that smoking and COPD is associated with the mitochondrial damage and dysfunction, altering the process of cellular metabolism and function 12,13. Similarly, telomere dysfunction was also associated with smoking and COPD 14,15. Mitochondrial and telomere dysfunction also play a causative role in the progression of IPF 16-18.
Several molecular mechanisms were identified and reported in relation to CS causing COPD and associated complications 19. Cigarette smoke alters several key functions in the cells, among them the crucial genes related to mitochondrial function, cellular senescence and telomeric length were selected in the current study to observe for any differential changes among young and old age groups categorized as non-smokers, smokers and COPD groups. Our previous studies showed independent contributions of these canonical signaling pathways and how they contribute towards the development of premature lung aging in chronic lung diseases, such as COPD/emphysema 20-24. Accumulating evidence suggest the close relationship of all these three pathways in influencing lung ageing and disease 11,25,26. Senescent cells are found in many age-related/chronic diseases 27. Studies from our laboratory showed that mice from different age group when exposed to chronic air and CS influence the process of lung inflammation and senescence. Chronic CS exposure in lung epithelial cells and mice increases several markers of cellular senescence 9,28. Recently, it was reported that serum from COPD patients can induce senescence in lung epithelial cells 29, giving strength to the importance of this area that needs to be explored further. Several DNA damaging agents occur in smoke, which may activate of DNA damage response by influencing telomere function over the time and tend to accumulate senescent cells. However, recent meta-analysis suggests that even though smokers are associated with shorter telomere length, the study implicates that smoking does not accelerate the telomere attrition in leucocytes 30.
Smoking and COPD conditions altered the expressions of many key genes involved in the above three crucial pathways of cellular maintenance. The current study was undertaken to determine the changes in genes related to mitochondrial biogenesis and function, telomere function and cellular senescence with respect to their age in human lungs. This study is important in unraveling some of the potential biomarkers to differentiate and follow the course of ageing in COPD. Keeping in view the importance of these genes in both COPD and IPF, we have also made comparisons between the similar age grouped COPD and IPF subjects, using the same set of gene panels.
Aging is one of the key components, which decides the subject’s susceptibility to various diseases and infections 31. Inflammaging/immunosenescence defines the condition to combat various illness. Recent pandemic of Severe Acute Respiratory Syndrome (ARDS) Coronavirus 2 (SARS-CoV-2) was thought to affect more elderly people especially men 32-34. Men with age-related comorbidities has higher mortality rate during Coronavirus Disease 2019 (COVID-19)35. Comorbidities like cardiovascular disease, diabetes and chronic respiratory diseases present a high mortality rate36. The studies involving the role of lung aging and senescence play an important role in understanding the role of multiple players involved in combating against SARS-CoV-2 infection and can lead to potential druggable targets to combat it. Considering the important role played by some of the crucial receptors and targets in COVID-19 37, we determined the protein expression of SARS-CoV-2 receptor Angiotensin Converting Enzyme 2 (ACE2) and aiding proteases like transmembrane protease serine protease-2 (TMPRSS-2) and spike protein furin in the lung homogenates of non-smokers, smokers, COPD and IPF subjects. We have also determined for the levels of dipeptidyl peptidase 4 (DPP4), which acts as a receptor for the similar class of coronavirus which is Middle East Respiratory Syndrome Coronavirus, (MERS-CoV)38.
Materials and Methods
Scientific rigor and reproducibility
Rigorous and unbiased approach were used to ensure full and detailed reporting of both methods and analyzed data.
Ethical approval: Institutional biosafety an review board approvals
Ethics statement
The current study was approved for the procurement of the human lung tissues as de-identified tissues by the Materials Transfer Agreement and Procurement (Institutional Review Board), and laboratory protocols by the Institutional Biosafety Committee of the University of Rochester Medical Center, Rochester, NY. Patients’ data or patients are not directly involved in this study as the lung tissues were procured from several agencies (see below). All patients/subjects were of age 21 and above. All methods were carried out in accordance with relevant guidelines and regulations of the University of Rochester, Rochester, NY.
Human lung tissues
The human peripheral lung tissues from non-smokers, smokers and COPD were procured/obtained from the NDRI (National Disease Research Interchange; the samples were collected from patients with various cause of deaths reported such as cardiac arrest or trauma/accidents, for most of the samples lower peripheral lung lobes were used or as supplied), LTRC (Lung Tissue Research Consortium of the NHLBI) and Department of Medicine and Pathology, and of the University of Helsinki Hospital, Finland as described in our previous reports 39. The clinical characteristics of the subjects used in the current study are given in Table 1. The subjects were broadly classified into two age groups: young age (≤55 years) group, old age (>55 years) group, as per previous studies 40. Although there were some co-morbid conditions reported for the specimens from COPD patients (which were on various medications), the tissues were assigned to different groups based on the age, smoking status and lung disease status (normal vs smoker’s normal vs COPD) reported during the procurements of the specimens. Even though the definitions across various reports vary, the samples in the non-smokers group were either reported to have no smoking history or doesn’t fit under current smokers with a long past tobacco use specially in old non-smokers, but have normal lungs, smokers have current smoking history with lung function, and COPD have diseased lungs classified as COPD. Additional comparisons were made between COPD (8 additional samples from above groups were added in addition to the 24 samples mentioned in Table 1) and IPF lung samples (3 in young IPF, 13 in old IPF) based on their age to check for the changes in the same custom gene panel (Tables 2 and 3).
RNA isolation from lung tissues
Total RNA was extracted from the human lung tissues stored at -80°C or in the RNA later, using Direct-Zol RNA miniprep plus kit (Zymo research, R2071) according to the manufacturer’s instructions. RNA concentration was measured using a Nanodrop 1000 (Thermo fisher scientific, USA). Various genes involved in Mitochondrial Biogenesis and Function, Telomere Replication and Maintenance and Cellular Senescence pathways were included in the custom code sets (Supplementary Table S1). The code set contained a total of 112 genes including 6 reference genes (Abcf1, Hprt, Polr1b, Rplp0, Ldha, Gusb) for gene normalization (Supplementary Table S2). The samples were processed through the NanoString nCounter system (NanoString Technologies Seattle, WA, USA). A total of 400ng RNA was submitted after adjusting the samples to a minimum of 20µg/µl as per the requirements.
Validation of gene targets using quantitative real-time PCR
Selected mRNA, which were found to be significantly and differentially altered which were further validated for their expression using qPCR in samples (three samples per group) used in the study as described earlier 41. The primers were obtained from Bio-Rad, except for FEN1 (F-CACCTGATGGGCATGTTCTAC, R-CTCGCCTGACTTGAGCTGT) and 18S (F-GTAACCCGTTGAACCCCATT, R-CCATCCAATCGGTAGTAGCG), used as internal control were obtained from integrated DNA Technologies. Results were represented as pairwise comparisons to compliment the observations made by NanoString analysis. Student t-test was used to determine the level of significance between two groups, while ANOVA was used for multiple group comparisons.
Western blot analysis in human lung homogenates
Total protein isolated from the lung homogenates of non-smokers, smokers, COPD and IPF were reduced and separated using the pre-made polyacrylamide gels (Bio-Rad) 41. The transferred membranes were probed for some of the important protein involved in the COVID-19, like tmprss2 (ab92323), furin (ab183495), DPP4 (ab28340) and ACE2 (ab108252). All the antibodies used in the current study was procured from the Abcam. The relative expression and equal loading as assessed using the Ponceau S staining or β-actin after stripping of the blots.
Data processing and statistical analysis
Nanostring mRNA counts were first normalized using the NanoStringNorm function in the statistical analysis software R on log2 transformed data. Geometric mean of reference genes was used to remove the technique variation and background expression. Differential analysis was conducted using linear models in the limma package (R/Bioconductor) after adjusting for the gender difference. Comparison among different experimental groups were performed using linear contrasts within the linear model framework; moderated t statistics was used to determine the differences in the gene expression levels between groups with empirical Bayes approach. The Benjamini-Hochberg procedure was used to adjust the p values to control the false discovery rates at 5%. The analyzed data was represented in the graphs as y-axis showing the negative log10 P-value and x-axis representing log2 fold change across each pairwise comparisons as described previously 42. The significantly altered gene data were shown as dot plot representation in supplementary Figures S1-9. Four samples from each age group were used for comparisons among non-smokers, smokers and COPD groups (as given in Supplementary Table S2). Comparisons with IPF includes all the samples as mentioned in the Table 1.
Results
Overall the study consisted of 24 lung tissue samples from different sources as mentioned above. The collected tissues were classified into six different groups based on age, the smoking and disease status. Further, comparative gene analysis was also done based on the smoking and disease status irrespective of the age. There were no genes in common that were changed in any of the comparisons involving all the three groups i.e. non-smokers, smokers and COPD. However, the individual group wise comparisons were reported here.
Differentially expressed genes in young non-smokers versus young smokers versus young COPD groups
First, we analyzed differentially expressed transcript levels among young non-smokers vs. young smokers, young smokers vs. young COPD and young non-smokers vs. young COPD (Figure 1). We found 5 genes were differentially expressed in young non-smokers vs. young smokers’ pairwise comparison. Out of 5 genes, the transcript levels of 4 genes (NFATC1, NFATC2, GADD45A, and CDKN1A) were decreased and 1 gene (PARP1) was increased in the young smokers as compared to young non-smokers group (Figures 2 and 3 A). Next, we compared genes differentially expressed in young smokers vs. young COPD pairwise comparison. Out of 5 genes, transcript levels of 1 gene (SIRT6) was decreased and the remaining 4 genes (RAD17, CDKN1C, COX10 and KLOTHO) were significantly increased in young smokers as compared to young COPD group (Figures 2A and 3B). Finally, we found 6 genes differentially expressed among young non-smokers vs. young COPD pairwise comparison. Out of 6 genes, the transcript levels of 2 genes CDKN1C and KLOTHO that belong to cellular senescence panel were decreased in the young COPD as compared to young non-smokers group. While the transcript levels of remaining 4 genes PARP1, SIRT6, TERT and SLX4 were increased in young COPD as compared to young non-smokers group (Figures 2A and 3C). Overall, 4 genes PARP1, SIRT6, KLOTHO and CDKN1C were among the common target genes that were differentially expressed in young COPD as compared to young non-smokers and young smokers groups.
Differentially expressed genes in old non-smokers versus old smokers versus old COPD groups
Here we analyzed differentially expressed transcript levels among old non-smokers vs. old smokers, old smokers vs. old COPD and old non-smokers vs. old COPD groups. Out of 7 genes, we found 3 genes IGF1, COX18 and RIF1 were decreased and remaining 4 genes NFATC1, NFATC2, RAD17 and PCNA were increased in old smokers as compared to old non-smokers group (Figures 2B and 4A). The transcript levels of IGF1, PARP1, PTEN, NBN, HSPD1 and RIF1 were decreased and GAR1 was increased in old smokers as compared to old COPD group (Figures 2B and 4B). Only 2 genes were affected among old non-smokers and old COPD group; RPA2 and PCNA were increased in old COPD as compared to old non-smokers group (Figures 2B and 4C). Overall, a total of 3 genes as mentioned above (IGF1, RIF1 and PCNA) were among the common targets that was found differentially expressed in old smokers as compared to old non-smokers and old COPD groups.
Altered gene expression levels in young and old non-smokers versus smokers versus COPD groups
We then analyzed differentially expressed genes among young non-smokers vs. old non-smokers, young smokers vs. old smokers and young COPD vs. old COPD pairwise comparisons. Transcript levels across different age groups were performed to better understand, whether age factor influences the measured outcomes in the current study. Accordingly, we found that 9 genes were significantly elevated in young non-smokers as compared to old non-smokers group. The following are the genes that were increased among young non-smokers (PCNA, NFATC2, ACD, GSK3β, HAT1, UCP2, CDKN1A, CDKN1C and SIRT1) as compared to old non-smokers group (Figures 2C and 5A). Although, young smokers show increased transcript levels of PARP1, UCP3 and E2F1 genes but decreased levels of NFATC1, NFATC2, MYC and GADD45A as compared to old smokers group as seen in Figures 2C and 5B. Additionally, 2 genes TNSK2 and PTEN were decreased in young COPD as compared to old COPD group. The transcript levels of GAR1, TERT, H2AX and FEN1 tend to increase in younger COPD as compared to old COPD group (Figures 2C and 5C). Interestingly, we noted that transcript levels of NFATC2 was decreased in lungs of old non-smokers, but increased in lungs of old smokers.
Combined analysis of differentially expressed genes among non-smokers, smokers and COPD groups
We performed grouped analysis of differentially expressed transcripts (comparisons without considering the age factor (young or old) among non-smokers, smokers and patients with COPD groups. Data from young and old that belong to same experimental groups were combined (young nonsmokers with old non-smokers group, young smokers with old smokers group and young COPD and old COPD group; n=8/group, as given in Figures 6 and 7A). Results indicated that smokers show decreased FOXO1 and increased RAD17 levels as compared to non-smokers group (Figure 8A). While decreased PARP1 and increased RAD17 levels were observed in smokers as compared to COPD group (Figure 8B). In patients with COPD KLOTHO gene was decreased and PARP1 and SLX4 genes were increased as compared to non-smokers group (Figure 8C). Furthermore, these comparisons revealed that smokers have significantly higher levels of RAD17 expression compared to both non-smokers and COPD groups. Whereas, COPD patients showed significantly higher levels of PARP1 as compared to both non-smokers and smokers groups.
Altered gene expression levels in young and old non-smokers versus IPF groups
Pairwise analysis of young non-smokers vs young IPF showed 25 significantly altered genes, which were given in Table 2. Comparisons between old non-smokers and old IPF showed 13 significantly altered genes, as listed in table 3 along with their observed level of significance. The gene comparisons among these groups were given in Figure 7B.
Altered gene expression levels in young and old COPD versus IPF groups
As detailed above both COPD and IPF are chronic age related diseases that severely alter lung function and share certain common features for the disease occurrence and progression. Here, we compared the altered gene levels related to the same pathways among COPD and IPF subjects as detailed above. There was no change in any of the genes analyzed in comparisons between young IPF (n=3) and old IPF (n=13). While, 16 genes were found to be altered in the comparisons between young (COPD vs IPF), as indicated in Table 4. A total of 6 genes were altered in comparisons between old COPD vs old IPF) as shown in Table 5.
Gene expression analysis of differentially expressed targets
Some of the differentially expressed mRNA targets predicted using Nanostring were selected for qPCR study. As given in the Figure 9 the expression trends of these selected genes matches with the trend observed using the NanoString mRNA analysis with varied levels of fold changes. Combined gene analysis for PARP1 was also given, and matches with the trend are given in Figure 8. The results clearly indicate that the genes validated using qPCR are in agreement and significant across all the pairwise comparisons made.
Protein expression levels of crucial SARS-CoV-2 targets in the lung homogenates
Western blots analysis (Figures 10 and 11) revealed a significant increase in the protein levels of TMPRSS2 protease (which plays a crucial role in the processing of the SARS-CoV-2 proteins) in COPD subjects as compared to smokers and non-smokers. Similarly, the levels of another important protease furin a spike protein, was also increased in smokers and COPD subjects, with a significant expression in COPD subjects. ACE2, which is considered crucial for the SARS-CoV2 binding as a receptor, was found to increase in smokers as compared to the rest of the groups. The samples were also further probed for the expression of DPP4, another crucial protein which plays an important role in MERS-CoV binding. The DPP4 expression was found to be significantly higher in smokers as compared to COPD and non-smokers. These results indicate that the levels of proteases, which aid in the processing and binding of the viral spike proteins were highly expressed in COPD and smokers. The expression results showed varied protein intensities in smokers and COPD for TMPRSS2 and DPP4 in the lungs, which may suggest a varied effect of virus entry/susceptibility based on specific cells in smokers and COPD/IPF subjects.
Discussion
Aging is a major contributor for the decline in lung function and as most of the COPD are old aged, so there is a need to define the role of aging influence in COPD. On the other hand, smoking is one of the major contributing factors for the development of COPD. COPD, is commonly observed in the older subjects compared to younger ones. Nonetheless there are growing evidences of young subjects with COPD, which needs a thorough and careful phenotypic characterization for potential markers to understand the cause and progression of disease. The current study examined the changes in the gene expression in the normal lungs of non-smokers and smokers and patients with COPD/IPF. The study included age as an influencing factor apart from the lung disease status. The extracted RNA was processed and analyzed using the sophisticated NanoString nCounter analysis platform. The NanoString has several advantages in simultaneous estimation of the several gene levels in a single sample with low amounts of sample input 43,44. We have successfully used this platform to report the changes in gene levels using different gene set panels in our earlier studies 17,29. In the current study, the custom designed panel included genes from three different pathways addressing the mitochondrial biogenesis and function, telomere function and cellular senescence, which play a major role in the lung inflammation and COPD development.
Most of the chronic diseases like COPD are associated with the mitochondrial dysfunctions. Several reports claim the causal role of mitochondrial dysfunctions in the initiation and progression of smoking associated COPD 45-47. We have previously reported the presence of mitochondrial dysfunctions in CS-induced lung damage models and in human lungs 11. Further, we along with others have reported that telomere dysfunction is also seen in the COPD patients and smoking plays as crucial role in influencing the telomere genes 10,14,21. Assessment of all these three pathway related genes in the same subjects throws light on the involvement and coordination of complex processes in smoking-related chronic lung disease such as COPD.
Pairwise comparisons revealed that a total of 21 genes were altered among the young and old subjects. Among them some of the genes were earlier reported to be altered in smoker and COPD patients. CDKN1A (p21) which is a cell dependent kinase (CDK), plays a vital role in the cellular senescence and proliferation was reported to be increased in smokers and COPD subjects 48. Further, p21 disruption attenuated CS-induced lung inflammation in mice 49. In the current study, there was a reduced expression in the p21 levels in the young smokers compared to NS in accordance with the previous reports, where CDKN1C (p57), along with p21 was decreased in aging lungs of mice 50.
Among the other important targets KLOTHO, SIRT6, PARP1 and PCNA were altered in the current study. COPD patients has reduced KLOTHO expression as compared to non-smokers, in accordance with previous reports 51. In similar lines PARP1 was also elevated in the COPD subjects compared to smokers and non-smokers, as reported earlier 52,53. We have reported that TERT levels were altered in mice (young and old) exposed to chronic CS 20. Accordingly, the current results demonstrates that TERT levels were significantly increased in COPD patients, but this gene may be influenced in a different way in humans, as younger COPD patients has higher TERT levels compared to older ones. Nuclear factor of activated T cells c2 (NFATC2) levels were significantly higher in aged smokers as compared to younger ones. Earlier reports claim that NFATC2 levels were increased by nicotine/smokers 54. It was also reported that NFATC2 enhances tumor-initiating phenotypes in lung adenocarcinoma 55. These observations were opposite in non-smokers, as they age the levels of NFATC2 were decreased. This suggests that NFATC2 may be used as a potential marker related to CS-induced lung damage. Among the other important genes that were altered and are crucial in the maintenance of these three pathways are ACD, which is related to TPP1 gene and coordinates with its function was elevated in COPD 21,56. FEN1 could be a novel biomarker for COPD which was increased in young COPD as compared to old. Prior studies showed mutation in FEN1 linking lung cancer progression in an age-dependent manner in mice exposed to benzo[α]pyrene which is present in tobacco smoke 21,57.
Pairwise comparisons between non-smokers and IPF subjects revealed changes in some of the important genes like PARP1, PCNA, FEN1, CDKN1B, NFATC2 and GADD45B, as discussed in above comparisons. However, the directionality of the changes varied between groups, which may be attributed to the small sample size in the study and the heterogeneity in the samples used. In view of the growing interests, comparisons were also made between the expression profiles of the young IPF and old IPF with their age matched COPD subjects. Interestingly some of the well characterized genes in IPF like NOX4, TNKS2 was increased in the young IPF as compared to the young COPD patients 58,59. Mitophagy is a well-known phenomenon occurring in both COPD and IPF and regulates the mitochondrial related damage response towards the disease 45,60. Genes participating in the mitochondrial dynamics and other quality control mechanisms like FIS1 and RHOT2 were found to be decreased in young IPF compared to their age matched COPD subjects. ERCC1 (Excision Repair Cross-Complementation Group 1) was also found to be high in young COPD as compared to IPF. Earlier reports claim that ERCC1 gene is strongly associated with COPD subjects 61. Some of the common gene targets like GADD45B needs further attention and characterization especially in the chronic lung diseases like COPD and IPF. Recent biomarker identification study indicated GADD45B in their list of genes that can be targeted in chronic diseases like asthma, IPF and COPD 6.
Several studies indicate that the patients with underlying chronic disease conditions like diabetes, hypertension and COPD are more prone to COVID19 infections and have higher chances of the hospitalization rates and mortality 33,35,36. Several crucial mechanisms and targets were reported to increases these susceptibility in these patients 37. In the current study, we have determined the expression of four such important protein targets reported to play an important role in SARS-CoV-2 COVID-19. As anticipated, COPD and IPF patients in our study have higher levels of TMPRSS2 proteins in the lungs, suggesting the ideal condition for the processing of the viral protein and attachment to its receptor ACE2. ACE2 receptor abundance expression was slightly increased, but was tend to be on higher side in smokers, COPD and IPF subjects 62. Furin, another crucial protease in the COVID-19 infection was also found to be higher in smokers, COPD and IPF as compared to the non-smokers. We have also assessed the levels of DPP4 in the same samples and found that DPP4 was significantly higher in the smokers as compared to non-smokers and COPD. DPP4 was found to play an important role in the entry of MERS-CoV and acts as its receptor (belonging to the similar class of beta coronaviruses) in humans, was also suggested to play an important role in COVID-19 infections 63. In agreement with recent reports 38, suggesting that smokers and COPD have higher DPP4, our findings suggest increase in DPP4 in smokers, but to our surprise we didn’t find any increase in COPD subjects. This may suggest a different mechanism in smokers and COPD, which required further studies. Nevertheless, the varied abundance of different SARS-CoV-2 COVID-19 proteins suggest cell-specific effects for viral entry in smokers and COPD/IPF subjects.
In conclusion, our study provides the novel directions in conducting the studies involving the crucial and interdependent mitochondrial, telomere and cellular senescence pathways in the same subjects in association with SARS-CoV-2 COVID-19 proteins. Whilst the study provides several differential gene expression patterns in non-smokers, smokers and COPD/IPF, there is a limitation regarding the small sample size and smoking history. Nonetheless the study is valuable in terms of approach and identifying a good number targets, which needs thorough characterizations in future studies.
Data Availability
All data are available within the manuscript.
Funding
This study was supported by the NIH R01 HL1377380, R01 HL135613, and R01 ES 029177 (all IR). DL is supported in part by the University of Rochester CTSA UL1 TR002001 of the NIH.
Author contributions
KPM, IKS and IR designed and interpreted the experimental results. KPM performed all the experiments, DL analyzed the data and edited the manuscript, and prepared Figures 1-8. KPM prepared and edited Figures 1-11. KPM, IKS and IR edited, improvised and approved the final version of the manuscript.
Additional Information
Supplementary information
accompanies this paper (Supplementary Table S1, Supplementary Table S2, and Supplementary Figures S1-11).
Competing Conflict of Interests Statement
The authors declare no competing interests.
Data availability statement
We declare that we have provided all the data in the manuscript.
Supplementary Tables
Supplementary Table S1: Pre-selected genes belonging to the mitochondrial biogenesis and function pathway with their annotations
Supplementary Table S2: Determined genes with raw and normalized counts based on Nano String analysis.
Supplementary Figures S1-9. Dot plot representation of the significantly altered genes across various groups indicated. The data were represented as the normalized counts among all the observed groups. Supplementary Figure 10-11. Full unedited gels/blots for Figures 10 and 11 (original and unprocessed) are shown.