Outcome of array CGH analysis for 255 subjects with intellectual disability and search for candidate genes using bioinformatics

Array CGH enables the detection of pathogenic copy number variants (CNVs) in 5–15% of individuals with intellectual disability (ID), making it a promising tool for uncovering ID candidate genes. However, most CNVs encompass multiple genes, making it difficult to identify key disease gene(s) underlying ID etiology. Using array CGH we identified 47 previously unreported unique CNVs in 45/255 probands. We prioritized ID candidate genes using five bioinformatic gene prioritization web tools. Gene priority lists were created by comparing integral genes from each CNV from our ID cohort with sets of training genes specific either to ID or randomly selected. Our findings suggest that different training sets alter gene prioritization only moderately; however, only the ID gene training set resulted in significant enrichment of genes with nervous system function (19%) in prioritized versus non-prioritized genes from the same de novo CNVs (7%, p < 0.05). This enrichment further increased to 31% when the five web tools were used in concert and included genes within mitogen-activated protein kinase (MAPK) and neuroactive ligand-receptor interaction pathways. Gene prioritization web tools enrich for genes with relevant function in ID and more readily facilitate the selection of ID candidate genes for functional studies, particularly for large CNVs.     

Microdeletion and microduplication syndromes

The widespread use of whole genome analysis based on array comparative genomic hybridization in diagnostics and research has led to a continuously growing number of microdeletion and microduplication syndromes (MMSs) connected to certain phenotypes. These MMSs also include increasing instances in which the critical region can be reciprocally deleted or duplicated. This review catalogues the currently known MMSs and the corresponding critical regions including phenotypic consequences. Besides the pathogenic pathways leading to such rearrangements, the different detection methods and their limitations are discussed. Finally, the databases available for distinguishing between reported benign or pathogenic copy number alterations are highlighted. Overall, a review of MMSs that previously were also denoted "genomic disorders" or "contiguous gene syndromes" is given.     

Reference-unbiased copy number variant analysis using CGH microarrays

Comparative genomic hybridization (CGH) microarrays have been used to determine copy number variations (CNVs) and their effects on complex diseases. Detection of absolute CNVs independent of genomic variants of an arbitrary reference sample has been a critical issue in CGH array experiments. Whole genome analysis using massively parallel sequencing with multiple ultra-high resolution CGH arrays provides an opportunity to catalog highly accurate genomic variants of the reference DNA (NA10851). Using information on variants, we developed a new method, the CGH array reference-free algorithm (CARA), which can determine reference-unbiased absolute CNVs from any CGH array platform. The algorithm enables the removal and rescue of false positive and false negative CNVs, respectively, which appear due to the effects of genomic variants of the reference sample in raw CGH array experiments. We found that the CARA remarkably enhanced the accuracy of CGH array in determining absolute CNVs. Our method thus provides a new approach to interpret CGH array data for personalized medicine.     

Appropriateness of array‐CGH in the ADHD clinics: A comparative study

Attention deficit hyperactivity disorder (ADHD) is one of the most common neurodevelopmental disorder with a worldwide prevalence of about 5%. The disorder is characterized by inattentive, hyperactive and impulsive behavior and is often comorbid with other neuropsychiatric conditions. Array comparative genomic hybridization (array‐CGH) testing has been proved to be useful to detect chromosomal aberrations in several neuropsychiatric conditions including autism spectrum disorders (ASD) and intellectual disability (ID). The usefulness of array‐CGH in the ADHD clinics is still debated and no conclusive evidence has been reached to date. We performed array‐CGH in 98 children and adolescents divided in two similarly sized groups according to the clinical diagnosis: (a) one group diagnosed with ADHD as primary diagnosis; (b) the other group in which ADHD was co‐morbid with ASD and/or ID. We detected pathogenetic and likely pathogenetic copy number variants (CNVs) in 12% subjects in which ADHD was co‐morbid with autism and/or intellectual disability and in 8.5% subjects diagnosed with ADHD as primary diagnosis. Detection of CNVs of unknown clinical significance was similar in the two groups being 27% and 32%, respectively. Benign and likely benign CNVs accounted for 61% and 59.5% in the first and second group, respectively. Differences in the diagnostic yield were not statistically significant between the two groups (P > .05). Our data strongly suggest that array‐CGH (a) is a valuable diagnostic tool to detect clinically significant CNVs in individuals with ADHD even in the absence of comorbidity with ASD and/or ID and (b) should be implemented routinely in the ADHD clinics.     

