Copy Number Variations Detection Methods Research Articles

A systematic benchmark of copy number variation detection tools for high density SNP genotyping arrays.

Copy Number Variations (CNVs) are crucial in various diseases, especially cancer, but detecting them accurately from SNP genotyping arrays remains challenging. Therefore, this study benchmarked five CNV detection tools-PennCNV, QuantiSNP, iPattern, EnsembleCNV, and R-GADA-using SNP array and WGS data from 2002 individuals of the DRAGEN re-analysis of the 1000 Genomes project. Results showed significant variability in tool performance. R-GADA had the highest recall but low precision, while PennCNV was the most reliable in terms of precision and F1 score. EnsembleCNV improved recall by combining multiple callers but increased false positives. Overall, current tools, including new methods, do not outperform PennCNV in precise CNV detection. Improved reference data and consensus on true positive CNV calls are necessary. This study provides valuable insights and scalable workflows for researchers selecting CNV detection methods in future studies.

CNV_MCD: Detection of copy number variations based on minimum covariance determinant using next-generation sequencing data

Copy number variation (CNV), a pivotal form of genomic structural variation, plays a critical role in the genetic diversity of cancer genomes. In numerous studies, the identification of CNVs is commonly approached as an issue of outlier detection. To address this, read depth (RD) signals for genomic segments are first extracted from next-generation sequencing (NGS) data. CNVs are detected by assigning outlier scores to genomic segments based on the distance between their RD signals and those of adjacent segments. However, the mean and covariance estimators of the global distribution commonly utilized for calculating distance are susceptible to the effect of CNVs, resulting in inaccurate CNV detection. To solve this problem, we introduce a new method, CNV_MCD, for detecting CNVs based on the minimum covariance determinant (MCD). CNV_MCD employs the MCD method to estimate the mean and covariance of the RD profile, circumventing the need for direct computation of these parameters and ensuring the minimization of the determinant of the covariance matrix. This approach enables the calculation of a robust distance for each genomic segment, which serves as an outlier score. Furthermore, we implement a fast median filtering to correct for baseline drift in the outlier scores and use a chi-squared approximation to determine the cutoff distance for CNVs. These enhancements facilitate the detection of small CNVs, establishing CNV_MCD as a complementary method for CNV detection in low-coverage sequencing data. Extensive experiments on both simulated and real datasets demonstrate that CNV_MCD outperforms other popular CNV detection methods. Overall, our method offers a more robust and reliable technique for CNV detection, playing a crucial role in elucidating the genetic mechanisms underlying complex diseases such as cancer.

Copy Number Variations in Hereditary Spastic Paraplegia-Related Genes: Evaluation of an Iranian Hereditary Spastic Paraplegia Cohort and Literature Review

Introduction: In human genetic disorders, copy number variations (CNVs) are considered a considerable underlying cause. CNVs are generally detected by array-based methods but can also be discovered by read-depth analysis of whole-exome sequencing (WES) data. We performed WES-based CNV identification in a cohort of 35 Iranian families with hereditary spastic paraplegia (HSP) patients. Methods: Thirty-five patients whose routine single-nucleotide variants (SNVs) and insertion/deletion analyses from exome data were unrevealing underwent a pipeline of CNV analysis using the read-depth detection method. Subsequently, a comprehensive search about the existence of CNVs in all 84 known HSP-causing genes was carried out in all reported HSP cases, so far. Results and Discussion: CNV analysis of exome data indicated that 1 patient harbored a heterozygous deletion in exon 17 of the SPAST gene. Multiplex ligation-dependent probe amplification analysis confirmed this deletion in the proband and his affected father. Literature review demonstrated that, to date, pathogenic CNVs have been identified in 30 out of 84 HSP-causing genes (∼36%). However, CNVs in only 17 of these genes were specifically associated with the HSP phenotype. Among them, CNVs were more common in L1CAM, PLP1, SPAST, SPG7, SPG11, and REEP1 genes. The identification of the CNV in 1 of our patients suggests that WES allows the detection of both SNVs and CNVs from a single method without additional costs and execution time. However, because of intrinsic issues of WES in the detection of large rearrangements, it may not yet be exploited to replace the CNV detection methods in standard clinical practice.

