Next Article in Journal
Premature MicroRNA-Based Therapeutic: A “One-Two Punch” against Cancers
Previous Article in Journal
Clinical and Genetic Tumor Characteristics of Responding and Non-Responding Patients to PD-1 Inhibition in Hepatocellular Carcinoma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 12
Issue 12
10.3390/cancers12123834
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Germline Genetic Variants of the Renin-Angiotensin System, Hypoxia and Angiogenesis in Non-Small Cell Lung Cancer Progression: Discovery and Validation Studies

by
Maria Joana Catarata
1,2,3,4,5,*,
Rui Medeiros
4,5,
Maria José Oliveira
1,2,4,
Alice Pêgo
3,
João Gonçalo Frade
6,
Maria Fátima Martins
6,7,
Carlos Robalo Cordeiro
3,7,
Felix J F Herth
8,9,
Michael Thomas
9,10,
Mark Kriegsmann
9,11,
Michael Meister
9,12,
Marc A Schneider
9,12,
Thomas Muley
9,12,† and
Ricardo Ribeiro
1,2,6,13,†
1
i3S-Institute for Research & Innovation in Health, University of Porto, 4200-135 Porto, Portugal
2
INEB- Institute of Biomedical Engineering, University of Porto, 4200-135 Porto, Portugal
3
Department of Pulmonology, Coimbra Hospital and University Center, 3000-075 Coimbra, Portugal
4
Faculty of Medicine, University of Porto, 4200-319 Porto, Portugal
5
Molecular Oncology and Viral Pathology Group-Research Centre, Portuguese Institute of Oncology, 4200-072 Porto, Portugal
6
Department of Clinical Pathology, Coimbra Hospital and University Center, 3000-075 Coimbra, Portugal
7
Faculty of Medicine, University of Coimbra, 3004-504 Coimbra, Portugal
8
Department of Pulmonology and Critical Care Medicine, Thoraxklinik, University of Heidelberg, 69126 Heidelberg, Germany
9
Translational Lung Research Centre (TLRC), member of the German Centre for Lung Research (DZL), 69120 Heidelberg, Germany
10
Department of Thoracic Oncology, Thoraxklinik, University of Heidelberg, 69126 Heidelberg, Germany
11
Institute of Pathology, University Hospital Heidelberg, 69120 Heidelberg, Germany
12
Translational Research Unit, Thoraxklinik, University of Heidelberg, 69126 Heidelberg, Germany
13
Laboratory of Genetics and Institute of Environmental Health, Faculty of Medicine, University of Lisbon, 1649-026 Lisboa, Portugal
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this paper.
Cancers 2020, 12(12), 3834; https://doi.org/10.3390/cancers12123834
Submission received: 6 November 2020 / Revised: 10 December 2020 / Accepted: 14 December 2020 / Published: 18 December 2020
Versions Notes

Abstract

:

Simple Summary

The presence of polymorphic gene variants in the human genome provides extensive genetic (and eventually phenotypic) variation affecting both normal physiological mechanisms and cancer pathogenesis. Functional genetic polymorphisms might have predictive and/or prognostic value in lung cancer, opening novel opportunities to improve prediction and guide clinical reasoning and therapeutics in lung cancer patients. Recent knowledge pinpoints a pleiotropic role for renin-angiotensin system, particularly in the lung and mainly through locally regulated alternative molecules and secondary pathways. Dysregulation of this system play a role in cell proliferation, hypoxia and angiogenesis, which processes are involved in lung cancer progression. Here we suggest that polymorphic variants in genes coding for renin-angiotensin system might play a role in Non-Small Cell Lung Cancer progression.

Abstract

Introduction: The renin–angiotensin system (RAS) is involved in cell proliferation, immunoinflammatory response, hypoxia and angiogenesis, which are critical biological processes in lung cancer. Our aim was to study the association of putatively functional genetic polymorphisms in genes coding for proteins involved in RAS, hypoxia and angiogenesis with non-small cell lung cancer (NSCLC) prognosis. Methods: Genotyping of 52 germline variants from genes of the RAS and hypoxic/angiogenic factors/receptors was performed using MassARRAY iPLEX Gold in a retrospective cohort (n = 167) of advanced NSCLC patients. Validation of the resulting genetic markers was conducted in an independent group (n = 190), matched by clinicopathological characteristics. Results: Multivariate analysis on the discovery set revealed that MME rs701109 C carriers were protected from disease progression in comparison with homozygous T (hazard ratio (HR) = 0.5, 95% confidence interval (CI) = 0.2–0.8, p = 0.010). Homozygous A and T genotypes for KDR rs1870377 were at increased risk for disease progression and death compared to heterozygous (HR = 1.7, 95% CI = 1.2–2.5, p = 0.005 and HR = 2.1, 95% CI = 1.2–3.4, p = 0.006, respectively). Carriers of homozygous genotypes for ACE2 rs908004 presented increased risk for disease progression, only in the subgroup of patients without tumour actionable driver mutations (HR = 2.9, 95% CI = 1.3–6.3, p = 0.010). Importantly, the association of homozygous genotypes in MME rs701109 with risk for disease progression was confirmed after multivariate analysis in the validation set. Conclusion: This study provides evidence that MME polymorphism, which encodes neprilysin, may modulate progression-free survival in advanced NSCLC. Present genetic variation findings will foster basic, translational, and clinical research on their role in NSCLC.

1. Introduction

Lung cancer is one of the most common malignancies worldwide and the most common cause of cancer death, with over one million people diagnosed worldwide each year [1].
The renin–angiotensin system (RAS) is an established primary regulator of blood pressure, homeostasis, and natriuresis [2]. Recent evidences indicate that angiotensin peptides might have a role in tumour cell proliferation, inflammation, immune surveillance, hypoxia and angiogenesis, which have been recognized as hallmarks of cancer, and thereby actively participating in lung tumour microenvironment regulation [3,4]. The enzymatic cascade leading to Angiotensin II (Ang II) production includes renin, which cleaves angiotensinogen (AGT) to form angiotensin I (Ang I). Angiotensin-converting enzyme (ACE) cleaves Ang I to produce Ang II [2]. The chymase, which might catalyse the conversion from Ang I to Ang II, has been associated with activation of proliferation, adhesion, regulation of E-cadherin expression, and modulation of the immunosuppressive microenvironment in lung carcinoma [5,6]. Angiotensin-(1-7) (Ang-(1-7)) has opposite effects to Ang II through binding to Mas receptor (MasR), and is obtained by catalysis from the enzymes ACE type 2 (ACE2) or neprilysin (NEP) [7]. On its turn, NEP has been associated with some malignancies, such as prostate, renal and lung cancer [8].
In vitro and in vivo studies showed that activation of the ACE/Ang II/ angiotensin II type 1 receptor (ATR1) axis promotes cancer cell proliferation, invasion, angiogenesis and reduced immunosurveillance [9,10,11]. Conversely, anti-tumour effects were demonstrated for the ACE2/Ang-(1-7)/MasR axis [12,13,14]. Upregulated expression of ATR1 transcripts in lung tumours suggests that Ang II, acting upon the ATR1, plays an important role in cell growth and invasion [15], whereas decreased expression of genes encoding for ACE2/Ang-(1–7)/MasR axis proteins in lung tumours halt MasR signalling, which is contrary to Ang II/ATR1 [16]. Notably, clinical studies revealed that RAS inhibitors combined with chemotherapy and Epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors prolonged progression-free survival in non-small cell lung cancer (NSCLC) patients [17,18,19].
Activity and modulation of the RAS are dependable of germline genetic variations [20], somatic mutations [21], environmental stimulus and from hypoxic lung tumour microenvironment [22]. Additionally, the RAS is causally and consequently associated with hypoxia and angiogenesis, considered in continuum as branches of RAS [4]. Notwithstanding, the RAS remains largely unexplored at translational and clinical levels in cancer, particularly NSCLC. In this study, we sought to evaluate whether functional genetic variants coding for proteins of RAS and hypoxic/angiogenic factors/receptors have predictive and/or prognostic value in NSCLC using progression-free survival (PFS) and overall survival (OS) as clinical endpoints.

2. Material and Methods

2.1. Study Population: Discovery and Validation Sets