Molekulare Karyotypisierung in der klinischen Anwendung

Microarray technology for the detection of putative pathological submicroscopic copy number variants (CNV) has become a standard tool in the field of molecular cytogenetics in recent years. In addition to the identification of somatic CNVs in tumour genetics this technology is increasingly used for the analysis of constitutional CNVs in patients with developmental delay. Array-based genomic hybridisation increases sensitivity in comparison to more conventional technologies such as comparative genomic hybridisation (CGH). Recent developments now allow a genome-wide detection of submicroscopic chromosomal alterations, deletions and duplications smaller than 100 Kb, thus significantly increasing the detection rate of chromosomal aberrations in patients suffering from idiopathic mental retardation. Several centers are already using array technology in their routine setting in the diagnostic approach to syndromes. Therefore, this overview focuses on the similarities, as well as the differences, of several basic array techniques.     

Interpretation of array comparative genome hybridization data: a major challenge

The advent and application of high-resolution array-based comparative genome hybridization (array CGH) has led to the detection of large numbers of copy number variants (CNVs) in patients with developmental delay and/or multiple congenital anomalies as well as in healthy individuals. The notion that CNVs are also abundantly present in the normal population challenges the interpretation of the clinical significance of detected CNVs in patients. In this review we will illustrate a general clinical workflow based on our own experience that can be used in routine diagnostics for the interpretation of CNVs.     

Medical applications of array CGH and the transformation of clinical cytogenetics

Microarray-based comparative genomic hybridization (array CGH) merges molecular diagnostics with traditional chromosome analysis and is transforming the field of cytogenetics. Prospective studies of individuals with developmental delay and dysmorphic features have demonstrated that array CGH has the ability to detect any genomic imbalance including deletions, duplications, aneuploidies and amplifications. Detection rates for chromosome abnormalities with array CGH range from 5-17% in individuals with normal results from prior routine cytogenetic testing. In addition, copy number variants (CNVs) were identified in all studies. These CNVs may include large-scale variation and can confound the diagnostic interpretations. Although cytogeneticists will require additional training and laboratories must become appropriately equipped, array CGH holds the promise of being the initial diagnostic tool in the identification of visible and submicroscopic chromosome abnormalities in mental retardation and other developmental disabilities.     

High resolution array in the clinical approach to chromosomal phenotypes

Array genomic hybridization (AGH) has recently been implemented as a diagnostic tool for the detection of submicroscopic copy number variants (CNVs) in patients with developmental disorders. However, there is no consensus regarding the choice of the platform, the minimal resolution needed and systematic interpretation of CNVs. We report our experience in the clinical diagnostic use of high resolution AGH up to 100 kb on 131 patients with chromosomal phenotypes but previously normal karyotype. We evaluated the usefulness in our clinics and laboratories by the detection rate of causal CNVs and CNVs of unknown clinical significance and to what extent their interpretation would challenge the systematic use of high-resolution arrays in clinical application. Prioritizing phenotype-genotype correlation in our interpretation strategy to criteria previously described, we identified 33 (25.2%) potentially pathogenic aberrations. 16 aberrations were confirmed pathogenic (16.4% syndromic, 8.5% non-syndromic patients); 9 were new and individual aberrations, 3 of them were pathogenic although inherited and one is as small as approx 200 kb. 13 of 16 further CNVs of unknown significance were classified likely benign, for 3 the significance remained unclear. High resolution array allows the detection of up to 12.2% of pathogenic aberrations in a diagnostic clinical setting. Although the majority of aberrations are larger, the detection of small causal aberrations may be relevant for family counseling. The number of remaining unclear CNVs is limited. Careful phenotype-genotype correlations of the individual CNVs and clinical features are challenging but remain a hallmark for CNV interpretation.     