CNVABNN: An AdaBoost algorithm and neural networks-based detection of copy number variations from NGS data

Copy number variation (CNV) is a non-negligible structural variation on the genome. And next-generation sequencing (NGS) technology is widely used to detect CNVs due to the feature of high throughput and low cost on the whole genome. Based on the original MFCNV method, this paper proposes an improved CNV detection method, which is called CNVABNN. In comparison to the MFCNV method, CNVABNN has three advantages: (1) It adds detectable categories, and refines the categories of loss into hemi_loss and homo_loss. (2) It utilizes the idea of integrated learning. The AdaBoost algorithm is used as the core framework and neural networks are used as weak classifiers, then CNVABNN combines all of the weak classifiers into a strong classifier. The overall performance of CNV detection is improved by using the strong classifier. (3) The detection is optimized by predicting CNVs twice through neural networks and voting mechanisms. To evaluate the performance of CNVABNN, six existing detection methods are used for comparison. The experimental results show that CNVABNN achieves better results in terms of precision, sensitivity, and F1-score for both simulated and real samples.

Clinical Selection of Prenatal Diagnostic Techniques Following Positive Noninvasive Prenatal Screening Results in Southwest China.

Background: This study aims to evaluate prenatal diagnosis methods following positive noninvasive prenatal screening (NIPS) results. Methods: According to the positive noninvasive prenatal screening results, 926 pregnant women were divided into three groups: main target disease group (high risk for trisomy 21, trisomy 18, or trisomy 13), sex chromosome aneuploidy (SCA) group, and other chromosomal abnormalities group [abnormal Z-scores for chromosomes other than trisomy (T)21/T18/T13 or SCAs]. The verification methods and results were then retrospectively analysed. Results: In the main target disease group, the positive rate of chromosomal abnormalities confirmed by quantitative fluorescence polymerase chain reaction (QF-PCR) was 75.18% (212/282), which was not significantly different from that by karyotyping (79.36%, 173/218) and copy number variation (CNV) detection methods (71.43%, 65/91). The positive rate of additional findings confirmed by karyotyping and copy number variation detection methods in main target disease group was 0.46% (1/218) and 8.79% (8/91), respectively. The positive rate of chromosomal abnormalities confirmed by karyotyping and CNV detection methods were 27.11% (45/166) and 38.46% (95/247) in the SCA group and 4.17% (1/24) and 20% (36/180) in the other chromosomal abnormalities group, respectively. Fetal sex chromosome mosaicism was detected in 16.13% (20/124) of the confirmed SCA cases. There were no significant differences in the detection rates of chromosomal microarray analysis (CMA) and CNV sequencing (CNVseq) among the three groups (p > 0.05). Conclusion: QF-PCR can quickly and accurately identify aneuploidies following NIPS-positive results for common aneuploidy, and in combination with karyotyping and CNV detection techniques can provide more comprehensive results. With the NIPS-positive results for SCA or other abnormalities, CMA and CNVseq may have the same effect on increasing the detection rate. The addition of fluorescence in situ hybridization assay may help to identify true fetal mosaicism.

CNV-MEANN: A Neural Network and Mind Evolutionary Algorithm-Based Detection of Copy Number Variations From Next-Generation Sequencing Data.

Copy number variation (CNV), is defined as repetitions or deletions of genomic segments of 1 Kb to 5 Mb, and is a major trigger for human disease. The high-throughput and low-cost characteristics of next-generation sequencing technology provide the possibility of the detection of CNVs in the whole genome, and also greatly improve the clinical practicability of next-generation sequencing (NGS) testing. However, current methods for the detection of CNVs are easily affected by sequencing and mapping errors, and uneven distribution of reads. In this paper, we propose an improved approach, CNV-MEANN, for the detection of CNVs, involving changing the structure of the neural network used in the MFCNV method. This method has three differences relative to the MFCNV method: (1) it utilizes a new feature, mapping quality, to replace two features in MFCNV, (2) it considers the influence of the loss categories of CNV on disease prediction, and refines the output structure, and (3) it uses a mind evolutionary algorithm to optimize the backpropagation (neural network) neural network model, and calculates individual scores for each genome bin to predict CNVs. Using both simulated and real datasets, we tested the performance of CNV-MEANN and compared its performance with those of seven widely used CNV detection methods. Experimental results demonstrated that the CNV-MEANN approach outperformed other methods with respect to sensitivity, precision, and F1-score. The proposed method was able to detect many CNVs that other approaches could not, and it reduced the boundary bias. CNV-MEANN is expected to be an effective method for the analysis of changes in CNVs in the genome.

