Pharmacogenetics in Italy: current landscape and future prospects.

Pharmacogenetics investigates sequence of genes that affect drug response, enabling personalized medication. This approach reduces drug-induced adverse reactions and improves clinical effectiveness, making it a crucial consideration for personalized medical care. Numerous guidelines, drawn by global consortia and scientific organizations, codify genotype-driven administration for over 120 active substances. As the scientific community acknowledges the benefits of genotype-tailored therapy over traditionally agnostic drug administration, the push for its implementation into Italian healthcare system is gaining momentum. This evolution is influenced by several factors, including the improved access to patient genotypes, the sequencing costs decrease, the growing of large-scale genetic studies, the rising popularity of direct-to-consumer pharmacogenetic tests, and the continuous improvement of pharmacogenetic guidelines. Since EMA (European Medicines Agency) and AIFA (Italian Medicines Agency) provide genotype information on drug leaflet without clear and explicit clinical indications for gene testing, the regulation of pharmacogenetic testing is a pressing matter in Italy. In this manuscript, we have reviewed how to overcome the obstacles in implementing pharmacogenetic testing in the clinical practice of the Italian healthcare system. Our particular emphasis has been on germline testing, given the absence of well-defined national directives in contrast to somatic pharmacogenetics.

Functional annotation and investigation of the 10q24.33 melanoma risk locus identifies a common variant that influences transcriptional regulation of OBFC1.

The 10q24.33 locus is known to be associated with susceptibility to cutaneous malignant melanoma (CMM), but the mechanisms underlying this association have been not extensively investigated. We carried out an integrative genomic analysis of 10q24.33 using epigenomic annotations and in vitro reporter gene assays to identify regulatory variants. We found two putative functional single nucleotide polymorphisms (SNPs) in an enhancer and in the promoter of OBFC1, respectively, in neural crest and CMM cells, one, rs2995264, altering enhancer activity. The minor allele G of rs2995264 correlated with lower OBFC1 expression in 470 CMM tumors and was confirmed to increase the CMM risk in a cohort of 484 CMM cases and 1801 controls of Italian origin. Hi-C and chromosome conformation capture (3C) experiments showed the interaction between the enhancer-SNP region and the promoter of OBFC1 and an isogenic model characterized by CRISPR-Cas9 deletion of the enhancer-SNP region confirmed the potential regulatory effect of rs2995264 on OBFC1 transcription. Moreover, the presence of G-rs2995264 risk allele reduced the binding affinity of the transcription factor MEOX2. Biologic investigations showed significant cell viability upon depletion of OBFC1, specifically in CMM cells that were homozygous for the protective allele. Clinically, high levels of OBFC1 expression associated with histologically favorable CMM tumors. Finally, preliminary results suggested the potential effect of decreased OBFC1 expression on telomerase activity in tumorigenic conditions. Our results support the hypothesis that reduced expression of OBFC1 gene through functional heritable DNA variation can contribute to malignant transformation of normal melanocytes.

Integrative genomic analyses identify neuroblastoma risk genes involved in neuronal differentiation.

Genome-Wide Association Studies (GWAS) have been decisive in elucidating the genetic predisposition of neuroblastoma (NB). The majority of genetic variants identified in GWAS are found in non-coding regions, suggesting that they can be causative of pathogenic dysregulations of gene expression. Nonetheless, pinpointing the potential causal genes within implicated genetic loci remains a major challenge. In this study, we integrated NB GWAS and expression Quantitative Trait Loci (eQTL) data from adrenal gland to identify candidate genes impacting NB susceptibility. We found that ZMYM1, CBL, GSKIP and WDR81 expression was dysregulated by NB predisposing variants. We further investigated the functional role of the identified genes through computational analysis of RNA sequencing (RNA-seq) data from single-cell and whole-tissue samples of NB, neural crest, and adrenal gland tissues, as well as through in vitro differentiation assays in NB cell cultures. Our results indicate that dysregulation of ZMYM1, CBL, GSKIP, WDR81 may lead to malignant transformation by affecting early and late stages of normal program of neuronal differentiation. Our findings enhance the understanding of how specific genes contribute to NB pathogenesis by highlighting their influence on neuronal differentiation and emphasizing the impact of genetic risk variants on the regulation of genes involved in critical biological processes.

FGFR1 is a potential therapeutic target in neuroblastoma.

