Next Article in Journal
Genetic Alterations and Therapeutic Targeting of Philadelphia-Like Acute Lymphoblastic Leukemia
Previous Article in Journal
miRNA Profiling in the Chicken Liver under the Influence of Early Microbiota Stimulation with Probiotic, Prebiotic, and Synbiotic
Previous Article in Special Issue
Unveiling Sex-Based Differences in the Effects of Alcohol Abuse: A Comprehensive Functional Meta-Analysis of Transcriptomic Studies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Analysis of Sex Disparities in the Genetic Architecture of Lung and Colorectal Cancers

by
Alireza Nazarian
and
Alexander M. Kulminski
*
Biodemography of Aging Research Unit, Social Science Research Institute, Duke University, Durham, NC 27705, USA
*
Author to whom correspondence should be addressed.
Genes 2021, 12(5), 686; https://doi.org/10.3390/genes12050686
Submission received: 17 March 2021 / Revised: 28 April 2021 / Accepted: 29 April 2021 / Published: 1 May 2021
(This article belongs to the Special Issue Sex/Gender Differences in Health from Omics Approaches)

Abstract

:
Almost all complex disorders have manifested epidemiological and clinical sex disparities which might partially arise from sex-specific genetic mechanisms. Addressing such differences can be important from a precision medicine perspective which aims to make medical interventions more personalized and effective. We investigated sex-specific genetic associations with colorectal (CRCa) and lung (LCa) cancers using genome-wide single-nucleotide polymorphisms (SNPs) data from three independent datasets. The genome-wide association analyses revealed that 33 SNPs were associated with CRCa/LCa at P < 5.0 × 10−6 neither males or females. Of these, 26 SNPs had sex-specific effects as their effect sizes were statistically different between the two sexes at a Bonferroni-adjusted significance level of 0.0015. None had proxy SNPs within their ±1 Mb regions and the closest genes to 32 SNPs were not previously associated with the corresponding cancers. The pathway enrichment analyses demonstrated the associations of 35 pathways with CRCa or LCa which were mostly implicated in immune system responses, cell cycle, and chromosome stability. The significant pathways were mostly enriched in either males or females. Our findings provided novel insights into the potential sex-specific genetic heterogeneity of CRCa and LCa at SNP and pathway levels.

1. Introduction

Sex disparities have been long reported in various malignancies, with most cancers predominantly affecting males and having better survival and lower mortality rates in females [1,2,3,4,5]. Lung (LCa) and colorectal (CRCa) cancers are among the top three common malignancies in both males and females. They jointly comprised around 25.1% of new cancer cases in men and 17.6% in women in 2018. They were also among the leading causes of cancer-related deaths in 2018, accounting for around 31% and 23% of such deaths in males and females [6]. A study of the relative risks of different cancer types revealed that LCa and CRCa were among 32 other cancers that had significantly higher incidence rates in men across various geographical regions and gross domestic product (GDP) groups, with average male-to-female incidence rate ratios of 2.08 and 1.33, respectively [3]. Sex has also been suggested as a potential favorable prognostic factor for these cancers conferring better survival to female patients [2,7,8,9]. LCa and CRCa were reported to have male-to-female mortality ratios of 2.31 and 1.42, respectively, and worse survival in males with significant male-to-female hazards ratios of 1.17 and 1.08 after adjusting models for the age of subjects and stage of tumors [2].
In addition to the sex-dependent differences in incidence, survival, and mortality rates, LCa and CRCa have displayed some other clinical and histopathological sex disparities in tumor topology, clinical manifestations, aggression potentials, and responses to therapy [10,11,12]. For instance, several studies reported that the female-to-male ratios were >1 and <1 in right- and left-sided CRCa, respectively, which have different clinical manifestations and are genetically heterogeneous [12,13,14]. As another example, while adenocarcinoma was found to be the most common LCa subtype in both sexes in most populations, the proportion of adenocarcinomas to squamous cell carcinomas was different in the two sexes [15]. It was also suggested that women might be more susceptible than men to tobacco carcinogenesis [11]. Sex disparities have been reported in the metastasis patterns of LCa and CRCa as well. The major distant metastatic sites of LCa include brain, bone, liver, and adrenal glands; and CRCa mainly metastasizes into the liver, lung, and bone [16]. It was suggested that bone metastases from LCa affects females more frequently than males possibly due to a feminized bone microenvironment [17], however, LCa-related brain metastases mostly occur in male patients [18]. Additionally, while the proportions of advanced-stage (i.e., Duke C and D) right- and left-sided CRCa were not statistically different in females, male patients with right-sided CRCa were more likely to present at advanced-stages compared to those with left-sided tumors [19]. Another aspect of sex differences was highlighted in previous reports indicating that female patients better responded to surgical treatment and chemotherapy of non-small cell LCa and CRCa than male patients [4]. In addition, females were found more susceptible than males to develop chemotherapy-related side effects (e.g., stomatitis, hematologic, and gastrointestinal toxicity) [8,20].
These disparities tend to be attributed to different life expectancies of the two sexes, hormonal effects, environmental exposures, and lifestyle risk factors. In addition, they might be in part mediated by the genetic and epigenetic mechanisms [5,10,21,22]. LCa and CRCa have evident genetic bases and in most cases are caused sporadically as multi-factorial disorders [23]. The genomes of patients with LCa and microsatellite stable CRCa may harbor ~10 to ~200 somatic mutations, mostly single-base substitutions, and CRCa cases with microsatellite instability (MSI) may accumulate >500 somatic mutations in their genomes [24]. The mutation density in the genomes of cancer patients demonstrated sex disparities as the numbers of single nucleotide mutations were in general higher in males than female [4]. It has been reported that LCa-linked mutations in P53, K-RAS, and EGFR genes are more commonly found in women than in men [11]. Additionally, while P53, APC, and K-RAS mutations are more commonly found in left-sided CRCa (more common in males), MSI, BRAF mutations, and eNOS and EPHB4 overexpression are more common in right-sided tumors (more common in females) [12,25]. More prevalent MSI-high (MSI-H) tumors in women with right-sided CRCa may confer them a decreased metastasis propensity [19]. In addition, genome-wide association studies (GWAS) have discovered several sex-dependent LCa/CRCa-associated single-nucleotide polymorphisms (SNPs) and haplotypes [22,26,27].
Genetic factors were also suggested to moderate the gender-specific recurrence and survival rates [10,12,28,29,30,31]. For instance, in a study of stage-III CRCa patients who were treated with adjuvant chemotherapy, female patients with non-zinc-binding mutations in TP53 DNA-binding domain had the best 5-year survival compared to females with zinc-binding mutations or wildtype genotypes. The same impacts from TP53 mutations on 5-year survival were not observed in men [30]. In addition, a polymorphism in PLS3 gene was found to be a predictor of tumor recurrence time in female patients with stage-II/III CRCa receiving adjuvant chemotherapy [31]. In another study, it was reported that among inoperable non-small-cell LCa cases who received carboplatin and gemcitabine therapy, the lack of ERCC1 gene expression conferred a survival advantage to male patients but not to females [29]. Additionally, women with advanced-staged non-small-cell LCa were more responsive to erlotinib, an EGFR-receptor inhibitor, compared to men which was partly attributed to higher mutations in EGFR in female patients [10,28]. Genetic factors may differentially impact chemotherapy-related toxicity as well. For instance, while women are more prone to fluorouracil (5-FU)-based chemotherapy toxicity [20], a splice site mutation in DPYD gene was strongly associated with sever 5-FU toxicity in men [32].
Despite the well-recognized sex differences in cancers, the underlying mechanisms have not been fully discovered and such disparities have not been consistently addressed in cancer research [5,22]. The evident contributions of genetic mechanisms to such sex disparities warrants further investigations into the sex-specific genetic architecture of LCa and CRCa, in particular due to their potential genetic heterogeneity. Exploring sex-specific genetic contributors to LCa and CRCa may provide more comprehensive insights into their underlying biological processes which in turn may help implementing more effective personalized and sex-specific medical interventions [5,10,22,33]. Searching genome-wide associations databases [34,35] shows that the genetic analysis of LCa’s and CRCa’s sex disparities has not received proper attention in previous GWAS. In this study, we performed sex-stratified genome-wide analyses of LCa and CRCa using phenotype and genotype data from three independent datasets to investigate potential sex disparities in the genetic predisposition to these common cancers.