CNproScan: Hybrid CNV detection for bacterial genomes

Discovering copy number variation (CNV) in bacteria is not in the spotlight compared to the attention focused on CNV detection in eukaryotes. However, challenges arising from bacterial drug resistance bring further interest to the topic of CNV and its role in drug resistance. General CNV detection methods do not consider bacteria's features and there is space to improve detection accuracy. Here, we present a CNV detection method called CNproScan focused on bacterial genomes. CNproScan implements a hybrid approach and other bacteria-focused features and depends only on NGS data. We benchmarked our method and compared it to the previously published methods and we can resolve to achieve a higher detection rate together with providing other beneficial features, such as CNV classification. Compared with other methods, CNproScan can detect much shorter CNV events.

A Novel Computational Framework to Predict Disease-Related Copy Number Variations by Integrating Multiple Data Sources.

Copy number variation (CNV) may contribute to the development of complex diseases. However, due to the complex mechanism of path association and the lack of sufficient samples, understanding the relationship between CNV and cancer remains a major challenge. The unprecedented abundance of CNV, gene, and disease label data provides us with an opportunity to design a new machine learning framework to predict potential disease-related CNVs. In this paper, we developed a novel machine learning approach, namely, IHI-BMLLR (Integrating Heterogeneous Information sources with Biweight Mid-correlation and L1-regularized Logistic Regression under stability selection), to predict the CNV-disease path associations by using a data set containing CNV, disease state labels, and gene data. CNVs, genes, and diseases are connected through edges and then constitute a biological association network. To construct a biological network, we first used a self-adaptive biweight mid-correlation (BM) formula to calculate correlation coefficients between CNVs and genes. Then, we used logistic regression with L1 penalty (LLR) function to detect genes related to disease. We added stability selection strategy, which can effectively reduce false positives, when using self-adaptive BM and LLR. Finally, a weighted path search algorithm was applied to find top D path associations and important CNVs. The experimental results on both simulation and prostate cancer data show that IHI-BMLLR is significantly better than two state-of-the-art CNV detection methods (i.e., CCRET and DPtest) under false-positive control. Furthermore, we applied IHI-BMLLR to prostate cancer data and found significant path associations. Three new cancer-related genes were discovered in the paths, and these genes need to be verified by biological research in the future.

A Cluster-Based Approach for the Discovery of Copy Number Variations From Next-Generation Sequencing Data.

The next-generation sequencing technology offers a wealth of data resources for the detection of copy number variations (CNVs) at a high resolution. However, it is still challenging to correctly detect CNVs of different lengths. It is necessary to develop new CNV detection tools to meet this demand. In this work, we propose a new CNV detection method, called CBCNV, for the detection of CNVs of different lengths from whole genome sequencing data. CBCNV uses a clustering algorithm to divide the read depth segment profile, and assigns an abnormal score to each read depth segment. Based on the abnormal score profile, Tukey’s fences method is adopted in CBCNV to forecast CNVs. The performance of the proposed method is evaluated on simulated data sets, and is compared with those of several existing methods. The experimental results prove that the performance of CBCNV is better than those of several existing methods. The proposed method is further tested and verified on real data sets, and the experimental results are found to be consistent with the simulation results. Therefore, the proposed method can be expected to become a routine tool in the analysis of CNVs from tumor-normal matched samples.

Copy Number Variant Detection with Low-Coverage Whole-Genome Sequencing Represents a Viable Alternative to the Conventional Array-CGH.

