Skip to main content
Download PDF

Exploring quantitative traits-associated copy number deletions through reanalysis of UK10K consortium whole genome sequencing cohorts

BMC Genomics volume 24, Article number: 787 (2023) Cite this article

Abstract

Objectives

We performed comprehensive association analyses of common high-confidence gnomAD-reported copy number deletions (CNDs) with 60 quantitative traits from UK10K consortium WGS data.

Methods

The study made use of data generated by the UK10K Consortium. UK10K consortium WGS data consist of TwinsUK (n = 1754, middle-aged females) and ALSPAC (n = 1867, birth to adolescence) cohorts. UK10K consortium called 18,739 CNDs (hg19) with GenomeSTRiP software. After filtering out variants with minor allele frequency < 0.05 or HWE P < 1.0 × 10− 6, 1222 (TwinsUK) and 1211 (ALSPAC) CNDs remained for association analyses with 60 normalized quantitative traits.

Results

We identified 23 genome-wide significant associations at 13 loci, among which 2 associations reached experiment-wide significance. We found that two common deletions in chromosome 4, located between WDR1 and ZNF518B (23.3 kb, dbVar ID:nssv15888957, 4:10211262–10,234,569 and 9.8 kb, dbVar ID:nssv15888975, 4:10392422–10,402,191), were associated with uric acid levels (P = 5.23 × 10− 11 and 2.29 × 10− 8, respectively). We also discovered a novel deletion spanning chromosome 18 (823 bp, dbVar ID: nssv15841628, 8:74347187–74,348,010) associated with low HDL cholesterol levels (P = 4.15 × 10− 7). Additionally, we observed two red blood cell traits-associated loci with genome-wide significance, a 13.2 kb deletion in 7q22.1 (nssv15922542) and a 3.7 kb deletion in 12q24.12 (nssv15813226), both of which were located in regions previously reported to be associated with red blood cell traits. Two deletions in 11q11 (nssv15803200 and nssv15802240), where clusters of multiple olfactory receptor genes exist, and a deletion (nssv15929560) upstream to DOCK5 were associated with childhood obesity. Finally, when defining Trait-Associated copy number Deletions (TADs) as CNDs with phenotype associations at sub-threshold significance (P < 10− 3), we identified 157 (97.5%) out of 161 TADs in non-coding regions, with a mean size of 4 kb (range: 209 − 47,942 bp).

Conclusion

We conducted a reanalysis of the UK10K Whole Genome Sequencing cohort, which led to the identification of multiple high confidence copy number deletions associated with quantitative traits. These deletions have standard dbVar IDs and replicate previous findings, as well as reveal novel loci that require further replication studies.

Peer Review reports

Introduction

Copy number variations (CNVs) are unbalanced rearrangements of DNA segments that involve at least 50 nucleotides with gain (duplication) or loss (deletion) of DNA contents [1, 2]. It is estimated that 4.8–9.5% of the genome contributes to CNVs and encompass more polymorphic base pairs than single nucleotide variations (SNPs) [3]. CNVs are reported to have a global effect on the transcriptome [4, 5] and contribute to various diseases and traits through numerous mechanisms disrupting protein-coding genes and cis-regulatory regions [6].

Microarrays including comparative genomic hybridization (CGH) and SNP-based arrays have been the primary approaches for CNV detection. However, with the advent of next-generation sequencing (NGS) technologies, new tools can detect CNVs from NGS data [7]. Compared to microarrays, NGS can detect smaller variants and have both higher sensitivity and the ability to provide accurate sequence-level breakpoint resolution, whereas array-based approaches have a limited resolution capacity [8].

Recently, a population reference for medical or population genetic studies, the Genome Aggregation Database (gnomAD), started to provide the reference of sequence-resolved structural variations like CNVs constructed from 14,891 whole genomes across diverse global populations [9]. Although whole genome sequencing (WGS) is rapidly becoming the predominant technology in large-scale human disease studies or biobank-based population cohort studies [10], few studies exist that investigated the association of high-confidence SVs from WGS with diseases or traits [11].

The UK10K consortium is one of the earliest large-scale WGS-based cohorts [12]. WGS was performed on a total of 3781 participants in two UK cohorts, the TwinsUK cohort and the ALSPAC cohort. TwinsUK is the UK’s largest twin cohort, comprising more than 14,000 twins who were predominantly female (82%) and middle-aged with the average of 59 years [13] and the ALSPAC cohort is a birth cohort in which more than 14,000 pregnant women and their children were followed for about 20 years from birth and health information from childhood to adolescence was collected [14].

In the present study, we performed a comprehensive association study of common high confidence copy number deletions (CNDs) reported in gnomAD with 60 distinct quantitative traits based on the UK10K consortium WGS cohort data. The traits were grouped into several broad classes such as lipid and glucose metabolism, body composition, cardiovascular health, serum chemistry, lung function, etc. We were able to replicate previously reported associations of CNDs with quantitative traits and identified new CNDs that showed associations with traits. Our findings suggest that incorporating structural variations in association studies could broaden our understanding of the connections between genotypes and phenotypes.

