Search

Real-time quantitative PCR is currently the most widely used method for the human pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) identification. Due to the rapid evolution of the SARS-CoV-2 genome, novel mutations on the primer binding sites will cause the failure of PCR. Therefore, in addition to a well-designed primer set, these primers need to be updated and evaluated regularly to ensure that the rapidly evolving genome primers can be amplified. In this protocol, (1) we firstly use assembled genome sequences in the SARS-CoV-2 database to identify and characterize indels and point mutations; (2) design primers skipping the sites of mutations; (3) check the coverage of the primers with the daily update SARS-CoV-2 database; (4) redesign them if novel mutations found in the primer binding sites. Although this protocol takes SARS-CoV-2 as an example, it is suitable for other species that have genomes accumulating mutations over time.

Quality control (QC) for lab-designed primers is crucial for the success of a polymerase chain reaction (PCR). Here, we present MFEprimer-3.0, a functional primer quality control program for checking non-specific amplicons, dimers, hairpins and other parameters. The new features of the current version include: (i) more sensitive binding site search using the updated k-mer algorithm that allows mismatches within the k-mer, except for the first base at the 3′ end. The binding sites of each primer with a stable 3′ end are listed in the output; (ii) new algorithms for rapidly identifying self-dimers, cross-dimers and hairpins; (iii) the command-line version, which has an added option of JSON output to enhance the versatility of MFEprimer by acting as a QC step in the ‘primer design → quality control → redesign’ pipeline; (iv) a function for checking whether the binding sites contain single nucleotide polymorphisms (SNPs), which will affect the consistency of binding efficiency among different samples. In summary, MFEprimer-3.0 is updated with the well-tested PCR primer QC program and it can be integrated into various PCR primer design applications as a QC module. The MFEprimer-3.0 server is freely accessible without any login requirement at: https://mfeprimer3.igenetech.com/ and https://www.mfeprimer.com/. The source code for the command-line version is available upon request.

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common malignancy worldwide. Somatic copy number alterations (CNAs) play a significant role in the development of this lethal cancer. In this study, we present a meta-analysis of CNAs for a total of 1,395 HNSCC samples. Publicly available R packages and in-house scripts were used for genomic array data processing, including normalization, segmentation and CNA calling. We detected 125 regions of significant gains or losses using GISTIC algorithm and found several potential driver genes in these regions. The incidence of chromothripsis in HNSCC was estimated to be 6%, and the chromosome pulverization hotspot regions were detected. We determined 323 genomic locations significantly enriched for breakpoints, which indicate HNSCC-specific genomic instability regions. Unsupervised clustering of genome-wide CNA data revealed a sub-cluster predominantly composed of nasopharynx tumors and presented a large proportion of HPV-positive samples. These results will facilitate the discovery of therapeutic candidates and extend our molecular understanding of HNSCC.

Background The yak is a species of livestock which is crucial for local communities of the Qinghai-Tibet Plateau and adjacent regions and naturally owns one more thoracic vertebra than cattle. Recently, a sub-population of yak termed as the Jinchuan yak has been identified with over half its members own a thoracolumbar vertebral formula of T15L5 instead of the natural T14L5 arrangement. The novel T15L5 positioning is a preferred genetic trait leading to enhanced meat and milk production. Selective breeding of this trait would have great agricultural value and exploration of the molecular mechanisms underlying this trait would both accelerate this process and provide us insight into the development and regulation of somitogenesis. Results Here we investigated the genetic background of the Jinchuan yak through resequencing fifteen individuals, comprising five T15L5 individuals and ten T14L5 individuals with an average sequencing depth of > 10X, whose thoracolumbar vertebral formulae were confirmed by anatomical observation. Principal component analysis, linkage disequilibrium analysis, phylogenetic analysis, and selective sweep analysis were carried out to explore Jinchuan yak’s genetic background. Three hundred and thirty candidate markers were identified as associated with the additional thoracic vertebrae and target sequencing was used to validate seven carefully selected markers in an additional 51 Jinchuan yaks. The accuracies of predicting 15 thoracic vertebrae and 20 thoracolumbar vertebrae with these 7 markers were 100.00 and 33.33% despite they both could only represent 20% of all possible genetic diversity. Two genes, PPP2R2B and TBLR1, were found to harbour the most candidate markers associated with the trait and likely contribute to the unique somitic number and identity according to their reported roles in the mechanism of somitogenesis. Conclusions Our findings provide a clear depiction of the Jinchuan yak’s genetic background and a solid foundation for marker-assistant selection. Further exploitation of this unique population and trait could be promoted with the aid of our genomic resource.

