Microarray analysis of copy number variation in single cells

We present a protocol for reliably detecting DNA copy number aberrations in a single human cell. Multiple displacement-amplified DNAs of a cell are hybridized to a 3,000-bacterial artificial chromosome (BAC) array and to an Affymetrix 250,000 (250K)-SNP array. Subsequent copy number calling is based on the integration of BAC probe-specific copy number probabilities that are estimated by comparing probe intensities with a single-cell whole-genome amplification (WGA) reference model for diploid chromosomes, as well as SNP copy number and loss-of-heterozygosity states estimated by hidden Markov models (HMM). All methods for detecting DNA copy number aberrations in single human cells have difficulty in confidently discriminating WGA artifacts from true genetic variants. Furthermore, some methods lack thorough validation for segmental DNA imbalance detection. Our protocol minimizes false-positive variant calling and enables uniparental isodisomy detection in single cells. Additionally, it provides quality assessment, allowing the exclusion of uninterpretable single-cell WGA samples. The protocol takes 5-7 d.     

The Performance of Whole Genome Amplification Methods and Next-Generation Sequencing for Pre-Implantation Genetic Diagnosis of Chromosomal Abnormalities

Reliable and accurate pre-implantation genetic diagnosis(PGD) of patient’s embryos by next-generation sequencing(NGS) is dependent on efficient whole genome amplification(WGA) of a representative biopsy sample. However, the performance of the current state of the art WGA methods has not been evaluated for sequencing. Using low template DNA(15 pg) and single cells, we showed that the two PCR-based WGA systems Sure Plex and MALBAC are superior to the REPLI-g WGA multiple displacement amplification(MDA) system in terms of consistent and reproducible genome coverage and sequence bias across the 24 chromosomes, allowing better normalization of test to reference sequencing data. When copy number variation sequencing(CNV-Seq) was applied to single cell WGA products derived by either Sure Plex or MALBAC amplification, we showed that known disease CNVs in the range of 3e15 Mb could be reliably and accurately detected at the correct genomic positions. These findings indicate that our CNV-Seq pipeline incorporating either Sure Plex or MALBAC as the key initial WGA step is a powerful methodology for clinical PGD to identify euploid embryos in a patient’s cohort for uterine transplantation.     

不同单细胞全基因组扩增方法的比较及MALBAC在辅助生殖中的应用

单细胞全基因组扩增(whole genome amplification, WGA)是指在单细胞水平对全基因组进行扩增的新技术,其原理是将分离的单个细胞的微量全基因组DNA进行扩增,获得高覆盖率的完整的基因组后进行高通量测序,用于揭示细胞异质性。目前,WGA方法主要包括引物延伸预扩增(primer extension preamplification PCR, PEP-PCR)、简并寡核苷酸引物PCR (degenerate oligonucleotide primed PCR, DOP-PCR)、多重置换扩增(multiple displacement amplification, MDA)、多次退火环状循环扩增(multiple annealing and looping-based amplification cycles, MALBAC)等。本文对不同的单细胞WGA方法的原理及应用情况分别进行了阐述,并对其扩增效率进行评价和比较,包括基因组覆盖度、均一性、重现性、SNV (single-nucleotide variants)和CNV (copy number variants)检测力等。综合对比不同单细胞WGA方法后发现,MALBAC的扩增均一性最高、等位基因脱扣率最低、重现性最好,且对于CNV和SNV的检测效果最好。本文还阐述了MALBAC技术在人类单精子减数重组、非整倍体分析以及人类卵细胞基因组研究中的应用。     

Single-cell paired-end genome sequencing reveals structural variation per cell cycle

The nature and pace of genome mutation is largely unknown. Because standard methods sequence DNA from populations of cells, the genetic composition of individual cells is lost, de novo mutations in cells are concealed within the bulk signal and per cell cycle mutation rates and mechanisms remain elusive. Although single-cell genome analyses could resolve these problems, such analyses are error-prone because of whole-genome amplification (WGA) artefacts and are limited in the types of DNA mutation that can be discerned. We developed methods for paired-end sequence analysis of single-cell WGA products that enable (i) detecting multiple classes of DNA mutation, (ii) distinguishing DNA copy number changes from allelic WGA-amplification artefacts by the discovery of matching aberrantly mapping read pairs among the surfeit of paired-end WGA and mapping artefacts and (iii) delineating the break points and architecture of structural variants. By applying the methods, we capture DNA copy number changes acquired over one cell cycle in breast cancer cells and in blastomeres derived from a human zygote after in vitro fertilization. Furthermore, we were able to discover and fine-map a heritable inter-chromosomal rearrangement t(1;16)(p36;p12) by sequencing a single blastomere. The methods will expedite applications in basic genome research and provide a stepping stone to novel approaches for clinical genetic diagnosis.     

