CopyMap: localization and calling of copy number variation by joint analysis of hybridization data from multiple individuals.
Ontology highlight
ABSTRACT: Summary
The program package CopyMap identifies copy number variation from oligo-hybridization and CGH data. Using a time-dependent hidden Markov model to combine evidence of copy number variants (CNVs) across multiple carriers, CopyMap is substantially more accurate than standard hidden Markov methods in identifying CNVs and calling CNV-carriers. Moreover, CopyMap provides more precise estimates of CNV-boundaries.Availability
The C-source code and detailed documentation for the program CopyMap is available on the Internet at http://www.sph.umich.edu/csg/szoellner/
SUBMITTER: Zollner S
PROVIDER: S-EPMC2981497 | biostudies-literature | 2010 Nov
REPOSITORIES: biostudies-literature
Publications
CopyMap: localization and calling of copy number variation by joint analysis of hybridization data from multiple individuals.
Bioinformatics (Oxford, England) 20101005 21
<h4>Summary</h4>The program package CopyMap identifies copy number variation from oligo-hybridization and CGH data. Using a time-dependent hidden Markov model to combine evidence of copy number variants (CNVs) across multiple carriers, CopyMap is substantially more accurate than standard hidden Markov methods in identifying CNVs and calling CNV-carriers. Moreover, CopyMap provides more precise estimates of CNV-boundaries.<h4>Availability</h4>The C-source code and detailed documentation for the p ...[more]
Similar Datasets
Project description:Copy number variants (CNVs) play an important role in a number of human diseases, but the accurate calling of CNVs remains challenging. Most current approaches to CNV detection use raw read alignments, which are computationally intensive to process. We use a regression tree-based approach to call germline CNVs from whole-genome sequencing (WGS, >18x) variant call sets in 6,898 samples across four European cohorts, and describe a rich large variation landscape comprising 1,320 CNVs. Eighty-one percent of detected events have been previously reported in the Database of Genomic Variants. Twenty-three percent of high-quality deletions affect entire genes, and we recapitulate known events such as the GSTM1 and RHD gene deletions. We test for association between the detected deletions and 275 protein levels in 1,457 individuals to assess the potential clinical impact of the detected CNVs. We describe complex CNV patterns underlying an association with levels of the CCL3 protein (MAF = 0.15, p = 3.6x10-12 ) at the CCL3L3 locus, and a novel cis-association between a low-frequency NOMO1 deletion and NOMO1 protein levels (MAF = 0.02, p = 2.2x10-7 ). This study demonstrates that existing population-wide WGS call sets can be mined for germline CNVs with minimal computational overhead, delivering insight into a less well-studied, yet potentially impactful class of genetic variant.
| S-EPMC8653900 | biostudies-literature
Project description:MotivationThe complexity of a large number of recently discovered copy number polymorphisms is much higher than initially thought, thus making it more difficult to detect them in the presence of significant measurement noise. In this scenario, separate normalization and segmentation is prone to lead to many false detections of changes in copy number. New approaches capable of jointly modeling the copy number and the non-copy number (noise) hybridization effects across multiple samples will potentially lead to more accurate results.MethodsIn this article, the genome alteration detection analysis (GADA) approach introduced in our previous work is extended to a multiple sample model. The copy number component is independent for each sample and uses a sparse Bayesian prior, while the reference hybridization level is not necessarily sparse but identical on all samples. The expectation maximization (EM) algorithm used to fit the model iteratively determines whether the observed hybridization levels are more likely due to a copy number variation or to a shared hybridization bias.ResultsThe new proposed approach is compared with the currently used strategy of separate normalization followed by independent segmentation of each array. Real microarray data obtained from HapMap samples are randomly partitioned to create different reference sets. Using the new approach, copy number and reference intensity estimates are significantly less variable if the reference set changes; and a higher consistency on copy numbers detected within HapMap family trios is obtained. Finally, the running time to fit the model grows linearly in the number samples and probes.Availabilityhttp://biron.usc.edu/~piquereg/GADA.
| S-EPMC2732310 | biostudies-literature
Project description:BackgroundPediatric onset multiple sclerosis (MS) accounts for 2-4% of all MS. It is unknown whether the disease shares the same underlying pathophysiology found in adult patients or an extreme early onset phenotype triggered by distinct biological mechanisms. It has been hypothesized that copy number variations (CNVs) may result in extreme early onset diseases because CNVs can have major effects on many genes in large genomic regions.Objectives and methodsThe objective of the current research was to identify CNVs, with a specific focus on de novo CNVs, potentially causing early onset MS by competitively hybridizing 30 white non-Hispanic pediatric MS patients with each of their parents via comparative genomic hybridization (CGH) analysis on the Agilent 1M CGH array.Results and discussionWe identified 10 CNVs not overlapping with any CNV regions currently reported in the Database of Genomic Variants (DGV). Fifty-five putatively de novo CNVs were also identified: all but one common in the DGV. We found the single rare CNV was a private variation harboring the SACS gene. SACS mutations cause autosomal-recessive spastic ataxia of Charlevoix-Saguenay (ARSACS) disease. Additional clinical review revealed that the patient with the SACS gene CNV shared some features of both MS and ARSACS.ConclusionsThis is the first reported study analyzing pediatric MS CNVs. While not yielding causal variation in our initial pediatric dataset, our approach confirmed diagnosis of an ARSACS-like disease in addition to MS in the affected individual, which led to a more complete understanding of the patient's disease course and prognosis.
| S-EPMC5507577 | biostudies-literature
Project description:The usage of next-generation sequencing (NGS) to detect copy number variation (CNV) is widely accepted in cancer research. Based on an AluScanCNV software developed by us previously, an AluScanCNV2 software has been developed in the present study as an R package that performs CNV detection from NGS data obtained through AluScan, whole-genome sequencing or other targeted NGS platforms. Its applications would include the expedited usage of somatic CNVs for cancer subtyping, and usage of recurrent germline CNVs to perform machine learning-assisted prediction of a test subject's susceptibility to cancer.
| S-EPMC6411622 | biostudies-literature
Project description:MotivationDNA copy number variants (CNVs) are gains and losses of segments of chromosomes, and comprise an important class of genetic variation. Recently, various microarray hybridization-based techniques have been developed for high-throughput measurement of DNA copy number. In many studies, multiple technical platforms or different versions of the same platform were used to interrogate the same samples; and it became necessary to pool information across these multiple sources to derive a consensus molecular profile for each sample. An integrated analysis is expected to maximize resolution and accuracy, yet currently there is no well-formulated statistical method to address the between-platform differences in probe coverage, assay methods, sensitivity and analytical complexity.ResultsThe conventional approach is to apply one of the CNV detection ('segmentation') algorithms to search for DNA segments of altered signal intensity. The results from multiple platforms are combined after segmentation. Here we propose a new method, Multi-Platform Circular Binary Segmentation (MPCBS), which pools statistical evidence across platforms during segmentation, and does not require pre-standardization of different data sources. It involves a weighted sum of t-statistics, which arises naturally from the generalized log-likelihood ratio of a multi-platform model. We show by comparing the integrated analysis of Affymetrix and Illumina SNP array data with Agilent and fosmid clone end-sequencing results on eight HapMap samples that MPCBS achieves improved spatial resolution, detection power and provides a natural consensus across platforms. We also apply the new method to analyze multi-platform data for tumor samples.AvailabilityThe R package for MPCBS is registered on R-Forge (http://r-forge.r-project.org/) under project name MPCBS.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC2852203 | biostudies-literature
Project description:BackgroundVariations in DNA copy number carry information on the modalities of genome evolution and mis-regulation of DNA replication in cancer cells. Their study can help localize tumor suppressor genes, distinguish different populations of cancerous cells, and identify genomic variations responsible for disease phenotypes. A number of different high throughput technologies can be used to identify copy number variable sites, and the literature documents multiple effective algorithms. We focus here on the specific problem of detecting regions where variation in copy number is relatively common in the sample at hand. This problem encompasses the cases of copy number polymorphisms, related samples, technical replicates, and cancerous sub-populations from the same individual.ResultsWe present a segmentation method named generalized fused lasso (GFL) to reconstruct copy number variant regions. GFL is based on penalized estimation and is capable of processing multiple signals jointly. Our approach is computationally very attractive and leads to sensitivity and specificity levels comparable to those of state-of-the-art specialized methodologies. We illustrate its applicability with simulated and real data sets.ConclusionsThe flexibility of our framework makes it applicable to data obtained with a wide range of technology. Its versatility and speed make GFL particularly useful in the initial screening stages of large data sets.
| S-EPMC3534631 | biostudies-literature
Project description:PHACE syndrome is the association of large segmental facial hemangiomas and congenital anomalies, such as posterior fossa malformations, cerebral arterial anomalies, coarctation of the aorta, eye anomalies, and sternal defects. To date, the reported cases of PHACE syndrome have been sporadic, suggesting that PHACE may have a complex pathogenesis. We report here genomic copy number variation (CNV) analysis of 98 individuals with PHACE syndrome as a first step in deciphering a potential genetic basis of PHACE syndrome. A total of 3,772 CNVs (2,507 duplications and 1,265 deletions) were detected in 98 individuals with PHACE syndrome. CNVs were then eliminated if they failed to meet established criteria for quality, spanned centromeres, or did not contain genes. CNVs were defined as "rare" if not documented in the database of genomic variants. Ten rare CNVs were discovered (size range: 134-406 ?kb), located at 1q32.1, 1q43, 3q26.32-3q26.33, 3p11.1, 7q33, 10q24.32, 12q24.13, 17q11.2, 18p11.31, and Xq28. There were no rare CNV events that occurred in more than one subject. Therefore, further study is needed to determine the significance of these CNVs in the pathogenesis of PHACE syndrome.
| S-EPMC3971866 | biostudies-literature
Project description:BACKGROUND:Copy number variants (CNVs) are known to play an important role in the development and progression of several diseases. However, detection of CNVs with whole-exome sequencing (WES) experiments is challenging. Usually, additional experiments have to be performed. FINDINGS:We developed a novel algorithm for somatic CNV calling in matched WES data called "CopyDetective". Different from other approaches, CNV calling with CopyDetective consists of a 2-step procedure: first, quality analysis is performed, determining individual detection thresholds for every sample. Second, actual CNV calling on the basis of the previously determined thresholds is performed. Our algorithm evaluates the change in variant allele frequency of polymorphisms and reports the fraction of affected cells for every CNV. Analyzing 4 WES data sets (n = 100) we observed superior performance of CopyDetective compared with ExomeCNV, VarScan2, ControlFREEC, ExomeDepth, and CNV-seq. CONCLUSIONS:Individual detection thresholds reveal that not every WES data set is equally apt for CNV calling. Initial quality analyses, determining individual detection thresholds-as realized by CopyDetective-can and should be performed prior to actual variant calling.
| S-EPMC7604644 | biostudies-literature
Project description:Hidden Markov model, based on Li and Stephens model that takes into account chromosome sharing of multiple individuals, results in mainstream haplotype phasing algorithms for genotyping arrays and next-generation sequencing (NGS) data. However, existing methods based on this model assume that the allele count data are independently observed at individual sites and do not consider haplotype informative reads, i.e. reads that cover multiple heterozygous sites, which carry useful haplotype information. In our previous work, we developed a new hidden Markov model to incorporate a two-site joint emission term that captures the haplotype information across two adjacent sites. Although our model improves the accuracy of genotype calling and haplotype phasing, haplotype information in reads covering non-adjacent sites and/or more than two adjacent sites is not used because of the severe computational burden.We develop a new probabilistic model for genotype calling and haplotype phasing from NGS data that incorporates haplotype information of multiple adjacent and/or non-adjacent sites covered by a read over an arbitrary distance. We develop a new hybrid Markov Chain Monte Carlo algorithm that combines the Gibbs sampling algorithm of HapSeq and Metropolis-Hastings algorithm and is computationally feasible. We show by simulation and real data from the 1000 Genomes Project that our model offers superior performance for haplotype phasing and genotype calling for population NGS data over existing methods.HapSeq2 is available at www.ssg.uab.edu/hapseq/.
| S-EPMC3777110 | biostudies-literature
Project description:Genomic variation includes single-nucleotide variants, small insertions or deletions (indels), and copy number variants (CNVs). CNVs affect gene expression by altering the genome structure and transposable elements within a region. CNVs are greater than 1 kb in size; hence, CNVs can produce more variation than can individual single-nucleotide variations that are detected by next-generation sequencing. Multiple system atrophy (MSA) is an ?-synucleinopathy adult-onset disorder. Pathologically, it is characterized by insoluble aggregation of filamentous ?-synuclein in brain oligodendrocytes. Generally, MSA is sporadic, although there are rare cases of familial MSA. In addition, the frequencies of the clinical phenotypes differ considerably among countries. Reports indicate that genetic factors play roles in the mechanisms involved in the pathology and onset of MSA. To evaluate the genetic background of this disorder, we attempted to determine whether there are differences in CNVs between patients with MSA and normal control subjects. We found that the number of CNVs on chromosomes 5, 22, and 4 was increased in MSA; 3 CNVs in non-coding regions were considered risk factors for MSA. Our results show that CNVs in non-coding regions influence the expression of genes through transcription-related mechanisms and potentially increase subsequent structural alterations of chromosomes. Therefore, these CNVs likely play roles in the molecular mechanisms underlying MSA.
| S-EPMC5708077 | biostudies-literature