Next Article in Journal
NPM-ALK: A Driver of Lymphoma Pathogenesis and a Therapeutic Target
Next Article in Special Issue
Systematic Review and Metanalysis of Oncomarkers in IPF Patients and Serial Changes of Oncomarkers in a Prospective Italian Real-Life Case Series
Previous Article in Journal
Global Burden, Risk Factors, and Trends of Esophageal Cancer: An Analysis of Cancer Registries from 48 Countries
Previous Article in Special Issue
RNA-Based Assay for Next-Generation Sequencing of Clinically Relevant Gene Fusions in Non-Small Cell Lung Cancer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 1
10.3390/cancers13010143
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:
Systematic Review

Functional Signatures in Non-Small-Cell Lung Cancer: A Systematic Review and Meta-Analysis of Sex-Based Differences in Transcriptomic Studies

by
Irene Pérez-Díez
1,2,
Marta R. Hidalgo
1,
Pablo Malmierca-Merlo
1,3,
Zoraida Andreu
1,
Sergio Romera-Giner
1,3,
Rosa Farràs
4,
María de la Iglesia-Vayá
2,
Mariano Provencio
5,
Atocha Romero
5 and
Francisco García-García
1,6,*
1
Bioinformatics and Biostatistics Unit, Principe Felipe Research Center (CIPF), 46012 Valencia, Spain
2
Biomedical Imaging Unit FISABIO-CIPF, Fundación para el Fomento de la Investigación Sanitaria y Biomédica de la Comunidad Valenciana, 46012 Valencia, Spain
3
Atos Research Innovation (ARI), 28037 Madrid, Spain
4
Department of Oncogenic Signalling, Principe Felipe Research Center (CIPF), 46012 Valencia, Spain
5
Medical Oncology Department, Hospital Universitario Puerta de Hierro-Majadahonda, 28222 Madrid, Spain
6
Spanish National Bioinformatics Institute, ELIXIR-Spain (INB, ELIXIR-ES), C/Eduardo Primo Yúfera, 3, 46012 Valencia, Spain
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(1), 143; https://doi.org/10.3390/cancers13010143
Submission received: 16 November 2020 / Revised: 29 December 2020 / Accepted: 31 December 2020 / Published: 5 January 2021
(This article belongs to the Special Issue Diagnostic and Predictive Biomarkers in Lung Cancer)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Differences in epidemiological and clinical patterns of lung adenocarcinoma have been described between male and female patients. Our study aimed to assess the molecular mechanisms underlying those differences through functional profiling and meta-analysis of lung adenocarcinoma expression datasets. We discovered an overrepresentation of terms related to acute inflammatory responses in female patients, suggesting an increase in certain cell populations, such as T CD8+ cells, and in interleukin production. We also identified purinergic signaling and lipid metabolism as relevant upregulated functional groups in female lung adenocarcinoma samples. Overall, we identified molecular processes and pathways that will aid the understanding of the differential clinical patterns described in lung adenocarcinoma in male and female patients and provide new pathways for the discovery of novel biomarkers and therapeutic targets.

Abstract

While studies have established the existence of differences in the epidemiological and clinical patterns of lung adenocarcinoma between male and female patients, we know relatively little regarding the molecular mechanisms underlying such sex-based differences. In this study, we explore said differences through a meta-analysis of transcriptomic data. We performed a meta-analysis of the functional profiling of nine public datasets that included 1366 samples from Gene Expression Omnibus and The Cancer Genome Atlas databases. Meta-analysis results from data merged, normalized, and corrected for batch effect show an enrichment for Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways related to the immune response, nucleic acid metabolism, and purinergic signaling. We discovered the overrepresentation of terms associated with the immune response, particularly with the acute inflammatory response, and purinergic signaling in female lung adenocarcinoma patients, which could influence reported clinical differences. Further evaluations of the identified differential biological processes and pathways could lead to the discovery of new biomarkers and therapeutic targets. Our findings also emphasize the relevance of sex-specific analyses in biomedicine, which represents a crucial aspect influencing biological variability in disease.

Graphical Abstract

1. Introduction

Lung cancer is the most frequently diagnosed cancer and the leading cause of cancer-related death worldwide, representing 18.4% of all cancer deaths [1]. Exposure to tobacco, domestic biomass fuels, asbestos, and radon represent the most relevant lung cancer risk factors [1,2,3]; however, as has become evident from studies of other cancer types, sex-based differences (sexual dimorphisms) may also have significant relevance in lung cancer [1,4,5].
Lung cancer exhibits sex-based disparities in clinical characteristics and outcomes, with better survival observed in women [3,6,7]. While lung cancer incidence worldwide is higher in men, there exists an increasing trend in women that cannot be solely explained by tobacco consumption [1,2]. Furthermore, studies have reported sex-dependent differences in estrogen receptors and their impact on lung cancer [8,9,10]; however, conflicting results have attributed lung cancer susceptibility in women to genetic variants, hormonal factors, molecular abnormalities, and oncogenic viruses [3,11,12,13].
Adenocarcinoma represents the most frequent non-small cell lung cancer (NSCLC) subtype in both sexes [14], with a higher predominance in women compared to men (41% of cases in women versus 34% in men) [3,7,13]. Interestingly, Wheatley-Price et al. demonstrated a more pronounced survival rate difference between male and female lung adenocarcinoma patients when compared to other tumor types [7]. The molecular causes underlying such sex-biased patterns remain largely unknown, as limited efforts have been made for lung adenocarcinoma, with few studies considering this differential perspective [15,16,17,18]. These limitations can be partially addressed through meta-analysis, a robust methodology that combines information from related but independent studies to derive results with increased statistical power and precision [19,20]. As current treatment strategies do not cure most lung cancer patients, and invasive diagnostic techniques (e.g., via biopsies and bronchoscopies) often induce discomfort in patients, meta-analyses that improve our understanding of sex-specific molecular mechanisms in lung adenocarcinoma may facilitate the discovery of non-invasive prognostic and diagnostic biomarkers.
We performed a meta-analysis based on functional profiles of transcriptomic studies to explore the molecular mechanisms underlying sex-based differences in early-stage lung adenocarcinoma. We carried out exhaustive review and selection steps to guarantee the homogeneity of the selected studies and the subsequent comparison and integration of the data in the meta-analysis with an appreciation of this strategy’s specific limitations. After the systematic review, we retrieved and analyzed nine studies from Gene Expression Omnibus (GEO) [21] and The Cancer Genome Atlas (TCGA) [22], and then combined the results in a random-effects meta-analysis. This approach allowed the identification of functional alterations caused by lung adenocarcinoma in both male and female patients, comparing tumor samples with adjacent non-tumor tissue. In this study, we identified immune responses, purinergic signaling, and lipid metabolism as the main biological processes that display differences between male and female lung adenocarcinoma patients, with the acute immune response increased in female patients. Overall, our findings provide evidence that sex-based differences influence cancer biology and may impact response to treatment. Furthermore, underlying sex-based differences may contribute to the discovery of sex-specific prognostic and diagnostic biomarkers and the improvement of personalized therapies.

2. Results

We organized our findings into three sections: the first describes the studies evaluated and selected in the systematic review; the second section reports on the results of the bioinformatic analysis of each of these selected studies as follows: (i) exploratory analysis, (ii) differential expression, and (iii) functional characterization; while the third section presents the results of the differential functional profiling by sex.

2.1. Study Search and Selection

The systematic review identified 207 non-duplicated studies, of which 48.8% included both male and female patients (Figure S1). We applied inclusion and exclusion criteria (see Figure 1) to select a set of homogeneous and comparable studies to ensure the reliability of the subsequent analyses. To ensure study homogeneity (and in the hope of contributing to the early diagnosis of disease), we focused on those studies of early-stage disease and selected nine transcriptomic studies for further analysis (Table 1). The selected studies represented a population of 1366 early-stage samples (369 controls and 997 cases), of which 44% were from men, and 56% from women (Figure 2), with a median age of 65.54 years old. Table 1, Figure 2, and Supplementary Table S1 contain further information regarding the selected studies and the clinicopathological characteristics of the study population.

2.2. Individual Analysis

