Probe-free allele-specific copy number detection and analysis of tumors

Cancer development and progression frequently involve nucleotide mutations as well as amplifications and deletions of genomic segments. Quantification of allele-specific copy number is an important step in characterizing tumor genomes for precision medicine. Despite advances in approaches to high-throughput genomic DNA analysis, inexpensive and simple methods for analyzing complex nucleotide and copy number variants are still needed. Real-time polymerase chain reaction (PCR) methods for discovering and genotyping single nucleotide polymorphisms are becoming increasingly important in genetic analysis. In this study, we describe a simple, single-tube, probe-free method that combines SYBR Green I-based quantitative real-time PCR and quantitative melting curve analysis both to detect specific nucleotide variants and to quantify allele-specific copy number variants of tumors. The approach is based on the quantification of the targets of interest and the relative abundance of two alleles in a single tube. The specificity, sensitivity, and utility of the assay were demonstrated in detecting allele-specific copy number changes critical for carcinogenesis and therapeutic intervention. Our approach would be useful for allele-specific copy number analysis or precise genotyping.     

Two Distinct Categories of Focal Deletions in Cancer Genomes

One of the key questions about genomic alterations in cancer is whether they are functional in the sense of contributing to the selective advantage of tumor cells. The frequency with which an alteration occurs might reflect its ability to increase cancer cell growth, or alternatively, enhanced instability of a locus may increase the frequency with which it is found to be aberrant in tumors, regardless of oncogenic impact. Here we’ve addressed this on a genome-wide scale for cancer-associated focal deletions, which are known to pinpoint both tumor suppressor genes (tumor suppressors) and unstable loci. Based on DNA copy number analysis of over one-thousand human cancers representing ten different tumor types, we observed five loci with focal deletion frequencies above 5%, including the A2BP1 gene at 16p13.3 and the MACROD2 gene at 20p12.1. However, neither RNA expression nor functional studies support a tumor suppressor role for either gene. Further analyses suggest instead that these are sites of increased genomic instability and that they resemble common fragile sites (CFS). Genome-wide analysis revealed properties of CFS-like recurrent deletions that distinguish them from deletions affecting tumor suppressor genes, including their isolation at specific loci away from other genomic deletion sites, a considerably smaller deletion size, and dispersal throughout the affected locus rather than assembly at a common site of overlap. Additionally, CFS-like deletions have less impact on gene expression and are enriched in cell lines compared to primary tumors. We show that loci affected by CFS-like deletions are often distinct from known common fragile sites. Indeed, we find that each tumor tissue type has its own spectrum of CFS-like deletions, and that colon cancers have many more CFS-like deletions than other tumor types. We present simple rules that can pinpoint focal deletions that are not CFS-like and more likely to affect functional tumor suppressors.     

Sequence organisation in nuclear DNA from Physarum polycephalum: Genomic organisation of DNA segments containing foldback sequences

DNA clones containing foldback sequences, derived from Physarum polycephalum nuclear DNA, can be classified according to their pattern of hydridisation to Southern blots of genomic DNA. One group of DNA clones map to unique DNA loci when used as a probe to restriction digests of Physarum nuclear DNA. These cloned segments appear to contain dispersed repetitive sequence elements located at many hundreds of sites in the genome. Similar patterns of hybridisation are generated when these cloned DNA probes are annealed to DNA restriction fragments of genomic DNA obtained from a number of different Physarum strains, indicating that no detectable alteration has occurred at these genomic loci subsequent to the divergence of the strains as a result of the introduction or deletion of mobile genetic elements. However, deletion of segments of some cloned DNA fragments occurs following their propagation in Escherichia coli. A second, distinct group of clones are shown to be derived from highly methylated segments of Physarum DNA which contain very abundant repetitive sequences with regular, though complex, arrangements of restriction sites at their various genomic locations. It is suggested that these DNA segments contain clustered repetitive sequence elements. The results lead to the conclusion that foldback elements in Physarum DNA are located in segments of the genome which display markedly different patterns of sequence organisation and degree of DNA methylation.     

A quantitative assay to detect α-thalassemia deletions and triplications using multiplex nested real-time quantitative polymerase chain reaction

