A comprehensive profile of DNA copy number variations in a Korean population: identification of copy number invariant regions among Koreans

Copy Number Variation

Differences between genomes can be due to single nucleotide variants (SNPs), translocations, inversions and copy number variants (CNVs, gain or loss of DNA). The latter can range from sub-microscopic events to complete chromosomal aneuploidies. Small CNVs are often benign but those larger than 250 kb are strongly associated with morbid consequences such as developmental disorders and cancer. Detecting CNVs within and between populations is essential to better understand the plasticity of our genome and to elucidate its possible contribution to disease or phenotypic traits.While the link between SNPs and disease susceptibility has been well studied, to date there are still very few published CNV genome-wide association studies; probably owing to the fact that CNV analysis remains a slightly more complex task than SNP analysis (both in term of bioinformatics workflow and uncertainty in the CNV calling leading to high false positive rates and unknown false negative rates). This chapter aims at explaining computational methods for the analysis of CNVs, ranging from study design, data processing and quality control, up to genome-wide association study with clinical traits.

Copy number variants in the sheep genome detected using multiple approaches

Background

Copy number variants (CNVs) are a type of polymorphism found to underlie phenotypic variation, both in humans and livestock. Most surveys of CNV in livestock have been conducted in the cattle genome, and often utilise only a single approach for the detection of copy number differences. Here we performed a study of CNV in sheep, using multiple methods to identify and characterise copy number changes. Comprehensive information from small pedigrees (trios) was collected using multiple platforms (array CGH, SNP chip and whole genome sequence data), with these data then analysed via multiple approaches to identify and verify CNVs.

Results

In total, 3,488 autosomal CNV regions (CNVRs) were identified in this study, which substantially builds on an initial survey of the sheep genome that identified 135 CNVRs. The average length of the identified CNVRs was 19 kb (range of 1 kb to 3.6 Mb), with shorter CNVRs being more frequent than longer CNVRs. The total length of all CNVRs was 67.6Mbps, which equates to 2.7 % of the sheep autosomes. For individuals this value ranged from 0.24 to 0.55 %, and the majority of CNVRs were identified in single animals. Rather than being uniformly distributed throughout the genome, CNVRs tended to be clustered. Application of three independent approaches for CNVR detection facilitated a comparison of validation rates. CNVs identified on the Roche-NimbleGen 2.1M CGH array generally had low validation rates with lower density arrays, while whole genome sequence data had the highest validation rate (>60 %).

Conclusions

This study represents the first comprehensive survey of the distribution, prevalence and characteristics of CNVR in sheep. Multiple approaches were used to detect CNV regions and it appears that the best method for verifying CNVR on a large scale involves using a combination of detection methodologies. The characteristics of the 3,488 autosomal CNV regions identified in this study are comparable to other CNV regions reported in the literature and provide a valuable and sizeable addition to the small subset of published sheep CNVs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2754-7) contains supplementary material, which is available to authorized users.

The Growing Importance of CNVs: New Insights for Detection and Clinical Interpretation

Differences between genomes can be due to single nucleotide variants, translocations, inversions, and copy number variants (CNVs, gain or loss of DNA). The latter can range from sub-microscopic events to complete chromosomal aneuploidies. Small CNVs are often benign but those larger than 500 kb are strongly associated with morbid consequences such as developmental disorders and cancer. Detecting CNVs within and between populations is essential to better understand the plasticity of our genome and to elucidate its possible contribution to disease. Hence there is a need for better-tailored and more robust tools for the detection and genome-wide analyses of CNVs. While a link between a given CNV and a disease may have often been established, the relative CNV contribution to disease progression and impact on drug response is not necessarily understood. In this review we discuss the progress, challenges, and limitations that occur at different stages of CNV analysis from the detection (using DNA microarrays and next-generation sequencing) and identification of recurrent CNVs to the association with phenotypes. We emphasize the importance of germline CNVs and propose strategies to aid clinicians to better interpret structural variations and assess their clinical implications.

Copy number variations in normal karyotype acute myeloid leukaemia and their association with treatment response

Copy number variation (CNV) has been reported to be associated with chemotherapy response, which affects disease prognosis. Here, we determined the frequency of genome-wide cytogenetic CNV aberrations in Korean patients with normal karyotype (NK) acute myeloid leukaemia (AML) and tested whether these genomic variations contribute to differences in Ara-C and anthracycline-based chemotherapy responses. Bone marrow aspirates and blood from 30 previously untreated de novo NK-AML patients were provided at the time of diagnosis for copy number analysis. Possible associations between cytogenetic aberrations and clinical parameters were analysed. CNVs were identified in 23 (76.7%) of the 30 cases tested. Multivariate analyses controlled for other clinical co-variates showed that patients having copy number loss had a decreased probability of complete remission (OR, 0.015 [95% CI, 0-0.737], p = 0.035). Patients who had a copy number gain of more than four regions tended to have shorter event-free survival (EFS) (p = 0.083) with multivariate analysis showing that CNV increase is an independent predictive factor for shorter EFS (HR, 22.104 [95% CI, 1.644-297.157], p = 0.020). In addition, we identified candidate genes that may be involved in Ara-C and anthracycline drug response in Korean patients with NK-AML. These results suggest that CNVs may affect the success of Ara-C and anthracycline-based chemotherapy in Korean patients with NK-AML.