2. Materials and Methods

2.1. Study Participants

Data from three independent studies were used including: Cardiovascular Health Study (CHS) [36], Framingham Heart Study (FHS) [37,38], and Health and Retirement Study (HRS) [39]. In each dataset, the genetic analyses were performed separately in females (i.e., CRCa-F and LCa-F) and males (i.e., CRCa-M and LCa-M). The cases comprised of 211 and 237 females as well as 186 and 220 males with CRCa and LCa, respectively. Also, 8382 and 8354 unaffected females and 6312 and 6278 unaffected males were included as controls in the CRCa-F, LCa-F, CRCa-M, and LCa-M analyses, respectively. The cases and controls were identified either by the study researchers (FHS) or by decoding medical diagnoses (CHS) or Medicare claims (HRS) using the International Classification of Disease codes, Ninth revision (ICD-9). Our genetic analyses were performed on subjects of Caucasian ancestry as there were not sufficient samples from other ethnicities. Table S1 (Supplementary File 1) provides summary demographic information for these three datasets. Figure 1 displays an overview of the analysis steps and main findings of our study.

2.2. Genotype Data and Quality Control (QC)

Our study made use of ~2 million genotyped and imputed SNPs. The imputation process has been detailed in [40]. Low-quality data were first filtered out including: (1) SNPs with imputation r2 < 0.7, (2) SNPs with minor allele frequencies (MAF) <5%, (3) SNPs/subjects with missing rates >5%, (4) SNPs deviated from Hardy–Weinberg equilibrium at P < 1.0 × 10−6, and (5) SNPs and subjects/families with Mendel error rates >2% in the case of FHS which is a family-based study. QC was performed using PLINK package [41]. This resulted in ~1.3–1.7 million SNPs in the datasets under consideration (Supplementary File 1: Table S2).

2.3. GWAS

2.3.1. Genetic Models

Additive genetic models were fitted using PLINK package [41] to identify the association between SNPs and cancers of interest after adjustment for birth year, smoking history, and body mass index (BMI) of subjects, and the top 3–4 principal components of genotype data obtained by GENESIS R package [42]. To address the risk of inflation of type-I errors due to ignoring family structure [43], SNPs nominally (i.e., P < 0.05) associated with CRCa/LCa in FHS were reanalyzed by fitting generalized linear mixed models (using lme4 R package [44]) which contained family IDs as a random-effects covariate in addition to the fixed-effects covariates stated above [40,45]. The results of GWAS of each cancer from the three datasets under consideration were then combined through an inverse-variance meta-analysis after adjustment for genomic inflation (i.e., λ values). Meta-analysis was performed using GWAMA package [46].

2.3.2. Discovery and Replication Analyses

We followed a commonly used discovery-replication strategy considering each of CHS, FHS, and HRS as a discovery set and the other two datasets as its counterpart replication sets. An association signal was considered replicated if a SNP had P < 5.0 × 10−8 (i.e., genome-wide significance) or 5.0 × 10−8P < 5.0 × 10−6 (i.e., suggestive significance) [40] in GWAS of one dataset and P < 0.05 in another dataset, and had consist directions of associations in the discovery and replication sets. The SNPs that were not among the replicated set of SNPs but had significant P-values at genome-wide or suggestive significance levels in conducted meta-analyses constituted the meta-analysis set of significant SNPs.

2.3.3. Novel Associations

CRCa/LCa-associated SNPs were considered as newly detected cancer variants if they were not associated with CRCa/LCa at P < 5.0 × 10−6 by previous GWAS available at databases such as GRASP [34] and NHGRI-EBI GWAS catalog [35]. LDlink web-tool [47] was then used to search possible proxy variants for the newly detected SNPs in the CEU population (i.e., Utah Residents with Northern and Western European Ancestry). A proxy variant was defined as a SNP that was located within ±1 Mb of a newly detected CRCa/LCa-associated SNP, was in LD with it (i.e., significant 𝒳2 in LD test) and was previously associated with the same cancer at P < 5.0 × 10−6.

2.3.4. Sex-specific Associations

SNPs disparately associated with CRCa/LCa in males and females were further analyzed by contrasting SNPs effects between males and females to determine if their effects were sex-specific [48]:
χ 2 = ( b f b m ) 2 s e f 2 + s e m 2
where χ2 is the Wald’s Chi-square statistics, bf and bm are the SNP effects (i.e., the natural logarithm of odds ratios) in females and males, and sef and sem are their standard errors.

2.4. Pathway Enrichment Analysis

Pathway enrichment analyses were performed by the GSA-SNP2 package [49] using compound gene-based P-values, obtained according to the fastBAT method [50,51], to identify potential biological processes associated with the studied cancers in males and females. The canonical pathways from the Broad Institute gene set enrichment analysis (GSEA) [52] were considered as the reference pathways [53,54,55,56]. The significant pathways were determined at false discovery rates (FDR) [57] of 0.025 (CRCa-F and CRCa-M) and 0.05 (LCa-F and LCa-M) to keep the numbers of possible false-positive findings below one in each analyzed cancer.

3. Results

3.1. Fixed-Effects Covariates

Smoking history, birth year, and BMI were included as fixed-effects covariates in our GWAS models to address their potential confounding effects on SNPs effects estimates, particularly due to their different distributions between males and females (Supplementary File 1: Table S1). Our meta-analyses revealed that smoking history and birth year were associated with CRCa and LCa in both males and females (P < 1.69 × 10−2) and BMI was associated with CRCa in both sexes (P < 7.74 ×10−3). However, their effects were not statistically different when their odds ratios were compared between males and females (Supplementary File 1: Table S3).

3.2. GWAS

Figures S1–S8 (Supplementary File 1) display the Manhattan and QQ plots from our GWAS. The λ values were smaller than 1.036 in these analyses (Supplementary File 1: Table S2), indicating the adequacy of population structure control [58]. Table 1 and Table S4 (Supplementary File 1) contain summary and detailed information regarding significant associations detected in our GWAS. We found that five SNPs (i.e., rs7593032, rs11000463, and rs11000467 in LCa-F; and rs9579517 and rs56357430 in CRCa-M) were associated with cancers of interest in a discovery dataset at suggestive significance level (i.e., 5.0 × 10−8P < 5.0 × 10−6) and were replicated at P < 0.05 in a replication dataset with the same directions of effects. In addition, there were 28 SNPs which were associated with CRCa or LCa in conducted meta-analyses (PMETA = 3.21 × 10−7 to 4.98 × 10−6; PQ = 1.72× 10−1 to 9.75 × 10−1; and i2 values between 0 and 0.432). As seen in Table 1, there were several genes (i.e., GLRX3 (CRCa-F), PRKG1 (LCa-F), MPHOSPH8 (CRCa-M), LINC02039, MAP7, and GRIK1 (LCa-M)) to which multiple significant associations signals were mapped. SNPs mapped to each of these genes were in high LD (0.855 ≤ r2 ≤ 1 and D’=1) with each other in the CEU population [47] (Supplementary File 1: Table S5).
None of the 33 detected SNPs and their corresponding chromosomal regions had significant association signals in both sexes. Of these, 26 SNPs were sex-specific as their effect sizes (i.e., the natural logarithm of odds ratios) were statistically different between males and females at a Bonferroni-adjusted significance level of 0.0015 (i.e., 0.05/33) when compared by a Wald’s Chi-square test (Table 2).

3.3. Pathway Enrichment Analysis

Our analyses (Table 3) revealed that 11 and 13 pathways were significantly associated with LCa-F and LCa-M, respectively, at an FDR of 0.05. They were mainly involved in meiosis, chromosome maintenance and telomere/centromere organization, and DNA transcription. Of these, eight pathways were significant in both sexes, while three pathways in females and five pathways in males were specifically enriched in one sex. We also found that 19 pathways were associated with CRCa-M at an FDR of 0.025. They were mainly involved in immune system responses and signal transduction. No pathway was enriched in CRCa-F analyses at FDRs of 0.025 or 0.05, however, 2 pathways were significant at FDR of 0.2. These 2 pathways were among the 19 CRCa-M-associated pathways. The other 17 pathways detected in CRCa-M were male-specific. None of the detected pathways were associated with both CRCa and LCa.

