SCNVSim: somatic copy number variation and structure variation simulator

Background

Somatically acquired structure variations (SVs) and copy number variations (CNVs) can induce genetic changes that are directly related to tumor genesis. Somatic SV/CNV detection using next-generation sequencing (NGS) data still faces major challenges introduced by tumor sample characteristics, such as ploidy, heterogeneity, and purity. A simulated cancer genome with known SVs and CNVs can serve as a benchmark for evaluating the performance of existing somatic SV/CNV detection tools and developing new methods.

Results

SCNVSim is a tool for simulating somatic CNVs and structure variations SVs. Other than multiple types of SV and CNV events, the tool is capable of simulating important features related to tumor samples including aneuploidy, heterogeneity and purity.

Conclusions

SCNVSim generates the genomes of a cancer cell population with detailed information of copy number status, loss of heterozygosity (LOH), and event break points, which is essential for developing and evaluating somatic CNV and SV detection methods in cancer genomics studies.     

Comparison of Sequencing Based CNV Discovery Methods Using Monozygotic Twin Quartets

Background

The advent of high throughput sequencing methods breeds an important amount of technical challenges. Among those is the one raised by the discovery of copy-number variations (CNVs) using whole-genome sequencing data. CNVs are genomic structural variations defined as a variation in the number of copies of a large genomic fragment, usually more than one kilobase. Here, we aim to compare different CNV calling methods in order to assess their ability to consistently identify CNVs by comparison of the calls in 9 quartets of identical twin pairs. The use of monozygotic twins provides a means of estimating the error rate of each algorithm by observing CNVs that are inconsistently called when considering the rules of Mendelian inheritance and the assumption of an identical genome between twins. The similarity between the calls from the different tools and the advantage of combining call sets were also considered.

Results

ERDS and CNVnator obtained the best performance when considering the inherited CNV rate with a mean of 0.74 and 0.70, respectively. Venn diagrams were generated to show the agreement between the different algorithms, before and after filtering out familial inconsistencies. This filtering revealed a high number of false positives for CNVer and Breakdancer. A low overall agreement between the methods suggested a high complementarity of the different tools when calling CNVs. The breakpoint sensitivity analysis indicated that CNVnator and ERDS achieved better resolution of CNV borders than the other tools. The highest inherited CNV rate was achieved through the intersection of these two tools (81%).

Conclusions

This study showed that ERDS and CNVnator provide good performance on whole genome sequencing data with respect to CNV consistency across families, CNV breakpoint resolution and CNV call specificity. The intersection of the calls from the two tools would be valuable for CNV genotyping pipelines.     

拷贝数变异的全基因组关联分析

基因组拷贝数变异(copy number variations,CNVs)是指与基因组参考序列相比,基因组中≥1 kb的DNA片段插入、缺失和/或扩增,及其互相组合衍生出的复杂变异．由于其具有分布范围广、可遗传、相对稳定和高度异质性等特点,目前认为,CNVs是一种新的可以作为疾病易感标志的基因组DNA多态性,其变异引起的基因剂量改变可以导致表型改变．最近,一种基于CNVs的新的疾病易感基因鉴定策略——CNV全基因组关联分析开始出现,这一策略和传统的基于单核苷酸多态性的关联分析具有互补性,通过认识基因组结构变异可以认识复杂疾病的分子机制和遗传基础．     

Genome-Wide Identification of Somatic Aberrations from Paired Normal-Tumor Samples