What a difference copy number variation makes

DNA copy number variation (CNV) represents a considerable source of human genetic diversity. Recently,1 a global map of copy number variation in the human genome has been drawn up which reveals not only the ubiquity but also the complexity of this type of variation. Thus, two human genomes may differ by more than 20 Mb and it is likely that the full extent of CNV still remains to be discovered. Nearly 3000 genes are associated with CNV. This high degree of variability with regard to gene copy number between two individuals challenges definitions of normality. Many CNVs are located in regions of complex genomic structure and this currently limits the extent to which these variants can be genotyped by using tagging SNPs. However, some CNVs are already amenable to genome-wide association studies so that their influence on human phenotypic diversity and disease susceptibility may soon be determined.     

Use of whole genome amplification and comparative genomic hybridisation to detect chromosomal copy number alterations in cell line material and tumour tissue

We have established that whole genome amplification (WGA), in conjunction with genomic DNA array comparative genomic hybridisation (gaCGH) allows for the identification of genome-wide copy number abnormalities (CNAs) in DNA extracted from both cell line and patient material. To determine the fidelity and reproducibility of WGA to detect copy number imbalances using gaCGH, well characterized cell line genomic DNA was analysed. The gaCGH data obtained from non-amplified DNA and amplified DNA for the neuroblastoma cell line NUB7 and a paediatric medulloblastoma patient was almost identical. In addition, laser capture microdissection (LCM) of prostate tumour cells and subsequent WGA allowed for the detection of a number of CNAs that may not have been identified if DNA had been extracted in bulk from heterogeneous tissue. The results presented here demonstrate the use of WGA for generating sufficient DNA for gaCGH analysis without the introduction of significant sequence representation bias. The combination of amplification and gaCGH using DNA extracted from archival patient material has the potential for permitting the studying of DNA from small cancerous or pre-cancerous foci, which may help to identify potential genomic markers for early diagnosis.     

A Simple Method for Encapsulating Single Cells in Alginate Microspheres Allows for Direct PCR and Whole Genome Amplification

Microdroplets are an effective platform for segregating individual cells and amplifying DNA. However, a key challenge is to recover the contents of individual droplets for downstream analysis. This paper offers a method for embedding cells in alginate microspheres and performing multiple serial operations on the isolated cells. Rhodobacter sphaeroides cells were diluted in alginate polymer and sprayed into microdroplets using a fingertip aerosol sprayer. The encapsulated cells were lysed and subjected either to conventional PCR, or whole genome amplification using either multiple displacement amplification (MDA) or a two-step PCR protocol. Microscopic examination after PCR showed that the lumen of the occupied microspheres contained fluorescently stained DNA product, but multiple displacement amplification with phi29 produced only a small number of polymerase colonies. The 2-step WGA protocol was successful in generating fluorescent material, and quantitative PCR from DNA extracted from aliquots of microspheres suggested that the copy number inside the microspheres was amplified up to 3 orders of magnitude. Microspheres containing fluorescent material were sorted by a dilution series and screened with a fluorescent plate reader to identify single microspheres. The DNA was extracted from individual isolates, re-amplified with full-length sequencing adapters, and then a single isolate was sequenced using the Illumina MiSeq platform. After filtering the reads, the only sequences that collectively matched a genome in the NCBI nucleotide database belonged to R. sphaeroides. This demonstrated that sequencing-ready DNA could be generated from the contents of a single microsphere without culturing. However, the 2-step WGA strategy showed limitations in terms of low genome coverage and an uneven frequency distribution of reads across the genome. This paper offers a simple method for embedding cells in alginate microspheres and performing PCR on isolated cells in common bulk reactions, although further work must be done to improve the amplification coverage of single genomes.     