Detection of copy number variations and their effects in Chinese bulls

Background

Copy number variations (CNVs) are a main source of genomic structural variations underlying animal evolution and production traits. Here, with one pure-blooded Angus bull as reference, we describe a genome-wide analysis of CNVs based on comparative genomic hybridization arrays in 29 Chinese domesticated bulls and examined their effects on gene expression and cattle growth traits.

Results

We identified 486 copy number variable regions (CNVRs), covering 2.45% of the bovine genome, in 24 taurine (Bos taurus), together with 161 ones in 2 yaks (Bos grunniens) and 163 ones in 3 buffaloes (Bubalus bubalis). Totally, we discovered 605 integrated CNVRs, with more “loss” events than both “gain” and “both” ones, and clearly clustered them into three cattle groups. Interestingly, we confirmed their uneven distributions across chromosomes, and the differences of mitochondrion DNA copy number (gain: taurine, loss: yak & buffalo). Furthermore, we confirmed approximately 41.8% (253/605) and 70.6% (427/605) CNVRs span cattle genes and quantitative trait loci (QTLs), respectively. Finally, we confirmed 6 CNVRs in 9 chosen ones by using quantitative PCR, and further demonstrated that CNVR22 had significantly negative effects on expression of PLA2G2D gene, and both CNVR22 and CNVR310 were associated with body measurements in Chinese cattle, suggesting their key effects on gene expression and cattle traits.

Conclusions

The results advanced our understanding of CNV as an important genomic structural variation in taurine, yak and buffalo. This study provides a highly valuable resource for Chinese cattle’s evolution and breeding researches.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-480) contains supplementary material, which is available to authorized users.     

Genomic copy number variation in Mus musculus

Background

Copy number variation is an important dimension of genetic diversity and has implications in development and disease. As an important model organism, the mouse is a prime candidate for copy number variant (CNV) characterization, but this has yet to be completed for a large sample size. Here we report CNV analysis of publicly available, high-density microarray data files for 351 mouse tail samples, including 290 mice that had not been characterized for CNVs previously.

Results

We found 9634 putative autosomal CNVs across the samples affecting 6.87 % of the mouse reference genome. We find significant differences in the degree of CNV uniqueness (single sample occurrence) and the nature of CNV-gene overlap between wild-caught mice and classical laboratory strains. CNV-gene overlap was associated with lipid metabolism, pheromone response and olfaction compared to immunity, carbohydrate metabolism and amino-acid metabolism for wild-caught mice and classical laboratory strains, respectively. Using two subspecies of wild-caught Mus musculus, we identified putative CNVs unique to those subspecies and show this diversity is better captured by wild-derived laboratory strains than by the classical laboratory strains. A total of 9 genic copy number variable regions (CNVRs) were selected for experimental confirmation by droplet digital PCR (ddPCR).

Conclusion

The analysis we present is a comprehensive, genome-wide analysis of CNVs in Mus musculus, which increases the number of known variants in the species and will accelerate the identification of novel variants in future studies.

Electronic supplementary material

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

Chromosome 15q13.3 microduplications are associated with treatment refractory major depressive disorder