Increasing evidence indicates that copy number variants (CNVs) have great relevance to common human diseases. In α-thalassemia, clinical phenotypes are related to genotypes, specifically copy number changes in the human α-globin gene cluster. Assays are available for high-throughput screening of unknown CNVs genome-wide and also for targeted CNV genotyping at loci associated with genetic disorders. Here we describe a universal quantitative approach based on nested real-time quantitative polymerase chain reaction for accurate determination of copy numbers at multiple particular gene loci. We used the α-globin gene as a model system, obtaining the reproducibility and sensitivity to analyze different gene copies and testing 95 DNA samples with 16 different known genotypes. Our results showed that this approach has high sensitivity and low standard deviations for correctly genotyping DNA samples containing different copy numbers of the α1 and α2 globin genes. Our method is rapid, simple, and reliable, and it could be used to simultaneously screen for α-thalassemia deletions or triplications. Moreover, it has potential as a versatile technology for the rapid genotyping of known CNVs in a targeted region.     

Detailed chromosomal and molecular genetic analysis of single cells by whole genome amplification and comparative genomic hybridisation.

Molecular genetic analysis of isolated single cells and other minute DNA samples is limited because there is insufficient DNA to perform more than one independent PCR amplification. One solution to this problem is to first amplify the entire genome, thus providing enough DNA for numerous subsequent PCRs. In this study we have investigated four different methods of whole genome amplification performed on single cells, and have identified a protocol that generates sufficient quantities of DNA for comparative genomic hybridisation (CGH) as well as more than 90 independent amplification reactions. Thus, numerous specific loci and the copy number of every chromosome can be assessed in a single cell. We report here the first reliable application of CGH to single cells from human preimplantation embryos (blastomeres) and to single fibroblasts, buccal cells and amniocytes.     

A segmental maximum a posteriori approach to genome-wide copy number profiling

MOTIVATION: Copy number profiling methods aim at assigning DNA copy numbers to chromosomal regions using measurements from microarray-based comparative genomic hybridizations. Among the proposed methods to this end, Hidden Markov Model (HMM)-based approaches seem promising since DNA copy number transitions are naturally captured in the model. Current discrete-index HMM-based approaches do not, however, take into account heterogeneous information regarding the genomic overlap between clones. Moreover, the majority of existing methods are restricted to chromosome-wise analysis. RESULTS: We introduce a novel Segmental Maximum A Posteriori approach, SMAP, for DNA copy number profiling. Our method is based on discrete-index Hidden Markov Modeling and incorporates genomic distance and overlap between clones. We exploit a priori information through user-controllable parameterization that enables the identification of copy number deviations of various lengths and amplitudes. The model parameters may be inferred at a genome-wide scale to avoid overfitting of model parameters often resulting from chromosome-wise model inference. We report superior performances of SMAP on synthetic data when compared with two recent methods. When applied on our new experimental data, SMAP readily recognizes already known genetic aberrations including both large-scale regions with aberrant DNA copy number and changes affecting only single features on the array. We highlight the differences between the prediction of SMAP and the compared methods and show that SMAP accurately determines copy number changes and benefits from overlap consideration.     

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.     

Digital karyotyping: an update of its applications in cancer

DNA copy number alterations, including entire chromosomal changes and small interstitial DNA amplifications and deletions, characterize the development of cancer. These changes usually affect the expression of target genes and subsequently the function of the target proteins. Since the completion of the human genome project, the capacity to comprehensively analyze the human cancer genome has expanded significantly. Techniques such as digital karyotyping have been developed to allow for the detection of DNA copy number alterations in cancer at the whole-genome scale. When compared with conventional methods such as spectral karyotyping, representational difference analysis, comparative genomic hybridization (CGH), or the more recent array CGH; digital karyotyping provides an evaluation of copy number of genetic material at higher resolution. Digital karyotyping has therefore promised to enhance our understanding of the cancer genome. This article provides an overview of digital karyotyping including the principle of the technology and its applications in identifying potential oncogenes and tumor suppressor genes.     

Accurate, high-throughput typing of copy number variation using paralogue ratios from dispersed repeats

Recent work has demonstrated an unexpected prevalence of copy number variation in the human genome, and has highlighted the part this variation may play in predisposition to common phenotypes. Some important genes vary in number over a high range (e.g. DEFB4, which commonly varies between two and seven copies), and have posed formidable technical challenges for accurate copy number typing, so that there are no simple, cheap, high-throughput approaches suitable for large-scale screening. We have developed a simple comparative PCR method based on dispersed repeat sequences, using a single pair of precisely designed primers to amplify products simultaneously from both test and reference loci, which are subsequently distinguished and quantified via internal sequence differences. We have validated the method for the measurement of copy number at DEFB4 by comparison of results from >800 DNA samples with copy number measurements by MAPH/REDVR, MLPA and array-CGH. The new Paralogue Ratio Test (PRT) method can require as little as 10 ng genomic DNA, appears to be comparable in accuracy to the other methods, and for the first time provides a rapid, simple and inexpensive method for copy number analysis, suitable for application to typing thousands of samples in large case-control association studies.     