BACKGROUND: FGFR1 regulates cell-cell adhesion and extracellular matrix architecture and acts as oncogene in several cancers. Potential cancer driver mutations of FGFR1 occur in neuroblastoma (NB), a neural crest-derived pediatric tumor arising in sympathetic nervous system, but so far they have not been studied experimentally. We investigated the driver-oncogene role of FGFR1 and the implication of N546K mutation in therapy-resistance in NB cells. METHODS: Public datasets were used to predict the correlation of FGFR1 expression with NB clinical outcomes. Whole genome sequencing data of 19 paired diagnostic and relapse NB samples were used to find somatic mutations. In NB cell lines, silencing by short hairpin RNA and transient overexpression of FGFR1 were performed to evaluate the effect of the identified mutation by cell growth, invasion and cologenicity assays. HEK293, SHSY5Y and SKNBE2 were selected to investigate subcellular wild-type and mutated protein localization. FGFR1 inhibitor (AZD4547), alone or in combination with PI3K inhibitor (GDC0941), was used to rescue malignant phenotypes induced by overexpression of FGFR1 wild-type and mutated protein. RESULTS: High FGFR1 expression correlated with low relapse-free survival in two independent NB gene expression datasets. In addition, we found the somatic mutation N546K, the most recurrent point mutation of FGFR1 in all cancers and already reported in NB, in one out of 19 matched primary and recurrent tumors. Loss of FGFR1 function attenuated invasion and cologenicity in NB cells, whereas FGFR1 overexpression enhanced oncogenicity. The overexpression of FGFR1N546K protein showed a higher nuclear localization compared to wild-type protein and increased cellular invasion and cologenicity. Moreover, N546K mutation caused the failure in response to treatment with FGFR1 inhibitor by activation of ERK, STAT3 and AKT pathways. The combination of FGFR1 and PI3K pathway inhibitors was effective in reducing the invasive and colonigenic ability of cells overexpressing FGFR1 mutated protein. CONCLUSIONS: FGFR1 is an actionable driver oncogene in NB and a promising therapy may consist in targeting FGFR1 mutations in patients with therapy-resistant NB.

Neural crest-derived tumor neuroblastoma and melanoma share 1p13.2 as susceptibility locus that shows a long-range interaction with the SLC16A1 gene.

Neuroblastoma (NB) and malignant cutaneous melanoma (CMM) are neural crest cells (NCC)-derived tumors and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association studies (GWAS). We took a three-staged approach to conduct cross-disease meta-analysis of GWAS for NB and CMM (2101 NB cases and 4202 controls; 12 874 CMM cases and 23 203 controls) to identify shared loci. Findings were replicated in 1403 NB cases and 1403 controls of European ancestry and in 636 NB, 508 CMM cases and 2066 controls of Italian origin. We found a cross-association at locus 1p13.2 (rs2153977, odds ratio = 0.91, P = 5.36 × 10-8). We also detected a suggestive (P < 10-7) NB-CMM cross-association at 2q37.1 with opposite effect on cancer risk. Pathway analysis of 110 NB-CMM risk loci with P < 10-4 demonstrated enrichment of biological processes such as cell migration, cell cycle, metabolism and immune response, which are essential of human NCC development, underlying both tumors. In vitro and in silico analyses indicated that the rs2153977-T protective allele, located in an NB and CMM enhancer, decreased expression of SLC16A1 via long-range loop formation and altered a T-box protein binding site. Upon depletion of SLC16A1, we observed a decrease of cellular proliferation and invasion in both NB and CMM cell lines, suggesting its role as oncogene. This is the largest study to date examining pleiotropy across two NC cell-derived tumors identifying 1p13.2 as common susceptibility locus for NB and CMM risk. We demonstrate that combining genome-wide association studies results across cancers with same origins can identify new loci common to neuroblastoma and melanoma arising from tissues which originate from neural crest cells. Our results also show 1p13.2 confer risk to neuroblastoma and melanoma by regulating SLC16A1.

Association of PARP1 polymorphisms with response to chemotherapy in patients with high-risk neuroblastoma.

The genetic aetiology and the molecular mechanisms that characterize high-risk neuroblastoma are still little understood. The majority of high-risk neuroblastoma patients do not take advantage of current induction therapy. So far, one of the main reasons liable for cancer therapeutic failure is the acquisition of resistance to cytotoxic anticancer drugs, because of the DNA repair system of tumour cells. PARP1 is one of the main DNA damage sensors involved in the DNA repair system and genomic stability. We observed that high PARP1 mRNA level is associated with unfavourable prognosis in 3 public gene expression NB patients' datasets and in 20 neuroblastomas analysed by qRT-PCR. Among 4983 SNPs in PARP1, we selected two potential functional SNPs. We investigated the association of rs907187, in PARP1 promoter, and rs2048426 in non-coding region with response chemotherapy in 121 Italian patients with high-risk NB. Results showed that minor G allele of rs907187 associated with induction response of patients (P = .02) and with decrease PARP1 mRNA levels in NB cell line (P = .003). Furthermore, rs907187 was predicted to alter the binding site of E2F1 transcription factor. Specifically, allele G had low binding affinity with E2F1 whose expression positively correlates with PARP1 expression and associated with poor prognosis of patients with NB. By contrast, we did not find genetic association for the SNP rs2048426. These data reveal rs907187 as a novel potential risk variant associated with the failure of induction therapy for high-risk NB.