Deoxyribonucleic Acid copy number aberrations in vasospastic angina patients using an array comparative genomic hybridization

BACKGROUND AND OBJECTIVES

Vasospastic angina (VA) is a specific type of coronary artery disease and develops as a result of coronary artery spasm. Recently, a few studies have revealed that VA caused by coronary artery spasm is related to genetic traits. The objective of this study was to use the recently developed technique of array comparative genomic hybridization (CGH) to screen the genetic aberrations of VA.

SUBJECTS AND METHODS

To identify candidate genes that might be causally involved in the pathogenesis of VA, genomic deoxyribonucleic acids were extracted from whole blood of 28 patients with VA who presented at Department of Cardiology at Seoul St. Mary's Hospital, Seoul, Korea. The copy number profiles of these patients was then analyzed using array CGH and reverse transcriptase (RT) quantitative polymerase chain reaction (PCR).

RESULTS

Array CGH revealed gains in 31 different regions, with losses in the 4q35.2, 7q22.1, 10q26.3, 15q11.2, 16p13.11, 17p11.2 and 19q13.3 regions (more than 32% aberration in these regions). Several loci were found to be frequently including gains of 5p and 11q (50% of samples). The most common losses were found in 7q (54% of samples). Copy number aberrations in chromosomal regions were detected and corresponding genes were confirmed by RT quantitative PCR. The fold change levels were highest in the CTDP1 (18q23), HDAC10 (22q13.33), KCNQ1 (11p15.5-p15.4), NINJ2 (12p13.33), NOTCH2 (1p12-p11.2), PCSK6 (15q26.3), SDHA (5p15.33), and MUC17 (7q22.1) genes.

CONCLUSION

Many candidate chromosomal regions that might be related to the pathogenesis of VA were detected by array CGH and should be systematically investigated to establish the causative and specific genes for VA.

Human lymphoblastoid cell lines: a goldmine for the biobankomics era

Biobanking became a necessity for translating genetic discoveries into clinical practice. Approaches to personalized medicine require a new model system for functional and pharmacogenomic studies of a variety of accumulating genetic variations, as well as new research environments such as biobankomics. Human lymphoblastoid cell lines (LCLs) will provide a valuable tool to meet such new demands in the biobankomics era. The National Biobank of Korea (NBK), which is leading the Korea Biobank Project, has a large collection of LCLs derived mostly from population-based cohort samples. Using a special long-term subculture collection of NBK LCLs, biological characteristics of early passage LCLs and terminally immortalized LCLs have been investigated to promote the utilization of LCLs and provide well quality-controlled LCLs for genetic and pharmacogenomic studies. As LCLs have been successfully phenotyped for cytotoxicity in response to various stimulators, including chemotherapeutic agents, environmental chemicals and irradiation, the utility of LCLs will increase in the future. Here, we discuss current and future applications of NBK LCLs for the biobankomics era.

Copy number variation at leptin receptor gene locus associated with metabolic traits and the risk of type 2 diabetes mellitus

Background

Recent efforts have been made to link complex human traits and disease susceptibility to DNA copy numbers. The leptin receptor (LEPR) has been implicated in obesity and diabetes. Mutations and genetic variations of LEPR gene have been discovered in rodents and humans. However, the association of DNA copy number variations at the LEPR gene locus with human complex diseases has not been reported. In an attempt to study DNA copy number variations associated with metabolic traits and type 2 diabetes mellitus (T2DM), we targeted the LEPR gene locus in DNA copy number analyses.

Results

We identified DNA copy number variations at the LEPR gene locus among a Korean population using genome-wide SNP chip data, and then quantified copy numbers of the E2 DNA sequence in the first two exons overlapped between LEPR and LEPROT genes by the quantitative multiplex PCR of short fluorescent fragment (QMPSF) method. Among the non-diabetic subjects (n = 1,067), lower E2 DNA copy numbers were associated with higher fasting glucose levels in men (p = 1.24 × 10^-7) and women (p = 9.45 × 10^-5), as well as higher total cholesterol levels in men (p = 9.96 × 10^-7). In addition, the significant association between lower E2 DNA copy numbers and lower level of postprandial 2hr insulin was evident only in non-diabetic women, whereas some obesity-related phenotypes and total cholesterol level exhibited significant associations only in non-diabetic men. Logistic regression analysis indicated that lower E2 DNA copy numbers were associated with T2DM (odds ratio, 1.92; 95% CI, 1.26~2.96; p < 0.003) in our nested case-control study. Interestingly, the E2 DNA copy number exhibited a negative correlation with LEPR gene expression, but a positive correlation with LEPROT gene expression.

Conclusions

This work suggests that a structural variation at the LEPR gene locus is functionally associated with complex metabolic traits and the risk of T2DM.