4. Discussion

Almost all complex diseases (including many cancer types) have manifested sex disparities in epidemiological and clinical studies (e.g., in incidence/prevalence rates or disease severity) [21] which may be due to the hormonal effects, lifestyle risk factors, and genetic mechanisms, among others [5,10,21,22,40]. Investigating sex disparities in the genetic mechanisms underlying complex disorders may have translational impacts on medical interventions and has been stressed by the National Institutes of Health (NIH) [5,10,22,33].
In this study, we analyzed potential sex disparities in the genetic architectures of CRCa and LCa in three independent datasets which, to the best of our knowledge, were not used previously for the study of sex-specific genetic contributions to the cancer phenotypes of interest. Our GWAS revealed replicated association signals of five SNPs at P < 5.0 × 10−6 that were associated with CRCa in males or LCa in females. In addition, 28 other SNPs were significantly associated with CRCa or LCa at suggestive significance level in conducted meta-analyses (Table 1 and Table S4). None of the detected SNPs attained genome-wide significance in our study. This might be due to the insufficient sample sizes of this study or the heterogeneity of SNPs effects in the studied cohorts (i.e., CHS, FHS, and HRS). All these 33 SNPs were potentially novel markers for the studied cancers as their association signals were not reported by previous GWAS and, in addition, there were no proxy CRCa/LCa-associated SNPs within their ±1 Mb flanking regions. It should be noted that the significant associations detected in GWAS do not imply causality. Functional studies are needed to investigate whether the identified SNPs themselves or other variants in nearby chromosomal regions that are in high LD with these index SNPs contribute to the genetic architecture of the studied cancers. A literature review further delineated potential implications of these SNPs, their closest genes, and variants in nearby regions in CRCa and LCa. We found that the closest genes to these SNPs were not associated with the same cancer at P < 5.0 × 10−6 by previous GWAS [34,35], except for CSMD1 gene (corresponding to rs13261356 detected in LCa-M) that was previously associated with LCa at genome-wide significance level [34,59]. Therefore, they can be considered as potentially novel genes for the studied cancers. However, nine of these genes (i.e., KHDRBS3, PRKG1, ZRANB1, THRB, FSTL4, HTR1E, LINC01377, MIR4675, and GRIK1) were previously implicated in cancers at other sites (i.e., other than those in our study) at PGWAS < 5.0 ×10−6 [34,35], and 7 genes (i.e., SNX31, KHDRBS3, GLRX3, THRB, MAP7, ATP8B1, and GRIK1) were prognostically linked to other cancers at P < 0.001 (The Human Protein Atlas [60]: www.proteinatlas.org accessed on 2019–2020) (Table 1). In addition, 21 of 33 detected SNPs were located within nine chromosomal regions which were not previously associated with the same cancers at PGWAS < 5.0 × 10−6 (i.e., KHDRBS3/8q24.23 and GLRX3/10q26.3 (CRCa-F); PRKG1/10q21.1 (LCa-F); THRB/3p24.2, HTR1E/6q14.3, and MPHOSPH8/13q12.11 (CRCa-M); MAP7/6q23.3, ATP8B1/18q21.31, and GRIK1/21q21.3 (LCa-M)).
Notably, 26 SNPs had sex-specific effects as they were significant only in males or females and their odds ratios were statistically different between the two sexes (Table 2). None of these SNPs were among or in LD with previously reported sex-linked SNPs [61,62]. The sex-specific SNPs associations with CRCa and LCa may advance the understanding of the underlying mechanisms of these common cancers in the two sexes by guiding functional studies in the detected chromosomal regions. Such sex-specific genetic factors may have translational implications in the era of personalized medicine by providing more efficient and cost-effective sex-specific medical interventions.
Our pathway enrichment analyses (Table 3) revealed that several pathways were significantly associated with CRCa and LCa (19 and 16 pathways, respectively). Most of the CRCa-associated pathways were related to the immune system functions. The intact/dysfunctional immune system responses were previously implicated in preventing/promoting tumorigenesis of CRCa [63,64]. The LCa-associated pathways were mostly involved in DNA replication/transcription and chromosome stability whose potential roles in LCa were previously highlighted [65,66,67]. Sex disparities were also noticed at the pathway level as most of the significant pathways (i.e., 17 pathways in CRCa and 8 pathways in LCa) were sex-specifically associated with the cancers of interest.
Limitations. Our analyses were generally underpowered for detecting association signals of SNPs with very small effect sizes and/or low MAFs (e.g., <0.05). Analyzing datasets with larger sample sizes would provide more statistical power and may replicate some of the detected disparate associations at genome-wide significance level and discover additional sex-specific associations. In addition, investigating potential sex disparities that may exist in the genetic architecture of different stages and/or histopathologic subtypes of CRCa and LCa may increase our knowledge about the genetic heterogeneity of these common cancers, although this requires sufficiently large sample sizes and availability of the staging and histopathologic data for the analyzed patients.

5. Conclusions

Our genome-wide analyses revealed associations of 33 SNPs (mapped to 19 genes) with CRCa or LCa at suggestive significance levels which were significant in either males or females. None of these associations were reported by previous GWAS, and there were no proxy SNPs within ±1 Mb regions of the identified SNPs. Of these, 26 SNPs had sex-specific effects evidenced by significantly different effect sizes (i.e., the natural logarithm of odds ratios) between the two sexes. Our pathway enrichment analyses revealed 35 pathways, mainly involved in immune system functions, DNA replication/transcription, and chromosome stability, were associated with the studied cancers. Twenty-five of these pathways were significant in either males or females. The potential sex-specific contributions to the genetic architecture of CRCa and LCa identified in our study provided novel insights into the genetic heterogeneity of these common cancers, although they did not imply causality. Such sex-specific associations, if replicated in independent genome-wide studies and/or corroborated in functional studies, may have translational impacts on the medical interventions in CRCa and LCa.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12050686/s1, Supplementary File 1: Supporting Acknowledgment, Table S1. Demographic information about the analyzed datasets; Table S2. Numbers (and percentages) of analyzed SNPs along with the genomic inflation factors (λ) resulted from our genome-wide association analyses; Table S3. The odds ratios of fixed-effects covariates in females and males from conducted meta-analyses; Table S4. Cancer-associated SNPs from genome-wide association analyses; Table S5. Linkage disequilibrium measures among SNPs with significant association signals that were mapped to the same gene; Figure S1. Manhattan plot of the genome-wide association analyses of colorectal cancer in females (CRCa-F); Figure S2. QQ plot of the genome-wide association analyses of colorectal cancer in females (CRCa-F); Figure S3. Manhattan plot of the genome-wide association analyses of lung cancer in females (LCa-F); Figure S4. QQ plot of the genome-wide association analyses of lung cancer in females (LCa-F); Figure S5. Manhattan plot of the genome-wide association analyses of colorectal cancer in males (CRCa-M); Figure S6. QQ plot of the genome-wide association analyses of colorectal cancer in males (CRCa-M); Figure S7. Manhattan plot of the genome-wide association analyses of lung cancer in males (LCa-M); Figure S8. QQ plot of the genome-wide association analyses of lung cancer in males (LCa-M).

Author Contributions

A.N. and A.M.K. designed the study, A.N. analyzed data, A.N. and A.M.K. wrote the manuscript. Both authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Grants from the National Institute on Aging (P01AG043352, R01AG047310, R01AG061853, R01AG065477, and R01AG070488). The funders had no role in study design, data collection and analysis, decision to publish, or manuscript preparation. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Institutional Review Board Statement

The CHS, FHS, and HRS studies, whose data were analyzed here, were approved by the corresponding institutional review boards (IRBs) and were conducted after obtaining written informed consent from participants. The authors accessed data upon approval by the Duke University IRB.

Data Availability Statement