Major depressive disorder (MDD) affects approximately 15 million Americans. Approximately 2 million of these are classified as being refractory to treatment (TR‐MDD). Because of the lack of available therapies for TR‐MDD, and the high risk of suicide, there is interest in identifying new treatment modalities and diagnostic methods. Understanding of the impact of genomic copy number variation in the etiology of a variety of neuropsychiatric phenotypes is increasing. Low copy repeat elements at 15q13.3 facilitate non‐allelic homologous recombination, resulting in recurrent copy number variants (CNVs). Numerous reports have described association between microdeletions in this region and a variety of neuropsychiatric phenotypes, with CHRNA7 implicated as a candidate gene. However, the pathogenicity of 15q13.3 duplications is less clear. As part of an ongoing study, in which we have identified a number of metabolomic anomalies in spinal fluid from TR‐MDD patients, we also evaluated genomic copy number variation in patients (n = 125) and controls (n = 26) via array‐based copy number genomic hybridization (CGH); the case frequency was compared with frequencies reported in a prior study as well as a larger population‐sized cohort. We identified five TR‐MDD patients with microduplications involving CHRNA7. CHRNA7 duplications are the most common CNVs identified by clinical CGH in this cohort. Therefore, this study provides insight into the potential involvement of CHRNA7 duplications in the etiology of TR‐MDD and informs those involved with care of affected individuals.     

Allele-specific disparity in breast cancer

Background

Molecular alterations critical to development of cancer include mutations, copy number alterations (amplifications and deletions) as well as genomic rearrangements resulting in gene fusions. Massively parallel next generation sequencing, which enables the discovery of such changes, uses considerable quantities of genomic DNA (> 5 ug), a serious limitation in ever smaller clinical samples. However, a commonly available microarray platforms such as array comparative genomic hybridization (array CGH) allows the characterization of gene copy number at a single gene resolution using much smaller amounts of genomic DNA. In this study we evaluate the sensitivity of ultra-dense array CGH platforms developed by Agilent, especially that of the 1 million probe array (1 M array), and their application when whole genome amplification is required because of limited sample quantities.

Methods

We performed array CGH on whole genome amplified and not amplified genomic DNA from MCF-7 breast cancer cells, using 244 K and 1 M Agilent arrays. The ADM-2 algorithm was used to identify micro-copy number alterations that measured less than 1 Mb in genomic length.

Results

DNA from MCF-7 breast cancer cells was analyzed for micro-copy number alterations, defined as measuring less than 1 Mb in genomic length. The 4-fold extra resolution of the 1 M array platform relative to the less dense 244 K array platform, led to the improved detection of copy number variations (CNVs) and micro-CNAs. The identification of intra-genic breakpoints in areas of DNA copy number gain signaled the possible presence of gene fusion events. However, the ultra-dense platforms, especially the densest 1 M array, detect artifacts inherent to whole genome amplification and should be used only with non-amplified DNA samples.

Conclusions

This is a first report using 1 M array CGH for the discovery of cancer genes and biomarkers. We show the remarkable capacity of this technology to discover CNVs, micro-copy number alterations and even gene fusions. However, these platforms require excellent genomic DNA quality and do not tolerate relatively small imperfections related to the whole genome amplification.     

Identification of copy number variation of CAPN10 in Thais with type 2 diabetes by multiplex PCR and denaturing high performance liquid chromatography (DHPLC)

Copy number variations (CNVs) have been shown to be associated with several diseases. They can cause deviation of genotypes from Hardy-Weinberg Equilibrium (HWE). Genetic case-control association studies in Thais revealed that genotype distribution of CAPN10 Indel19 was deviated from HWE after correction of genotyping error. Therefore, we aim to identify CNVs within CAPN10 Indel19 region. The semi-quantitative denaturating high performance liquid chromatography (DHPLC) method was used to detect CNVs in the region of CAPN10 Indel19 marker in cohort of 305 patients with type 2 diabetes and 250 control subjects without diabetes. CNVs in the region of CAPN10 Indel19 was successfully detected by DHPLC. After correction of genotype calling based on the status of identified CNVs, CAPN10 Indel19 genotypes were well-fitted for HWE (p>0.05). However, we did not find association between CNV genotypes and risk of type 2 diabetes in our population. CNVs in CAPN10 have been identified in Thais. These CNVs lead to deviation from HWE of CAPN10 Indel19 genotypes. After excluding identified CNVs from the analysis, CAPN10 Indel19 was associated with type 2 diabetes. The information obtained from our study would be helpful for genotyping accuracies of SNPs residing in the CNVs region.     

