DNA methylation alterations at RE1-silencing transcription factor binding sites and their flanking regions in cancer.

BACKGROUND: DNA methylation changes, frequent early events in cancer, can modulate the binding of transcription factors. RE1-silencing transcription factor (REST) plays a fundamental role in regulating the expression of neuronal genes, and in particular their silencing in non-neuronal tissues, by inducing chromatin modifications, including DNA methylation changes, not only in the proximity of its binding sites but also in the flanking regions. REST has been found aberrantly expressed in brain cancer and other cancer types. In this work, we investigated DNA methylation alterations at REST binding sites and their flanking regions in a brain cancer (pilocytic astrocytoma), two gastrointestinal tumours (colorectal cancer and biliary tract cancer) and a blood cancer (chronic lymphocytic leukemia). RESULTS: Differential methylation analyses focused on REST binding sites and their flanking regions were conducted between tumour and normal samples from our experimental datasets analysed by Illumina microarrays and the identified alterations were validated using publicly available datasets. We discovered distinct DNA methylation patterns between pilocytic astrocytoma and the other cancer types in agreement with the opposite oncogenic and tumour suppressive role of REST in glioma and non-brain tumours. CONCLUSIONS: Our results suggest that these DNA methylation alterations in cancer may be associated with REST dysfunction opening the enthusiastic possibility to develop novel therapeutic interventions based on the modulation of this master regulator in order to restore the aberrant methylation of its target regions into a normal status.

HLA-C dysregulation as a possible mechanism of immune evasion in SARS-CoV-2 and other RNA-virus infections.

One of the mechanisms by which viruses can evade the host's immune system is to modify the host's DNA methylation pattern. This work aims to investigate the DNA methylation and gene expression profile of COVID-19 patients, divided into symptomatic and asymptomatic, and healthy controls, focusing on genes involved in the immune response. In this study, changes in the methylome of COVID-19 patients' upper airways cells, the first barrier against respiratory infections and the first cells presenting viral antigens, are shown for the first time. Our results showed alterations in the methylation pattern of genes encoding proteins implicated in the response against pathogens, in particular the HLA-C gene, also important for the T-cell mediated memory response. HLA-C expression significantly decreases in COVID-19 patients, especially in those with a more severe prognosis and without other possibly confounding co-morbidities. Moreover, our bionformatic analysis revealed that the identified methylation alteration overlaps with enhancers regulating HLA-C expression, suggesting an additional mechanism exploited by SARS-CoV-2 to inhibit this fundamental player in the host's immune response. HLA-C could therefore represent both a prognostic marker and an excellent therapeutic target, also suggesting a preventive intervention that conjugate a virus-specific antigenic stimulation with an adjuvant increasing the T-cell mediated memory response.

HOXD8 hypermethylation as a fully sensitive and specific biomarker for biliary tract cancer detectable in tissue and bile samples.

BACKGROUND: Biliary tract cancers (BTC) are rare but highly aggressive tumours with poor prognosis, usually detected at advanced stages. Herein, we aimed at identifying BTC-specific DNA methylation alterations. METHODS: Study design included statistical power and sample size estimation. A genome-wide methylation study of an explorative cohort (50 BTC and ten matched non-tumoral tissue samples) has been performed. BTC-specific altered CpG islands were validated in over 180 samples (174 BTCs and 13 non-tumoral controls). The final biomarkers, selected by a machine-learning approach, were validated in independent tissue (18 BTCs, 14 matched non-tumoral samples) and bile (24 BTCs, five non-tumoral samples) replication series, using droplet digital PCR. RESULTS: We identified and successfully validated BTC-specific DNA methylation alterations in over 200 BTC samples. The two-biomarker panel, selected by an in-house algorithm, showed an AUC > 0.97. The best-performing biomarker (chr2:176993479-176995557), associated with HOXD8, a pivotal gene in cancer-related pathways, achieved 100% sensitivity and specificity in a new series of tissue and bile samples. CONCLUSIONS: We identified a novel fully efficient BTC biomarker, associated with HOXD8 gene, detectable both in tissue and bile by a standardised assay ready-to-use in clinical trials also including samples from non-invasive matrices.