Methods

Whole genome sequencing and phenotype data of UK10K project

The access to whole genome sequencing and phenotype data of UK10K project (https://www.uk10k.org) was granted on 14th April 2018 (request ID:6410). Data were downloaded for 1754 subjects in TwinsUK cohort and 1867 subjects in ALSPAC (Avon Longitudinal Study of Parents and Children) cohort.

In the TwinsUK and ALSPAC cohorts, 48 and 42 unique phenotypes were identified, with 30 overlapping phenotypes between them (see Supplementary table S1). These phenotypes or traits were linked to various biological functions, including lipid and glucose metabolism, body composition, pulmonary function, cardiovascular health, blood cell composition, and serum biochemistry.

The TwinsUK cohort primarily consisted of adult female twins with phenotypes measured at 56.0 (standard deviation = 12.3) years of age on average. On the other hand, in ALSPAC cohort, phenotypes were longitudinally collected from childhood to adolescence. Specifically, glucose metabolic traits such as fasting glucose, insulin levels, insulin secretory function (HOMA-Beta), and insulin resistance (HOMA-IR) [15] were measured at 15 years of age. Lipid traits like total cholesterol, VLDL, HDL, LDL, triglycerides, APO-A1, and APO-B were measured at 9 years of age. Leptin, adiponectin, CRP, homocysteine, and hemoglobin were profiled at age 9. Pulmonary function tests were carried out at 8 years of age and blood pressures were measured at 9 years of age and heart rate was measured at 13 years of age. Body mass index (BMI) was serially measured at age 7, 9, 13, and 15 years of age and other body compositions like total fat mass and total lean mass by dual-energy x-ray absorptiometry (DEXA), and waist-to-hip ratio were measured at age 9.

Copy number deletion dataset and annotation

Using whole genome sequencing data at a mean depth ×7, deletions ranging in size from 100 bp to 1 Mb were identified by UK10K consortium using the GenomeSTRiP(version1.04.1068) software [16] as described in the UK10K consortium flagship paper [12]. A total of 18,739 copy number deletions profiled for each cohort were downloaded in multi-sample vcf files (Supplementary table S2). Chromosome positions followed hg19 version of human reference genome. We annotated copy number deletions with VEP [17] and AnnotSV softwares [18, 19].

Association analysis

We excluded copy number deletions on sex chromosomes and retained 18,056 variants. Further, we filtered out variants with a minor allele frequency < 0.05 or Hardy-Weinberg Equilibrium P < 1.0 × 10 − 6 to focus our analysis on common variants. This resulted in a total of 1222 and 1211 common copy number deletions for association analysis in the TwinsUK and ALSPAC cohorts, respectively (Supplementary table S3). We applied the same filtering criteria for the association study of single nucleotide variations (SNV) and small insertions/deletions (Indel) with more than 5 million SNVs and Indels in each cohort. Prior to conducting association analysis, continuous phenotype variables were transformed by rank-based inverse normal transformation. We performed GWAS across all normalized traits adjusting for age and sex and assumed additive, dominant, or recessive genetic models (see Supplementary table S4 for GWAS summary statistics for CNDs). Hail version 0.2 software (https://github.com/hail-is/hail) was used for variant filtering and association analyses. After association analysis, we restricted analysis to CNDs with confirmed and comparable population allele frequencies reported in gnomAD SV database [9] and reported in dbVar [20]. In the meta-analysis of the two cohorts for 30 common phenotypes, METAL was used to conduct a fixed-effects meta-analysis and p values for association were combined [21]. Local association plot was generated with LocusZoom [22]. Other statistical analyses were conducted with R v4.0.2 (https://www.r-project.org/). Network diagrams were drawn with Cytoscape software [23].

Results

Association analysis of common copy number deletions

The genome-wide significance threshold was set as 4.09 × 10− 5 and 4.13 × 10− 5 after Bonferroni correction at level α = 0.05 accounting for 1222 and 1211 variables in TwinsUK and ALSPAC cohorts, respectively. Furthermore, the experiment-wide significance threshold was set as 2.84 × 10− 7 and 3.28 × 10− 7 correcting for 48 and 42 independent variables within TwinsUK and ALSPAC cohorts, respectively, and three assumed models of genetic inheritance.

In the TwinsUK cohort, two loci associated with uric acid levels surpassed the experiment-wide significance threshold (Table 1). A common 23.3 kb deletion in chromosome 4 (dbVar ID:nssv15888957, 4:10211262–10,234,569) with deletion allele frequency 0.69 (0.68 among gnomAD Europeans) was associated with reduced uric acid levels, whereas adjacent common 9.8 kb deletion in chromosome 4 (dbVar ID:nssv15888975, 4:10392422–10,402,191) with deletion allele frequency 0.24 (0.25 among gnomAD Europeans) was associated with elevated uric acid levels (Fig. 1). A GWAS including SNVs and Indels with uric acid levels in TwinsUK cohort confirmed the significant peak on the chromosome 4 (Supplementary Figure S1). These two copy number deletions exhibited moderate linkage disequilibrium (R2 = 0.69) and placed between WDR1(WD Repeat Domain 1) and ZNF518B (Zinc Finger Protein 518B) genes near the urate transporter SLC2A9 (Solute Carrier Family 2 Member 9) gene which is reported to be associated with uric acid levels and gout. Among common SNVs in the region, the A allele of rs874432, an intron variant of SLC2A9, was associated with lower uric acid levels with the most statistical significance (P = 1.98 × 10− 26). After adjusting for age and rs874432, nssv15888957 (adjusted P = 3.69 × 10− 4) and nssv15888975 (adjusted P = 2.26 × 10− 3) remained associated with serum uric acid levels at nominal significance. At combined analysis with adjacent SNVs in the region, nssv15888957 formed a 254 kb-ranging haplotype block in strong linkage disequilibrium with 224 SNVs (R2 > 0.8)(Fig. 1). The haplotype block is placed in active regulatory element region that harbors two enhancers, GH04J010179 (4:10180825–10,188,676) and GH04J010344(4:10346574–10,348,589), curated in GeneHancer database, that target SLC2A9, WDR1, and ZNF518B in many cell types including kidney and urothelial cells [24].

Table 1 Copy number deletions associated with traits at experiment-wide or genome-wide significance
Full size table

In addition to experiment-wide significant associations, there were 21 genome-wide significant associations for 11 loci in two cohorts in different genetic models (Table 1). At TwinsUK cohort, there were 8 genome-wide significance associations for 6 copy number deletions.

A common (deletion allele frequency 0.29) 823 bp deletion that spans on chromosome 18 (dbVar ID: nssv15841628, 18:74347187–74,348,010) was newly associated with low HDL cholesterol and ApoA1 levels at genome-wide significance, assuming recessive inheritance (Table 1). This deletion is located in intron 3 of lncRNA, LINC01927 (Fig. 2) and GWAS of SNVs and Indels with HDL levels also showed sub-threshold peak in the region (Supplementary Figure S2). Low HDL cholesterol < 40 mg/dL (1.034mmol/L) is a well-known major cardiovascular risk factor, and we categorized the subjects with HDL < 1.034mmol/L according to the number of deletion alleles. Subjects with two deletion alleles had higher proportion of subjects with HDL < 1.034mmol/L versus the rest with odd ratio = 1.90.

A 5 kb deletion in STPG2 and STPG2-AS1 (dbVar ID: nssv15893432, 4:98354274–98,359,233) was associated with lower total lean mass (Table 1 and supplementary figure  S3). nssv15893432 was not only associated with total lean mass but also lower BMI (P = 2.73 × 10− 3, β=-0.16), height (P = 1.04 × 10− 3, β=-0.17), FEV1(P = 4.09 × 10− 3, β=-0.12) and FVC(P = 4.24 × 10− 3, beta=-0.13) suggesting the link between poor lung function and sarcopenia [25].

Two red blood cell traits-associated copy number deletions were found near SNV loci with previously reported associations with red cell traits [26]. A 13.2 kb deletion in 7q22.1 with loss of function on ZAN (Zonadhesin) gene and downstream of EPO gene was associated with MCH. Another 3.7 kb deletion in 12q24.12 associated with MCV is found in intron 5 of ATXN2 (Ataxin 2) gene which is associated with multiple red cell traits (Table 1 and supplementary figures S4-S5).

It is of note that a common deletion (dbVar ID: nssv15913056, 6:134589734–134,594,596) in intron 1 of SGK1 (Serum/Glucocorticoid Regulated Kinase 1) which is involved in sodium regulation was associated with higher sodium levels (P = 1.74 × 10− 5, Supplementary figure S6).

In the ALSPAC cohort, 13 genome-wide significance associations were found for 5 copy number deletions (Table 1). A 6.96 kb copy number deletion on chromosome 11 (dbVar ID: nssv15802240, 11:55031591–55,038,554) with loss of E3 ubiquitin-protein ligase TRIM48 (tripartite motif containing 48) gene was associated with childhood obesity (Table 1). nssv15802240 was in strong linkage disequilibrium with adjacent SNVs and formed a long haplotype block harboring multiple genes in olfactory receptor gene cluster. The block contains a copy number deletion involving OR4P4, OR4S2, and OR4C6 genes (Fig. 3) which is associated with early extreme obesity in previous studies using genotyping arrays and suggesting the link between olfactory dysfunction and obesity [27]. In the present study with whole genome sequencing, the region is a 79.6 kb common copy number deletion (dbVar ID:nssv15802291, 11:55364273–55,443,871, gnomAD deletion allele frequency 0.27 in Europeans) that spans olfactory receptor family 4 genes (OR4C11, OR4P4, OR4S2, and OR4C6) and was also nominally associated with early obesity (BMI at 15 years of age, β = 0.11, dominant, P = 3.54 × 10− 2) in the ALSPAC cohort. We also found a 1.59 kb copy number deletion (dbVar ID: nssv15803200, 11:49759283–49,760,872) within intron 2 of GRM5P1(GRM5 pseudogene 1) associated with early obesity. It also formed a long-ranging haplotype block with adjacent SNVs (Fig. 3) that involves olfactory receptor family 4 genes (OR4C12 and OR4C13).

Another copy number deletion associated with early central obesity (truncal fat mass) was an 18.5 kb deletion in intergenic region on chromosome 8 (dbVar ID: nssv15929560) that was placed upstream of DOCK5 (Dedicator Of Cytokinesis 5) gene (Fig. 3). DOCK5 has been reported associated with severe obesity and regulate hepatic insulin resistance [28, 29].

A 1.1 kb deletion on chromosome 1 (dbVar ID: nssv15849193, 1:234318630–234,319,762) that disrupts the thiamine transporter SLC35F3 (Solute Carrier Family 35 Member F3) was associated with higher LDL cholesterol in childhood, suggesting the link between thiamine deficiency and dyslipidemia (supplementary figure S7).

Meta-analysis and sub-threshold associations replicated in the other cohort

Although phenotypes in children and adult women may have different characteristics, we performed meta-analysis for 30 quantitative traits shared by the two cohorts from GWAS summary statistics. Two loci were discovered surpassing the Bonferroni-corrected genome-wide significance threshold at combined meta-analysis p value (Table 2). A 3 kb deletion in GUSBP1 (GUSB Pseudogene 1) enhancer on chromosome 5 (dbVar ID: nssv15898820, 5: 21611834–21614796) was associated with waist circumference and total fat mass and a 0.8 kb deletion on chromosome 3 (dbVar ID: nssv15883024, 3:141572125–141572896) in enhancer region targeting multiple genes was associated with anemia. Additionally, we also sought to find sub-threshold associations (P < 10− 3) that are replicated in the other cohort as traits can persist regardless of age and found 4 sub-threshold associations for 4 copy number deletions replicated in the other cohort (Table 2). The association of 4.7 kb deletion in 3’-UTR region of FAM149A gene (dbVar ID: nssv15897230, 4:187093507–187,098,168) with reduced lung function (FVC and FEV1) in adult women was replicated in childhood cohort. It is interesting that the deletion was associated with lower total lean mass (P = 3.47 × 10− 3) and body mass indices (BMI_7, P = 1.04 × 10− 3; BMI_9, P = 7.13 × 10− 3; BMI_13, P = 2.48 × 10− 2; BMI_15, P = 7.30 × 10− 3) in children, suggesting the genetic link between childhood sarcopenia and adult poor lung function (supplementary figure S8) [30].

Table 2 Copy number deletions with sub-threshold associations with traits replicated at the other cohort
Full size table

Trait-associated copy number deletions (TAD)

We defined Trait-Associated copy number Deletions (TADs) as copy number deletions associated with any phenotypes with at least sub-threshold significance (P < 10− 3). Out of 239 associations (Supplementary table S5), we discovered 161 Trait-Associated copy number Deletions (TADs) in the two cohorts (Supplementary table S6). The mean size of TADs was 4 kb (range: 209 − 47,942 bp). Overall, 157 (97.5%) out of 161 TADs were found in non-coding regions (Fig. 4). There were 60 TADs annotated in gene regions with 56 TADs in introns and only 4 TADs resulting in loss of function deletions of coding regions. Out of 56 TADs in intron regions, 28 TADs were annotated to have enhancer activity. Additionally, among 101 TADs in intergenic regions, 34 deletions were located in enhancer regions.

Pleiotropy of trait-associated copy number deletions (TADs)

The comprehensive exploration of the associations between multiple phenotypes from childhood to adult and CNDs enabled us to investigate the pleiotropy of TADs by network analysis. The full network (Supplementary Fig. S4) illustrates that adult and childhood phenotypes are generally distinct, suggesting that adult and childhood phenotypes are distinct with respect to TADs. The largest connected core of the full network is shown in Fig. 5. Adult and childhood lung functions are connected by 719 bp deletion on chromosome 4 (1,602,940–1,603,659) associated with obstructive lung functions (FEV1/FVC ratio in childhood and FEV1 in adult). It is concordant with the previous epidemiologic study reporting that being in the lowest quartile for lung function at age 7 might have long-term consequences for the development of COPD by middle age [31]. A 1.2 kb deletion on chromosome 18 (75,267,000–75,268,164) was associated with both elevated adult serum uric acid levels and stronger grip strength in children. Epidemiologic data suggests the association of hand grip strength, a marker of sarcopenia, and serum uric acid levels in middle aged adults [32], and it warrants further study if muscle mass in children predicts serum uric acid levels in adults.

Discussion

We reanalyzed UK10K WGS cohorts’ data to uncover copy number deletions associated with 60 quantitative traits and found 23 genome-wide significant associations at 13 copy number deletions. Some of the associations replicate previously suggested links between copy number deletions and quantitative traits: the uric acid association signals around the SLC2A9 locus or childhood obesity related 11q11 region. Others such as red blood cell trait associations are found in regions tagged by red cell traits associated SNVs from cohorts with huge number of subjects with microarray data. The others are novel associations of copy number deletions with traits such as sodium, HDL cholesterol, and lean body mass that require further validation in separate cohorts. We also observed pleiotropy of common copy number deletions that might underlie complex epidemiological associations.

One of the strengths of the study is that we annotated copy number deletions with dbVar reference IDs and population allele frequencies reported in gnomAD structural variation (SV) database with definite breakpoint positions. Since the first descriptions of the genome-wide presence of copy number variations(CNV) in the human genome in year 2004 [1, 2], attempts have been made to build public CNV databases [3]. However, compared to SNPs or Indels that are well curated in dbSNP database with non-redundant and unique Reference SNP (rs) numbers, the mutational spectra of SVs are only recently beginning to be reported based on high-depth whole genome sequencing [9, 10]. As WGS is rapidly becoming the predominant technology in large-scale human studies with progressively decreasing sequencing cost, we expect that trait or disease mapping studies involving structural variations will report reference IDs such as dbVar IDs, enabling curated databases of phenotype related structural variations.

Another strength of the study is that we analyzed association signals by combining SNVs and copy number deletions. The combined analysis found long-ranging haplotype blocks that overlap with regulatory region or gene clusters. For example, around the SLC2A9 gene on chromosome 4 with very high association signal with serum uric acid [33], nssv15888957 formed a 254 kb-ranging haplotype block in strong linkage disequilibrium with 224 SNVs (R2 > 0.8) spanning putative enhancer region. With respect to childhood obesity, we found several haplotype blocks of CNDs and SNVs that span multiple olfactory receptor gene clusters. Considering studies that elevated BMI is known to be associated with olfactory dysfunction and mouse experiment suggesting that ectopic olfactory receptor activation reverses obesity, the existence of early obesity related CNDs and SNVs in strong LD in 11q11 olfactory gene cluster region might corroborate the role of olfactory dysfunction in obesity [34, 35].

We also found novel trait-associated CNDs in functionally relevant genes which had no previous reports of associations with SNVs or Indels. For example, CND in Intron 1 of the SGK1 gene (nssv15913056) was associated with elevated sodium levels. SGK1 is a gene that plays an important role in regulating sodium concentration, and it is involved in the regulation of sodium reabsorption by modulating the activity of the epithelial sodium channel (ENaC) in the distal tubules of the kidney [36]. However, no sodium level associated SNVs or Indels at or near SGK1 gene have been registered in the GWAS catalogue (https://www.ebi.ac.uk/gwas/) database despite the definite involvement of SGK1 in sodium metabolism. Further human genetic epidemiologic studies that validate the novel associations in the present study are warranted.

The majority (97.5%) out of 161 Trait-Associated copy number Deletions (TADs) was placed in non-coding regions such as introns or intergenic areas. This genomic impact of common copy number deletions, or the bias away from genes, is consistent with study by Conrad et al. that reported 95.3% of copy number deletions were placed in non-coding regions [37]. Out of 60 TADs annotated in gene regions, 56 TADs were in introns and structural variation in intron might be associated with gene expression (CNV-eQTL) [38] or regulate genes with tissue-specific activity [39]. It is interesting that 42% of TADs were found in non-enhancer intergenic regions and what further genomic features are associated with quantitative traits is the subject of further study.

The study is unique in that it analyzed both adult women and children’s cohorts and we could study the pleiotropy of CNDs with respect to adult and childhood traits. For example, nssv15897230 was associated with childhood sarcopenia and poor lung function in adults that is found in previous epidemiologic studies [30]. In the pleiotropy analysis, 719 bp deletion on chromosome 4 (dbVar ID: nssv15887243, 4:1602940–1,603,659) was associated with obstructive lung functions in childhood and adults (FEV1/FVC ratio in childhood and FEV1 in adult). The longitudinal associations between childhood traits and their impacts on disease in later adult life are important epidemiological questions and further mendelian randomization studies using SNVs and copy number deletions profiled in UK10K WGS cohorts may give further insight.

There are several limitations of the study: (1) The study included relatively moderate number of subjects, (2) only the associations with traits were analyzed because more clinically relevant phenotypes like disease were not available, (3) the replication study was limited as the two cohorts, TwinsUK and ALSPAC, represented two clinically distinct subject groups, middle-aged women versus children, and (4) genome-wide or experiment-wide significant associations need further replicate validation analyses and the pleiotropy analysis could contain false positive findings that require attention in interpretation. With the prospect that millions of genomes will be sequenced in the coming years from national initiatives we expect there will be available replicate WGS cohorts with high-confidence SV calls.

In conclusion, the reanalysis of UK10K WGS cohorts with combined analysis of SNVs and copy number deletions uncovered multiple novel and replicating associations with multiple quantitative traits. The recent explosive growth in whole genome sequence data from multiple biobanks that call both SNVs and structural variations might enable a comprehensive understanding into the genotype to phenotype relationships.

Fig. 1
figure 1

(A) QQ plot and (B) Manhattan plot of uric acid associated copy number deletions. Two deletions passing experiment-wide significance are marked with red circles. (C) nssv15888957 and nssv15888975 copy number deletions are associated with lower and higher serum uric acid levels, respectively. (D) Fine mapping analysis of deletions and SNVs in linkage disequilibrium with nssv15888957. Variants were colored by the category of linkage disequilibrium R2 values with respect to nssv15888957. It demonstrates that nssv15888957 forms a 254 kb-ranging haplotype block in strong linkage disequilibrium with 224 SNVs. Upper panel shows LocusZoom plot and lower panel is from UCSC genome browser annotated with histone marks and GeneHancer reported regulatory elements

Full size image
Fig. 2
figure 2

(A) QQ plot and (B) Manhattan plot of HDL cholesterol associated copy number deletions. nssv15841628 is marked with red circle. (C) nssv15841628 deletion is associated with lower HDL and ApoA1 levels. (D) Fine mapping analysis of deletions and SNVs in linkage disequilibrium with nssv15841628. Variants were colored by the category of linkage disequilibrium R2 values with respect to nssv15841628. The start position of nssv15841628 is marked with diamond

Full size image
Fig. 3
figure 3

Fine mapping analysis of childhood obesity associated deletions. Variants were colored by the category of linkage disequilibrium R2 values with respect to copy number deletions. Haplotypes including SNVs in linkage disequilibrium with nssv15802240 associated with BMI (A) and nssv15803200 associated with total fat mass (B) formed a long haplotype block, encompassing multiple olfactory receptor genes in 11q11 region. (C) nssv15929560 associated with truncal fat mass is in upstream of DOCK5 gene. The start positions of deletions are marked with diamonds

Full size image
Fig. 4
figure 4

Genomic impact of 161 Trait-Associated copy number Deletions (TADs). Most (97.5%) TADs were found in non-coding regions. Out of 157 deletions in non-exon regions, 62 (39.5%) TADs were annotated to have enhancer activity

Full size image
Fig. 5
figure 5

The comprehensive associations between multiple childhood to adult phenotypes based on pleiotropy of CNDs. Phenotypes (blue nodes) and copy number deletions (orange nodes) are connected by edges if they are associated with at least sub-threshold significance (P < 10− 3) at TwinsUK(purple edges) and ALSPAC(green edges) cohorts with positive (continuous line) or negative (dotted line) β values. Two CND nodes that connect categories of phenotypes are colored with red

Full size image

Data availability

The datasets supporting the conclusions of this article are included within the article and the supplementary materials. Any additional data are available on request to the corresponding author.

References

  1. Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y, Scherer SW, Lee C. Detection of large-scale variation in the human genome. Nat Genet. 2004;36(9):949–51.

    Article  CAS  PubMed  Google Scholar 

  2. Sebat J, Lakshmi B, Troge J, Alexander J, Young J, Lundin P, Maner S, Massa H, Walker M, Chi M, et al. Large-scale copy number polymorphism in the human genome. Science. 2004;305(5683):525–8.

    Article  CAS  PubMed  Google Scholar 

  3. Zarrei M, MacDonald JR, Merico D, Scherer SW. A copy number variation map of the human genome. Nat Rev Genet. 2015;16(3):172–83.

    Article  CAS  PubMed  Google Scholar 

  4. Haraksingh RR, Snyder MP. Impacts of variation in the human genome on gene regulation. J Mol Biol. 2013;425(21):3970–7.

    Article  CAS  PubMed  Google Scholar 

  5. Henrichsen CN, Chaignat E, Reymond A. Copy number variants, Diseases and gene expression. Hum Mol Genet. 2009;18(R1):R1–8.

    Article  CAS  PubMed  Google Scholar 

  6. Hurles ME, Dermitzakis ET, Tyler-Smith C. The functional impact of structural variation in humans. Trends Genet. 2008;24(5):238–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Pirooznia M, Goes FS, Zandi PP. Whole-genome CNV analysis: advances in computational approaches. Front Genet. 2015;6:138.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Pang AW, MacDonald JR, Pinto D, Wei J, Rafiq MA, Conrad DF, Park H, Hurles ME, Lee C, Venter JC, et al. Towards a comprehensive structural variation map of an individual human genome. Genome Biol. 2010;11(5):R52.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Collins RL, Brand H, Karczewski KJ, Zhao X, Alfoldi J, Francioli LC, Khera AV, Lowther C, Gauthier LD, Wang H, et al. A structural variation reference for medical and population genetics. Nature. 2020;581(7809):444–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Halldorsson BV, Eggertsson HP, Moore KHS, Hauswedell H, Eiriksson O, Ulfarsson MO, Palsson G, Hardarson MT, Oddsson A, Jensson BO, et al. The sequences of 150,119 genomes in the UK Biobank. Nature. 2022;607(7920):732–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Chen L, Abel HJ, Das I, Larson DE, Ganel L, Kanchi KL, Regier AA, Young EP, Kang CJ, Scott AJ, et al. Association of structural variation with cardiometabolic traits in finns. Am J Hum Genet. 2021;108(4):583–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Consortium UK, Walter K, Min JL, Huang J, Crooks L, Memari Y, McCarthy S, Perry JR, Xu C, Futema M, et al. The UK10K project identifies rare variants in health and Disease. Nature. 2015;526(7571):82–90.

    Article  Google Scholar 

  13. Moayyeri A, Hammond CJ, Hart DJ, Spector TD. The UK Adult Twin Registry (TwinsUK Resource). Twin Res Hum Genet. 2013;16(1):144–9.

    Article  PubMed  Google Scholar 

  14. Boyd A, Golding J, Macleod J, Lawlor DA, Fraser A, Henderson J, Molloy L, Ness A, Ring S, Davey Smith G. Cohort Profile: the ‘children of the 90s’--the index offspring of the Avon Longitudinal Study of parents and children. Int J Epidemiol. 2013;42(1):111–27.

    Article  PubMed  Google Scholar 

  15. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985;28(7):412–9.

    Article  CAS  PubMed  Google Scholar 

  16. Handsaker RE, Korn JM, Nemesh J, McCarroll SA. Discovery and genotyping of genome structural polymorphism by sequencing on a population scale. Nat Genet. 2011;43(3):269–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GR, Thormann A, Flicek P, Cunningham F. The Ensembl variant effect predictor. Genome Biol. 2016;17(1):122.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Geoffroy V, Herenger Y, Kress A, Stoetzel C, Piton A, Dollfus H, Muller J. AnnotSV: an integrated tool for structural variations annotation. Bioinformatics. 2018;34(20):3572–4.

    Article  CAS  PubMed  Google Scholar 

  19. Geoffroy V, Guignard T, Kress A, Gaillard JB, Solli-Nowlan T, Schalk A, Gatinois V, Dollfus H, Scheidecker S, Muller J. AnnotSV and knotAnnotSV: a web server for human structural variations annotations, ranking and analysis. Nucleic Acids Res. 2021;49(W1):W21–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Lappalainen I, Lopez J, Skipper L, Hefferon T, Spalding JD, Garner J, Chen C, Maguire M, Corbett M, Zhou G, et al. DbVar and DGVa: public archives for genomic structural variation. Nucleic Acids Res. 2013;41(Database issue):D936–941.

    CAS  PubMed  Google Scholar 

  21. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Boughton AP, Welch RP, Flickinger M, VandeHaar P, Taliun D, Abecasis GR, Boehnke M. LocusZoom.js: interactive and embeddable visualization of genetic association study results. Bioinformatics. 2021;37(18):3017–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Fishilevich S, Nudel R, Rappaport N, Hadar R, Plaschkes I, Iny Stein T, Rosen N, Kohn A, Twik M, Safran M et al. GeneHancer: genome-wide integration of enhancers and target genes in GeneCards. Database (Oxford) : 2017, 2017.

  25. Moon JH, Kong MH, Kim HJ. Implication of Sarcopenia and Sarcopenic obesity on lung function in healthy Elderly: using Korean National Health and Nutrition Examination Survey. J Korean Med Sci. 2015;30(11):1682–8.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Ganesh SK, Zakai NA, van Rooij FJ, Soranzo N, Smith AV, Nalls MA, Chen MH, Kottgen A, Glazer NL, Dehghan A, et al. Multiple loci influence erythrocyte phenotypes in the CHARGE Consortium. Nat Genet. 2009;41(11):1191–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Jarick I, Vogel CI, Scherag S, Schafer H, Hebebrand J, Hinney A, Scherag A. Novel common copy number variation for early onset extreme obesity on chromosome 11q11 identified by a genome-wide analysis. Hum Mol Genet. 2011;20(4):840–52.

    Article  CAS  PubMed  Google Scholar 

  28. El-Sayed Moustafa JS, Eleftherohorinou H, de Smith AJ, Andersson-Assarsson JC, Alves AC, Hadjigeorgiou E, Walters RG, Asher JE, Bottolo L, Buxton JL, et al. Novel association approach for variable number tandem repeats (VNTRs) identifies DOCK5 as a susceptibility gene for severe obesity. Hum Mol Genet. 2012;21(16):3727–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Lai Y, Zhao A, Tan M, Yang M, Lin Y, Li S, Song J, Zheng H, Zhu Z, Liu D, et al. DOCK5 regulates energy balance and hepatic insulin sensitivity by targeting mTORC1 signaling. EMBO Rep. 2020;21(2):e49473.

    Article  CAS  PubMed  Google Scholar 

  30. Peralta GP, Fuertes E, Granell R, Mahmoud O, Roda C, Serra I, Jarvis D, Henderson J, Garcia-Aymerich J. Childhood body composition trajectories and adolescent lung function. Findings from the ALSPAC study. Am J Respir Crit Care Med. 2019;200(1):75–83.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Bui DS, Burgess JA, Lowe AJ, Perret JL, Lodge CJ, Bui M, Morrison S, Thompson BR, Thomas PS, Giles GG, et al. Childhood lung function predicts adult Chronic Obstructive Pulmonary Disease and Asthma-Chronic Obstructive Pulmonary Disease Overlap Syndrome. Am J Respir Crit Care Med. 2017;196(1):39–46.

    Article  PubMed  Google Scholar 

  32. Kawamoto R, Ninomiya D, Kasai Y, Kusunoki T, Ohtsuka N, Kumagi T, Abe M. Serum uric acid is positively Associated with Handgrip Strength among Japanese Community-Dwelling Elderly women. PLoS ONE. 2016;11(4):e0151044.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Merriman TR. An update on the genetic architecture of hyperuricemia and gout. Arthritis Res Ther. 2015;17(1):98.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Faour M, Magnan C, Gurden H, Martin C. Olfaction in the context of obesity and Diabetes: insights from animal models to humans. Neuropharmacology. 2022;206:108923.

    Article  CAS  PubMed  Google Scholar 

  35. Wu C, Hwang SH, Jia Y, Choi J, Kim YJ, Choi D, Pathiraja D, Choi IG, Koo SH, Lee SJ. Olfactory receptor 544 reduces adiposity by steering fuel preference toward fats. J Clin Invest. 2017;127(11):4118–23.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Chen SY, Bhargava A, Mastroberardino L, Meijer OC, Wang J, Buse P, Firestone GL, Verrey F, Pearce D. Epithelial sodium channel regulated by aldosterone-induced protein sgk. Proc Natl Acad Sci U S A. 1999;96(5):2514–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, et al. Origins and functional impact of copy number variation in the human genome. Nature. 2010;464(7289):704–12.

    Article  CAS  PubMed  Google Scholar 

  38. Rigau M, Juan D, Valencia A, Rico D. Intronic CNVs and gene expression variation in human populations. PLoS Genet. 2019;15(1):e1007902.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Borsari B, Villegas-Miron P, Perez-Lluch S, Turpin I, Laayouni H, Segarra-Casas A, Bertranpetit J, Guigo R, Acosta S. Enhancers with tissue-specific activity are enriched in intronic regions. Genome Res. 2021;31(8):1325–36.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The study makes use of data generated by the UK10K Consortium, derived from samples from UK10K_COHORTS_TWINSUK (EGA Study ID: EGAS00001000108) and UK10K_COHORT_ALSPAC (EGA Study ID: EGAS00001000090). A full list of the investigators who contributed to the generation of the data is available from www.UK10K.org. Funding for UK10K was provided by the Wellcome Trust under award WT091310.

Funding

This work was supported by a grant (No. 18-2018-0006) from the Seoul National University Bundang Hospital Research Fund, and by a grant of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number: HR22C1363).

Author information

Authors and Affiliations

  1. Precision Medicine Center, Future Innovation Research Division, Seoul National University Bundang Hospital, 173-82, Gumi-ro, Bundang-gu, Seongnam, Gyeonggi-do, 13620, South Korea

    Sejoon Lee, Jinho Kim & Jung Hun Ohn

  2. Department of Pathology, Seoul National University Bundang Hospital, 173-82, Gumi-ro, Bundang-gu, Seongnam, Gyeonggi-do, 13620, South Korea

    Sejoon Lee

  3. Department of Laboratory Medicine, Seoul National University Bundang Hospital, 173-82, Gumi-ro, Bundang-gu, Seongnam, Gyeonggi-do, 13620, South Korea

    Jinho Kim

  4. Department of Internal Medicine, Seoul National University Bundang Hospital, 173-82, Gumi-ro, Bundang-gu, Seongnam, Gyeonggi-do, 13620, South Korea

    Jung Hun Ohn

  5. Department of Internal Medicine, College of Medicine, Seoul National University, 103, Daehak-ro, Jongno-gu, Seoul, 03080, South Korea

    Jung Hun Ohn

Authors
  1. Sejoon Lee
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Jinho Kim
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jung Hun Ohn
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

J.H.O. thought of the study and S.L., J.K., J.H.O. wrote the main manuscript text. S.L. and J.H.O. analyzed UK10K data. All authors reviewed the manuscript.

Corresponding author

Correspondence to Jung Hun Ohn.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

Not applicable as the study is a re-analysis of UK10K WGS dataset.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lee, S., Kim, J. & Ohn, J.H. Exploring quantitative traits-associated copy number deletions through reanalysis of UK10K consortium whole genome sequencing cohorts. BMC Genomics 24, 787 (2023). https://doi.org/10.1186/s12864-023-09903-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09903-3

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us