labelSeg: segment annotation for tumor copy number alteration profiles.

Somatic copy number alterations (SCNAs) are a predominant type of oncogenomic alterations that affect a large proportion of the genome in the majority of cancer samples. Current technologies allow high-throughput measurement of such copy number aberrations, generating results consisting of frequently large sets of SCNA segments. However, the automated annotation and integration of such data are particularly challenging because the measured signals reflect biased, relative copy number ratios. In this study, we introduce labelSeg, an algorithm designed for rapid and accurate annotation of CNA segments, with the aim of enhancing the interpretation of tumor SCNA profiles. Leveraging density-based clustering and exploiting the length-amplitude relationships of SCNA, our algorithm proficiently identifies distinct relative copy number states from individual segment profiles. Its compatibility with most CNA measurement platforms makes it suitable for large-scale integrative data analysis. We confirmed its performance on both simulated and sample-derived data from The Cancer Genome Atlas reference dataset, and we demonstrated its utility in integrating heterogeneous segment profiles from different data sources and measurement platforms. Our comparative and integrative analysis revealed common SCNA patterns in cancer and protein-coding genes with a strong correlation between SCNA and messenger RNA expression, promoting the investigation into the role of SCNA in cancer development.

Leveraging European infrastructures to access 1 million human genomes by 2022.

Human genomics is undergoing a step change from being a predominantly research-driven activity to one driven through health care as many countries in Europe now have nascent precision medicine programmes. To maximize the value of the genomic data generated, these data will need to be shared between institutions and across countries. In recognition of this challenge, 21 European countries recently signed a declaration to transnationally share data on at least 1 million human genomes by 2022. In this Roadmap, we identify the challenges of data sharing across borders and demonstrate that European research infrastructures are well-positioned to support the rapid implementation of widespread genomic data access.

Author Correction: Leveraging European infrastructures to access 1 million human genomes by 2022.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Beacon v2 and Beacon networks: A "lingua franca" for federated data discovery in biomedical genomics, and beyond.

Beacon is a basic data discovery protocol issued by the Global Alliance for Genomics and Health (GA4GH). The main goal addressed by version 1 of the Beacon protocol was to test the feasibility of broadly sharing human genomic data, through providing simple "yes" or "no" responses to queries about the presence of a given variant in datasets hosted by Beacon providers. The popularity of this concept has fostered the design of a version 2, that better serves real-world requirements and addresses the needs of clinical genomics research and healthcare, as assessed by several contributing projects and organizations. Particularly, rare disease genetics and cancer research will benefit from new case level and genomic variant level requests and the enabling of richer phenotype and clinical queries as well as support for fuzzy searches. Beacon is designed as a "lingua franca" to bridge data collections hosted in software solutions with different and rich interfaces. Beacon version 2 works alongside popular standards like Phenopackets, OMOP, or FHIR, allowing implementing consortia to return matches in beacon responses and provide a handover to their preferred data exchange format. The protocol is being explored by other research domains and is being tested in several international projects.

Minimum error calibration and normalization for genomic copy number analysis.

BACKGROUND: Copy number variations (CNV) are regional deviations from the normal autosomal bi-allelic DNA content. While germline CNVs are a major contributor to genomic syndromes and inherited diseases, the majority of cancers accumulate extensive "somatic" CNV (sCNV or CNA) during the process of oncogenetic transformation and progression. While specific sCNV have closely been associated with tumorigenesis, intriguingly many neoplasias exhibit recurrent sCNV patterns beyond the involvement of a few cancer driver genes. Currently, CNV profiles of tumor samples are generated using genomic micro-arrays or high-throughput DNA sequencing. Regardless of the underlying technology, genomic copy number data is derived from the relative assessment and integration of multiple signals, with the data generation process being prone to contamination from several sources. Estimated copy number values have no absolute or strictly linear correlation to their corresponding DNA levels, and the extent of deviation differs between sample profiles, which poses a great challenge for data integration and comparison in large scale genome analysis. RESULTS: In this study, we present a novel method named "Minimum Error Calibration and Normalization for Copy Numbers Analysis" (Mecan4CNA). It only requires CNV segmentation files as input, is platform independent, and has a high performance with limited hardware requirements. For a given multi-sample copy number dataset, Mecan4CNA can batch-normalize all samples to the corresponding true copy number levels of the main tumor clones. Experiments of Mecan4CNA on simulated data showed an overall accuracy of 93% and 91% in determining the normal level and single copy alteration (i.e. duplication or loss of one allele), respectively. Comparison of estimated normal levels and single copy alternations with existing methods and karyotyping data on the NCI-60 tumor cell line produced coherent results. To estimate the method's impact on downstream analyses, we performed GISTIC analyses on the original and Mecan4CNA normalized data from the Cancer Genome Atlas (TCGA) where the normalized data showed prominent improvements of both sensitivity and specificity in detecting focal regions. CONCLUSIONS: Mecan4CNA provides an advanced method for CNA data normalization, especially in meta-analyses involving large profile numbers and heterogeneous source data quality. With its informative output and visualization options, Mecan4CNA also can improve the interpretation of individual CNA profiles. Mecan4CNA is freely available as a Python package and through its code repository on Github.