Colorectal cancer promoter methylation alteration affects the expression of glutamate ionotropic receptor AMPA type subunit 4 alternative isoforms potentially relevant in colon tissue.

DNA methylation alterations are early events during tumourigenesis, affecting genes involved in the crosstalk between cells and surroundings in colorectal cancer (CRC). Among these genes, GRIA4, Glutamate Ionotropic Receptor AMPA Type Subunit 4, displays hypermethylation in the promoter region, and is an early diagnostic biomarker. It is well known that methylation can also affect alternative transcription. The purpose of this study is to evaluate the expression, at transcript and protein level, of GRIA4 main isoforms (the canonical one and a short variant) in 23 CRC and matched normal samples, of which we previously verified the methylation status. We further predicted miRNA/transcript target interactions as a possible post-transcriptional regulation using bioinformatics tools. As expected, downregulation of both variants has been observed in tumours. Interestingly, in contrast to what observed at transcriptional level, the GluR4 protein short isoform displayed higher expression than the canonical one either in normal or tumoural tissues. This may be explained by miRNA specifically targeting the canonical isoform. Our study is the first one that shows the expression of both isoforms in colon tissues. To note, the evident expression of the short isoform suggests a functional role in intestinal cell biology.

Colorectal Cancer Early Detection in Stool Samples Tracing CpG Islands Methylation Alterations Affecting Gene Expression.

Colorectal cancer (CRC) is a major cause of cancer mortality. Early diagnosis is relevant for its prevention and treatment. Since DNA methylation alterations are early events in tumourigenesis and can be detected in cell-free DNA, they represent promising biomarkers for early CRC diagnosis through non-invasive methods. In our previous work, we identified 74 early altered CpG islands (CGIs) associated with genes involved in cell cross-talking and cell signalling pathways. The aim of this work was to test whether methylation-based biomarkers could be detected in non-invasive matrices. Our results confirmed methylation alterations of GRIA4 and VIPR2 in CRC tissues, using MethyLight, as well as in stool samples, using a much more sensitive technique as droplet digital PCR. Furthermore, we analysed expression levels of selected genes whose promoter CGIs were hypermethylated in CRC, detecting downregulation at mRNA and protein levels in CRC tissue for GRIA4, VIPR2, SPOCK1 and SLC6A3. Most of these genes were already lowly expressed in colon normal tissues supporting the idea that cancer DNA methylation targets genes already barely expressed in the matched normal tissues. Our study suggests GRIA4 and VIPR2 as biomarkers for early CRC diagnosis using stool samples and confirms downregulation of genes hypermethylated in CRC.

ELMOD3-SH2D6 gene fusion as a possible co-star actor in autism spectrum disorder scenario.

Autism spectrum disorder (ASD) is a group of neurodevelopmental disorders characterized by high heritability. It is known that genetic factors contribute to ASD pathogenesis. In particular, copy number variants (CNVs) are involved in ASD susceptibility and can affect gene expression regulation. 2p11.2 microdeletions encompassing ELMOD3, CAPG and SH2D6 genes have been described in four unrelated ASD families. The present study revealed that this microdeletion is responsible for the production of a chimeric transcript generated from the fusion between ELMOD3 and SH2D6. The identified transcript showed significantly higher expression levels in subjects carrying the deletion compared to control subjects, suggesting that it is not subjected to nonsense-mediated decay and might encode for a chimeric protein. In conclusion, this study suggests the possible involvement of this gene fusion, together with the other previously identified variants, in ASD.

Methylation alteration of SHANK1 as a predictive, diagnostic and prognostic biomarker for chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is a clinically heterogeneous disease characterized by the clonal expansion of malignant B cells. To predict the clinical course of the disease, the identification of diagnostic biomarkers is urgently needed. Aberrant methylation patterns may predict CLL development and its course, being very early changes during carcinogenesis. Our aim was to identify CLL specific methylation patterns and to evaluate whether methylation aberrations in selected genes are associated with changes in gene expression. Here, by performing a genome-wide methylation analysis, we identified several CLL-specific methylation alterations. We focused on the most altered one, at a CpG island located in the body of SHANK1 gene, in our CLL cases compared to healthy controls. This methylation alteration was successfully validated in a larger cohort including 139 CLL and 20 control in silico samples. We also found a positive correlation between SHANK1 methylation level and absolute lymphocyte count, in particular CD19+ B cells, in CLL patients. Moreover, we were able to detect gains of methylation at SHANK1 in blood samples collected years prior to diagnosis. Overall, our results suggest methylation alteration at this SHANK1 CpG island as a biomarker for risk and diagnosis of CLL, and also in the personalized quantification of tumor aggressiveness.

