Next Article in Journal
A Precision Medicine Approach to Optimize Modulator Therapy for Rare CFTR Folding Mutants
Next Article in Special Issue
Challenges to Building a Gene Variant Commons to Assess Hereditary Cancer Risk: Results of a Modified Policy Delphi Panel Deliberation
Previous Article in Journal
Evolution of CT Findings and Lung Residue in Patients with COVID-19 Pneumonia: Quantitative Analysis of the Disease with a Computer Automatic Tool
Previous Article in Special Issue
Somatic Tumor Profile Analysis in a Patient with Germline PMS2 Mutation and Synchronous Ovarian and Uterine Carcinomas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calcium-Sensing Receptor Polymorphisms at rs1801725 Are Associated with Increased Risk of Secondary Malignancies

by
Ky’Era V. Actkins
1,†,
Heather K. Beasley
2,†,
Annika B. Faucon
3,†,
Lea K. Davis
2,4 and
Amos M. Sakwe
2,*
1
Department of Microbiology, Immunology and Physiology, School of Graduate Studies and Research, Meharry Medical College, Nashville, TN 37208, USA
2
Department of Biochemistry, Cancer Biology, Neuroscience and Pharmacology, School of Graduate Studies and Research, Meharry Medical College, Nashville, TN 37208, USA
3
Vanderbilt University Medical Center, Vanderbilt Genetics Institute, Nashville, TN 37232, USA
4
Department of Medicine, Division of Genetic Medicine, Vanderbilt University Medical Center, Nashville, TN 37232, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally.
J. Pers. Med. 2021, 11(7), 642; https://doi.org/10.3390/jpm11070642
Submission received: 28 April 2021 / Revised: 7 June 2021 / Accepted: 2 July 2021 / Published: 6 July 2021
(This article belongs to the Special Issue Cancer Genetics and Genomics Research)

Abstract

:
Dysregulation of systemic calcium homeostasis during malignancy is common in most patients with high-grade tumors. However, it remains unclear whether single nucleotide polymorphisms (SNPs) that alter the sensitivity of the calcium-sensing receptor (CaSR) to circulating calcium are associated with primary and/or secondary neoplasms at specific pathological sites in patients of European and African ancestry. Multivariable logistic regression models were used to analyze the association of CASR SNPs with circulating calcium, parathyroid hormone, vitamin D, and primary and secondary neoplasms. Circulating calcium is associated with an increased risk for breast, prostate, and skin cancers. In patients of European descent, the rs1801725 CASR SNP is associated with bone-related cancer phenotypes, deficiency of humoral immunity, and a higher risk of secondary neoplasms in the lungs and bone. Interestingly, circulating calcium levels are higher in homozygous patients for the inactivating CASR variant at rs1801725 (TT genotype), and this is associated with a higher risk of secondary malignancies. Our data suggest that expression of CaSR variants at rs1801725 is associated with a higher risk of developing secondary neoplastic lesions in the lungs and bone, due in part to cancer-induced hypercalcemia and/or tumor immune suppression. Screening of patients for CASR variants at this locus may lead to improved management of high calcium associated tumor progression.

1. Introduction

The calcium-sensing receptor (CaSR) plays an essential role in systemic calcium homeostasis by sensing slight increases in circulating calcium levels. Besides its role in calcium homeostasis, the CaSR plays additional roles in several tissues, including several tumor types. The receptor is a tumor suppressor in certain cancers, e.g., neuroblastoma, gastric, colon, and parathyroid cancers [1,2], a tumor promoter in other cancers such as breast, prostate, ovarian, and kidney cancers [2,3,4], and plays a significant role in chemosensitivity of breast and other neoplasms [5,6,7]. For normal systemic calcium homeostasis, activation of the receptor by high calcium alters intracellular signaling pathways that reduce both parathyroid hormone (PTH) secretion by parathyroid chief cells and renal calcium reabsorption [8,9], which together restores systemic calcium to near normal levels.
During malignancy, calcium homeostasis is progressively disrupted by the secretion of osteolytic factors such as parathyroid hormone-related protein (PTHrP) by tumor cells, which like PTH promotes osteolysis. With tumor progression, the increased circulating PTHrP stimulates the biosynthesis and secretion of receptor activator of nuclear factor κB ligand (RANKL) by osteoblasts that further activate preosteoclasts osteoclast driven osteolysis [10]. This bone resorption signaling is inhibited by osteoprotegerin (OPG), which binds to RANKL and inhibits RANK signaling and bone resorption [11]. However, the release of calcium and other growth factors from bone resorption contributes to tumor growth and a vicious osteolytic cycle that ultimately results in cancer-induced hypercalcemia (CIH) or humoral hypercalcemia of malignancy [3,12,13,14]. CIH is a benign paraneoplastic syndrome common in patients with solid tumors and is associated with high mortality rates, with a median patient survival of 1–3 months from initial diagnosis [15,16]. Interestingly, high circulating calcium levels are associated with aggressive breast tumors in premenopausal women and larger breast tumors in postmenopausal women [17,18]. Although this suggests the regulation of tumor progression and/or metastasis by high calcium, whether distinct variants of the CaSR mediate the effects of high calcium remains poorly understood.
Several variants and single nucleotide polymorphisms (SNPs) have been reported in the CASR gene [2,19,20]. Some of these SNPs are inactivating (loss-of-function) and are associated with hypercalcemia syndromes, while others are gain-of-function or activating variants associated with hypocalcemia [21]. Specifically, missense polymorphisms at rs1801725 (e.g., A986S) and rs1801726 (e.g., Q1011E) in exon 7 of the receptor are differentially associated with calcium in various chronic diseases [22,23,24], some cancers [25,26], and various heart diseases [27]. In several other studies, the A986S CaSR variant at rs1801725 is common among subjects of European ancestry, while the Q1011E variant at rs1801726 is common in individuals of African descent [22,25,26,28,29,30]. In these studies, the A986S CaSR variant is consistently associated with calcium, while the association of the Q1011E CaSR variant with calcium has been contradictory. However, whether the expression of polymorphic CASR variants at these loci, which presumably affects the sensitivity of the receptor to calcium, influences cancer diagnosis at specific pathological sites in patients of European versus African descent remains unclear.
In this study, we proposed that the expression of polymorphic CASR variants at rs1801725 and at rs1801726 loci are associated with circulating calcium, PTH, Vitamin D and that these factors independently influence cancer diagnosis at specific pathological sites in patients of European versus African descent. Our analyses reveal that circulating calcium, PTH, and vitamin D levels are associated with neoplasms at multiple pathological sites, including breast, prostate, and skin cancers. In cancer patients, CASR polymorphisms at rs1801725, but not at rs1801726, are associated with circulating calcium levels, deficiency of humoral immunity, and an increased risk of secondary neoplastic lesions in the lungs and bone. Our data suggest that the CaSR variant at rs1801725, is associated with a higher risk of developing secondary neoplastic lesions in the lungs and bone due in part to cancer-induced hypercalcemia and/or tumor immune suppression. Screening of patients for CASR variants at this locus may lead to improved management of high calcium associated tumor progression.

2. Materials and Methods

2.1. Patient Samples