Copy number variations (CNVs) represent a type of structural variant involving alterations in the number of copies of specific regions of DNA that can either be deleted or duplicated. CNVs contribute substantially to normal population variability, however, abnormal CNVs cause numerous genetic disorders. At present, several methods for CNV detection are applied, ranging from the conventional cytogenetic analysis, through microarray-based methods (aCGH), to next-generation sequencing (NGS). In this paper, we present GenomeScreen, an NGS-based CNV detection method for low-coverage, whole-genome sequencing. We determined the theoretical limits of its accuracy and obtained confirmation in an extensive in silico study and in real patient samples with known genotypes. In theory, at least 6 M uniquely mapped reads are required to detect a CNV with the length of 100 kilobases (kb) or more with high confidence (Z-score > 7). In practice, the in silico analysis required at least 8 M to obtain >99% accuracy (for 100 kb deviations). We compared GenomeScreen with one of the currently used aCGH methods in diagnostic laboratories, which has mean resolution of 200 kb. GenomeScreen and aCGH both detected 59 deviations, while GenomeScreen furthermore detected 134 other (usually) smaller variations. When compared to aCGH, overall performance of the proposed GenemoScreen tool is comparable or superior in terms of accuracy, turn-around time, and cost-effectiveness, thus providing reasonable benefits, particularly in a prenatal diagnosis setting.

Copy-number variation contributes 9% of pathogenicity in the inherited retinal degenerations

PurposeCurrent sequencing strategies can genetically solve 55–60% of inherited retinal degeneration (IRD) cases, despite recent progress in sequencing. This can partially be attributed to elusive pathogenic variants (PVs) in known IRD genes, including copy-number variations (CNVs), which have been shown as major contributors to unsolved IRD cases. MethodsFive hundred IRD patients were analyzed with targeted next-generation sequencing (NGS). The NGS data were used to detect CNVs with ExomeDepth and gCNV and the results were compared with CNV detection with a single-nucleotide polymorphism (SNP) array. Likely causal CNV predictions were validated by quantitative polymerase chain reaction (qPCR). ResultsLikely disease-causing single-nucleotide variants (SNVs) and small indels were found in 55.6% of subjects. PVs in USH2A (11.6%), RPGR (4%), and EYS (4%) were the most common. Likely causal CNVs were found in an additional 8.8% of patients. Of the three CNV detection methods, gCNV showed the highest accuracy. Approximately 30% of unsolved subjects had a single likely PV in a recessive IRD gene. ConclusionCNV detection using NGS-based algorithms is a reliable method that greatly increases the genetic diagnostic rate of IRDs. Experimentally validating CNVs helps estimate the rate at which IRDs might be solved by a CNV plus a more elusive variant.

CNV Radar: an improved method for somatic copy number alteration characterization in oncology

BackgroundCancer associated copy number variation (CNV) events provide important information for identifying patient subgroups and suggesting treatment strategies. Technical and logistical issues, however, make it challenging to accurately detect abnormal copy number events in a cost-effective manner in clinical studies.ResultsHere we present CNV Radar, a software tool that utilizes next-generation sequencing read depth information and variant allele frequency patterns, to infer the true copy number status of genes and genomic regions from whole exome sequencing data. Evaluation of CNV Radar in a public multiple myeloma dataset demonstrated that CNV Radar was able to detect a variety of CNVs associated with risk of progression, and we observed > 70% concordance with fluorescence in situ hybridization (FISH) results. Compared to other CNV callers, CNV Radar showed high sensitivity and specificity. Similar results were observed when comparing CNV Radar calls to single nucleotide polymorphism array results from acute myeloid leukemia and prostate cancer datasets available on TCGA. More importantly, CNV Radar demonstrated its utility in the clinical trial setting: in POLLUX and CASTOR, two phase 3 studies in patients with relapsed or refractory multiple myeloma, we observed a high concordance rate with FISH for del17p, a risk defining CNV event (88% in POLLUX and 90% in CASTOR), therefore allowing for efficacy assessments in clinically relevant disease subgroups. Our case studies also showed that CNV Radar is capable of detecting abnormalities such as copy-neutral loss of heterozygosity that elude other approaches.ConclusionsWe demonstrated that CNV Radar is more sensitive than other CNV detection methods, accurately detects clinically important cytogenetic events, and allows for further interrogation of novel disease biology. Overall, CNV Radar exhibited high concordance with standard methods such as FISH, and its success in the POLLUX and CASTOR clinical trials demonstrated its potential utility for informing clinical and therapeutic decisions.

