regioneReloaded: evaluating the association of multiple genomic region sets.

MOTIVATION: Next-generation sequencing methods continue improving the annotation of genomes in part by determining the distribution of features such as epigenetic marks. Evaluating and interpreting the association between genomic regions and their features has become a common and challenging analysis in genomic and epigenomic studies. RESULTS: With regioneR we provided an R package allowing to assess the statistical significance of pairwise associations between genomic region sets using permutation tests. We now present the R package regioneReloaded that builds upon regioneR's statistical foundation and extends the functionality for the simultaneous analysis and visualization of the associations between multiple genomic region sets. Thus, we provide a novel discovery tool for the identification of significant associations that warrant to be tested for functional interdependence. AVAILABILITY AND IMPLEMENTATION: regioneReloaded is an R package released under an Artistic-2.0 License. The source code and documentation are freely available through Bioconductor: http://www.bioconductor.org/packages/regioneReloaded.

Screening of CNVs using NGS data improves mutation detection yield and decreases costs in genetic testing for hereditary cancer.

INTRODUCTION: Germline CNVs are important contributors to hereditary cancer. In genetic diagnostics, multiplex ligation-dependent probe amplification (MLPA) is commonly used to identify them. However, MLPA is time-consuming and expensive if applied to many genes, hence many routine laboratories test only a subset of genes of interest. METHODS AND RESULTS: We evaluated a next-generation sequencing (NGS)-based CNV detection tool (DECoN) as first-tier screening to decrease costs and turnaround time and expand CNV analysis to all genes of clinical interest in our diagnostics routine. We used DECoN in a retrospective cohort of 1860 patients where a limited number of genes were previously analysed by MLPA, and in a prospective cohort of 2041 patients, without MLPA analysis. In the retrospective cohort, 6 new CNVs were identified and confirmed by MLPA. In the prospective cohort, 19 CNVs were identified and confirmed by MLPA, 8 of these would have been lost in our previous MLPA-restricted detection strategy. Also, the number of genes tested by MLPA across all samples decreased by 93.0% in the prospective cohort. CONCLUSION: Including an in silico germline NGS CNV detection tool improved our genetic diagnostics strategy in hereditary cancer, both increasing the number of CNVs detected and reducing turnaround time and costs.

CNVfilteR: an R/Bioconductor package to identify false positives produced by germline NGS CNV detection tools.