Stability of ultramer as copy number standards in real-time PCR

Since commercial copy number standards are not always available for real-time PCR, alternative sources of DNA are used. Unfortunately, stored genomic DNA or PCR amplicon has been shown to be unstable, resulting in variable copy number. More recently, the use of ultramer as copy number standard has been reported. However, there is little information on the stability of ultramer under different storage conditions. Thus the aim of this study was to determine the stability of ultramer as copy number standard under different storage conditions using different mixing methods. We found that ultramer copy number was not affected by storage at either 4 °C or − 20 °C over a period of 30 days. Furthermore, the method of mixing the ultramer did not appear to contribute to variability in results. Irrespective of storage temperature or mixing method, there was less than 5% variance in Ct value over a period of 30 days. A duplicate set of standards costs approximately $0.01. Therefore, the use of ultramer as copy number standards in real-time PCR, is cost effective and convenient.     

Gene amplification profiling of esophageal squamous cell carcinomas by DNA array CGH

Gene amplification is one of the basic mechanisms that lead to overexpression of oncogenes. DNA array comparative genomic hybridization (CGH) has great potential for comprehensive analysis of both a relative gene-copy number and altered chromosomal regions in cancers, which enables us to identify new amplified genes and unstable chromosomal loci. We examined the amplification status in 32 esophageal squamous cell carcinomas (ESCCs) and 13 ESCC cell lines on 51 frequently amplified loci in a variety of cancers by both DNA array CGH and Southern blot analyses. The 1p34 locus containing MYCL1, 2p24 (MYCN), 7p12 (EGFR), and 12q14 (MDM2) were amplified in one of the 32 cases (3%), and the 17q12 locus (ERBB2) and 8p11 (FGFR1) in two of the 32 cases (6%), while only the 11q13 locus (Cyclin D1, FGF4, and EMS1) was frequently amplified (28%, 9/32), demonstrating this locus to be a major target in ESCCs. One locus, 8q24 (c-MYC) was found to be amplified only in the cell lines. Eight out of 51 loci (15.7%) were found to be amplified in at least one of the 32 primary ESCCs or the 13 ESCC cell lines, suggesting that chromosomal loci frequently amplified in a type of human cancer may also be amplified in other types of cancers. This paper is the first report of an application of DNA array CGH to ESCCs.     

An idic(15) associated with POF (premature ovarian failure): molecular cytogenetic definition of a case and review of the literature

We report on a 36-year-old infertile woman, presenting a premature ovarian failure with an otherwise normal female phenotype. Cytogenetic analyses showed the presence of a supernumerary marker chromosome, that was characterized by FISH (fluorescent in situ hybridization) and array CGH (comparative genomic hybridization). This marker chromosome was derived from chromosome 15, and contained only heterochromatic material. The Prader Willi/Angelman region was not present. No duplications of the 15q regions were detected by array CGH. Supernumerary markers of chromosome 15 have been reported in cases of infertility and amenorrhea, that is also described in cases with marker derived by other acrocentric chromosomes. The case here presented constitutes a further example that etiology of POF is not always associated with a defective gene, but in some cases oocytes atresia can be the consequence of the abnormal meiotic pairing of chromosomes.     

Spectrum of Cytogenomic Abnormalities Revealed by Array Comparative Genomic Hybridization on Products of Conception Culture Failure and Normal Karyotype Samples