The CHS and FHS datasets as well as the HRS genetic data were obtained from dbGaP (https://www.ncbi.nlm.nih.gov/gap, accession numbers: phs000287.v5.p1 (CHS), phs000007.v28.p10 (FHS), and phs000428.v1.p1 (HRS)). Phenotypic HRS data were obtained from the University of Michigan and the Centers for Medicare and Medical Services (CMS) in 2015–2016 and downloaded from the HRS public website (http://hrsonline.isr.umich.edu/index.php?p=data).

Acknowledgments

This manuscript was prepared using limited access datasets obtained from dbGaP and the University of Michigan. Please also see the Supporting Acknowledgments in S File 1 regarding these datasets. The authors thank Konstantin G. Arbeev and Arseniy P. Yashkin for help with downloading and preparing data.

Conflicts of Interest

The authors declare no conflicts of Interest.

References

  1. 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]
  2. Cook, M.B.; McGlynn, K.A.; Devesa, S.S.; Freedman, N.D.; Anderson, W.F. Sex Disparities in Cancer Mortality and Survival. Cancer Epidemiol. Biomark. Prev. 2011, 20, 1629–1637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Edgren, G.; Liang, L.; Adami, H.-O.; Chang, E.T. Enigmatic Sex Disparities in Cancer Incidence. Eur. J. Epidemiol. 2012, 27, 187–196. [Google Scholar] [CrossRef] [PubMed]
  4. Li, C.H.; Haider, S.; Shiah, Y.-J.; Thai, K.; Boutros, P.C. Sex Differences in Cancer Driver Genes and Biomarkers. Cancer Res. 2018, 78, 5527–5537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Rubin, J.B.; Lagas, J.S.; Broestl, L.; Sponagel, J.; Rockwell, N.; Rhee, G.; Rosen, S.F.; Chen, S.; Klein, R.S.; Imoukhuede, P.; et al. Sex Differences in Cancer Mechanisms. Biol. Sex Differ. 2020, 11, 17. [Google Scholar] [CrossRef] [Green Version]
  6. 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 A Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [Green Version]
  7. Ferguson, M.K.; Skosey, C.; Hoffman, P.C.; Golomb, H.M. Sex-Associated Differences in Presentation and Survival in Patients with Lung Cancer. J. Clin. Oncol 1990, 8, 1402–1407. [Google Scholar] [CrossRef]
  8. Singh, S.; Parulekar, W.; Murray, N.; Feld, R.; Evans, W.K.; Tu, D.; Shepherd, F.A. Influence of Sex on Toxicity and Treatment Outcome in Small-Cell Lung Cancer. J. Clin. Oncol. 2005, 23, 850–856. [Google Scholar] [CrossRef]
  9. Yang, Y.; Wang, G.; He, J.; Ren, S.; Wu, F.; Zhang, J.; Wang, F. Gender Differences in Colorectal Cancer Survival: A Meta-Analysis. Int. J. Cancer 2017, 141, 1942–1949. [Google Scholar] [CrossRef]
  10. 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] [PubMed]
  11. Kiyohara, C.; Ohno, Y. Sex Differences in Lung Cancer Susceptibility: A Review. Gend. Med. 2010, 7, 381–401. [Google Scholar] [CrossRef] [PubMed]
  12. Kim, S.-E.; Paik, H.Y.; Yoon, H.; Lee, J.E.; Kim, N.; Sung, M.-K. Sex- and Gender-Specific Disparities in Colorectal Cancer Risk. World J. Gastroenterol. 2015, 21, 5167–5175. [Google Scholar] [CrossRef]
  13. Benedix, F.; Kube, R.; Meyer, F.; Schmidt, U.; Gastinger, I.; Lippert, H. Colon/Rectum Carcinomas (Primary Tumor) Study Group Comparison of 17,641 Patients with Right- and Left-Sided Colon Cancer: Differences in Epidemiology, Perioperative Course, Histology, and Survival. Dis. Colon Rectum 2010, 53, 57–64. [Google Scholar] [CrossRef]
  14. Hansen, I.O.; Jess, P. Possible Better Long-Term Survival in Left versus Right-Sided Colon Cancer—A Systematic Review. Dan Med. J. 2012, 59, A4444. [Google Scholar]
  15. Cheng, T.-Y.D.; Cramb, S.M.; Baade, P.D.; Youlden, D.R.; Nwogu, C.; Reid, M.E. The International Epidemiology of Lung Cancer: Latest Trends, Disparities, and Tumor Characteristics. J. Thorac. Oncol. 2016, 11, 1653–1671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Wilkinson, I.B.; Raine, T.; Wiles, K.; Goodhart, A.; Hall, C.; O’Neill, H. Oxford Handbook of Clinical Medicine, 10th ed.; Sinharay, R., Furmedge, D., Eds.; Oxford University Press: Oxford UK, 2019; ISBN 978-0-19-883490-8. [Google Scholar]
  17. Farach-Carson, M.C.; Lin, S.-H.; Nalty, T.; Satcher, R.L. Sex Differences and Bone Metastases of Breast, Lung, and Prostate Cancers: Do Bone Homing Cancers Favor Feminized Bone Marrow? Front. Oncol. 2017, 7, 163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Sun, T.; Warrington, N.M.; Rubin, J.B. Why Does Jack, and Not Jill, Break His Crown? Sex Disparity in Brain Tumors. Biol. Sex Differ. 2012, 3, 3. [Google Scholar] [CrossRef] [Green Version]
  19. Koo, J.H.; Jalaludin, B.; Wong, S.K.C.; Kneebone, A.; Connor, S.J.; Leong, R.W.L. Improved Survival in Young Women with Colorectal Cancer. Off. J. Am. Coll. Gastroenterol. 2008, 103, 1488–1495. [Google Scholar] [CrossRef]
  20. Sloan, J.A.; Goldberg, R.M.; Sargent, D.J.; Vargas-Chanes, D.; Nair, S.; Cha, S.S.; Novotny, P.J.; Poon, M.A.; O’Connell, M.J.; Loprinzi, C.L. Women Experience Greater Toxicity with Fluorouracil-Based Chemotherapy for Colorectal Cancer. J. Clin. Oncol. 2002, 20, 1491–1498. [Google Scholar] [CrossRef]
  21. Ober, C.; Loisel, D.A.; Gilad, Y. Sex-Specific Genetic Architecture of Human Disease. Nat. Rev. Genet. 2008, 9, 911. [Google Scholar] [CrossRef] [Green Version]
  22. Dorak, M.T.; Karpuzoglu, E. Gender Differences in Cancer Susceptibility: An Inadequately Addressed Issue. Front. Genet. 2012, 3, 268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Galvan, A.; Ioannidis, J.P.A.; Dragani, T.A. Beyond Genome-Wide Association Studies: Genetic Heterogeneity and Individual Predisposition to Cancer. Trends Genet. 2010, 26, 132–141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Vogelstein, B.; Papadopoulos, N.; Velculescu, V.E.; Zhou, S.; Diaz, L.A.; Kinzler, K.W. Cancer Genome Landscapes. Science 2013, 339, 1546–1558. [Google Scholar] [CrossRef] [PubMed]
  25. Ulivi, P.; Scarpi, E.; Chiadini, E.; Marisi, G.; Valgiusti, M.; Capelli, L.; Casadei Gardini, A.; Monti, M.; Ruscelli, S.; Frassineti, G.L.; et al. Right- vs. Left-Sided Metastatic Colorectal Cancer: Differences in Tumor Biology and Bevacizumab Efficacy. Int. J. Mol. Sci. 2017, 18, 1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Gaj, P.; Maryan, N.; Hennig, E.E.; Ledwon, J.K.; Paziewska, A.; Majewska, A.; Karczmarski, J.; Nesteruk, M.; Wolski, J.; Antoniewicz, A.A.; et al. Pooled Sample-Based GWAS: A Cost-Effective Alternative for Identifying Colorectal and Prostate Cancer Risk Variants in the Polish Population. PLoS ONE 2012, 7, e35307. [Google Scholar] [CrossRef] [Green Version]
  27. Timofeeva, M.N.; Hung, R.J.; Rafnar, T.; Christiani, D.C.; Field, J.K.; Bickeböller, H.; Risch, A.; McKay, J.D.; Wang, Y.; Dai, J.; et al. Influence of Common Genetic Variation on Lung Cancer Risk: Meta-Analysis of 14 900 Cases and 29 485 Controls. Hum. Mol. Genet. 2012, 21, 4980–4995. [Google Scholar] [CrossRef] [PubMed]
  28. Shepherd, F.A.; Rodrigues Pereira, J.; Ciuleanu, T.; Tan, E.H.; Hirsh, V.; Thongprasert, S.; Campos, D.; Maoleekoonpiroj, S.; Smylie, M.; Martins, R.; et al. Erlotinib in Previously Treated Non-Small-Cell Lung Cancer. N. Engl. J. Med. 2005, 353, 123–132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Holm, B.; Mellemgaard, A.; Skov, T.; Skov, B.G. Different Impact of Excision Repair Cross-Complementation Group 1 on Survival in Male and Female Patients with Inoperable Non-Small-Cell Lung Cancer Treated with Carboplatin and Gemcitabine. J. Clin. Oncol. 2009, 27, 4254–4259. [Google Scholar] [CrossRef]
  30. Warren, R.S.; Atreya, C.E.; Niedzwiecki, D.; Weinberg, V.K.; Donner, D.B.; Mayer, R.J.; Goldberg, R.M.; Compton, C.C.; Zuraek, M.B.; Ye, C.; et al. Association of TP53 Mutational Status and Gender with Survival after Adjuvant Treatment for Stage III Colon Cancer: Results of CALGB 89803. Clin. Cancer Res. 2013, 19, 5777–5787. [Google Scholar] [CrossRef] [Green Version]
  31. Ning, Y.; Gerger, A.; Zhang, W.; Hanna, D.L.; Yang, D.; Winder, T.; Wakatsuki, T.; Labonte, M.J.; Stintzing, S.; Volz, N.; et al. Plastin Polymorphisms Predict Gender- and Stage-Specific Colon Cancer Recurrence after Adjuvant Chemotherapy. Mol. Cancer Ther. 2014, 13, 528–539. [Google Scholar] [CrossRef] [Green Version]
  32. Schwab, M.; Zanger, U.M.; Marx, C.; Schaeffeler, E.; Klein, K.; Dippon, J.; Kerb, R.; Blievernicht, J.; Fischer, J.; Hofmann, U.; et al. Role of Genetic and Nongenetic Factors for Fluorouracil Treatment-Related Severe Toxicity: A Prospective Clinical Trial by the German 5-FU Toxicity Study Group. J. Clin. Oncol. 2008, 26, 2131–2138. [Google Scholar] [CrossRef] [PubMed]
  33. Clayton, J.A.; Collins, F.S. Policy: NIH to Balance Sex in Cell and Animal Studies. Nat. News 2014, 509, 282–283. [Google Scholar] [CrossRef] [PubMed]
  34. Leslie, R.; O’Donnell, C.J.; Johnson, A.D. GRASP: Analysis of Genotype-Phenotype Results from 1390 Genome-Wide Association Studies and Corresponding Open Access Database. Bioinformatics 2014, 30, i185–i194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. MacArthur, J.; Bowler, E.; Cerezo, M.; Gil, L.; Hall, P.; Hastings, E.; Junkins, H.; McMahon, A.; Milano, A.; Morales, J.; et al. The New NHGRI-EBI Catalog of Published Genome-Wide Association Studies (GWAS Catalog). Nucleic Acids Res. 2017, 45, D896–D901. [Google Scholar] [CrossRef] [PubMed]
  36. Fried, L.P.; Borhani, N.O.; Enright, P.; Furberg, C.D.; Gardin, J.M.; Kronmal, R.A.; Kuller, L.H.; Manolio, T.A.; Mittelmark, M.B.; Newman, A. The Cardiovascular Health Study: Design and Rationale. Ann. Epidemiol. 1991, 1, 263–276. [Google Scholar] [CrossRef]
  37. Dawber, T.R.; Meadors, G.F.; Moore, F.E. Epidemiological Approaches to Heart Disease: The Framingham Study. Am. J. Public Health Nations Health 1951, 41, 279–286. [Google Scholar] [CrossRef]
  38. Feinleib, M.; Kannel, W.B.; Garrison, R.J.; McNamara, P.M.; Castelli, W.P. The Framingham Offspring Study: Design and Preliminary Data. Prev. Med. 1975, 4, 518–525. [Google Scholar] [CrossRef]
  39. Sonnega, A.; Faul, J.D.; Ofstedal, M.B.; Langa, K.M.; Phillips, J.W.; Weir, D.R. Cohort Profile: The Health and Retirement Study (HRS). Int. J. Epidemiol. 2014, 43, 576–585. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Nazarian, A.; Yashin, A.I.; Kulminski, A.M. Genome-Wide Analysis of Genetic Predisposition to Alzheimer’s Disease and Related Sex Disparities. Alzheimer’s Res. Ther. 2019, 11, 5. [Google Scholar] [CrossRef] [Green Version]
  41. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  42. Conomos, M.P.; Miller, M.B.; Thornton, T.A. Robust Inference of Population Structure for Ancestry Prediction and Correction of Stratification in the Presence of Relatedness. Genet. Epidemiol. 2015, 39, 276–293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. McArdle, P.F.; O’Connell, J.R.; Pollin, T.I.; Baumgarten, M.; Shuldiner, A.R.; Peyser, P.A.; Mitchell, B.D. Accounting for Relatedness in Family Based Genetic Association Studies. Hum. Hered. 2007, 64, 234–242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models Using Lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  45. Nazarian, A.; Arbeev, K.G.; Kulminski, A.M. The Impact of Disregarding Family Structure on Genome-Wide Association Analysis of Complex Diseases in Cohorts with Simple Pedigrees. J. Appl. Genet. 2020, 61, 75–86. [Google Scholar] [CrossRef]
  46. Mägi, R.; Morris, A.P. GWAMA: Software for Genome-Wide Association Meta-Analysis. BMC Bioinform. 2010, 11, 288. [Google Scholar] [CrossRef] [Green Version]
  47. Machiela, M.J.; Chanock, S.J. LDlink: A Web-Based Application for Exploring Population-Specific Haplotype Structure and Linking Correlated Alleles of Possible Functional Variants. Bioinformatics 2015, 31, 3555–3557. [Google Scholar] [CrossRef]
  48. Allison, P.D. Comparing Logit and Probit Coefficients across Groups. Soc. Methods Res. 1999, 28, 186–208. [Google Scholar] [CrossRef]
  49. Yoon, S.; Nguyen, H.C.T.; Yoo, Y.J.; Kim, J.; Baik, B.; Kim, S.; Kim, J.; Kim, S.; Nam, D. Efficient Pathway Enrichment and Network Analysis of GWAS Summary Data Using GSA-SNP2. Nucleic Acids Res. 2018, 46, e60. [Google Scholar] [CrossRef] [Green Version]
  50. Yang, J.; Lee, S.H.; Goddard, M.E.; Visscher, P.M. GCTA: A Tool for Genome-Wide Complex Trait Analysis. Am. J. Hum. Genet. 2011, 88, 76–82. [Google Scholar] [CrossRef] [Green Version]
  51. Bakshi, A.; Zhu, Z.; Vinkhuyzen, A.A.E.; Hill, W.D.; McRae, A.F.; Visscher, P.M.; Yang, J. Fast Set-Based Association Analysis Using Summary Data from GWAS Identifies Novel Gene Loci for Human Complex Traits. Sci. Rep. 2016, 6, 32894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  54. Nishimura, D. BioCarta. Biotech Softw. Internet Rep. 2001, 2, 117–120. [Google Scholar] [CrossRef]
  55. Schaefer, C.F.; Anthony, K.; Krupa, S.; Buchoff, J.; Day, M.; Hannay, T.; Buetow, K.H. PID: The Pathway Interaction Database. Nucleic Acids Res. 2009, 37, D674–D679. [Google Scholar] [CrossRef]
  56. Fabregat, A.; Jupe, S.; Matthews, L.; Sidiropoulos, K.; Gillespie, M.; Garapati, P.; Haw, R.; Jassal, B.; Korninger, F.; May, B.; et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2018, 46, D649–D655. [Google Scholar] [CrossRef]
  57. 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]
  58. Wellcome Trust Case Control Consortium Genome-Wide Association Study of 14,000 Cases of Seven Common Diseases and 3000 Shared Controls. Nature 2007, 447, 661–678. [CrossRef] [Green Version]
  59. Luo, W.; Schork, N.J.; Marschke, K.B.; Ng, S.-C.; Hermann, T.W.; Zhang, J.; Sanders, J.M.; Tooker, P.; Malo, N.; Zapala, M.A.; et al. Identification of Polymorphisms Associated with Hypertriglyceridemia and Prolonged Survival Induced by Bexarotene in Treating Non-Small Cell Lung Cancer. Anticancer Res. 2011, 31, 2303–2311. [Google Scholar]
  60. Uhlen, M.; Zhang, C.; Lee, S.; Sjöstedt, E.; Fagerberg, L.; Bidkhori, G.; Benfeitas, R.; Arif, M.; Liu, Z.; Edfors, F.; et al. A Pathology Atlas of the Human Cancer Transcriptome. Science 2017, 357, 2507. [Google Scholar] [CrossRef] [Green Version]
  61. Galichon, P.; Mesnard, L.; Hertig, A.; Stengel, B.; Rondeau, E. Unrecognized Sequence Homologies May Confound Genome-Wide Association Studies. Nucleic Acids Res. 2012, 40, 4774–4782. [Google Scholar] [CrossRef] [Green Version]
  62. Boraska, V.; Jerončić, A.; Colonna, V.; Southam, L.; Nyholt, D.R.; Rayner, N.W.; Perry, J.R.B.; Toniolo, D.; Albrecht, E.; Ang, W.; et al. Genome-Wide Meta-Analysis of Common Variant Differences between Men and Women. Hum. Mol. Genet. 2012, 21, 4805–4815. [Google Scholar] [CrossRef]
  63. Colussi, D.; Brandi, G.; Bazzoli, F.; Ricciardiello, L. Molecular Pathways Involved in Colorectal Cancer: Implications for Disease Behavior and Prevention. Int. J. Mol. Sci. 2013, 14, 16365–16385. [Google Scholar] [CrossRef] [Green Version]
  64. Markman, J.L.; Shiao, S.L. Impact of the Immune System and Immunotherapy in Colorectal Cancer. J. Gastrointest. Oncol. 2015, 6, 208–223. [Google Scholar] [CrossRef]
  65. Eymin, B.; Gazzeri, S. Role of Cell Cycle Regulators in Lung Carcinogenesis. Cell Adh. Migr. 2010, 4, 114–123. [Google Scholar] [CrossRef]
  66. Jafri, M.A.; Ansari, S.A.; Alqahtani, M.H.; Shay, J.W. Roles of Telomeres and Telomerase in Cancer, and Advances in Telomerase-Targeted Therapies. Genome Med. 2016, 8, 69. [Google Scholar] [CrossRef] [Green Version]
  67. Zhang, W.; Mao, J.-H.; Zhu, W.; Jain, A.K.; Liu, K.; Brown, J.B.; Karpen, G.H. Centromere and Kinetochore Gene Misexpression Predicts Cancer Patient Survival and Response to Radiotherapy and Chemotherapy. Nat. Commun. 2016, 7, 12619. [Google Scholar] [CrossRef]