Genomic copy number alteration and allelic imbalance are distinct features of cancer cells, and recent advances in the genotyping technology have greatly boosted the research in the cancer genome. However, the complicated nature of tumor usually hampers the dissection of the SNP arrays. In this study, we describe a bioinformatic tool, named GIANT, for genome-wide identification of somatic aberrations from paired normal-tumor samples measured with SNP arrays. By efficiently incorporating genotype information of matched normal sample, it accurately detects different types of aberrations in cancer genome, even for aneuploid tumor samples with severe normal cell contamination. Furthermore, it allows for discovery of recurrent aberrations with critical biological properties in tumorigenesis by using statistical significance test. We demonstrate the superior performance of the proposed method on various datasets including tumor replicate pairs, simulated SNP arrays and dilution series of normal-cancer cell lines. Results show that GIANT has the potential to detect the genomic aberration even when the cancer cell proportion is as low as 5∼10%. Application on a large number of paired tumor samples delivers a genome-wide profile of the statistical significance of the various aberrations, including amplification, deletion and LOH. We believe that GIANT represents a powerful bioinformatic tool for interpreting the complex genomic aberration, and thus assisting both academic study and the clinical treatment of cancer.     

Genome-wide mapping of copy number variation in humans: comparative analysis of high resolution array platforms

Accurate and efficient genome-wide detection of copy number variants (CNVs) is essential for understanding human genomic variation, genome-wide CNV association type studies, cytogenetics research and diagnostics, and independent validation of CNVs identified from sequencing based technologies. Numerous, array-based platforms for CNV detection exist utilizing array Comparative Genome Hybridization (aCGH), Single Nucleotide Polymorphism (SNP) genotyping or both. We have quantitatively assessed the abilities of twelve leading genome-wide CNV detection platforms to accurately detect Gold Standard sets of CNVs in the genome of HapMap CEU sample NA12878, and found significant differences in performance. The technologies analyzed were the NimbleGen 4.2 M, 2.1 M and 3×720 K Whole Genome and CNV focused arrays, the Agilent 1×1 M CGH and High Resolution and 2×400 K CNV and SNP+CGH arrays, the Illumina Human Omni1Quad array and the Affymetrix SNP 6.0 array. The Gold Standards used were a 1000 Genomes Project sequencing-based set of 3997 validated CNVs and an ultra high-resolution aCGH-based set of 756 validated CNVs. We found that sensitivity, total number, size range and breakpoint resolution of CNV calls were highest for CNV focused arrays. Our results are important for cost effective CNV detection and validation for both basic and clinical applications.     

基因组拷贝数变异研究进展

基因组结构变异分为两个层次:显微水平(microscopic)和亚显微水平(submicroscopic)。显微水平的基因组结构变异主要是指显微镜下可见的染色体畸变,包括整倍体或非整倍体、缺失、插入、倒位、易位、脆性位点等结构变异。亚显微水平的基因组结构变异是指DNA片段长度在1Kb-3Mb的基因组结构变异,包括缺失、插入、重复、重排、倒位、DNA拷贝数目变化(copy numbervariation,CNV),这些统称为CNV或者CNP(copy number polymorphisms,CNP)。对CNV的研究能够帮助研究者建立遗传检测假说,进而发现疾病易感基因,同时加深对表型变异的理解,为今后研究人类生物功能、进化、疾病奠定基础。本文主要从CNV的研究历史、分子机制、研究方法、研究意义等四个方面进行综述.。     

SNP rs10248565 in HDAC9 as a novel genomic aberration biomarker of lung adenocarcinoma in non-smoking women

Background

Numerous efforts have been made to elucidate the etiology and improve the treatment of lung cancer, but the overall five-year survival rate is still only 15%. Although cigarette smoking is the primary risk factor for lung cancer, only 7% of female lung cancer patients in Taiwan have a history of smoking. Since cancer results from progressive accumulation of genetic aberrations, genomic rearrangements may be early events in carcinogenesis.

Results

In order to identify biomarkers of early-stage adenocarcinoma, the genome-wide DNA aberrations of 60 pairs of lung adenocarcinoma and adjacent normal lung tissue in non-smoking women were examined using Affymetrix Genome-Wide Human SNP 6.0 arrays. Common copy number variation (CNV) regions were identified by ≥30% of patients with copy number beyond 2 ± 0.5 of copy numbers for each single nucleotide polymorphism (SNP) and at least 100 continuous SNP variant loci. SNPs associated with lung adenocarcinoma were identified by McNemar’s test. Loss of heterozygosity (LOH) SNPs were identified in ≥18% of patients with LOH in the locus. Aberration of SNP rs10248565 at HDAC9 in chromosome 7p21.1 was identified from concurrent analyses of CNVs, SNPs, and LOH.