This study comprised a retrospective cohort of histologically confirmed NSCLC patients (n = 167) with clinical advanced stage at diagnosis (IIIA-IVB), between 2012 and 2018. Participants were consecutively recruited during routine follow-up visit, since August 2017 to October 2018 from Coimbra University Hospital (discovery set). The validation set comprised an independent cohort of NSCLC patients (n = 190) selected from the Lung Biobank Heidelberg, member of the Biomaterialbank Heidelberg (BMBH) and the biobank platform of the German Centre for Lung Research (DZL). All patients were recruited at Thoraxklinik (University Hospital of Heidelberg) between 2015–2018, using the same matched clinical and pathological criteria from the discovery set.
Histolopathogical diagnoses in both sets were made by experienced pathologists following the 2009 and 2015 WHO classification criteria for lung cancer, accordingly to the year of patient’s diagnosis. The existence of concomitant primary tumours in another organ was considered as exclusion criteria for both cohorts. Clinical data was retrieved from clinical charts on pathological background, previous medication, Eastern Cooperative Oncology Group performance status (ECOG PS) at diagnosis, TNM staging according to 8th classification, tumour mutational status, type of cancer treatment, disease progression/death by clinical authors. Informed consent was obtained from each participant and all procedures are in agreement with Helsinki Declaration.
In the discovery set, somatic mutations in EGFR, ALK and ROS1 were determined using fluorescence in situ hybridization (FISH), in patients with adenocarcinoma stage IVA-IVB, whereas for the validation set, somatic mutations were determined through panel-based next generation sequencing in advanced NSCLC patients [23]. In both cohorts, targeted therapies were administered to carriers of actionable driver mutations, whereas checkpoint inhibitors were used as second-, or third-line therapy for patients with PD-L1 expression on tumour cells below 50%.
The primary endpoint was PFS and the time-to-disease progression was calculated in months, from the date of first line chemotherapy/target therapy, until the date of progression according to RECIST criteria. OS was included as secondary endpoint, and the time-to-death was computed in months from the date of first line chemotherapy/target therapy until the date of death/date of last visit.

2.2. Selection of Genetic Polymorphisms

Fifty-two single nucleotide polymorphisms (SNP) included in the present study from the genes ACE, ACE2, AGT, AGTR1, AGTR2, MME, CMA1, MAS1, HIF1A, VEGFA, KDR, PGF and FLT1 (Table S1) were selected from search in genomic public databases, in silico analysis and review of scientific literature, to identify putatively functional polymorphisms [24,25,26,27,28,29,30,31,32,33,34,35,36,37]. Polymorphisms with minor allele frequency below 1% were excluded [37].

2.3. Genotyping

2.3.1. Discovery Set

Each patient donated a sample of blood (~4 mL) for research, collected to EDTA-Vacutainer tubes, at the same time of blood collection for routine analytic follow-up. The collected blood was centrifuged to obtain the buffy coat, which was used to isolate and purify DNA by the EZ1 DNA Blood kit in the EZ1 BioRobot (QIAgen, Hilden, Germany). The selected SNPs were genotyped using the Sequenom Mass ARRAY matrix-assisted laser desorption/ionization time-of-flight mass spectrometry platform (Sequenom, San Diego, CA, USA). Primers were designed using semi-automated Assay Design 3.1 Software (Sequenom). The global genotyping success rate was 99.85% with 100% replication agreement among one third of the samples.

2.3.2. Validation Set

Buffy coat samples were provided by Lung Biobank Heidelberg from patients diagnosed with advanced NSCLC (clinical stage IIIA-IVB). DNA extraction and purification were conducted using the FlexiGene DNA Kit (QIAgen). The resulting genetic polymorphisms, associated with disease progression or overall survival from multivariate analysis in the discovery set, were genotyped in the validation set using Real-Time PCR. Allelic discrimination was obtained using TaqMan SNP genotyping assays: ACE2 rs908004 (Assay ID: C___8816953_20), MME rs701109 (Assay ID: C____964850_10) and KDR rs1870377 (Assay ID: C__11895315_20) (all from Thermo Fisher Scientific, Waltham, MA, USA). Quality control included non-template controls in all runs and blind replicate assessment in 5% of the samples. Genotype was assigned for 96.3% of all samples tested.

2.4. Statistical Analysis

Descriptives were depicted as average ± standard deviation or median (interquartile range) according to departure from normality (Shapiro-Wilk test). Central tendency data was compared using the t-test, Mann-Whitney or Kruskal-Wallis tests according to normality and stratification of independent variables. Categorical variables were expressed as frequency (percentage) and proportions tested using the chi-square test. The difference in median months of PFS and OS from both cohorts was calculated through Kaplan-Meier curve with log rank test.
Four genetic models were considered in this study: each SNP was stratified according to the wild type allele into additive (“AA” genotype versus “Aa” versus “aa”), dominant (“AA” genotype versus “Aa” and “aa” genotypes), recessive (“Aa” and “AA” versus “aa” genotype) and heterosis (“Aa” versus “AA” and “aa” genotypes). The wild-type allele for each SNP was designated in agreement with Ensembl [37]. Empirical analyses were conducted on clinicopathological and genetic variables to test the association with time-to-outcome, using Cox regression analyses. After univariate analyses, p < 0.10 and an effect > 30% in hazard ratio (HR) were used as criteria for inclusion in a multivariated Cox with stepwise regression model. Genetic models with under-represented variant genotype (n < 4) were excluded from the multivariate analysis, and only the genetic model with the highest likelihood ratio for each SNP was included in the multivariate model to avoid collinearity. The resulting multivariated model was further tested with bootstrapping (MonteCarlo simulations, 1000 iterations), to elect the resulting relevant markers to be tested in the validation set. To investigate the possible predictive value of each SNP in treatment response for first line therapy, and time-to-progression or time-to-death, exploratory analyses were conducted examining the effect of genotype in a subgroup of patients treated with first line platinum doublet chemotherapy and in another subgroup of wild-type EGFR (epidermal growth factor receptor) patients, using the same survival analysis methods. Sample size, power, effect size and number of events for survival analyses (assuming alpha = 0.05 and power > 0.8) were estimated in both cohorts. Statistical analyses were performed on STATA software version 16.0 (StataCorp, College Station, TX, USA).

2.5. Ethics Approval

This project has been reviewed and approved by Coimbra University Hospital’s Ethical Committee (ref. no. 0111/CES) and by the Portuguese National Committee for data protection (ref. no. 2588/2017). The biobank protocol of Lungbiobank Heidelberg was approved by the ethics committee of the University of Heidelberg (S-270/2001).

3. Results

