A genome wide association study between copy number variation （CNV） and human height in Chinese population

Copy number variation (CNV) is a type of genetic variation which may have important roles in phenotypic variability and disease susceptibility. To hunt for genetic variants underlying human height variation, we performed a genome wide CNV association study for human height in 618 Chinese unrelated subjects using Affymetrix 500K array set. After adjusting for age and sex, we found that four CNVs at 6p21.3, 8p23.3-23.2, 9p23 and 16p12.1 were associated with human height (with borderline significant p value: 0.013, 0.011, 0.024, 0.049; respectively). However, after multiple tests correction, none of them was associated with human height. We observed that the gain of copy number (more than 2 copies) at 8p23.3-23.2 was associated with lower height (normal copy number vs. gain of copy number; 161.2 cm vs. 153.7 cm, p = 0.011), which accounted for 0.9% of height variation. Loss of copy number (less than 2 copies) at 6p21.3 was associated with 0.8% lower height (loss of copy number vs. normal copy number: 154.5 cm vs. 161.1 cm, p = 0.013). Since no important genes influencing height located in CNVs at loci of 8p23.3-23.2 and 6p21.3, the two CNVs may cause the structural rearrangements of neighbored important candidate genes, thus regulates the variation of height. Our results expand our knowledge of the genetic factors underlying height variation and the biological regulation of human height.     

Diversity of copy number variation in a worldwide population of sheep

Copy number variation (CNV) represents a major source of genomic variation. We investigated the diversity of CNV distribution using SNP array data collected from a comprehensive collection of geographically dispersed sheep breeds. We identified 24,558 putative CNVs, which can be merged into 619 CNV regions, spanning 197 Mb of total length and corresponding to ~ 6.9% of the sheep genome. Our results reveal a population differentiation in CNV between different geographical areas, including Africa, America, Asia, Southwestern Asia, Central Europe, Northern Europe and Southwestern Europe. We observed clear distinctions in CNV prevalence between diverse groups, possibly reflecting the population history of different sheep breeds. We sought to determine the gene content of CNV, and found several important CNV-overlapping genes (BTG3, PTGS1 and PSPH) which were involved in fetal muscle development, prostaglandin (PG) synthesis, and bone color. Our study generates a comprehensive CNV map, which may contribute to genome annotation in sheep.     

Germ-line DNA copy number variation frequencies in a large North American population

Genomic copy number variation (CNV) is a recently identified form of global genetic variation in the human genome. The Affymetrix GeneChip 100 and 500 K SNP genotyping platforms were used to perform a large-scale population-based study of CNV frequency. We constructed a genomic map of 578 CNV regions, covering approximately 220 Mb (7.3%) of the human genome, identifying 183 previously unknown intervals. Copy number changes were observed to occur infrequently (<1%) in the majority (>93%) of these genomic regions, but encompass hundreds of genes and disease loci. This North American population-based map will be a useful resource for future genetic studies. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Deletion 17q12 is a recurrent copy number variant that confers high risk of autism and schizophrenia

Autism spectrum disorders (ASD) and schizophrenia are neurodevelopmental disorders for which recent evidence indicates an important etiologic role for rare copy number variants (CNVs) and suggests common genetic mechanisms. We performed cytogenomic array analysis in a discovery sample of patients with neurodevelopmental disorders referred for clinical testing. We detected a recurrent 1.4 Mb deletion at 17q12, which harbors HNF1B, the gene responsible for renal cysts and diabetes syndrome (RCAD), in 18/15,749 patients, including several with ASD, but 0/4,519 controls. We identified additional shared phenotypic features among nine patients available for clinical assessment, including macrocephaly, characteristic facial features, renal anomalies, and neurocognitive impairments. In a large follow-up sample, the same deletion was identified in 2/1,182 ASD/neurocognitive impairment and in 4/6,340 schizophrenia patients, but in 0/47,929 controls (corrected p = 7.37 × 10⁻⁵). These data demonstrate that deletion 17q12 is a recurrent, pathogenic CNV that confers a very high risk for ASD and schizophrenia and show that one or more of the 15 genes in the deleted interval is dosage sensitive and essential for normal brain development and function. In addition, the phenotypic features of patients with this CNV are consistent with a contiguous gene syndrome that extends beyond RCAD, which is caused by HNF1B mutations only.     