Copy Number Variation Detection Using Total Variation.

Next-generation sequencing (NGS) technologies offer new opportunities for precise and accurate identification of genomic aberrations, including copy number variations (CNVs). For high-throughput NGS data, using depth of coverage has become a major approach to identify CNVs, especially for whole exome sequencing (WES) data. Due to the high level of noise and biases of read-count data and complexity of the WES data, existing CNV detection tools identify many false CNV segments. Besides, NGS generates a huge amount of data, requiring to use effective and efficient methods. In this work, we propose a novel segmentation algorithm based on the total variation approach to detect CNVs more precisely and efficiently using WES data. The proposed method also filters out outlier read-counts and identifies significant change points to reduce false positives. We used real and simulated data to evaluate the performance of the proposed method and compare its performance with those of other commonly used CNV detection methods. Using simulated and real data, we show that the proposed method outperforms the existing CNV detection methods in terms of accuracy and false discovery rate and has a faster runtime compared to the circular binary segmentation method.

Vient de paraître

Current sequencing strategies can genetically solve 55–60% of inherited retinal degeneration (IRD) cases, despite recent progress in sequencing. This can partially be attributed to elusive pathogenic variants (PVs) in known IRD genes, including copy-number variations (CNVs), which have been shown as major contributors to unsolved IRD cases.Five hundred IRD patients were analyzed with targeted next-generation sequencing (NGS). The NGS data were used to detect CNVs with ExomeDepth and gCNV and the results were compared with CNV detection with a single-nucleotide polymorphism (SNP) array. Likely causal CNV predictions were validated by quantitative polymerase chain reaction (qPCR).Likely disease-causing single-nucleotide variants (SNVs) and small indels were found in 55.6% of subjects. PVs in USH2A (11.6%), RPGR (4%), and EYS (4%) were the most common. Likely causal CNVs were found in an additional 8.8% of patients. Of the three CNV detection methods, gCNV showed the highest accuracy. Approximately 30% of unsolved subjects had a single likely PV in a recessive IRD gene.CNV detection using NGS-based algorithms is a reliable method that greatly increases the genetic diagnostic rate of IRDs. Experimentally validating CNVs helps estimate the rate at which IRDs might be solved by a CNV plus a more elusive variant.

BagGMM: Calling copy number variation by bagging multiple Gaussian mixture models from tumor and matched normal next-generation sequencing data

Copy number variations (CNVs) contribute significantly to human genomic variability, some of which lead to diseases. However, effective detection of CNVs from whole genome next generation sequencing data (NGS) remains challenging. Here, we present BagGMM, a new method to call CNVs using tumor-normal matched samples from NGS data. BagGMM extracts read depth ratios of tumor samples to normal samples, divides the genomic sequences into segments by sliding windows to count the average coverage ratio of each segment, filters candidate deletions and duplications based on a coarse criterion of coverage ratio, and then builds Gaussian Mixture Model (GMM) for remaining ratios to identify the remaining ambiguous copy number states after filtration. Bagging multiple GMMs makes false positive calls descent instead of using one GMM, thus enhancing the detection power of BagGMM. Considering the computation speed of GMMs and false positive calls, we employ a segmentation procedure “large window and then small windows”, which is also helpful to determine boundary of CNV regions. We apply BagGMM to three simulation datasets and two groups of human whole genome sequencing (WGS) data for breast cancer patients and ovarian cancer patients to identify CNVs, respectively. All performed experiments demonstrate that BagGMM has the capability of robustly identification of CNVs with different sizes and states. The performance of this tool is compared to four peer existing CNV detection methods. BagGMM shows a significant improvement in both sensitivity and specificity for detecting both copy number gains and losses.