A Robust Method to Analyze Copy Number Alterations of Less than 100 kb in Single Cells Using Oligonucleotide Array CGH

Comprehensive genome wide analyses of single cells became increasingly important in cancer research, but remain to be a technically challenging task. Here, we provide a protocol for array comparative genomic hybridization (aCGH) of single cells. The protocol is based on an established adapter-linker PCR (WGAM) and allowed us to detect copy number alterations as small as 56 kb in single cells. In addition we report on factors influencing the success of single cell aCGH downstream of the amplification method, including the characteristics of the reference DNA, the labeling technique, the amount of input DNA, reamplification, the aCGH resolution, and data analysis. In comparison with two other commercially available non-linear single cell amplification methods, WGAM showed a very good performance in aCGH experiments. Finally, we demonstrate that cancer cells that were processed and identified by the CellSearch® System and that were subsequently isolated from the CellSearch® cartridge as single cells by fluorescence activated cell sorting (FACS) could be successfully analyzed using our WGAM-aCGH protocol. We believe that even in the era of next-generation sequencing, our single cell aCGH protocol will be a useful and (cost-) effective approach to study copy number alterations in single cells at resolution comparable to those reported currently for single cell digital karyotyping based on next generation sequencing data.     

Global copy number analyses by next generation sequencing provide insight into pig genome variation

Background

Copy number variations (CNVs) confer significant effects on genetic innovation and phenotypic variation. Previous CNV studies in swine seldom focused on in-depth characterization of global CNVs.

Results

Using whole-genome assembly comparison (WGAC) and whole-genome shotgun sequence detection (WSSD) approaches by next generation sequencing (NGS), we probed formation signatures of both segmental duplications (SDs) and individualized CNVs in an integrated fashion, building the finest resolution CNV and SD maps of pigs so far. We obtained copy number estimates of all protein-coding genes with copy number variation carried by individuals, and further confirmed two genes with high copy numbers in Meishan pigs through an enlarged population. We determined genome-wide CNV hotspots, which were significantly enriched in SD regions, suggesting evolution of CNV hotspots may be affected by ancestral SDs. Through systematically enrichment analyses based on simulations and bioinformatics analyses, we revealed CNV-related genes undergo a different selective constraint from those CNV-unrelated regions, and CNVs may be associated with or affect pig health and production performance under recent selection.

Conclusions

Our studies lay out one way for characterization of CNVs in the pig genome, provide insight into the pig genome variation and prompt CNV mechanisms studies when using pigs as biomedical models for human diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-593) contains supplementary material, which is available to authorized users.     

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.     

Impact of whole genome amplification on analysis of copy number variants

Large-scale copy number variants (CNVs) have recently been recognized to play a role in human genome variation and disease. Approaches for analysis of CNVs in small samples such as microdissected tissues can be confounded by limited amounts of material. To facilitate analyses of such samples, whole genome amplification (WGA) techniques were developed. In this study, we explored the impact of Phi29 multiple-strand displacement amplification on detection of CNVs using oligonucleotide arrays. We extracted DNA from fresh frozen lymph node samples and used this for amplification and analysis on the Affymetrix Mapping 500k SNP array platform. We demonstrated that the WGA procedure introduces hundreds of potentially confounding CNV artifacts that can obscure detection of bona fide variants. Our analysis indicates that many artifacts are reproducible, and may correlate with proximity to chromosome ends and GC content. Pair-wise comparison of amplified products considerably reduced the number of apparent artifacts and partially restored the ability to detect real CNVs. Our results suggest WGA material may be appropriate for copy number analysis when amplified samples are compared to similarly amplified samples and that only the CNVs with the greatest significance values detected by such comparisons are likely to be representative of the unamplified samples.     

Rapid visualisation of microarray copy number data for the detection of structural variations linked to a disease phenotype