The clinicopathological characteristics of participating subjects from discovery and validation sets are described in Table 1. Regarding mutational status from discovery set (n = 167), we observed that 10.8% of patients (n = 18) had EGFR mutation (exon 19 deletions or exon 21 mutation), whereas 3.0% (n = 5) had rearrangements in the gene encoding anaplastic lymphocyte kinase (ALK). Mutations of ROS1 were not detected. From the validation set (n = 190), 15.3% (n = 29) of NSCLC patients had EGFR mutation (exon 19 deletions or exon 21 mutation), 13.7% (n = 26) were KRAS mutated, 4.2% (n = 8) had EML4-ALK fusion, 2.6% (n = 5) had BRAF V600E mutation. No other targetable mutations were identified. Concerning respiratory co-morbidities from discovery set, twenty patients (12.0%) had moderate-severe COPD with lung emphysema and one patient (0.6%) had sarcoidosis, whereas in the validation set there were fifty patients (26.3%) with moderate-severe COPD and lung emphysema, and five patients (2.6%) had idiopathic pulmonary fibrosis.
The description of clinical-pathological variables is depicted in Table S2 with report of univariate analyses for PFS and OS. Genotypic distribution of functional SNPs in genes of the renin-angiotensin pathway, hypoxic and angiogenic factors using additive, dominant, recessive and heterosis models from discovery set are described in Table S3, together with hazard and survival univariate analyses. The univariate analysis showed that carriers of allele T (n = 112, 67.1%) for the SNP rs4316 in ACE gene had a significant protection to disease progression (Hazard ratio (HR) = 0.6, 95% confidence interval (CI) (0.4–0.9), p = 0.011). Homozygous GG/CC for ACE2 rs908004 (n = 151, 90.4%) had increased risk for disease progression (HR = 2.1, 95% CI (1.1–3.9), p = 0.021). Another SNP from MME gene that encodes for neprilysin, showed a tendency between the dominant genetic model of MME rs701109 and PFS, having the heterozygous CT and homozygous CC protection to disease progression (HR = 0.6, 95% CI (0.3–1.0), p = 0.064). Univariate analysis also showed a significant association between the functional genetic variant of VEGFA rs25648 with time to tumour progression, being the homozygous C and T (n = 127, 78.9%) associated with worse prognosis (HR = 1.7, 95% CI (1.1–2.6), p = 0.030). The data showed a tendency between the heterosis genetic model of KDR rs1870377 and time to tumour progression, having the homozygous A and T increased risk to disease progression (HR = 1.4, 95% CI (1.0–2.0), p = 0.054) and to death (HR = 1.7, 95% CI (1.0–2.7), p = 0.034).
Multivariate analysis included the relevant variables (both clinical-pathological and genetic) from univariate empirical analysis (Table 2, Table 3). The C carriers (n = 149, 89.7%) for MME rs701109 had significantly decreased risk for disease progression than homozygotes T (n = 17, 10.3%), regardless of relevant clinical-pathological variables, namely tumour size, distant metastasis, type of systemic therapy, treatment modality or from other relevant SNPs (Table 2). Homozygous A and T for KDR rs1870377 (n = 99, 59.3%) had significantly lesser PFS and OS compared to heterozygous (n = 68, 40.7%), regardless of relevant clinical and pathological variables (Table 2 and Table 3). Notably, these results maintained its significance after boostrapping, as a statistical method for internal validation.
Exploratory subgroup analyses were conducted among NSCLC patients from discovery set that included only patients under first line platinum-based doublet chemotherapy without actionable driver mutations (n = 145) or EGFR wild-type (n = 149). Within the MME (rs701109) and KDR (rs1870377) polymorphisms, we found evidence for an independent effect at PFS and OS respectively, across all subgroups. Notably, CT and CC patients for MME (rs701109) under first line platinum-based doublet chemotherapy or EGFR wild-type had protection for disease progression (HR = 0.4, 95% CI (0.2–0.7), p = 0.005 and HR = 0.4, 95% CI (0.2–0.8), p = 0.006, respectively). Homozygous A and T for KDR rs1870377 had increased risk to death also in patients under first line platinum doublet chemotherapy or EGFR wild-type (HR = 1.9, 95% CI (1.1–3.2), p = 0.019 and HR = 2.0, 95% CI (1.2–3.4), p = 0.011, respectively). Moreover, from the multivariate analysis conducted in the subgroup of patients treated with first line platinum-based doublet chemotherapy without actionable driver mutations (n = 145), we also found that homozygous for ACE2 rs908004 had increased risk for disease progression, regardless of relevant clinical-pathological covariates (HR = 2.9, 95% CI (1.3–6.3), p = 0.010).
The validation set included an independent matched cohort of NSCLC patients (Table 1) to validate the resulting genetic SNPs from the discovery set. Results from univariate analysis of clinicopathological variables in association with PFS and OS are described in Table S4. The genotype distribution and univariate analysis according to additive, dominant, recessive or heterosis models of MME rs701109, KDR rs1870377 and ACE2 rs908004 for hazard and survival analysis are depicted in Table S5. The genetic distribution of genotypes in MME rs701109, KDR rs1870377 and ACE2 rs908004 are similar among discovery and validation sets (p = 0.204, p = 0.846 and p = 0.251, respectively) (Pearson Chi-square test).
The multivariate model for PFS in the validation cohort is presented in Table 4. Homozygous (n = 102, 53.9%) for MME rs701109 had significantly increased risk to disease progression compared to CT genotype (n = 87, 46.0%), regardless of histology, type of systemic therapy or treatment modality (chemotherapy or chemo and radiotherapy). There was no significant effect on overall survival for MME rs701109, KDR rs1870377 and ACE2 rs908004.
Subgroup analyses were also conducted in this set, for first line platinum-based doublet chemotherapy (n = 155) and EGFR wild-type patients (n = 161). The MME rs701109 polymorphism was associated with PFS only in the subgroup of patients with wild-type EGFR, being the homozygous at increased risk to disease progression (HR = 1.7, 95% CI (1.1–2.5), p = 0.009).

4. Discussion