A preliminary study of copy number variation in tibetans

Genetic features of Tibetans have been broadly investigated, but the properties of copy number variation (CNV) have not been well examined. To get a preliminary view of CNV in Tibetans, we scanned 29 Tibetan genomes with the Illumina Human-1 M high-resolution genotyping microarray and identified 139 putative copy number variable regions (CNVRs), consisting of 70 deletions, 61 duplications, and 8 multi-allelic loci. Thirty-four of the 139 CNVRs showed differential allele frequencies versus other East-Asian populations, with P values <0.0001. These results indicated a distinct pattern of CNVR allele frequency distribution in Tibetans. The Tibetan CNVRs are enriched for genes in the disease class of human reproduction (such as genes from the DAZ, BPY2, CDY, and HLA-DQ and -DR gene clusters) and biological process categories of "response to DNA damage stimulus" and "DNA repair" (such as RAD51, RAD52, and MRE11A). These genes are related to the adaptive traits of high infant birth weight and darker skin tone of Tibetans, and may be attributed to recent local adaptation. Our results provide a different view of genetic diversity in Tibetans and new insights into their high-altitude adaptation.     

Analysis of copy number variation in the normal human population within a region containing complex segmental duplications on 22q11 using high-resolution array-CGH

A previously detected copy number polymorphism (Ep CNP) in patients affected with neuroectodermal tumors led us to investigate its frequency and length in the normal population. For this purpose, a program called Sequence Allocator was developed and applied for the construction of an array that consisted of unique and duplicated fragments, allowing the assessment of copy number variation within regions of segmental duplications. The average resolution of this array was 11 kb and we determined the size of the Ep CNP to be 290 kb. Analysis of normal controls identified 7.7 and 7.1% gains in peripheral blood and lymphoblastoid cell line (LCL) DNA, respectively, while deletions were found only in the LCL group (7.1%). This array platform allows the detection of DNA copy number variation within regions of pronounced genomic complexity, which constitutes an improvement over available technologies.     

The genetic effect of copy number variations on the risk of type 2 diabetes in a Korean population

Background

Unlike Caucasian populations, genetic factors contributing to the risk of type 2 diabetes mellitus (T2DM) are not well studied in Asian populations. In light of this, and the fact that copy number variation (CNV) is emerging as a new way to understand human genomic variation, the objective of this study was to identify type 2 diabetes–associated CNV in a Korean cohort.

Methodology/Principal Findings

Using the Illumina HumanHap300 BeadChip (317,503 markers), genome-wide genotyping was performed to obtain signal and allelic intensities from 275 patients with type 2 diabetes mellitus (T2DM) and 496 nondiabetic subjects (Total n = 771). To increase the sensitivity of CNV identification, we incorporated multiple factors using PennCNV, a program that is based on the hidden Markov model (HMM). To assess the genetic effect of CNV on T2DM, a multivariate logistic regression model controlling for age and gender was used. We identified a total of 7,478 CNVs (average of 9.7 CNVs per individual) and 2,554 CNV regions (CNVRs; 164 common CNVRs for frequency>1%) in this study. Although we failed to demonstrate robust associations between CNVs and the risk of T2DM, our results revealed a putative association between several CNVRs including chr15:45994758–45999227 (P = 8.6E-04, P^corr = 0.01) and the risk of T2DM. The identified CNVs in this study were validated using overlapping analysis with the Database of Genomic Variants (DGV; 71.7% overlap), and quantitative PCR (qPCR). The identified variations, which encompassed functional genes, were significantly enriched in the cellular part, in the membrane-bound organelle, in the development process, in cell communication, in signal transduction, and in biological regulation.

Conclusion/Significance

We expect that the methods and findings in this study will contribute in particular to genome studies of Asian populations.     

Copy number variation of age-related macular degeneration relevant genes in the Korean population