M-CGH: Analysing microarray-based CGH experiments

Background  

Microarray-based comparative genomic hybridisation (array CGH) is a technique by which variation in relative copy numbers between two genomes can be analysed by competitive hybridisation to DNA microarrays. This technology has most commonly been used to detect chromosomal amplifications and deletions in cancer. Dedicated tools are needed to analyse the results of such experiments, which include appropriate visualisation, and to take into consideration the physical relation in the genome between the probes on the array.     

Chromosome‐specific physical localisation of expressed sequence tag loci in Corchorus olitorius L.

Jute (Corchorus spp.), as a natural fibre‐producing species, ranks next only to cotton. Inadequate understanding of its genetic architecture is a major lacuna for genetic improvement of this crop in terms of yield and quality. Establishment of a physical map provides a genomic tool that helps in positional cloning of valuable genes. In this report, an attempt was initiated to study association and localisation of single copy expressed sequence tag (EST) loci in the genome of Corchorus olitorius. The chromosome‐specific association of EST was determined based on the appearance of an extra signal for a single copy cDNA probe in mitotic interphase nuclei of specific trisomic(s) for fluorescence in situ hybridisation, and validated using a cDNA fragment of the 26S rRNA gene (600 bp) as molecular probe. The probe exhibited three signals in meiotic interphase nuclei of trisomic 5, instead of two as observed in diploids and other trisomics, indicating its association with chromosome 5. Subsequent hybridisation of the same probe on the pachytene chromosomes of diploids confirmed that 26S rRNA occupies the terminal end of the short arm of chromosome 5 in C. olitorius. Subsequently, chromosome‐specific association of 63 single copy EST and their physical localisation were determined on chromosomes 2, 4, 5 and 7. The study describes chromosome‐specific physical localisation of genes in jute. The approach used here could be a step towards construction of genome‐wide physical maps for any recalcitrant plant species like jute.     

Analysis of copy number variation using quantitative interspecies competitive PCR

Over recent years small submicroscopic DNA copy-number variants (CNVs) have been highlighted as an important source of variation in the human genome, human phenotypic diversity and disease susceptibility. Consequently, there is a pressing need for the development of methods that allow the efficient, accurate and cheap measurement of genomic copy number polymorphisms in clinical cohorts. We have developed a simple competitive PCR based method to determine DNA copy number which uses the entire genome of a single chimpanzee as a competitor thus eliminating the requirement for competitive sequences to be synthesized for each assay. This results in the requirement for only a single reference sample for all assays and dramatically increases the potential for large numbers of loci to be analysed in multiplex. In this study we establish proof of concept by accurately detecting previously characterized mutations at the PARK2 locus and then demonstrating the potential of quantitative interspecies competitive PCR (qicPCR) to accurately genotype CNVs in association studies by analysing chromosome 22q11 deletions in a sample of previously characterized patients and normal controls.     

Methylation-specific MLPA (MS-MLPA): simultaneous detection of CpG methylation and copy number changes of up to 40 sequences

Copy number changes and CpG methylation of various genes are hallmarks of tumor development but are not yet widely used in diagnostic settings. The recently developed multiplex ligation-dependent probe amplification (MLPA) method has increased the possibilities for multiplex detection of gene copy number aberrations in a routine laboratory. Here we describe a novel robust method: the methylation-specific MLPA (MS-MLPA) that can detect changes in both CpG methylation as well as copy number of up to 40 chromosomal sequences in a simple reaction. In MS-MLPA, the ligation of MLPA probe oligonucleotides is combined with digestion of the genomic DNA–probe hybrid complexes with methylation-sensitive endonucleases. Digestion of the genomic DNA–probe complex, rather than double-stranded genomic DNA, allowed the use of DNA derived from the formalin treated paraffin-embedded tissue samples, enabling retrospective studies. To validate this novel method, we used MS-MLPA to detect aberrant methylation in DNA samples of patients with Prader–Willy syndrome, Angelman syndrome or acute myeloid leukemia.     

