Influence of single nucleotide polymorphisms in the MMP1 promoter region on cutaneous melanoma progression

The relationship between single nucleotide polymorphisms and skin cancer susceptibility: A systematic review and network meta-analysis

Background

Single nucleotide polymorphisms (SNPs) interfere with the function of certain genes and thus may influence the probability of skin cancer. The correlation between SNPs and skin cancer (SC) lacks statistical power, however. Therefore, the purpose of this study was to identify the gene polymorphisms involved in skin cancer susceptibility using network meta-analysis and to determine the relationship between SNPs and SC risk.

Methods

PubMed, Embase, and Web of Science were searched for articles including "SNP" and different types of SC as keywords between January 2005 and May 2022. The Newcastle-Ottawa Scale was used to assess bias judgments. The odds ratio (ORs) and their 95% confidence intervals (CIs) were determined to estimate heterogeneity within and between studies. Meta-analysis and network meta-analysis were carried out to identify the SNPs associated with SC. The P-score of each SNP was compared to obtain the rank of probability. Subgroup analyses were performed by cancer type.

Results

A total of 275 SNPs from 59 studies were included in the study. Two subgroup SNP networks using the allele model and dominant model were analyzed. The alternative alleles of rs2228570 (FokI) and rs13181 (ERCC2) were the first-ranking SNPs in both subgroups one and two of the allele model, respectively. The homozygous dominant genotype and heterozygous genotype of rs475007 in subgroup one and the homozygous recessive genotype of rs238406 in subgroup two were most likely to be associated with skin cancer based on the dominant model.

Conclusions

According to the allele model, SNPs FokI rs2228570 and ERCC2 rs13181 and, according to the dominant model, SNPs MMP1 rs475007 and ERCC2 rs238406 are closely linked to SC risk.

Ulcerated Cutaneous Melanoma: A Review of the Clinical, Histologic, and Molecular Features Associated with a Clinically Aggressive Histologic Phenotype

The presence of ulceration in melanoma is associated with poor clinical outcomes and is the third most powerful predictor of survival in the AJCC Melanoma Staging System after tumor thickness and mitotic activity. The aggressive biological behavior associated with ulceration has been hypothesized to be the result of an intrinsic biological attribute that favors dissemination and presents locally with the loss of epidermal integrity. Among the features of ulcerated melanoma, many show promise as potential prognostic tools, markers of differential immunogenicity and indicators of oncogenic drivers of invasion and metastasis. The incidence of ulcerated melanoma is greater in males, increases with age and with systemic inflammatory risk factors (diabetes, smoking, low vitamin D, elevated body mass index). Patients with ulcerated primary tumors seem to exclusively benefit from adjuvant interferon (IFN) therapy, which is likely the consequence of an altered tumor microenvironment. When ulceration is present, there is a higher density of macrophages and dendritic cells and enhanced expression of pro-inflammatory cytokines, such as IL-6. There is also an increased expression of proteins involved in tumor antigen presentation in ulcerated melanomas. Histologically, vascular density, vasculogenic mimicry and angiotropism are all significantly correlated with ulceration in melanoma. The presence of ulceration is associated with reduced protein expression of E-cadherin and PTEN and elevated levels of N-cadherin and the matrix metalloproteinases. Differential microRNA expression also holds promise as a potential prognostic biomarker of malignancy and disease spread within the setting of ulceration. However, the molecular and cellular differences associated with the ulcerated state are complex and further study will aid in determining how these differences can be harnessed to improve care for patients with melanoma.

Association of Matrix Metalloproteinase (MMP) Gene Polymorphisms With Knee Osteoarthritis: A Review of the Literature

Progressive matrix metalloproteinase (MMP)-induced degradation of the extracellular matrix (ECM) of the articular cartilage is one of the major pathogenic osteoarthritis (OA) events. Several single nucleotide polymorphisms (SNPs) in genes encoding MMPs have been identified as affecting MMP expression, production, and enzymatic activity. This study systematically reviews the literature regarding the association between the SNPs of genes encoding MMPs and the risk of knee OA. An electronic search in the PubMed and Web of Science databases from conception to January 2021 was performed addressing studies relating MMPs genetic polymorphisms with the risk of knee OA. We included case-control studies that used validated genotyping methods to detect the SNPs’ association in MMP genes with primary knee OA risk. Ten studies were finally included in this systematic review, evaluating different SNPs in six MMP genes in terms of knee OA pathogenesis: MMP-1 (3 SNPs), MMP-2 (1 SNP), MMP-3 (9 SNPs), MMP-8 (10 SNPs), MMP-9 (6 SNPs), and MMP-13 (1 SNP). Among them, nine SNPs of four MMP genes have been associated with knee OA: (a) MMP-1 -1607 1G/2G (Turkish, Chinese), (b) MMP-3 rs650108, rs650108, rs520540, rs602128, rs679620 (Chinese), (c) MMP-8 rs1940475 and rs376520 (Finnish), and (d) MMP-13 77A/ (rs2252070) (Chinese). The present review summarizes all known SNPs of MMP genes related to a higher risk of knee OA. There are at least nine SNPs in four MMP genes associated with knee OA. No solid correlation between MMP genotype and knee OA phenotype exists. More high-quality studies and modern genetic testing methods are needed to fully elucidate the role of polymorphisms of MMP genes in knee OA pathogenesis.

Germline variants are associated with increased primary melanoma tumor thickness at diagnosis

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma.

Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness.

To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness.

The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry.

We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10⁻⁸) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

Impact of phosphoinositide-3-kinase and vitamin D3 nuclear receptor single-nucleotide polymorphisms on the outcome of malignant melanoma patients

BACKGROUND

Several studies associating single nucleotide polymorphisms (SNPs) frequencies with tumors outcome have been conducted, nevertheless malignant melanoma literature data are inconclusive.Therefore we evaluate the impact of different genotypes for phosphoinositide-3-kinase (PI3K) and vitamin D3 nuclear receptor (VDR) SNPs on melanoma patients' outcome.

MATERIALS AND METHODS

Genomic DNA of 88 patients was extracted from blood and tumor samples. SNPs were determined by PCR using TaqMan assays. We selected polymorphisms of the regulatory and catalytic subunit of PI3K (PIK3R1 and PIK3CA genes, respectively), analyzing rs2699887C>T of PIK3CA and rs3730089G>A of PIK3R1 SNPs. Furthermore we considered the following VDR SNPs: rs2228570A>G (Fok1), rs731236A>G (Taq1) and rs1544410C>T (Bsm1).Progression free survival (PFS) and overall survival (OS) were estimated with the Kaplan-Meier method and with Mantel-Haenszel log-rank test.

RESULTS

The statistical analysis for Fok1 of VDR showed a significant difference in PFS after the first line therapy (median PFS= 21.2 months in the homozygous recessive genotype group vs. 3.3 months of homozygous dominant and heterozygous ones, p= 0.03). In particular, in homozygous recessive patients for Fok1 SNPs of VDR a high rate of histological regression and BRAF (B- Rapidly Accelerated Fibrosarcoma gene) mutation were observed. Furthermore, more efficacy of BRAF +/- MEK (MAPK-ERK-Kinase) inhibitors therapies in homozygous recessive patients vs. homozygous dominant and heterozygous ones was shown.

CONCLUSIONS

Our study showed a significant correlation between homozygous recessive genotype of Fok1 SNPs of VDR gene and an increased PFS in patients who underwent a first line therapy with BRAF inhibitors.

Component-wise gradient boosting and false discovery control in survival analysis with high-dimensional covariates

MOTIVATION

Technological advances that allow routine identification of high-dimensional risk factors have led to high demand for statistical techniques that enable full utilization of these rich sources of information for genetics studies. Variable selection for censored outcome data as well as control of false discoveries (i.e. inclusion of irrelevant variables) in the presence of high-dimensional predictors present serious challenges. This article develops a computationally feasible method based on boosting and stability selection. Specifically, we modified the component-wise gradient boosting to improve the computational feasibility and introduced random permutation in stability selection for controlling false discoveries.

RESULTS

We have proposed a high-dimensional variable selection method by incorporating stability selection to control false discovery. Comparisons between the proposed method and the commonly used univariate and Lasso approaches for variable selection reveal that the proposed method yields fewer false discoveries. The proposed method is applied to study the associations of 2339 common single-nucleotide polymorphisms (SNPs) with overall survival among cutaneous melanoma (CM) patients. The results have confirmed that BRCA2 pathway SNPs are likely to be associated with overall survival, as reported by previous literature. Moreover, we have identified several new Fanconi anemia (FA) pathway SNPs that are likely to modulate survival of CM patients.

AVAILABILITY AND IMPLEMENTATION

The related source code and documents are freely available at https://sites.google.com/site/bestumich/issues.

CONTACT

yili@umich.edu.

Germline determinants of clinical outcome of cutaneous melanoma

Cutaneous melanoma (CM) is the most lethal form of skin cancer. Despite the constant increase in melanoma incidence, which is in part due to incremental advances in early diagnostic modalities, mortality rates have not improved over the last decade and for advanced stages remain steadily high. While conventional prognostic biomarkers currently in use find significant utility for predicting overall general survival probabilities, they are not sensitive enough for a more personalized clinical assessment on an individual level. In recent years, the advent of genomic technologies has brought the promise of identification of germline DNA alterations that may associate with CM outcomes and hence represent novel biomarkers for clinical utilization. This review attempts to summarize the current state of knowledge of germline genetic factors studied for their impact on melanoma clinical outcomes. We also discuss ongoing problems and hurdles in validating such surrogates, and we also project future directions in discovery of more powerful germline genetic factors with clinical utility in melanoma prognostication.

Vitamin D receptor polymorphisms and survival in patients with cutaneous melanoma: a population-based study

Factors known to affect melanoma survival include age at presentation, sex and tumor characteristics. Polymorphisms also appear to modulate survival following diagnosis. Result from other studies suggest that vitamin D receptor (VDR) polymorphisms (SNPs) impact survival in patients with glioma, renal cell carcinoma, lung, breast, prostate and other cancers; however, a comprehensive study of VDR polymorphisms and melanoma-specific survival is lacking. We aimed to investigate whether VDR genetic variation influences survival in patients with cutaneous melanoma. The analysis involved 3566 incident single and multiple primary melanoma cases enrolled in the international population-based Genes, Environment, and Melanoma Study. Melanoma-specific survival outcomes were calculated for each of 38 VDR SNPs using a competing risk analysis after adjustment for covariates. There were 254 (7.1%) deaths due to melanoma during the median 7.6 years follow-up period. VDR SNPs rs7299460, rs3782905, rs2239182, rs12370156, rs2238140, rs7305032, rs1544410 (BsmI) and rs731236 (TaqI) each had a statistically significant (trend P values < 0.05) association with melanoma-specific survival in multivariate analysis. One functional SNP (rs2239182) remained significant after adjustment for multiple testing using the Monte Carlo method. None of the SNPs associated with survival were significantly associated with Breslow thickness, ulceration or mitosis. These results suggest that the VDR gene may influence survival from melanoma, although the mechanism by which VDR exerts its effect does not seem driven by tumor aggressiveness. Further investigations are needed to confirm our results and to understand the relationship between VDR and survival in the combined context of tumor and host characteristics.

Genetic variants in fanconi anemia pathway genes BRCA2 and FANCA predict melanoma survival

Cutaneous melanoma (CM) is the most lethal skin cancer. The Fanconi Anemia (FA) pathway involved in DNA crosslinks repair may affect CM susceptibility and prognosis. Using data derived from published genome-wide association study, we comprehensively analyzed the associations of 2339 common single nucleotide polymorphisms (SNPs) in 14 autosomal FA genes with overall survival (OS) in 858 CM patients. By performing false-positive report probability corrections and stepwise Cox proportional hazards regression analyses, we identified significant associations between CM OS and four putatively functional SNPs: BRCA2 rs10492396 [AG vs. GG: adjusted hazard ratio (adjHR)=1.85, 95% confident interval (CI)=1.16-2.95, P=0.010], rs206118 (CC vs. TT+TC: adjHR=2.44, 95% CI=1.27-4.67, P=0.007), rs3752447 (CC vs. TT+TC: adjHR=2.10, 95% CI=1.38-3.18, P=0.0005), and FANCA rs62068372 (TT vs. CC+CT: adjHR=1.85, 95% CI=1.27-2.69, P=0.001). Moreover, patients with an increasing number of unfavorable genotypes (NUG) of these loci had markedly reduced OS and melanoma-specific survival (MSS). The final model incorporating with NUG, tumor stage and Breslow thickness showed an improved discriminatory ability to classify both 5-year OS and 5-year MSS. Additional investigations, preferably prospective studies, are needed to validate our findings.

Metalloproteinases in melanoma

Tumour cell adhesion, motility, proteolytic activities and cell receptors have important roles in cancer invasion. These processes are involved from early development of melanoma within the epidermis, to tumour cell invasion of the underlying tissue until intravasation of lymphatic or blood vessels, and thereafter, dissemination into distant organs occur. The activity of several proteolytic enzymes was shown to be pivotal in promoting melanoma cell invasion. These enzymes not only remodel the extracellular matrix, but also release active factors and shed cell surface receptors thereby mediating melanoma cross-communication with their microenvironment. This leads to the generation of a favourable environment for melanoma growth. Several proteases are involved in melanoma invasion and include serine, cysteine proteases, matrix metalloproteases (MMPs) and the disintegrin and metalloproteases (ADAMs). This study summarises the current knowledge on the role of metalloproteinases, MMPs and ADAMs, in melanoma.

Melanoma risk loci as determinants of melanoma recurrence and survival

Background

Steadily high melanoma mortality rates urge for the availability of novel biomarkers with a more personalized ability to predict melanoma clinical outcomes. Germline risk variants are promising candidates for this purpose; however, their prognostic potential in melanoma has never been systematically tested.

Methods

We examined the effect of 108 melanoma susceptibility single nucleotide polymorphisms (SNPs), associated in recent GWAS with melanoma and melanoma-related phenotypes, on recurrence-free survival (RFS) and overall survival (OS), in 891 prospectively accrued melanoma patients. Cox proportional hazards models (Cox PH) were used to test the associations between 108 melanoma risk SNPs and RFS and OS adjusted by age at diagnosis, gender, tumor stage, histological subtype and other primary tumor characteristics.

Results

We identified significant associations for rs7538876 (RCC2) with RFS (HR = 1.48, 95% CI = 1.20-1.83, p = 0.0005) and rs9960018 (DLGAP1) with both RFS and OS (HR = 1.43, 95% CI = 1.07-1.91, p = 0.01, HR = 1.52, 95% CI = 1.09-2.12, p = 0.01, respectively) using multivariable Cox PH models. In addition, we developed a logistic regression model that incorporates rs7538876, rs9960018, primary tumor histological type and stage at diagnosis that has an improved discriminatory ability to classify 3-year recurrence (AUC = 82%) compared to histological type and stage alone (AUC = 78%).

Conclusions

We identified associations between melanoma risk variants and melanoma outcomes. The significant associations observed for rs7538876 and rs9960018 suggest a biological implication of these loci in melanoma progression. The observed predictive patterns of associated variants with clinical end-points suggest for the first time the potential for utilization of genetic risk markers in melanoma prognostication.

Gene variants in angiogenesis and lymphangiogenesis and cutaneous melanoma progression

BACKGROUND

Angiogenesis and lymphangiogenesis are important in the progression of melanoma. We investigated associations between genetic variants in these pathways with sentinel lymph node (SLN) metastasis and mortality in 2 independent series of patients with melanoma.

METHODS

Participants at Moffitt Cancer Center were 552 patients, all Caucasian, with primary cutaneous melanoma referred for SLN biopsy. A total of 177 patients had SLN metastasis, among whom 60 died from melanoma. Associations between 238 single-nucleotide polymorphisms (SNP) in 26 genes and SLN metastasis were estimated as ORs and 95% confidence intervals (CI) using logistic regression. Competing risk regression was used to estimate HRs and 95% CI for each SNP and melanoma-specific mortality. We attempted to replicate significant findings using data from a genome-wide association study comprising 1,115 patients with melanoma who were referred for SLN biopsy from MD Anderson Cancer Center (MDACC), among whom 189 patients had SLN metastasis and 92 patients died from melanoma.

RESULTS

In the Moffitt dataset, we observed significant associations in 18 SNPs with SLN metastasis and 17 SNPs with mortality. Multiple SNPs in COL18A1, EGF receptor (EGFR), FLT1, interleukin (IL)-10, platelet-derived growth factor D (PDGFD), PIK3CA, and toll-like receptor (TLR)-3 were associated with the risk of SLN metastasis and/or patient mortality. The MDACC data set replicated an association between mortality and rs2220377 in PDGFD. Furthermore, in a meta-analysis, 3 additional SNPs were significantly associated with SLN metastasis (EGFR rs723526 and TLR3 rs3775292) and melanoma-specific death (TLR3 rs7668666).

CONCLUSIONS

These findings suggest that genetic variation in angiogenesis and lymphangiogenesis contributes to regional nodal metastasis and progression of melanoma.

IMPACT

Additional research attempting to replicate these results is warranted.

Interplay of microRNAs, transcription factors and target genes: linking dynamic expression changes to function

MicroRNAs (miRNAs) are ubiquitously expressed small non-coding RNAs that, in most cases, negatively regulate gene expression at the post-transcriptional level. miRNAs are involved in fine-tuning fundamental cellular processes such as proliferation, cell death and cell cycle control and are believed to confer robustness to biological responses. Here, we investigated simultaneously the transcriptional changes of miRNA and mRNA expression levels over time after activation of the Janus kinase/Signal transducer and activator of transcription (Jak/STAT) pathway by interferon-γ stimulation of melanoma cells. To examine global miRNA and mRNA expression patterns, time-series microarray data were analysed. We observed delayed responses of miRNAs (after 24-48 h) with respect to mRNAs (12-24 h) and identified biological functions involved at each step of the cellular response. Inference of the upstream regulators allowed for identification of transcriptional regulators involved in cellular reactions to interferon-γ stimulation. Linking expression profiles of transcriptional regulators and miRNAs with their annotated functions, we demonstrate the dynamic interplay of miRNAs and upstream regulators with biological functions. Finally, our data revealed network motifs in the form of feed-forward loops involving transcriptional regulators, mRNAs and miRNAs. Additional information obtained from integrating time-series mRNA and miRNA data may represent an important step towards understanding the regulatory principles of gene expression.