Evaluation of oxidative stress mechanisms and the effects of phytotherapic extracts on Parkinson's disease Drosophila PINK1B9 model.

Oxidative stress is commonly observed in both idiopathic and genetic cases of Parkinson's disease (PD). It plays an important role in the degeneration of dopaminergic neurons, and it has been associated with altered telomere length (TL). There is currently no cure for PD, and extracts of antioxidative plant, such as Mucuna pruriens and Withania somnifera, are commonly used in Ayurveda to treat patients with PD. In this study, we evaluated 2 enzymatic markers of oxidative stress, glutathione (GSH) system and superoxide dismutase (SOD), and TL in a Drosophila melanogaster model for PD [phosphatase and tensin homolog-induced putative kinase 1 (PINK1)B9]. This evaluation was also performed after treatment with the phytoextracts. PINK1B9 mutants showed a decrease in GSH amount and SOD activity and unexpected longer telomeres compared with wild-type flies. M. pruriens treatment seemed to have a beneficial effect on the oxidative stress conditions. On the other hand, W. somnifera treatment did not show any improvements in the studied oxidative stress mechanisms and even seemed to favor the selection of flies with longer telomeres. In summary, our study suggests the importance of testing antioxidant phytoextracts in a PINK1B9 model to identify beneficial effects for PD.-Baroli, B., Loi, E., Solari, P., Kasture, A., Moi, L., Muroni, P., Kasture, S., Setzu, M. D., Liscia, A., Zavattari, P. Evaluation of oxidative stress mechanisms and the effects of phytotherapic extracts on Parkinson's disease Drosophila PINK1B9 model.

Clustered protocadherins methylation alterations in cancer.

BACKGROUND: Clustered protocadherins (PCDHs) map in tandem at human chromosome 5q31 and comprise three multi-genes clusters: α-, ß- and γ-PCDH. The expression of this cluster consists of a complex mechanism involving DNA hub formation through DNA-CCTC binding factor (CTCF) interaction. Methylation alterations can affect this interaction, leading to transcriptional dysregulation. In cancer, clustered PCDHs undergo a mechanism of long-range epigenetic silencing by hypermethylation. RESULTS: In this study, we detected frequent methylation alterations at CpG islands associated to these clustered PCDHs in all the solid tumours analysed (colorectal, gastric and biliary tract cancers, pilocytic astrocytoma), but not hematologic neoplasms such as chronic lymphocytic leukemia. Importantly, several altered CpG islands were associated with CTCF binding sites. Interestingly, our analysis revealed a hypomethylation event in pilocytic astrocytoma, suggesting that in neuronal tissue, where PCDHs are highly expressed, these genes become hypomethylated in this type of cancer. On the other hand, in tissues where PCDHs are lowly expressed, these CpG islands are targeted by DNA methylation. In fact, PCDH-associated CpG islands resulted hypermethylated in gastrointestinal tumours. CONCLUSIONS: Our study highlighted a strong alteration of the clustered PCDHs methylation pattern in the analysed solid cancers and suggested these methylation aberrations in the CpG islands associated with PCDH genes as powerful diagnostic biomarkers.

Analysis of a Sardinian Multiplex Family with Autism Spectrum Disorder Points to Post-Synaptic Density Gene Variants and Identifies CAPG as a Functionally Relevant Candidate Gene.