Triple‑negative breast cancer (TNBC) accounts for ~15% of all breast cancer diagnoses each year. Patients with TNBC tend to have a higher risk for early relapse and a worse prognosis. TNBC is characterized by extensive somatic copy number alterations (CNAs). However, the DNA CNA profile of TNBC remains to be extensively investigated. The present study assessed the genomic profile of CNAs in 201 TNBC samples, aiming to identify recurrent CNAs that may drive the pathogenesis of TNBC. In total, 123 regions of significant amplification and deletion were detected using the Genomic Identification of Significant Targets in Cancer algorithm, and potential driver genes for TNBC were identified. A total of 31 samples exhibited signs of chromothripsis and revealed chromosome pulverization hotspot regions. The present study further determined 199 genomic locations that were significantly enriched for breakpoints, which indicated TNBC‑specific genomic instability regions. Unsupervised hierarchical clustering of tumors resulted in three main subgroups that exhibited distinct CNA profiles, which may reveal the heterogeneity of molecular mechanisms in TNBC subgroups. These results will extend the molecular understanding of TNBC and will facilitate the discovery of therapeutic and diagnostic target candidates.

Comprehensive genomic analyses of cancers have revealed substantial intrapatient molecular heterogeneities that may explain some instances of drug resistance and treatment failures. Examination of the clonal composition of an individual tumor and its evolution through disease progression and treatment may enable identification of precise therapeutic targets for drug design. Multi-region and single-cell sequencing are powerful tools that can be used to capture intratumor heterogeneity. Here, we present a database we’ve named CancerTracer (http://cailab.labshare.cn/cancertracer): a manually curated database designed to track and characterize the evolutionary trajectories of tumor growth in individual patients. We collected over 6000 tumor samples from 1548 patients corresponding to 45 different types of cancer. Patient-specific tumor phylogenetic trees were constructed based on somatic mutations or copy number alterations identified in multiple biopsies. Using the structured heterogeneity data, researchers can identify common driver events shared by all tumor regions, and the heterogeneous somatic events present in different regions of a tumor of interest. The database can also be used to investigate the phylogenetic relationships between primary and metastatic tumors. It is our hope that CancerTracer will significantly improve our understanding of the evolutionary histories of tumors, and may facilitate the identification of predictive biomarkers for personalized cancer therapies.

Chromothripsis is a recently observed phenomenon in cancer cells in which one or several chromosomes shatter into pieces with subsequent inaccurate reassembly and clonal propagation. This type of event generates a potentially vast number of mutations within a relatively short-time period, and has been considered as a new paradigm in cancer development. Despite recent advances, much work is still required to better understand the molecular mechanisms of this phenomenon, and thus an easy-to-use tool is in urgent need for automatically detecting and annotating chromothripsis. Here we present CTLPScanner, a web server for detection of chromothripsis-like pattern (CTLP) in genomic array data. The output interface presents intuitive graphical representations of detected chromosome pulverization region, as well as detailed results in table format. CTLPScanner also provides additional information for associated genes in chromothripsis region to help identify the potential candidates involved in tumorigenesis. To assist in performing meta-data analysis, we integrated over 50 000 pre-processed genomic arrays from The Cancer Genome Atlas and Gene Expression Omnibus into CTLPScanner. The server allows users to explore the presence of chromothripsis signatures from public data resources, without carrying out any local data processing. CTLPScanner is freely available at http://cgma.scu.edu.cn/CTLPScanner/.

Background The integration of DNA methylation and copy number alteration data promises to provide valuable insight into the underlying molecular mechanisms responsible for cancer initiation and progression. However, the generation and processing of these datasets are costly and time-consuming if carried out separately. The Illumina Infinium HumanMethylation450 BeadChip, initially designed for the evaluation of DNA methylation levels, allows copy number variant calling using bioinformatics tools. Results A substantial amount of Infinium HumanMethylation450 data across various cancer types has been accumulated in recent years and is a valuable resource for large-scale data analysis. Here we present MethCNA, a comprehensive database for genomic and epigenomic data integration in human cancer. In the current release, MethCNA contains about 10,000 tumor samples representing 37 cancer types. All raw array data were collected from The Cancer Genome Atlas and NCBI Gene Expression Omnibus database and analyzed using a pipeline that integrated multiple computational resources and tools. The normalized copy number aberration data and DNA methylation alterations were obtained. We provide a user-friendly web-interface for data mining and visualization. Conclusions The Illumina Infinium HumanMethylation450 BeadChip enables the interrogation and integration of both genomic and epigenomic data from exactly the same DNA specimen, and thus can aid in distinguishing driver from passenger mutations in cancer. We expect MethCNA will enable researchers to explore DNA methylation and copy number alteration patterns, identify key oncogenic drivers in cancer, and assist in the development of targeted therapies. MethCNA is publicly available online at http://cgma.scu.edu.cn/MethCNA. Electronic supplementary material The online version of this article (10.1186/s12864-018-4525-0) contains supplementary material, which is available to authorized users.

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 (CGH) 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. Additionally, 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., Accepted at Nucleic Acid Research (NAR 2014 database issue)

At the onset of lactation, there are three distinct stages of mammary tissue development and function including mammogenesis, colostrogenesis, and lactogenesis. The mechanism of the transition from colostrogenesis to lactogenesis of Maiwa Yak is still unknown. In this study, mammary tissues from three Maiwa Yaks were collected at 1 and 30 d after parturition for transcriptome exploration using Affymetrix Bovine Genome Arrays. Comparing 1 and 30 d results, a total of 517 annotated differentially expressed genes (DEGs) were identified at the criteria of a P ≤ 0.05. The ratio of up-regulated genes to the down-regulated ones was around 1:2 (more specifically, 164:353). To depict the profile of DEGs, a dynamic impact approach (DIA) was used to analyse the microarray data based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. GO terms “fatty acid transport” and “monocarboxylic acid transport” were significantly induced during the colostrum period. The strongly impacted KEGG pa...

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. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3074-7) contains supplementary material, which is available to authorized users.