HIF-1 transcription activity: HIF1A driven response in normoxia and in hypoxia.

BACKGROUND: HIF1A (Hypoxia-Inducible-Factor 1A) expression in solid tumors is relevant to establish resistance to therapeutic approaches. The use of compounds direct against hypoxia signaling and HIF1A does not show clinical efficiency because of changeable oxygen concentrations in solid tumor areas. The identification of HIF1A targets expressed in both normoxia and hypoxia and of HIF1A/hypoxia signatures might meliorate the prognostic stratification and therapeutic successes in patients with high-risk solid tumors. METHODS: In this study, we conducted a combined analysis of RNA expression and DNA methylation of neuroblastoma cells silenced or unsilenced for HIF1A expression, grown in normoxia and hypoxia conditions. RESULTS: The analysis of pathways highlights HIF-1 (heterodimeric transcription factor 1) activity in normoxia in metabolic process and HIF-1 activity in hypoxia in neuronal differentiation process. HIF1A driven transcriptional response in hypoxia depends on epigenetic control at DNA methylation status of gene regulatory regions. Furthermore, low oxygen levels generate HIF1A-dependent or HIF1A-independent signatures, able to stratify patients according to risk categories. CONCLUSIONS: These findings may help to understand the molecular mechanisms by which low oxygen levels reshape gene signatures and provide new direction for hypoxia targeting in solid tumor.

Improving single nucleotide polymorphisms genotyping accuracy for dihydropyrimidine dehydrogenase testing in pharmacogenetics.

Fluoropyrimidines, crucial in cancer treatment, often cause toxicity concerns even at standard doses. Toxic accumulation of fluoropyrimidine metabolites, culminating in adverse effects, can stem from impaired dihydropyrimidine dehydrogenase (DPYD) enzymatic function. Emerging evidence underscores the role of single nucleotide polymorphisms (SNPs) in DPYD gene, capable of inducing DPYD activity deficiency. Consequently, DPYD genotyping's importance is on the rise in clinical practice before initiating fluoropyrimidine treatment. Although polymerase chain reaction (PCR) followed by Sanger sequencing (SS; PCR-SS) is a prevalent method for DPYD genotyping, it may encounter limitations. In this context, there is reported a case in which a routine PCR-SS approach for genotyping DPYD SNP rs55886062 failed in a proband of African descent. The Clinical Pharmacogenetics Implementation Consortium (CPIC) categorizes the guanine (G) allele of this SNP as non-functional. The enforcement of whole genome sequencing (WGS) approach led to the identification of two adenine (A) insertions near the PCR primers annealing regions in the proband, responsible for a sequence frameshift and a genotyping error for rs55886062. These SNPs (rs145228578, 1-97981199-T-TA and rs141050810, 1-97981622-G-GA) were extremely rare in non-Finnish Europeans (0.05%) but prevalent in African populations (16%). Although limited evidence was available for these SNPs, they were catalogued as benign variants in public databases. Notably, these two SNPs exhibited a high linkage disequilibrium [LD; squared correlation coefficient (R2) = 0.98]. These findings highlighted the importance to consider the prevalence of genetic variants within diverse ethnic populations when designing primers and probes for SNP genotyping in pharmacogenetic testing. This preventive measure is essential to avoid sequence frameshifts or primer misalignments arising from SNP occurrences in the genome, which can compromise PCR-SS and lead to genotyping failures. Furthermore, this case highlights the significance of exploring alternative genotyping approaches, like WGS, when confronted with challenges associated with conventional techniques.

Label-free liquid biopsy through the identification of tumor cells by machine learning-powered tomographic phase imaging flow cytometry.