Variable copy number of macronuclear DNA molecules in Tetrahymena

In Tetrahymena, the DNA of the macronucleus exists as very large (100 to 4,000-kb) linear molecules that are randomly partitioned to the daughter cells during cell division. This genetic system leads directly to an assortment of alleles such that all loci become homozygous during vegetative growth. Apparently, there is a copy number control mechanism operative that adjusts the number of each macronuclear DNA molecule so that macronuclear DNA molecules (with their loci) are not lost and aneuploid death is a rare event. In comparing Southern analyses of the DNA from various species of Tetrahymena using histone H4 genes as a probe, we find different band intensities in many species. These differences in band intensities primarily reflect differences in the copy number of macronuclear DNA molecules. The variation in copy number of macronuclear DNA molecules in some species is greater than an order of magnitude. These observations are consistent with a developmental control mechanism that operates by increasing the macronuclear copy number of specific DNA molecules (and the genes located on these molecules) to provide the relatively high gene copy number required for highly expressed proteins. © 1992 Wiley-Liss, Inc.     

Methylation-specific multiplex ligation-dependent probe amplification analysis of subjects with chromosome 15 abnormalities

Prader-Willi syndrome (PWS) and Angelman syndrome (AS) are neurodevelopmental disorders caused by loss of expression of imprinted genes from the 15q11-q13 region. They arise from similar defects in the region but differ in parent of origin. There are two recognized typical 15q11-q13 deletions depending on size and several diagnostic assays are available but each has limitations. We evaluated the usefulness of a methylation-specific multiplex ligation-dependent probe amplification (MLPA) kit consisting of 43 probes to detect copy number changes and methylation status in the region. We used the MLPA kit to genotype 82 subjects with chromosome 15 abnormalities (62 PWS, 10 AS and 10 individuals with other chromosome 15 abnormalities) and 13 with normal cytogenetic findings. We developed an algorithm for MLPA probe analysis which correctly identified methylation abnormalities associated with PWS and AS and accurately determined copy number in previously assigned genetic subtypes including microdeletions of the imprinting center. Furthermore, MLPA analysis identified copy number changes in those with distal 15q deletions and ring 15s. MLPA is a relatively simple, cost-effective technique found to be useful and accurate for methylation status, copy number and analysis of genetic subtype in PWS and AS, as well as other chromosome 15 abnormalities.     

Construction and characterization of deletions with defined end points in Drosophila using P elements in trans

We used P-element transposase-mediated "male recombination" between two P elements in trans to create genetic deletions that removed a number of loci, including the gene encoding the neuropeptide crustacean cardioactive peptide (CCAP). Two classes of recombinant chromosomes were produced. Approximately one-quarter were viable when homozygous or hemizygous, whereas the remaining lines caused homozygous and hemizygous lethality. Preliminary analyses using PCR and CCAP immunohistochemistry suggested that, whereas the DNA of the viable lines was largely intact, most lethal lines contained chromosomal deletions that were roughly bounded by the insertion sites of the two P elements used. Southern blot analyses of select lethal lines showed that the DNA flanking the deletion was indeed grossly intact whereas the intervening DNA could not be detected. Sequencing across the deletion in three of these lethal lines identified a single line bearing intact genomic DNA on either side of the deletion separated by 30 bp of P-element DNA. The method described here suggests a simple procedure for creating deletions with defined end points. Importantly, it can use preexisting P-element insertion strains and does not rely on the use of transposable elements that are engineered to cause specific DNA rearrangements.     

Somatic Mutation Allelic Ratio Test Using ddPCR (SMART-ddPCR): An Accurate Method for Assessment of Preferential Allelic Imbalance in Tumor DNA

