Gene copy number variations in adaptive evolution: The genomic distribution of gene copy number variations revealed by genetic mapping and their adaptive role in an undomesticated species,white spruce (Picea glauca)

Gene copy number variation (CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays (aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution.     

Gene copy number variations as signatures of adaptive evolution in the parthenogenetic,plant‐parasitic nematode Meloidogyne incognita

Adaptation to changing environmental conditions represents a challenge to parthenogenetic organisms, and until now, how phenotypic variants are generated in clones in response to the selection pressure of their environment remains poorly known. The obligatory parthenogenetic root‐knot nematode species Meloidogyne incognita has a worldwide distribution and is the most devastating plant‐parasitic nematode. Despite its asexual reproduction, this species exhibits an unexpected capacity of adaptation to environmental constraints, for example, resistant hosts. Here, we used a genomewide comparative hybridization strategy to evaluate variations in gene copy numbers between genotypes of M. incognita resulting from two parallel experimental evolution assays on a susceptible vs. resistant host plant. We detected gene copy number variations (CNVs) associated with the ability of the nematodes to overcome resistance of the host plant, and this genetic variation may reflect an adaptive response to host resistance in this parthenogenetic species. The CNV distribution throughout the nematode genome is not random and suggests the occurrence of genomic regions more prone to undergo duplications and losses in response to the selection pressure of the host resistance. Furthermore, our analysis revealed an outstanding level of gene loss events in nematode genotypes that have overcome the resistance. Overall, our results support the view that gene loss could be a common class of adaptive genetic mechanism in response to a challenging new biotic environment in clonal animals.     

Parallel evolution? Microsatellite variation of recently isolated marine and freshwater three-spined stickleback

Variation at nine microsatellite loci, four of which are linked to phenotypic traits (spine length and lateral plate morphology) in Canadian three-spined stickleback Gasterosteus aculeatus , are used to test for selection on marine and freshwater three-spined stickleback morphs in Iceland. There are indications of strong selection on loci linked to dorsal spine length, providing another potential example of parallel divergence at the genomic level in three-spined stickleback.     

Temporal variation in divergent selection on spine number in threespine stickleback

Short-term temporal cycles in ecological pressures, such as shifts in predation regime, are widespread in nature yet estimates of temporal variation in the direction and intensity of natural selection are few. Previous work on threespine stickleback (Gasterosteus aculeatus) has revealed that dorsal and pelvic spines are a defense against gape-limited predators but may be detrimental against grappling insect predators. In this study, we examined a 15-year database from an endemic population of threespine stickleback to look for evidence of temporal shifts in exposure to these divergent predation regimes and correlated shifts in selection on spine number. For juveniles, we detected selection for increased spine number during winter when gape-limited avian piscivores were most common but selection for decreased spine number during summer when odonate predation was more common. For subadults and adults, which are taken primarily by avian piscivores, we predicted selection should generally be for increased spine number in all seasons. Among 59 comparisons, four selection differentials were significant (Bonferroni corrected) and in the predicted direction. However, there was also substantial variability in remaining differentials, including two examples with strong selection for spine reduction. These reversals were associated with increased tendency of the fish to shift to a benthic niche, as determined from examination of stomach contents. These dietary data suggest that increased encounter rates with odonate predation select for spine reduction. Strong selection on spine number was followed by changes in mean spine number during subsequent years and a standard quantitative genetic formula revealed that spine number has a heritable component. Our results provide evidence of rapid morphological responses to selection from predators and suggest that temporal variation in selection may help maintain variation within populations. Furthermore, our findings indicate that variable selection can be predicted if the agents of selection are known.     

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

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

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.     

Copy number variations of HLA-DRB5 is associated with systemic lupus erythematosus risk in Chinese Han population

Systemic lupus erythematosus （SLE） is a polygenic, systemic, autoimmune disease. Copy number variants （CNVS） have been discovered to be associated with a number of complex disorders. We undertook the current study to explore the potential associations between genomic CNVS and SLE in Chinese Han population. In the discovery stage, seven SLE patients were examined with the high-density comparative genomic hybridization microarrays in the screening test for SLE associated CNVS. Then, in the validation stage, 135 SLE patients and 219 matched healthy subjects were investigated for the CNVS of gene HLA-DRB5 by AccuCopyTM technol- ogy. Quantitative polymerase chain reaction was carried out to determine the copy number （CN） and mRNA level of HLA- DRB5 in SLE patients. Although the mRNA level of HLA- DRB5 between the CN deletion group and the CN normal group in SLE patients was not statistically positive （P = 0.46）, our results still showed more CN of HLA-DRB5 in SLE patients than in healthy controls （P = 3.98×10^-6）. Odds ratio for CN deletion was 0.38 （95% confidence interval （C1）, 0.23-0.61, P = 7.79×10^-5） and for CN duplication was 1.89 （95% CI, 0.56-7.66, P = 0.37）, respectively. These findings indicated that CNVS of HLA-DRB5 was associated with the risk of SLE, and CN deletion appeared to be protective for SLE.     

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.     

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.     

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.     

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.     

Segmental duplications are hot spots of copy number variants affecting barley gene content

Copy number variants (CNVs) are pervasive in several animal and plant genomes and contribute to shaping genetic diversity. In barley, there is evidence that changes in gene copy number underlie important agronomic traits. The recently released reference sequence of barley represents a valuable genomic resource for unveiling the incidence of CNVs that affect gene content and for identifying sequence features associated with CNV formation. Using exome sequencing and read count data, we detected 16 605 deletions and duplications that affect barley gene content by surveying a diverse panel of 172 cultivars, 171 landraces, 22 wild relatives and other 32 uncategorized domesticated accessions. The quest for segmental duplications (SDs) in the reference sequence revealed many low‐copy repeats, most of which overlap predicted coding sequences. Statistical analyses revealed that the incidence of CNVs increases significantly in SD‐rich regions, indicating that these sequence elements act as hot spots for the formation of CNVs. The present study delivers a comprehensive genome‐wide study of CNVs affecting barley gene content and implicates SDs in the molecular mechanisms that lead to the formation of this class of CNVs.