Figure 1. Schematic summary of the analysis steps and main findings of this study. Abbreviations: CHS = Cardiovascular Health Study; FHS = Framingham Heart Study; HRS = Health and Retirement Study; CRCa = colorectal cancer; LCa = lung cancer; SNP = single-nucleotide polymorphism; MAF = minor allele frequency; GWAS = genome-wide association study; LD = linkage disequilibrium.
Figure 1. Schematic summary of the analysis steps and main findings of this study. Abbreviations: CHS = Cardiovascular Health Study; FHS = Framingham Heart Study; HRS = Health and Retirement Study; CRCa = colorectal cancer; LCa = lung cancer; SNP = single-nucleotide polymorphism; MAF = minor allele frequency; GWAS = genome-wide association study; LD = linkage disequilibrium.
Genes 12 00686 g001
Table 1. Significant findings from genome-wide association analyses.
Table 1. Significant findings from genome-wide association analyses.
ChrGeneSNPPosA1PHRSPFHSPCHSEffectsFreqOR (se)PMETAPQi2N1Mb?Proxy?Region?Gene?Prognostic?
CRCa-F
8q22.3SNX31rs1078186100613251C4.19 × 10−49.48 × 10−42.21 × 10−1---0.5740.615 (0.056)1.60 × 10−64.49 × 10−108590NNSNO
8q24.23KHDRBS3rs118174020135843271CNA3.92 × 10−48.73 × 10−4?--0.9490.357 (0.062)1.41 × 10−68.49 × 10−104752NNNOO
10q26.3GLRX3rs3118492130435775A2.74 × 10−38.05 × 10−31.69 × 10−2+++0.2291.672 (0.168)4.15 × 10−69.71 × 10−108544NNNNO
10q26.3GLRX3rs372813130442933A2.89 × 10−36.08 × 10−31.68 × 10−2+++0.2331.661 (0.163)3.25 × 10−69.75 × 10−108592NNNNO
10q26.3GLRX3rs359063130452326T4.64 × 10−33.89 × 10−31.39 × 10−2+++0.2301.669 (0.164)2.86 × 10−69.14 × 10−108545NNNNO
10q26.3GLRX3rs307385130452402T3.61 × 10−33.50 × 10−31.39 × 10−2+++0.2301.684 (0.167)2.16 × 10−69.23 × 10−108543NNNNO
LCa-F
2p21LOC102723824rs7593032 *42122050T4.63 × 10−63.32 × 10−13.68 × 10−2---0.7130.631 (0.057)3.66 × 10−69.14 × 10−20.5828577NNSNN
5q13.2LINC02056rs4277572584687T2.43 × 10−41.97 × 10−28.54 × 10−2+++0.3251.561 (0.138)4.98 × 10−66.74 × 10−108516SNSNN
10q21.1PRKG1rs11000463 *52014359T4.60 × 10−63.77 × 10−27.65 × 10−1--+0.9030.558 (0.067)1.75 × 10−53.31 × 10−20.7078591NNNON
10q21.1PRKG1rs11000467 *52015312G4.67 × 10−62.87 × 10−27.65 × 10−1--+0.9050.548 (0.066)1.37 × 10−53.41 × 10−20.7048570NNNON
10q26.13ZRANB1rs76972397124961380T1.70 × 10−4NA1.43 × 10−3-?-0.9300.432 (0.063)9.62 × 10−76.72 × 10−105770NNGON
18q11.2HRH4rs48296224532843A2.41 × 10−22.04 × 10−34.52 × 10−3+++0.2621.607 (0.151)4.90 × 10−64.75 × 10−108585NNSNN
CRCa-M
3p24.2THRBrs5775157824486090T5.16 × 10−45.80 × 10−24.39 × 10−3---0.9110.483 (0.064)1.89 × 10−66.98 × 10−106491NNNOO
5q31.1FSTL4rs4631227133152748A1.81 × 10−45.77 × 10−21.49 × 10−2+++0.2931.663 (0.163)3.20 × 10−64.76 × 10−106481NNGON
6q14.3HTR1Ers7290725186515650C7.49 × 10−51.55 × 10−34.64 × 10−1---0.9370.446 (0.063)1.06 × 10−63.73 × 10−106496NNNON
13q12.11MPHOSPH8rs9579517 *19603453C2.43 × 10−68.91 × 10−12.62 × 10−2---0.9000.519 (0.066)8.88 × 10−63.31 × 10−20.7076490NNNNN
13q12.11MPHOSPH8rs56357430 *19632421C2.43 × 10−69.20 × 10−13.69 × 10−2-+-0.9000.528 (0.068)1.77 × 10−52.25 × 10−20.7366498NNNNN
LCa-M
5p15.33LINC01377rs126577423123661A5.98 × 10−32.66 × 10−32.27 × 10−2---0.8640.555 (0.062)3.67 × 10−69.63 × 10−106490NNGON
5q23.2LINC02039rs77914729126003571C6.07 × 10−31.33 × 10−32.67 × 10−2---0.9290.475 (0.065)2.67 × 10−68.09 × 10−106495NNSNN
5q23.2LINC02039rs75874914126003878C6.25 × 10−31.64 × 10−32.66 × 10−2---0.9290.478 (0.065)3.24 × 10−68.29 × 10−106496NNSNN
5q23.2LINC02039rs11954381126006435G5.32 × 10−31.73 × 10−32.66 × 10−2---0.9290.477 (0.065)2.81 × 10−68.63 × 10−106491NNSNN
6q23.3MAP7rs9399183136496350A1.03 × 10−57.56 × 10−21.33 × 10−1---0.8330.585 (0.060)3.33 × 10−62.19 × 10−10.3416494NNNNO
6q23.3MAP7rs3799451136508794A4.91 × 10−65.78 × 10−21.61 × 10−1---0.8520.565 (0.060)1.62 × 10−62.01 × 10−10.3776495NNNNO
6q23.3MAP7rs3799453136518050G4.70 × 10−65.78 × 10−22.42 × 10−1---0.8530.570 (0.061)2.46 × 10−61.72 × 10−10.4326495NNNNO
6q23.3MAP7rs3799454136518078T6.90 × 10−66.85 × 10−22.16 × 10−1---0.8350.586 (0.060)3.58 × 10−61.76 × 10−10.4256498NNNNO
6q23.3MAP7rs3799462136524611G1.10 × 10−55.78 × 10−22.42 × 10−1---0.8530.577 (0.062)4.49 × 10−62.12 × 10−10.3566498NNNNO
8p23.2CSMD1rs132613564605119C1.04 × 10−38.05 × 10−4NA--?0.9000.482 (0.065)3.23 × 10−67.69 × 10−105209GNGSON
10p12.31MIR4675rs1101212920545511T2.15 × 10−31.68 × 10−5NA--?0.9450.375 (0.060)3.21 × 10−72.14 × 10−10.3525207NNSON
18q21.31ATP8B1rs243703757795758A2.39 × 10−4NA4.50 × 10−3-?-0.9090.449 (0.066)3.95 × 10−68.71 × 10−104216NNNNO
21q21.3GRIK1rs36343329588811C8.08 × 10−41.03 × 10−21.20 × 10−2---0.8630.539 (0.061)1.14 × 10−69.30 × 10−106490NNNOO
21q21.3GRIK1rs36343229591270G3.23 × 10−31.03 × 10−21.12 × 10−2---0.8630.554 (0.063)4.10 × 10−69.03 × 10−106495NNNOO
21q21.3GRIK1rs283240529591425T5.58 × 10−41.24 × 10−21.98 × 10−2---0.8540.548 (0.061)1.50 × 10−69.16 × 10−106497NNNOO
21q21.3GRIK1rs36347229595793G1.77 × 10−31.29 × 10−21.78 × 10−2---0.8560.562 (0.062)4.23 × 10−69.26 × 10−106490NNNOO
Abbreviations: SNP = single-nucleotide polymorphism; CRCa-F = colorectal cancer in females; LCa-F = lung cancer in females; CRCa-M = colorectal cancer in males; LCa-M = lung cancer in males; Chr = chromosomal region (i.e., cytogenetic band); Gene = the gene within which the SNP is located or else the closest gene to the SNP; Pos = SNP position based on Human Genome version 38 (hg38); A1 = effect allele; PHRS, PFHS, and PCHS = P-values in HRS (Health and Retirement Study), FHS (Framingham Heart Study), and CHS (Cardiovascular Health Study) datasets, respectively; Effects = directions of SNP’s effects in the aforementioned datasets (Positive, Negative, Missing); Freq = frequency of effect allele in conducted meta-analysis; OR (se) = odds ratio and its standard error in conducted meta-analysis; PMETA = P-value in conducted meta-analysis; PQ = P-value of Q-statistics (Cochran’s heterogeneity test); i2 = i-squared inconsistency metric; N = Number of subjects with non-missing genotype data in conducted meta-analysis; 1Mb? = whether any polymorphism associated with the same cancer was previously discovered within the ±1 Mb of the SNP detected here (G: polymorphism with P < 5.0 ×10−8, and S: polymorphism with 5.0 × 10−8P < 5.0 × 10−6, N: none); Proxy? = whether any of the 1Mb?-SNPs is in linkage disequilibrium (LD) with the SNP detected here and its P-value is less than the one detected in this study (Y: yes, N: no); Region? = whether any SNP within the chromosomal region corresponding to the SNP detected here was previously associated with the same cancer at P < 5.0 × 10−6 (Y: yes, N: no); Gene? = whether any SNP mapped to the gene corresponding to the SNP detected here was previously associated with cancers at P < 5.0 × 10−6 (S: same cancer, O: other cancers, N: none); Prognostic? = whether the closest gene to the SNP detected here was previously reported as a cancer prognostic factor at P < 0.001 (S: same cancer, O: other cancers, N: none); * = Replicated SNPs.
Table 2. Wald’s Chi-square test to compare odds ratios of detected cancer-associated SNPs between males and females.
Table 2. Wald’s Chi-square test to compare odds ratios of detected cancer-associated SNPs between males and females.
ChrGeneSNPPosA1MalesFemalesComparison
FreqORseP-ValueFreqORseP-ValueChi-SquareP-Value
CRCa-F
8q22.3SNX31rs1078186100613251C0.5680.9890.0979.19 × 10−10.5740.6150.0561.60 × 10−618.0332.17 × 10−5
8q24.23KHDRBS3rs118174020 +135843271CNANANANA0.9490.3570.0621.41 × 10−6NANA
10q26.3GLRX3rs3118492 +130435775A0.2351.0320.1198.07 × 10−10.2291.6720.1684.15 × 10−65.5071.89 × 10−2
10q26.3GLRX3rs372813 +130442933A0.2381.0370.1187.79 × 10−10.2331.6610.1633.25 × 10−65.4971.90 × 10−2
10q26.3GLRX3rs359063 +130452326T0.2351.0610.1216.43 × 10−10.2301.6690.1642.86 × 10−64.9352.63 × 10−2
10q26.3GLRX3rs307385 +130452402T0.2341.0660.1226.21 × 10−10.2301.6840.1672.16 × 10−64.9232.65 × 10−2
LCa-F
2p21LOC102723824rs759303242122050T0.7111.2060.1229.73 × 10−20.7130.6310.0573.66 × 10−623.0961.54 × 10−6
5q13.2LOC102477328rs42775 +72584687T0.3260.9580.0936.92 × 10−10.3251.5610.1384.98 × 10−68.5423.47 × 10−3
10q21.1PRKG1rs1100046352014359T0.9031.2660.1992.09 × 10−10.9030.5580.0671.75 × 10−515.2459.44 × 10−5
10q21.1PRKG1rs1100046752015312G0.9041.2520.2012.43 × 10−10.9050.5480.0661.37 × 10−515.2839.25 × 10−5
10q26.13ZRANB1rs76972397124961380T0.9350.9450.1938.27 × 10−10.9300.4320.0639.62 × 10−714.9281.12 × 10−4
18q11.2HRH4rs482962+24532843A0.2521.1320.1162.76 × 10−10.2621.6070.1514.90 × 10−63.3906.56 × 10−2
CRCa-M
3p24.2THRBrs5775157824486090T0.9110.4830.0641.89 × 10−60.9091.1310.1755.03 × 10−120.8474.97 × 10−6
5q31.1FSTL4rs4631227133152748A0.2931.6630.1633.20 × 10−60.2930.8920.0903.09 × 10−111.1498.41 × 10−4
6q14.3HTR1Ers7290725186515650C0.9370.4460.0631.06 × 10−60.9350.7740.1221.74 × 10−116.1335.91 × 10−5
13q12.11MPHOSPH8rs957951719603453C0.9000.5190.0668.88 × 10−60.9010.9190.136.11 × 10−115.3908.74 × 10−5
13q12.11MPHOSPH8rs5635743019632421C0.9000.5280.0681.77 × 10−50.9010.9150.1285.89 × 10−114.2941.56 × 10−4
LCa-M
5p15.33LOC102467074rs126577423123661A0.8640.5550.0623.67 × 10−60.8680.9750.1188.56 × 10−117.8092.44 × 10−5
5q23.2LOC102546228rs77914729126003571C0.9290.4750.0652.67 × 10−60.9321.0470.1738.16 × 10−118.3961.79 × 10−5
5q23.2LOC102546228rs75874914126003878C0.9290.4780.0653.24 × 10−60.9321.0950.1826.51 × 10−118.3071.88 × 10−5
5q23.2LOC102546228rs11954381126006435G0.9290.4770.0652.81 × 10−60.9321.1330.1915.41 × 10−118.4081.78 × 10−5
6q23.3MAP7rs9399183136496350A0.8330.5850.0603.33 × 10−60.8380.8630.0972.43 × 10−111.6696.36 × 10−4
6q23.3MAP7rs3799451136508794A0.8520.5650.0601.62 × 10−60.8540.8330.0951.57 × 10−111.8775.68 × 10−4
6q23.3MAP7rs3799453136518050G0.8530.5700.0612.46 × 10−60.8540.8320.0951.54 × 10−111.2877.81 × 10−4
6q23.3MAP7rs3799454136518078T0.8350.5860.0603.58 × 10−60.8400.8520.0952.03 × 10−111.0109.06 × 10−4
6q23.3MAP7rs3799462136524611G0.8530.5770.0624.49 × 10−60.8550.8290.0951.47 × 10−110.2911.34 × 10−3
8p23.2CSMD1rs132613564605119C0.9000.4820.0653.23 × 10−60.9030.8520.1283.68 × 10−115.8396.90 × 10−5
10p12.31MIR4675rs1101212920545511T0.9450.3750.0603.21 × 10−70.9451.4650.2881.23 × 10−121.5423.46 × 10−6
18q21.31ATP8B1rs243703757795758A0.9090.4490.0663.95 × 10−60.9101.0750.1857.30 × 10−119.8048.58 × 10−6
21q21.3GRIK1rs36343329588811C0.8630.5390.0611.14 × 10−60.8610.8820.1023.40 × 10−117.2063.35 × 10−5
21q21.3GRIK1rs36343229591270G0.8630.5540.0634.10 × 10−60.8610.8790.1023.26 × 10−114.8821.14 × 10−4
21q21.3GRIK1rs283240529591425T0.8540.5480.0611.50 × 10−60.8530.9240.1065.45 × 10−118.1971.99 × 10−5
21q21.3GRIK1rs36347229595793G0.8560.5620.0624.23 × 10−60.8550.8950.1023.89 × 10−115.0711.04 × 10−4
Please see the description provided below Table 1. + denotes the SNP did not have sex-specific effects (i.e., P-value ≥ 0.0015 in Chi-square test comparing SNP effects in males and females).
Table 3. Significant findings from pathway enrichment analyses.
Table 3. Significant findings from pathway enrichment analyses.
PathwayPathway SourceGSEA IDSizeMalesFemales
Countz-Scorep-Valueq-ValueCountz-Scorep-Valueq-Value
CRCa
Regulation of IFNA signalingREACTOMEM98224248.3300243.641.36 × 10−41.72 × 10−1
Antigen processing and presentationKEGGM1600489847.696.83 × 10−154.32 × 10−12NSNSNSNS
Natural killer cell mediated cytotoxicityKEGGM56691371356.719.61 × 10−124.05 × 10−9NSNSNSNS
Regulation of autophagyKEGGM638235355.677.09 × 10−92.24 × 10−6NSNSNSNS
TRAF6 mediated IRF7 activationREACTOMEM93630305.639.24 × 10−92.34 × 10−6NSNSNSNS
Autoimmune thyroid diseaseKEGGM1310353505.629.64 × 10−92.34 × 10−6NSNSNSNS
CD8/TCR downstream signalingPIDM27265654.973.27 × 10−75.91 × 10−5NSNSNSNS
Interferon alpha/beta signalingREACTOMEM97364624.788.82 × 10−71.40 × 10−4NSNSNSNS
Toll-like receptor signaling pathwayKEGGM32611021004.542.76 × 10−63.88 × 10−41003.512.26 × 10−41.72 × 10−1
RIG-I-like receptor signaling pathwayKEGGM1591371694.298.86 × 10−61.12 × 10−3NSNSNSNS
Translocation of ZAP-70 to immunological synapseREACTOMEM72214134.151.66 × 10−51.91 × 10−3NSNSNSNS
JAK-STAT signaling pathwayKEGGM174111551464.072.32 × 10−52.44 × 10−3NSNSNSNS
Immunoregulatory interactions between a lymphoid and a non-lymphoid cellREACTOMEM824070643.914.68 × 10−54.56 × 10−3NSNSNSNS
Generation of second messenger moleculesREACTOMEM1652327263.846.05 × 10−55.46 × 10−3NSNSNSNS
Cytosolic DNA-sensing pathwayKEGGM1184456553.691.10 × 10−49.28 × 10−3NSNSNSNS
Phosphorylation of CD3 and TCR zeta chainsREACTOMEM1249416153.681.18 × 10−49.33 × 10−3NSNSNSNS
Interferon signalingREACTOMEM9831591523.671.21 × 10−49.33 × 10−3NSNSNSNS
P2Y receptorsREACTOMEM1096012123.591.66 × 10−41.17 × 10−2NSNSNSNS
PD-1 signalingREACTOMEM1881018173.413.19 × 10−42.13 × 10−2NSNSNSNS
LCa
RNA polymerase I promoter openingREACTOMEM88462595.521.67 × 10−82.11 × 10−5595.783.83 × 10−94.84 × 10−6
Meiotic recombinationREACTOMEM101186824.836.77 × 10−74.28 × 10−4824.191.38 × 10−52.49 × 10−3
Amyloid fiber formationREACTOMEM107683794.681.42 × 10−64.97 × 10−4794.631.80 × 10−65.69 × 10−4
RNA polymerase I transcriptionREACTOMEM72889834.721.18 × 10−64.97 × 10−4834.251.08 × 10−52.28 × 10−3
Interleukin-7 signalingREACTOMEM54211114.464.04 × 10−61.02 × 10−3NSNSNSNS
MeiosisREACTOMEM5291161104.261.02 × 10−52.16 × 10−31103.472.61 × 10−43.31 × 10−2
Packaging of telomere endsREACTOMEM1769548484.191.42 × 10−52.57 × 10−3485.472.20 × 10−81.39 × 10−5
RNA polymerase I, RNA polymerase III, and mitochondrial transcriptionREACTOMEM8581221163.993.36 × 10−55.31 × 10−3NSNSNSNS
Meiotic synapsisREACTOMEM106173713.885.20 × 10−57.32 × 10−3714.318.14 × 10−62.06 × 10−3
A6B1 and A6B4 integrin signalingPIDM23946463.836.29 × 10−57.95 × 10−3NSNSNSNS
Interleukin-7 signalingBIOCARTAM129617173.621.49 × 10−41.71 × 10−2NSNSNSNS
Interleukin-2 family signalingREACTOMEM101241393.581.70 × 10−41.79 × 10−2NSNSNSNS
Deposition of new CENPA-containing nucleosomes at the centromereREACTOMEM87164613.442.88 × 10−42.81 × 10−2613.797.53 × 10−51.12 × 10−2
Telomere maintenanceREACTOMEM405275NSNSNSNS755.042.37 × 10−71.00 × 10−4
Systemic lupus erythematosusKEGGM4741140NSNSNSNS1333.817.09 × 10−51.12 × 10−2
Chromosome maintenanceREACTOMEM868122NSNSNSNS1173.43.42 × 10−43.94 × 10−2
Abbreviations: CRCa = colorectal cancer; LCa = lung cancer; GSEA = Gene Set Enrichment Analysis Platform; KEGG = Kyoto Encyclopedia of Genes and Genomes [53]; BIOCARTA = BIOCARTA pathways [54]; PID = Pathway Interaction Database [55]; REACTOME = REACTOME pathway knowledgebase [56]; Size = number of genes in the pathway; Count = number of enriched genes in the pathway; NS = non-significant.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nazarian, A.; Kulminski, A.M. Genome-Wide Analysis of Sex Disparities in the Genetic Architecture of Lung and Colorectal Cancers. Genes 2021, 12, 686. https://doi.org/10.3390/genes12050686

AMA Style

Nazarian A, Kulminski AM. Genome-Wide Analysis of Sex Disparities in the Genetic Architecture of Lung and Colorectal Cancers. Genes. 2021; 12(5):686. https://doi.org/10.3390/genes12050686

Chicago/Turabian Style

Nazarian, Alireza, and Alexander M. Kulminski. 2021. "Genome-Wide Analysis of Sex Disparities in the Genetic Architecture of Lung and Colorectal Cancers" Genes 12, no. 5: 686. https://doi.org/10.3390/genes12050686

APA Style

Nazarian, A., & Kulminski, A. M. (2021). Genome-Wide Analysis of Sex Disparities in the Genetic Architecture of Lung and Colorectal Cancers. Genes, 12(5), 686. https://doi.org/10.3390/genes12050686

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