The extent to which heritable genetic variants can affect tumor development has yet to be fully elucidated. Tumor selection of single nucleotide polymorphism (SNP) risk alleles, a phenomenon called preferential allelic imbalance (PAI), has been demonstrated in some cancer types. We developed a novel application of digital PCR termed Somatic Mutation Allelic Ratio Test using Droplet Digital PCR (SMART-ddPCR) for accurate assessment of tumor PAI, and have applied this method to test the hypothesis that heritable SNPs associated with childhood acute lymphoblastic leukemia (ALL) may demonstrate tumor PAI. These SNPs are located at CDKN2A (rs3731217) and IKZF1 (rs4132601), genes frequently lost in ALL, and at CEBPE (rs2239633), ARID5B (rs7089424), PIP4K2A (rs10764338), and GATA3 (rs3824662), genes located on chromosomes gained in high-hyperdiploid ALL. We established thresholds of AI using constitutional DNA from SNP heterozygotes, and subsequently measured allelic copy number in tumor DNA from 19–142 heterozygote samples per SNP locus. We did not find significant tumor PAI at these loci, though CDKN2A and IKZF1 SNPs showed a trend towards preferential selection of the risk allele (p = 0.17 and p = 0.23, respectively). Using a genomic copy number control ddPCR assay, we investigated somatic copy number alterations (SCNA) underlying AI at CDKN2A and IKZF1, revealing a complex range of alterations including homozygous and hemizygous deletions and copy-neutral loss of heterozygosity, with varying degrees of clonality. Copy number estimates from ddPCR showed high agreement with those from multiplex ligation-dependent probe amplification (MLPA) assays. We demonstrate that SMART-ddPCR is a highly accurate method for investigation of tumor PAI and for assessment of the somatic alterations underlying AI. Furthermore, analysis of publicly available data from The Cancer Genome Atlas identified 16 recurrent SCNA loci that contain heritable cancer risk SNPs associated with a matching tumor type, and which represent candidate PAI regions warranting further investigation.     

Evolution of the ovine MHC DQA region

Southern hybridisation was used to define an apparent gene duplication event at the ovine DQA2 locus. Approximately 500 sheep from five different breeds were genotyped at their DQA1 and DQA2 loci. A subset of these were selected for further characterisation. Southern hybridisation of TaqI digested DNA revealed no DQA1 region in some sheep. It was also noted in these DQA1 null animals the DQA2 specific probe hybridised to two bands. An EcoRV-RFLP designed to distinguish copy number confirmed this duplication of the DQA2 region. The results showed that the duplication was exclusively associated with the DQA1 null haplotype and occurred only in alleles DQA2-F, -G, -I and -J. Comparison with bovine MHC genes revealed that they also contained a DQA1 null haplotype and that this haplotype was associated with a putative DQA3 gene. The potential for an ovine DQA3 locus is discussed.     

Variable copy number of macronuclear DNA molecules in Tetrahymena.

In Tetrahymena, the DNA of the macronucleus exists as very large (100 to 4,000-kb) linear molecules that are randomly partitioned to the daughter cells during cell division. This genetic system leads directly to an assortment of alleles such that all loci become homozygous during vegetative growth. Apparently, there is a copy number control mechanism operative that adjusts the number of each macronuclear DNA molecule so that macronuclear DNA molecules (with their loci) are not lost and aneuploid death is a rare event. In comparing Southern analyses of the DNA from various species of Tetrahymena using histone H4 genes as a probe, we find different band intensities in many species. These differences in band intensities primarily reflect differences in the copy number of macronuclear DNA molecules. The variation in copy number of macronuclear DNA molecules in some species is greater than an order of magnitude. These observations are consistent with a developmental control mechanism that operates by increasing the macronuclear copy number of specific DNA molecules (and the genes located on these molecules) to provide the relatively high gene copy number required for highly expressed proteins.     

SSR allelic variation in almond (Prunus dulcis Mill.)

Sixteen SSR markers including eight EST-SSR and eight genomic SSRs were used for genetic diversity analysis of 23 Chinese and 15 international almond cultivars. EST- and genomic SSR markers previously reported in species of Prunus, mainly peach, proved to be useful for almond genetic analysis. DNA sequences of 117 alleles of six of the 16 SSR loci were analysed to reveal sequence variation among the 38 almond accessions. For the four SSR loci with AG/CT repeats, no insertions or deletions were observed in the flanking regions of the 98 alleles sequenced. Allelic size variation of these loci resulted exclusively from differences in the structures of repeat motifs, which involved interruptions or occurrences of new motif repeats in addition to varying number of AG/CT repeats. Some alleles had a high number of uninterrupted repeat motifs, indicating that SSR mutational patterns differ among alleles at a given SSR locus within the almond species. Allelic homoplasy was observed in the SSR loci because of base substitutions, interruptions or compound repeat motifs. Substitutions in the repeat regions were found at two SSR loci, suggesting that point mutations operate on SSRs and hinder the further SSR expansion by introducing repeat interruptions to stabilize SSR loci. Furthermore, it was shown that some potential point mutations in the flanking regions are linked with new SSR repeat motif variation in almond and peach. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.