- Primary research
- Open Access
- Published:
A comprehensive analysis of copy number variation in a Turkish dementia cohort
Human Genomics volume 15, Article number: 48 (2021)
Abstract
Background
Copy number variants (CNVs) include deletions or multiplications spanning genomic regions. These regions vary in size and may span genes known to play a role in human diseases. As examples, duplications and triplications of SNCA have been shown to cause forms of Parkinson’s disease, while duplications of APP cause early onset Alzheimer’s disease (AD).
Results
Here, we performed a systematic analysis of CNVs in a Turkish dementia cohort in order to further characterize the genetic causes of dementia in this population. One hundred twenty-four Turkish individuals, either at risk of dementia due to family history, diagnosed with mild cognitive impairment, AD, or frontotemporal dementia, were whole-genome genotyped and CNVs were detected. We integrated family analysis with a comprehensive assessment of potentially disease-associated CNVs in this Turkish dementia cohort. We also utilized both dementia and non-dementia individuals from the UK Biobank in order to further elucidate the potential role of the identified CNVs in neurodegenerative diseases.
We report CNVs overlapping the previously implicated genes ZNF804A, SNORA70B, USP34, XPO1, and a locus on chromosome 9 which includes a cluster of olfactory receptors and ABCA1. Additionally, we also describe novel CNVs potentially associated with dementia, overlapping the genes AFG1L, SNX3, VWDE, and BC039545.
Conclusions
Genotyping data from understudied populations can be utilized to identify copy number variation which may contribute to dementia.
Background
Dementia is an umbrella term for symptoms presenting in different diseases such as Alzheimer’s disease (AD), frontotemporal dementia (FTD), and dementia with Lewy bodies (DLB). All of these diseases have significant genetic components. Mendelian early-onset Alzheimer’s disease (EOAD), for example, has been considered to be almost completely genetically determined [1]. However, causative mutations in known genes account for only a small proportion of all cases (~10% of EOAD cases), indicating that other genetic factors, still to be identified, play a substantial role in this form of AD. This missing heritability may be partly explained by genetic elements that are to-date poorly investigated and characterized, such as copy number variants (CNVs).
CNVs are segments of DNA that can have different sizes (ranging from 1 kb to several megabases) for which differences in the number of copies have been identified. These can be copy number gains (duplications or insertional transpositions), losses (deletions), gains or losses of the same locus, or multiallelic or complex rearrangements. The large size of many CNVs causes them to span or partially overlap with various genes, which may lead to unbalanced dosages of those genes. This phenomenon can lead to the development of various diseases. CNVs have been implicated in different neuropsychiatric disorders, such as autism and schizophrenia, and in neurodegenerative diseases [2]. Examples of CNVs that have been proven to cause neurodegenerative diseases include duplications of APP leading to EOAD [3], multiplications of alpha-synuclein (SNCA) causing forms of Parkinson’s disease (PD) and DLB [4], as well as both deletions and duplications of PRKN that are known to cause forms of PD [5].
Only a limited number of studies have assessed the role of CNVs and other structural variants in different dementias. Due to a lack of evidence for segregation, it is often not clear if these are causative, or contributing to the risk of disease. While it is understood that CNVs in known disease-associated genes play a role in the development of diseases, the effects of CNVs in other genes, with no known disease-causing mutations, can be difficult to relate to disease. A genome-wide scan of CNVs in 331 dementia cases (in which >80% of patients had a clinical diagnosis of AD), and 368 controls revealed no CNVs associated with disease using a genome-wide threshold of significance. However, the authors were able to identify a duplication of potential interest encompassing CHRNA7, after screening loci previously reported to be risk regions for neuropsychiatric diseases [6]. The study of 261 non-Mendelian AD families with at least one EOAD case revealed 5 deletions and 5 duplications that segregated with the phenotype. Two of these CNVs overlapped FTD genes (deletion of CHMP2B and duplication of MAPT) [7]. CNV assessments in FTD patients have revealed duplications encompassing MAPT [8] and a partial GRN deletion [9]. In two loci genome-wide associated with AD, structural variants have also been reported to be associated with disease risk: an 8-kb insertion in CR1 that creates the CR1-S isoform harboring an extra set of C3b/C4b binding sites [10], and a 44 base-pair frameshift deletion in ABCA7 that increases AD risk in African Americans and Caribbean Hispanics [11]. Interestingly, this latter association could be ethnic-specific as the frequency of the deletion in non-Hispanic white cohorts is much lower and has so far not resulted in significant associations with risk of disease.
There is increased awareness about the need to represent diverse populations in genetic studies (for a recent review related to AD genetics see [12]). Such contributions include our previous studies of a Turkish dementia cohort where we described mutations in the known Mendelian AD genes [13], a NOTCH3 variant, known to be causative for cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy, in an AD case [14], three homozygous TREM2 variants and a compound heterozygous variant in FTD-like families [15, 16], and no variants or CNVs in TYROBP/DAP12, which encodes a TREM2 ligand [17]. Additionally, this cohort has been screened for variants in genes causative for FTD [18]. We previously reported a preliminary analysis of copy number variation in a subset of this cohort as a poster presentation [19]. Here, we complement this work by assessing the potential role of CNVs in the etiology of dementia in this population.
Results
Family analysis
CNVs were called in trios or quartets for seven Turkish families to analyze segregation of CNVs; pedigree trees are illustrated in Fig. 1.
Pedigree trees for the seven families studied. Squares represent males, circles represent females, and diamonds represent any sex, with numbers representing more than one individual. Diagonal lines indicate that the individual is deceased. Double horizontal bars represent consanguineous marriages and a dotted line suggests consanguinity without support from the full pedigree tree. Black fill indicates affected status, white indicates unaffected and gray indicates risk. Family members genotyped in this study are marked with an asterisk (*). Arrows indicate probands
The sibling pair from family C were both affected with dementia, with the proband receiving a diagnosis of FTD. The proband and three offspring carry a GRN variant ENST00000053867.3:c.415T>C, which co-segregates with a GRN variant ENST00000053867.3:c.-22C>T; both of these variants were previously reported with unclear pathogenicity [18]. The proband died aged 77 years, and the age at last examination was 74 years for the affected sibling genotyped. One duplication overlapping AFG1L and SNX3 was found to be shared by the siblings (Supplementary Fig. 1). Subsequent assessment of this candidate CNV in the entire Turkish dementia cohort revealed no other CNVs in the same locus. There was also no overlap with CNVs called in individuals from the UK Biobank (UKBB).
In family G, where both father and son were diagnosed with FTD (ages at onset 77 and 45 years, respectively), a shared duplication was detected overlapping VWDE. This CNV intersected with a duplication in 1 non-dementia and zero dementia individuals from the UKBB.
Following our filtering criteria, no candidate CNVs were identified in families A, B, D, E, and F. Of these, only family A was a real quartet, with both parents and two offspring genotyped. A probably pathogenic NOTCH3 variant was previously reported in both the father diagnosed with AD (age at onset 62 years) and one offspring who was not genotyped [14]. Ages at last examination were 74, 73, 51, and 50 years for father, mother, son, and daughter, respectively.
From family B, we genotyped a trio with consanguineous parents including the father diagnosed with FTD (age at onset 82 years), and both the mother and son at risk of developing disease, with ages at last examination for the latter being 79 and 49 years, respectively. Based on relatedness between parents, as indicated on the pedigree tree (identity-by-descent 0.1048), mother-father segregation of CNVs was possible and taken into account in the analysis.
In family D, the proband was diagnosed with FTD aged 39 years, the mother was unaffected at last examination aged 79 years, and the additionally genotyped siblings were 37 and 34 years at last examination.
From family E, we genotyped the mother diagnosed with FTD (age at onset 51 years) and daughter at risk (36 years at last examination).
The two siblings in family F were previously reported to carry a PSEN2 variant ENST00000366783.3:c.520A>G reported as probably pathogenic [13], although another report suggested this is likely non-pathogenic [20]. The proband was diagnosed with FTD at 52 years, whilst at last examination, the additionally genotyped sibling was 36 years.
Table 1 summarizes the candidate CNVs identified from this analysis in families C and G.
CNVs in neurogenes or previously reported in neurodegenerative diseases
We next sought to determine if CNVs overlapping well-established neurodegenerative disease genes, or genes implicated in similar analyses were present in this Turkish dementia cohort. Some previously reported findings, such as CR1 duplications, could not be confirmed with our data due to a lack of coverage or resolution at the low copy repeat regions.
As summarized in Table 2, there were CNVs which were both not common in the Database of Genomic Variants (DGV) and had no evidence to refute their association with dementia in the UKBB. Such CNVs included a duplication overlapping ZNF804A detected in one Turkish AD patient; heterozygous deletions in a region overlapping this gene have been previously reported in 2 AD cases, 4 mild cognitive impairment (MCI) cases, and one control [21].
We also detected a duplication overlapping SNORA70B, USP34, and XPO1 in one AD patient; a deletion has been previously reported in an MCI patient [21].
We observed a duplication overlapping GALNTL6 in one AD patient. A previous study reported an association between GALNTL6 CNVs and AD AAO; the significance of this association was lost when corrected for APOE; however, the authors mentioned that GALNTL6 had been previously associated with lipid metabolism, body mass index, and hypertension, suggesting a potential association with AD through a vascular mechanism [22]. Heterozygous deletions intersecting our identified CNV were detected in 1 dementia and 2 non-dementia individuals from the UKBB.
We detected a duplication in one Turkish MCI case overlapping MIR3691; a duplication overlapping this microRNA was reported in one French early-onset AD case [23]. In the UKBB, 1 dementia and 1 non-dementia individual carried heterozygous deletions intersecting this region.
We also observed a heterozygous deletion overlapping ABCA1, LOC286367, NIPSNAP3A, NIPSNAP3B, OR13C2, OR13C5, OR13C8, OR13C9, and OR13D1 in one MCI patient. Heterozygous deletions in a CNV region overlapping these genes have been previously reported in 1 MCI and 4 AD cases [21]. CNVs of unreported type overlapping ABCA1 were described in 3 AD participants and no controls by Swaminathan and colleagues [24]. Although olfactory receptor genes frequently harbor CNVs [25], the heterozygous deletion spanning this region was not considered common in the genome aggregation database structural variants (gnomAD-SV). From the UKBB, only one dementia individual harbored a duplication intersecting this region. Additionally, one deletion spanning this region was reported in DGV, in a Thai infectious disease cohort [26]; these CNVs are collectively illustrated in Fig. 2. We observed instances of copy number discordance for CNVs overlapping the same loci. Potential considerations when interpreting such results include differences in cohort sizes, populations, and diagnoses.
Plot for CNVs overlapping ABCA1 and a cluster of olfactory receptor genes on chromosome 9. The Turkish CNV region is highlighted in green, with the longest transcripts for the overlapping genes plotted below. CNVs extending beyond the Turkish CNV are plotted above the chromosome. AD= Alzheimer’s disease, ADNI= Alzheimer’s Disease Neuroimaging Initiative, gnomAD-SV= genome aggregation database structural variants, MCI= mild cognitive impairment, NCRAD= National Centralized Repository for Alzheimer’s Disease and Related Dementias, NIA-LOAD= National Institute on Aging- late onset Alzheimer’s disease, UKBB= UK Biobank. *This CNV region encompasses heterozygous deletions reported in 1 MCI patient and 4 AD patients. **This CNV was reported in 3 AD patients.
Moreover, we observed CNVs that we defined as common on DGV or report evidence suggesting a lack of association with dementia in the UKBB. One such CNVs was a duplication overlapping HNRNPCL1 detected in an AD patient; CNVs of unreported type overlapping this gene were previously described in 6 AD cases, 4 MCI cases, and zero healthy controls [27].
We detected a duplication overlapping both CYFIP1 and NIPA1 in one AD and one FTD patient. A similar duplication was previously reported in 10 Caribbean Hispanic AD cases and 3 controls, with dosage increase of CYFIP1 and NIPA1 genes confirmed by quantitative PCR [28].
A duplication overlapping DGCR2 was identified in only the affected father of family G. Although this was not common on DGV, there was no significant association with dementia in the UKBB. CNVs of unreported type overlapping this gene have been previously reported in at least 4 AD cases and no controls [29]. This CNV also overlaps PRODH; this encodes proline dehydrogenase which is involved in a key step before the production of glutamate. Elevated proline levels have been previously associated with schizophrenia and bipolar disorder [30], and this has been proposed to be investigated in AD.
We observed a duplication overlapping CNTN6 in one MCI patient; two duplications overlapping this gene were previously reported in a Caribbean Hispanic AD study [28].
We report one heterozygous deletion and two duplications spanning MICA; to note, the heterozygous deletion overlapping HCG26, HCP5, MICA, and PMSP was detected in the to-date unaffected sibling in family F (no CNV overlapping this gene was detected in the affected sibling). This was not reported in the family analysis since deletions (and for that matter, duplications) spanning this region were common on DGV. Moreover, CNVs intersecting this region were not significantly associated with dementia in the UKBB. CNVs of unreported type overlapping this gene were previously reported in four AD patients (and no controls) [24].
CNVs shared within the cohort
Although many candidate CNVs were identified in previous analyses in our study, our approach to test for shared CNVs enabled us to identify additional, potentially novel, CNVs associated with dementia in our cohort, as summarized in Table 3. There were no additional reports in the literature on BC039545, where we observed overlapping deletions in two Turkish FTD patients.
Discussion
This is the first study reporting CNVs in a Turkish dementia cohort. We report novel CNVs that may be causing or be associated with disease, replicate some CNVs previously reported in other dementia cohorts, and utilize a larger dataset of both dementia and non-dementia individuals from the UKBB to further investigate the potential association of candidate CNVs with the studied phenotypes. When comparing CNVs observed in this Turkish cohort to those detected in the UKBB dementia individuals, we observed both CNVs with the same and opposite dosage. The latter is perhaps a more frequent phenomenon than expected, where both deletions and duplications are observed at over 99% of copy number variant regions [31].
From our family analysis, we report potentially novel duplications associated with disease overlapping the genes AFG1L, SNX3, and VWDE. Deletions overlapping BC039545 were additional findings from our analysis of shared genes. As expected, due to lack of previous reports, most CNVs are of uncertain significance. AFG1L encodes a mitochondrial membrane protein and increased protein expression has been shown to induce apoptosis [32]. Downregulated AFG1L expression (also known as LACE1) was reported in AD patients compared to controls [33]. Variants in genes that regulate retromer function, including SNX3, have been associated with AD [34]. Moreover, there is evidence to suggest that SNX3 regulates amyloid-beta production [35]. Retromer dysfunction has been implicated in neurodegenerative diseases including FTD [36]. This duplication, together with the previously reported GRN variants, may increase the genetic risk of FTD. VWDE encodes Von Willebrand Factor D and EGF Domains. This locus has been previously identified as a potential modifier of disease risk in a genome-wide association study of FTD patients with GRN mutations [37].
Our analysis of CNVs overlapping neurogenes and genes implicated in previous CNV studies of similar cohorts contributes further evidence for understanding the potential involvement of such events in disease. One example of a region with conflicting reports overlaps CYFIP1 and NIPA1 where we observed duplications in two Turkish dementia patients. Since previous reports identified CYFIP1 and NIPA1 duplications in three unaffected individuals, CNVs in these genes are most likely not causative for AD; however, potential associations with increased risk of AD were considered [28]. CYFIP1 is linked with the APP translation repressor, FMRP, suggesting that this CNV may alter APP turnover. CYFIP1 deletions have been reported in CNV analyses of schizophrenia [38]. Our candidate CNVs overlapping both CYFIP1 and NIPA1 were predicted as likely benign by ClassifyCNV, and we did not observe a significant association with such CNVs and dementia in the UKBB. However, these duplications in Turkish individuals had a positive likelihood ratio above 1 (Supplementary Table 2).
Deletions spanning a cluster of olfactory receptors on chromosome 9 may increase susceptibility to dementia since there are no structural variants in gnomAD-SV spanning all genes in this cluster; this is interesting considering that olfactory dysfunction is a common prodromal symptom of neurodegenerative diseases [39]. Increased copy number encompassing a different cluster of olfactory receptors was previously associated with earlier AD age at onset; we did not detect any overlap between these genes and CNVs from our broader dementia cohort [40]. The ABCA1 gene, also impacted by these deletions, is a cholesterol transporter from the ATP-binding cassette transporter superfamily. Mutations in this gene cause Tangier disease which is characterized by reduced levels of plasma high density lipoproteins (HDL) [41]. Intriguingly, the long noncoding RNA LOC286367 has been implicated in mediation of ABCA1 expression, at least in the context of propofol treatment [42]. In a recent AD GWAS, rs1800978 in ABCA1 was identified as the lead SNP in a new genome-wide significant locus [43]. Together with the vasoprotective functions of HDL and our observed deletion spanning this gene in an MCI patient, there is increasing evidence suggesting a role of ABCA1 in AD [44]. Also located in this CNV region, NIPSNAP3A gene expression in white matter has been significantly negatively correlated with neurocognitive impairment [45].
We identified a duplication overlapping ZNF804A in one AD patient. Heterozygous deletions overlapping this gene were previously reported in 2 AD cases, 4 MCI cases, and one control [21]. The gene encoding zinc finger binding protein ZNF804A is a well-recognized schizophrenia risk gene knockdown of which has been shown to decrease TYROBP protein expression in a variant CNS catecholaminergic cell line [46, 47]. Variants in TYROBP, encoding an adaptor protein for TREM2, have been reported in EOAD patients [48].
We also detected a duplication overlapping SNORA70B, USP34, and XPO1 in one AD patient; a deletion has been previously reported in an MCI patient [21]. USP34 encodes a ubiquitin-specific protease contributing to the prevention of proteasomal degradation of proteins. USP34 expression was downregulated in hippocampi of patients with AD compared to controls [49]. USP34 has been identified as a YAP1-regulated gene; YAP1 has been proposed as an upstream regulator of AD development [50]. More recently, Usp34 was identified as a differentially expressed gene between lipid-droplet-low and lipid-droplet-rich microglia from aged mice with such changes perhaps contributing to impaired phagocytosis [51]. XPO1 encodes exportin 1 which is involved in the nuclear export of proteins. XPO1 has been reported to modify amyloid-beta toxicity; however, reports for the association of XPO1 variants with AD were not conclusive [52, 53]. Inhibition of XPO1 is not sufficient to promote nuclear localization of TDP-43 and FUS; two proteins which accumulate in the cytoplasm in ALS and FTD [54,55,56]. Potential mechanisms by which XPO1 inhibition may prevent neurodegeneration is by modulating autophagy [57]. Haploinsufficiency of USP34 and XPO1 has been postulated as a mechanism of disease in a patient with mild intellectual disability and cranio-facial dysmorphisms [58].
Limitations of this study include both low sample size and the absence of Turkish controls to cross-reference with and propose new CNVs. However, we utilize data from both dementia and non-dementia UKBB individuals in order to conduct this comparison. Additionally, novel small CNVs are difficult to identify using genotyping data, particularly when the intensity of the probes frequently deviates from the normalized value. To minimize this issue, we used best practices to make the CNV calls; still, there is the caveat that small CNVs may be missed.
There is a clear lack of studies utilizing whole genome genotyping data to detect CNVs in Alzheimer’s disease. This type of data is widely available in the AD field, so there is the ability to conduct similar CNV analyses in much larger datasets to unequivocally replicate existing candidate CNVs and identify new ones. Additionally, large datasets with raw genotyping data available for download, such as the UKBB, are particularly powerful resources, although it should be noted that this dataset is predominantly composed of individuals with European ancestry.
Conclusions
Our analysis of CNVs in a Turkish dementia cohort extends previous reports of CNVs overlapping the genes ZNF804A, SNORA70B, USP34, XPO1, and a locus on chromosome 9 which includes a cluster of olfactory receptors and ABCA1. With regard to ZNF804A, both duplications and deletions have now been reported in AD. The duplication spanning SNORA70B, USP34, and XPO1 in a Turkish AD patient contrasts the deletion previously reported in MCI and previous molecular studies implicating these genes in neurodegeneration. Our report of a Turkish MCI patient with a deletion spanning ABCA1 replicates previous reports in both AD and MCI and increases evidence for a role of such CNVs in disease. We also report novel duplications, potentially associated with FTD, overlapping the genes AFG1L, SNX3, and VWDE from our family analyses. Deletions overlapping BC039545 in two Turkish FTD patients were additional findings from our analysis of shared genes. Future studies in both larger cohorts (including extended family members and new families demonstrating segregation of these CNVs) and additional populations are necessary to replicate these findings for better understanding of their contributions to disease. The detailed genetic study of this cohort has led to identifying several genetic factors associated with dementia in the Turkish population. The current study adds CNVs to the variants with potential roles in different forms of dementia in this population. In addition, given the degree of consanguinity observed in the cohort, we are also systematically assessing the role of variants within extended runs of homozygosity to characterize the genetics of dementia in the Turkish population completely.
Methods
Cohort
We studied a Turkish dementia cohort of one hundred twenty-four individuals with either MCI, AD, FTD, or at risk of such diseases due to family history and age below the age of disease onset. All participants were recruited at the Behavioral Neurology and Movement Disorders Unit outpatient clinic of Istanbul Faculty of Medicine, Istanbul University, and underwent detailed clinical and neuropsychological examination and, if possible, cerebral magnetic resonance imaging (MRI) or positron emission tomography (PET) imaging. Diagnosis of dementia was based on the National Institute of Neurological and Communicative Disorders and Stroke and the Alzheimer’s Disease and Related Disorders Association (NINCDS/ADRDA) [59], FTD was defined following the Lund and Manchester criteria [60], and MCI patients were diagnosed as described by Frank and Petersen [61]. The individuals selected to be genotyped were commonly probands of affected families. Descriptions of cohort phenotypes are in Table 4. Since CNVs can also be involved in phenotype modulation, we included individuals who have had variants previously reported, pathogenic, or otherwise, as listed in Supplementary Table 3. DNA was extracted from whole blood using standard procedures. The study was approved by the Ethics Committee of Istanbul Faculty of Medicine, Istanbul University. A neurologist took the necessary clinical information after obtaining informed written consent from the patients and their participating family members. Consent was provided by the legally authorized representative for subjects unable to consent.
Genotyping
One hundred twenty-four Turkish individuals were genotyped using Illumina’s HumanOmniExpress arrays (Illumina, Inc., CA, USA). Intensity files were analyzed using GenomeStudio (GS) v2.0.4 software (Illumina, Inc., CA, USA) along with the respective manufacturer’s cluster files. Quality control (QC) procedures were performed in GS prior to CNV analysis as described by Jarick and colleagues, including exclusion of SNPs with a GenTrain score below 0.7 [62]. Five genotyped individuals (2 male MCI, 1 female AD, 2 FTD (1 male, 1 female)) were not analyzed for CNVs based on call rate <0.97 or low-quality raw data upon visualization in GS. Therefore, CNVs called from 119 individuals were used in downstream analyses.
Calling CNVs
Copy number variants in autosomes were called using signal intensity files exported by GS, with PennCNV-1.0.5 software [63] (positions in hg19). PennCNV analysis was performed using a PFB file for Human OmniExpress arrays from www.penncnv.openbioinformatics.org, updated on August 16, 2011, and a GC-model adjustment, as per the recommended protocol. Automatic quality control of CNV calls was not used with PennCNV since individuals with low-quality raw data were manually excluded before the analysis. CNVs were filtered to consider those with at least 10 SNPs and at least 50kb in length. Adjacent CNVs were merged based on the default parameters, and CNVs were excluded if they were overlapping telomeres, centromeres, known segmental duplications, and the immunoglobulin, or T cell receptor loci. The final step of the PennCNV analysis involved identifying genes overlapping the CNV calls. Results were further filtered based on complete overlap with CNVs of the same copy number reported in the genome aggregation database structural variants (gnomAD-SV) [64] and part of the curated NCBI Common Structural Variants track which uses the following criteria to define common structural variants: >50bp, variants occurring in at least 100 individuals, allele frequency ≥0.01. The 2020-02-25 release of the Database of Genomic Variants (DGV) was utilized in a similar manner; Turkish CNVs with complete overlap by CNVs of the same type reported in at least 100 individuals in DGV were filtered [65]. Candidate CNVs were visually inspected with GS using cnvPartition v3.2.1 plug-in (Illumina, Inc). In order to determine the presence and copy number of a CNV, measurements of the log R ratio and the B allele frequency were utilized.
CNVs segregating within families
This cohort included seven families; therefore, we first conducted a genome-wide scan for CNVs that segregated with disease within each of these families. We identified families based on information from pedigree trees (Fig. 1). We also calculated identity-by-descent to confirm relatedness within families; based on the same measure, there was no relationship between the 7 families described.
We specifically re-called CNVs in family trios or quartets using PennCNV. The CNV had to be present in at least 1 affected individual and absent in unaffected family members, or only present in unaffected family members which would suggest a protective effect. Family analyses were not conducted on two families with definitely pathogenic variants reported in TREM2. These analyses were only conducted in families where either no variants were identified to be the cause of disease, or in families where variants were identified, but these were of unknown significance for the disease. Candidate CNVs resulting from these family analyses were subsequently assessed in the entire Turkish dementia cohort.
Overlap with “neurogenes” and genes previously reported in CNV analyses of similar dementia cohorts
We assessed our CNV data from PennCNV for CNVs overlapping specific genes known to be related to neurodegenerative diseases [66]; we refer to these as “neurogenes” (Supplementary Table 4).
Additionally, we conducted a literature search for CNVs previously reported in dementia, MCI, AD, or FTD and compiled the respective list of genes (Supplementary Table 5). We combined this list with the “neurogenes” list to obtain a final set of genes (all in autosomes) for which overlapping CNVs were specifically assessed.
CNVs shared within the cohort
We searched for CNVs of the same copy number overlapping the same genes, present in two or more individuals of the studied cohort. Results were further filtered to remove CNVs that were completely localized within 2Mb of autosome ends, based on reported UCSC positions (hg19).
Intersection with dementia and non-dementia individuals from the UK Biobank (UKBB)
We used data from the UK Biobank to validate results. For this, CEL files from 2594 and 4570 UK Biobank dementia or non-dementia individuals, respectively, were downloaded and CNVs were called using both Affymetrix Power Tools [67], and PennCNV, as described above. There were two quality control (QC) steps. The first QC step utilized dish QC (DQC) values which are based on intensities of probe sequences for non-polymorphic genome locations. A DQC value of zero indicates no resolution between the distributions of AT and GC probe contrast values, and a DQC value of 1 indicates perfect resolution. Individuals with a DQC value less than the default threshold of 0.82 were filtered out. The second QC step filtered out individuals with a QC call rate value less than the default threshold of 97%. Following QC, CNVs were detected from 2323 and 4567 UK Biobank individuals with dementia or non-dementia, respectively (hg19). All individuals were actively enrolled in the UKBB study. Exact positions of candidate CNVs from the three analyses were intersected with the CNVs called for UKBB individuals; for the analysis of shared CNVs, the largest encompassing region was used for the intersect. Fisher’s exact test (two-sided) was performed to compare the detection of CNVs overlapping candidate loci between UKBB dementia and non-dementia individuals. When reporting novel dementia CNVs, results were considered significant if p<0.05.
Prediction of CNV pathogenicity and likelihood ratios
We assessed CNV pathogenicity using ClassifyCNV [68]. This tool implements the 2019 American College of Medical Genetics and Genomics classification guidelines [69]. We also report values from likelihood ratio tests for both CNVs reported in Turkish individuals, and overlapping UKBB CNVs with the same copy number. Positive likelihood ratio = \( \frac{\frac{TP}{TP+ FN}}{\frac{FP}{FP+ TN}} \) and negative likelihood ratio = \( \frac{\frac{FN}{TP+ FN}}{\frac{TN}{FP+ TN}} \), where TP = true positive, FN = false negative, TN = true negative, and FP = false positive.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- CNVs:
-
Copy number variants
- AD:
-
Alzheimer’s disease
- FTD:
-
Frontotemporal dementia
- DLB:
-
Dementia with Lewy bodies
- EOAD:
-
Early-onset Alzheimer’s disease
- SNCA:
-
Alpha-synuclein
- PD:
-
Parkinson’s disease
- MCI:
-
Mild cognitive impairment
- MRI:
-
Magnetic resonance imaging
- PET:
-
Positron emission tomography
- NINCDS/ADRDA:
-
National Institute of Neurological and Communicative Disorders and Stroke and the Alzheimer’s Disease and Related Disorders Association
- GS:
-
GenomeStudio
- QC:
-
Quality control
- gnomAD-SV:
-
Genome aggregation database structural variants
- DGV:
-
Database of Genomic Variants
- UKBB:
-
UK Biobank
- DQC:
-
Dish QC
- HDL:
-
High-density lipoproteins
References
Wingo TS, Lah JJ, Levey AI, Cutler DJ. Autosomal recessive causes likely in early-onset Alzheimer disease. Arch Neurol. 2012;69:59–64.
Cook EH Jr, Scherer SW. Copy-number variations associated with neuropsychiatric conditions. Nature. 2008;455:919–23.
Sleegers K, Brouwers N, Gijselinck I, Theuns J, Goossens D, Wauters J, et al. APP duplication is sufficient to cause early onset Alzheimer’s dementia with cerebral amyloid angiopathy. Brain. 2006;129:2977–83.
Singleton AB, Farrer M, Johnson J, Singleton A, Hague S, Kachergus J, et al. alpha-Synuclein locus triplication causes Parkinson’s disease. Science. 2003;302:841.
Simon-Sanchez J, Scholz S, Matarin MM, Fung H-C, Hernandez D, Gibbs JR, et al. Genomewide SNP assay reveals mutations underlying Parkinson disease. Hum Mutat. 2008;29:315–22.
Heinzen EL, Need AC, Hayden KM, Chiba-Falek O, Roses AD, Strittmatter WJ, et al. Genome-wide scan of copy number variation in late-onset Alzheimer’s disease. J Alzheimers Dis. 2010;19:69–77.
Hooli BV, Kovacs-Vajna ZM, Mullin K, Blumenthal MA, Mattheisen M, Zhang C, et al. Rare autosomal copy number variations in early-onset familial Alzheimer’s disease. Mol Psychiatry. 2014;19:676–81.
Rovelet-Lecrux A, Hannequin D, Guillin O, Legallic S, Jurici S, Wallon D, et al. Frontotemporal dementia phenotype associated with MAPT gene duplication. J Alzheimers Dis. 2010;21:897–902.
Blauwendraat C, Wilke C, Simón-Sánchez J, Jansen IE, Reifschneider A, Capell A, et al. The wide genetic landscape of clinical frontotemporal dementia: systematic combined sequencing of 121 consecutive subjects. Genet Med. 2018;20:240–9.
Brouwers N, Van Cauwenberghe C, Engelborghs S, Lambert J-C, Bettens K, Le Bastard N, et al. Alzheimer risk associated with a copy number variation in the complement receptor 1 increasing C3b/C4b binding sites. Mol Psychiatry. 2012;17:223–33.
Cukier HN, Kunkle BW, Vardarajan BN, Rolati S, Hamilton-Nelson KL, Kohli MA, et al. ABCA7 frameshift deletion associated with Alzheimer disease in African Americans. Neurol Genet. 2016;2:e79.
Dehghani N, Bras J, Guerreiro R. How understudied populations have contributed to our understanding of Alzheimer’s disease genetics. Brain. 2021. https://doi.org/10.1093/brain/awab028.
Lohmann E, Guerreiro RJ, Erginel-Unaltuna N, Gurunlian N, Bilgic B, Gurvit H, et al. Identification of PSEN1 and PSEN2 gene mutations and variants in Turkish dementia patients. Neurobiol Aging. 1850;2012(33):e17–27.
Guerreiro RJ, Lohmann E, Kinsella E, Brás JM, Luu N, Gurunlian N, et al. Exome sequencing reveals an unexpected genetic cause of disease: NOTCH3 mutation in a Turkish family with Alzheimer’s disease. Neurobiol Aging. 2012;33:1008 e17–23.
Guerreiro RJ, Lohmann E, Brás JM, Gibbs JR, Rohrer JD, Gurunlian N, et al. Using exome sequencing to reveal mutations in TREM2 presenting as a frontotemporal dementia-like syndrome without bone involvement. JAMA Neurol. 2013;70:78–84.
Guerreiro R, Bilgic B, Guven G, Brás J, Rohrer J, Lohmann E, et al. Novel compound heterozygous mutation in TREM2 found in a Turkish frontotemporal dementia-like family. Neurobiol Aging. 2013;34:2890 e1–5.
Darwent L, Carmona S, Lohmann E, Guven G, Kun-Rodrigues C, Bilgic B, et al. Mutations in TYROBP are not a common cause of dementia in a Turkish cohort. Neurobiol Aging. 2017;58:240 e1–240.e3.
Guven G, Lohmann E, Bras J, Gibbs JR, Gurvit H, Bilgic B, et al. Mutation frequency of the major frontotemporal dementia genes, MAPT, GRN and C9ORF72 in a Turkish cohort of dementia patients. PLoS One. 2016;11:e0162592.
Dehghani N, Kun-Rodrigues C, Foster K, Hanagasi H, Lohmann E, Guven G, et al. Analysis of copy number variation in a Turkish dementia cohort. Alzheimers Dement, Wiley. 2020:16 Available from: https://onlinelibrary.wiley.com/doi/10.1002/alz.044868.
Cruchaga C, Haller G, Chakraverty S, Mayo K, Vallania FLM, Mitra RD, et al. Rare variants in APP, PSEN1 and PSEN2 increase risk for AD in late-onset Alzheimer’s disease families. PLoS One. 2012;7:e31039.
Guffanti G, Torri F, Rasmussen J, Clark AP, Lakatos A, Turner JA, et al. Increased CNV-region deletions in mild cognitive impairment (MCI) and Alzheimer’s disease (AD) subjects in the ADNI sample. Genomics. 2013;102:112–22.
Szigeti K, Kellermayer B, Lentini JM, Trummer B, Lal D, Doody RS, et al. Ordered subset analysis of copy number variation association with age at onset of Alzheimer’s disease. J Alzheimers Dis. 2014;41:1063–71.
Boscher E, Husson T, Quenez O, Laquerrière A, Marguet F, Cassinari K, et al. Copy number variants in miR-138 as a potential risk factor for early-onset Alzheimer’s disease. J Alzheimers Dis. 2019;68:1243–55.
Swaminathan S, Shen L, Kim S, Inlow M, West JD, Faber KM, et al. Analysis of copy number variation in Alzheimer’s disease: the NIALOAD/ NCRAD Family Study. Curr Alzheimer Res. 2012;9:801–14.
Young JM, Endicott RM, Parghi SS, Walker M, Kidd JM, Trask BJ. Extensive copy-number variation of the human olfactory receptor gene family. Am J Hum Genet. 2008;83:228–42.
Suktitipat B, Naktang C, Mhuantong W, Tularak T, Artiwet P, Pasomsap E, et al. Copy number variation in Thai population. PLoS One. 2014;9:e104355.
Swaminathan S, Kim S, Shen L, Risacher SL, Foroud T, Pankratz N, et al. Genomic copy number analysis in Alzheimer’s disease and mild cognitive impairment: an ADNI study. Int J Alzheimers Dis. 2011;2011:729478.
Ghani M, Pinto D, Lee JH, Grinberg Y, Sato C, Moreno D, et al. Genome-wide survey of large rare copy number variants in Alzheimer’s disease among Caribbean hispanics. G3. 2012;2:71–8.
Swaminathan S, Huentelman MJ, Corneveaux JJ, Myers AJ, Faber KM, Foroud T, et al. Analysis of copy number variation in Alzheimer’s disease in a cohort of clinically characterized and neuropathologically verified individuals. PLoS One. 2012;7:e50640.
Clelland CL, Drouet V, Rilett KC, Smeed JA, Nadrich RH, Rajparia A, et al. Evidence that COMT genotype and proline interact on negative-symptom outcomes in schizophrenia and bipolar disorder. Transl Psychiatry. 2016;6:e891.
Li YR, Glessner JT, Coe BP, Li J, Mohebnasab M, Chang X, et al. Rare copy number variants in over 100,000 European ancestry subjects reveal multiple disease associations. Nat Commun. 2020;11:255.
Cesnekova J, Spacilova J, Hansikova H, Houstek J, Zeman J, Stiburek L. LACE1 interacts with p53 and mediates its mitochondrial translocation and apoptosis. Oncotarget. 2016;7:47687–98.
Jiang Z, Tan G, Wang Z. Comprehensive analysis reveals A six-gene signature and associated drugs in Alzheimer disease [Internet]. Research Square. 2020; Available from: https://www.researchsquare.com/article/rs-30210/v1.
Vardarajan BN, Bruesegem SY, Harbour ME, Inzelberg R, Friedland R, St George-Hyslop P, et al. Identification of Alzheimer disease-associated variants in genes that regulate retromer function. Neurobiol Aging. 2012;33:2231 e15–2231.e30.
Xu S, Nigam SM, Brodin L. Overexpression of SNX3 decreases amyloid-β peptide production by reducing internalization of amyloid precursor protein. Neurodegener Dis. 2018;18:26–37.
Carosi JM, Denton D, Kumar S, Sargeant TJ. Retromer dysfunction at the nexus of tauopathies. Cell Death Differ. 2021. https://doi.org/10.1038/s41418-020-00727-2.
Pottier C, Zhou X, Perkerson RB 3rd, Baker M, Jenkins GD, Serie DJ, et al. Potential genetic modifiers of disease risk and age at onset in patients with frontotemporal lobar degeneration and GRN mutations: a genome-wide association study. Lancet Neurol. 2018;17:548–58.
Lew AR, Kellermayer TR, Sule BP, Szigeti K. Copy number variations in adult-onset neuropsychiatric diseases. Curr Genomics. 2018;19:420–30.
Zou Y-M, Lu D, Liu L-P, Zhang H-H, Zhou Y-Y. Olfactory dysfunction in Alzheimer’s disease. Neuropsychiatr Dis Treat. 2016;12:869–75.
Shaw CA, Li Y, Wiszniewska J, Chasse S, Zaidi SNY, Jin W, et al. Olfactory copy number association with age at onset of Alzheimer disease. Neurology. 2011;76:1302–9.
Remaley AT, Rust S, Rosier M, Knapper C, Naudin L, Broccardo C, et al. Human ATP-binding cassette transporter 1 (ABC1): genomic organization and identification of the genetic defect in the original Tangier disease kindred. Proc Natl Acad Sci U S A. 1999;96:12685–90.
Ma X, Wang T, Zhao Z-L, Jiang Y, Ye S. Propofol suppresses proinflammatory cytokine production by increasing ABCA1 expression via mediation by the long noncoding RNA LOC286367. Mediators Inflamm. 2018;2018:8907143.
Bellenguez C, Küçükali F, Jansen I, Andrade V, Moreno-Grau S, Amin N, et al. New insights on the genetic etiology of Alzheimer’s and related dementia [Internet]. bioRxiv. medRxiv. 2020; Available from: http://medrxiv.org/lookup/doi/10.1101/2020.10.01.20200659.
Button EB, Robert J, Caffrey TM, Fan J, Zhao W, Wellington CL. HDL from an Alzheimer’s disease perspective. Curr Opin Lipidol. 2019;30:224–34.
Levine AJ, Miller JA, Shapshak P, Gelman B, Singer EJ, Hinkin CH, et al. Systems analysis of human brain gene expression: mechanisms for HIV-associated neurocognitive impairment and common pathways with Alzheimer’s disease. BMC Med Genomics. 2013;6:4.
Chang H, Xiao X, Li M. The schizophrenia risk gene ZNF804A: clinical associations, biological mechanisms and neuronal functions. Mol Psychiatry. 2017;22:944–53.
Zhou Y, Dong F, Lanz TA, Reinhart V, Li M, Liu L, et al. Interactome analysis reveals ZNF804A, a schizophrenia risk gene, as a novel component of protein translational machinery critical for embryonic neurodevelopment. Mol Psychiatry. 2018;23:952–62.
Pottier C, Ravenscroft TA, Brown PH, Finch NA, Baker M, Parsons M, et al. TYROBP genetic variants in early-onset Alzheimer’s disease. Neurobiol Aging. 2016;48:222 e9–222.e15.
Ho Kim J, Franck J, Kang T, Heinsen H, Ravid R, Ferrer I, et al. Proteome-wide characterization of signalling interactions in the hippocampal CA4/DG subfield of patients with Alzheimer’s disease. Sci Rep. 2015;5:11138.
Xu M, Zhang D-F, Luo R, Wu Y, Zhou H, Kong L-L, et al. A systematic integrated analysis of brain expression profiles reveals YAP1 and other prioritized hub genes as important upstream regulators in Alzheimer’s disease. Alzheimers Dement. 2018;14:215–29.
Marschallinger J, Iram T, Zardeneta M, Lee SE, Lehallier B, Haney MS, et al. Lipid-droplet-accumulating microglia represent a dysfunctional and proinflammatory state in the aging brain. Nat Neurosci. 2020;23:194–208.
Treusch S, Hamamichi S, Goodman JL, Matlack KES, Chung CY, Baru V, et al. Functional links between Aβ toxicity, endocytic trafficking, and Alzheimer’s disease risk factors in yeast. Science. 2011;334:1241–5.
Rosenthal SL, Wang X, Demirci FY, Barmada MM, Ganguli M, Lopez OL, et al. Beta-amyloid toxicity modifier genes and the risk of Alzheimer’s disease. Am J Neurodegener Dis. 2012;1:191–8.
Archbold HC, Jackson KL, Arora A, Weskamp K, Tank EM-H, Li X, et al. TDP43 nuclear export and neurodegeneration in models of amyotrophic lateral sclerosis and frontotemporal dementia. Sci Rep. 2018;8:4606.
Ederle H, Funk C, Abou-Ajram C, Hutten S, Funk EBE, Kehlenbach RH, et al. Nuclear egress of TDP-43 and FUS occurs independently of Exportin-1/CRM1. Sci Rep. 2018;8:7084.
Pinarbasi ES, Cağatay T, Fung HYJ, Li YC, Chook YM, Thomas PJ. Active nuclear import and passive nuclear export are the primary determinants of TDP-43 localization. Sci Rep. 2018;8:7083.
Silvestrini MJ, Johnson JR, Kumar AV, Thakurta TG, Blais K, Neill ZA, et al. Nuclear export inhibition enhances HLH-30/TFEB activity, autophagy, and lifespan. Cell Rep. 2018;23:1915–21.
Fannemel M, Barøy T, Holmgren A, Rødningen OK, Haugsand TM, Hansen B, et al. Haploinsufficiency of XPO1 and USP34 by a de novo 230 kb deletion in 2p15, in a patient with mild intellectual disability and cranio-facial dysmorphisms. Eur J Med Genet. 2014;57:513–9.
McKhann G, Drachman D, Folstein M, Katzman R, Price D, Stadlan EM. Clinical diagnosis of Alzheimer’s disease: report of the NINCDS-ADRDA Work Group under the auspices of Department of Health and Human Services Task Force on Alzheimer's Disease. Neurology. 1984;34:939–44.
Clinical and neuropathological criteria for frontotemporal dementia. The Lund and Manchester Groups. J Neurol Neurosurg Psychiatry. 1994;57:416–8.
Frank AR, Petersen RC. Mild cognitive impairment. Handb Clin Neurol. 2008;89:217–21.
Jarick I, Volckmar A-L, Pütter C, Pechlivanis S, Nguyen TT, Dauvermann MR, et al. Genome-wide analysis of rare copy number variations reveals PARK2 as a candidate gene for attention-deficit/hyperactivity disorder. Mol Psychiatry. 2014;19:115–21.
Wang K, Li M, Hadley D, Liu R, Glessner J, Grant SFA, et al. PennCNV: an integrated hidden Markov model designed for high-resolution copy number variation detection in whole-genome SNP genotyping data. Genome Res. 2007;17:1665–74.
Collins RL, Brand H, Karczewski KJ, Zhao X, Alföldi J, Francioli LC, et al. A structural variation reference for medical and population genetics. Nature. 2020;581:444–51.
MacDonald JR, Ziman R, Yuen RKC, Feuk L, Scherer SW. The Database of Genomic Variants: a curated collection of structural variation in the human genome. Nucleic Acids Res. 2014;42:D986–92.
Kun-Rodrigues C, Orme T, Carmona S, Hernandez DG, Ross OA, Eicher JD, et al. A comprehensive screening of copy number variability in dementia with Lewy bodies. Neurobiol Aging. 2019;75:223 e1–223.e10.
Kendall KM, Rees E, Escott-Price V, Einon M, Thomas R, Hewitt J, et al. Cognitive performance among carriers of pathogenic copy number variants: analysis of 152,000 UK Biobank Subjects. Biol Psychiatry. 2017;82:103–10.
Gurbich TA, Ilinsky VV. ClassifyCNV: a tool for clinical annotation of copy-number variants [Internet]. Cold Spring Harbor Laboratory. 2020 [cited 2020 Nov 13]. p. 2020.07.20.213215. Available from: https://www.biorxiv.org/content/10.1101/2020.07.20.213215v1.full-text
Riggs ER, Andersen EF, Cherry AM, Kantarci S, Kearney H, Patel A, et al. Technical standards for the interpretation and reporting of constitutional copy-number variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics (ACMG) and the Clinical Genome Resource (ClinGen). Genet Med. 2020;22:245–57.
Acknowledgements
This research has been conducted using data from UK Biobank, a major biomedical database (www.ukbiobank.ac.uk).
We thank the Van Andel Institute Bioinformatics and Biostatistics Core, especially Zach Madaj, for their statistical expertise.
Funding
There was no specific funding for this study.
Author information
Authors and Affiliations
Contributions
ND, RG, and JB designed the study and wrote the manuscript. ND, CKR, CG, KF, RG, and JB performed data analyses and interpreted the data. GG, HH, EL, BS, HG, and BB collected and characterized samples for inclusion in the study. All authors provided critical feedback and helped shape the research, analysis and manuscript. The authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The study was approved by the Ethics Committee of Istanbul Faculty of Medicine, Istanbul University. A neurologist took the necessary clinical information after obtaining informed written consent from the patients and their participating family members. Consent was provided by the legally authorized representative for subjects unable to consent.
Consent for publication
Details relating to individuals are entirely unidentifiable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Supplementary Figure 1
. GS plot illustrating duplications spanning AFGL1/LACE1 and SNX3 in a sibling pair from family C.
Additional file 2: Supplementary Table 1
. Shared CNVs of the same copy number spanning the same genes in two or more individuals.
Additional file 3: Supplementary Table 2
. Likelihood ratios for CNVs reported from analysis of both Turkish dementia cohort and UK Biobank.
Additional file 4: Supplementary Table 3
. Previously reported definitely and potentially pathogenic variants in individuals included in this cohort.
Additional file 5: Supplementary Table 4
. List of neurogenes.
Additional file 6: Supplementary Table 5
. Genes reported to overlap CNVs associated with dementia.
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.
About this article
Cite this article
Dehghani, N., Guven, G., Kun-Rodrigues, C. et al. A comprehensive analysis of copy number variation in a Turkish dementia cohort. Hum Genomics 15, 48 (2021). https://doi.org/10.1186/s40246-021-00346-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40246-021-00346-z
Keywords
- Copy number variants
- Dementia
- Genotyping