Noise cancellation using total variation for copy number variation detection

BackgroundDue to recent advances in sequencing technologies, sequence-based analysis has been widely applied to detecting copy number variations (CNVs). There are several techniques for identifying CNVs using next generation sequencing (NGS) data, however methods employing depth of coverage or read depth (RD) have recently become a main technique to identify CNVs. The main assumption of the RD-based CNV detection methods is that the readcount value at a specific genomic location is correlated with the copy number at that location. However, readcount data’s noise and biases distort the association between the readcounts and copy numbers. For more accurate CNV identification, these biases and noise need to be mitigated. In this work, to detect CNVs more precisely and efficiently we propose a novel denoising method based on the total variation approach and the Taut String algorithm.ResultsTo investigate the performance of the proposed denoising method, we computed sensitivities, false discovery rates and specificities of CNV detection when employing denoising, using both simulated and real data. We also compared the performance of the proposed denoising method, Taut String, with that of the commonly used approaches such as moving average (MA) and discrete wavelet transforms (DWT) in terms of sensitivity of detecting true CNVs and time complexity. The results show that Taut String works better than DWT and MA and has a better power to identify very narrow CNVs. The ability of Taut String denoising in preserving CNV segments’ breakpoints and narrow CNVs increases the detection accuracy of segmentation algorithms, resulting in higher sensitivities and lower false discovery rates.ConclusionsIn this study, we proposed a new denoising method for sequence-based CNV detection based on a signal processing technique. Existing CNV detection algorithms identify many false CNV segments and fail in detecting short CNV segments due to noise and biases. Employing an effective and efficient denoising method can significantly enhance the detection accuracy of the CNV segmentation algorithms. Advanced denoising methods from the signal processing field can be employed to implement such algorithms. We showed that non-linear denoising methods that consider sparsity and piecewise constant characteristics of CNV data result in better performance in CNV detection.

A remark on copy number variation detection methods.

Copy number variations (CNVs) are gain and loss of DNA sequence of a genome. High throughput platforms such as microarrays and next generation sequencing technologies (NGS) have been applied for genome wide copy number losses. Although progress has been made in both approaches, the accuracy and consistency of CNV calling from the two platforms remain in dispute. In this study, we perform a deep analysis on copy number losses on 254 human DNA samples, which have both SNP microarray data and NGS data publicly available from Hapmap Project and 1000 Genomes Project respectively. We show that the copy number losses reported from Hapmap Project and 1000 Genome Project only have < 30% overlap, while these reports are required to have cross-platform (e.g. PCR, microarray and high-throughput sequencing) experimental supporting by their corresponding projects, even though state-of-art calling methods were employed. On the other hand, copy number losses are found directly from HapMap microarray data by an accurate algorithm, i.e. CNVhac, almost all of which have lower read mapping depth in NGS data; furthermore, 88% of which can be supported by the sequences with breakpoint in NGS data. Our results suggest the ability of microarray calling CNVs and the possible introduction of false negatives from the unessential requirement of the additional cross-platform supporting. The inconsistency of CNV reports from Hapmap Project and 1000 Genomes Project might result from the inadequate information containing in microarray data, the inconsistent detection criteria, or the filtration effect of cross-platform supporting. The statistical test on CNVs called from CNVhac show that the microarray data can offer reliable CNV reports, and majority of CNV candidates can be confirmed by raw sequences. Therefore, the CNV candidates given by a good caller could be highly reliable without cross-platform supporting, so additional experimental information should be applied in need instead of necessarily.

Detection of copy number variations by pair analysis using next-generation sequencing data in inherited kidney diseases.