PurposeStudies that analyzed single nucleotide polymorphisms (SNP) in various genes have shown that genetic factors are strongly associated with age-related macular degeneration (AMD) susceptibility. Copy number variation (CNV) may be an additional type of genetic variation that contributes to AMD pathogenesis. This study investigated CNV in 4 AMD-relevant genes in Korean AMD patients and control subjects.MethodsFour CNV candidate regions located in AMD-relevant genes (VEGFA, ARMS2/HTRA1, CFH and VLDLR), were selected based on the outcomes of our previous study which elucidated common CNVs in the Asian populations. Real-time PCR based TaqMan Copy Number Assays were performed on CNV candidates in 273 AMD patients and 257 control subjects.ResultsThe predicted copy number (PCN, 0, 1, 2 or 3+) of each region was called using the CopyCaller program. All candidate genes except ARMS2/HTRA1 showed CNV in at least one individual, in which losses of VEGFA and VLDLR represent novel findings in the Asian population. When the frequencies of PCN were compared, only the gain in VLDLR showed significant differences between AMD patients and control subjects (p = 0.025). Comparisons of the raw copy values (RCV) revealed that 3 of 4 candidate genes showed significant differences (2.03 vs. 1.92 for VEGFA, p<0.01; 2.01 vs. 1.97 for CFH, p<0.01; 1.97 vs. 2.01, p<0.01 for ARMS2/HTRA1).ConclusionCNVs located in AMD-relevant genes may be associated with AMD susceptibility. Further investigations encompassing larger patient cohorts are needed to elucidate the role of CNV in AMD pathogenesis.     

Genetic variants at 4q21, 4q23 and 12q24 are associated with esophageal squamous cell carcinoma risk in a Chinese population

A recently published genome-wide association study (GWAS) in European populations identified several loci at 4q21, 4q23 and 12q24 that were associated with risk of upper aerodigestive tract (UADT) cancers, including esophageal squamous cell carcinoma (ESCC). In the current study, we conducted a case–control study in a Chinese population including 2,139 ESCC cases and 2,273 controls to evaluate the associations of six reported single nucleotide polymorphisms (SNPs) (rs1494961, rs1229984, rs1789924, rs971074, rs671 and rs4767364) with risk of ESCC. We found significant association with risk of ESCC for four SNPs, including rs1494961 in HEL308 at 4q21 [odds ratio (OR) = 1.15, 95 % confidence interval (CI) = 1.05–1.26], rs1229984 in ADH1B at 4q23 (OR = 1.24, 95 % CI = 1.13–1.36) and rs1789924 near ADH1C at 4q23 (OR = 1.20, 95 % CI = 1.03-1.39), and rs671 in ALDH2 at 12q24 (OR = 0.83, 95 % CI = 0.75–0.91). Combined analysis of these four SNPs showed a significant allele-dosage effect on ESCC risk for individuals with different number of risk alleles (P trend = 2.23 × 10^?11). Compared with individuals with “0–2” risk allele, those carrying “3”, “4” or “5 or more” risk alleles had 1.42-, 1.66-, or 1.76-fold risk of ESCC, respectively. Thus, our findings indicate that rs1494961 at 4q21, rs1229984 and rs1789924 at 4q23, and rs671 at 12q24 may be used as genetic biomarkers for ESCC susceptibility in Chinese population.     

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.     

Increased paternal age and the influence on burden of genomic copy number variation in the general population

Genomic copy number variations (CNVs) and increased parental age are both associated with the risk to develop a variety of clinical neuropsychiatric disorders such as autism, schizophrenia and bipolar disorder. At the same time, it has been shown that the rate of transmitted de novo single nucleotide mutations is increased with paternal age. To address whether paternal age also affects the burden of structural genomic deletions and duplications, we examined various types of CNV burden in a large population sample from the Netherlands. Healthy participants with parental age information (n = 6,773) were collected at different University Medical Centers. CNVs were called with the PennCNV algorithm using Illumina genome-wide SNP array data. We observed no evidence in support of a paternal age effect on CNV load in the offspring. Our results were negative for global measures as well as several proxies for de novo CNV events in this unique sample. While recent studies suggest de novo single nucleotide mutation rate to be dominated by the age of the father at conception, our results strongly suggest that at the level of global CNV burden there is no influence of increased paternal age. While it remains possible that local genomic effects may exist for specific phenotypes, this study indicates that global CNV burden and increased father’s age may be independent disease risk factors.     

Rapid loss of MHC class II variation in a bottlenecked population is explained by drift and loss of copy number variation

Population bottlenecks may reduce genetic variation and potentially increase the risk of extinction. Here, we present the first study to use historic samples to analyse loss of variation at the major histocompatibility complex (MHC), which plays a central role in vertebrate disease resistance. Balancing selection acts on the MHC and could moderate the loss of variation expected from drift; however, in a Wisconsin population of greater prairie-chickens (Tympanuchus cupido), the number of MHC class II B alleles per individual declined by 44% following a population bottleneck, compared to a loss of only 8% at microsatellites. Simulations indicate that drift likely reduced MHC variation at the population level, as well as within individuals by reducing the number of gene copies per individual or by fixing the same alleles across multiple loci. These multiple effects of genetic drift on MHC variation could have important implications for immunity and fitness.     

Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers

Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on a tRNA ligase gene (Migut.N02091; RLG1a) exhibiting unprecedented, and fitness‐relevant, CNV within an annual population of the yellow monkeyflower Mimulus guttatus. RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate‐frequency three‐copy variants of RLG1a (trip+; 5/35 = 14%), and trip+ lines exhibited elevated RLG1a expression under multiple conditions. trip+ carriers, in addition to being over‐represented in late‐flowering and large‐flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rare RLG1a variant (high+) that carries 250–300 copies of RLG1a totalling ~5.7 Mb (20–40% of a chromosome). In the progeny of a high+ carrier, Mendelian segregation of diagnostic alleles and qPCR‐based copy counts indicate that high+ is a single tandem array unlinked to the single‐copy RLG1a locus. In the wild, high+ carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; both p < 0.01), while single‐copy individuals were twice as fecund as either CNV type in a lush year (2016: p < 0.005). Our results demonstrate fluctuating selection on CNVs affecting phenological traits in a wild population, suggest that plant tRNA ligases mediate stress‐responsive life‐history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.     

Plasmid pGA1 from Corynebacterium glutamicum codes for a gene product that positively influences plasmid copy number.

The complete nucleotide sequence (4,826 bp) of the cryptic plasmid pGA1 from Corynebacterium glutamicum was determined. DNA sequence analysis revealed four putative coding regions (open reading frame A [ORFA], ORFA2, ORFB, and ORFC). ORFC was identified as a rep gene coding for an initiator of plasmid replication (Rep) according to the high level of homology of its deduced amino acid sequence with the Rep proteins of plasmids pSR1 (from C. glutamicum) and pNG2 (from Corynebacterium diphtheriae). This function was confirmed by deletion mapping of the minimal replicon of pGA1 (1.7 kb) which contains only ORFC. Deletion derivatives of pGA1 devoid of ORFA exhibited significant decreases in the copy number in C. glutamicum cells and displayed segregational instability. Introduction of ORFA in trans into the cells harboring these deletion plasmids dramatically increased their copy number and segregational stability. The ORFA gene product thus positively influences plasmid copy number. This is the first report on such activity associated with a nonintegrating bacterial plasmid. The related plasmids pGA1, pSR1, and pNG2 lacking significant homology with any other plasmid seem to be representatives of a new group of plasmids replicating in the rolling-circle mode.     

Novel origins of copy number variation in the dog genome

ABSTRACT: BACKGROUND: Copy number variants (CNVs) account for substantial variation between genomes and are a major source of normal and pathogenic phenotypic differences. The dog is an ideal model to investigate mutational mechanisms that generate CNVs as its genome lacks a functional ortholog of the PRDM9 gene implicated in recombination and CNV formation in humans. Here we comprehensively assay CNVs using high-density array comparative genomic hybridization in 50 dogs from 17 dog breeds and 3 gray wolves. RESULTS: We use a stringent new method to identify a total of 430 high-confidence CNV loci, that range in size from 9 kb to 1.6 Mb and span 26.4 Mb, or 1.08%, of the assayed dog genome, overlapping 413 annotated genes. 98% of CNVs observed in each breed are also observed in multiple breeds. CNVs predicted to disrupt gene function are significantly less common than expected by chance. We identify a significant overrepresentation of peaks of GC content, previously shown to be enriched in dog recombination hotspots, in the vicinity of CNV breakpoints. CONCLUSIONS: A number of the CNVs identified by this study are candidates for generating breed-specific phenotypes. Purifying selection seems to be a major factor shaping structural variation in the dog genome, suggesting that many CNVs are deleterious. Localized peaks of GC content appear to be novel sites of CNV formation in the dog genome by non-allelic homologous recombination, potentially activated by the loss of PRDM9. These sequence features may have driven genome instability and chromosomal rearrangements throughout canid evolution.     

Identification of copy number variation and population analysis of the sacred lotus (Nelumbo nucifera)

ABSTRACT

The sacred lotus (Nelumbo nucifera) is widely cultured in East Asia for its horticultural, agricultural, and medicinal values. Although many molecular markers had been used to extrapolate population genetics of the sacred lotus, a study of large variations, such as copy number variation (CNV), are absent up to now. In this study, we applied whole-genome re-sequencing to 24 lotus accessions, and use read depth information to genotype and filter original CNV call. Totally 448 duplications and 4,267 deletions were identified in the final CNV set. Further analysis of population structure revealed that the population structure patterns revealed by CNV and SNP are largely consistent with each other. Our result indicated that deep sequencing followed by genotyping is a quick and straightforward way to mine out CNV from the population, and the CNV along with SNP could enable us to better comprehend the biology of the plant.     

Insertion-deletions burden in copy number polymorphisms of the Tibetan population

BACKGROUND:

Many studies have been conducted to identify either insertions-deletions (inDels) or copy number variations (CNVs) in humans, but few studies have been conducted to identify both of these forms coexisting in the same region.

AIMS AND OBJECTIVES:

To map the functionally significant sites within human genes that are likely to influence human traits and diseases.

MATERIALS AND METHODS:

In this report, we describe an inDel map in the 1051 Tibetan CNV regions obtained through CNV genotyping using Affymetrix Genome-wide single nucleotide polymorphism 6.0 chip. InDel polymorphisms in these copy number polymorphism regions were identified with a computational approach using the 2500 deoxyribonucleic acid sequences obtained from the 1000 Genome Project.

RESULTS:

The study identified a total of 95935 inDels that range from 1 bp to several bps in length which were found scattered across regulatory regions, exons and in introns of genes underlying the CNVs. A study on the distribution of inDels revealed that the majority of inDels were found in coding regions of the genome than the noncoding, while within the genes, inDels in intron regions were more followed by exonic regions and finally the regulatory regions.

CONCLUSION:

Study of inDels in CNV regions contribute to the enhanced understanding of the role played by the two variations and their collective influence on the genome. Further, a collection of these inDel genetic markers will aid in genetic mapping, further understanding of the phenotypic variability, identification of disease genes and in detecting novel CNVs.     

F-box基因拷贝数目变异的机制研究:以12种果蝇为例

拷贝数目变异是一种对表型变异和生物进化具有重要意义的基因组结构变异.以前的研究表明不同物种中F-box基因的拷贝数目差异较大.为了深入探索拷贝数目变异的式样和机制.我们以12个果蝇近缘种为研究对象,分析了F-box基因的系统发育关系、进化式样以及它们在染色体上的位置.结果发现,虽然各个物种中F-box基因的拷贝数目差别不大(42-47个),但是仍然存在着很多引起拷贝数目变异的基因获得和丢失事件.这说明表面上变化不大的拷贝数目在一定程度上掩盖了频繁发生的基因获得和丢失事件.通过比较这些基因在染色体上的位置,发现只有在亲缘关系很近的物种之间才能鉴定出有明显微共线性关系的基因组区段.我们还发现,造成F-box基因拷贝数目增加的主要机制是散在重复和串联重复,而反转录转座和新基因的非编码区起源也是两种值得注意的机制.此外,序列变异导致的外显子边界变化以及外显子丢失是引起拷贝数目减少的两种机制.在12种果蝇的最近共同祖先中,F-box基因的拷贝数目与现存物种基本相似,但是基因的获得和丢失事件使得现存物种中的F-box基因在构成上已经有了明显的差别.对数目变异的式样及其与基因功能的关系的研究表明,拷贝数目变异是F-box基因家族"生与死"的进化在基因组层面的系统反映,并有可能为表型变异提供了原始材料.