Objectives MicroRNAs (miRNAs) are endogenous ~23 nucleotides (nt) RNAs, regulating gene expression by pairing to the mRNAs of protein-coding genes to direct their post-transcriptional repression. Both in normal and aberrant activities, miRNAs contribute to a recurring paradigm of cellular behaviors in pathological settings, especially in gerontology. Autophagy, a multi-step lysosomal degradation process with function to degrade long-lived proteins and damaged organelles, has significant impact on gerontology. Thus, elucidating how miRNAs participate in autophagy may enlarge the scope of miRNA in autophagy and facilitate researches in gerontology. Materials and methods Herein, based upon the published studies, predicted targets and gerontology-related diseases, we constructed a web resource named Gerontology-Autophagic-MicroRNA Database (GAMDB) (http://gamdb.liu-lab.com/index.php), which contained 836 autophagy-related miRNAs, 197 targeted genes/proteins and 56 aging-related diseases such as Parkinson' disease, Alzheimer's disease and Huntington's disease. Results and Conclusion We made use of large amounts of data to elucidate the intricate relationships between microRNA-regulated autophagic mechanisms and gerontology. This database will facilitate better understanding of autophagy regulation network in gerontology and thus promoting gerontology-related therapy in the future.

Protein–protein interactions (PPIs) play crucial roles in a number of biological processes. Recently, protein interaction networks (PINs) for several model organisms and humans have been generated, but few large-scale researches for mice have ever been made neither experimentally nor computationally. In the work, we undertook an effort to map a mouse PIN, in which protein interactions are hidden in enormous amount of biomedical literatures. Following a co-occurrence-based text-mining approach, a probabilistic model—naive Bayesian was used to filter false-positive interactions by integrating heterogeneous kinds of evidence from genomic and proteomic datasets. A support vector machine algorithm was further used to choose protein pairs with physical interactions. By comparing with the currently available PPI datasets from several model organisms and humans, it showed that the derived mouse PINs have similar topological properties at the global level, but a high local divergence. The mouse protein interaction dataset is stored in the Mouse protein–protein interaction DataBase (MppDB) that is useful source of information for system-level understanding of gene function and biological processes in mammals. Access to the MppDB database is public available at http://bio.scu.edu.cn/mppi .

Predicting protein-coding genes still remains a signif icant challenge. Although a variety of computational programs that use commonly machine learning methods have emerged, the accuracy of predictions remains a low level when implementing in large genomic sequences. Moreover, computational gene f inding in newly sequenced genomes is especially a dif f icult task due to the absence of a training set of abundant validated genes. Here we present a new gene-f inding program, SCGPred, to improve the accuracy of prediction by combining multiple sources of evidence. SCGPred can perform both supervised method in previously well-studied genomes and unsupervised one in novel genomes. By testing with datasets composed of large DNA sequences from human and a novel genome of Ustilago maydi ,S CGPred gains a signif icant improvement in comparison to the popular ab initio gene predictors. We also demonstrate that SCGPred can signif icantly improve prediction in novel genomes by combining several foreign gene f inders with similarity alignments, which is superior to other unsupervised methods. Therefore, SCGPred can serve as an alternative gene-f inding tool for newly sequenced eukaryotic genomes. The program is freely available at http://bio.scu.edu.cn/SCGPred/.

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.

Objectives Autophagy, a highly conserved lysosomal degradation process in eukaryotic cells, can digest long-lived proteins and damaged organelles through vesicular trafficking pathways. Nowadays, mechanisms of autophagy have been gradually elucidated and thus the discovery of small-molecule drugs targeting autophagy has always been drawing much attention. So far, some autophagy-related web servers have been available online to facilitate scientists to obtain the information relevant to autophagy conveniently, such as HADb, CTLPScanner, iLIR server and ncRDeathDB. However, to the best of our knowledge, there is not any web server available about the autophagy-modulating compounds. Methods According to published articles, all the compounds and their relations with autophagy were anatomized. Subsequently, an online Autophagic Compound Database (ACDB) (http://www.acdbliulab.com/) was constructed, which contained information of 357 compounds with 164 corresponding signalling pathways and potential targets in different diseases. Results We achieved a great deal of information of autophagy-modulating compounds, including compounds, targets/pathways and diseases. ACDB is a valuable resource for users to access to more than 300 curated small-molecule compounds correlated with autophagy. Conclusions Autophagic compound database will facilitate to the discovery of more novel therapeutic drugs in the near future.

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.

Zusammenfassung Array-basierte Methoden zum Nachweis chromsomaler Imbalancen haben in der vergangenen Dekade zunehmende Bedeutung in der tumorgenetischen Analytik gewonnen. Dabei werden im Wesentlichen auf Array-CGH („comparative genomic hybridization“) und SNP(Single-nucleotide-polymorphism)-Array basierte Technologien unterschieden, die je nach Fragestellung, Ausgangsmaterial und gewünschter Auflösung Vor- und Nachteile haben. So erlauben SNP-basierte Methoden im Gegensatz zum klassischen Array-CGH-Ansatz den gleichzeitigen Nachweis von chromosomalen Imbalancen und von Verlust der Heterozygotie ohne Veränderung der Kopienzahl („copy-neutral loss of heterozygosity“, CN-LOH). Bei allen Array-basierten Analysen von Tumoren ist zu beachten, dass im Gegensatz zu den Analysen zum Nachweis konstitutioneller Veränderungen zumeist nicht alle untersuchten Zellen dem neoplastischen Klon oder einem chromosomal aberranten Subklon angehören. Einsatzgebiete von Array-basierten Technologien bei Tumoren sind z. B. die Charakterisierung pathogenetisch relevanter Imbalancen, die Definition von molekularen und klinischen Subgruppen von Tumoren oder die Identifizierung von Targets für eine individualisierte Therapie.

Little is known about frequency, association with clinical characteristics, and prognostic impact of DNA copy number alterations (CNA) on survival in central primitive neuroectodermal tumors (CNS-PNET) and tumors of the pineal region. Searches of MEDLINE, Pubmed, and EMBASE--after the original description of comparative genomic hybridization in 1992 and July 2010--identified 15 case series of patients with CNS-PNET and tumors of the pineal region whose tumors were investigated for genome-wide CNA. One additional case study was identified from contact with experts. Individual patient data were extracted from publications or obtained from investigators, and CNAs were converted to a digitized format suitable for data mining and subgroup identification. Summary profiles for genomic imbalances were generated from case-specific data. Overall survival (OS) was estimated using the Kaplan-Meier method, and by univariable and multivariable Cox regression models. In their overall CNA profiles, low grade tumors of the pineal region clearly diverged from CNS-PNET and pineoblastoma. At a median follow-up of 89 months, 7-year OS rates of CNS-PNET, pineoblastoma, and low grade tumors of the pineal region were 22.9 ± 6, 0 ± 0, and 87.5 ± 12 %, respectively. Multivariable analysis revealed that histology (CNS-PNET), age (≤2.5 years), and possibly recurrent CNAs were associated with unfavorable OS. DNA copy number profiling suggests a close relationship between CNS-PNET and pineoblastoma. Low grade tumors of the pineal region differed from CNS-PNET and pineoblastoma. Due to their high biological and clinical variability, a coordinated prospective validation in future studies is necessary to establish robust risk factors.

Background: The delineation of genomic copy number abnormalities (CNAs) from cancer samples has been instrumental for identification of tumor suppressor genes and oncogenes and proven useful for clinical marker detection. An increasing number of projects have mapped CNAs using high-resolution microarray based techniques. So far, no single resource does provide a global collection of readily accessible oncoge- nomic array data. Methodology/Principal Findings: We here present arrayMap, a curated reference database and bioinformatics resource targeting copy number profiling data in human cancer. The arrayMap database provides a platform for meta-analysis and systems level data integration of high-resolution oncogenomic CNA data. To date, the resource incorporates more than 40,000 arrays in 224 cancer types extracted from several resources, including the NCBI's Gene Expression Omnibus (GEO), EBIs ArrayExpress (AE), The Cancer Genome Atlas (TCGA), publication supplements and direct submissions. For the majority of the included datasets, probe level and integrated visualization facilitate gene level and genome wide data re- view. Results from multi-case selections can be connected to downstream data analysis and visualization tools. Conclusions/Significance: To our knowledge, currently no data source provides an extensive collection of high resolution oncogenomic CNA data which readily could be used for genomic feature mining, across a representative range of cancer entities. arrayMap represents our effort for providing a long term platform for oncogenomic CNA data independent of specific platform considerations or specific project dependence. The online database can be accessed at http://www.arraymap.org., Comment: 17 pages, 5 inline figures, 3 tables, supplementary figures/tables split into 4 PDF files; manuscript submitted to PLoS ONE

BACKGROUND: Regional genomic copy number alterations (CNA) are observed in the vast majority of cancers. Besides specifically targeting well-known, canonical oncogenes, CNAs may also play more subtle roles in terms of modulating genetic potential and broad gene expression patterns of developing tumors. Any significant differences in the overall CNA patterns between different cancer types may thus point towards specific biological mechanisms acting in those cancers. In addition, differences among CNA profiles may prove valuable for cancer classifications beyond existing annotation systems. PRINCIPAL FINDINGS: We have analyzed molecular-cytogenetic data from 25579 tumors samples, which were classified into 160 cancer types according to the International Classification of Disease (ICD) coding system. When correcting for differences in the overall CNA frequencies between cancer types, related cancers were often found to cluster together according to similarities in their CNA profiles. Based on a randomization approach, distance measures from the cluster dendrograms were used to identify those specific genomic regions that contributed significantly to this signal. This approach identified 43 non-neutral genomic regions whose propensity for the occurrence of copy number alterations varied with the type of cancer at hand. Only a subset of these identified loci overlapped with previously implied, highly recurrent (hot-spot) cytogenetic imbalance regions. CONCLUSIONS: Thus, for many genomic regions, a simple null-hypothesis of independence between cancer type and relative copy number alteration frequency can be rejected. Since a subset of these regions display relatively low overall CNA frequencies, they may point towards second-tier genomic targets that are adaptively relevant but not necessarily essential for cancer development.

Background: Copy number alterations (CNA) play a key role in cancer development and progression. Since more than one CNA can be detected in most tumors, frequently co-occurring genetic CNA may point to cooperating cancer related genes. Existing methods for co-occurrence evaluation so far have not considered the overall heterogeneity of CNA per tumor, resulting in a preferential detection of frequent changes with limited specificity for each association due to the high genetic instability of many samples. Method: We hypothesize that in cancer some linkage-independent CNA may display a non-random co-occurrence, and that these CNA could be of pathogenetic relevance for the respective cancer. We also hypothesize that the statistical relevance of co-occurring CNA may depend on the sample specific CNA complexity. We verify our hypotheses with a simulation based algorithm CDCOCA (complexity dependence of co-occurring chromosomal aberrations). Results: Application of CDCOCA to example data sets identified co-occurring CNA from low complex background which otherwise went unnoticed. Identification of cancer associated genes in these co-occurring changes can provide insights of cooperative genes involved in oncogenesis. Conclusions: We have developed a method to detect associations of regional copy number abnormalities in cancer data. Along with finding statistically relevant CNA co-occurrences, our algorithm points towards a generally low specificity for co-occurrence of regional imbalances in CNA rich samples, which may have negative impact on pathway modeling approaches relying on frequent CNA events., BMC Medical Genomics, 4, ISSN:1755-8794