Approximately 30% of pregnancies after implantation end up in spontaneous abortions, and 50% of them are caused by chromosomal abnormalities. However, the spectrum of genomic copy number variants(CNVs) in products of conception(POC) and the underlying genedosage-sensitive mechanisms causing spontaneous abortions remain largely unknown. In this study, array comparative genomic hybridization(a CGH) analysis was performed as a salvage procedure for 128 POC culture failure(POC-CF) samples and as a supplemental procedure for106 POC normal karyotype(POC-NK) samples. Chromosomal abnormalities were detected in 10% of POC-CF and pathogenic CNVs were detected in 3.9% of POC-CF and 5.7% of POC-NK samples. Compiled results from this study and relevant case series through a literature review demonstrated an abnormality detection rate(ADR) of 35% for chromosomal abnormalities in POC-CF samples, 3.7% for pathogenic CNVs in POC-CF samples, and 4.6% for pathogenic CNVs in POC-NK samples. Ingenuity Pathway Analysis(IPA) was performed on the genes from pathogenic CNVs found in POC samples. The denoted primary gene networks suggested that apoptosis and cell proliferation pathways are involved in miscarriage. In summary, a similar spectrum of cytogenomic abnormalities was observed in POC culture success and POC-CF samples. A threshold effect correlating the number of dosage-sensitive genes in a chromosome with the observed frequency of autosomal trisomy is proposed. A rationalized approach using firstly fluorescence in situ hybridization(FISH) testing with probes of chromosomes X/Y/18, 13/21, and 15/16/22 for common aneuploidies and polyploidies and secondly a CGH for other cytogenomic abnormalities is recommended for POC-CF samples.     

Considering specific clinical features as evidence of pathogenic copy number variants

Since the introduction of high-resolution microarray technologies, it has become apparent that structural chromosomal rearrangements can lead to a wide variety of clinical manifestations, including developmental delay/intellectual disability (DD/ID). It has been shown previously that the diagnostic yield of genome-wide array-based identification of submicroscopic alterations in patients with ID varies widely and depends on the patient selection criteria. More attempts have recently been made to define the phenotypic clues of pathogenic copy number variants (CNVs). The aim of this study was to investigate a well-phenotyped cohort of patients with DD/ID and determine whether certain clinical features may serve as indicators for pathogenic CNVs. A retrospective analysis was conducted for patients with DD/ID (n?=?211) who were tested using genome-wide chromosomal microarray technologies and a review of the clinical data was performed. Pathogenic CNVs were detected in 29 patients. In comparison with individuals who had normal molecular karyotyping results (n?=?182), malformations of the musculoskeletal system; congenital malformations of the CNS (particularly hydrocephalus and congenital malformations of the corpus callosum); minor anomalies of the eye, face, and neck subgroup (particularly downward-slanting palpebral fissures, minor anomalies of the ear, and micrognathia); brachydactyly; and umbilical hernia were more common in patients with chromosomal alterations. A multivariate logistic regression analysis allowed the identification of three independent pathogenic CNV predictors: congenital malformations of the corpus callosum, minor anomalies of the ear, and brachydactyly. Insights into the chromosomal phenotype may help to increase the diagnostic yield of microarray technologies and sharpen the distinction between chromosomal alterations and other conditions.     

Novel URAT1 mutations caused acute renal failure after exercise in two Chinese families with renal hypouricemia

Renal hypouricemia (RHUC), as an infrequent hereditary disease, is associated with severe complications such as exercise-induced acute renal failure (EIARF). Loss-of-function mutations in urate transporter gene URAT1 (Type 1) and in glucose transporter gene GLUT9 (Type 2) are major causes of this disorder. In this study, URAT1 and GLUT9 were screened in two uncorrelated families from mainland China and a total of five mutations were identified in exons, including two novel heterozygous URAT1 mutations. In four members of the first family, c.151delG (p.A51fsX64) in exon 1 was detected, which resulted in a frameshift and truncated the original 553-residue-protein to 63 amino acid protein. A missense mutation c.C1546A (p.P516T) in exon 9 in GLUT9 was revealed in the second family, which caused a functional protein substitution at codon 516. These two novel mutations were neither identified in the subsequent scanning of 200 ethnically matched healthy control subjects with normal serum UA level nor in a 1000 genome project database. Thus our report identifies two novel loss-of-function mutations (c.151delG in URAT1 and p.P516T in GLUT9) which cause RHUC and renal dysfunction in two independent RHUC pedigrees.