The work developed here sought to test the hypothesis that variants in genes coding for the renin–angiotensin system and for interacting hypoxia/angiogenesis key molecules might have predictive and prognostic value, respectively for NSCLC progression and survival. SNPs are the most common type of polymorphisms, accounting for 90% of genetic variations [38]. Nucleotide substitution, according to its location in germline DNA, may elicit functional modifications, ultimately at protein level, contributing to inter-individual variability, and correlating with relevant phenotypic alterations. The association between genetic polymorphisms and cancer susceptibility, progression, aggressiveness and death has been thoroughly described in the literature [38].
RAS was first discovered and related with the physiological regulation of systemic arterial pressure and involved in hypertension pathogenesis [39]. As a gatekeeper system, RAS is responsible for more than one function, including cell cycle regulation, inflammation, angiogenesis and fibrosis [3]. The capillary blood vessels in the lung are one of the major sites of ACE expression and Ang II production in the human body [40]. In fact, the increased ACE expression observed in several interstitial lung diseases supports pulmonary RAS and establishes a putative role for Ang II in the response to lung injury and fibrosis [4]. A previous report showed that the G-6A polymorphism of the AGT gene is associated with increased angiotensin production and idiopathic pulmonary fibrosis progression [41], which is associated with increased risk of lung cancer [42]. Furthermore, the insertion/deletion (I/D) polymorphism of ACE gene has been linked to the pathogenesis and progression of human cancers [43,44] and the lowest serum ACE levels were correlated with poorer prognosis and higher relapse rates in lung cancer [15]. We found that the association between ACE rs4316 genetic variant and PFS, despite being significant in univariate, the finding was not confirmed on multivariate analysis. It is an I/D polymorphism and is localized at intron 16 of ACE gene. This SNP has been associated with plasma ACE levels, with DD genotype carriers independently associated to higher expression [45], albeit no effect was noted on AngII [45]. Thus, it seems unlikely this SNP might impact the ACE/AngII/ATR1 axis, which has been associated with tumour cell growth and angiogenesis, partially justifying the absence of clinical impact in NSCLC patients with advanced stage disease. There are other possible reasons for the failure to replicate the prior findings. Many of the initial reports were present only in retrospectively defined subsets of the subjects studied and the main hypotheses investigated were related to cancer risk. Another possibility is that varied ethnic background of study populations may result in a different genetic background which confers a different biologic significance to a given functional polymorphism.
ACE2 is responsible for converting AngII to Ang-(1-7), an agonist of Mas receptor (MasR). In a mouse lung xenograft model, it was previously shown that ACE2 overexpression inhibited both lung cancer cell proliferation In vitro and reduced tumour growth in vivo [46]. The ACE2 rs908004 is an intronic variant at the X-linked ACE2 gene [47]. Findings from the discovery cohort revealed that ACE2 rs908004 heterozygous carriers presented shorter time-to-progression, only in the subgroup of patients treated with first line platinum-based chemotherapy without actionable driver mutations. We hypothesize that heterozygous produces reduced levels of ACE2, contributing to the up-regulation of ACE/Ang II/AT1 receptor axis, or to linkage with another genomic region that might be relevant for this clinical phenotype. These hypotheses need to be confirmed in further studies. Furthermore, as with all gene association studies, subgroup analyses can lead to spurious associations and might justify our result. Nevertheless, there is lack of studies on genotype–phenotype relationship with this SNP in lung cancer, opening new opportunities for further investigations.
Ang-(1-7) has opposite effects to Ang II upon binding to MasR and might be produced through hydroxylation of Ang I by neprilysin (NEP), an enzyme encoded by the gene MME [48]. NEP, also called CD10, has been described as a tumour-specific antigen for lymphocytic leukemia [49]. In the later years, it was demonstrated that NEP functions are implicated in oncogenesis and tumour microenvironment [49]. Moreover, alterations in NEP protein expression have also been reported in solid tumours, including lung cancer [50]. Both in discovery and validation cohorts, NSCLC subjects carrying the heterozygous TC genotype were at decreased risk for disease progression, albeit in the discovery set the protection was conferred to carriers of the TC/CC combined genotypes (dominant model), whereas in the validation cohort it was conferred only to TC (heterosis model). Heterosis for a single gene refers to a heterozygote that has a different phenotypic trait than for either allele homozygotes [51], and heterozygosity has been associated with human healthy aging [52]. Despite insufficient knowledge on the molecular basis of heterosis, an accepted mechanism is related to environmentally-driven epigenetic modifications, namely methylation and microRNAs (miRNAs) [53]. MicroRNAs regulate gene expression by binding to target sequences mostly at the 3’-UTR of mRNAs, and evidence exist for their functional effect in candidate SNPs of various types of cancer [54]. Noteworthy, recently reported data demonstrated such a mechanism with functional repercussion, the SNP rs6665, located downstream of rs701109 in the 3’-UTR region sequence of MME [55].
We conducted an in-silico proof-of-principle analysis using the MicroSNiPer database [56,57] to predict MME rs701109 SNP effects on putative microRNA targets, followed by study of the candidate miRNAs using the RNAhybrid program [58,59] to calculate the minimum free energy (MFE) of hybridization between miRNAs and target sequence (3′-UTR sequence of MME extracted from Ensembl [37]. The strongest binding affinity was found for hsa-miR-362-3p, which overlaps the locus of rs701109 at 3′-UTR region. A stronger predicted paring affinity between the hsa-miR-362-3p miRNA and the target mRNA of MME was observed for C-allele (−26.1 kcal/mol) in comparison to T-allele (−23.2 kcal/mol), stronger interaction with the miR-362-3p, resulting in lower expression of MME, and consecutively decreased risk for disease progression. Notably, previous data confirmed that hsa-miR-362-3p is expressed in lung tissue, aberrantly expressed in lung adenocarcinoma, and a promoter of metastasis [60] The differences among cohorts, particularly due to distinct therapeutics, genotype distribution and environmental stimuli, might account for different epigenetic regulation. Experimental studies in lung cancer tissues and In vitro with lung cancer organoids are warranted to further validate current findings.
The MME rs701109 variant consists in a nucleotide substitution of T-by-C, that has been shown to associate with lower NEP protein expression in the brain of patients with Alzheimer’s disease [61]. Notably, previous studies revealed that higher NEP expression associated with poorer prognosis in NSCLC and that tumour NEP was an independent predictor of recurrence in lung adenocarcinomas [50,62]. Taking into consideration that increased NEP expression in lung tumours associates with worst prognosis, and that TC and CC genotype carriers secrete less protein, our findings showing an independent protection for disease progression in TC carriers likely reflect the influence of this SNP in decreasing NEP availability. This result strengthens the involvement of NEP in NSCLC and supports the putative interest of including this genetic biomarker in a multilocus model for predicting relapse. Nevertheless, despite the known ability of NEP to degrade mitogenic peptides, the lack of Ang II degradation, as well as reduced production of Ang-(1-7) by NEP, might account for the growth of lung cancer cells [15]. Taken together, accumulated data and current findings support the need to foster research in neprilysin actions on cancer, particularly in lung cancer.
Hypoxia is a hallmark of solid tumours, and consequently related to enhanced aggressiveness, risk of disease progression and therapy resistance, effects known to be largely mediated by hypoxia-inducible factors [63]. The downstream upregulated vascular endothelial growth factor (VEGF) further contributes to cancer growth and metastasis, directly targeting tumour cells. Increased VEGF expression and/or circulating VEGF levels have been consistently reported in lung cancer [64]. Noteworthy, the activation of RAS in the tumour microenvironment usually involves Ang II, which has been described to induce hypoxia through reactive oxygen species production, and subsequent activation of pro-inflammatory and pro-angiogenic signals [65]. Albeit no clinical impact was demonstrated for polymorphisms in HIFA, PGF and FLT1 genes, associations were observed in the discovery set for VEGFA rs25648 with PFS (univariate analysis) and KDR rs1870377 with PFS and OS (multivariate analysis). Koukourakis et al. [66], described that inherited VEGF sequence variations were strong determinants of the molecular VEGF and VEGF-downstream phenotype of NSCLC, being the substitution C > A at locus −2578 associated with higher levels of VEGF. Another study showed that VEGFA rs25648 (−73C > T) was a prognostic biomarker for PFS and OS in patients with hepatocellular carcinoma, with C-allele conferring improved prognosis, despite not confirmed in multivariate analysis [67]. In the discovery cohort, lung cancer patients had a different outcome since the homozygous C for VEGF rs25648 had lower time until disease progression in univariate analysis. The apparently conflicting results might be due to particularities related with different cancer models. Notably, studies investigating the potential association between VEGF polymorphisms and circulating VEGF levels also led to controversial results [68,69,70], opening new opportunities for further investigations to clarify the role of germline VEGF variations in NSCLC progression and its repercussion in VEGF levels.
The KDR rs1870377 was significantly associated with PFS and OS in the discovery set. Homozygous T and C had increased risk for tumour progression PFS and death, independent of clinicopathological variables and after bootstrap analysis. Findings were not replicated in the validation cohort. KDR encodes VEGFR-2, the key receptor mediating endothelial cell proliferation and sprouting, leading to angiogenesis and tumour promotion upon binding and signalling through the VEGFA pathway [71]. The KDR rs1870377 (c.1719A > T) variant induces a missense alteration in exon 11 that results in the substitution A > T and in the amino acid modification Gln472His. This substitution was previously described to alter the function of VEGF-KDR signalling pathway, enhancing ligand/receptor binding efficacy and upregulated signal transduction, ultimately leading to increased microvessel density [72]. This phenotypic correlation with the KDR rs1870377 might support T homozygous findings, hypothetically by promoting tumour vascular density for homozygous carriers. Previous findings reported that microvessel count, reflecting tumour neoangiogenesis, appeared to be a poor prognostic factor for survival in NSCLC patients [73], despite some conflicting results [74].
This study provides evidence that SNPs coding for molecules involved in the RAS (ACE2 rs908004, MME rs701109, KDR rs1870377) might influence PFS and OS in advanced NSCLC patients, although only MME rs701109 was validated in an independent cohort. Candidate SNP studies searching associations with clinically relevant endpoints, have been reported to yield false-positive findings. Here, the multivariate model was first developed in a cohort of consecutively recruited advanced NSCLC patients. Furthermore, we included a robust data analyses setting that included: (i) sample size and number of events in both cohorts were within adequated estimations, (ii) initial empirical analyses to include only the most relevant variables in the multivariate models, (iii) bootstrapping using Montecarlo simulations with 1000 iterations to strengthen internal validation in each cohort, and (iv) validation of important SNPs from the discovery set in an independent cohort of lung cancer patients matched for demographic and pathological characteristics. Nevertheless, despite matching for clinicopathological characteristics, differences between discovery and validation cohorts existed, and this could partially explain the validation of a single polymorphism. Particularly, the validation set was selected from a larger lung cancer database and somatic mutational analyses were determined using distinct methods, revealing different mutational status for both cohorts, which likely influenced the course of disease and the subsequent therapeutic decisions. In addition, differences in respiratory co-morbidities, which may have conducted to mismatched treatment modalities (surgery with adjuvant chemotherapy versus chemo with radiotherapy for stages IIIA–IIIB), combined with the genetic background may also add to explain failure to validate the other genetic polymorphisms.

5. Conclusions

Germline genetic variants are common genomic events, putatively inducing phenotypic modifications affecting both normal physiology and cancer physiopathology. Data presented here supports the rationale for NEP and RAS mechanistic involvement in NSCLC, whereas the polymorphism MME rs701109 might prove to be a useful biomarker for NSCLC progression. These findings will foster further and larger studies to confirm the role of RAS, and particularly NEP in NSCLC, as well as to determine the usefulness of MME rs701109 as a molecular marker for disease progression. Taken together, these results can lead to a better understanding of NEP influence in the NSCLC biological mechanisms, contributing towards the definition of a genetic profile for this disease and eventually serve as guidance for therapeutic decisions.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/12/3834/s1, Table S1: Representative single nucleotide polymorphisms (SNPs) of genes involved in RAS (ACE, ACE2, AGT, AGTR1, AGTR2, MME, CMA1, MAS1) and SNPs of genes that express hypoxia-inducible factor 1-alpha, hypoxia-inducible factor 1-beta, hypoxia-inducible factor 2-alpha, hypoxia-inducible factor 2-beta, vascular endothelial growth factor and its receptor, placental growth factor and its receptor, Table S2. Clinical-pathological variables and progression-free survival (PFS) and overall survival (OS) in discovery set (n = 167)—Univariate analysis, Table S3: Discovery Set (n = 167): progression-free survival (PFS) and overall survival (OS)—Univariate analysis, Table S4. Clinical-pathological variables and progression-free survival (PFS) and overall survival (OS) in validation set (n = 190)—Univariate analysis, Table S5. Validation Set (n = 190)—Univariate analysis.

Author Contributions

Conceptualization, M.J.C., R.M. and R.R.; Data curation, M.J.C., A.P., J.G.F. and T.M.; Formal analysis, M.J.C., J.G.F., M.F.M., M.T., M.K., M.M., M.A.S., T.M. and R.R.; Funding acquisition, M.J.C. and C.R.C.; Investigation, M.J.C., C.R.C., F.J.F.H., M.T., M.M. and R.R.; Methodology, M.J.C., M.K., M.M., T.M. and R.R.; Project administration, M.A.S. and T.M.; Resources, M.J.O., A.P., J.G.F., M.F.M., F.J.F.H., M.K., M.M., M.A.S. and T.M.; Supervision, R.M., M.J.O., C.R.C., M.M., T.M. and R.R.; Validation, T.M.; Writing—original draft, M.J.C.; Writing—review & editing, M.J.O., T.M. and R.R. All authors have read and agreed to the published version of the manuscript.

Funding

M.J.C. was supported by the Associação de Estudos Respiratórios and the Portuguese Pulmonology Society.

Acknowledgments

M.J.C. would like to dedicate this work in the lovely memory of Maria João Matos, eminent senior pulmonologist, who was a continuous support and inspiration throughout this study and during her Pulmonology residency. Authors acknowledge Margarida Lourenço and technical support from Elisabete Camilo, Isabel Marques and Andreia Coelho in DNA extraction. Authors also acknowledge to Liz Meister and Martin Fallenbüchel for technical assistance and Maria Torres for genotyping using MASS ARRAY iPLEX Gold.

Conflicts of Interest

All authors declare that they have no conflict of interest.

References

  1. Ferlay, J.; Colombet, M.; Soerjomataram, I.; Mathers, C.; Parkin, D.M.; Piñeros, M.; Znaor, A.; Bray, F. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int. J. Cancer 2019, 144, 1941–1953. [Google Scholar] [CrossRef] [Green Version]
  2. Lavoie, J.L.; Sigmund, C.D. Minireview: Overview of the Renin-Angiotensin System—An Endocrine and Paracrine System. Endocrinology 2003, 144, 2179–2183. [Google Scholar] [CrossRef] [Green Version]
  3. Wegman-Ostrosky, T.; Soto-Reyes, E.; Vidal-Millán, S.; Sánchez-Corona, J. The renin-angiotensin system meets the hallmarks of cancer. J. Renin-Angiotensin-Aldosterone Syst. 2015, 16, 227–233. [Google Scholar] [CrossRef] [Green Version]
  4. Catarata, M.J.P.; Ribeiro, R.; Oliveira, M.J.; Robalo-Cordeiro, C.; Medeiros, R. Renin-Angiotensin System in Lung Tumor and Microenvironment Interactions. Cancers 2020, 12, 1457. [Google Scholar] [CrossRef]
  5. Takai, S.; Jin, D.; Miyazaki, M. New approaches to blockade of the renin-angiotensin-aldosterone system: Chymase as an important target to prevent organ damage. J. Pharmacol. Sci. 2010, 113, 301–309. [Google Scholar] [CrossRef] [Green Version]
  6. Jiang, Y.; Wu, Y.; Hardie, W.J.; Zhou, X. Mast cell chymase affects the proliferation and metastasis of lung carcinoma cells in vitro. Oncol. Lett. 2017, 14, 3193–3198. [Google Scholar] [CrossRef] [Green Version]
  7. Ferrario, C.M.; Trask, A.J.; Jessup, J.A. Advances in biochemical and functional roles of angiotensin-converting enzyme 2 and angiotensin-(1–7) in regulation of cardiovascular function. Am. J. Physiol. Circ. Physiol. 2005, 289, H2281–H2290. [Google Scholar] [CrossRef] [Green Version]
  8. Papandreou, C.N.; Usmani, B.; Geng, Y.; Bogenrieder, T.; Freeman, R.; Wilk, S.; Finstad, C.L.; Reuter, V.E.; Powell, C.T.; Scheinberg, D.; et al. Neutral endopeptidase 24.11 loss in metastatic human prostate cancer contributes to androgen-independent progression. Nat. Med. 1998, 4, 50–57. [Google Scholar] [CrossRef]
  9. Imai, N.; Hashimoto, T.; Kihara, M.; Yoshida, S.-I.; Kawana, I.; Yazawa, T.; Kitamura, H.; Umemura, S. Roles for host and tumor angiotensin II type 1 receptor in tumor growth and tumor-associated angiogenesis. Lab. Investig. 2006, 87, 189–198. [Google Scholar] [CrossRef]
  10. Jethon, A.; Puła, B.; Piotrowska, A.; Wojnar, A.; Ryś, J.; Dzięgiel, P.; Podhorska-Okolow, M. Angiotensin II Type 1 Receptor (AT-1R) Expression Correlates with VEGF-A and VEGF-D Expression in Invasive Ductal Breast Cancer. Pathol. Oncol. Res. 2012, 18, 867–873. [Google Scholar] [CrossRef] [Green Version]
  11. Lewandowska, U.; Lachowicz-Ochędalska, A.; Domińska, K.; Kaszewska, D.; Rębas, E. Angiotensin II as a factor modulating protein tyrosine kinase activity in two breast cancer lines—MCF-7 and MDA-MB-231. Endokrynol. Polska 2011, 62, 151–158. [Google Scholar]
  12. Li, H.; Qi, Y.; Li, C.; Braseth, L.N.; Gao, Y.; Shabashvili, A.E.; Katovich, M.J.; Sumners, C. Angiotensin type 2 receptor-mediated apoptosis of human prostate cancer cells. Mol. Cancer Ther. 2009, 8, 3255–3265. [Google Scholar] [CrossRef] [Green Version]
  13. Krishnan, B.; Torti, F.M.; Gallagher, P.E.; Tallant, E.A. Angiotensin-(1-7) reduces proliferation and angiogenesis of human prostate cancer xenografts with a decrease in angiogenic factors and an increase in sFlt-1. Prostate 2013, 73, 60–70. [Google Scholar] [CrossRef]
  14. Menon, J.; Soto-Pantoja, D.R.; Callahan, M.F.; Cline, J.M.; Ferrario, C.M.; Tallant, E.A.; Gallagher, P.E. Angiotensin-(1-7) Inhibits Growth of Human Lung Adenocarcinoma Xenografts in Nude Mice through a Reduction in Cyclooxygenase-2. Cancer Res. 2007, 67, 2809–2815. [Google Scholar] [CrossRef] [Green Version]
  15. Gallagher, P.E.; Cook, K.L.; Soto-Pantoja, D.; Menon, J.; Tallant, E.A. Angiotensin peptides and lung cancer. Curr. Cancer Drug Targets 2011, 11, 394–404. [Google Scholar] [CrossRef] [Green Version]
  16. Goldstein, B.; Trivedi, M.; Speth, R.C. Alterations in Gene Expression of Components of the Renin-Angiotensin System and Its Related Enzymes in Lung Cancer. Lung Cancer Int. 2017, 2017, 1–8. [Google Scholar] [CrossRef]
  17. Miao, L.; Chen, W.; Zhou, L.; Wan, H.; Gao, B.; Feng, Y. Impact of Angiotensin I-converting Enzyme Inhibitors and Angiotensin II Type-1 Receptor Blockers on Survival of Patients with NSCLC. Sci. Rep. 2016, 6, 21359. [Google Scholar] [CrossRef]
  18. Aydiner, A.; Ciftci, R.; Sen, F. Renin-Angiotensin System Blockers May Prolong Survival of Metastatic Non-Small Cell Lung Cancer Patients Receiving Erlotinib. Medicine 2015, 94, e887. [Google Scholar] [CrossRef]
  19. Zhang, W.; Liang, Z.; Li, J.; Cai, S. Angiotensin receptor blockers use and the risk of lung cancer: A meta-analysis. J. Renin-Angiotensin-Aldosterone Syst. 2015, 16, 768–773. [Google Scholar] [CrossRef] [Green Version]
  20. Tiret, L.; Rigat, B.; Visvikis, S.; Breda, C.; Corvol, P.; Cambien, F.; Soubrier, F. Evidence, from combined segregation and linkage analysis, that a variant of the angiotensin I-converting enzyme (ACE) gene controls plasma ACE levels. Am. J. Hum. Genet. 1992, 51, 197–205. [Google Scholar]
  21. Chen, J.; Guo, L.; Peiffer, D.A.; Zhou, L.; Chan, O.T.M.; Bibikova, M.; Wickham-Garcia, E.; Lu, S.-H.; Zhan, Q.; Wang-Rodriguez, J.; et al. Genomic profiling of 766 cancer-related genes in archived esophageal normal and carcinoma tissues. Int. J. Cancer 2008, 122, 2249–2254. [Google Scholar] [CrossRef] [PubMed]
  22. Fan, L.; Feng, Y.; Wan, H.; Ni, L.; Qian, Y.; Guo, Y.; Xiang, Y.; Li, Q. Hypoxia induces dysregulation of local renin-angiotensin system in mouse Lewis lung carcinoma cells. Genet. Mol. Res. 2014, 13, 10562–10573. [Google Scholar] [CrossRef] [PubMed]
  23. Christopoulos, P.; Kirchner, M.; Roeper, J.; Saalfeld, F.; Janning, M.; Bozorgmehr, F.; Magios, N.; Kazdal, D.; Volckmar, A.; Brückner, L.; et al. Risk stratification of EGFR+ lung cancer diagnosed with panel-based next-generation sequencing. Lung Cancer 2020, 148, 105–112. [Google Scholar] [CrossRef] [PubMed]
  24. Srivastava, K.; Srivastava, A.; Mittal, B. Angiotensin I-Converting Enzyme Insertion/Deletion Polymorphism and Increased Risk of Gall Bladder Cancer in Women. DNA Cell Biol. 2010, 29, 417–422. [Google Scholar] [CrossRef]
  25. Huang, J.; Chen, S.; Lu, X.; Zhao, Q.; Rao, D.C.; Jaquish, C.E.; Hixson, J.E.; Chen, J.; Wang, L.; Cao, J.; et al. Polymorphisms of ACE2 are Associated with Blood Pressure Response to Cold Pressor Test: The GenSalt Study. Am. J. Hypertens. 2012, 25, 937–942. [Google Scholar] [CrossRef] [Green Version]
  26. Hajjar, I.; Sorond, F.; Hsu, Y.-H.; Galica, A.; Cupples, L.A.; Lipsitz, L.A. Renin angiotensin system gene polymorphisms and cerebral blood flow regulation: The MOBILIZE Boston study. Stroke 2010, 41, 635–640. [Google Scholar] [CrossRef]
  27. Deckers, I.A.G.; Brandt, P.V.D.; Van Engeland, M.; Van Schooten, F.-J.; Godschalk, R.W.; Keszei, A.P.; Schouten, L.J. Polymorphisms in genes of the renin-angiotensin-aldosterone system and renal cell cancer risk: Interplay with hypertension and intakes of sodium, potassium and fluid. Int. J. Cancer 2014, 136, 1104–1116. [Google Scholar] [CrossRef]
  28. Ouyang, J.; Wu, Z.; Xing, J.; Yan, Y.; Zhang, G.; Wang, B.; Li, H.; Ma, X.; Zhang, X. Association of polymorphisms in angiotensin II receptor genes with aldosterone-producing adenoma. J. Huazhong Univ. Sci. Technol. Med. Sci. 2011, 31, 301. [Google Scholar] [CrossRef]
  29. Helisalmi, S.; Hiltunen, M.; Vepsalainen, S.; Iivonen, S.; Mannermaa, A.; Lehtovirta, M.; Koivisto, A.M.; Alafuzoff, I.; Soininen, H. Polymorphisms in neprilysin gene affect the risk of Alzheimer’s disease in Finnish patients. J. Neurol. Neurosurg. Psychiatry 2004, 75, 1746–1748. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, B.; Liu, Q.; Song, Y.; Li, X.; Wang, Y.; Wan, S.; Zhang, Z.; Su, H. Polymorphisms of HIF1A gene are associated with prognosis of early stage non-small-cell lung cancer patients after surgery. Med. Oncol. 2014, 31, 877. [Google Scholar] [CrossRef] [Green Version]
  31. Masago, K.; Fujita, S.; Kim, Y.H.; Hatachi, Y.; Fukuhara, A.; Nagai, H.; Irisa, K.; Ichikawa, M.; Mio, T.; Mishima, M. Effect of vascular endothelial growth factor polymorphisms on survival in advanced-stage non-small-cell lung cancer. Cancer Sci. 2009, 100, 1917–1922. [Google Scholar] [CrossRef] [PubMed]
  32. Sorice, R.; Ruggiero, D.; Nutile, T.; Aversano, M.; Husemoen, L.; Linneberg, A.; Bourgain, C.; Leutenegger, A.-L.; Ciullo, M. Genetic and Environmental Factors Influencing the Placental Growth Factor (PGF) Variation in Two Populations. PLoS ONE 2012, 7, e42537. [Google Scholar] [CrossRef] [PubMed]
  33. Zhou, D.; Chen, Y.; Wu, J.; Shen, J.; Shang, Y.; Zheng, L.; Xie, X. Association between chymase gene polymorphisms and atrial fibrillation in Chinese Han population. BMC Cardiovasc. Disord. 2019, 19, 1–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Kim, D.H.; Xu, W.; Kamel-Reid, S.; Liu, X.; Jung, C.W.; Kim, S.; Lipton, J.H. Clinical relevance of vascular endothelial growth factor (VEGFA) and VEGF receptor (VEGFR2) gene polymorphism on the treatment outcome following imatinib therapy. Ann. Oncol. 2010, 21, 1179–1188. [Google Scholar] [CrossRef]
  35. Luo, Y.; Tanabe, E.; Kitayoshi, M.; Nishiguchi, Y.; Fujiwara, R.; Matsushima, S.; Sasaki, T.; Sasahira, T.; Chihara, Y.; Nakae, D.; et al. Expression of MAS 1 in breast cancer. Cancer Sci. 2015, 106, 1240–1248. [Google Scholar] [CrossRef]
  36. Paradowska-Gorycka, A.; Sowińska, A.; Pawlik, A.; Malinowski, D.; Stypińska, B.; Haladyj, E.; Romanowska-Próchnicka, K.; Olesinska, M. FLT-1 gene polymorphisms and protein expression profile in rheumatoid arthritis. PLoS ONE 2017, 12, e0172018. [Google Scholar] [CrossRef]
  37. Ensembl. Available online: https://www.ensembl.org/index.html (accessed on 3 July 2017).
  38. Brookes, A.J. The essence of SNPs. Gene 1999, 234, 177–186. [Google Scholar] [CrossRef]
  39. Putnam, K.; Shoemaker, R.; Yiannikouris, F.; Cassis, L.A. The renin-angiotensin system: A target of and contributor to dyslipidemias, altered glucose homeostasis, and hypertension of the metabolic syndrome. Am. J. Physiol. Circ. Physiol. 2012, 302, H1219–H1230. [Google Scholar] [CrossRef] [Green Version]
  40. Orte, C.; Polak, J.M.; Haworth, S.G.; Yacoub, M.H.; Morrell, N.W. Expression of pulmonary vascular angiotensin-converting enzyme in primary and secondary plexiform pulmonary hypertension. J. Pathol. 2000, 192, 379–384. [Google Scholar] [CrossRef]
  41. Molina, M.; Xaubet, A.; Li, X.; Abdul-Hafez, A.; Friderici, K.; Jernigan, K.; Fu, W.; Ding, Q.; Pereda, J.; Serrano-Mollar, A.; et al. Angiotensinogen gene G-6A polymorphism influences idiopathic pulmonary fibrosis disease progression. Eur. Respir. J. 2008, 32, 1004–1008. [Google Scholar] [CrossRef] [Green Version]
  42. Tzouvelekis, A.; Gomatou, G.; Bouros, E.; Trigidou, R.; Tzilas, V.; Bouros, D. Common Pathogenic Mechanisms Between Idiopathic Pulmonary Fibrosis and Lung Cancer. Chest 2019, 156, 383–391. [Google Scholar] [CrossRef] [PubMed]
  43. Medeiros, R.; Vasconcelos, A.; Costa, S.; Pinto, D.; Lobo, F.; Morais, A.; Oliveira, J.; Lopes, C. Linkage of angiotensin I-converting enzyme gene insertion/deletion polymorphism to the progression of human prostate cancer. J. Pathol. 2004, 202, 330–335. [Google Scholar] [CrossRef] [PubMed]
  44. Singh, A.; Srivastava, N.; Amit, S.; Prasad, S.; Misra, M.; Ateeq, B. Association of AGTR1 (A1166C) and ACE (I/D) Polymorphisms with Breast Cancer Risk in North Indian Population. Transl. Oncol. 2018, 11, 233–242. [Google Scholar] [CrossRef] [PubMed]
  45. Dai, S.; Ding, M.; Liang, N.; Li, Z.; Li, D.; Guan, L.; Liu, H. Associations of ACE I/D polymorphism with the levels of ACE, kallikrein, angiotensin II and interleukin-6 in STEMI patients. Sci. Rep. 2019, 9, 1–8. [Google Scholar] [CrossRef] [Green Version]
  46. Wan, H.; Feng, Y.; Ni, L.; Fan, L.; Fei, X.; Ma, Q.; Gao, B.; Xiang, Y.; Che, J.; Li, Q.-Y. Overexpression of ACE2 produces antitumor effects via inhibition of angiogenesis and tumor cell invasion in vivo and in vitro. Oncol. Rep. 2011, 26, 1157–1164. [Google Scholar] [CrossRef]
  47. Lopera Maya, E.A.; van der Graaf, A.; Lanting, P.; van der Geest, M.; Fu, J.; Swertz, M.; Franke, L.; Wijmenga, C.; Deelen, P.; Zhernakova, A.; et al. Lack of Association Between Genetic Variants at ACE2 and TMPRSS2 Genes Involved in SARS-CoV-2 Infection and Human Quantitative Phenotypes. Front Genet. 2020, 11, 613. [Google Scholar] [CrossRef]
  48. Santos, R.A.S.; Sampaio, W.O.; Alzamora, A.C.; Motta-Santos, D.; Alenina, N.; Bader, M.; Campagnole-Santos, M.J. The ACE2/Angiotensin-(1-7)/MAS Axis of the Renin-Angiotensin System: Focus on Angiotensin-(1-7). Physiol. Rev. 2018, 98, 505–553. [Google Scholar] [CrossRef] [Green Version]
  49. Pavo, N.; Arfsten, H.; Cho, A.; Goliasch, G.; Bartko, P.E.; Wurm, R.; Freitag, C.; Gisslinger, H.; Kornek, G.; Strunk, G.; et al. The circulating form of neprilysin is not a general biomarker for overall survival in treatment-naïve cancer patients. Sci. Rep. 2019, 9, 1–8. [Google Scholar] [CrossRef] [Green Version]
  50. Kadota, K.; Buitrago, D.H.; Lee, M.-C.; Villena-Vargas, J.; Sima, C.S.; Jones, D.R.; Travis, W.D.; Adusumilli, P.S. Tumoral CD10 expression correlates with high-grade histology and increases risk of recurrence in patients with stage I lung adenocarcinoma. Lung Cancer 2015, 89, 329–336. [Google Scholar] [CrossRef] [Green Version]
  51. Comings, D.E.; MacMurray, J.P. Molecular Heterosis: A Review. Mol. Genet. Metab. 2000, 71, 19–31. [Google Scholar] [CrossRef]
  52. Xu, K.; Kosoy, R.; Shameer, K.; Kumar, S.; Liu, L.; Readhead, B.; Belbin, G.M.; Lee, H.-C.; Chen, R.; Dudley, J.T. Genome-wide analysis indicates association between heterozygote advantage and healthy aging in humans. BMC Genet. 2019, 20, 1–14. [Google Scholar] [CrossRef]
  53. Lippman, Z.B.; Zamir, D. Heterosis: Revisiting the magic. Trends Genet. 2007, 23, 60–66. [Google Scholar] [CrossRef] [PubMed]
  54. Manikandan, M.; Munirajan, A.K. Single Nucleotide Polymorphisms in MicroRNA Binding Sites of Oncogenes: Implications in Cancer and Pharmacogenomics. OMICS J. Integr. Biol. 2014, 18, 142–154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Liu, C.X.; Tan, L.; Sun, F.R.; Zhang, W.; Miao, D.; Tan, M.S.; Wan, Y.; Tan, C.-C.; Yu, J.T.; Tan, L. The Association of MME microRNA Binding Site Polymorphism with the Risk of Late Onset Alzheimer’s Disease in Northern Han Chinese. Curr. Neurovasc. Res. 2017, 14, 90–95. [Google Scholar] [CrossRef]
  56. Barenboim, M.; Zoltick, B.J.; Guo, Y.; Weinberger, D.R. MicroSNiPer: A web tool for prediction of SNP effects on putative microRNA targets. Hum. Mutat. 2010, 31, 1223–1232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. MicroSNiPer. Available online: http://vm24141.virt.gwdg.de/services/microsniper/ (accessed on 20 October 2020).
  58. Krüger, J.; Rehmsmeier, M. RNAhybrid: MicroRNA target prediction easy, fast and flexible. Nucleic Acids Res. 2006, 34, W451–W454. [Google Scholar] [CrossRef]
  59. BiBiServ. Available online: http://bibiserv.techfak.uni-bielefeld.de/rnahybrid (accessed on 20 October 2020).
  60. Ludwig, N.; Leidinger, P.; Becker, K.; Backes, C.; Fehlmann, T.; Pallasch, C.; Rheinheimer, S.; Meder, B.; Stähler, C.; Meese, E.; et al. Distribution of miRNA expression across human tissues. Nucleic Acids Res. 2016, 44, 3865–3877. [Google Scholar] [CrossRef]
  61. Miners, S.; Van Helmond, Z.; Barker, R.; A Passmore, P.; A Johnston, J.; Todd, S.; McGuinness, B.M.; Panza, F.; Seripa, D.; Solfrizzi, V.; et al. Genetic variation in MME in relation to neprilysin protein and enzyme activity, Aβ levels, and Alzheimer’s disease risk. Int. J. Mol. Epidemiol. Genet. 2012, 3, 30–38. [Google Scholar]
  62. Leithner, A.; Wohlkoenig, C.; Stacher, E.; Lindenmann, J.; Hofmann, N.A.; Gallé, B.; Guelly, C.; Quehenberger, F.; Stiegler, P.; Smolle-Jüttner, F.M.; et al. Hypoxia increases membrane metallo-endopeptidase expression in a novel lung cancer ex vivo model—Role of tumor stroma cells. BMC Cancer 2014, 14, 40. [Google Scholar] [CrossRef] [Green Version]
  63. Muz, B.; De La Puente, P.; Azab, F.; Azab, A.K. The role of hypoxia in cancer progression, angiogenesis, metastasis, and resistance to therapy. Hypoxia 2015, 3, 83–92. [Google Scholar] [CrossRef] [Green Version]
  64. Guo, S.; Martin, M.G.; Tian, C.; Cui, J.; Wang, L.; Wu, S.; Gu, W. Evaluation of Detection Methods and Values of Circulating Vascular Endothelial Growth Factor in Lung Cancer. J. Cancer 2018, 9, 1287–1300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. George, A.J.; Thomas, W.G.; Hannan, R.D. The renin–angiotensin system and cancer: Old dog, new tricks. Nat. Rev. Cancer 2010, 10, 745–759. [Google Scholar] [CrossRef] [PubMed]
  66. Koukourakis, M.I.; Papazoglou, D.; Giatromanolaki, A.; Bougioukas, G.; Maltezos, E.; Sivridis, E. VEGF gene sequence variation defines VEGF gene expression status and angiogenic activity in non-small cell lung cancer. Lung Cancer 2004, 46, 293–298. [Google Scholar] [CrossRef] [PubMed]
  67. Scartozzi, M.; Faloppi, L.; Svegliati-Baroni, G.; Loretelli, C.; Piscaglia, F.; Iavarone, M.; Toniutto, P.; Fava, G.; De Minicis, S.; Mandolesi, A.; et al. VEGF and VEGFR genotyping in the prediction of clinical outcome for HCC patients receiving sorafenib: The ALICE-1 study. Int. J. Cancer 2014, 135, 1247–1256. [Google Scholar] [CrossRef] [Green Version]
  68. Zhai, R.; Gong, M.N.; Zhou, W.; Thompson, T.B.; Kraft, P.; Su, L.; Christiani, D.C. Genotypes and haplotypes of the VEGF gene are associated with higher mortality and lower VEGF plasma levels in patients with ARDS. Thorax 2007, 62, 718–722. [Google Scholar] [CrossRef] [Green Version]
  69. Krippl, P.; Langsenlehner, U.; Renner, W.; Yazdani-Biuki, B.; Wolf, G.; Wascher, T.C.; Paulweber, B.; Haas, J.; Samonigg, H. A common 936 C/T gene polymorphism of vascular endothelial growth factor is associated with decreased breast cancer risk. Int. J. Cancer 2003, 106, 468–471. [Google Scholar] [CrossRef]
  70. Petrovic, M.G.; Korosec, P.; Kosnik, M.; Osredkar, J.; Hawlina, M.; Peterlin, B.; Petrovič, D. Local and genetic determinants of vascular endothelial growth factor expression in advanced proliferative diabetic retinopathy. Mol. Vis. 2008, 14, 1382–1387. [Google Scholar]
  71. Glubb, D.M.; Cerri, E.; Giese, A.; Zhang, W.; Mirza, O.; Thompson, E.E.; Chen, P.; Das, S.; Jassem, J.; Rzyman, W.; et al. Novel Functional Germline Variants in the VEGF Receptor 2 Gene and Their Effect on Gene Expression and Microvessel Density in Lung Cancer. Clin. Cancer Res. 2011, 17, 5257–5267. [Google Scholar] [CrossRef] [Green Version]
  72. Keshavarz, L.; Yavarian, M. The association of Q472H variant in the KDR gene with recurrent pregnancy loss in Southern Iran: A case-control study. Int. J. Reprod. Biomed. 2019, 17, 473–480. [Google Scholar] [CrossRef]
  73. Meert, A.-P.; Paesmans, M.; Martin, B.; Delmotte, P.; Berghmans, T.; Verdebout, J.-M.; Lafitte, J.-J.; Mascaux, C.; Sculier, J.-P. The role of microvessel density on the survival of patients with lung cancer: A systematic review of the literature with meta-analysis. Br. J. Cancer 2002, 87, 694–701. [Google Scholar] [CrossRef] [Green Version]
  74. Trivella, M.; Pezzella, F.; Pastorino, U.; Harris, A.L.; Altman, D.G.; Prognosis In Lung Cancer Collaborative Study Group. Microvessel density as a prognostic factor in non-small-cell lung carcinoma: A meta-analysis of individual patient data. Lancet Oncol. 2007, 8, 488–499. [Google Scholar] [CrossRef]
Table
Table 1. Patient′s clinical characteristics.
Table 1. Patient′s clinical characteristics.
Clinical Variables.Discovery Set (N = 167)Validation Set (N = 190)p
Gender, n
Males123 (73.7%)133 (70.0%)
Females44 (26.3%)57 (30.0%)0.445 *
Age (years)
Median (IQR)64.0 (58.0–72.0)65.0 (58.0–73.0)0.773 **
Smoking Status, n (%)
Never smokers31 (27.7%)20 (10.5%)
Previously smokers68 (60.7%)97 (51.1%)
Smokers13 (11.6%)73 (38.4%)<0.0001 *
Hypertension, n (%)
Yes67 (48.6%)91 (47.9%)
No71 (51.4%)99 (52.1%)0.911 *
Anti-hypertension drugs, n (%)
No86 (65.2%)126 (67.0%)
iACE/ARB46 (34.8%)62 (33.0%)0.728 *
ECOG PS, n (%)
0–2163 (97.6%)188 (98.9%)
>24 (2.4%)2 (1.1%)0.153 *
Histology, n (%)
Adenocarcinoma116 (69.5%)126 (66.3%)
Squamous cell carcinoma42 (25.1%)62 (32.6%)
Adeno + Squamous cell9 (5.4%)2 (1.1%)0.066 *
cTNM Stage, n (%)
IIIA20 (12.0%)27 (14.2%)
IIIB 32 (19.2%)36 (18.9%)
IIIC21 (12.6%)18 (9.5%)
IVA 65 (38.9%)66 (34.7%)
IVB29 (17.4%)43 (22.6%)0.597 *
Treatment modality, n (%)
Surgery + ChT12 (7.2%)0 (0.0%)
ChT + RT31 (18.6%)81 (42.6%)
ChT124 (74.3%) 109 (57.4%)<0.0001 *
First-line treatment, n (%)
Platinum based ChT145 (86.8%)155 (81.6%)
TKI22 (12.8%)35 (18.4%)0.194 *
PFS (months)
Median, 95% CI7.5 (6.1–9.0)8.1 (6.8–9.4)0.415 ***
OS (months)
median, 95% CI30.0 (16.9–43.2)20.6 (16.3–24.9)0.036 ***
IQR—interquartil range; iACE—inhibitors of angiotensin-converting enzyme, ARB—AngiotensinII receptor blocker, ECOG PS—Eastern cooperative Oncology Group Performance Status; ChT—chemotherapy; RT—radiotherapy; TKI—tyrosine kinase inhibitors; 95% CI—confidence interval. PFS—progression-free survival, OS—overall survival. * Pearson Chi Square test. ** Mann-Whitney U test. *** Kaplan Meier curve with log rank test.
Table
Table 2. Multivariate analysis for progression-free survival (PFS) in the discovery cohort (n = 167).
Table 2. Multivariate analysis for progression-free survival (PFS) in the discovery cohort (n = 167).
VariablesPFS
(Multivariate Stepwise Cox Proportional Hazards Model)		Bootstrap
HR (95% CI)pHR (95% CI)p
T
 T1Referent Referent
 T2
 T3-41.7 (1.3–2.2)<0.0001 *1.7 (1.3–2.3)<0.0001 *
Distant metastasis
 NoReferent Referent
 Yes1.7 (1.2–2.4)0.0051.7 (1.1–2.5)0.009
Type of Therapy
 Surgery + ChT Referent Referent
 ChT
 ChT + RT1.5 (1.1–2.1)0.012 *1.5 (1.0–2.3)0.053 *
Type of Systemic Therapy, 1st line
 Platinum based ChTReferent Referent
 TKi0.4 (0.2–0.8)0.0030.4 (0.2–0.7)0.002
KDR rs1870377, Heterosis Model
 TAReferent Referent
 TT/AA1.7 (1.2–2.5)0.0051.8 (1.2–7)0.003
MME rs701109, Dominant Model
 TTReferent Referent
 TC/CC0.5 (0.3–0.9)0.0180.5 (0.2–0.9)0.030
Analyses were conducted using stepwise Cox regression, followed by Bootstrap analysis (Montecarlo simulations with 1000 iterations). HR, hazard ratio, 95% CI, confident interval, ChT, chemotherapy, RT, radiotherapy, TKi, tyrosine kinase inhibitors. * p trend value.
Table
Table 3. Multivariate analysis of relevant variables for overall survival (OS) in the discovery set (n = 167) using stepwise Cox regression, followed by Bootstrap analysis (Montecarlo simulations with 1000 iterations).
Table 3. Multivariate analysis of relevant variables for overall survival (OS) in the discovery set (n = 167) using stepwise Cox regression, followed by Bootstrap analysis (Montecarlo simulations with 1000 iterations).
VariablesOS
(Multivariate Stepwise Cox Proportional Hazards Model)		Bootstrap
HR (95% CI)pHR (95% CI)p
T
 T1Referent Referent
 T2
 T3-41.5 (1.1–2.2)0.018 *1.5 (1.0–2.2)0.066 *
Distant metastasis
 NoReferent Referent
 Yes2.5 (1.5–4.3)0.0012.6 (1.5–4.5)0.001
Age median
 ≤64.0 Referent Referent
 >64.01.7 (1.0–2.7)0.0311.6 (1.0–2.8)0.069
KDR rs1870377, Heterosis Model
 TAReferent Referent
 TT/AA2.1 (1.2–3.4)0.0062.0 (1.2–3.3)0.012
HR, hazard ratio, 95% CI, confident interval. * p trend value.
Table
Table 4. Multivariate analysis of relevant variables for disease progression using stepwise Cox regression, followed by Bootstrap analysis with 1000 iterations from validation set (n = 190).
Table 4. Multivariate analysis of relevant variables for disease progression using stepwise Cox regression, followed by Bootstrap analysis with 1000 iterations from validation set (n = 190).
VariablesPFS
(Multivariate Stepwise Cox Proportional Hazards Model)		Bootstrap
HR (95%CI)pHR (95%CI)p
Histology
 AdenocarcinomaReferent Referent
 Squamous cell1.7 (1.2–2.4)0.003 **1.7 (1.2–2.5)0.003 **
 Other *- -
Type of Therapy
 ChTReferent Referent
 ChT + RT0.4 (0.3–0.6)<0.00010.4 (0.3–0.6)<0.0001
Type of Systemic Therapy, 1st line
 Platinum based ChTReferent Referent
 TKi0.4 (0.3–0.7)0.0010.4 (0.3–0.7)0.002
MME rs701109, Heterosis Model
 TCReferent Referent
 TT/CC1.6 (1.2–2.3)0.0051.6 (1.1–2.4)0.007
PFS—progression-free survival, HR—hazard ratio, 95% CI—confident interval, ChT—chemotherapy, RT—radiotherapy, TKi—inhibitors of tyrosine kinase. * Adenocarcinoma+squamous cell. ** p trend value.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Catarata, M.J.; Medeiros, R.; Oliveira, M.J.; Pêgo, A.; Frade, J.G.; Martins, M.F.; Robalo Cordeiro, C.; Herth, F.J.F.; Thomas, M.; Kriegsmann, M.; et al. Germline Genetic Variants of the Renin-Angiotensin System, Hypoxia and Angiogenesis in Non-Small Cell Lung Cancer Progression: Discovery and Validation Studies. Cancers 2020, 12, 3834. https://doi.org/10.3390/cancers12123834

AMA Style

Catarata MJ, Medeiros R, Oliveira MJ, Pêgo A, Frade JG, Martins MF, Robalo Cordeiro C, Herth FJF, Thomas M, Kriegsmann M, et al. Germline Genetic Variants of the Renin-Angiotensin System, Hypoxia and Angiogenesis in Non-Small Cell Lung Cancer Progression: Discovery and Validation Studies. Cancers. 2020; 12(12):3834. https://doi.org/10.3390/cancers12123834

Chicago/Turabian Style

Catarata, Maria Joana, Rui Medeiros, Maria José Oliveira, Alice Pêgo, João Gonçalo Frade, Maria Fátima Martins, Carlos Robalo Cordeiro, Felix J F Herth, Michael Thomas, Mark Kriegsmann, and et al. 2020. "Germline Genetic Variants of the Renin-Angiotensin System, Hypoxia and Angiogenesis in Non-Small Cell Lung Cancer Progression: Discovery and Validation Studies" Cancers 12, no. 12: 3834. https://doi.org/10.3390/cancers12123834

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop