Delineation of the frequency and boundary of chromosomal copy number variations in paediatric neuroblastoma

Neuroblastoma, the most common solid tumour in early childhood, is characterized by very frequent chromosomal copy number variations (CNVs). While chromosome 2p amplification, 17q gain, 1p and 11q deletion in human neuroblastoma tissues are well-known, the exact frequencies and boundaries of the chromosomal CNVs have not been delineated. We analysed the publicly available single nucleotide polymorphism (SNP) array data which were originally generated by the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative, defined the frequencies and boundaries of chromosomes 2p11.2 – 2p25.3 amplification, 17q11.1-17q25.3 gain, 1p13.3-1p36.33 deletion and 11q13.3-11q25 deletion in neuroblastoma tissues, and identified chromosome 7q14.1 (Chr7:38254795-38346971) and chromosome 14q11.2 (Chr14:21637401-22024617) deletion in blood and bone marrow samples from neuroblastoma patients, but not in tumour tissues. Kaplan Meier analysis showed that double deletion of Chr7q14.1 and Chr14q11.2 correlated with poor prognosis in MYCN gene amplified neuroblastoma patients. In conclusion, the oncogenes amplified or gained and tumour suppressor genes deleted within the boundaries of chromosomal CNVs in tumour tissues should be studied for their roles in tumourigenesis and as therapeutic targets. Focal deletions of Chr7q14.1 and Chr14q11.2 together in blood and bone marrow samples from neuroblastoma patients can be used as a marker for poorer prognosis and more aggressive therapies.     

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.     

Optimization protein productivity of human interleukin-2 through codon usage,gene copy number and intracellular tRNA concentration in CHO cells

Transfer RNA (tRNA) abundance is one of the critical factors for the enhancement of protein productivity in prokaryotic and eukaryotic hosts. Gene copy number of tRNA and tRNA codon usage bias are generally used to match tRNA abundance of protein-expressing hosts and to optimize the codons of recombinant proteins. Because sufficient concentration of intracellular tRNA and optimized codons of recombinant proteins enhanced translation efficiency, we hypothesized that sufficient supplement of host’s tRNA improved protein productivity in mammalian cells. First, the small tRNA sequencing results of CHO-K1 cells showed moderate positive correlation with gene copy number and codon usage bias. Modification of human interleukin-2 (IL-2) through codons with high gene copy number and high codon usage bias (IL-2 HH, modified on Leu, Thr, Glu) significantly increased protein productivity in CHO-K1 cells. In contrast, modification through codons with relatively high gene copy number and low codon usage bias (IL-2 HL, modified on Ala, Thr, Val), or relatively low gene copy number and low codon usage bias (IL-2 LH, modified on Ala, Thr, Val) did not increase IL-2 productivity significantly. Furthermore, supplement of the alanine tRNA or threonine tRNA increased IL-2 productivity of IL-2 HL. In summary, we revealed a potential strategy to enhance productivity of recombinant proteins, which may be applied in production of protein drug or design of DNA vaccine.     

质粒拷贝数的测定方法

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

Determination of 16S ribosomal RNA gene copy number in Streptococcus uberis, S. agalactiae, S. dysgalactiae and S. parauberis

Abstract Species-specific oligonucleotide probes and a universal oligonucleotide probe derived from sequences of 16S rRNA were hybridised to chromosomal DNA from Streptococcus agalactiae, S. dysgalactiae, S. parauberis and S. uberis following digestion with Eco RI. Due to the presence of a unique Eco RI site in each 16S rRNA gene, the number of hybridised fragments was indicative of the number of 16S rRNA genes. Southern hybridisation indicated six 16S rRNA genes in ten isolates of S. agalactiae , five genes in ten isolates of S. uberis , five genes in six isolates and six in another isolate of S. dysgalactiae , and six genes in four isolates of S. parauberis . For a fifth isolate of S. parauberis , six 16S rRNA genes were indicated by the universal probe but only five when hybridised to the species-specific probe, indicating sequence variation (microheterogeneity) within the probe target region.     

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.     

Detection of copy number variations and their effects in Chinese bulls

Background

Copy number variations (CNVs) are a main source of genomic structural variations underlying animal evolution and production traits. Here, with one pure-blooded Angus bull as reference, we describe a genome-wide analysis of CNVs based on comparative genomic hybridization arrays in 29 Chinese domesticated bulls and examined their effects on gene expression and cattle growth traits.

Results

We identified 486 copy number variable regions (CNVRs), covering 2.45% of the bovine genome, in 24 taurine (Bos taurus), together with 161 ones in 2 yaks (Bos grunniens) and 163 ones in 3 buffaloes (Bubalus bubalis). Totally, we discovered 605 integrated CNVRs, with more “loss” events than both “gain” and “both” ones, and clearly clustered them into three cattle groups. Interestingly, we confirmed their uneven distributions across chromosomes, and the differences of mitochondrion DNA copy number (gain: taurine, loss: yak & buffalo). Furthermore, we confirmed approximately 41.8% (253/605) and 70.6% (427/605) CNVRs span cattle genes and quantitative trait loci (QTLs), respectively. Finally, we confirmed 6 CNVRs in 9 chosen ones by using quantitative PCR, and further demonstrated that CNVR22 had significantly negative effects on expression of PLA2G2D gene, and both CNVR22 and CNVR310 were associated with body measurements in Chinese cattle, suggesting their key effects on gene expression and cattle traits.

Conclusions

The results advanced our understanding of CNV as an important genomic structural variation in taurine, yak and buffalo. This study provides a highly valuable resource for Chinese cattle’s evolution and breeding researches.

Electronic supplementary material

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

Enhancement of bovine growth hormone gene expression by increasing the plasmid copy number

Effect of copy number on the expression of bovine growth hormone gene (bGH) was investigated using the copy number mutants such as pKBJ10, pBJ( tet)10, pUBJ10-1, and pUBJ10 plasmids. The cells harboring plasmids below 84 copies/cell did not produced detectable levels of bGH. When the ColE1 replicon was replaced with the mutated ColE1 replicon originated from pUC19 plasmid, the copy number was increased to about 300 copies/cell and bGH production was enhanced by 11.5% (pUBJ10-1) and 12.3% of total cell protein (pUBJ10). A large amount of mRNA caused by increment of copy number would be needed to overcome some inhibitory threshold and might be an important factor for regulating bGH expression.     

tRNA gene copy number variation in humans

The human tRNAome consists of more than 500 interspersed tRNA genes comprising 51 anticodon families of largely unequal copy number. We examined tRNA gene copy number variation (tgCNV) in six individuals; two kindreds of two parents and a child, using high coverage whole genome sequence data. Such differences may be important because translation of some mRNAs is sensitive to the relative amounts of tRNAs and because tRNA competition determines translational efficiency vs. fidelity and production of native vs. misfolded proteins. We identified several tRNA gene clusters with CNV, which in some cases were part of larger iterations. In addition there was an isolated tRNALysCUU gene that was absent as a homozygous deletion in one of the parents. When assessed by semiquantitative PCR in 98 DNA samples representing a wide variety of ethnicities, this allele was found deleted in hetero- or homozygosity in all groups at ~ 50% frequency. This is the first report of copy number variation of human tRNA genes. We conclude that tgCNV exists at significant levels among individual humans and discuss the results in terms of genetic diversity and prior genome wide association studies (GWAS) that suggest the importance of the ratio of tRNALys isoacceptors in Type-2 diabetes.     

Germline DNA copy number variation in individuals with Argyrophilic grain disease reveals CTNS as a plausible candidate gene

Argyrophilic grain disease (AGD) is a progressive neurodegenerative disease of the human brain that has never been associated to a particular gene locus. In the present study, we report the results of a CNV investigation in 29 individuals whose anatomopathologic investigation of the brain showed AGD. Rare CNVs were identified in six patients (21%), in particular a 40 kb deletion at 17p13.2 encompassing the CTNS gene. Homozygote mutations in CTNS are known to cause cystinosis, a disorder characterized by the intralysosomal accumulation of cystine in all tissues. We present the first CNV results in individuals presenting AGD and a possible candidate gene implicated in the disorder.     

Transgene copy number comparison in recombinant mammalian cell lines: critical reflection of quantitative real-time PCR evaluation

Nucleic acid quantification is a relevant issue for the characterization of mammalian recombinant cell lines and also for the registration of producer clones. Quantitative real-time PCR is a powerful tool to investigate nucleic acid levels but numerous different quantification strategies exist, which sometimes lead to misinterpretation of obtained qPCR data. In contrast to absolute quantification using amplicon- or plasmid standard curves, relative quantification strategies relate the gene of interest to an endogenous reference gene. The relative quantification methods also consider the amplification efficiency for the calculation of the gene copy number and thus more accurate results compared to absolute quantification methods are generated. In this study two recombinant Chinese hamster ovary cell lines were analysed for their transgene copy number using different relative quantification strategies. The individual calculation methods resulted in differences of relative gene copy numbers because efficiency calculations have strong impact on gene copy numbers. However, in context of comparing transgene copy numbers of two individual clones the influence of the calculation method is marginal. Therefore especially for the comparison of two cell lines with the identical transgene any of the relative qPCR methods was proven as powerful tool.     

Theoretical framework of population genetics with somatic mutations taken into account: application to copy number variations in humans

Traditionally, population genetics focuses on the dynamics of frequencies of alleles acquired by mutations on germ-lines, because only such mutations are heritable. Typical genotyping experiments, however, use DNA from some somatic tissues such as blood, which harbors somatic mutations at the current generation in addition to germ-line mutations accumulated since the most recent common ancestor of the sample. This common practice may sometimes cause erroneous interpretations of polymorphism data, unless we properly understand the role of somatic mutations in population genetics. We here introduce a very basic theoretical framework of population genetics with somatic mutations taken into account. It is easy to imagine that somatic mutations at the current generation simply add individual-specific variations, as errors in mutation detection do. Our theory quantifies this increment under various conditions. We find that the major contribution of somatic mutations plus errors is to very rare variants, particularly to singletons. The relative contribution is markedly large when mutations are deleterious. Because negative selection also increases rare variants, it is important to distinguish the roles of these mutually confounding factors when we interpret the data, even after correcting for demography. We apply this theory to human copy number variations (CNVs), for which the composite effect of somatic mutations and errors may not be negligible. Using genome-wide CNV data, we demonstrate how the joint action of the two factors, selection and somatic mutations plus errors, shapes the observed pattern of polymorphism.     

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.     

Integrated analyses of copy number variations and gene differential expression in lung squamous-cell carcinoma

Background

Although numerous efforts have been made, the pathogenesis underlying lung squamous-cell carcinoma (SCC) remains unclear. This study aimed to identify the CNV-driven genes by an integrated analysis of both the gene differential expression and copy number variation (CNV).

Results

A higher burden of the CNVs was found in 10–50 kb length. The 16 CNV-driven genes mainly located in chr 1 and chr 3 were enriched in immune response [e.g. complement factor H (CFH) and Fc fragment of IgG, low affinity IIIa, receptor (FCGR3A)], starch and sucrose metabolism [e.g. amylase alpha 2A (AMY2A)]. Furthermore, 38 TFs were screened for the 9 CNV-driven genes and then the regulatory network was constructed, in which the GATA-binding factor 1, 2, and 3 (GATA1, GATA2, GATA3) jointly regulated the expression of TP63.

Conclusions

The above CNV-driven genes might be potential contributors to the development of lung SCC.     

Modelling the evolution of genomes with integrated external and internal functions

The genomes that organisms transmit between generations contain information about different kinds of functions. The genome with the "best" mix and number of genes for these functions is the one that natural selection favours. Here I introduce a new way to model simple organisms with genes for external and internal functions, and use it to study the evolution of genome size. The external functions are exemplified by resource use and the internal functions by mutation control (repair). It is shown that even with a suitable proportion of genes for mutation control, the genomes in the organisms do not forever incorporate genes that increase resource use. Instead they evolve towards an optimal genome of limited size. The optimal proportion of genes for mutation control is shown to have an upper limit given by the ease with which transmission accuracy is improved by adding extra genes for this purpose to the genome. The model illustrates how natural selection on genomes integrates systems for the transmission of genetic information with systems relating to the external adaptation of the organism. It also opens up for other, more detailed theoretical investigations of genome functions.     

An evaluation of models for the effect of plasmid copy number on bacterial growth rate

The rates of growth of recombinant bacteria depend on their plasmid content. This is modelled by expressing the specific growth rate in terms of the number of copies of the plasmid per cell. Three models in common use have been tested with different Escherichia coli strains and one strain of Bacillus stearothermophilus containing different plasmids. While no particular model was decisively better than others for all data, that of Bentley & Quiroga (Biotechnol. Bioeng. 1993, 42: 222–234) was the best for specific growth rates which vary inversely with the plasmid copy number, and a modified form of the model of Satyagal & Agarwal (Biotechnol. Bioeng. 1989, 33: 1135–1144) was the best for growth rates which increase with the copy number. The differences appear to be linked to the plasmid replication mechanisms. Contrary to some claims, no model portrayed the experimentally observed inflection points.     

Parallel evolution of IDH2 gene in cetaceans,primates and bats

Cetaceans and primates both have large brains that require large amounts of aerobic energy metabolism. In bats, the cost of flight makes locomotion energetically demanding. These mammalian groups may represent three independent evolutionary origins of an energy-demanding lifestyle in mammals. IDH2 encodes an enzyme in the tricarboxylic acid cycle in the mitochondrion, which plays a key role in aerobic energy metabolism. In this study, we cloned and sequenced this gene in two cetaceans, and 19 bat species, and compared the data with available primate sequences to test its evolution. We found significant signals of parallel evolution in this gene among these three groups. Parallel evolution of this gene may reflect their parallel evolution towards a higher demand for energy.     

Genome-wide integrative analysis revealed a correlation between lengths of copy number segments and corresponding gene expression profile

Microarray analysis has been applied to comprehensively reveal the abnormalities of DNA copy number (CN) and gene expression in human cancer research during the last decade. These analyses have individually contributed to identify the genes associated with carcinogenesis, progression, metastasis of tumor cells and poor prognosis of cancer patients. However, it is known that the correlation between profiles of CN and gene expression does not highly correlate. Factors which determine the degree of correlation remain largely unexplained. To investigate one such factor, we performed trend analyses between the lengths of CN segments and corresponding gene expression profiles from microarray data in hepatocellular carcinoma (HCC) and colorectal carcinoma (CRC). Significant correlations were observed in CN gain of HCC and CRC (p<0.05). The trend of the CN loss showed a significant correlation in HCC although there was no correlation between the length of CN loss segments and gene expression in CRC. Our findings suggest that the influence of CN on gene expression highly depends on the length of CN region, especially in the case of CN gain. To the best of our knowledge, this is the first study describing the correlation between lengths of CNA segments and expression profiles of corresponding genes.