Autism spectrum disorders (ASDs) are a group of neurodevelopmental disorders with high heritability, although their underlying genetic factors are still largely unknown. Here we present a comprehensive genetic characterization of two ASD siblings from Sardinia by genome-wide copy number variation analysis and whole exome sequencing (WES), to identify novel genetic alterations associated with this disorder. Single nucleotide polymorphism (SNP) array data revealed a rare microdeletion involving CAPG, ELMOD3, and SH2D6 genes, in both siblings. CAPG encodes for a postsynaptic density (PSD) protein known to regulate spine morphogenesis and synaptic formation. The reduced CAPG mRNA and protein expression levels in ASD patients, in the presence of hemizygosity or a particular genetic and/or epigenetic background, highlighted the functional relevance of CAPG as a candidate gene for ASD. WES analysis led to the identification in both affected siblings of a rare frameshift mutation in VDAC3, a gene intolerant to loss of function mutation, encoding for a voltage-dependent anion channel localized on PSD. Moreover, four missense damaging variants were identified in genes intolerant to loss of function variation encoding for PSD proteins: PLXNA2, KCTD16, ARHGAP21, and SLC4A1. This study identifies CAPG and VDAC3 as candidate genes and provides additional support for genes encoding PSD proteins in ASD susceptibility.

Integrated DNA methylation analysis identifies topographical and tumoral biomarkers in pilocytic astrocytomas.

Pilocytic astrocytoma (PA) is the most common glioma in pediatric patients and occurs in different locations. Chromosomal alterations are mostly located at chromosome 7q34 comprising the BRAF oncogene with consequent activation of the mitogen-activated protein kinase pathway. Although genetic and epigenetic alterations characterizing PA from different localizations have been reported, the role of epigenetic alterations in PA development is still not clear. The aim of this study was to investigate whether distinctive methylation patterns may define biologically relevant groups of PAs. Integrated DNA methylation analysis was performed on 20 PAs and 4 normal brain samples by Illumina Infinium HumanMethylation27 BeadChips. We identified distinct methylation profiles characterizing PAs from different locations (infratentorial vs supratentorial) and tumors with onset before and after 3 years of age. These results suggest that PA may be related to the specific brain site where the tumor arises from region-specific cells of origin. We identified and validated in silico the methylation alterations of some CpG islands. Furthermore, we evaluated the expression levels of selected differentially methylated genes and identified two biomarkers, one, IRX2, related to the tumor localization and the other, TOX2, as tumoral biomarker.

Colorectal cancer early methylation alterations affect the crosstalk between cell and surrounding environment, tracing a biomarker signature specific for this tumor.

Colorectal cancer (CRC) develops through the accumulation of both genetic and epigenetic alterations. However, while the former are already used as prognostic and predictive biomarkers, the latter are less well characterized. Here, performing global methylation analysis on both CRCs and adenomas by Illumina Infinium HumanMethylation450 Bead Chips, we identified a panel of 74 altered CpG islands, demonstrating that the earliest methylation alterations affect genes coding for proteins involved in the crosstalk between cell and surrounding environment. The panel discriminates CRCs and adenomas from peritumoral and normal mucosa with very high specificity (100%) and sensitivity (99.9%). Interestingly, over 70% of the hypermethylated islands resulted in downregulation of gene expression. To establish the possible usefulness of these non-invasive markers for detection of colon cancer, we selected three biomarkers and identified the presence of altered methylation in stool DNA and plasma cell-free circulating DNA from CRC patients.

Estimation of a significance threshold for epigenome-wide association studies.

Epigenome-wide association studies (EWAS) are designed to characterise population-level epigenetic differences across the genome and link them to disease. Most commonly, they assess DNA-methylation status at cytosine-guanine dinucleotide (CpG) sites, using platforms such as the Illumina 450k array that profile a subset of CpGs genome wide. An important challenge in the context of EWAS is determining a significance threshold for declaring a CpG site as differentially methylated, taking multiple testing into account. We used a permutation method to estimate a significance threshold specifically for the 450k array and a simulation extrapolation approach to estimate a genome-wide threshold. These methods were applied to five different EWAS datasets derived from a variety of populations and tissue types. We obtained an estimate of α=2.4×10-7 for the 450k array, and a genome-wide estimate of α=3.6×10-8. We further demonstrate the importance of these results by showing that previously recommended sample sizes for EWAS should be adjusted upwards, requiring samples between ∼10% and ∼20% larger in order to maintain type-1 errors at the desired level.

KIAA1549:BRAF fusion gene in pediatric brain tumors of various histogenesis.