SUMMARY: Germline copy-number variants (CNVs) are relevant mutations for multiple genetics fields, such as the study of hereditary diseases. However, available benchmarks show that all next-generation sequencing (NGS) CNV calling tools produce false positives. We developed CNVfilteR, an R package that uses the single-nucleotide variant calls usually obtained in germline NGS pipelines to identify those false positives. The package can detect both false deletions and false duplications. We evaluated CNVfilteR performance on callsets generated by 13 CNV calling tools on three whole-genome sequencing and 541 panel samples, showing a decrease of up to 44.8% in false positives and consistent F1-score increase. Using CNVfilteR to detect false-positive calls can improve the overall performance of existing CNV calling pipelines. AVAILABILITY AND IMPLEMENTATION: CNVfilteR is released under Artistic-2.0 License. Source code and documentation are freely available at Bioconductor (http://www.bioconductor.org/packages/CNVfilteR). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Chromosomal translocations inactivating CDKN2A support a single path for malignant peripheral nerve sheath tumor initiation.

Malignant peripheral nerve sheath tumors (MPNST) are aggressive soft tissue sarcomas with poor prognosis, developing either sporadically or in persons with neurofibromatosis type 1 (NF1). Loss of CDKN2A/B is an important early event in MPNST progression. However, many reported MPNSTs exhibit partial or no inactivation of CDKN2A/B, raising the question of whether there is more than one molecular path for MPNST initiation. We present here a comprehensive genomic analysis of MPNST cell lines and tumors to explore in depth the status of CDKN2A. After accounting for CDKN2A deletions and point mutations, we uncovered a previously unnoticed high frequency of chromosomal translocations involving CDKN2A in both MPNST cell lines and primary tumors. Most identified translocation breakpoints were validated by PCR amplification and Sanger sequencing. Many breakpoints clustered in an intronic 500 bp hotspot region adjacent to CDKN2A exon 2. We demonstrate the bi-allelic inactivation of CDKN2A in all tumors (n = 15) and cell lines (n = 8) analyzed, supporting a single molecular path for MPNST initiation in both sporadic and NF1-related MPNSTs. This general CDKN2A inactivation in MPNSTs has implications for MPNST diagnostics and treatment. Our findings might be relevant for other tumor types with high frequencies of CDKN2A inactivation.

Mutational spectrum by phenotype: panel-based NGS testing of patients with clinical suspicion of RASopathy and children with multiple café-au-lait macules.

Children with neurofibromatosis type 1 (NF1) may exhibit an incomplete clinical presentation, making difficult to reach a clinical diagnosis. A phenotypic overlap may exist in children with other RASopathies or with other genetic conditions if only multiple café-au-lait macules (CALMs) are present. The syndromes that can converge in these inconclusive phenotypes have different clinical courses. In this context, an early genetic testing has been proposed to be clinically useful to manage these patients. We present the validation and implementation into diagnostics of a custom NGS panel (I2HCP, ICO-IMPPC Hereditary Cancer Panel) for testing patients with a clinical suspicion of a RASopathy (n = 48) and children presenting multiple CALMs (n = 102). We describe the mutational spectrum and the detection rates identified in these two groups of individuals. We identified pathogenic variants in 21 out of 48 patients with clinical suspicion of RASopathy, with mutations in NF1 accounting for 10% of cases. Furthermore, we identified pathogenic mutations mainly in the NF1 gene, but also in SPRED1, in more than 50% of children with multiple CALMs, exhibiting an NF1 mutational spectrum different from a group of clinically diagnosed NF1 patients (n = 80). An NGS panel strategy for the genetic testing of these two phenotype-defined groups outperforms previous strategies.

Analysis of intratumor heterogeneity in Neurofibromatosis type 1 plexiform neurofibromas and neurofibromas with atypical features: Correlating histological and genomic findings.

Plexiform neurofibromas (PNFs) are benign peripheral nerve sheath tumors involving large nerves present in 30%-50% Neurofibromatosis type 1 (NF1) patients. Atypical neurofibromas (ANF) are distinct nodular lesions with atypical features on histology that arise from PNFs. The risk and timeline of malignant transformation in ANF is difficult to assess. A recent NIH workshop has stratified ANFs and separated a subgroup with multiple atypical features and higher risk of malignant transformation termed atypical neurofibromatous neoplasms with uncertain biological potential (ANNUBP). We performed an analysis of intratumor heterogeneity on eight PNFs to link histological and genomic findings. Tumors were homogeneous although histological and molecular heterogeneity was identified. All tumors were 2n, almost mutation-free and had a clonal NF1(-/-) origin. Two ANFs from the same patient showed atypical features on histology and deletions of CDKN2A/B. One of the ANFs exhibited different areas in which the degree of histological atypia correlated with the heterozygous or homozygous loss of the CDKN2A/B loci. CDKN2A/B deletions in different areas originated independently. Results may indicate that loss of a single CDKN2A/B copy in NF1(-/-) cells is sufficient to start ANF development and that total inactivation of both copies of CDKN2A/B is necessary to form an ANNUBP.

karyoploteR: an R/Bioconductor package to plot customizable genomes displaying arbitrary data.

MOTIVATION: Data visualization is a crucial tool for data exploration, analysis and interpretation. For the visualization of genomic data there lacks a tool to create customizable non-circular plots of whole genomes from any species. RESULTS: We have developed karyoploteR, an R/Bioconductor package to create linear chromosomal representations of any genome with genomic annotations and experimental data plotted along them. Plot creation process is inspired in R base graphics, with a main function creating karyoplots with no data and multiple additional functions, including custom functions written by the end-user, adding data and other graphical elements. This approach allows the creation of highly customizable plots from arbitrary data with complete freedom on data positioning and representation. AVAILABILITY AND IMPLEMENTATION: karyoploteR is released under Artistic-2.0 License. Source code and documentation are freely available through Bioconductor (http://www.bioconductor.org/packages/karyoploteR) and at the examples and tutorial page at https://bernatgel.github.io/karyoploter_tutorial. CONTACT: bgel@igtp.cat.

regioneR: an R/Bioconductor package for the association analysis of genomic regions based on permutation tests.

MOTIVATION: Statistically assessing the relation between a set of genomic regions and other genomic features is a common challenging task in genomic and epigenomic analyses. Randomization based approaches implicitly take into account the complexity of the genome without the need of assuming an underlying statistical model. SUMMARY: regioneR is an R package that implements a permutation test framework specifically designed to work with genomic regions. In addition to the predefined randomization and evaluation strategies, regioneR is fully customizable allowing the use of custom strategies to adapt it to specific questions. Finally, it also implements a novel function to evaluate the local specificity of the detected association. AVAILABILITY AND IMPLEMENTATION: regioneR is an R package released under Artistic-2.0 License. The source code and documents are freely available through Bioconductor (http://www.bioconductor.org/packages/regioneR). CONTACT: rmalinverni@carrerasresearch.org.

Unbalancing cAMP and Ras/MAPK pathways as a therapeutic strategy for cutaneous neurofibromas.

Cutaneous neurofibromas (cNFs) are benign Schwann cell (SC) tumors arising from subepidermal glia. Individuals with neurofibromatosis type 1 (NF1) may develop thousands of cNFs, which greatly affect their quality of life. cNF growth is driven by the proliferation of NF1-/- SCs and their interaction with the NF1+/- microenvironment. We analyzed the crosstalk between human cNF-derived SCs and fibroblasts (FBs), identifying an expression signature specific to the SC-FB interaction. We validated the secretion of proteins involved in immune cell migration, suggesting a role of SC-FB crosstalk in immune cell recruitment. The signature also captured components of developmental signaling pathways, including the cAMP elevator G protein-coupled receptor 68 (GPR68). Activation of Gpr68 by ogerin in combination with the MEK inhibitor (MEKi) selumetinib reduced viability and induced differentiation and death of human cNF-derived primary SCs, a result corroborated using an induced pluripotent stem cell-derived 3D neurofibromasphere model. Similar results were obtained using other Gpr68 activators or cAMP analogs/adenylyl cyclase activators in combination with selumetinib. Interestingly, whereas primary SC cultures restarted their proliferation after treatment with selumetinib alone was stopped, the combination of ogerin-selumetinib elicited a permanent halt on SC expansion that persisted after drug removal. These results indicate that unbalancing the Ras and cAMP pathways by combining MEKi and cAMP elevators could be used as a potential treatment for cNFs.

Expanding a precision medicine platform for malignant peripheral nerve sheath tumors: New patient-derived orthotopic xenografts, cell lines and tumor entities.

Malignant peripheral nerve sheath tumors (MPNSTs) are aggressive soft-tissue sarcomas with a poor survival rate, presenting either sporadically or in the context of neurofibromatosis type 1 (NF1). The histological diagnosis of MPNSTs can be challenging, with different tumors exhibiting great histological and marker expression overlap. This heterogeneity could be partly responsible for the observed disparity in treatment response due to the inherent diversity of the preclinical models used. For several years, our group has been generating a large patient-derived orthotopic xenograft (PDOX) MPNST platform for identifying new precision medicine treatments. Herein, we describe the expansion of this platform using six primary tumors clinically diagnosed as MPNSTs, from which we obtained six additional PDOX mouse models and three cell lines, thus generating three pairs of in vitro-in vivo models. We extensively characterized these tumors and derived preclinical models, including genomic, epigenomic, and histological analyses. Tumors were reclassified after these analyses: three remained as MPNSTs (two being classic MPNSTs), one was a melanoma, another was a neurotrophic tyrosine receptor kinase (NTRK)-rearranged spindle cell neoplasm, and, finally, the last was an unclassifiable tumor bearing neurofibromin-2 (NF2) inactivation, a neuroblastoma RAS viral oncogene homolog (NRAS) oncogenic mutation, and a SWI/SNF-related matrix-associated actin-dependent regulator of chromatin (SMARCA4) heterozygous truncated variant. New cell lines and PDOXs faithfully recapitulated histology, marker expression, and genomic characteristics of the primary tumors. The diversity in tumor identity and their specific associated genomic alterations impacted treatment responses obtained when we used the new cell lines for testing compounds against known altered pathways in MPNSTs. In summary, we present here an extension of our MPNST precision medicine platform, with new PDOXs and cell lines, including tumor entities confounded as MPNSTs in a real clinical scenario. This platform may constitute a useful tool for obtaining correct preclinical information to guide MPNST clinical trials.

SpadaHC: a database to improve the classification of variants in hereditary cancer genes in the Spanish population.

Accurate classification of genetic variants is crucial for clinical decision-making in hereditary cancer. In Spain, genetic diagnostic laboratories have traditionally approached this task independently due to the lack of a dedicated resource. Here we present SpadaHC, a web-based database for sharing variants in hereditary cancer genes in the Spanish population. SpadaHC is implemented using a three-tier architecture consisting of a relational database, a web tool and a bioinformatics pipeline. Contributing laboratories can share variant classifications and variants from individuals in Variant Calling Format (VCF) format. The platform supports open and restricted access, flexible dataset submissions, automatic pseudo-anonymization, VCF quality control, variant normalization and liftover between genome builds. Users can flexibly explore and search data, receive automatic discrepancy notifications and access SpadaHC population frequencies based on many criteria. In February 2024, SpadaHC included 18 laboratory members, storing 1.17 million variants from 4306 patients and 16 343 laboratory classifications. In the first analysis of the shared data, we identified 84 genetic variants with clinically relevant discrepancies in their classifications and addressed them through a three-phase resolution strategy. This work highlights the importance of data sharing to promote consistency in variant classifications among laboratories, so patients and family members can benefit from more accurate clinical management. Database URL: https://spadahc.ciberisciii.es/.

Generation of human iPSC-derived neurofibromaspheres for in vitro and in vivo uses.

Neurofibromas are benign peripheral nervous system tumors associated with neurofibromatosis type 1, which originate from NF1(-/-) Schwann cell precursors. We describe a protocol to generate neurofibromaspheres by differentiating NF1(-/-) Schwann cells from induced pluripotent stem cells and combining them with neurofibroma primary fibroblasts. We also describe the development of neurofibroma-like tumors when neurofibromaspheres are engrafted in the sciatic nerve of nude mice. This model constitutes a versatile platform for drug screening and the study of neurofibroma biology. For complete details on the use and execution of this protocol, please refer to Mazuelas et al. (2022).1.

Deep genomic analysis of malignant peripheral nerve sheath tumor cell lines challenges current malignant peripheral nerve sheath tumor diagnosis.

Malignant peripheral nerve sheath tumors (MPNSTs) are soft-tissue sarcomas of the peripheral nervous system that develop either sporadically or in the context of neurofibromatosis type 1 (NF1). MPNST diagnosis can be challenging and treatment outcomes are poor. We present here a resource consisting of the genomic characterization of 9 widely used human MPNST cell lines for their use in translational research. NF1-related cell lines recapitulated primary MPNST copy number profiles, exhibited NF1, CDKN2A, and SUZ12/EED tumor suppressor gene (TSG) inactivation, and presented no gain-of-function mutations. In contrast, sporadic cell lines collectively displayed different TSG inactivation patterns and presented kinase-activating mutations, fusion genes, altered mutational frequencies and COSMIC signatures, and different methylome-based classifications. Cell lines re-classified as melanomas and other sarcomas exhibited a different drug-treatment response. Deep genomic analysis, methylome-based classification, and cell-identity marker expression, challenged the identity of common MPNST cell lines, opening an opportunity to revise MPNST differential diagnosis.

Dasty3, a WEB framework for DAS.

MOTIVATION: Dasty3 is a highly interactive and extensible Web-based framework. It provides a rich Application Programming Interface upon which it is possible to develop specialized clients capable of retrieving information from DAS sources as well as from data providers not using the DAS protocol. Dasty3 provides significant improvements on previous Web-based frameworks and is implemented using the 1.6 DAS specification. AVAILABILITY: Dasty3 is an open-source tool freely available at http://www.ebi.ac.uk/dasty/ under the terms of the GNU General public license. Source and documentation can be found at http://code.google.com/p/dasty/. CONTACT: hhe@ebi.ac.uk.

Modeling iPSC-derived human neurofibroma-like tumors in mice uncovers the heterogeneity of Schwann cells within plexiform neurofibromas.

Plexiform neurofibromas (pNFs) are developmental tumors that appear in neurofibromatosis type 1 individuals, constituting a major source of morbidity and potentially transforming into a highly metastatic sarcoma (MPNST). pNFs arise after NF1 inactivation in a cell of the neural crest (NC)-Schwann cell (SC) lineage. Here, we develop an iPSC-based NC-SC in vitro differentiation system and construct a lineage expression roadmap for the analysis of different 2D and 3D NF models. The best model consists of generating heterotypic spheroids (neurofibromaspheres) composed of iPSC-derived differentiating NF1(-/-) SCs and NF1(+/-) pNF-derived fibroblasts (Fbs). Neurofibromaspheres form by maintaining highly proliferative NF1(-/-) cells committed to the NC-SC axis due to SC-SC and SC-Fb interactions, resulting in SC linage cells at different maturation points. Upon engraftment on the mouse sciatic nerve, neurofibromaspheres consistently generate human NF-like tumors. Analysis of expression roadmap genes in human pNF single-cell RNA-seq data uncovers the presence of SC subpopulations at distinct differentiation states.

Regulation of JAK2 by miR-135a: prognostic impact in classic Hodgkin lymphoma.

The behavior of classic Hodgkin lymphoma (cHL) is determined by both the intrinsic features of the tumor cells and the characteristics of the microenvironment, making the analysis of entire lymph nodes an effective approach to understanding the disease. We examined the influence of our previously reported 25-microRNA signature for cHL on clinical outcome in 89 homogeneously treated cHL patients with a median follow-up of 80 months. Patients with low miR-135a expression had a higher probability of relapse (P = .04) and a shorter disease-free survival (P = .02). Functional analysis of cHL cell lines showed that mature miR-135a levels increased after pre-miR-135a transfection, causing apoptosis and decreased cell growth. Target analysis showed a direct regulation by miR-135a of JAK2, a cytoplasmic tyrosine kinase involved in a specific subset of cytokine receptor signaling pathways. miR-135a-mediated JAK2 down-regulation led to decreased mRNA and protein levels of the antiapoptotic gene Bcl-xL, suggesting a role for Bcl-xL in miR-135a/JAK2-mediated apoptosis. Our findings confirm the critical role of miR-135a in the survival of cHL cells and in the prognosis of cHL patients, indicating that novel treatment approaches targeting miR-135a may potentially benefit these patients.

Single nucleotide polymorphisms in tobacco metabolism and DNA repair genes and prognosis in resected non-small-cell lung cancer.

BACKGROUND: If tobacco-related carcinogens are not inactivated or extruded from the cell, they can damage the DNA. Single nucleotide polymorphisms (SNPs) in genes involved in tobacco metabolism, DNA repair, and multidrug resistance have been related to lung cancer susceptibility. We examined 13 SNPs in 10 of these genes and correlated the results with time to progression (TTP) and overall survival (OS) in 71 smoker or former smoker patients with resected non-small-cell lung cancer (NSCLC). MATERIALS AND METHODS: DNA was obtained from paraffin-embedded tumor. SNP analysis of the candidate genes was performed by allelic discrimination assay. Log-rank test, Kaplan-Meier plots, and Cox multivariate analysis were used to evaluate the association of TTP and survival with the SNPs evaluated. RESULTS: Patients with wild-type (wt) XPC rs2228001, wt CYP2C8 rs10509681, or non-wt NAT2 rs1799930 had a longer TTP. Patients with wt ERCC1 showed a nonsignificant trend towards longer TTP. No other relation between SNPs and TTP were observed. Patients harboring at least two unfavorable genotypes in these four genes had a shorter TTP and OS than patients with either one or no unfavorable genotypes. In the multivariate analysis, non-wt XPC rs2228001 and the presence of at least two unfavorable genotypes emerged as independent markers for shorter TTP. CONCLUSIONS: SNPs in tobacco metabolism and DNA repair genes may influence the clinical outcome of resected NSCLC.

Prognostic implications of miR-16 expression levels in resected non-small-cell lung cancer.

BACKGROUND: MicroRNAs are novel regulators of gene expression that are linked to the main oncogene networks, including the p53 pathway. p53 regulates the maturation process of miR-16 and miR-143. We analyzed the role as prognostic markers of miR-16 and miR-143 in 70 non-small-cell lung cancer (NSCLC) patients. METHODS: MicroRNAs were analyzed by TaqMan MicroRNA assays. Disease-free survival (DFS) and overall survival (OS) were examined using Kaplan-Meier curves with log-rank tests and the Cox proportional hazard model. RESULTS: When patients were classified in three groups according to their miR-16 expression levels, those with normal levels had the best outcome while those with high levels had the worst. DFS was 22.4 months for patients with high levels, 71.8 months for those with normal levels, and 55.8 months for those with low levels (P = 0.05). OS was 23.9 months for patients with high levels, 97.6 months for those with normal levels, and 63.5 months for those with low levels (P < 0.001). In the multivariate analyses, high miR-16 levels emerged as an independent prognostic factor for poor DFS (P = 0.001) and OS (<0.001). CONCLUSIONS: Our results provide the first hints that miR-16 levels in tumor samples may be a prognostic marker in NSCLC.

Using antisense oligonucleotides for the physiological modulation of the alternative splicing of NF1 exon 23a during PC12 neuronal differentiation.

Neurofibromatosis Type 1 (NF1) is a genetic condition affecting approximately 1:3500 persons worldwide. The NF1 gene codes for neurofibromin protein, a GTPase activating protein (GAP) and a negative regulator of RAS. The NF1 gene undergoes alternative splicing of exon 23a (E23a) that codes for 21 amino acids placed at the center of the GAP related domain (GRD). E23a-containing type II neurofibromin exhibits a weaker Ras-GAP activity compared to E23a-less type I isoform. Exon E23a has been related with the cognitive impairment present in NF1 individuals. We designed antisense Phosphorodiamidate Morpholino Oligomers (PMOs) to modulate E23a alternative splicing at physiological conditions of gene expression and tested their impact during PC12 cell line neuronal differentiation. Results show that any dynamic modification of the natural ratio between type I and type II isoforms disturbed neuronal differentiation, altering the proper formation of neurites and deregulating both the MAPK/ERK and cAMP/PKA signaling pathways. Our results suggest an opposite regulation of these pathways by neurofibromin and the possible existence of a feedback loop sensing neurofibromin-related signaling. The present work illustrates the utility of PMOs to study alternative splicing that could be applied to other alternatively spliced genes in vitro and in vivo.

Common variants in NLRP2 and NLRP3 genes are strong prognostic factors for the outcome of HLA-identical sibling allogeneic stem cell transplantation.

The inflammasomes are macromolecular cytosolic complexes involved in the production of interleukin-1beta (IL-1beta) and IL-18 in response to several pathogen-derived stimuli. Such interleukins have been implicated in the origin of severe allogeneic stem cell transplant (allo-SCT) complications. We analyzed the relationship between the interindividual variability in inflammasome protein-encoding genes in donors and patients and clinical outcome after allo-SCT. Fourteen common genetic variants in 5 genes of the inflammasome, namely, NLRP1, NLRP2, NLRP3, CARD8, and CASP5, were genotyped in 133 human leukocyte antigen-identical sibling pairs undergoing allo-SCT. In the multivariate analysis, donor variants in NLRP2 and NLRP3 were the most important prognostic factors for the clinical outcome after allo-SCT. Thus, donor TT genotype at rs10925027 in NLRP3 was associated with disease relapse (odds ratio (OR) = 6.3, P = 1 x 10(-7)), and donor GG genotype at rs1043684 in NLRP2 was associated with nonrelapse mortality (OR = 4.4, P = 6 x 10(-4)) and overall survival (OR = 3.1, P = .001). In addition, patient AA genotype at rs5862 in NLRP1 was associated with nonrelapse mortality (OR = 2.8, P = .005) and overall survival (OR = 2.0, P = .009). These results suggest that inflammasome genetic variants are important prognostic factors for the outcome of allo-SCT. If validated in larger studies, including unrelated allo-SCT, NLRPs genotype would become an important factor in donor selection.