Tissue-Specific Methylation Biosignatures for Monitoring Diseases: An In Silico Approach

Tissue-specific gene methylation events are key to the pathogenesis of several diseases and can be utilized for diagnosis and monitoring. Here, we established an in silico pipeline to analyze high-throughput methylome datasets to identify specific methylation fingerprints in three pathological entities of major burden, i.e., breast cancer (BrCa), osteoarthritis (OA) and diabetes mellitus (DM). Differential methylation analysis was conducted to compare tissues/cells related to the pathology and different types of healthy tissues, revealing Differentially Methylated Genes (DMGs). Highly performing and low feature number biosignatures were built with automated machine learning, including: (1) a five-gene biosignature discriminating BrCa tissue from healthy tissues (AUC 0.987 and precision 0.987), (2) three equivalent OA cartilage-specific biosignatures containing four genes each (AUC 0.978 and precision 0.986) and (3) a four-gene pancreatic β-cell-specific biosignature (AUC 0.984 and precision 0.995). Next, the BrCa biosignature was validated using an independent ccfDNA dataset showing an AUC and precision of 1.000, verifying the biosignature’s applicability in liquid biopsy. Functional and protein interaction prediction analysis revealed that most DMGs identified are involved in pathways known to be related to the studied diseases or pointed to new ones. Overall, our data-driven approach contributes to the maximum exploitation of high-throughput methylome readings, helping to establish specific disease profiles to be applied in clinical practice and to understand human pathology.     

Development of Novel Ecto-Nucleotide Pyrophosphatase/Phosphodiesterase 1 (ENPP1) Inhibitors for Tumor Immunotherapy

The cyclic guanosine monophosphate–adenosine monophosphate synthase–stimulator of interferon genes–TANK-binding kinase 1–interferon regulating factor 3 (cGAS-STING-TBK1-IRF3) axis is now acknowledged as the major signaling pathway in innate immune responses. However, 2′,3′-cGAMP as a STING stimulator is easily recognized and degraded by ecto-nucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1), which reduces the effect of tumor immunotherapy and promotes metastatic progression. In this investigation, the structure-based virtual screening strategy was adopted to discover eight candidate compounds containing zinc-binding quinazolin-4(3H)-one scaffold as ENPP1 inhibitors. Subsequently, these novel inhibitors targeting ENPP1 were synthesized and characterized by NMR and high-resolution mass spectra (HRMS). In bioassays, 7-fluoro-2-(((5-methoxy-1H-imidazo[4,5-b]pyridin-2-yl)thio)methyl)quina-zolin-4(3H)-one(compound 4e) showed excellent activity against the ENPP1 at the molecular and cellular levels, with IC₅₀ values of 0.188 μM and 0.732 μM, respectively. Additionally, compound 4e had superior selectivity towards metastatic breast cancer cells (4T1) than towards normal cells (LO2 and 293T) in comparison with cisplatin, indicating that compound 4e can potentially be used in metastatic breast cancer therapy. On the other hand, compound 4e upgraded the expression levels of IFN-β in vivo by preventing the ENPP1 from hydrolyzing the cGAMP to stimulate a more potent innate immune response. Therefore, this compound might be applied to boost antitumor immunity for cancer immunotherapy. Overall, our work provides a strategy for the development of a promising drug candidate targeting ENPP1 for tumor immunotherapy.     

The Impact of [C16Pyr][Amp] on the Aggressiveness in Breast and Prostate Cancer Cell Lines

Breast (BrCa) and prostate (PCa) cancers are the most common malignancies in women and men, respectively. The available therapeutic options for these tumors are still not curative and have severe side effects. Therefore, there is an urgent need for more effective antineoplastic agents. Herein, BrCa, PCa, and benign cell lines were treated with two ionic liquids and two quinoxalines and functional experiments were performed—namely cell viability, apoptosis, cytotoxicity, and colony formation assays. At the molecular level, an array of gene expressions encompassing several molecular pathways were used to explore the impact of treatment on gene expression. Although both quinoxalines and the ionic liquid [C2OHMIM][Amp] did not show any effect on the BrCa and PCa cell lines, [C16Pyr][Amp] significantly decreased cell viability and colony formation ability, while it increased the apoptosis levels of all cell lines. Importantly, [C16Pyr][Amp] was found to be more selective for cancer cells and less toxic than cisplatin. At the molecular level, this ionic liquid was also associated with reduced expression levels of CPT2, LDHA, MCM2, and SKP2, in both BrCa and PCa cell lines. Hence, [C16Pyr][Amp] was shown to be a promising anticancer therapeutic agent for BrCa and PCa cell lines.     

SOCS3 Methylation Predicts a Poor Prognosis in HBV Infection-Related Hepatocellular Carcinoma

Suppressor of cytokine signaling 3 (SOCS3) plays crucial roles in JAK/STAT signaling pathway inhibition in hepatocellular carcinoma (HCC). However, the methylation status of SOCS3 in HBV infection-related HCC and the relationship between SOCS3 methylation and the clinical outcome remain unknown. Here, we reported that in HCC tumor tissues, two regions of the CpG island (CGI) in the SOCS3 promoter were subjected to methylation analysis and only the region close to the translational start site of SOCS3 was hypermethylated. In HCC tumor tissues, SOCS3 showed an increased methylation frequency and intensity compared with that in the adjacent non-tumor tissues. Moreover, SOCS3 expression was significantly down-regulated in HCC cell lines and tumor tissues, and this was inversely correlated with methylation. Kaplan–Meier curve analysis revealed that in patients with an hepatitis B virus (HBV) infection background, SOCS3 hypermethylation was significantly correlated with a poor clinical outcome of HCC patients. Our findings indicated that SOCS3 hypermethylation has already happened in non-tumor tissues and increased in both frequency and intensity in tumor tissues. This suggests that the methylation of SOCS3 could predict a poor prognosis in HBV infection-related HCC patients.     

ESR1 Methylation Measured in Cell-Free DNA to Evaluate Endocrine Resistance in Metastatic Breast Cancer Patients

ESR1 methylation was proposed as mechanism for endocrine resistance in metastatic breast cancer patients. To evaluate its potential as a minimally invasive biomarker, we investigated the feasibility of measuring ESR1 methylation in cell-free DNA (cfDNA) and its association with endocrine resistance. First, we provided evidence that demethylation in vitro restores ER expression. Subsequently, we found that ESR1 methylation in cfDNA was not enriched in endocrine-resistant versus endocrine-sensitive patients. Interestingly, we found a correlation between ESR1 methylation and age. Publicly available data confirm an age-related increase in ESR1 methylation in leukocytes, confounding the determination of the ESR1 methylation status of tumors using cfDNA.     

Increased Autotaxin Levels in Severe COVID-19, Correlating with IL-6 Levels,Endothelial Dysfunction Biomarkers,and Impaired Functions of Dendritic Cells

Autotaxin (ATX; ENPP2) is a secreted lysophospholipase D catalyzing the extracellular production of lysophosphatidic acid (LPA), a pleiotropic signaling phospholipid. Genetic and pharmacologic studies have previously established a pathologic role for ATX and LPA signaling in pulmonary injury, inflammation, and fibrosis. Here, increased ENPP2 mRNA levels were detected in immune cells from nasopharyngeal swab samples of COVID-19 patients, and increased ATX serum levels were found in severe COVID-19 patients. ATX serum levels correlated with the corresponding increased serum levels of IL-6 and endothelial damage biomarkers, suggesting an interplay of the ATX/LPA axis with hyperinflammation and the associated vascular dysfunction in COVID-19. Accordingly, dexamethasone (Dex) treatment of mechanically ventilated patients reduced ATX levels, as shown in two independent cohorts, indicating that the therapeutic benefits of Dex include the suppression of ATX. Moreover, large scale analysis of multiple single cell RNA sequencing datasets revealed the expression landscape of ENPP2 in COVID-19 and further suggested a role for ATX in the homeostasis of dendritic cells, which exhibit both numerical and functional deficits in COVID-19. Therefore, ATX has likely a multifunctional role in COVID-19 pathogenesis, suggesting that its pharmacological targeting might represent an additional therapeutic option, both during and after hospitalization.     

The Value of EphB2 Receptor and Cognate Ephrin Ligands in Prognostic and Predictive Assessments of Human Breast Cancer

Cell–cell communication proteins Eph and ephrin constitute the largest family of receptor tyrosine kinases (RTKs). They are distinguished by the fact that both receptors and ligands are membrane-bound, and both can drive intracellular signaling in their respective cells. Ever since these RTKs have been found to be involved in cancer development, strategies to target them therapeutically have been actively pursued. However, before this goal can be rationally achieved, the contributions of either Eph receptors or their ephrin ligands to cancer development and progression should be scrutinized in depth. To assess the clinical pertinence of this concern, we performed a systematic review and meta-analysis of the prognostic/predictive value of EphB2 and its multiple cognate ephrin ligands in breast cancer. We found that EphB2 has prognostic value, as indicated by the association of higher EphB2 expression levels with lower distant metastasis-free survival (DMFS), and the association of lower EphB2 expression levels with poorer relapse-free survival (RFS). We also found that higher EphB2 expression could be a prognostic factor for distant metastasis, specifically in the luminal subtypes of breast cancer. EFNB2 showed a marked correlation between higher expression levels and shorter DMFS. EFNA5 or EFNB1 overexpression is correlated with longer RFS. Increased EFNB1 expression is correlated with longer OS in lymph node (LN)-negative patients and the luminal B subtype. Higher levels of EFNB2 or EFNA5 are significantly correlated with shorter RFS, regardless of LN status. However, while this correlation with shorter RFS is true for EFNB2 in all subtypes except basal, it is also true for EFNA5 in all subtypes except HER2+. The analysis also points to possible predictive value for EphB2. In systemically treated patients who have undergone either endocrine therapy or chemotherapy, we found that higher expression of EphB2 is correlated with better rates of RFS. Bearing in mind the limitations inherent to any mRNA-based profiling method, we complemented our analysis with an immunohistochemical assessment of expression levels of both the EphB2 receptor and cognate ephrin ligands. We found that the latter are significantly more expressed in cancers than in normal tissues, and even more so in invasive and metastatic samples than in ductal carcinoma in situ (DCIS). Finally, in an in vitro cellular model of breast cancer progression, based on H-Ras-transformation of the MCF10A benign mammary cell line, we observed dramatic increases in the mRNA expression of EphB2 receptor and EFNB1 and EFNB2 ligands in transformed and invasive cells in comparison with their benign counterparts. Taken together, these data show the clinical validity of a model whereby EphB2, along with its cognate ephrin ligands, have dual anti- and pro-tumor progression effects. In so doing, they reinforce the necessity of further biological investigations into Ephs and ephrins, prior to using them in targeted therapies.     

Activation of VCAM-1 and Its Associated Molecule CD44 Leads to Increased Malignant Potential of Breast Cancer Cells

VCAM-1 (CD106), a transmembrane glycoprotein, was first reported to play an important role in leukocyte adhesion, leukocyte transendothelial migration and cell activation by binding to integrin VLA-1 (α4β1). In the present study, we observed that VCAM-1 expression can be induced in many breast cancer epithelial cells by cytokine stimulation in vitro and its up-regulation directly correlated with advanced clinical breast cancer stage. We found that VCAM-1 over-expression in the NMuMG breast epithelial cells controls the epithelial and mesenchymal transition (EMT) program to increase cell motility rates and promote chemoresistance to doxorubicin and cisplatin in vitro. Conversely, in the established MDAMB231 metastatic breast cancer cell line, we confirmed that knockdown of endogenous VCAM-1 expression reduced cell proliferation and inhibited TGFβ1 or IL-6 mediated cell migration, and increased chemosensitivity. Furthermore, we demonstrated that knockdown of endogenous VCAM-1 expression in MDAMB231 cells reduced tumor formation in a SCID xenograft mouse model. Signaling studies showed that VCAM-1 physically associates with CD44 and enhances CD44 and ABCG2 expression. Our findings uncover the possible mechanism of VCAM-1 activation facilitating breast cancer progression, and suggest that targeting VCAM-1 is an attractive strategy for therapeutic intervention.     

Promoter Methylation Changes in KRT17: A Novel Epigenetic Marker for Wool Production in Angora Rabbit

Wool production is an important economic trait of Angora rabbits. Exploring molecular markers related to wool production is one of the essentials of Angora rabbits’ breeding. KRT17 (Keratin 17) is an important gene of hair follicle development, which must be explored for genetic/epigenetic variation to assess its effect on wool production. Based on the effective wool production data of 217 Angora rabbits, the high and low yield groups were screened with 1.5 standard deviations of the population mean. The full-length sequence of KRT17 was obtained by rapid amplification of cDNA ends technology, and the polymorphism was analyzed in the promoter, exon, and intron regions by direct sequencing. KRT17, SP1 over-expression plasmids, and siRNA were constructed and transfected into dermal papilla cells. The mRNA expressions of relevant genes were analyzed by RT-qPCR. The methylation level of the KRT17 promoter was determined by Bisulfite Sequencing PCR. Dual-luciferase system, site-directed mutagenesis, and electrophoretic mobility shift assays were used to analyze the binding relationship between SP1 and the promoter of KRT17. The structure map of KRT17 was drawn, and no SNPs were found in the promoter, exon, and intron, indicating a relatively conserved structure of KRT17. Expression of KRT17 was significantly higher in cutaneous tissues than in other tissues and was significantly upregulated in the high-yield group compared to the low-yield group (p < 0.05). Furthermore, the overall high methylation levels of KRT17 CpG I and CpG III showed significant association with low wool yield; the methylation levels of 5 CpG locus (CpG I site 4 and CpG III site 2–5) were significantly different between the high and low yield groups (p < 0.05). The methylation levels of 3 CpG locus (CpG I site 4 and CpG III site 4, 14) showed a significant correlation with KRT17 expression (p < 0.05). Overall, CpG III site 4 significantly affects wool production and KRT17 expressions (p < 0.05). This site promotes SP1 binding to the KRT17 promoter region (CGCTACGCCC) to positively regulate the KRT17 expression. KRT17 CpG III site 4 can be used as candidate epigenetic markers for the breeding of high wool-producing Angora rabbits.     

KIT Expression Is Regulated by DNA Methylation in Uveal Melanoma Tumors

Uveal melanoma (UM) is an ocular tumor with a dismal prognosis. Despite the availability of precise molecular and cytogenetic techniques, clinicopathologic features with limited accuracy are widely used to predict metastatic potential. In 51 UM tissues, we assessed a correlation between the expression of nine proteins evaluated by immunohistochemistry (IHC) (Melan-A, S100, HMB45, Cyclin D1, Ki-67, p53, KIT, BCL2, and AIFM1) and the presence of UM-specific chromosomal rearrangements measured by multiplex ligation-dependent probe amplification (MLPA), to find IHC markers with increased prognostic information. Furthermore, mRNA expression and DNA methylation values were extracted from the whole-genome data, achieved by analyzing 22 fresh frozen UM tissues. KIT positivity was associated with monosomy 3, increasing the risk of poor prognosis more than 17-fold (95% CI 1.53–198.69, p = 0.021). A strong negative correlation was identified between mRNA expression and DNA methylation values for 12 of 20 analyzed positions, five located in regulatory regions of the KIT gene (r = −0.658, p = 0.001; r = −0.662, p = 0.001; r = −0.816; p < 0.001; r = −0.689, p = 0.001; r = −0.809, p < 0.001, respectively). DNA methylation β values were also inversely associated with KIT protein expression (p = 0.001; p = 0.001; p = 0.015; p = 0.025; p = 0.002). Our findings, showing epigenetic deregulation of KIT expression, may contribute to understanding the past failure to therapeutically target KIT in UM.     

CES2, ABCG2, TS and Topo-I Primary and Synchronous Metastasis Expression and Clinical Outcome in Metastatic Colorectal Cancer Patients Treated with First-Line FOLFIRI Regimen

Enzymatic activation of irinotecan (CPT-11) is due to carboxylesterase (CES), and its pharmacological behavior is influenced by drug resistance-related proteins. We previously reported that the clinical response and prognosis of metastatic colorectal cancer (mCRC) patients did not differ in tumors with different thymidylate synthase (TS) or topoisomerase-I (Topo-I) expression. Using immunohistochemistry (IHC), we evaluated the biological role of CES2 and the expression of breast cancer resistance protein (BCRP/ABCG2) in 58 consecutive mCRC patients, who had undergone a first-line CPT-11/5-FU/leucovirin (FOLFIRI) regimen. The expression of these proteins was also examined in a group of synchronous lymph nodes and liver metastases. Furthermore, all samples were revaluated for TS and Topo-I expression. High expression of CES2, ABCG2, TS and Topo-I was observed in 55%, 56%, 38% and 49% of patients, respectively. There was a significant association between high TS and high ABCG2 expression (p = 0.049). Univariate analysis showed that only TS expression significantly impacted on time to progression (p = 0.005). Moreover, Cox’ multivariate analysis revealed that TS expression was significantly associated with overall survival (p = 0.01). No significant correlation was found between investigated markers expression and clinical response. Topo-I expression resulted in being significantly higher in liver metastases with respect to the corresponding primary tumors (p < 0.0001), emphasizing the role of Topo-I expression in metastatic cancer biology. In primary tumor tissues, CES2 expression tended to be higher than that observed in liver metastasis tissues (p = 0.05). These preliminary data may suggest CES2 over-expression as a potential marker of malignant phenotype. In light of these findings, we suggest that Topo-I expression together with TS expression could be associated with metastatic progression of CRC. Further studies are warranted with the aim of evaluating the potential predictive and prognostic role of CES2 and ABCG2 in larger series of patients.     

Analysis of Intrinsic Breast Cancer Subtypes: The Clinical Utility of Epigenetic Biomarkers and TP53 Mutation Status in Triple-Negative Cases

This study aimed at analyzing the DNA methylation pattern and TP53 mutation status of intrinsic breast cancer (BC) subtypes for improved characterization and survival prediction. DNA methylation of 17 genes was tested by methylation-specific PCR in 116 non-familial BRCA mutation-negative BC and 29 control noncancerous cases. At least one gene methylation was detected in all BC specimens and a 10-gene panel statistically significantly separated tumors from noncancerous breast tissues. Methylation of FILIP1L and MT1E was predominant in triple-negative (TN) BC, while other BC subtypes were characterized by RASSF1, PRKCB, MT1G, APC, and RUNX3 hypermethylation. TP53 mutation (TP53-mut) was found in 38% of sequenced samples and mainly affected TN BC cases (87%). Cox analysis revealed that TN status, age at diagnosis, and RUNX3 methylation are independent prognostic factors for overall survival (OS) in BC. The combinations of methylated biomarkers, RUNX3 with MT1E or FILIP1L, were also predictive for shorter OS, whereas methylated FILIP1L was predictive of a poor outcome in the TP53-mut subgroup. Therefore, DNA methylation patterns of specific genes significantly separate BC from noncancerous breast tissues and distinguishes TN cases from non-TN BC, whereas the combination of two-to-three epigenetic biomarkers can be an informative tool for BC outcome predictions.     

Expression and Function of StAR in Cancerous and Non-Cancerous Human and Mouse Breast Tissues: New Insights into Diagnosis and Treatment of Hormone-Sensitive Breast Cancer

Breast cancer (BC) is primarily triggered by estrogens, especially 17β-estradiol (E2), which are synthesized by the aromatase enzyme. While all steroid hormones are derived from cholesterol, the rate-limiting step in steroid biosynthesis is mediated by the steroidogenic acute regulatory (StAR) protein. Herein, we demonstrate that StAR mRNA expression was aberrantly high in human hormone-dependent BC (MCF7, MDA-MB-361, and T-47D), modest in hormone-independent triple negative BC (TNBC; MDA-MB-468, BT-549, and MDA-MB-231), and had little to none in non-cancerous mammary epithelial (HMEC, MCF10A, and MCF12F) cells. In contrast, these cell lines showed abundant expression of aromatase (CYP19A1) mRNA. Immunofluorescence displayed qualitatively similar patterns of both StAR and aromatase expression in various breast cells. Additionally, three different transgenic (Tg) mouse models of spontaneous breast tumors, i.e., MMTV-Neu, MMTV-HRAS, and MMTV-PyMT, demonstrated markedly higher expression of StAR mRNA/protein in breast tumors than in normal mammary tissue. While breast tumors in these mouse models exhibited higher expression of ERα, ERβ, and PR mRNAs, their levels were undetected in TNBC tumors. Accumulation of E2 in plasma and breast tissues, from MMTV-PyMT and non-cancerous Tg mice, correlated with StAR, but not with aromatase, signifying the importance of StAR in governing E2 biosynthesis in mammary tissue. Treatment with a variety of histone deacetylase inhibitors (HDACIs) in primary cultures of enriched breast tumor epithelial cells, from MMTV-PyMT mice, resulted in suppression of StAR and E2 levels. Importantly, inhibition of StAR, concomitant with E2 synthesis, by various HDACIs, at clinical and preclinical doses, in MCF7 cells, indicated therapeutic relevance of StAR in hormone-dependent BCs. These findings provide insights into the molecular events underlying the differential expression of StAR in human and mouse cancerous and non-cancerous breast cells/tissues, highlighting StAR could serve not only as a novel diagnostic maker but also as a therapeutic target for the most prevalent hormone-sensitive BCs.