Whilst the majority of inherited diseases have been found to be caused by single base substitutions, small insertions or deletions (<1Kb), a significant proportion of genetic variability is due to copy number variation (CNV). The possible role of CNV in monogenic and complex diseases has recently attracted considerable interest. However, until the development of whole genome, oligonucleotide micro-arrays, designed specifically to detect the presence of copy number variation, it was not easy to screen an individual for the presence of unknown deletions or duplications with sizes below the level of sensitivity of optical microscopy (3-5 Mb). Now that currently available oligonucleotide micro-arrays have in excess of a million probes, the problem of copy number analysis has moved from one of data production to that of data analysis. We have developed CNViewer, to identify copy number variation that co-segregates with a disease phenotype in small nuclear families, from genome-wide oligonucleotide micro-array data. This freely available program should constitute a useful addition to the diagnostic armamentarium of clinical geneticists.     

Concurrent Whole-Genome Haplotyping and Copy-Number Profiling of Single Cells

Methods for haplotyping and DNA copy-number typing of single cells are paramount for studying genomic heterogeneity and enabling genetic diagnosis. Before analyzing the DNA of a single cell by microarray or next-generation sequencing, a whole-genome amplification (WGA) process is required, but it substantially distorts the frequency and composition of the cell’s alleles. As a consequence, haplotyping methods suffer from error-prone discrete SNP genotypes (AA, AB, BB) and DNA copy-number profiling remains difficult because true DNA copy-number aberrations have to be discriminated from WGA artifacts. Here, we developed a single-cell genome analysis method that reconstructs genome-wide haplotype architectures as well as the copy-number and segregational origin of those haplotypes by employing phased parental genotypes and deciphering WGA-distorted SNP B-allele fractions via a process we coin haplarithmisis. We demonstrate that the method can be applied as a generic method for preimplantation genetic diagnosis on single cells biopsied from human embryos, enabling diagnosis of disease alleles genome wide as well as numerical and structural chromosomal anomalies. Moreover, meiotic segregation errors can be distinguished from mitotic ones.     

Copy Number Variations in Tilapia Genomes

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R ² > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.     

The Role of Copy Number Variation in Susceptibility to Amyotrophic Lateral Sclerosis: Genome-Wide Association Study and Comparison with Published Loci

Background

The genetic contribution to sporadic amyotrophic lateral sclerosis (ALS) has not been fully elucidated. There are increasing efforts to characterise the role of copy number variants (CNVs) in human diseases; two previous studies concluded that CNVs may influence risk of sporadic ALS, with multiple rare CNVs more important than common CNVs. A little-explored issue surrounding genome-wide CNV association studies is that of post-calling filtering and merging of raw CNV calls. We undertook simulations to define filter thresholds and considered optimal ways of merging overlapping CNV calls for association testing, taking into consideration possibly overlapping or nested, but distinct, CNVs and boundary estimation uncertainty.

Methodology and Principal Findings

In this study we screened Illumina 300K SNP genotyping data from 730 ALS cases and 789 controls for copy number variation. Following quality control filters using thresholds defined by simulation, a total of 11321 CNV calls were made across 575 cases and 621 controls. Using region-based and gene-based association analyses, we identified several loci showing nominally significant association. However, the choice of criteria for combining calls for association testing has an impact on the ranking of the results by their significance. Several loci which were previously reported as being associated with ALS were identified here. However, of another 15 genes previously reported as exhibiting ALS-specific copy number variation, only four exhibited copy number variation in this study. Potentially interesting novel loci, including EEF1D, a translation elongation factor involved in the delivery of aminoacyl tRNAs to the ribosome (a process which has previously been implicated in genetic studies of spinal muscular atrophy) were identified but must be treated with caution due to concerns surrounding genomic location and platform suitability.

Conclusions and Significance

Interpretation of CNV association findings must take into account the effects of filtering and combining CNV calls when based on early genome-wide genotyping platforms and modest study sizes.     

Genome-wide copy number variation in Hanwoo, Black Angus, and Holstein cattle