Mountains and Chasms: Surveying the Oncogenomic Publication Landscape.

Cancers arise from the accumulation of somatic genome mutations, with varying contributions of intrinsic (i.e., genetic predisposition) and extrinsic (i.e., environmental) factors. For the understanding of malignant clones, precise information about their genomic composition has to be correlated with morphological, clinical, and individual features in the context of the available medical knowledge. Rapid improvements in molecular profiling techniques, the accumulation of a large amount of data in genomic alterations in human malignancies, and the expansion of bioinformatic tools and methodologies have facilitated the understanding of the molecular changes during oncogenesis, and their correlation with clinicopathological phenotypes. Far beyond a limited set of "driver" genes, oncogenomic profiling has identified a large variety of somatic mutations, and whole-genome sequencing studies of healthy individuals have improved the knowledge of heritable genome variation. Nevertheless, the main challenges arise from the skewed representation of individuals from varying population backgrounds in biomedical studies, and also through the limited extent in which some cancer entities are represented in the scientific literature. Content analyses of oncogenomic publications could provide guidance for the planning and support of future studies aiming at filling prominent knowledge gaps.

Oncology Informatics: Status Quo and Outlook.

Oncology has undergone rapid progress, with emerging developments in areas including cancer stem cells, molecularly targeted therapies, genomic analyses, and individually tailored immunotherapy. These advances have expanded the tools available in the fight against cancer. Some of these have seen broad media coverage resulting in justified public attention. However, these achievements have only been possible due to rapid developments in the expanding field of biomedical informatics and information technology (IT). Artificial intelligence, radiomics, electronic health records, and electronic patient-reported outcome measures (ePROMS) are only a few of the developments enabling further progress in oncology. The promising impact of IT in oncology will only become reality through a multidisciplinary approach to the complex challenges ahead.

PKCα and HMGB1 antagonistically control hydrogen peroxide-induced poly-ADP-ribose formation.

Harmful oxidation of proteins, lipids and nucleic acids is observed when reactive oxygen species (ROS) are produced excessively and/or the antioxidant capacity is reduced, causing 'oxidative stress'. Nuclear poly-ADP-ribose (PAR) formation is thought to be induced in response to oxidative DNA damage and to promote cell death under sustained oxidative stress conditions. However, what exactly triggers PAR induction in response to oxidative stress is incompletely understood. Using reverse phase protein array (RPPA) and in-depth analysis of key stress signaling components, we observed that PAR formation induced by H2O2 was mediated by the PLC/IP3R/Ca(2+)/PKCα signaling axis. Mechanistically, H2O2-induced PAR formation correlated with Ca(2+)-dependent DNA damage, which, however, was PKCα-independent. In contrast, PAR formation was completely lost upon knockdown of PKCα, suggesting that DNA damage alone was not sufficient for inducing PAR formation, but required a PKCα-dependent process. Intriguingly, the loss of PAR formation observed upon PKCα depletion was overcome when the chromatin structure-modifying protein HMGB1 was co-depleted with PKCα, suggesting that activation and nuclear translocation of PKCα releases the inhibitory effect of HMGB1 on PAR formation. Together, these results identify PKCα and HMGB1 as important co-regulators involved in H2O2-induced PAR formation, a finding that may have important relevance for oxidative stress-associated pathophysiological conditions.

arrayMap 2014: an updated cancer genome resource.

Somatic copy number aberrations (CNA) represent a mutation type encountered in the majority of cancer genomes. Here, we present the 2014 edition of arrayMap (http://www.arraymap.org), a publicly accessible collection of pre-processed oncogenomic array data sets and CNA profiles, representing a vast range of human malignancies. Since the initial release, we have enhanced this resource both in content and especially with regard to data mining support. The 2014 release of arrayMap contains more than 64,000 genomic array data sets, representing about 250 tumor diagnoses. Data sets included in arrayMap have been assembled from public repositories as well as additional resources, and integrated by applying custom processing pipelines. Online tools have been upgraded for a more flexible array data visualization, including options for processing user provided, non-public data sets. Data integration has been improved by mapping to multiple editions of the human reference genome, with the majority of the data now being available for the UCSC hg18 as well as GRCh37 versions. The large amount of tumor CNA data in arrayMap can be freely downloaded by users to promote data mining projects, and to explore special events such as chromothripsis-like genome patterns.

CNARA: reliability assessment for genomic copy number profiles.

BACKGROUND: DNA copy number profiles from microarray and sequencing experiments sometimes contain wave artefacts which may be introduced during sample preparation and cannot be removed completely by existing preprocessing methods. Besides, large derivative log ratio spread (DLRS) of the probes correlating with poor DNA quality is sometimes observed in genome screening experiments and may lead to unreliable copy number profiles. Depending on the extent of these artefacts and the resulting misidentification of copy number alterations/variations (CNA/CNV), it may be desirable to exclude such samples from analyses or to adapt the downstream data analysis strategy accordingly. RESULTS: Here, we propose a method to distinguish reliable genomic copy number profiles from those containing heavy wave artefacts and/or large DLRS. We define four features that adequately summarize the copy number profiles for reliability assessment, and train a classifier on a dataset of 1522 copy number profiles from various microarray platforms. The method can be applied to predict the reliability of copy number profiles irrespective of the underlying microarray platform and may be adapted for those sequencing platforms from which copy number estimates could be computed as a piecewise constant signal. Further details can be found at https://github.com/baudisgroup/CNARA . CONCLUSIONS: We have developed a method for the assessment of genomic copy number profiling data, and suggest to apply the method in addition to and after other state-of-the-art noise correction and quality control procedures. CNARA could be instrumental in improving the assessment of data used for genomic data mining experiments and support the reliable functional attribution of copy number aberrations especially in cancer research.

Progenetix: 12 years of oncogenomic data curation.

DNA copy number aberrations (CNAs) can be found in the majority of cancer genomes and are crucial for understanding the potential mechanisms underlying tumor initiation and progression. Since the first release in 2001, the Progenetix project (http://www.progenetix.org) has provided a reference resource dedicated to provide the most comprehensive collection of genome-wide CNA profiles. Reflecting the application of comparative genomic hybridization techniques to tens of thousands of cancer genomes, over the past 12 years our data curation efforts have resulted in a more than 60-fold increase in the number of cancer samples presented through Progenetix. In addition, new data exploration tools and visualization options have been added. In particular, the gene-specific CNA frequency analysis should facilitate the assignment of cancer genes to related cancer types. In addition, the new user file processing interface allows users to take advantage of the online tools, including various data representation options for proprietary data pre-publication. In this update article, we report recent improvements of the database in terms of content, user interface and online tools.

Chromothripsis-like patterns are recurring but heterogeneously distributed features in a survey of 22,347 cancer genome screens.

BACKGROUND: Chromothripsis is a recently discovered phenomenon of genomic rearrangement, possibly arising during a single genome-shattering event. This could provide an alternative paradigm in cancer development, replacing the gradual accumulation of genomic changes with a "one-off" catastrophic event. However, the term has been used with varying operational definitions, with the minimal consensus being a large number of locally clustered copy number aberrations. The mechanisms underlying these chromothripsis-like patterns (CTLP) and their specific impact on tumorigenesis are still poorly understood. RESULTS: Here, we identified CTLP in 918 cancer samples, from a dataset of more than 22,000 oncogenomic arrays covering 132 cancer types. Fragmentation hotspots were found to be located on chromosome 8, 11, 12 and 17. Among the various cancer types, soft-tissue tumors exhibited particularly high CTLP frequencies. Genomic context analysis revealed that CTLP rearrangements frequently occurred in genomes that additionally harbored multiple copy number aberrations (CNAs). An investigation into the affected chromosomal regions showed a large proportion of arm-level pulverization and telomere related events, which would be compatible to a number of underlying mechanisms. We also report evidence that these genomic events may be correlated with patient age, stage and survival rate. CONCLUSIONS: Through a large-scale analysis of oncogenomic array data sets, this study characterized features associated with genomic aberrations patterns, compatible to the spectrum of "chromothripsis"-definitions as previously used. While quantifying clustered genomic copy number aberrations in cancer samples, our data indicates an underlying biological heterogeneity behind these chromothripsis-like patterns, beyond a well defined "chromthripsis" phenomenon.

SIL1 mutations and clinical spectrum in patients with Marinesco-Sjogren syndrome.

Marinesco-Sjögren syndrome is a rare autosomal recessive multisystem disorder featuring cerebellar ataxia, early-onset cataracts, chronic myopathy, variable intellectual disability and delayed motor development. More recently, mutations in the SIL1 gene, which encodes an endoplasmic reticulum resident co-chaperone, were identified as the main cause of Marinesco-Sjögren syndrome. Here we describe the results of SIL1 mutation analysis in 62 patients presenting with early-onset ataxia, cataracts and myopathy or combinations of at least two of these. We obtained a mutation detection rate of 60% (15/25) among patients with the characteristic Marinesco-Sjögren syndrome triad (ataxia, cataracts, myopathy) whereas the detection rate in the group of patients with more variable phenotypic presentation was below 3% (1/37). We report 16 unrelated families with a total of 19 different SIL1 mutations. Among these mutations are 15 previously unreported changes, including single- and multi-exon deletions. Based on data from our screening cohort and data compiled from the literature we found that SIL1 mutations are invariably associated with the combination of a cerebellar syndrome and chronic myopathy. Cataracts were observed in all patients beyond the age of 7 years, but might be missing in infants. Six patients with SIL1 mutations had no intellectual disability, extending the known wide range of cognitive capabilities in Marinesco-Sjögren syndrome to include normal intelligence. Modestly constant features were somatic growth retardation, skeletal abnormalities and pyramidal tract signs. Examination of mutant SIL1 expression in cultured patient lymphoblasts suggested that SIL1 mutations result in severely reduced SIL1 protein levels irrespective of the type and position of mutations. Our data broaden the SIL1 mutation spectrum and confirm that SIL1 is the major Marinesco-Sjögren syndrome gene. SIL1 patients usually present with the characteristic triad but cataracts might be missing in young children. As cognitive impairment is not obligatory, patients without intellectual disability but a Marinesco-Sjögren syndrome-compatible phenotype should receive SIL1 mutation analysis. Despite allelic heterogeneity and many families with private mutations, the phenotype related to SIL1 mutations is relatively homogenous. Based on SIL1 expression studies we speculate that this may arise from a uniform effect of different mutations on protein expression.

High resolution copy number analysis of IRF4 translocation-positive diffuse large B-cell and follicular lymphomas.

Translocations affecting chromosome subband 6p25.3 containing the IRF4 gene have been recently described as characteristic alterations in a molecularly distinct subset of germinal center B-cell-derived lymphomas. Secondary changes have yet only been described in few of these lymphomas. Here, we performed array-comparative genomic hybridization and molecular inversion probe microarray analyses on DNA from 12 formalin-fixed paraffin-embedded and two fresh-frozen IRF4 translocation-positive lymphomas, which together with the previously published data on nine cases allowed the extension of copy number analyses to a total of 23 of these lymphomas. All except one case carried chromosomal imbalances, most frequently gains in Xq28, 11q22.3-qter, and 7q32.1-qter and losses in 6q13-16.1, 15q14-22.31, and 17p. No recurrent copy-neutral losses of heterozygosity were observed. TP53 point mutations were detected in three of six cases with loss of 17p. Overall this study unravels a recurrent pattern of secondary genetic alterations in IRF4 translocation-positive lymphomas.

cancercelllines.org-a novel resource for genomic variants in cancer cell lines.

Cancer cell lines are an important component in biological and medical research, enabling studies of cellular mechanisms as well as the development and testing of pharmaceuticals. Genomic alterations in cancer cell lines are widely studied as models for oncogenetic events and are represented in a wide range of primary resources. We have created a comprehensive, curated knowledge resource-cancercelllines.org-with the aim to enable easy access to genomic profiling data in cancer cell lines, curated from a variety of resources and integrating both copy number and single nucleotide variants data. We have gathered over 5600 copy number profiles as well as single nucleotide variant annotations for 16 000 cell lines and provide these data with mappings to the GRCh38 reference genome. Both genomic variations and associated curated metadata can be queried through the GA4GH Beacon v2 Application Programming Interface (API) and a graphical user interface with extensive data retrieval enabled using GA4GH data schemas under a permissive licensing scheme. Database URL: https://cancercelllines.org.

Data-driven information extraction and enrichment of molecular profiling data for cancer cell lines.

Motivation: With the proliferation of research means and computational methodologies, published biomedical literature is growing exponentially in numbers and volume. Cancer cell lines are frequently used models in biological and medical research that are currently applied for a wide range of purposes, from studies of cellular mechanisms to drug development, which has led to a wealth of related data and publications. Sifting through large quantities of text to gather relevant information on cell lines of interest is tedious and extremely slow when performed by humans. Hence, novel computational information extraction and correlation mechanisms are required to boost meaningful knowledge extraction. Results: In this work, we present the design, implementation, and application of a novel data extraction and exploration system. This system extracts deep semantic relations between textual entities from scientific literature to enrich existing structured clinical data concerning cancer cell lines. We introduce a new public data exploration portal, which enables automatic linking of genomic copy number variants plots with ranked, related entities such as affected genes. Each relation is accompanied by literature-derived evidences, allowing for deep, yet rapid, literature search, using existing structured data as a springboard. Availability and implementation: Our system is publicly available on the web at https://cancercelllines.org.

Short tandem repeat mutations regulate gene expression in colorectal cancer.

Short tandem repeat (STR) mutations are prevalent in colorectal cancer (CRC), especially in tumours with the microsatellite instability (MSI) phenotype. While STR length variations are known to regulate gene expression under physiological conditions, the functional impact of STR mutations in CRC remains unclear. Here, we integrate STR mutation data with clinical information and gene expression data to study the gene regulatory effects of STR mutations in CRC. We confirm that STR mutability in CRC highly depends on the MSI status, repeat unit size, and repeat length. Furthermore, we present a set of 1244 putative expression STRs (eSTRs) for which the STR length is associated with gene expression levels in CRC tumours. The length of 73 eSTRs is associated with expression levels of cancer-related genes, nine of which are CRC-specific genes. We show that linear models describing eSTR-gene expression relationships allow for predictions of gene expression changes in response to eSTR mutations. Moreover, we found an increased mutability of eSTRs in MSI tumours. Our evidence of gene regulatory roles for eSTRs in CRC highlights a mostly overlooked way through which tumours may modulate their phenotypes. Future extensions of these findings could uncover new STR-based targets in the treatment of cancer.

Recurrent loss of heterozygosity in 1p36 associated with TNFRSF14 mutations in IRF4 translocation negative pediatric follicular lymphomas.

Pediatric follicular lymphoma is a rare disease that differs genetically and clinically from its adult counterpart. With the exception of pediatric follicular lymphoma with IRF4-translocation, the genetic events associated with these lymphomas have not yet been defined. We applied array-comparative genomic hybridization and molecular inversion probe assay analyses to formalin-fixed paraffin-embedded tissues from 18 patients aged 18 years and under with IRF4 translocation negative follicular lymphoma. All evaluable cases lacked t(14;18). Only 6 of 16 evaluable cases displayed chromosomal imbalances with gains or amplifications of 6pter-p24.3 (including IRF4) and deletion and copy number neutral-loss of heterozygosity in 1p36 (including TNFRSF14) being most frequent. Sequencing of TNFRSF14 located in the minimal region of loss in 1p36.32 showed nine mutations in 7 cases from our series. Two subsets of pediatric follicular lymphoma were delineated according to the presence of molecular alterations, one with genomic aberrations associated with higher grade and/or diffuse large B-cell lymphoma component and more widespread disease, and another one lacking genetic alterations associated with more limited disease.

Improved multiplex ligation-dependent probe amplification analysis identifies a deleterious PMS2 allele generated by recombination with crossover between PMS2 and PMS2CL.

Heterozygous PMS2 germline mutations are associated with Lynch syndrome. Up to one third of these mutations are genomic deletions. Their detection is complicated by a pseudogene (PMS2CL), which--owing to extensive interparalog sequence exchange--closely resembles PMS2 downstream of exon 12. A recently redesigned multiplex ligation-dependent probe amplification (MLPA) assay identifies PMS2 copy number alterations with improved reliability when used with reference DNAs containing equal numbers of PMS2- and PMS2CL-specific sequences. We selected eight such reference samples--all publicly available--and used them with this assay to study 13 patients with PMS2-defective colorectal tumors. Three presented deleterious alterations: an Alu-mediated exon deletion; a 125-kb deletion encompassing PMS2 and four additional genes (two with tumor-suppressing functions); and a novel deleterious hybrid PMS2 allele produced by recombination with crossover between PMS2 and PMS2CL, with the breakpoint in intron 10 (the most 5' breakpoint of its kind reported thus far). We discuss mechanisms that might generate this allele in different chromosomal configurations (and their diagnostic implications) and describe an allele-specific PCR assay that facilitates its detection. Our data indicate that the redesigned PMS2 MLPA assay is a valid first-line option. In our series, it identified roughly a quarter of all PMS2 mutations.