As the normalized data derives from different platforms, we performed exploratory and processing steps for the data set to ensure the comparability and integration of subsequent analyses. The exploratory analysis found a lack of abnormal behavior except for three samples in the principal component analysis (PCA) and unsupervised clustering; therefore, we excluded the GSM47570 and GSM47578 samples in study GSE19188, and the GMS773784 sample in study GSE31210 from further analysis.
The differential expression results for each study demonstrated a large number of differentially expressed genes when comparing female lung adenocarcinoma samples to female control samples and male adenocarcinoma samples to male control samples (Supplementary Table S2). However, the evaluation of sex-based differences in lung adenocarcinoma patients provided a small number of significantly affected genes (see Supplementary Table S3), with no intersecting genes.
We performed an individual functional enrichment analysis of Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) to identify the possible implications of these sex-specific differentially-expressed genes in pathways relevant to lung adenocarcinoma. The identified pathways revealed a diversity of significant results among datasets, which we have summarized in Table 2. When analyzing intersections, UpSet plots (Figure 3, equivalent to a Venn diagram) illustrate the degree of intersection between studies, demonstrating that most significant results are exclusive of each study. This data highlights the need for integrated strategies, such as meta-analyses, to increase the statistical power of any findings.

2.3. Meta-Analysis

We performed a functional meta-analysis for each of the 8672 GO functions and KEGG pathways, including every term found in at least two studies. Results with a false discovery rate (FDR) of <0.05 included 106 GO biological processes (BP), 3 GO molecular functions (MF), and 20 KEGG pathways, which were associated with 21 wider functional groups. We rejected potential bias on the significant results after the inspection of funnel plots; furthermore, sensitivity analyses failed to indicate alterations in the results due to the inclusion of any study. The results for the 129 significant GO terms and KEGG pathways are further detailed in Supplementary Table S4, including FDR, the log odds ratio (LOR), and its 95% confidence interval (CI), and the standard error (SE) of the LOR. In addition, Supplementary Table S5 details the number of functions involved in each of the genes of this functional signature, and Supplementary Figure S2 displays the prognostic description of this set of significant functions.

2.3.1. Upregulated Functions

We discovered that 43.88% of detected functions related to the immune response (Supplementary Table S4 and Figure 4), which all displayed upregulation in female lung adenocarcinoma patients. This finding agrees with studies that report more robust innate and adaptive immune responses in women who, for example, present with more efficient antigen-presenting cells (APCs) than males [33,34]. The results provided evidence for the positive regulation of an acute inflammatory response, with CD8+ alpha-beta T cell differentiation and activation, B cell proliferation and activation, and an increase of interleukin (IL) biosynthesis, including IL-2, 6, 8, 10, and 17. Several immune-related signaling pathways also displayed differences between female and male lung adenocarcinoma patients. We discovered the upregulation of the MyD88-independent toll-like receptor signaling pathway, NIK/NF-kappa β signaling pathway, FC-epsilon receptor signaling pathway, B-cell receptor signaling pathway, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, and RIG-I-like receptor signaling pathways in female lung adenocarcinoma patients. Overall, these findings suggest the relevance of sex-based differences in immune responses to lung cancer, which may represent a significant contributor to the marked differences observed during disease progression in male and female patients. Furthermore, such findings may help to define novel therapeutic targets and may have important implications for immunotherapy.
We also uncovered evident sex-based differences in cell metabolism. Studies have shown the significant upregulation of lipid metabolism (ceramide and sphingolipid biosynthetic processes) in females and the higher utilization of carbohydrates by males [35,36]. Furthermore, lipid metabolism and signaling are widely accepted as major players in cancer biology [37]. “Metabolism- Nucleic acids metabolism and signaling” was the second most abundant functional group upregulated in female lung adenocarcinoma patients, comprising 8.63% of the altered functions. These GO terms and KEGG pathways are mainly related to purinergic signaling through G protein-coupled receptors and cytidine metabolism. Other functional groups upregulated in female lung adenocarcinoma patients include cell migration and homeostasis. Overall, further explorations of sex differences in metabolic pathways may provide new perspectives for treatment approaches with sex-specific effects.

2.3.2. Downregulated Functions

It should be noted that 23.26% of the significant functions exhibited lower activity in female compared to male lung adenocarcinoma patients. Downregulated functional groups include those related to cell cycle progression, cell junctions, DNA repair and telomere protection, mitochondrial processes, neural development, post-translational changes, post-transcriptional changes, protein degradation, and transcription regulation (Supplementary Table S4 and Figure 4). Taken together, the downregulated pathways in female patients suggest the existence of lower levels of oxidative stress compared to male patients [38], which may contribute to the existence of a less permissive tumorigenic environment.

2.4. Metafun-NSCLC Web Tool