The KIAA1549:BRAF fusion gene is considered a driver genetic event in pilocytic astrocytoma. We investigated a series of 69 pediatric brain neoplasms of diverse histogenesis and grade using the RT-PCR and sequencing. We detected the KIAA1549:BRAF fusion gene in five of 34 non-PA tumors (14.7%), that is, one glioblastoma, one anaplastic astrocytoma, one anaplastic pleomorphic xanthoastrocytoma, 1 ependymoma, and 1 Atypical Teratoid Rhabdoid Tumor. Our study showed that the K-B, although uncommon, it can be detected in non-PA tumors of various histogenesis and grading.

Response of recurrent BRAFV600E mutated ganglioglioma to Vemurafenib as single agent.

BACKGROUND: Ganglioglioma (GG) and pilocytic astrocytoma (PA) represent the most frequent low-grade gliomas (LGG) occurring in paediatric age. LGGs not amenable of complete resection (CR) represent a challenging subgroup where traditional treatments often fail. Activation of the MAP Kinase (MAPK) pathway caused by the BRAFV600E mutation or the KIAA1549-BRAF fusion has been reported in pediatric GG and PA, respectively. CASE PRESENTATION: We report on a case of BRAFV600E mutated cervicomedullary GG treated with standard chemotherapy and surgery. After multiple relapse, BRAF status was analyzed by immunohistochemistry and sequencing showing a BRAFV600E mutation. Treatment with Vemurafenib as single agent was started. For the first time, a radiological and clinical response was obtained after 3 months of treatment and sustained after 6 months. CONCLUSION: Our experience underline the importance of understanding the driver molecular alterations of LGG and suggests a role for Vemurafenib in the treatment of pediatric GG not amenable of complete surgical resection.

KIAA1549-BRAF fusions and IDH mutations can coexist in diffuse gliomas of adults.

KIAA1549-BRAF fusion gene and isocitrate dehydrogenase (IDH) mutations are considered two mutually exclusive genetic events in pilocytic astrocytomas and diffuse gliomas, respectively. We investigated the presence of the KIAA1549-BRAF fusion gene in conjunction with IDH mutations and 1p/19q loss in 185 adult diffuse gliomas. Moreover BRAF(v600E) mutation was also screened. The KIAA1549-BRAF fusion gene was evaluated by reverse-transcription polymerase chain reaction (RT-PCR) and sequencing. We found IDH mutations in 125 out 175 cases (71.4%). There were KIAA1549-BRAF fusion gene in 17 out of 180 (9.4%) cases and BRAF(v600E) in 2 out of 133 (1.5%) cases. In 11 of these 17 cases, both IDH mutations and the KIAA1549-BRAF fusion were present, as independent molecular events. Moreover, 6 of 17 cases showed co-presence of 1p/19q loss, IDH mutations and KIAA1549-BRAF fusion. Among the 17 cases with KIAA1549-BRAF fusion gene 15 (88.2%) were oligodendroglial neoplasms. Similarly, the two cases with BRAF(v600E) mutation were both oligodendroglioma and one had IDH mutations and 1p/19q co-deletion. Our results suggest that in a small fraction of diffuse gliomas, KIAA1549-BRAF fusion gene and BRAF(v600E) mutation may be responsible for deregulation of the Ras-RAF-ERK signaling pathway. Such alterations are more frequent in oligodendroglial neoplasm and may be co-present with IDH mutations and 1p/19q loss.

rs9939609 in the FTO gene is associated with obesity but not with several biochemical parameters in Sardinian obese children.

Several studies have reported an association of the intronic single nucleotide polymorphism (SNP) rs9939609 of the fat mass and obesity-associated (FTO) gene with obesity and with a number of obesity-related features. We studied the association of rs9939609 with obesity in 912 obese children and adolescents (426 males and 486 females, mean ± SD age 10.5 ± 3.3 years) and in 543 normal weight subjects. A number of biochemical and clinical parameters was also evaluated in 700 of these patients. In the obese group, mean body mass index standard deviation score (BMI-SDS) was similar between the three genotypes. The A allele was present in 55% of the patients' and in 43% of controls' chromosomes. The distribution of heterozygotes was similar between patients and controls (47%), while the distribution of AA homozygotes was significantly higher in patients (31% vs. 20%). Logistic regression analysis on the genotypes yielded a χ(2) of 35.5 with an odds ratio of 1.6 (CI = 1.3-1.8), P < 1 × 10(-5) . None of the clinical and metabolic parameters tested was associated with the genotype. In conclusion, we have confirmed the strong association between FTO and obesity, and shown that only AA homozygotes are predisposed to develop obesity while TT homozygotes might be protected. Finally, we found no association between rs9939609 and a number of obesity-related abnormalities.