Comprehensive genetic approaches for diagnosing inherited kidney diseases using next-generation sequencing (NGS) have recently been established. However, even with these approaches, we are still failing to detect gene defects in some patients who appear to suffer from genetic diseases. One of the reasons for this is the difficulty of detecting copy number variations (CNVs) using our current approaches. For such cases, we can apply methods of array-based comparative genomic hybridization (aCGH) or multiplex ligation and probe amplification (MLPA); however, these are expensive and laborious and also often fail to identify CNVs. Here, we report seven cases with CNVs in various inherited kidney diseases screened by NGS pair analysis. Targeted sequencing analysis for causative genes was conducted for cases with suspected inherited kidney diseases, for some of which a definitive genetic diagnosis had not been achieved. We conducted pair analysis using NGS data for those cases. When CNVs were detected by pair analysis, they were confirmed by aCGH and/or MLPA. In seven cases, CNVs in various causative genes of inherited kidney diseases were detected by pair analysis. With aCGH and/or MLPA, pathogenic CNV variants were confirmed: COL4A5 or HNF1B in two cases each, and EYA1, CLCNKB, or PAX2 in one each. We presented seven cases with CNVs in various genes that were screened by pair analysis. The NGS-based CNV detection method is useful for comprehensive screening of CNVs, and our results revealed that, for a certain proportion of cases, CNV analysis is necessary for accurate genetic diagnosis.

CNVcaller: highly efficient and widely applicable software for detecting copy number variations in large populations.

BackgroundThe increasing amount of sequencing data available for a wide variety of species can be theoretically used for detecting copy number variations (CNVs) at the population level. However, the growing sample sizes and the divergent complexity of nonhuman genomes challenge the efficiency and robustness of current human-oriented CNV detection methods.ResultsHere, we present CNVcaller, a read-depth method for discovering CNVs in population sequencing data. The computational speed of CNVcaller was 1–2 orders of magnitude faster than CNVnator and Genome STRiP for complex genomes with thousands of unmapped scaffolds. CNV detection of 232 goats required only 1.4 days on a single compute node. Additionally, the Mendelian consistency of sheep trios indicated that CNVcaller mitigated the influence of high proportions of gaps and misassembled duplications in the nonhuman reference genome assembly. Furthermore, multiple evaluations using real sheep and human data indicated that CNVcaller achieved the best accuracy and sensitivity for detecting duplications.ConclusionsThe fast generalized detection algorithms included in CNVcaller overcome prior computational barriers for detecting CNVs in large-scale sequencing data with complex genomic structures. Therefore, CNVcaller promotes population genetic analyses of functional CNVs in more species.

An evaluation of copy number variation detection tools for cancer using whole exome sequencing data

BackgroundRecently copy number variation (CNV) has gained considerable interest as a type of genomic/genetic variation that plays an important role in disease susceptibility. Advances in sequencing technology have created an opportunity for detecting CNVs more accurately. Recently whole exome sequencing (WES) has become primary strategy for sequencing patient samples and study their genomics aberrations. However, compared to whole genome sequencing, WES introduces more biases and noise that make CNV detection very challenging. Additionally, tumors’ complexity makes the detection of cancer specific CNVs even more difficult. Although many CNV detection tools have been developed since introducing NGS data, there are few tools for somatic CNV detection for WES data in cancer.ResultsIn this study, we evaluated the performance of the most recent and commonly used CNV detection tools for WES data in cancer to address their limitations and provide guidelines for developing new ones. We focused on the tools that have been designed or have the ability to detect cancer somatic aberrations. We compared the performance of the tools in terms of sensitivity and false discovery rate (FDR) using real data and simulated data. Comparative analysis of the results of the tools showed that there is a low consensus among the tools in calling CNVs. Using real data, tools show moderate sensitivity (~50% - ~80%), fair specificity (~70% - ~94%) and poor FDRs (~27% - ~60%). Also, using simulated data we observed that increasing the coverage more than 10× in exonic regions does not improve the detection power of the tools significantly.ConclusionsThe limited performance of the current CNV detection tools for WES data in cancer indicates the need for developing more efficient and precise CNV detection methods. Due to the complexity of tumors and high level of noise and biases in WES data, employing advanced novel segmentation, normalization and de-noising techniques that are designed specifically for cancer data is necessary. Also, CNV detection development suffers from the lack of a gold standard for performance evaluation. Finally, developing tools with user-friendly user interfaces and visualization features can enhance CNV studies for a broader range of users.