Conclusion

The results elucidate the genetic etiology of lung adenocarcinoma by demonstrating that SNP rs10248565 may be a potential biomarker of cancer susceptibility.     

Progress in the genomics and genome-wide study of sake yeast

ABSTRACT

Completion of the whole genome sequence of a laboratory yeast strain Saccharomyces cerevisiae in 1996 ushered in the development of genome-wide experimental tools and accelerated subsequent genetic study of S. cerevisiae. The study of sake yeast also shared the benefit of such tools as DNA microarrays, gene disruption-mutant collections, and others. Moreover, whole genome analysis of representative sake yeast strain Kyokai no. 7 was performed in the late 2000s, and enabled comparative genomics between sake yeast and laboratory yeast, resulting in some notable finding for of sake yeast genetics. Development of next-generation DNA sequencing and bioinformatics also drastically changed the field of the genetics, including for sake yeast. Genomics and the genome-wide study of sake yeast have progressed under these circumstances during the last two decades, and are summarized in this article.

Abbreviations: AFLP: amplified fragment length polymorphism; CGH: comparative genomic hybridization; CNV: copy number variation; DMS: dimethyl succinate; DSW: deep sea water; LOH: loss of heterozygosity; NGS: next generation sequencer; QTL: quantitative trait loci; QTN: quantitative trait nucleotide; SAM: S-adenosyl methionine; SNV: single nucleotide variation     

Allelic genome structural variations in maize detected by array comparative genome hybridization

DNA polymorphisms such as insertion/deletions and duplications affecting genome segments larger than 1 kb are known as copy-number variations (CNVs) or structural variations (SVs). They have been recently studied in animals and humans by using array-comparative genome hybridization (aCGH), and have been associated with several human diseases. Their presence and phenotypic effects in plants have not been investigated on a genomic scale, although individual structural variations affecting traits have been described. We used aCGH to investigate the presence of CNVs in maize by comparing the genome of 13 maize inbred lines to B73. Analysis of hybridization signal ratios of 60,472 60-mer oligonucleotide probes between inbreds in relation to their location in the reference genome (B73) allowed us to identify clusters of probes that deviated from the ratio expected for equal copy-numbers. We found CNVs distributed along the maize genome in all chromosome arms. They occur with appreciable frequency in different germplasm subgroups, suggesting ancient origin. Validation of several CNV regions showed both insertion/deletions and copy-number differences. The nature of CNVs detected suggests CNVs might have a considerable impact on plant phenotypes, including disease response and heterosis.     

High-resolution genome-wide analysis of irradiated (UV and γ-rays) diploid yeast cells reveals a high frequency of genomic loss of heterozygosity (LOH) events

In diploid eukaryotes, repair of double-stranded DNA breaks by homologous recombination often leads to loss of heterozygosity (LOH). Most previous studies of mitotic recombination in Saccharomyces cerevisiae have focused on a single chromosome or a single region of one chromosome at which LOH events can be selected. In this study, we used two techniques (single-nucleotide polymorphism microarrays and high-throughput DNA sequencing) to examine genome-wide LOH in a diploid yeast strain at a resolution averaging 1 kb. We examined both selected LOH events on chromosome V and unselected events throughout the genome in untreated cells and in cells treated with either γ-radiation or ultraviolet (UV) radiation. Our analysis shows the following: (1) spontaneous and damage-induced mitotic gene conversion tracts are more than three times larger than meiotic conversion tracts, and conversion tracts associated with crossovers are usually longer and more complex than those unassociated with crossovers; (2) most of the crossovers and conversions reflect the repair of two sister chromatids broken at the same position; and (3) both UV and γ-radiation efficiently induce LOH at doses of radiation that cause no significant loss of viability. Using high-throughput DNA sequencing, we also detected new mutations induced by γ-rays and UV. To our knowledge, our study represents the first high-resolution genome-wide analysis of DNA damage-induced LOH events performed in any eukaryote.     