The Vanderbilt University Medical Center (VUMC, Nashville, TN, USA) maintains a clinical research warehouse known as the Synthetic Derivative (SD), which includes a de-identified mirror image of the VUMC EHR available for research purposes [31]. The SD contains over 3 million records with data on demographics, billing codes according to the International Classification of Diseases (ICD 9th and 10th revisions), procedural codes, clinical notes, medications, and laboratory measurements. In addition, over 250,000 DNA samples are linked to these records in a biorepository (BioVU).
Inclusion criteria for this study were EHRs from individuals with at least 5 ICD codes of any type on separate days over a 3-year period, which defined a “data floor” to select a population of VUMC patients enriched for primary care. ICD9 and ICD10 codes were then mapped to hierarchical phenotype classifications (i.e., “phecodes”) using the phecode map version 1.2 and 1.2b1 [32]. Phecodes were used to identify primary cancer phenotypes at specific pathological sites (e.g., “Cancer of mouth,” “Cancer of the esophagus,” etc.) and secondary malignancies (Table S1). Cases were identified by the presence of at least one cancer phecode (which required at least two corresponding ICD codes) and controls were those who did not have any cancer phecodes. Two ancestry-specific datasets were subsequently established (Table 1), one composed of individuals of European descent herein designated European dataset (n = 53,682), and another containing predominantly individuals of African descent designated African dataset (n = 10,777).

2.2. Clinical Laboratory Measurements

All clinical laboratory measurements (labs), including total calcium, ionized calcium, intact parathyrin (parathyroid hormone, PTH), and 25-hydroxy vitamin D (herein referred to as vitamin D), were extracted from the SD and cleaned using the recently published QualityLab pipeline [33]. In brief, we were restricted to lab tests for which at least 70% of the observations were measured in the same unit and required that each lab was measured on a minimum of 100 patients resulting in a minimum of 1000 numeric observations recorded. We then applied lab-specific quality control filters to remove infinite and non-numeric values, as well as observations outside of 4 standard deviations from the overall sample mean, indicative of biologically implausible values due to technical or recording errors, monogenic disorders, or extreme environmental influence. Labs were adjusted for the cubic splines of age (4-knots) and transformed by inverse normalization for analysis by linear and logistic regression models. Median labs in the European and African datasets were then stratified into pre-diagnosis and post-diagnosis values, based on the earliest recorded cancer diagnosis date.

2.3. Genotyping and Quality Control Pipeline

Genotyping was performed by using the Illumina MEGAEX array platform [34]. PLINK v1.9 was used to filter SNPs (<0.98) and individual subjects (<0.98) with low call rates, sex discrepancies, and excessive or significantly reduced heterozygosity (Fhet > 0.2). Principal component (PC) analysis was performed using flashPCA2 was used to select individuals of European ancestry and, separately, of African ancestry [35]. Minor allele frequency was calculated for CASR rs1801725 (c.2956G > T, p.Ala986Ser, NM_000388.4) and rs1801726 (c.3031C > G, p.Gln1011Glu, NM_000245.4) using PLINK v1.9.

2.4. Statistical Analysis

We first examined the relationship between the clinical labs (circulating calcium, vitamin D, and PTH levels) and cancer status (primary and secondary cancer diagnoses) after accounting for covariates including median age across the medical record, EHR-reported race, and sex. Each lab was then treated as an independent variable in separate multivariable logistic regression models of twenty-two primary and seven secondary cancer diagnoses for a total of 116 models tested.
Next, we tested the relationship between CASR SNPs (rs1801725 and rs1801726) and circulating calcium, ionized calcium, vitamin D, and PTH levels after adjusting for sex and the top 10 principal components (PCs) calculated from the genetic data (to adjust for population stratification) in separate multivariable linear regression models. The analysis was performed separately in European and African datasets. A total of 16 regression models were tested. In a subsequent sensitivity analysis, clinical labs were further stratified into pre- and post-cancer diagnoses according to the lab draw date.
Finally, for each SNP, we tested both an additive (rs1801725: GG vs. GT vs. TT, rs1801726: CC vs. CG vs. GG) and a recessive (rs1801725: TT vs. (GG + GT), rs1801726: GG vs. (CC + CG)) genetic model in a logistic regression framework to evaluate the relationship between the CASR variants and seven secondary malignant cancers in the European dataset resulting in a total of 14 models tested. Median age across the medical record, sex, and the top 10 PCs were included as covariates in genetic analyses of lab traits and cancer diagnoses.

2.5. Statistical Power and Correction for Multiple Comparisons

The primary analyses in this report included a total of 148 independent statistical tests, reflecting a priori derived hypotheses. A strict Bonferroni correction of p < 3.38 × 10−4 (0.05/148) was used to determine the statistical significance of primary results from four sets of regression analyses (i.e., testing associations between primary cancer and labs, CASR SNPs and labs, CASR SNPs, and primary cancer, and CASR SNPs and secondary malignancies). The exploratory phenome-wide association study (PheWAS), a hypothesis generating analysis, was not included in the multiple testing correction.

3. Results

3.1. Characteristics of European and African Datasets

The sample demographics of each of our datasets are shown in Table 1. Although the datasets for individuals of European and African descent differed in sample size and average age, both sets had a similar proportion of males and females. The datasets were defined by genetic ancestry clusters mapped onto principal components derived from genetic data. We analyzed the European descent and African descent populations separately in all genetic analyses to avoid population stratification confounders. In the EHR-based data without accompanying genetic data, we did not have access to ancestry information. Previous studies have reported differences in calcium and vitamin D levels by race [34,35]. Although race is not an appropriate proxy for genetic ancestry, we stratified our EHR-based phenotypic analyses by race to avoid confounding by race.

3.2. Allele Frequencies for rs1801725 and rs1801726

Consistent with previous reports [25], the A986S variant of the CaSR was more frequent in the European dataset compared to the African dataset. In the European dataset, the minor allele frequency (MAF) for rs1801725 was 0.14 and for rs1801726 was 0.08. In the African dataset, the MAF for rs1801725 was 0.04, while that for rs1801726 was 0.16. The European dataset had more recessive homozygotes for rs1801725 (EUR n = 1005; AFR n = 20) while the African dataset had more recessive homozygotes for rs1801726 (AFR n = 309; EUR n = 87) [28].
We then carried out a preliminary analysis of the association between rs1801725 and the rs1801726 SNPs and disease phenotypes in these datasets by using additive and recessive PheWAS models. Based on a strict Bonferroni correction (p < 3.65 × 10−5), these exploratory analyses revealed that polymorphisms at rs1801725, but not rs1801726, were strongly associated with deficiency of humoral immunity (Figure S1). Although we also observed that these CASR SNPs were associated with several other disease phenotypes, including hypercalcemia and secondary malignancy of respiratory organs, these phenotypes did not pass the strict correction threshold (Figures S2–S4).

3.3. Circulating Calcium, Vitamin D, and PTH Levels in Individuals of European and African Ancestries

In a previous study, total serum calcium was found to be higher in African Americans compared to those in European American subjects [25], but the underlying cause for this difference was not further studied. Here, we determined that both circulating total (p = <0.001) and ionized calcium (p = <0.001) levels are higher in individuals of African descent compared to individuals of European descent (Table 1). Interestingly, in individuals of African descent, serum vitamin D levels are lower (p = <0.001), while intact PTH levels are higher (p = <0.001) than those in individuals of European descent (Table 1). These differences in calciotropic hormone levels suggest that the relatively higher circulating calcium in individuals of African descent may be attributed to reduced vitamin D facilitated reabsorption of calcium in the kidneys and increased PTH mediated bone resorption.