Image-based identification of circulating tumor cells in microfluidic cytometry condition is one of the most challenging perspectives in the Liquid Biopsy scenario. Here we show a machine learning-powered tomographic phase imaging flow cytometry system capable to provide high-throughput 3D phase-contrast tomograms of each single cell. In fact, we show that discrimination of tumor cells against white blood cells is potentially achievable with the aid of artificial intelligence in a label-free flow-cyto-tomography method. We propose a hierarchical machine learning decision-maker, working on a set of features calculated from the 3D tomograms of the cells' refractive index. We prove that 3D morphological features are adequately distinctive to identify tumor cells versus the white blood cell background in the first stage and, moreover, in recognizing the tumor type at the second decision step. Proof-of-concept experiments are shown, in which two different tumor cell lines, namely neuroblastoma cancer cells and ovarian cancer cells, are used against monocytes. The reported results allow claiming the identification of tumor cells with a success rate higher than 97% and with an accuracy over 97% in discriminating between the two cancer cell types, thus opening in a near future the route to a new Liquid Biopsy tool for detecting and classifying circulating tumor cells in blood by stain-free method.

Phenotyping neuroblastoma cells through intelligent scrutiny of stain-free biomarkers in holographic flow cytometry.

To efficiently tackle certain tumor types, finding new biomarkers for rapid and complete phenotyping of cancer cells is highly demanded. This is especially the case for the most common pediatric solid tumor of the sympathetic nervous system, namely, neuroblastoma (NB). Liquid biopsy is in principle a very promising tool for this purpose, but usually enrichment and isolation of circulating tumor cells in such patients remain difficult due to the unavailability of universal NB cell-specific surface markers. Here, we show that rapid screening and phenotyping of NB cells through stain-free biomarkers supported by artificial intelligence is a viable route for liquid biopsy. We demonstrate the concept through a flow cytometry based on label-free holographic quantitative phase-contrast microscopy empowered by machine learning. In detail, we exploit a hierarchical decision scheme where at first level NB cells are classified from monocytes with 97.9% accuracy. Then we demonstrate that different phenotypes are discriminated within NB class. Indeed, for each cell classified as NB its belonging to one of four NB sub-populations (i.e., CHP212, SKNBE2, SHSY5Y, and SKNSH) is evaluated thus achieving accuracy in the range 73.6%-89.1%. The achieved results solve the realistic problem related to the identification circulating tumor cell, i.e., the possibility to recognize and detect tumor cells morphologically similar to blood cells, which is the core issue in liquid biopsy based on stain-free microscopy. The presented approach operates at lab-on-chip scale and emulates real-world scenarios, thus representing a future route for liquid biopsy by exploiting intelligent biomedical imaging.

Somatic Mutations Enriched in Cis-Regulatory Elements Affect Genes Involved in Embryonic Development and Immune System Response in Neuroblastoma.

Noncoding cis-regulatory variants have gained interest as cancer drivers, yet progress in understanding their significance is hindered by the numerous challenges and limitations of variant prioritization. To overcome these limitations, we focused on active cis-regulatory elements (aCRE) to design a customized panel for the deep sequencing of 56 neuroblastoma tumor and normal DNA sample pairs. To search for driver mutations, aCREs were defined by reanalysis of H3K27ac chromatin immunoprecipitation sequencing peaks in 25 neuroblastoma cell lines. These regulatory genomic regions were tested for an excess of somatic mutations and assessed for statistical significance using a global approach that accounted for chromatin accessibility and replication timing. Additional validation was provided by whole genome sequence analysis of 151 neuroblastomas. Analysis of HiC data determined the presence of candidate target genes interacting with mutated regions. An excess of somatic mutations in aCREs of diverse genes were identified, including IPO7, HAND2, and ARID3A. CRISPR-Cas9 editing was utilized to assess the functional consequences of mutations in the IPO7-aCRE. Patients with noncoding mutations in aCREs showed inferior overall and event-free survival independent of age at diagnosis, stage, risk stratification, and MYCN status. Expression of aCRE-interacting genes correlated strongly with negative prognostic markers and low survival rates. Moreover, a convergence between the biological functions of aCRE target genes and transcription factors with mutated binding motifs was associated with embryonic development and immune system response. Overall, this strategy enabled the identification of somatic mutations in regulatory elements that collectively promote neuroblastoma tumorigenesis. SIGNIFICANCE: Assessment of noncoding cis-regulatory variants and long-range interaction data highlight the combined effect of somatic mutations in regulatory elements in driving neuroblastoma.

Stain-free identification of cell nuclei using tomographic phase microscopy in flow cytometry.

Quantitative Phase Imaging (QPI) has gained popularity in bioimaging because it can avoid the need for cell staining, which in some cases is difficult or impossible. However, as a result, QPI does not provide labelling of various specific intracellular structures. Here we show a novel computational segmentation method based on statistical inference that makes it possible for QPI techniques to identify the cell nucleus. We demonstrate the approach with refractive index tomograms of stain-free cells reconstructed through the tomographic phase microscopy in flow cytometry mode. In particular, by means of numerical simulations and two cancer cell lines, we demonstrate that the nucleus can be accurately distinguished within the stain-free tomograms. We show that our experimental results are consistent with confocal fluorescence microscopy (FM) data and microfluidic cytofluorimeter outputs. This is a significant step towards extracting specific three-dimensional intracellular structures directly from the phase-contrast data in a typical flow cytometry configuration.

Genetic Predisposition to Solid Pediatric Cancers.

Progresses over the past years have extensively improved our capacity to use genome-scale analyses-including high-density genotyping and exome and genome sequencing-to identify the genetic basis of pediatric tumors. In particular, exome sequencing has contributed to the evidence that about 10% of children and adolescents with tumors have germline genetic variants associated with cancer predisposition. In this review, we provide an overview of genetic variations predisposing to solid pediatric tumors (medulloblastoma, ependymoma, astrocytoma, neuroblastoma, retinoblastoma, Wilms tumor, osteosarcoma, rhabdomyosarcoma, and Ewing sarcoma) and outline the biological processes affected by the involved mutated genes. A careful description of the genetic basis underlying a large number of syndromes associated with an increased risk of pediatric cancer is also reported. We place particular emphasis on the emerging view that interactions between germline and somatic alterations are a key determinant of cancer development. We propose future research directions, which focus on the biological function of pediatric risk alleles and on the potential links between the germline genome and somatic changes. Finally, the importance of developing new molecular diagnostic tests including all the identified risk germline mutations and of considering the genetic predisposition in screening tests and novel therapies is emphasized.

Functional characterization of full-length BARD1 strengthens its role as a tumor suppressor in neuroblastoma.

BARD1 is associated with the development of high-risk neuroblastoma patients. Particularly, the expression of full length (FL) isoform, FL BARD1, correlates to high-risk neuroblastoma development and its inhibition is sufficient to induce neuroblastoma cells towards a worst phenotype. Here we have investigated the mechanisms of FL BARD1 in neuroblastoma cell lines depleted for FL BARD1 expression. We have shown that FL BARD1 expression protects the cells from spontaneous DNA damage and from damage accumulated after irradiation. We demonstrated a role for FL BARD1 as tumor suppressor to prevent unscheduled mitotic entry of DNA damaged cells and to lead to death cells that have bypassed cell cycle checkpoints. FL BARD1-depleted cells that have survived to checkpoints acquire features of aggressiveness. Overall, our results show that FL BARD1 may defend cells against cancer and prevent malignant transformation of cells.

Transcription Factors Involved in Tumorigenesis Are Over-Represented in Mutated Active DNA-Binding Sites in Neuroblastoma.

The contribution of coding mutations to oncogenesis has been largely clarified, whereas little is known about somatic mutations in noncoding DNA and their role in driving tumors remains controversial. Here, we used an alternative approach to interpret the functional significance of noncoding somatic mutations in promoting tumorigenesis. Noncoding somatic mutations of 151 neuroblastomas were integrated with ENCODE data to locate somatic mutations in regulatory elements specifically active in neuroblastoma cells, nonspecifically active in neuroblastoma cells, and nonactive. Within these types of elements, transcription factors (TF) were identified whose binding sites were enriched or depleted in mutations. For these TFs, a gene expression signature was built to assess their implication in neuroblastoma. DNA- and RNA-sequencing data were integrated to assess the effects of those mutations on mRNA levels. The pathogenicity of mutations was significantly higher in transcription factor binding site (TFBS) of regulatory elements specifically active in neuroblastoma cells, as compared with the others. Within these elements, there were 18 over-represented TFs involved mainly in cell-cycle phase transitions and 15 under-represented TFs primarily regulating cell differentiation. A gene expression signature based on over-represented TFs correlated with poor survival and unfavorable prognostic markers. Moreover, recurrent mutations in TFBS of over-represented TFs such as EZH2 affected MCF2L and ADP-ribosylhydrolase like 1 expression, among the others. We propose a novel approach to study the involvement of regulatory variants in neuroblastoma that could be extended to other cancers and provide further evidence that alterations of gene expression may have relevant effects in neuroblastoma development. SIGNIFICANCE: These findings propose a novel approach to study regulatory variants in neuroblastoma and suggest that noncoding somatic mutations have relevant implications in neuroblastoma development.