An INSIG2 polymorphism affects glucose homeostasis in Sardinian obese children and adolescents.

Allelic variants of a single nucleotide polymorphism (SNP), rs7566605, located approximately 10 kb upstream of the INSIG2 gene have been found in association with body weight and with other clinical features related to obesity in some populations but not in others. Our objective was to test the association of this SNP in obese children and adolescents from the genetically isolated population of Sardinia. We tested the association of rs7566605 with body mass index (BMI) and with serum glucose and insulin concentrations and a surrogate measure of insulin resistance (HOMA-IR) in a cohort of 747 Sardinian obese children and adolescents. A case control analysis was performed using 548 ethnically-matched healthy controls. Allelic frequencies of the SNP were similar between patients and controls. Mean glucose and insulin concentration and mean HOMA-IR values were significantly higher in patients carrying the CC genotype than in the CG and GG carriers. In the patients with impaired fasting glucose (IFG) and/or impaired glucose tolerance (IGT), allele C was significantly more frequent than in controls. Although INSIG2 polymorphisms do not consistently associate with BMI, the observation of an association with glucose concentration would support a role for this gene in the metabolic complications of obesity.

Variants within the immunoregulatory CBLB gene are associated with multiple sclerosis.

A genome-wide association scan of approximately 6.6 million genotyped or imputed variants in 882 Sardinian individuals with multiple sclerosis (cases) and 872 controls suggested association of CBLB gene variants with disease, which was confirmed in 1,775 cases and 2,005 controls (rs9657904, overall P = 1.60 x 10(-10), OR = 1.40). CBLB encodes a negative regulator of adaptive immune responses, and mice lacking the ortholog are prone to experimental autoimmune encephalomyelitis, the animal model of multiple sclerosis.

Genetic loci linked to type 1 diabetes and multiple sclerosis families in Sardinia.

BACKGROUND: The Mediterranean island of Sardinia has a strikingly high incidence of the autoimmune disorders Type 1 Diabetes (T1D) and Multiple Sclerosis (MS). Furthermore, the two diseases tend to be co-inherited in the same individuals and in the same families. These observations suggest that some unknown autoimmunity variant with relevant effect size could be fairly common in this founder population and could be detected using linkage analysis. METHODS: To search for T1D and MS loci as well as any that predispose to both diseases, we performed a whole genome linkage scan, sequentially genotyping 593 microsatellite marker loci in 954 individuals distributed in 175 Sardinian families. In total, 413 patients were studied; 285 with T1D, 116 with MS and 12 with both disorders. Model-free linkage analysis was performed on the genotyped samples using the Kong and Cox logarithm of odds (LOD) score statistic. RESULTS: In T1D, aside from the HLA locus, we found four regions showing a lod-score > or =1; 1p31.1, 6q26, 10q21.2 and 22q11.22. In MS we found three regions showing a lod-score > or =1; 1q42.2, 18p11.21 and 20p12.3. In the combined T1D-MS scan for shared autoimmunity loci, four regions showed a LOD >1, including 6q26, 10q21.2, 20p12.3 and 22q11.22. When we typed more markers in these intervals we obtained suggestive evidence of linkage in the T1D scan at 10q21.2 (LOD = 2.1), in the MS scan at 1q42.2 (LOD = 2.5) and at 18p11.22 (LOD = 2.6). When all T1D and MS families were analysed jointly we obtained suggestive evidence in two regions: at 10q21.1 (LOD score = 2.3) and at 20p12.3 (LOD score = 2.5). CONCLUSION: This suggestive evidence of linkage with T1D, MS and both diseases indicates critical chromosome intervals to be followed up in downstream association studies.