Hanwoo, Korean native cattle, is indigenous to the Korean peninsula. They have been used mainly as draft animals for about 5,000 years; however, in the last 30 years, their main role has been changed to meat production by selective breeding which has led to substantial increases in their productivity. Massively parallel sequencing technology has recently made possible the systematic identification of structural variations in cattle genomes. In particular, copy number variation (CNV) has been recognized as an important genetic variation complementary to single-nucleotide polymorphisms that can be used to account for variations of economically important traits in cattle. Here we report genome-wide copy number variation regions (CNVRs) in Hanwoo cattle obtained by comparing the whole genome sequence of Hanwoo with Black Angus and Holstein sequence datasets. We identified 1,173 and 963 putative CNVRs representing 16.7 and 7.8 Mbp from comparisons between Black Angus and Hanwoo and between Holstein and Hanwoo, respectively. The potential functional roles of the CNVRs were assessed by Gene Ontology enrichment analysis. The results showed that response to stimulus, immune system process, and cellular component organization were highly enriched in the genic-CNVRs that overlapped with annotated cattle genes. Of the 11 CNVRs that were selected for validation by quantitative real-time PCR, 9 exhibited the expected copy number differences. The results reported in this study show that genome-wide CNVs were detected successfully using massively parallel sequencing technology. The CNVs may be a valuable resource for further studies to correlate CNVs and economically important traits in cattle.     

Genome‐based estimates of fungal rDNA copy number variation across phylogenetic scales and ecological lifestyles

Ribosomal DNA (rDNA) copy number variation (CNV) has major physiological implications for all organisms, but how it varies for fungi, an ecologically ubiquitous and important group of microorganisms, has yet to be systemically investigated. Here, we examine rDNA CNV using an in silico read depth approach for 91 fungal taxa with sequenced genomes and assess copy number conservation across phylogenetic scales and ecological lifestyles. rDNA copy number varied considerably across fungi, ranging from an estimated 14 to 1,442 copies (mean = 113, median = 82), and copy number similarity was inversely correlated with phylogenetic distance. No correlations were found between rDNA CNV and fungal trophic mode, ecological guild or genome size. Taken together, these results show that like other microorganisms, fungi exhibit substantial variation in rDNA copy number, which is linked to their phylogeny in a scale‐dependent manner.     

Variation of CNV distribution in five different ethnic populations

Recent studies have revealed a new type of variation in the human genome encompassing relatively large genomic segments ( approximately 100 kb-2.5 Mb), commonly referred to as copy number variation (CNV). The full nature and extent of CNV and its frequency in different ethnic populations is still largely unknown. In this study we surveyed a set of 12 CNVs previously detected by array-CGH. More than 300 individuals from five different ethnic populations, including three distinct European, one Asian and one African population, were tested for the occurrence of CNV using multiplex ligation-dependent probe amplification (MLPA). Seven of these loci indeed showed CNV, i.e., showed copy numbers that deviated from the population median. More precise estimations of the actual genomic copy numbers for (part of) the NSF gene locus, revealed copy numbers ranging from two to at least seven. Additionally, significant inter-population differences in the distribution of these copy numbers were observed. These data suggest that insight into absolute DNA copy numbers for loci exhibiting CNV is required to determine their potential contribution to normal phenotypic variation and, in addition, disease susceptibility.     

Genomic and genetic aspects of autism spectrum disorder

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with a strong genetic component. The past decade has witnessed tremendous progress in the genetic studies of ASD. In this article, we review the accumulating literatures on the monogenic forms of ASD and chromosomal abnormalities associated with ASD, the genome-wide linkage and association studies, the copy number variation (CNV) and the next generation sequencing (NGS) studies. With more than hundreds of mutations being implicated, the convergent biological pathways are emerging and the genetic landscape of ASD becomes clearer. The genetic studies provide a solid basis for future translational study for better diagnoses, intervention and treatment of ASD.     

Whole-genome multiple displacement amplification from single cells

Multiple displacement amplification (MDA) is a recently described method of whole-genome amplification (WGA) that has proven efficient in the amplification of small amounts of DNA, including DNA from single cells. Compared with PCR-based WGA methods, MDA generates DNA with a higher molecular weight and shows better genome coverage. This protocol was developed for preimplantation genetic diagnosis, and details a method for performing single-cell MDA using the phi29 DNA polymerase. It can also be useful for the amplification of other minute quantities of DNA, such as from forensic material or microdissected tissue. The protocol includes the collection and lysis of single cells, and all materials and steps involved in the MDA reaction. The whole procedure takes 3 h and generates 1-2 microg of DNA from a single cell, which is suitable for multiple downstream applications, such as sequencing, short tandem repeat analysis or array comparative genomic hybridization.