Patchwork: allele-specific copy number analysis of whole-genome sequenced tumor tissue

Whole-genome sequencing of tumor tissue has the potential to provide comprehensive characterization of genomic alterations in tumor samples. We present Patchwork, a new bioinformatic tool for allele-specific copy number analysis using whole-genome sequencing data. Patchwork can be used to determine the copy number of homologous sequences throughout the genome, even in aneuploid samples with moderate sequence coverage and tumor cell content. No prior knowledge of average ploidy or tumor cell content is required. Patchwork is freely available as an R package, installable via R-Forge (http://patchwork.r-forge.r-project.org/).     

Rapid PCR-based determination of transgene copy number in rice

We present a simple, rapid, and low-cost method to determine transgene copy number in rice. More than 100 first- and second-generation transgenic rice plants were tested. The plasmid (pRCopy) used for rice transformation contains the specific gene of interest and a partially deleted cytochrome c gene (cyc), a single-copy gene in rice. A 132-bp segment of the cloned ricecyc was shortened to 108 bp by deleting a 24-bp internal fragment. After PCR amplification of the genomic DNA from transgenic rice harboring pRCopy, the 2 expected bands were found. The 121-bp band corresponds to the endogenouscyc; the 97-bp band comes from the integrated pRCopy. Clear distinctions can be made between single and multiple copies of the transgene by comparing band densities.     

Methodological strategies for transgene copy number quantification in goats (Capra hircus) using real‐time PCR

Taking into account the importance of goats as transgenic models, as well as the rarity of copy number (CN) studies in farm animals, the present work aimed to evaluate methodological strategies for accurate and precise transgene CN quantification in goats using quantitative polymerase chain reaction (qPCR). Mouse and goat lines transgenic for human granulocyte‐colony stimulating factor were used. After selecting the best genomic DNA extraction method to be applied in mouse and goat samples, intra‐assay variations, accuracy and precision of CN quantifications were assessed. The optimized conditions were submitted to mathematical strategies and used to quantify CN in goat lines. The findings were as follows: validation of qPCR conditions is required, and amplification efficiency is the most important. Absolute and relative quantifications are able to produce similar results. For normalized absolute quantification, the same plasmid fragment used to generate goat lines must be mixed with wild‐type goat genomic DNA, allowing the choice of an endogenous reference gene for data normalization. For relative quantifications, a resin‐based genomic DNA extraction method is strongly recommended when using mouse tail tips as calibrators to avoid tissue‐specific inhibitors. Efficient qPCR amplifications (≥95%) allow reliable CN measurements with SYBR technology. TaqMan must be used with caution in goats if the nucleotide sequence of the endogenous reference gene is not yet well understood. Adhering to these general guidelines can result in more exact CN determination in goats. Even when working under nonoptimal circumstances, if assays are performed that respect the minimum qPCR requirements, good estimations of transgene CN can be achieved. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 30:1390–1400, 2014     

质粒拷贝数的测定方法

在研究生物工程重组蛋白产量的过程中 ,质粒的拷贝数是一个需要重点考虑的参数 ,对它进行精确定量至关重要。本文概述了几种主要的测定方法的原理和特点。     

High-throughput transgene copy number estimation by competitive PCR

Transgene copy number affects the level and stability of gene expression. Therefore, it is important to determine the copy number of each transgenic line. Polymerase chain reaction (PCR) is widely employed to quantify amounts of target sequences. Although PCR is not inherently quantitative, various means of overcoming this limitation have been devised. Recent real-time PCR methods are rapid; however, they typically lack a suitable internal standard, limit the size of the target sequence, and require expensive specialized equipment. Competitive PCR techniques avoid these problems, but traditional competitive methods are time consuming. Here we apply mathematical modeling to create a rapid, simple, and inexpensive copy number determination method that retains the robustness of competitive PCR.     

实时荧光定量PCR(TaqMan)法测定外源基因的拷贝数

实时荧光定量PCR是近年新兴的一项技术,因其快速、方便、便宜,需要DNA样品量少,无需放射性检测等优点被广泛应用于基因的定量分析。该文就实时荧光定量PCR(TaqMan)技术的发展、基本原理及测定外源基因拷贝数的技术流程做一介绍。     

Recurrent patterns of DNA copy number alterations in tumors reflect metabolic selection pressures

Copy number alteration (CNA) profiling of human tumors has revealed recurrent patterns of DNA amplifications and deletions across diverse cancer types. These patterns are suggestive of conserved selection pressures during tumor evolution but cannot be fully explained by known oncogenes and tumor suppressor genes. Using a pan‐cancer analysis of CNA data from patient tumors and experimental systems, here we show that principal component analysis‐defined CNA signatures are predictive of glycolytic phenotypes, including ¹⁸F‐fluorodeoxy‐glucose (FDG) avidity of patient tumors, and increased proliferation. The primary CNA signature is enriched for p53 mutations and is associated with glycolysis through coordinate amplification of glycolytic genes and other cancer‐linked metabolic enzymes. A pan‐cancer and cross‐species comparison of CNAs highlighted 26 consistently altered DNA regions, containing 11 enzymes in the glycolysis pathway in addition to known cancer‐driving genes. Furthermore, exogenous expression of hexokinase and enolase enzymes in an experimental immortalization system altered the subsequent copy number status of the corresponding endogenous loci, supporting the hypothesis that these metabolic genes act as drivers within the conserved CNA amplification regions. Taken together, these results demonstrate that metabolic stress acts as a selective pressure underlying the recurrent CNAs observed in human tumors, and further cast genomic instability as an enabling event in tumorigenesis and metabolic evolution.     

Discovering tumor suppressor genes through genome-wide copy number analysis

Classical tumor suppressor gene discovery has largely involved linkage analysis and loss-of-heterozygosity (LOH) screens, followed by detailed mapping of relatively large chromosomal regions. Subsequent efforts made use of genome-wide PCR-based methods to detect rare homozygous deletions. More recently, high-resolution genomic arrays have been applied to cancer gene discovery. However, accurate characterization of regions of genomic loss is particularly challenging due to sample heterogeneity, the small size of deleted regions and the high frequency of germline copy number polymorphisms. Here, we review the application of genome-wide copy number analysis to the specific problem of identifying tumor suppressor genes.     

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

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

Detection of differential gene copy number using denaturing high performance liquid chromatography

Genomic rearrangements leading to deletion or duplication of gene(s) resulting in alterations in gene copy number underlie the molecular lesion in several genetic disorders. Methods currently used to determine gene copy number including real time PCR, southern hybridization, fluorescence in situ hybridization, densitometric scanning of PCR product etc. have certain disadvantages and are also expensive and time consuming. Herein, we describe a simple and rapid method to assess gene copy number using denaturing high performance liquid chromatography (dHPLC). We used X chromosome genes as model to compare the gene copy numbers present on this chromosome in males and females. DNA from these samples were amplified by biplex PCR using primer pairs specific for X chromosome genes only (target gene) and for genes present on both X and Y chromosomes (internal control). Amplified products were analyzed using HPLC under non-denaturing conditions. The ratio of peak areas (target gene/internal control) of the amplified products was approximately twice in female samples than male samples (p < 0.001) demonstrating that the differential gene copy number can be easily detected using this method. This method can potentially be used for diagnostic purpose where the need is to distinguish samples based on the differential gene copy numbers.     

Detection of Stachybotrys chartarum using rRNA, tri5, and β-tubulin primers and determining their relative copy number by real-time PCR

Highly conserved regions are attractive targets for detection and quantitation by PCR, but designing species-specific primer sets can be difficult. Ultimately, almost all primer sets are designed based upon literature searches in public domain databases, such as the National Center for Biotechnology Information (NCBI). Prudence suggests that the researcher needs to evaluate as many sequences as available for designing species-specific PCR primers. In this report, we aligned 11, 9, and 16 DNA sequences entered for Stachybotrys spp. rRNA, tri5, and β-tubulin regions, respectively. Although we were able to align and determine consensus primer sets for the 9 tri5 and the 16 β-tubulin sequences, there was no consensus sequence that could be derived from alignment of the 11 rRNA sequences. However, by judicious clustering of the sequences that aligned well, we were able to design three sets of primers for the rRNA region of S. chartarum. The two primer sets for tri5 and β-tubulin produced satisfactory PCR results for all four strains of S. chartarum used in this study whereas only one rRNA primer set of three produced similar satisfactory results. Ultimately, we were able to show that rRNA copy number is approximately 2-log greater than for tri5 and β-tubulin in the four strains of S. chartarum tested.     

Using real-time PCR to determine transgene copy number in wheat

Transgene copy number is usually determined by means of Southern blot analysis which can be time consuming and laborious. In this study, quantitative real-time PCR was developed to determine transgene copy number in transgenic wheat. A conserved wheat housekeeping gene,puroindoline-b, was used as an internal control to calculate transgene copy number. Estimated copy number in transgenic lines using real-time quantitative PCR was correlated with actual copy number based on Southern blot analysis. Real-time PCR can analyze hundreds of samples in a day, making it an efficient method for estimating copy number in transgenic wheat.     

Plasmid copy number and stability determination in Clostridium perfringens transformants

The copy number of the streptococcal plasmid pAM beta 1 (26.5 kb), and its deletion derivatives, pVA1 (11 kb) and pVA677 (7.6 kb) contained in Clostridium perfringens 3624A transformants was determined by incorporation of [methyl-3H]thymidine (4 muCi/ml) into chromosomal and plasmid DNA and sizing of the C. perfringens genome using transverse alternating field electrophoresis. Plasmids pAM beta 1, pVA1, and pVA677 were found to be present at 1.0, 97, and 216 copies/cell, respectively. 10.2, 54, and 96% of the initial pAM beta 1-, pVA1- and pVA677-containing transformants, respectively, remained resistant to erythromycin over 220 generations of growth. The results indicate a size-dependent relationship between plasmid stability and plasmid copy number in C. perfringens.     

High-throughput detection of multiple genetic polymorphisms influencing drug metabolism with mismatch primers in allele-specific polymerase chain reaction

Because of genetic polymorphisms of drug-metabolizing enzyme genes, the activities of the enzymes in humans vary widely and alter the metabolism of commonly used clinical agents. Severe adverse effects or resistance to therapy may result. We have developed a rapid and high-throughput genotyping method for detecting polymorphisms of the drug-metabolizing enzyme genes CYP2C9*3, CYP2C19*2, *3, CYP2D6*2, *4, *10, *14, *21, NAT2*5, *6, *7, and TPMT*3 using allele-specific polymerase chain reaction (PCR) with mismatch primers (ASPCR-MP) and CYP2D6*5, *36, and CYP2D6xN using stepdown PCR with detection by SYBR Green I. We analyzed genomic DNA from 139 Japanese volunteers. Identical genotyping results were obtained by using ASPCR-MP, stepdown PCR, and conventional PCR. We found that the methods clearly differentiate three specific profiles with no overlap in the signals. Moreover, both ASPCR-MP and stepdown PCR for genotyping took less than 3-4h. To our knowledge, this is the first report of successful simultaneous detection of multiple genetic polymorphisms with point mutations using ASPCR-MP or multiple genetic polymorphisms with large structural alterations using stepdown PCR. In conclusion, ASPCR-MP and stepdown PCR appear to be suitable for large clinical and epidemiological studies as methods that enable highly sensitive genotyping and yield a high-throughput.     

Number matters: control of mammalian mitochondrial DNA copy number

Regulation of mitochondrial biogenesis is essential for proper cellular functioning. Mitochondrial DNA （mtDNA） depletion and the resulting mitochondrial malfunction have been implicated in cancer, neurodegeneration, diabetes, aging, and many other human diseases. Although it is known that the dynamics of the mammalian mitochondrial genome are not linked with that of the nuclear genome, very little is known about the mechanism of mtDNA propagation. Nevertheless, our understanding of the mode of mtDNA replication has ad- vanced in recent years, though not without some controversies. This review summarizes our current knowledge of mtDNA copy number control in mammalian cells, while focusing on both mtDNA replication and turnover. Although mtDNA copy number is seemingly in excess, we reason that mtDNA copy number control is an important aspect of mitochondrial genetics and biogenesis and is essential for normal cellular function.     

Rapid determination of transgene copy number in stably-transfected mammalian cells by competitive PCR

We describe here an application of the competitive PCR technique to the analysis of copy number of recombinant rat parathyroid hormone-related protein (rPTHrP) gene in stably-transfected murine erythroleukemia (MEL) cell lines. A single-copy reference gene (endogenous mouse PTHrP gene or mPTHrP) is used as an internal control. This control gene, present in the genome of MEL cells, shares the same primer binding sites as the rPTHrP cDNA but contains an internal PvuII site, which allows resolution of the amplified products after restriction enzyme digestion by polyacrylamide gel electrophoresis (PAGE). The transgene copy number is determined by the ratio of band intensity of the rPTHrP product to that of the mPTHrP product. Using this method, we have determined the copy number of the rPTHrP transgene from isolated genomic DNA, and compared the results with those obtained from Southern blot analysis. In addition, we have demonstrated that the procedure can be applied very simply to whole MEL cells without DNA extractions and that as few as 10⁴ cells are required for the analysis.     

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30–35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.     

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.