3.4. Association of Circulating Calcium, Vitamin D, and PTH with Primary Cancer Phenotypes

Next, we assessed whether total serum calcium, ionized calcium, vitamin D, and PTH levels are associated with specific primary cancer types (Figure 1, Tables S2–S5). This analysis revealed that circulating calcium, vitamin D, and PTH levels are independently but mostly inversely associated with neoplasms at several pathological sites. However, total serum calcium was positively associated with breast (OR = 1.26, 95% CI = 1.23 − 1.29, p = 3.83 × 10−9), prostate (OR = 1.18, 95% CI = 1.15 − 1.21, p = 1.58 × 10−40), and skin (OR = 1.23, 95% CI = 1.21 − 1.25, p = 4.79 × 10−136) cancers (Figure 1, Table S2). Similarly, ionized calcium was positively associated with skin cancer (OR = 1.07, 95% CI = 1.04 − 1.11, p = 1.47 × 10−5) and liver cancer (OR = 1.19, 95% CI = 1.11 − 1.28, p = 1.55 × 10−6) (Figure 1, Table S3), while circulating vitamin D was not only associated with breast (OR = 1.13, 95% CI = 1.10 − 1.17, p = 9.89 × 10−14) and skin cancers (OR = 1.06, 95% CI = 1.03 − 1.09, p = 2.17 × 10−5), but also thyroid cancer (OR = 1.17, 95% CI = 1.10 − 1.25, p = 1.92 × 10−7) (Figure 1, Table S4). Meanwhile, PTH was positively associated with cancer of urinary organs (OR = 1.25, 95% CI = 1.17 − 1.33, p = 1.16 × 10−11) (Figure 1, Table S5).
We further demonstrate that circulating total and ionized calcium are also significantly, but negatively associated with secondary malignancies of lymph nodes, respiratory organs, digestive system, liver, brain/spine, and bone (Table S6). Serum vitamin D levels are also significantly negatively associated with secondary malignancies of respiratory organs, digestive system, and liver, while circulating PTH levels are significantly associated with secondary malignancy of lymph nodes. Together, this suggests that low circulating calcium and calciotropic hormones are associated with secondary neoplasms at diverse pathological sites.

3.5. Association of CASR SNPS with Circulating Calcium, Vitamin D, and PTH

Several studies have reported the association of rs1801725 and rs1801726 CASR SNPs with calcium levels [25,26,27,28,29,30]. Prior to analyses at specific pathological sites, we assessed the relationship between these SNPs and the CaSR dependent variables, as well as primary cancer and secondary malignancy at all sites as single variables. The significance threshold for this analysis was set at p < 3.38 × 10−4 representing a strict Bonferroni correction (0.05/148) from a total of 148 independent statistical tests. In our European dataset, we confirmed that polymorphisms at rs1801725, but not rs1801726, were strongly associated with total and ionized calcium. The association of rs1801725 with secondary malignancy at all sites with a p-value of 0.03 did not pass the strict Bonferroni correction (Table 2). This analysis also revealed that Vitamin D, PTH, and primary cancers at all sites as a variable were not associated with these CaSR SNPs (Table 2). We also found that median calcium values were significantly (p < 0.05) different between all rs1801725 allele carriers (GG vs. GT vs. TT) and that individuals expressing the TT genotype had the highest median calcium levels compared to those expressing the GG genotype (Figure 2a, Table S7). Although this trend was also observed for ionized calcium, only homozygous dominant and heterozygous carriers (GG vs. GT) and homozygous dominant and homozygous recessive carriers (GG vs. TT) had significantly (p < 0.05) different ionized calcium levels (Figure 2b). The same analytic approach also revealed a similar trend for individuals of African descent regarding total and ionized serum calcium for carriers of polymorphisms at rs1801725 (GG vs. GT) (Figure 2c, Figure 2d, Table S8).
We also confirmed previous reports that the T-allele at rs1801725 is associated with increased total serum calcium and, interestingly, increased ionized calcium levels in the European population (Table S9) [25]. The biggest effect was observed with the assessment of total serum calcium levels with an OR of 1.15 (95% CI = 1.13 − 1.16, p = 8.40 × 10−73), followed by ionized calcium with an OR of 1.10 (95% CI = 1.07 − 1.14, p = 1.43 × 10−9). A similar linear regression analysis revealed that polymorphisms at rs1801726 in individuals of European descent are not associated with circulating calcium, vitamin D, or PTH levels. We next showed that CASR polymorphisms at the rs1801725 locus are associated with circulating total and ionized calcium levels prior to and after cancer diagnosis in individuals of European descent (Table S10).
For individuals of African descent, total serum calcium was modestly (p < 0.05) associated with rs1801725 (OR = 1.12, 95% CI = 1.06 − 1.20, p = 2.34 × 10−4). As in the European population, the rs1801726 locus was not significantly associated with calcium, vitamin D, or PTH (Table S9). Together, this suggests that polymorphisms at rs1801725 but not at rs1801726 are associated with circulating serum calcium in individuals of European or African descent.

3.6. Association of CASR SNPs with Secondary Cancer Types