The Metafun-NSCLC web tool (https://bioinfo.cipf.es/metafun-nsclc) contains information related to the nine studies and 1329 samples involved in this study. For each study, this resource includes fold-changes of genes and LOR of functions and pathways that users can explore to identify profiles of interest.
We carried out a total of 8672 meta-analyses. For each of the 129 significant functions and pathways, Metafun-NSCLC depicts the global activation level for all studies and each study’s specific contribution using statistical indicators (LOR, CI, and p-value) and graphical representations by function as forest and funnel plots. This open resource hopes to contribute to data sharing between researchers, the elaboration of innovative studies, and the discovery of new findings.
Here, we also highlight the importance of including/reporting sex-related data in the results of clinical studies given their general importance in tumor risk, treatment response, and outcomes in lung adenocarcinoma and other cancers/disorders. The integration of sex-based differences in this manner has the potential to significantly impact the cancer biology field.

3. Discussion

Despite the profuse evidence for the influence of sex on rates and patterns of metastasis, the expression of prognostic biomarkers, and therapeutic responses in several cancer types [39,40], sex-based differences have not been consistently considered when studying cancer, designing therapies, or constructing clinical trials. Cases of NSCLC, including adenocarcinoma, exhibit differences in incidence, prevalence, and severity in female and male patients [1,3,7,41]. Elucidating the molecular basis for this sex-based differential impact will have clinical relevance, as this information can guide/improve both diagnosis and treatment.
Biomedical research generally underrepresents female patients, with sex-based differences rarely considered [42,43]. Our systematic review of transcriptomic studies revealed that only 48.8% of lung adenocarcinoma-related datasets considered both sexes, a figure similar (49%) to that reported by Woitowich et al. [43]. Sex-based differences impact disease biomarkers, drug responses, and treatment [42], and, therefore, sex must represent a critical component of experimental design. Added to this problem, we faced a lack of standardization among studies and detailed clinical information (i.e., mutations, smoking status, stages) when searching for suitable datasets. The consideration of Findable, Accessible, Interoperable, and Reusable (FAIR) data principles [44], a requisite for quality science, would ensure that generated data can be of further use throughout the scientific community.
To the best of our knowledge, only four studies have attempted to address the functional alterations caused by lung adenocarcinoma in both male and female patients—those by Shi et al. [16], which considered female patients, Araujo et al. [15], Yuan et al. [17], and Li et al. [18]. Shi et al. [16] integrated samples from two datasets for a differential expression analysis followed by a functional enrichment analysis, whereas Araujo et al. [15] independently processed six datasets and jointly analyzed their results. Yuan et al. [17] compared male and female patients with various cancer types (including lung adenocarcinoma) using the TCGA dataset, but did not include control samples in the statistical comparison. Li et al. [18] evaluated the differences in lung adenocarcinoma by focusing only on metabolic pathways. By including both male and female patients and controls in our gene expression comparison, in contrast with the analysis performed by Shi et al. [16] and Yuan et al. [17], we have effectively unveiled sex-based differences occurring in lung adenocarcinoma. Of note, 88% of the results reported by Yuan et al. [17] relate to the sex chromosome, but not necessarily due to cancer. Shi et al. [16] described the consequences of cancer development in female patients, but the authors failed to compare said effects with male patients. Although Araujo et al. [15] do not describe the statistical comparisons performed (and do not perform a statistical integration such as meta-analysis), the authors describe the results obtained for each dataset. Our study addressed sex-based differences in male and female lung adenocarcinoma patients through meta-analysis to address previous limitations and improve on those approaches employed by others. Despite certain supposed limitations to our approach (the presence of studies with different sample sizes and types of platforms), meta-analyses can integrate selected studies by eliminating inconsistency in individual studies, thereby increasing the statistical power, and highlighting robust disease-associated functions. We performed a meta-analysis using a random-effects model on Gene Set Enrichment Analysis (GSEA) results independently obtained from each study to evaluate the functions differentially altered between male and female lung adenocarcinoma patients. To make results comparable and reduce biases in the type of analysis used, we applied the same bioinformatics strategy from normalized expression matrices to GSEA results. This robust methodology integrates groups of data and provides results with higher statistical power and precision [19,20] and reveals findings that cannot be obtained through the intersection or addition of results in individual studies. Nevertheless, we selected only samples from early-stage lung adenocarcinoma patients to reduce variability and included smoking status in the differential expression linear model. While the inclusion of the mutational status of relevant genes to lung adenocarcinoma (e.g., EGFR) could have revealed important insight, the lack of this information in the majority of the studies and the resultant limited sample size hampered this aim. We would support the inclusion of this type of data as a requirement for the publication of new data to repositories, such as the GEO.
The immune system plays a crucial role in the development of cancer [45], and several studies have reported sex-based differences in immune responses (reviewed by Klein and Flanagan [33]). Tumor cells evade the immune system using different strategies, including the modulation of antigen-presentation and the suppression of regulatory T cells. Therefore, sex differences in APCs and their downstream effector cells, among other components, may contribute to the sexual disparity observed in various aspects of cancer development and may significantly impact antitumor immunity and immunotherapy. Adult females generally present more robust innate and adaptive immune responses than males, as evidenced by increased phagocytic activity of neutrophils and macrophages, more efficient APCs, and differences in lymphocyte subsets (B cells, CD4+T cells, CD8+T cells) and cytokine production. Accordingly, our results demonstrate an enrichment of immune response-related terms in female lung adenocarcinoma patients, which agrees with the findings of Araujo et al. [15]. The analyzed functions suggest the positive regulation of CD8+ alpha-beta T cell activation and differentiation in female lung adenocarcinoma patients, which play an essential role in antitumor immunity [46,47]. Furthermore, Ye et al. [47] discovered a more abundant population of effector memory CD8+ T cells in female lung adenocarcinoma patients, which agrees with our results. A previous study described CD8+ lymphocyte levels as a prognostic biomarker in NSCLC [48], and specifically in lung adenocarcinoma [49], with a correlation between higher levels of CD8+ lymphocytes with higher survival rates and lower disease recurrence. Elevated levels of active CD8+ T cells in female lung adenocarcinoma patients could form part of the molecular mechanism underlying higher survival rates when compared to male lung adenocarcinoma patients. Activation of the Notch signaling pathway decreases CD8+ T lymphocyte activity in lung adenocarcinoma [50]; therefore, the downregulation of the Notch signaling pathway discovered in female lung adenocarcinoma patients could explain higher CD8+ T activity when compared to male lung adenocarcinoma patients.
Concerning the immune response, we also detected differences that supported the increased production of IL-2, which is known to stimulate T cell proliferation and the production of effector T cells, thereby amplifying the lymphocytic response [51]. Higher levels of IL-2 could also relate to increased activity of CD8+ T cells in female lung adenocarcinoma patients. Increased levels of IL-10 are also supported in female lung adenocarcinoma patients and, although IL-10 has anti-inflammatory and anti-immune activities [52,53], studies have suggested a dual role in cancer. In advanced lung adenocarcinoma, high expression of IL-10 receptor 1 correlates with worse prognosis [52], while IL-10 expression by T-regulatory cells inhibits apoptosis through Programmed death-ligand 1 inhibition [53]. Nevertheless, IL-10 correlates with better prognosis when expressed by CD8+ T cells in early-stage NSCLC [54], and it seems to activate the antitumor control of CD8+ T cells [55]. IL-2 and IL-10 could display increased activity in early-stage female patients, alongside a higher population of active CD8+ T cells than males, conferring women a survival advantage.
We also detected the positive regulation of IL-6 biosynthesis in female lung adenocarcinoma patients, with increased IL-6 levels correlating with worse prognosis in NSCLC patients in previous studies [56,57]. Network analysis in non-smoking female lung adenocarcinoma patients described IL-6 as one of the pathology’s central nodes [16], and these findings agree with our results, which provide evidence of the critical role of IL-6 in tumor progression in female lung adenocarcinoma patients. IL-8 and IL-17 exhibit increased production and biosynthesis in female lung adenocarcinoma patients, with said interleukins known to influence tumor growth and metastasis and correlate with worse prognosis [58,59,60].
Although altered immune responses can positively and negatively influence tumor progression, our findings have detected GO terms that point to an elevated acute immune response in female compared to male lung adenocarcinoma patients. Of note, sex-based immunological differences in lung adenocarcinoma might have an impact on immunotherapy response. Different studies have addressed the role of sex in immunotherapy [5,47,61,62] and established improved survival for female NSCLC patients. The discovered molecular pathways differentially activated between male and female lung adenocarcinoma patients may underlie phenotypic differences regarding immunotherapy response.
Sex-based differences in metabolism occur under physiological conditions and in the presence of cancer. Here, we detected an upregulation of purinergic signaling and nucleic acid metabolism in female lung adenocarcinoma patients, a finding not described in previous lung cancer studies. An NSLC-based study described an antitumor effect of the P2X4 receptor [63], which also exhibits sexual dimorphism in murine brain microglia [64]. Other P2 and A2 receptors play a role in NSCLC [63], but evidence of sex-based differences in receptor expression in human NSCLC patients has yet to be reported. Purinergic signaling and the role of purinergic receptors may also have relevance to innate and adaptive responses in different inflammatory and neurodegenerative diseases and several cancer types, including pancreatic ductal carcinoma (PDAC), hepatocellular, hepatobiliary carcinoma cells, and breast cancer [65,66,67,68,69]. Of note, studies have linked the upregulation of purinergic signaling with poor prognosis in PDAC [66]. We also discovered significant differences in lipid metabolism, with the positive regulation of ceramide and sphingolipid biosynthetic processes upregulated in female lung adenocarcinoma patients. The presence of lipids can promote tumorigenesis, while higher adipose tissue levels are associated with poorer outcomes in several cancers [37]. Thus, exploring the differential roles of purinergic signaling and lipid metabolism between male and female lung adenocarcinoma patients may represent an interesting proposition to improve sex-specific risk-stratification of patients, prevention, diagnosis, and treatment.
DNA damage and repair-related genes also presented sex-based differences in lung adenocarcinoma patients. In general, males present with a higher level of DNA damage, and females present with a lower DNA repair capacity [70], and, in agreement, we detected DNA repair and telomere protection as a downregulated functional group in female lung adenocarcinoma patients (with both mechanisms involved in tumor growth prevention [71]). Furthermore, we discovered the upregulation of DNA repair and an increase of DNA unwinding in male lung adenocarcinoma patients. DNA unwinding has emerged as a new target in cancer therapy with a primary focus on helicase inhibitors [72]. Besides, regarding DNA repair, the evaluation of poly (ADP-ribose) polymerase inhibitors in NSCLC cell lines has suggested potential therapeutic activity [73,74], and there may be value in exploring both treatment approaches in NSCLC patients, especially male patients.
With these results in mind, we propose future studies focused on DNA repair and lipid/purinergic metabolism in female lung adenocarcinoma patients and the immune response in male lung adenocarcinoma patients in the hope of developing enhanced therapeutic strategies.
Our study has characterized functional differences between the sexes in lung adenocarcinoma, shedding light on the functional basis behind this pathology in male and female patients. While our meta-analysis confirmed the conclusions of other studies, we also report previously undescribed alterations in biological processes that may broaden this field of study. Further knowledge regarding how those factors related to the functional mechanisms, described above, differentially impact male and female lung adenocarcinoma patients and may improve our understanding of the disease and improve treatment and diagnosis through biomarker identification.

4. Materials and Methods

Bioinformatics and statistical analysis employed R software v.3.5.3 [75]. Supplementary Table S6 details R packages and versions.

4.1. Study Search and Selection

Publicly available datasets were collected from GEO [21], ArrayExpress [76], and TCGA [22]. A systematic search of studies published in the period 2004–2018 was conducted in 2019 following the preferred reporting items for systematic reviews and meta-analyses (PRISMA) guidelines [32]. Two researchers involved in the study carried out the literature search, and the consistency of the review and selection procedures used was evaluated and confirmed. Several keywords were employed in the search, including lung adenocarcinoma (ADC), non-small-cell lung carcinoma (NSCLC), Homo sapiens, and excluding cell lines.
Eleven variables were considered for each study, including the clinical characteristics of the patients (e.g., sex and smoking habit) and experimental design (e.g., sample size and sample extraction source). The final inclusion criteria were:
  • Sex, disease stage, and smoking habit variables registered;
  • RNA extracted directly from human lung biopsies;
  • Both normal and lung adenocarcinoma samples available;
  • Patients who had not undergone treatment before biopsy;
  • Sample size of > 3 for case and control groups in both sexes.
Finally, normalized gene expression data of six array NSCLC datasets (GSE10072, GSE19188, GSE31210, GSE32863, GSE63459, and GSE75037) and counts matrix of three RNA-seq NSCLC datasets (GSE81089, GSE87340, and TCGA-LUAD) were retrieved.

4.2. Individual Transcriptomics Analysis

Individual transcriptomics analysis consisted of three steps: pre-processing, differential expression analysis, and functional enrichment analysis (Figure 5a).
Data pre-processing included the standardization of the nomenclature of the clinical variables included in each study, normalization of RNA-seq counts matrix, and exploratory analysis. RNA-seq counts were pre-processed with the edgeR [77] R package using the trimmed mean of m-values (TMM) method [78]. We assessed the normalization methods performed by the original authors for each dataset, and log2 transformed the matrices when necessary. Annotation from probe set to Entrez identifiers from the National Center for Biotechnology Information [79] database and gene symbol was carried out with the biomaRt [80] R package. When dealing with duplicated probe-to-Entrez mappings, the median of their expression values was calculated. The exploratory analysis included unsupervised clustering and PCA to detect patterns of expression between samples and genes and the presence of batch effects in each study (Figure 5b).
Differential expression analyses were performed using the limma [81] R package. To detect differentially expressed genes in male and female lung adenocarcinoma patients, the following contrast was applied:
(ADC.W − Control.W) − (ADC.M − Control.M)
where ADC.W, Control.W, ADC.M and Control.M correspond to lung adenocarcinoma affected women, control women, lung adenocarcinoma affected men, and control men, respectively. Paired samples design was implemented, and tobacco consumption was included as a batch effect on the limma linear model to reduce its impact on data. p-values were calculated and corrected for FDR [82]. This comparison allows the detection of genes and functions altered by the disease and that have higher or lower activity in women when compared to men. Significant functions and genes were considered when FDR < 0.05.
Functional enrichment analyses were performed using the Gene Set Enrichment Analysis (GSEA) implemented in the mdgsa [83] R package. p-values were, again, corrected for FDR. For functional annotation, two functional databases were used: the KEGG Pathway database [84] and GO [85]. GO terms were analyzed and propagated independently for each GO ontology: BPs, MFs, and cellular components (CC). Those annotations excessively specific or generic were filtered out, keeping functions with blocks of annotations between 10 and 500. Intersections within groups were analyzed with UpSet plots [86] (Figure 5c).

4.3. Functional Meta-Analysis

Functional GSEA results were integrated into a functional meta-analysis [87] implemented with mdgsa and metafor [88] R packages. Meta-analysis was applied under the DerSimonian and Laird random-effects model [89], taking into account individual study heterogeneity. This model considers the variability of individual studies by increasing the weights of studies with less variability when computing meta-analysis results. Thus, the most robust functions between studies are highlighted.
A total of 6467 GO BP terms, 785 GO CC terms, 1207 GO MF terms, and 213 KEGG pathways were evaluated. p-values, FDR corrected p-values, LOR, and 95% CIs of the LOR were calculated for each evaluated function. Functions and pathways with FDR < 0.05 were considered significant, and both funnel and forest plots were computed for each (Figure 5d). These representations were checked to assess for possible biased results, where LOR represents the effect size of a function, and the SE of the LOR serves as a study precision measure [90]. Sensitivity analysis (leave-one-out cross-validation [88]) was conducted for each significant function to verify possible alterations in the results due to the inclusion of any study.

4.4. Metafun-NSCLC Web Tool

All data and results generated in the different steps of the meta-analysis are available in the Metafun-NSCLC web tool (https://bioinfo.cipf.es/metafun-nsclc), which is freely accessible to any user and allows the confirmation of the results described in this manuscript and the exploration of other results of interest.
The front-end was developed using the Angular Framework. All graphics used in this web resource have been implemented with Plot.ly except for the exploratory analysis cluster plot, which was generated with ggplot2 [91].
This easy-to-use resource is divided into five sections: (1) Summary of analysis results in each phase. Then, for each of the studies, the detailed results of the (2) exploratory analysis, (3) differential expression, and (4) functional profiling. The user can interact with the web tool through graphics and tables and search for specific information for a gene or function. Finally, Section (5) provides several indicators for the significant functions identified in the meta-analysis that inform whether they are more active in men or women.

5. Conclusions

Sex-based molecular differences may influence the incidence and outcome of lung adenocarcinoma and, therefore, may have important clinical implications. We identified immune responses, purinergic signaling, and lipid-related processes as the main biological processes altered between male and female lung adenocarcinoma patients by a meta-analysis of transcriptomic datasets. Said processes exhibit increased activity in female lung adenocarcinoma patients, whereas other processes (such as DNA repair) are more active in male lung adenocarcinoma. Although further studies are required to verify and fully explore these findings, our results provide new clues to understand the molecular mechanisms of sex-based differences in lung adenocarcinoma patients and new perspectives regarding the identification of biomarkers and therapeutic targets.

Supplementary Materials

The following materials are available online at https://www.mdpi.com/2072-6694/13/1/143/s1, Figure S1: Distribution of sex information among the reviewed studies, Figure S2: Prognostic effect of transcriptional pathways, Table S1: Distribution of clinicopathological characteristics of each study population, Table S2: Summary of differential expression analysis results, Table S3: Significant genes differentially expressed between male and female lung adenocarcinoma patients, Table S4: All significant GO terms and KEGG pathways in the functional meta-analysis, Table S5: Genes of significant GO terms and KEGG pathways in functional meta-analysis, Table S6: Software and versions used in this work.

Author Contributions

Conceptualization, F.G.-G.; methodology, M.R.H. and F.G.-G.; software, P.M.-M., S.R.-G., I.P.-D., M.R.H. and F.G.-G.; validation, P.M.-M., A.R., R.F., M.R.H. and F.G.-G.; formal analysis, I.P.-D., M.R.H., S.R.-G. and F.G.-G.; investigation, I.P.-D., M.R.H., Z.A. and F.G.-G.; data curation, P.M.-M., I.P.-D.; writing—original draft preparation, I.P.-D., M.R.H. and F.G.-G.; writing—review and editing, I.P.-D., M.R.H., A.R., M.P., Z.A., R.F., M.d.l.I.-V. and F.G.-G.; visualization, P.M.-M., I.P.-D., S.R.-G., M.R.H. and F.G.-G.; supervision, M.R.H. and F.G.-G.; funding acquisition, F.G.-G.; project administration, F.G.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Fondo de Investigación Sanitaria (ISCIII PI15-00209), GV/2020/186, and ISCIII PT17/0009/0015 FEDER.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used for the analyses described in this work is publicly available at GEO (https://www.ncbi.nlm.nih.gov/geo/) and the Genomic Data Commons Data Portal (TCGA, https://portal.gdc.cancer.gov/). The accession numbers of the GEO datasets downloaded are: GSE10072, GSE19188, GSE31210, GSE32863, GSE63459, GSE75037, GSE81089, and GSE87340. Only the samples associated with the project TCGA-LUAD were downloaded from the TCGA.

Computer Code and Software

The code developed for the analyses described in this work is publicly available at GitLab (https://gitlab.com/ubb-cipf/metafunr). All software and versions used are detailed in Supplementary Table S6.

Acknowledgments

The authors thank the Principe Felipe Research Center (CIPF) for providing access to the cluster, which is co-funded by European Regional Development Funds (FEDER) in Valencian Community 2014–2020. The results published here are in whole or part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Jemal, A.; Miller, K.D.; Ma, J.; Siegel, R.L.; Fedewa, S.A.; Islami, F.; Devesa, S.S.; Thun, M.J. Higher Lung Cancer Incidence in Young Women Than Young Men in the United States. N. Engl. J. Med. 2018, 378, 1999–2009. [Google Scholar] [CrossRef]
  3. Barta, J.A.; Powell, C.A.; Wisnivesky, J.P. Global Epidemiology of Lung Cancer. Ann. Glob. Health 2019, 85, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2020. CA Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef] [PubMed]
  5. Conforti, F.; Pala, L.; Bagnardi, V.; De Pas, T.; Martinetti, M.; Viale, G.; Gelber, R.D.; Goldhirsch, A. Cancer immunotherapy efficacy and patients’ sex: A systematic review and meta-analysis. Lancet Oncol. 2018, 19, 737–746. [Google Scholar] [CrossRef]
  6. Barquín, M.; Calvo, V.; García-García, F.; Nuñez, B.; Sánchez-Herrero, E.; Serna-Blasco, R.; Auglytė, M.; Carcereny, E.; Rodriguez-Abreu, D.; López Castro, R.; et al. Sex is a strong prognostic factor in stage IV non-small-cell lung cancer patients and should be considered in survival rate estimation. Cancer Epidemiol. 2020, 67, 101737. [Google Scholar] [CrossRef]
  7. Wheatley-Price, P.; Blackhall, F.; Lee, S.-M.; Ma, C.; Ashcroft, L.; Jitlal, M.; Qian, W.; Hackshaw, A.; Rudd, R.; Booton, R.; et al. The influence of sex and histology on outcomes in non-small-cell lung cancer: A pooled analysis of five randomized trials. Ann. Oncol. 2010, 21, 2023–2028. [Google Scholar] [CrossRef]
  8. Fasco, M.J.; Hurteau, G.J.; Spivack, S.D. Gender-dependent expression of alpha and beta estrogen receptors in human nontumor and tumor lung tissue. Mol. Cell. Endocrinol. 2002, 188, 125–140. [Google Scholar] [CrossRef]
  9. Stabile, L.P.; Davis, A.L.G.; Gubish, C.T.; Hopkins, T.M.; Luketich, J.D.; Christie, N.; Finkelstein, S.; Siegfried, J.M. Human Non-Small Cell Lung Tumors and Cells Derived from Normal Lung Express Both Estrogen Receptor α and β and Show Biological Responses to Estrogen. Cancer Res. 2002, 62, 2141–2150. [Google Scholar]
  10. Kawai, H.; Ishii, A.; Washiya, K.; Konno, T.; Kon, H.; Yamaya, C.; Ono, I.; Minamiya, Y.; Ogawa, J. Estrogen Receptor α and β are Prognostic Factors in Non–Small Cell Lung Cancer. Clin. Cancer Res. 2005, 11, 5084–5089. [Google Scholar] [CrossRef] [Green Version]
  11. Kligerman, S.; White, C. Epidemiology of Lung Cancer in Women: Risk Factors, Survival, and Screening. Am. J. Roentgenol. 2011, 196, 287–295. [Google Scholar] [CrossRef] [PubMed]
  12. Bain, C.; Feskanich, D.; Speizer, F.E.; Thun, M.; Hertzmark, E.; Rosner, B.A.; Colditz, G.A. Lung Cancer Rates in Men and Women With Comparable Histories of Smoking. JNCI J. Natl. Cancer Inst. 2004, 96, 826–834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. North, C.M.; Christiani, D.C. Women and Lung Cancer: What is New? Semin. Thorac. Cardiovasc. Surg. 2013, 25, 87–94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Travis, W.D.; Brambilla, E.; Noguchi, M.; Nicholson, A.G.; Geisinger, K.R.; Yatabe, Y.; Beer, D.G.; Powell, C.A.; Riely, G.J.; Van Schil, P.E.; et al. International Association for the Study of Lung Cancer/American Thoracic Society/European Respiratory Society International Multidisciplinary Classification of Lung Adenocarcinoma. J. Thorac. Oncol. 2011, 6, 244–285. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Araujo, J.M.; Prado, A.; Cardenas, N.K.; Zaharia, M.; Dyer, R.; Doimi, F.; Bravo, L.; Pinillos, L.; Morante, Z.; Aguilar, A.; et al. Repeated observation of immune gene sets enrichment in women with non-small cell lung cancer. Oncotarget 2016, 7, 20282–20292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Shi, K.; Li, N.; Yang, M.; Li, W. Identification of Key Genes and Pathways in Female Lung Cancer Patients Who Never Smoked by a Bioinformatics Analysis. J. Cancer 2019, 10, 51–60. [Google Scholar] [CrossRef] [Green Version]
  17. Yuan, Y.; Liu, L.; Chen, H.; Wang, Y.; Xu, Y.; Mao, H.; Li, J.; Mills, G.B.; Shu, Y.; Li, L.; et al. Comprehensive Characterization of Molecular Differences in Cancer between Male and Female Patients. Cancer Cell 2016, 29, 711–722. [Google Scholar] [CrossRef] [Green Version]
  18. Li, Y.; He, C.-L.; Li, W.-X.; Zhang, R.-X.; Duan, Y. Transcriptome analysis reveals gender-specific differences in overall metabolic response of male and female patients in lung adenocarcinoma. PLoS ONE 2020, 15, e0230796. [Google Scholar] [CrossRef] [Green Version]
  19. Normand, S.-L.T. Meta-analysis: Formulating, evaluating, combining, and reporting. Stat. Med. 1999, 18, 321–359. [Google Scholar] [CrossRef]
  20. Higgins, J.; Thomas, J. Cochrane Handbook for Systematic Reviews of Interventions, 2nd ed.; Wiley-Blackwell: Chichester, UK, 2019. [Google Scholar]
  21. Edgar, R.; Domrachev, M.; Lash, A.E. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30, 207–210. [Google Scholar] [CrossRef] [Green Version]
  22. Weinstein, J.N.; Collisson, E.A.; Mills, G.B.; Shaw, K.R.M.; Ozenberger, B.A.; Ellrott, K.; Shmulevich, I.; Sander, C.; Stuart, J.M. The Cancer Genome Atlas Pan-Cancer analysis project. Nat. Genet. 2013, 45, 1113–1120. [Google Scholar] [CrossRef] [PubMed]
  23. Landi, M.T.; Dracheva, T.; Rotunno, M.; Figueroa, J.D.; Liu, H.; Dasgupta, A.; Mann, F.E.; Fukuoka, J.; Hames, M.; Bergen, A.W.; et al. Gene Expression Signature of Cigarette Smoking and Its Role in Lung Adenocarcinoma Development and Survival. PLoS ONE 2008, 3, e1651. [Google Scholar] [CrossRef] [PubMed]
  24. Hou, J.; Aerts, J.; den Hamer, B.; van IJcken, W.; den Bakker, M.; Riegman, P.; van der Leest, C.; van der Spek, P.; Foekens, J.A.; Hoogsteden, H.C.; et al. Gene Expression-Based Classification of Non-Small Cell Lung Carcinomas and Survival Prediction. PLoS ONE 2010, 5, e10312. [Google Scholar] [CrossRef] [PubMed]
  25. Okayama, H.; Kohno, T.; Ishii, Y.; Shimada, Y.; Shiraishi, K.; Iwakawa, R.; Furuta, K.; Tsuta, K.; Shibata, T.; Yamamoto, S.; et al. Identification of Genes Upregulated in ALK-Positive and EGFR/KRAS/ALK-Negative Lung Adenocarcinomas. Cancer Res. 2012, 72, 100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Yamauchi, M.; Yamaguchi, R.; Nakata, A.; Kohno, T.; Nagasaki, M.; Shimamura, T.; Imoto, S.; Saito, A.; Ueno, K.; Hatanaka, Y.; et al. Epidermal Growth Factor Receptor Tyrosine Kinase Defines Critical Prognostic Genes of Stage I Lung Adenocarcinoma. PLoS ONE 2012, 7, e43923. [Google Scholar] [CrossRef]
  27. Selamat, S.A.; Chung, B.S.; Girard, L.; Zhang, W.; Zhang, Y.; Campan, M.; Siegmund, K.D.; Koss, M.N.; Hagen, J.A.; Lam, W.L.; et al. Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression. Genome Res. 2012, 22, 1197–1211. [Google Scholar] [CrossRef] [Green Version]
  28. Robles, A.I.; Arai, E.; Mathé, E.A.; Okayama, H.; Schetter, A.J.; Brown, D.; Petersen, D.; Bowman, E.D.; Noro, R.; Welsh, J.A.; et al. An Integrated Prognostic Classifier for Stage I Lung Adenocarcinoma Based on mRNA, microRNA, and DNA Methylation Biomarkers. J. Thorac. Oncol. 2015, 10, 1037–1048. [Google Scholar] [CrossRef] [Green Version]
  29. Girard, L.; Rodriguez-Canales, J.; Behrens, C.; Thompson, D.M.; Botros, I.W.; Tang, H.; Xie, Y.; Rekhtman, N.; Travis, W.D.; Wistuba, I.I.; et al. An Expression Signature as an Aid to the Histologic Classification of Non–Small Cell Lung Cancer. Clin. Cancer Res. 2016, 22, 4880. [Google Scholar] [CrossRef] [Green Version]
  30. Mezheyeuski, A.; Bergsland, C.H.; Backman, M.; Djureinovic, D.; Sjöblom, T.; Bruun, J.; Micke, P. Multispectral imaging for quantitative and compartment-specific immune infiltrates reveals distinct immune profiles that classify lung cancer patients. J. Pathol. 2018, 244, 421–431. [Google Scholar] [CrossRef]
  31. Sun, Z.; Wang, L.; Eckloff, B.W.; Deng, B.; Wang, Y.; Wampfler, J.A.; Jang, J.; Wieben, E.D.; Jen, J.; You, M.; et al. Conserved recurrent gene mutations correlate with pathway deregulation and clinical outcomes of lung adenocarcinoma in never-smokers. BMC Med. Genom. 2014, 7, 486. [Google Scholar] [CrossRef] [Green Version]
  32. Moher, D.; Liberati, A.; Tetzlaff, J.; Altman, D.G.; Group, T.P. Preferred Reporting Items for Systematic Reviews and Meta-Analyses: The PRISMA Statement. PLoS Med. 2009, 6, e1000097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Klein, S.L.; Flanagan, K.L. Sex differences in immune responses. Nat. Rev. Immunol. 2016, 16, 626–638. [Google Scholar] [CrossRef] [PubMed]
  34. Cook, M.B.; Dawsey, S.M.; Freedman, N.D.; Inskip, P.D.; Wichner, S.M.; Quraishi, S.M.; Devesa, S.S.; McGlynn, K.A. Sex Disparities in Cancer Incidence by Period and Age. Cancer Epidemiol. Biomark. Prev. 2009, 18, 1174–1182. [Google Scholar] [CrossRef] [Green Version]
  35. Lukey, M.J.; Katt, W.P.; Cerione, R.A. Targeting amino acid metabolism for cancer therapy. Drug Discov. Today 2017, 22, 796–804. [Google Scholar] [CrossRef] [Green Version]
  36. Poff, A.; Koutnik, A.P.; Egan, K.M.; Sahebjam, S.; D’Agostino, D.; Kumar, N.B. Targeting the Warburg effect for cancer treatment: Ketogenic diets for management of glioma. Semin. Cancer Biol. 2019, 56, 135–148. [Google Scholar] [CrossRef]
  37. Butler, L.M.; Perone, Y.; Dehairs, J.; Lupien, L.E.; de Laat, V.; Talebi, A.; Loda, M.; Kinlaw, W.B.; Swinnen, J.V. Lipids and cancer: Emerging roles in pathogenesis, diagnosis and therapeutic intervention. Adv. Drug Deliv. Rev. 2020. [Google Scholar] [CrossRef] [PubMed]
  38. Gaignard, P.; Savouroux, S.; Liere, P.; Pianos, A.; Thérond, P.; Schumacher, M.; Slama, A.; Guennoun, R. Effect of Sex Differences on Brain Mitochondrial Function and Its Suppression by Ovariectomy and in Aged Mice. Endocrinology 2015, 156, 2893–2904. [Google Scholar] [CrossRef]
  39. Mervic, L. Time Course and Pattern of Metastasis of Cutaneous Melanoma Differ between Men and Women. PLoS ONE 2012, 7, e32955. [Google Scholar] [CrossRef] [Green Version]
  40. Pal, S.K.; Hurria, A. Impact of Age, Sex, and Comorbidity on Cancer Therapy and Disease Progression. J. Clin. Oncol. 2010, 28, 4086–4093. [Google Scholar] [CrossRef]
  41. Sun, S.; Schiller, J.H.; Gazdar, A.F. Lung cancer in never smokers—A different disease. Nat. Rev. Cancer 2007, 7, 778–790. [Google Scholar] [CrossRef]
  42. Kim, A.M.; Tingen, C.M.; Woodruff, T.K. Sex bias in trials and treatment must end. Nature 2010, 465, 688–689. [Google Scholar] [CrossRef] [PubMed]
  43. Woitowich, N.C.; Beery, A.; Woodruff, T. A 10-year follow-up study of sex inclusion in the biological sciences. eLife 2020, 9, e56344. [Google Scholar] [CrossRef] [PubMed]
  44. Wilkinson, M.D.; Dumontier, M.; Aalbersberg, I.J.; Appleton, G.; Axton, M.; Baak, A.; Blomberg, N.; Boiten, J.-W.; da Silva Santos, L.B.; Bourne, P.E.; et al. The FAIR Guiding Principles for scientific data management and stewardship. Sci. Data 2016, 3, 160018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Vinay, D.S.; Ryan, E.P.; Pawelec, G.; Talib, W.H.; Stagg, J.; Elkord, E.; Lichtor, T.; Decker, W.K.; Whelan, R.L.; Kumara, H.M.C.S.; et al. Immune evasion in cancer: Mechanistic basis and therapeutic strategies. Semin. Cancer Biol. 2015, 35, S185–S198. [Google Scholar] [CrossRef]
  46. Altorki, N.K.; Markowitz, G.J.; Gao, D.; Port, J.L.; Saxena, A.; Stiles, B.; McGraw, T.; Mittal, V. The lung microenvironment: An important regulator of tumour growth and metastasis. Nat. Rev. Cancer 2019, 19, 9–31. [Google Scholar] [CrossRef] [Green Version]
  47. Ye, Y.; Jing, Y.; Li, L.; Mills, G.B.; Diao, L.; Liu, H.; Han, L. Sex-associated molecular differences for cancer immunotherapy. Nat. Commun. 2020, 11, 1779. [Google Scholar] [CrossRef] [Green Version]
  48. Zeng, D.-Q.; Yu, Y.-F.; Ou, Q.-Y.; Li, X.-Y.; Zhong, R.-Z.; Xie, C.-M.; Hu, Q.-G. Prognostic and predictive value of tumor-infiltrating lymphocytes for clinical therapeutic research in patients with non-small cell lung cancer. Oncotarget 2016, 7, 13765–13781. [Google Scholar] [CrossRef] [Green Version]
  49. Ye, S.-L.; Li, X.-Y.; Zhao, K.; Feng, T. High expression of CD8 predicts favorable prognosis in patients with lung adenocarcinoma: A cohort study. Medicine 2017, 96, e6472. [Google Scholar] [CrossRef]
  50. Li, S.; Wang, Z.; Li, X.-J. Notch signaling pathway suppresses CD8+ T cells activity in patients with lung adenocarcinoma. Int. Immunopharmacol. 2018, 63, 129–136. [Google Scholar] [CrossRef]
  51. Abbas, A.K.; Trotta, E.; Simeonov, D.; Marson, A.; Bluestone, J.A. Revisiting IL-2: Biology and therapeutic prospects. Sci. Immunol. 2018, 3, eaat1482. [Google Scholar] [CrossRef] [Green Version]
  52. Gao, Y.; Lu, J.; Zeng, C.; Yang, J.; Huang, B.; Zhang, N.; Li, L.; Fu, X. IL-10 suppresses IFN-γ-mediated signaling in lung adenocarcinoma. Clin. Exp. Med. 2020, 20, 449–459. [Google Scholar] [CrossRef] [PubMed]
  53. Vahl, J.M.; Friedrich, J.; Mittler, S.; Trump, S.; Heim, L.; Kachler, K.; Balabko, L.; Fuhrich, N.; Geppert, C.-I.; Trufa, D.I.; et al. Interleukin-10-regulated tumour tolerance in non-small cell lung cancer. Br. J. Cancer 2017, 117, 1644–1655. [Google Scholar] [CrossRef] [Green Version]
  54. Miotto, D.; Cascio, N.L.; Stendardo, M.; Querzoli, P.; Pedriali, M.; De Rosa, E.; Fabbri, L.M.; Mapp, C.E.; Boschetto, P. CD8+ T cells expressing IL-10 are associated with a favourable prognosis in lung cancer. Lung Cancer 2010, 69, 355–360. [Google Scholar] [CrossRef]
  55. Emmerich, J.; Mumm, J.B.; Chan, I.H.; LaFace, D.; Truong, H.; McClanahan, T.; Gorman, D.M.; Oft, M. IL-10 Directly Activates and Expands Tumor-Resident CD8+ T Cells without De Novo Infiltration from Secondary Lymphoid Organs. Cancer Res. 2012, 72, 3570. [Google Scholar] [CrossRef] [Green Version]
  56. Shang, G.-S.; Liu, L.; Qin, Y.-W. IL-6 and TNF-α promote metastasis of lung cancer by inducing epithelial-mesenchymal transition. Oncol. Lett. 2017, 13, 4657–4660. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Silva, E.M.; Mariano, V.S.; Pastrez, P.R.A.; Pinto, M.C.; Castro, A.G.; Syrjanen, K.J.; Longatto-Filho, A. High systemic IL-6 is associated with worse prognosis in patients with non-small cell lung cancer. PLoS ONE 2017, 12, e0181125. [Google Scholar] [CrossRef] [PubMed]
  58. Lewis, A.M.; Varghese, S.; Xu, H.; Alexander, H.R. Interleukin-1 and cancer progression: The emerging role of interleukin-1 receptor antagonist as a novel therapeutic agent in cancer treatment. J. Transl. Med. 2006, 4, 48. [Google Scholar] [CrossRef] [Green Version]
  59. Gottschlich, A.; Endres, S.; Kobold, S. Can we use interleukin-1β blockade for lung cancer treatment? Transl. Lung Cancer Res. 2018, 7, S160–S164. [Google Scholar] [CrossRef] [Green Version]
  60. Wu, F.; Xu, J.; Huang, Q.; Han, J.; Duan, L.; Fan, J.; Lv, Z.; Guo, M.; Hu, G.; Chen, L.; et al. The Role of Interleukin-17 in Lung Cancer. Mediat. Inflamm. 2016, 2016, 8494079. [Google Scholar] [CrossRef] [Green Version]
  61. Wang, C.; Qiao, W.; Jiang, Y.; Zhu, M.; Shao, J.; Ren, P.; Liu, D.; Li, W. Effect of sex on the efficacy of patients receiving immune checkpoint inhibitors in advanced non-small cell lung cancer. Cancer Med. 2019, 8, 4023–4031. [Google Scholar] [CrossRef] [Green Version]
  62. Conforti, F.; Pala, L.; Bagnardi, V.; Viale, G.; De Pas, T.; Pagan, E.; Pennacchioli, E.; Cocorocchio, E.; Ferrucci, P.F.; De Marinis, F.; et al. Sex-Based Heterogeneity in Response to Lung Cancer Immunotherapy: A Systematic Review and Meta-Analysis. JNCI J. Natl. Cancer Inst. 2019, 111, 772–781. [Google Scholar] [CrossRef]
  63. Vijayan, D.; Smyth, M.J.; Teng, M.W.L. Purinergic Receptors: Novel Targets for Cancer Immunotherapy. In Oncoimmunology: A Practical Guide for Cancer Immunotherapy; Zitvogel, L., Kroemer, G., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 115–141. [Google Scholar]
  64. Crain, J.M.; Nikodemova, M.; Watters, J.J. Expression of P2 nucleotide receptors varies with age and sex in murine brain microglia. J. Neuroinflamm. 2009, 6, 24. [Google Scholar] [CrossRef] [Green Version]
  65. Savio, L.E.B.; de Andrade Mello, P.; da Silva, C.G.; Coutinho-Silva, R. The P2X7 Receptor in Inflammatory Diseases: Angel or Demon? Front. Pharmacol. 2018, 9, 52. [Google Scholar] [CrossRef] [Green Version]
  66. Hu, L.-P.; Zhang, X.-X.; Jiang, S.-H.; Tao, L.-Y.; Li, Q.; Zhu, L.-L.; Yang, M.-W.; Huo, Y.-M.; Jiang, Y.-S.; Tian, G.-A.; et al. Targeting Purinergic Receptor P2Y2 Prevents the Growth of Pancreatic Ductal Adenocarcinoma by Inhibiting Cancer Cell Glycolysis. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2019, 25, 1318–1330. [Google Scholar] [CrossRef] [Green Version]
  67. Tak, E.; Jun, D.Y.; Kim, S.-H.; Park, G.-C.; Lee, J.; Hwang, S.; Song, G.-W.; Lee, S.-G. Upregulation of P2Y2 nucleotide receptor in human hepatocellular carcinoma cells. J. Int. Med. Res. 2016, 44, 1234–1247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Asif, A.; Khalid, M.; Manzoor, S.; Ahmad, H.; Rehman, A.U. Role of purinergic receptors in hepatobiliary carcinoma in Pakistani population: An approach towards proinflammatory role of P2X4 and P2X7 receptors. Purinergic Signal. 2019, 15, 367–374. [Google Scholar] [CrossRef]
  69. Xu, M.; Almasi, S.; Yang, Y.; Yan, C.; Sterea, A.M.; Rizvi Syeda, A.K.; Shen, B.; Richard Derek, C.; Huang, P.; Gujar, S.; et al. The lysosomal TRPML1 channel regulates triple negative breast cancer development by promoting mTORC1 and purinergic signaling pathways. Cell Calcium 2019, 79, 80–88. [Google Scholar] [CrossRef]
  70. Zhu, Y.; Shao, X.; Wang, X.; Liu, L.; Liang, H. Sex disparities in cancer. Cancer Lett. 2019, 466, 35–38. [Google Scholar] [CrossRef]
  71. Shay, J.W.; Wright, W.E. Senescence and immortalization: Role of telomeres and telomerase. Carcinogenesis 2005, 26, 867–874. [Google Scholar] [CrossRef]
  72. Datta, A.; Brosh, R.M. New Insights Into DNA Helicases as Druggable Targets for Cancer Therapy. Front. Mol. Biosci. 2018, 5, 59. [Google Scholar] [CrossRef]
  73. Luo, J.; Dai, X.; Hu, H.; Chen, J.; Zhao, L.; Yang, C.; Sun, J.; Zhang, L.; Wang, Q.; Xu, S.; et al. Fluzoparib increases radiation sensitivity of non-small cell lung cancer (NSCLC) cells without BRCA1/2 mutation, a novel PARP1 inhibitor undergoing clinical trials. J. Cancer Res. Clin. Oncol. 2020, 146, 721–737. [Google Scholar] [CrossRef]
  74. Abbotts, R.; Topper, M.J.; Biondi, C.; Fontaine, D.; Goswami, R.; Stojanovic, L.; Choi, E.Y.; McLaughlin, L.; Kogan, A.A.; Xia, L.; et al. DNA methyltransferase inhibitors induce a BRCAness phenotype that sensitizes NSCLC to PARP inhibitor and ionizing radiation. Proc. Natl. Acad. Sci. USA 2019, 116, 22609. [Google Scholar] [CrossRef]
  75. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2014. [Google Scholar]
  76. Athar, A.; Füllgrabe, A.; George, N.; Iqbal, H.; Huerta, L.; Ali, A.; Snow, C.; Fonseca, N.A.; Petryszak, R.; Papatheodorou, I.; et al. ArrayExpress update-from bulk to single-cell expression data. Nucleic Acids Res. 2019, 47, D711–D715. [Google Scholar] [CrossRef]
  77. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Robinson, M.D.; Oshlack, A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010, 11, R25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Maglott, D.; Ostell, J.; Pruitt, K.D.; Tatusova, T. Entrez Gene: Gene-centered information at NCBI. Nucleic Acids Res. 2005, 33, D54–D58. [Google Scholar] [CrossRef]
  80. Durinck, S.; Moreau, Y.; Kasprzyk, A.; Davis, S.; De Moor, B.; Brazma, A.; Huber, W. BioMart and Bioconductor: A powerful link between biological databases and microarray data analysis. Bioinformatics 2005, 21, 3439–3440. [Google Scholar] [CrossRef] [Green Version]
  81. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef]
  82. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  83. Montaner, D.; Dopazo, J. Multidimensional Gene Set Analysis of Genomic Data. PLoS ONE 2010, 5, e10348. [Google Scholar] [CrossRef] [Green Version]
  84. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  85. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene Ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [Green Version]
  86. Lex, A.; Gehlenborg, N.; Strobelt, H.; Vuillemot, R.; Pfister, H. UpSet: Visualization of Intersecting Sets. IEEE Trans. Vis. Comput. Graph. 2014, 20, 1983–1992. [Google Scholar] [CrossRef]
  87. García-García, F. Methods of Functional Enrichment Analysis in Genomic Studies; Universitat de València: València, Spain, 2016. [Google Scholar]
  88. Viechtbauer, W. Conducting Meta-Analyses in R with the metafor Package. J. Stat. Softw. 2010, 36, 104408. [Google Scholar] [CrossRef] [Green Version]
  89. DerSimonian, R.; Laird, N. Meta-analysis in clinical trials. Control. Clin. Trials 1986, 7, 177–188. [Google Scholar] [CrossRef]
  90. Sterne, J.A.C.; Egger, M. Funnel plots for detecting bias in meta-analysis: Guidelines on choice of axis. J. Clin. Epidemiol. 2001, 54, 1046–1055. [Google Scholar] [CrossRef]
  91. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016; ISBN 978-3-319-24277-4. [Google Scholar]
Cancers 13 00143 g001 550
Figure 1. Flow of information through the different phases of the systematic review, following PRISMA Statement guidelines [32].
Figure 1. Flow of information through the different phases of the systematic review, following PRISMA Statement guidelines [32].
Cancers 13 00143 g001
Cancers 13 00143 g002 550
Figure 2. Number of samples per study, divided by sex and experimental group (ADC: lung adenocarcinoma samples).
Figure 2. Number of samples per study, divided by sex and experimental group (ADC: lung adenocarcinoma samples).
Cancers 13 00143 g002
Cancers 13 00143 g003 550
Figure 3. The intersection of significant functions between studies. UpSet plots for (a) Gene Ontology (GO) biological process, (b) GO molecular functions, (c) GO cellular components, and (d) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. UpSet plots detailing the number of common elements among GO terms in our functional enrichment analysis. Horizontal bars indicate the number of significant elements in each study. The vertical bars indicate the common elements in the sets, indicated with dots under each bar. The single points represent the number of unique elements in each group.
Figure 3. The intersection of significant functions between studies. UpSet plots for (a) Gene Ontology (GO) biological process, (b) GO molecular functions, (c) GO cellular components, and (d) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. UpSet plots detailing the number of common elements among GO terms in our functional enrichment analysis. Horizontal bars indicate the number of significant elements in each study. The vertical bars indicate the common elements in the sets, indicated with dots under each bar. The single points represent the number of unique elements in each group.
Cancers 13 00143 g003
Cancers 13 00143 g004 550
Figure 4. Summary dot plot of GO biological processes (BP) meta-analysis results. Only those significant terms with a log odds ratio (LOR) over 0.4 are shown.
Figure 4. Summary dot plot of GO biological processes (BP) meta-analysis results. Only those significant terms with a log odds ratio (LOR) over 0.4 are shown.
Cancers 13 00143 g004
Cancers 13 00143 g005 550
Figure 5. Workflow and analysis design. (a) Summary of the analysis design followed in this work, (b) PCA plot and HC analysis as an example of exploratory analysis (to explore possible batch effects and to assure expected data behavior) performed at the pre-processing stage to assess for the integrity of the data, (c) example of UpSet plot as an intersection analysis for functional enrichment analysis results, and (d) examples of forest and funnel plots to assess meta-analysis results.
Figure 5. Workflow and analysis design. (a) Summary of the analysis design followed in this work, (b) PCA plot and HC analysis as an example of exploratory analysis (to explore possible batch effects and to assure expected data behavior) performed at the pre-processing stage to assess for the integrity of the data, (c) example of UpSet plot as an intersection analysis for functional enrichment analysis results, and (d) examples of forest and funnel plots to assess meta-analysis results.
Cancers 13 00143 g005
Table
Table 1. Studies selected after the systematic review.
Table 1. Studies selected after the systematic review.
StudyPlatformPublication
GSE10072Affymetrix Human Genome U133A Array[23]
GSE19188Affymetrix Human Genome U133 Plus 2.0 Array[24]
GSE31210Affymetrix Human Genome U133 Plus 2.0 Array[25,26]
GSE32863Illumina HumanWG-6 v3.0 Expression BeadChip[27]
GSE63459Illumina HumanRef-8 v3.0 Expression BeadChip[28]
GSE75037Illumina HumanWG-6 v3.0 Expression BeadChip[29]
GSE81089Illumina HiSeq 2500[30]
GSE87340Illumina HiSeq 2000[31]
TCGAIllumina HiSeq 2000[22]
Table
Table 2. Summary of functional enrichment analysis results by Gene Ontology functions (BP: Biological Process, MF: Molecular Functions, CC: Cellular Component), and KEGG pathways. “Up” terms are overrepresented in female lung adenocarcinoma patients, while “Down” terms are overrepresented in male lung adenocarcinoma patients.
Table 2. Summary of functional enrichment analysis results by Gene Ontology functions (BP: Biological Process, MF: Molecular Functions, CC: Cellular Component), and KEGG pathways. “Up” terms are overrepresented in female lung adenocarcinoma patients, while “Down” terms are overrepresented in male lung adenocarcinoma patients.
StudySignificant GO BP TermsSignificant GO MF TermsSignificant GO CC TermsSignificant KEGG Pathways
UpDownTotalUpDownTotalUpDownTotalUpDownTotal
GSE100722615317929837164056156
GSE1918871219033112031145
GSE31210212230004048210
GSE32863428514792853340418127633
GSE634590262603382735145
GSE75037245352801441814213522527
GSE810892137187411112
GSE873401786224030348136126127
TCGA294228522283058217091301747
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pérez-Díez, I.; Hidalgo, M.R.; Malmierca-Merlo, P.; Andreu, Z.; Romera-Giner, S.; Farràs, R.; de la Iglesia-Vayá, M.; Provencio, M.; Romero, A.; García-García, F. Functional Signatures in Non-Small-Cell Lung Cancer: A Systematic Review and Meta-Analysis of Sex-Based Differences in Transcriptomic Studies. Cancers 2021, 13, 143. https://doi.org/10.3390/cancers13010143

AMA Style

Pérez-Díez I, Hidalgo MR, Malmierca-Merlo P, Andreu Z, Romera-Giner S, Farràs R, de la Iglesia-Vayá M, Provencio M, Romero A, García-García F. Functional Signatures in Non-Small-Cell Lung Cancer: A Systematic Review and Meta-Analysis of Sex-Based Differences in Transcriptomic Studies. Cancers. 2021; 13(1):143. https://doi.org/10.3390/cancers13010143

Chicago/Turabian Style

Pérez-Díez, Irene, Marta R. Hidalgo, Pablo Malmierca-Merlo, Zoraida Andreu, Sergio Romera-Giner, Rosa Farràs, María de la Iglesia-Vayá, Mariano Provencio, Atocha Romero, and Francisco García-García. 2021. "Functional Signatures in Non-Small-Cell Lung Cancer: A Systematic Review and Meta-Analysis of Sex-Based Differences in Transcriptomic Studies" Cancers 13, no. 1: 143. https://doi.org/10.3390/cancers13010143

APA Style

Pérez-Díez, I., Hidalgo, M. R., Malmierca-Merlo, P., Andreu, Z., Romera-Giner, S., Farràs, R., de la Iglesia-Vayá, M., Provencio, M., Romero, A., & García-García, F. (2021). Functional Signatures in Non-Small-Cell Lung Cancer: A Systematic Review and Meta-Analysis of Sex-Based Differences in Transcriptomic Studies. Cancers, 13(1), 143. https://doi.org/10.3390/cancers13010143

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