Increased paternal age and the influence on burden of genomic copy number variation in the general population

Genomic copy number variations (CNVs) and increased parental age are both associated with the risk to develop a variety of clinical neuropsychiatric disorders such as autism, schizophrenia and bipolar disorder. At the same time, it has been shown that the rate of transmitted de novo single nucleotide mutations is increased with paternal age. To address whether paternal age also affects the burden of structural genomic deletions and duplications, we examined various types of CNV burden in a large population sample from the Netherlands. Healthy participants with parental age information (n = 6,773) were collected at different University Medical Centers. CNVs were called with the PennCNV algorithm using Illumina genome-wide SNP array data. We observed no evidence in support of a paternal age effect on CNV load in the offspring. Our results were negative for global measures as well as several proxies for de novo CNV events in this unique sample. While recent studies suggest de novo single nucleotide mutation rate to be dominated by the age of the father at conception, our results strongly suggest that at the level of global CNV burden there is no influence of increased paternal age. While it remains possible that local genomic effects may exist for specific phenotypes, this study indicates that global CNV burden and increased father’s age may be independent disease risk factors.     

Copy number variation in the cattle genome

Copy number variations (CNVs) are gains and losses of genomic sequence greater than 50?bp between two individuals of a species. While single nucleotide polymorphisms (SNPs) are more frequent, CNVs impact a higher percentage of genomic sequence and have potentially greater effects, including the changing of gene structure and dosage, altering gene regulation and exposing recessive alleles. In particular, segmental duplications (SDs) were shown to be one of the catalysts and hotspots for CNV formation. Substantial progress has been made in understanding CNVs in mammals, especially in humans and rodents. CNVs have been shown to be important in both normal phenotypic variability and disease susceptibility. Recently, interest in CNV study has extended into domesticated animals, including cattle. Multiple genome-wide cattle CNV studies have been carried out using both microarray and next generation sequencing technologies. Integration of SD and CNV results with SNP and other datasets are beginning to reveal impacts of CNVs on cattle domestication, health, and production traits.     

Identification of copy number variation hotspots in human populations

Copy number variants (CNVs) in the human genome contribute to both Mendelian and complex traits as well as to genomic plasticity in evolution. The investigation of mutational rates of CNVs is critical to understanding genomic instability and the etiology of the copy number variation (CNV)-related traits. However, the evaluation of the CNV mutation rate at the genome level poses an insurmountable practical challenge that requires large samples and accurate typing. In this study, we show that an approximate estimation of the CNV mutation rate could be achieved by using the phylogeny information of flanking SNPs. This allows a genome-wide comparison of mutation rates between CNVs with the use of vast, readily available data of SNP genotyping. A total of 4187 CNV regions (CNVRs) previously identified in HapMap populations were investigated in this study. We showed that the mutation rates for the majority of these CNVRs are at the order of 10⁻⁵ per generation, consistent with experimental observations at individual loci. Notably, the mutation rates of 104 (2.5%) CNVRs were estimated at the order of 10⁻³ per generation; therefore, they were identified as potential hotspots. Additional analyses revealed that genome architecture at CNV loci has a potential role in inciting mutational hotspots in the human genome. Interestingly, 49 (47%) CNV hotspots include human genes, some of which are known to be functional CNV loci (e.g., CNVs of C4 and β-defensin causing autoimmune diseases and CNVs of HYDIN with implication in control of cerebral cortex size), implicating the important role of CNV in human health and evolution, especially in common and complex diseases.