Although assessment of tumor characteristics and/or progression patterns are not available in the EHRs, we hypothesized that the association of CaSR SNPs with tumors at secondary sites (secondary malignancies) is indicative of the progression of primary tumors to these sites. To test this, we first assessed whether the CaSR SNPs were associated with any secondary malignancy (i.e., cases with any secondary malignancy code). This analysis revealed that the polymorphism at rs1801725 (TT genotype), but not rs1801726, was significantly associated (p < 0.05) with secondary malignancies in the European dataset (Tables S11 and S12).
We next determined which secondary cancers were associated with the rs1801725 CASR SNP. Using an additive logistic regression model, we observed that the rs1801725-T allele was modestly associated with the presence of secondary malignancy of respiratory organs (OR = 1.21, 95% CI = 1.08 − 1.35, p = 9.61 × 10−4) and with secondary malignancy of bone (OR of 1.20 (95% CI = 1.06 − 1.36, p = 4.28 × 10−3) in the European dataset (Figure 3, Table S11). Using a recessive model, we confirmed that the TT genotype at rs1801725 was more common among patients with secondary malignancy of bone (OR = 1.82, 95% CI = 1.25 − 2.65, p = 1.80 × 10−3) (Table S12). These results, while supportive of the hypothesis, did not exceed our strict multiple correction threshold. The small sample size of TT allele carriers in the population of African ancestry for rs1801725 precluded further genotype analysis in this population (Table 1).

4. Discussion

In this study, we assembled larger, genotyped datasets consisting of individuals of European and African descent. We interrogated whether CASR polymorphisms at rs1801725 and rs1801726 are associated with circulating total and ionized calcium, vitamin D, and intact PTH levels. We also investigated whether these SNPs, circulating calcium PTH, and vitamin D are associated with primary cancers at multiple pathological sites as well as with and secondary malignancies for individuals with primary cancer diagnosis. We replicated previous studies showing that circulating calcium levels in individuals of African descent are higher than those of European descent [25]. We further demonstrate that in individuals of European descent, the CASR variant at rs1801725, but not rs1801726, are associated with calcium and bone-related secondary cancer phenotypes, a higher risk of developing secondary neoplastic lesions in the lungs and bone, and interestingly, deficiency of humoral immunity. Overall, our data suggest that dysregulation of calcium homeostasis via expression of potentially inactivating CASR variants at rs1801725 (e.g., A986S) may be an important driver for the progression of primary neoplasms at several pathological sites, including breast, prostate, and skin, as well as metastasis to lung and bone tissues. This study provides further evidence suggesting that dysregulation of calcium homeostasis is distinct in cancer patients of African and European descent. The distinction is more pronounced in patients expressing polymorphic CaSR variants at rs1801725 rather than at rs1801726. Data from this study suggest that the higher circulating calcium in individuals of African descent is presumably due to lower circulating 1,25-OH vitamin D and higher PTH levels. This is especially interesting as reduced vitamin D levels are associated with reduced reabsorption of calcium in the kidneys [36], while higher circulating PTH will lead to increased bone resorption [37] and, consequently, higher circulating calcium levels.
This study used a relatively large dataset of African American patients comprising 10,777 records from this cohort. We were not able to detect any significant associations between polymorphisms at rs1801725 and circulating calcium or cancer phenotypes in this cohort primarily due to the low minor allele frequency (MAF) at this locus. Similar to other studies, the rs1801726 SNP was not associated with any of the cancer phenotypes in either the European or the African ancestry datasets [38]. However, in individuals of African descent, we found modest associations between these CASR SNPs and several disease phenotypes in our exploratory PheWAS. This includes a positive association between rs1801725 and osteoporosis, whereas previous reports found no associations between osteoporosis and rs1801725 in individuals of European ancestry [39]. Finally, the lack of association between rs1801726 and circulating calcium in individuals of African ancestry may suggest a lack of effect of CASR variants at this locus on the activity of the receptor. Although some studies have linked the development of hypercalcemia to CASR variants at these loci [24], direct evidence on the activity of the A986S and the Q1011E receptor variants remains to be clearly demonstrated. Overall, the expression of CASR variants at rs1801725 in >20% of cancer patients may be linked to the reduced sensitivity of the receptor to calcium, and therefore, the potential for the development of CIH. However, whether the difference in circulating calcium levels and /or the expression of CASR variants at whether rs1801725 underlies the disparate aggressiveness of especially prostate and breast cancers in patients of African versus European ancestries remains to be clearly delineated.
As previously reported, the CaSR plays distinct roles in several tissues, including several tumor types. It is a tumor suppressor in certain cancers (e.g., neuroblastoma, gastric, colon, and parathyroid cancers) but a tumor promoter in other cancers such as breast, prostate, ovarian, and kidney cancers [2]. In many of these tissues, CaSR functions by activating several downstream effectors depending on the coupled G proteins [40]. This biased coupling of G proteins to the activated CaSR partly accounts for the diverse role of the receptor in several tissues. Our findings are consistent with the diverse functions of CaSR in that circulating calcium is positively associated with prostate, breast, and skin cancers in which CaSR expression stimulates tumor cell growth and motility [2,3,41] while most other malignancies, including colorectal cancer in which downregulation of CaSR stimulates tumor growth are inversely associated with calcium [42,43]. Although pending further investigation, it is also possible that the expression of distinct variants of the receptor at rs1801725 may lead to distinct CaSR signaling and/or tumor progression patterns in specific subtypes of highly heterogeneous cancers such as breast cancer.
Besides the association with circulating calcium, we also found a significant association between CASR rs1801725 polymorphisms and deficiency in humoral immunity in our exploratory PheWAS analysis. This finding is not only novel but also intriguing. Expression of CaSR in immune cells such as neutrophils, monocytes/macrophages, and T lymphocytes [44,45] is associated with proinflammatory cytokine secretion. Previous studies reported that CaSR regulates immune response by controlling the NLR family pyrin domain containing 3 (NLRP3) in macrophages [46] and in monocytes [47]. The NLRP3 inflammasome is an innate immune mediator responsible for active interleukin (IL)-1β biosynthesis and secretion, which is critical for tumor development and immunity [48]. This is supported by studies showing that activation of NLRP3 inflammasome/IL-1β pathway promotes head and neck squamous cell carcinoma tumorigenesis while inactivation of this pathway delayed tumor growth [49]. Conditional knockout of CaSR in mouse intestinal epithelium also resulted in several defects, including immunity that was skewed towards a pro-inflammatory response [50], consistent with the tumor suppressor role of the receptor in colorectal cancer. While these studies suggest the importance of CaSR in immune cells and innate immunity, others have demonstrated the presence of circulating autoantibodies against the receptor due to either aberrant immune responses against the parathyroid glands [51] or adverse immune events associated with immunosuppressive therapy such as immune checkpoint inhibitors [52,53]. While CaSR autoantibodies may play a limited role, it is possible that the association of CaSR mutants with deficiency in humoral immunity may at least in part be related to diminished activation of NLRP3 inflammasome/IL-1β pathway and anti-tumor immune responses due to decreased sensitivity of the mutant CaSR at rs1801725 to calcium. However, the true implication of exon 7 CaSR variants in humoral immunity and whether this is associated with tumorigenesis or tumor progression requires further investigations.
In addition to its role in calcium homeostasis, vitamin D also plays a central role in innate and adaptive immune function through the vitamin D receptor [54]. Individuals predisposed to higher circulating calcium due to vitamin D deficiency may have a weakened immune system that may further increase susceptibility to the more deadly secondary malignancies. This, coupled with increased PTH and hypercalcemia, may especially compromise high-risk populations such as those of African ancestry who are disproportionately diagnosed with more aggressive tumors and, consequently, higher cancer mortality rates [55].

Study Limitations

In spite of the robustness of our findings and opportunities for future studies, some limitations of this study are worth noting. First, we were unable to analyze PTHrP, the major osteolytic factor during malignancy. This is because PTHrP levels are not routinely measured in a clinical setting, and the small sample size in our EHR datasets warranted their exclusion during the implementation of the QualityLab pipeline. Second, although our dataset containing individuals of African descent was well powered for rs1801725 in the additive models, it was not possible for us to examine the effects of the T-allele in the recessive model because of the low MAF within the African descent population. Third, the skewed nature of the frequency of SNPs at rs1801725 or rs1801726 in individuals of African and European ancestries made it difficult to generalize the effects of the rs1801725 CASR variants in the African population. Fourth, due to differing linkage disequilibrium structures in European and African populations, it is possible that the associations we observed with the European descent do not tag the same causal allele associations in the African descent dataset. Finally, the EHR-based data from which we drew our study population are a clinically ascertained sample; thus, ascertainment biases could impact this study. More functional in vitro and in vivo studies are needed to understand causal links between CaSR genetic variants and cancer progression and/or metastasis.

5. Conclusions

In summary, this study provides robust population-level genetic evidence supporting the notion that expression of the inactivating CaSR rs1801725 SNP predisposes breast, prostate, and skin cancer patients to secondary neoplastic lesions in the lung and bone tissues. Whether the development of secondary malignancies at these sites is exclusively mediated by cancer-induced hypercalcemia via CaSR signaling or includes other complications such as CaSR-mediated stimulation of inflammation and/or tumor immune suppression remains to be fully elucidated.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jpm11070642/s1, Figure S1: Recessive phenome-wide association study of rs1801725 (A986S) and rs1801726 (Q1011E) in the European descent dataset; Figure S2: Additive phenome-wide association study of rs1801725 (A986S) and rs1801726 (Q1011E) in the European descent dataset; Figure S3: Recessive phenome-wide association study of rs1801725 (A986S) and rs1801726 (Q1011E) in the African descent dataset; Figure S4: Additive phenome-wide association study of rs1801725 (A986S) and rs1801726 (Q1011E) in the African descent dataset; Table S1: Description of cancer phecodes used in this study. Table S2: Relationship between circulating total calcium and primary cancer diagnosis at various pathological sites. Table S3: Relationship between circulating ionized calcium and primary cancer diagnosis at various pathological sites. Table S4: Relationship between circulating vitamin D and primary cancer diagnosis at various pathological sites. Table S5: Relationship between circulating intact parathyroid hormone and primary cancer diagnosis at various pathological sites. Table S6: Association of circulating calcium, vitamin D, and parathyroid hormone with secondary malignancies. Table S7: Associating of calcium, vitamin D, and parathyroid hormone levels with polymorphisms at rs1801725 and rs1801726 in individuals of European descent. Table S8: Associating of calcium, vitamin D, and parathyroid hormone levels with polymorphisms at rs1801725 and rs1801726 in individuals of African descent. Table S9: Relationship between CaSR SNPs and circulating calcium, vitamin D, and parathyroid hormone levels. Table S10: Relationship between circulating calcium, vitamin D, and parathyroid hormone levels in CaSR rs1801725 allele carriers before and after cancer diagnosis. Table S11: Association of CaSR SNPs with secondary malignancies in European descent individuals—Additive regression analysis. Table S12: Association of CaSR SNPs with secondary malignancies in European descent individuals—Recessive regression analysis

Author Contributions

Conceptualization: K.V.A., H.K.B., A.B.F., L.K.D., A.M.S.; data curation: K.V.A., A.B.F. formal analysis: K.V.A., A.B.F., L.K.D.; funding acquisition: L.K.D., A.M.S.; investigation: K.V.A., A.B.F.; methodology: K.V.A., A.B.F., L.K.D. project administration: K.V.A., H.K.B.; resources: L.K.D.; supervision: L.K.D., A.M.S.; visualization: K.V.A.; writing—original draft: K.V.A., H.K.B., A.B.F.; writing—reviewing and editing: K.V.A., H.K.B., A.B.F., L.K.D., A.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

KVA is funded by the NIH 5T32GM007628 training grant. HKB is funded by the NIH R25 GM059994 and U54 CA163069 grants. AMS is supported by SC1 CA211030, P50 CA095103, U54 MD007593. LKD is supported by 1R56MH120736. This study used resources from the Advanced Computing Center for Research and Education and the Synthetic Derivative at Vanderbilt University Medical Center, which are supported by Shared Instrumentation Grant S10RR025141 and CTSA grants UL1TR002243, UL1TR000445, and UL1RR024975. Genomic data repositories are supported by investigator-led projects that include U01HG004798, R01NS032830, RC2GM092618, P50GM115305, U01HG006378, U19HL065962, R01HD074711; and additional funding sources listed at https://victr.vumc.org/biovu-funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Vanderbilt University Medical Center (IRB #190418).

Informed Consent Statement

All individuals were de-identified and provided informed consent.

Data Availability Statement

The data that support the findings of this study are available from Vanderbilt University Medical Center, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are, however, available from the authors upon reasonable request and with permission of Vanderbilt University Medical Center.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Aggarwal, A.; Prinz-Wohlgenannt, M.; Tennakoon, S.; Höbaus, J.; Boudot, C.; Mentaverri, R.; Brown, E.M.; Baumgartner-Parzer, S.; Kállay, E. The Calcium-Sensing Receptor: A Promising Target for Prevention of Colorectal Cancer. Biochim. Biophys. Acta (BBA)-Mol. Cell Res. 2015, 1853, 2158–2167. [Google Scholar] [CrossRef] [Green Version]
  2. Tennakoon, S.; Aggarwal, A.; Kállay, E. The Calcium-Sensing Receptor and the Hallmarks of Cancer. Biochim. Biophys. Acta (BBA)-Mol. Cell Res. 2016, 1863, 1398–1407. [Google Scholar] [CrossRef] [Green Version]
  3. Kim, W.; Wysolmerski, J.J. Calcium-Sensing Receptor in Breast Physiology and Cancer. Front. Physiol. 2016, 7, 440. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Masvidal, L.; Iniesta, R.; García, M.; Casalà, C.; Lavarino, C.; Mora, J.; de Torres, C. Genetic Variants in the Promoter Region of the Calcium-Sensing Receptor Gene Are Associated with Its down-Regulation in Neuroblastic Tumors. Mol. Carcinog. 2017, 56, 1281–1289. [Google Scholar] [CrossRef] [PubMed]
  5. Hernández-Bedolla, M.A.; Carretero-Ortega, J.; Valadez-Sánchez, M.; Vázquez-Prado, J.; Reyes-Cruz, G. Chemotactic and Proangiogenic Role of Calcium Sensing Receptor Is Linked to Secretion of Multiple Cytokines and Growth Factors in Breast Cancer MDA-MB-231 Cells. Biochim. Biophys. Acta (BBA)-Mol. Cell Res. 2015, 1853, 166–182. [Google Scholar] [CrossRef] [Green Version]
  6. Liu, G.; Hu, X.; Chakrabarty, S. Calcium Sensing Receptor Down-Regulates Malignant Cell Behavior and Promotes Chemosensitivity in Human Breast Cancer Cells. Cell Calcium 2009, 45, 216–225. [Google Scholar] [CrossRef] [PubMed]
  7. Singh, N.; Aslam, M.N.; Varani, J.; Chakrabarty, S. Induction of Calcium Sensing Receptor in Human Colon Cancer Cells by Calcium, Vitamin D and Aquamin: Promotion of a More Differentiated, Less Malignant and Indolent Phenotype. Mol. Carcinog. 2015, 54, 543–553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Tfelt-Hansen, J.; Brown, E.M. The Calcium-Sensing Receptor in Normal Physiology and Pathophysiology: A Review. Crit Rev. Clin. Lab. Sci. 2005, 42, 35–70. [Google Scholar] [CrossRef] [PubMed]
  9. Brown, E.M.; Gamba, G.; Riccardi, D.; Lombardi, M.; Butters, R.; Kifor, O.; Sun, A.; Hediger, M.A.; Lytton, J.; Hebert, S.C. Cloning and Characterization of an Extracellular Ca(2+) Sensing Receptor from Bovine Parathyroid. Nature 1993, 366, 575–580. [Google Scholar] [CrossRef]
  10. Ming, J.; Cronin, S.J.F.; Penninger, J.M. Targeting the RANKL/RANK/OPG Axis for Cancer Therapy. Front. Oncol. 2020, 10, 1283. [Google Scholar] [CrossRef]
  11. Boyce, B.F.; Xing, L. Functions of RANKL/RANK/OPG in Bone Modeling and Remodeling. Arch. Biochem. Biophys. 2008, 473, 139–146. [Google Scholar] [CrossRef] [Green Version]
  12. Guise, T.A.; Kozlow, W.M.; Heras-Herzig, A.; Padalecki, S.S.; Yin, J.J.; Chirgwin, J.M. Molecular Mechanisms of Breast Cancer Metastases to Bone. Clin. Breast Cancer 2005, 5, S46–S53. [Google Scholar] [CrossRef]
  13. Mundy, G.R. Metastasis to Bone: Causes, Consequences and Therapeutic Opportunities. Nat. Rev. Cancer 2002, 2, 584–593. [Google Scholar] [CrossRef] [PubMed]
  14. Sanders, J.L.; Chattopadhyay, N.; Kifor, O.; Yamaguchi, T.; Butters, R.R.; Brown, E.M. Extracellular Calcium-Sensing Receptor Expression and Its Potential Role in Regulating Parathyroid Hormone-Related Peptide Secretion in Human Breast Cancer Cell Lines. Endocrinology 2000, 141, 4357–4364. [Google Scholar] [CrossRef]
  15. Ralston, S.H.; Gallacher, S.J.; Patel, U.; Campbell, J.; Boyle, I.T. Cancer-Associated Hypercalcemia: Morbidity and Mortality. Clinical Experience in 126 Treated Patients. Ann. Intern. Med. 1990, 112, 499–504. [Google Scholar] [CrossRef]
  16. Stewart, A.F. Hypercalcemia Associated with Cancer. N. Engl. J. Med. 2005, 352, 373–379. [Google Scholar] [CrossRef]
  17. Almquist, M.; Anagnostaki, L.; Bondeson, L.; Bondeson, A.-G.; Borgquist, S.; Landberg, G.; Malina, J.; Malm, J.; Manjer, J. Serum Calcium and Tumour Aggressiveness in Breast Cancer: A Prospective Study of 7847 Women. Eur. J. Cancer Prev. 2009, 18, 354–360. [Google Scholar] [CrossRef] [PubMed]
  18. Thaw, S.S.H.; Sahmoun, A.; Schwartz, G.G. Serum Calcium, Tumor Size, and Hormone Receptor Status in Women with Untreated Breast Cancer. Cancer Biol. Ther. 2012, 13, 467–471. [Google Scholar] [CrossRef]
  19. Hendy, G.N.; D’Souza-Li, L.; Yang, B.; Canaff, L.; Cole, D.E. Mutations of the Calcium-Sensing Receptor (CASR) in Familial Hypocalciuric Hypercalcemia, Neonatal Severe Hyperparathyroidism, and Autosomal Dominant Hypocalcemia. Hum. Mutat. 2000, 16, 281–296. [Google Scholar] [CrossRef]
  20. Pidasheva, S.; D’Souza-Li, L.; Canaff, L.; Cole, D.E.C.; Hendy, G.N. CASRdb: Calcium-Sensing Receptor Locus-Specific Database for Mutations Causing Familial (Benign) Hypocalciuric Hypercalcemia, Neonatal Severe Hyperparathyroidism, and Autosomal Dominant Hypocalcemia. Hum. Mutat. 2004, 24, 107–111. [Google Scholar] [CrossRef] [PubMed]
  21. Thakker, R.V. Diseases Associated with the Extracellular Calcium-Sensing Receptor. Cell Calcium 2004, 35, 275–282. [Google Scholar] [CrossRef]
  22. Pratt, J.H.; Ambrosius, W.T.; Wagner, M.A.; Maharry, K. Molecular Variations in the Calcium-Sensing Receptor in Relation to Sodium Balance and Presence of Hypertension in Blacks and Whites. Am. J. Hypertens. 2000, 13, 654–658. [Google Scholar] [CrossRef] [Green Version]
  23. Muddana, V.; Lamb, J.; Greer, J.B.; Elinoff, B.; Hawes, R.H.; Cotton, P.B.; Anderson, M.A.; Brand, R.E.; Slivka, A.; Whitcomb, D.C. Association between Calcium Sensing Receptor Gene Polymorphisms and Chronic Pancreatitis in a US Population: Role of Serine Protease Inhibitor Kazal 1type and Alcohol. World J. Gastroenterol. 2008, 14, 4486. [Google Scholar] [CrossRef]
  24. Lorch, G.; Viatchenko-Karpinski, S.; Ho, H.-T.; Dirksen, W.P.; Toribio, R.E.; Foley, J.; Györke, S.; Rosol, T.J. The Calcium-Sensing Receptor Is Necessary for the Rapid Development of Hypercalcemia in Human Lung Squamous Cell Carcinoma. Neoplasia 2011, 13, 428–438. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wang, L.; Widatalla, S.E.; Whalen, D.S.; Ochieng, J.; Sakwe, A.M. Association of Calcium Sensing Receptor Polymorphisms at Rs1801725 with Circulating Calcium in Breast Cancer Patients. BMC Cancer 2017, 17, 511. [Google Scholar] [CrossRef] [PubMed]
  26. Schwartz, G.G.; John, E.M.; Rowland, G.; Ingles, S.A. Prostate Cancer in African-American Men and Polymorphism in the Calcium-Sensing Receptor. Cancer Biol. Ther. 2010, 9, 994–999. [Google Scholar] [CrossRef] [Green Version]
  27. März, W.; Seelhorst, U.; Wellnitz, B.; Tiran, B.; Obermayer-Pietsch, B.; Renner, W.; Boehm, B.O.; Ritz, E.; Hoffmann, M.M. Alanine to Serine Polymorphism at Position 986 of the Calcium-Sensing Receptor Associated with Coronary Heart Disease, Myocardial Infarction, All-Cause, and Cardiovascular Mortality. J. Clin. Endocrinol. Metab. 2007, 92, 2363–2369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Scillitani, A.; Guarnieri, V.; De Geronimo, S.; Muscarella, L.A.; Battista, C.; D’Agruma, L.; Bertoldo, F.; Florio, C.; Minisola, S.; Hendy, G.N.; et al. Blood Ionized Calcium Is Associated with Clustered Polymorphisms in the Carboxyl-Terminal Tail of the Calcium-Sensing Receptor. J. Clin. Endocrinol. Metab. 2004, 89, 5634–5638. [Google Scholar] [CrossRef] [Green Version]
  29. Cole, D.E.; Vieth, R.; Trang, H.M.; Wong, B.Y.; Hendy, G.N.; Rubin, L.A. Association between Total Serum Calcium and the A986S Polymorphism of the Calcium-Sensing Receptor Gene. Mol. Genet. Metab 2001, 72, 168–174. [Google Scholar] [CrossRef]
  30. Vezzoli, G.; Scillitani, A.; Corbetta, S.; Terranegra, A.; Dogliotti, E.; Guarnieri, V.; Arcidiacono, T.; Macrina, L.; Mingione, A.; Brasacchio, C.; et al. Risk of Nephrolithiasis in Primary Hyperparathyroidism Is Associated with Two Polymorphisms of the Calcium-Sensing Receptor Gene. J. Nephrol. 2015, 28, 67–72. [Google Scholar] [CrossRef] [PubMed]
  31. Roden, D.M.; Pulley, J.M.; Basford, M.A.; Bernard, G.R.; Clayton, E.W.; Balser, J.R.; Masys, D.R. Development of a Large-Scale de-Identified DNA Biobank to Enable Personalized Medicine. Clin. Pharmacol. Ther. 2008, 84, 362–369. [Google Scholar] [CrossRef] [Green Version]
  32. Denny, J.C.; Bastarache, L.; Ritchie, M.D.; Carroll, R.J.; Zink, R.; Mosley, J.D.; Field, J.R.; Pulley, J.M.; Ramirez, A.H.; Bowton, E.; et al. Systematic Comparison of Phenome-Wide Association Study of Electronic Medical Record Data and Genome-Wide Association Study Data. Nat. Biotechnol. 2013, 31, 1102–1110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Dennis, J.K.; Sealock, J.M.; Straub, P.; Lee, Y.H.; Hucks, D.; Actkins, K.; Faucon, A.; Feng, Y.-C.A.; Ge, T.; Goleva, S.B.; et al. Clinical Laboratory Test-Wide Association Scan of Polygenic Scores Identifies Biomarkers of Complex Disease. Genome Med. 2021, 13, 6. [Google Scholar] [CrossRef] [PubMed]
  34. Bien, S.A.; Wojcik, G.L.; Zubair, N.; Gignoux, C.R.; Martin, A.R.; Kocarnik, J.M.; Martin, L.W.; Buyske, S.; Haessler, J.; Walker, R.W.; et al. Strategies for Enriching Variant Coverage in Candidate Disease Loci on a Multiethnic Genotyping Array. PLoS ONE 2016, 11, e0167758. [Google Scholar] [CrossRef] [PubMed]
  35. Abraham, G.; Qiu, Y.; Inouye, M. FlashPCA2: Principal Component Analysis of Biobank-Scale Genotype Datasets. Bioinformatics 2017, 33, 2776–2778. [Google Scholar] [CrossRef]
  36. Kumar, R.; Tebben, P.J.; Thompson, J.R. Vitamin D and the Kidney. Arch. Biochem. Biophys. 2012, 523, 77–86. [Google Scholar] [CrossRef] [Green Version]
  37. Silva, B.C.; Bilezikian, J.P. Parathyroid Hormone: Anabolic and Catabolic Actions on the Skeleton. Curr. Opin. Pharmacol. 2015, 22, 41–50. [Google Scholar] [CrossRef] [Green Version]
  38. Huang, H.; Li, T.; Liao, D.; Zhu, Z.; Dong, Y. Quantitative Assessment of the Clinical Susceptibility of Calcium-Sensing Receptor Polymorphisms in Cancer Patients. Cancer Manag. Res. 2018, 10, 755–763. [Google Scholar] [CrossRef] [Green Version]
  39. Kapur, K.; Johnson, T.; Beckmann, N.D.; Sehmi, J.; Tanaka, T.; Kutalik, Z.; Styrkarsdottir, U.; Zhang, W.; Marek, D.; Gudbjartsson, D.F.; et al. Genome-Wide Meta-Analysis for Serum Calcium Identifies Significantly Associated SNPs near the Calcium-Sensing Receptor (CASR) Gene. PLoS Genet. 2010, 6, e1001035. [Google Scholar] [CrossRef] [Green Version]
  40. Abid, H.A.; Inoue, A.; Gorvin, C.M. Heterogeneity of G-Protein Activation by the Calcium-Sensing Receptor. J. Mol. Endocrinol. 2021. [Google Scholar] [CrossRef]
  41. Saidak, Z.; Boudot, C.; Abdoune, R.; Petit, L.; Brazier, M.; Mentaverri, R.; Kamel, S. Extracellular Calcium Promotes the Migration of Breast Cancer Cells through the Activation of the Calcium Sensing Receptor. Exp. Cell Res. 2009, 315, 2072–2080. [Google Scholar] [CrossRef]
  42. Chakrabarty, S.; Wang, H.; Canaff, L.; Hendy, G.N.; Appelman, H.; Varani, J. Calcium Sensing Receptor in Human Colon Carcinoma: Interaction with Ca(2+) and 1,25-Dihydroxyvitamin D(3). Cancer Res. 2005, 65, 493–498. [Google Scholar] [PubMed]
  43. Whitfield, J.F. Calcium, Calcium-Sensing Receptor and Colon Cancer. Cancer Lett. 2009, 275, 9–16. [Google Scholar] [CrossRef]
  44. Iamartino, L.; Elajnaf, T.; Kallay, E.; Schepelmann, M. Calcium-Sensing Receptor in Colorectal Inflammation and Cancer: Current Insights and Future Perspectives. World J. Gastroenterol. 2018, 24, 4119–4131. [Google Scholar] [CrossRef]
  45. Li, F.; Tian, Z. The Liver Works as a School to Educate Regulatory Immune Cells. Cell Mol. Immunol. 2013, 10, 292–302. [Google Scholar] [CrossRef] [Green Version]
  46. Lee, G.-S.; Subramanian, N.; Kim, A.I.; Aksentijevich, I.; Goldbach-Mansky, R.; Sacks, D.B.; Germain, R.N.; Kastner, D.L.; Chae, J.J. The Calcium-Sensing Receptor Regulates the NLRP3 Inflammasome through Ca2+ and CAMP. Nature 2012, 492, 123–127. [Google Scholar] [CrossRef] [Green Version]
  47. Rossol, M.; Pierer, M.; Raulien, N.; Quandt, D.; Meusch, U.; Rothe, K.; Schubert, K.; Schöneberg, T.; Schaefer, M.; Krügel, U.; et al. Extracellular Ca2+ Is a Danger Signal Activating the NLRP3 Inflammasome through G Protein-Coupled Calcium Sensing Receptors. Nat. Commun. 2012, 3, 1329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Moossavi, M.; Parsamanesh, N.; Bahrami, A.; Atkin, S.L.; Sahebkar, A. Role of the NLRP3 Inflammasome in Cancer. Mol. Cancer 2018, 17, 158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Chen, L.; Huang, C.-F.; Li, Y.-C.; Deng, W.-W.; Mao, L.; Wu, L.; Zhang, W.-F.; Zhang, L.; Sun, Z.-J. Blockage of the NLRP3 Inflammasome by MCC950 Improves Anti-Tumor Immune Responses in Head and Neck Squamous Cell Carcinoma. Cell Mol. Life Sci. 2018, 75, 2045–2058. [Google Scholar] [CrossRef]
  50. Cheng, S.X.; Lightfoot, Y.L.; Yang, T.; Zadeh, M.; Tang, L.; Sahay, B.; Wang, G.P.; Owen, J.L.; Mohamadzadeh, M. Epithelial CaSR Deficiency Alters Intestinal Integrity and Promotes Proinflammatory Immune Responses. FEBS Lett. 2014, 588, 4158–4166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Piranavan, P.; Li, Y.; Brown, E.; Kemp, E.H.; Trivedi, N. Immune Checkpoint Inhibitor-Induced Hypoparathyroidism Associated With Calcium-Sensing Receptor-Activating Autoantibodies. J. Clin. Endocrinol. Metab. 2019, 104, 550–556. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Lupi, I.; Brancatella, A.; Cetani, F.; Latrofa, F.; Kemp, E.H.; Marcocci, C. Activating Antibodies to The Calcium-Sensing Receptor in Immunotherapy-Induced Hypoparathyroidism. J. Clin. Endocrinol. Metab. 2020, 105, 1581–1588. [Google Scholar] [CrossRef] [PubMed]
  53. Chamberlin, M.; Kemp, E.H.; Weetman, A.P.; Khadka, B.; Brown, E.M. Immunosuppressive Therapy of Autoimmune Hypoparathyroidism in a Patient with Activating Autoantibodies against the Calcium-Sensing Receptor. Clin. Endocrinol. 2019, 90, 214–221. [Google Scholar] [CrossRef] [Green Version]
  54. Aranow, C. Vitamin D and the Immune System. J. Investig. Med. 2011, 59, 881–886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Zagzag, J.; Hu, M.I.; Fisher, S.B.; Perrier, N.D. Hypercalcemia and Cancer: Differential Diagnosis and Treatment. CA Cancer J. Clin. 2018, 68, 377–386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Association of circulating calcium, parathyroid hormone, and vitamin D levels with cancer risk. Forest plots showing the odds associated with circulating levels of the indicated laboratory values and diagnosis of cancer at the indicated pathologic sites. The dashed line in each plot represents an odds ratio (OR) of 1. PTH = intact parathyroid hormone, VITD = vitamin D, IONCAL = ionized calcium.
Figure 1. Association of circulating calcium, parathyroid hormone, and vitamin D levels with cancer risk. Forest plots showing the odds associated with circulating levels of the indicated laboratory values and diagnosis of cancer at the indicated pathologic sites. The dashed line in each plot represents an odds ratio (OR) of 1. PTH = intact parathyroid hormone, VITD = vitamin D, IONCAL = ionized calcium.
Jpm 11 00642 g001
Figure 2. Circulating total calcium and ionized calcium in CaSR rs1801725 allele carriers. Top row (a,b): European dataset; bottom row (c,d): African dataset. The Wilcoxon rank sum test was used to compare the median circulating total (a) and ionized (b) calcium levels in individuals of European descent, and total (c) and ionized (d) calcium levels in individuals of African descent expressing the indicated polymorphisms at rs1801725 depicted by the genotype. p < 0.05 was considered statistically significant.
Figure 2. Circulating total calcium and ionized calcium in CaSR rs1801725 allele carriers. Top row (a,b): European dataset; bottom row (c,d): African dataset. The Wilcoxon rank sum test was used to compare the median circulating total (a) and ionized (b) calcium levels in individuals of European descent, and total (c) and ionized (d) calcium levels in individuals of African descent expressing the indicated polymorphisms at rs1801725 depicted by the genotype. p < 0.05 was considered statistically significant.
Jpm 11 00642 g002
Figure 3. Association of circulating calcium levels with secondary cancer phenotypes. Logistic regression models were performed on secondary cancer and the CASR variants at rs1801725. The additive (red) and recessive (green) models were adjusted for median age, sex, and the first ten principal components. * indicates p < 0.01 in the logistic regression model. OR = odds ratio.
Figure 3. Association of circulating calcium levels with secondary cancer phenotypes. Logistic regression models were performed on secondary cancer and the CASR variants at rs1801725. The additive (red) and recessive (green) models were adjusted for median age, sex, and the first ten principal components. * indicates p < 0.01 in the logistic regression model. OR = odds ratio.
Jpm 11 00642 g003
Table 1. Descriptive statistics for European and African descent datasets.
Table 1. Descriptive statistics for European and African descent datasets.
DemographicsEuropean DatasetAfrican Datasetp
N5368010777
Average Age (SD)48.74 (22.03)37.98 (21.57)<0.001
Average BMI (SD)27.75 (7.25)29.04 (8.65)<0.001
Sex (%) <0.001
Female30801 (57.4)6767 (62.8)
Male22879 (42.6)4010 (37.2)
rs1801725 Alleles (%) <0.001
GG39939 (74.4)9976 (92.6)
GT12736 (23.7) 773 (7.2)
TT1005 (1.9)20 (0.2)
rs1801726 Alleles (%) <0.001
CC49733 (92.7)7487 (69.5)
CG3858 (7.2)2972 (27.6)
GG87 (0.2)309 (2.9)
EHR-Reported Race (%) <0.001
White52461 (97.7)285 (2.6)
Black59 (0.1)10313 (95.7)
Calcium (mg/dL)
N511769878
Median [IQR]9.30 (9.05, 9.60)9.40 (9.10, 9.60)<0.001
Ionized Calcium (mg/dL)
N138992336
Median [IQR]4.52 (4.32, 4.73)4.57 (4.35, 4.80)<0.001
Vitamin D (ng/mL)
N188853341
Median [IQR]31 (24, 39)23 (17, 30.50)<0.001
Parathyroid Hormone (pg/mL)
N69371430
Median [IQR]56 (36, 91)82 (50, 160)<0.001
Categorical data and allele information were statistically analyzed by chi-square tests, lab data were analyzed by Wilcoxon rank sum tests, and age and BMI were analyzed by Student’s t-test. SD = standard deviation, EHR = electronic health record, IQR = interquartile range.
Table 2. Relationship between CASR exon 7 variants and calcium, PTH, vitamin D, and cancer in European descent individuals.
Table 2. Relationship between CASR exon 7 variants and calcium, PTH, vitamin D, and cancer in European descent individuals.
rs1801725rs1801726
OR (95% CI)pOR (95% CI)p
Calcium1.15 (1.13–1.16)8.40 × 10−73*0.98 (0.96–1.01)0.16
Ionized Calcium1.10 (1.07–1.14)1.43 × 10−9*1.03 (0.97–1.09)0.28
Vitamin D0.99 (0.96–1.02)0.500.98 (0.94–1.03)0.51
Parathyroid Hormone1.00 (0.95–1.05)0.981.03 (0.94–1.12)0.57
Primary Cancer1.02 (0.97–1.08)0.401.00 (0.91–1.09)0.95
Secondary Malignancy1.08 (1.01–1.15)0.031.00 (0.88–1.14)0.97
The additive regression model adjusted for median age, sex, and the first ten principal components was performed between CASR SNPs (rs180125 (G vs. T) and rs1801726 (C vs. G)) and the indicated laboratory values as well as cancer diagnosis. In this analysis, the phecodes used for the primary cancers and secondary malignancies were each collapsed into a single diagnostic code. * p < 3.38 × 10−4 was considered statistically 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

Actkins, K.V.; Beasley, H.K.; Faucon, A.B.; Davis, L.K.; Sakwe, A.M. Calcium-Sensing Receptor Polymorphisms at rs1801725 Are Associated with Increased Risk of Secondary Malignancies. J. Pers. Med. 2021, 11, 642. https://doi.org/10.3390/jpm11070642

AMA Style

Actkins KV, Beasley HK, Faucon AB, Davis LK, Sakwe AM. Calcium-Sensing Receptor Polymorphisms at rs1801725 Are Associated with Increased Risk of Secondary Malignancies. Journal of Personalized Medicine. 2021; 11(7):642. https://doi.org/10.3390/jpm11070642

Chicago/Turabian Style

Actkins, Ky’Era V., Heather K. Beasley, Annika B. Faucon, Lea K. Davis, and Amos M. Sakwe. 2021. "Calcium-Sensing Receptor Polymorphisms at rs1801725 Are Associated with Increased Risk of Secondary Malignancies" Journal of Personalized Medicine 11, no. 7: 642. https://doi.org/10.3390/jpm11070642

APA Style

Actkins, K. V., Beasley, H. K., Faucon, A. B., Davis, L. K., & Sakwe, A. M. (2021). Calcium-Sensing Receptor Polymorphisms at rs1801725 Are Associated with Increased Risk of Secondary Malignancies. Journal of Personalized Medicine, 11(7), 642. https://doi.org/10.3390/jpm11070642

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