|Home | About | Journals | Submit | Contact Us | Français|
Array genomic hybridization is being used clinically to detect pathogenic copy number variants in children with intellectual disability and other birth defects. However, there is no agreement regarding the kind of array, the distribution of probes across the genome, or the resolution that is most appropriate for clinical use.
We performed 500 K Affymetrix GeneChip® array genomic hybridization in 100 idiopathic intellectual disability trios, each comprised of a child with intellectual disability of unknown cause and both unaffected parents. We found pathogenic genomic imbalance in 16 of these 100 individuals with idiopathic intellectual disability. In comparison, we had found pathogenic genomic imbalance in 11 of 100 children with idiopathic intellectual disability in a previous cohort who had been studied by 100 K GeneChip® array genomic hybridization. Among 54 intellectual disability trios selected from the previous cohort who were re-tested with 500 K GeneChip® array genomic hybridization, we identified all 10 previously-detected pathogenic genomic alterations and at least one additional pathogenic copy number variant that had not been detected with 100 K GeneChip® array genomic hybridization. Many benign copy number variants, including one that was de novo, were also detected with 500 K array genomic hybridization, but it was possible to distinguish the benign and pathogenic copy number variants with confidence in all but 3 (1.9%) of the 154 intellectual disability trios studied.
Affymetrix GeneChip® 500 K array genomic hybridization detected pathogenic genomic imbalance in 10 of 10 patients with idiopathic developmental disability in whom 100 K GeneChip® array genomic hybridization had found genomic imbalance, 1 of 44 patients in whom 100 K GeneChip® array genomic hybridization had found no abnormality, and 16 of 100 patients who had not previously been tested. Effective clinical interpretation of these studies requires considerable skill and experience.
Chromosomal imbalance has been recognized as the most frequent cause of intellectual disability (ID) for 50 years [1-3]. Until recently, most of this genomic imbalance was diagnosed by cytogenetic analysis, but studies over the past few years have found that ID is caused by constitutional gains or losses of submicroscopic genomic segments even more often than by microscopically-apparent chromosomal aberrations [4-6]. Our ability to recognize these submicroscopic genomic changes, which are usually called copy number variants (CNVs), as the most frequent cause of ID depends on the use of Array Genomic Hybridization (AGH) (also known as array-comparative genomic hybridization, chromosomal microarray analysis, or copy number analysis). AGH can survey the entire genome for imbalance that is 1/100th the size of that detectable by conventional cytogenetic analysis.
Although AGH is now being used routinely as a clinical test for the identification of chromosomal imbalance in people with ID and other birth defects, controversy still exists regarding the most appropriate platform to use for this purpose. Initial clinical studies were done with arrays having a few thousand BACs distributed at 1-3 Mb intervals across the genome or with BACs targeted to regions where pathogenic submicroscopic deletions or duplications were known to occur. More recent studies have shown that arrays with higher resolution and genome-wide coverage provide better detection rates for pathogenic CNVs in children with ID and normal cytogenetic analysis [7-12]. Other methods have identified pathogenic CNVs that are too small to have been detected by the array platforms used in most AGH studies [13-15], so analysis at even higher resolution may be necessary to detect all pathogenic CNVs in children with ID.
Unfortunately, use of higher resolution AGH for detection of genomic imbalance is confounded because most CNVs are not pathogenic. Estimates of the mean number of apparently benign CNVs per person range from 6-824 in various studies, depending on the technology used to identify the variants and the size range used to define a CNV [16-25]. Sequencing of the complete diploid genomes of normal individuals has shown that the number of CNVs per person is actually even greater [26-29]. Distinguishing these benign CNVs from those that cause ID and other birth defects is the most serious challenge to routine clinical use of AGH.
We previously reported our findings using 100 K Affymetrix GeneChip® AGH [10,16] to perform a genome-wide survey of benign and pathogenic CNVs in 100 idiopathic ID trios, each comprised of an affected child and both unaffected parents. Here we describe the results of a study of 100 new idiopathic ID trios, as well as 54 of the trios tested previously, using 500 K Affymetrix GeneChip® AGH. We found that higher resolution AGH detected a larger number of apparently pathogenic CNVs in both groups. Many benign CNVs, including at least one that was de novo, were also detected with the 500 K AGH, but it was possible to distinguish benign and pathogenic CNVs with confidence in almost all cases.
We performed AGH with 500 K Affymetrix GeneChip® arrays on 154 children with idiopathic ID and on both parents of each affected child. Fifty-four of these trios (called "the 100 K cohort") had previously been studied with lower-resolution 100 K GeneChip® AGH [10,16]; the other 100 trios (called "the new cohort") were studied by AGH for the first time with the 500 K GeneChip® arrays. Data were analyzed to determine copy number along the length of all chromosomes except the Y (for which there are no probes on either the 100 K or 500 K GeneChip® arrays), and the findings for each child were compared with those for his or her parents. Autosomal CNVs seen in the child and in at least one parent were considered likely to be benign polymorphisms. Autosomal CNVs found in the child but not in either parent were evaluated by an independent method to confirm the presence of the CNV and its de novo occurrence. CNV calls on the X-chromosome in a female child were validated by an independent method if they appeared to have occurred de novo; CNV calls on the X-chromosome in a male child were validated by an independent method whether they appeared to be de novo or to have been inherited from the mother.
We found a total of 4577 hits (putative CNVs including at least 10 contiguous SNPs called by the SMD software with a p-value below 1 × 10-8) in the 462 samples (154 trios) analyzed by 500 K GeneChip® AGH. This is an average of about 10 hits per sample, which is an underestimate of the total number of benign CNVs present because of the stringent cutoff used to obtain a false discovery rate of less than 5% (see Methods).
Within the 54 trios who were studied with both platforms, we found four times as many hits with 500 K AGH as we did with 100 K AGH (Table (Table1).1). The ratio of the number of putative CNVs called with the 500 K platform to the number called with the 100 K platform was 11.0 for CNVs between 100 kb and 200 kb in size but was lower for both smaller and larger CNVs (4.1 for those < 100 kb and 3.4 for those between 200 kb and 500 kb).
The 4577 putative CNVs called in all 154 trios were subjected to further bioinformatic analysis (see Methods) to produce a final annotated list of 58 apparently de novo CNVs and two cases of mosaic trisomy that were called in 50 patients. These apparent genomic imbalances were subjected to validation by independent methods. Thus, an apparent de novo CNV call that required independent validation was made in about 1 child in 3. Thirty-three of these CNVs in 30 patients and both cases of mosaic trisomy were confirmed to be de novo by an independent method and are described in detail below. The other putative CNVs were found to be present in both of the parents as well as in the child (false negative AGH calls in the parents) in two instances or could not be confirmed to be present in the child (false positive AGH calls) in 22 instances. Altogether, false positive CNV calls were made in 21 (13.6%) of the 154 trios studied and false negative CNV calls were made in 2 (1.3%) of the 154 trios studied. In one other trio (Family 5202), an apparent de novo deletion of chromosome 14 called on AGH in the child was found by FISH to be a duplication of the region in both parents instead.
Nineteen of 100 children with ID in the new cohort were found by Affymetrix 500 K GeneChip® AGH to have de novo genomic imbalance that was confirmed by FISH, MLPA, AGH on an Agilent® 244 K platform or cytogenetic re-analysis (Table (Table2).2). One of these children (Patient 8056) had mosaic trisomy 9, and two were found to have de novo unbalanced reciprocal translocations - a der(10)t(2;10)(q37;q26.13) in Patient 873 and a der(4)t(4;8)(p16.1;p23.1) in Patient 5814 - each producing both a terminal duplication and a terminal deletion identified by AGH. We found de novo submicroscopic deletions in 13 other patients and de novo submicroscopic duplications in three other patients. The deletions ranged in size from 89 kb to 11.0 Mb; six were less than 1 Mb. The duplications ranged in size from 362 kb to 11.1 Mb; one was less than 1 Mb.
Using 500 K Affymetrix GeneChip® AGH, we also confirmed the genomic imbalance that had previously been identified in 10 of the 54 ID trios from the 100 K cohort (Table (Table33 and Additional File 1: Supplemental Table S1). In addition, we identified and confirmed by FISH two de novo CNVs that were not called on the 100 K assay - a 1.6 Mb duplication of 8q23.2-23.3 in Patient 3890 and a 1.5 Mb deletion of 4p16.3 in Patient 4840 (Table (Table33 and Figure Figure11).
We found two instances of uniparental disomy (UPD), diagnosed by the occurrence of mendelian inconsistency in a region with a normal copy number of 2 , among the 100 ID trios in the new cohort studied by 500 K Affymetrix GeneChip® AGH (Table (Table44 and Figures Figures2A2A and and2B).2B). Patient 6904 has mosaic paternal isodisomy of most of the short arm of chromosome 11. Patient 1658 has maternal UPD 16, being heterodisomic for the central portion of the chromosome and isodisomic for both ends. Both cases were confirmed to be disomic with microsatellite markers (Figures (Figures2C2C and and2D2D).
We judged the mosaic trisomy 9, both unbalanced reciprocal translocations, 11 of the other de novo deletions, and two of the other de novo duplications found in the new cohort to be pathogenic (Table (Table2).2). A mosaic 107 kb de novo deletion of chromosome 14 q11.2 (Patient 818) and a homozygous 89 kb deletion of the HLA-G region that resulted from transmission by parents who both carried heterozygous deletions of the same region (Patient 216) were judged to be benign variants. A 186 kb deletion of chromosome 21q22.11 (Patient 8619), a 362 kb duplication of chromosome 22 q11.21 (Patient 9979) and both cases of UPD (Patients 1658 and 6904) were of uncertain clinical significance.
Of the two de novo CNVs identified by 500 K Affymetrix GeneChip® AGH in the 44 children with ID whose 100 K GeneChip® AGH studies had been interpreted as normal (Table (Table3),3), one is likely to be pathogenic and the other is of uncertain clinical significance. A 1.5 Mb deletion of 4p16.3 in Patient 4840 is pathogenic because of its size, the inclusion of two genes that have been implicated in the Wolf-Hirschhorn syndrome (WHSC1 and WHSC2), overlap with known pathogenic CNVs, and a compatible clinical phenotype. The 1.6 Mb duplication of 8q23.2-23.3 in Patient 3890 is of unknown clinical significance. It includes no RefSeq genes but involves a large (1.6 Mb) genomic region, most of which has never been reported to be polymorphic in normal individuals.
Because its detection rate for pathogenic genomic imbalance is much higher than that of conventional cytogenetic analysis, a consensus has developed that AGH should be used clinically for the evaluation of patients with ID and other birth defects [31-38]. It is clear that AGH using "targeted" arrays that only include probes for genomic regions known to be involved in microdeletion or microduplication syndromes has substantially lower detection rates for CNVs that cause ID than AGH using arrays that provide genome-wide coverage [37,39-41]. Beyond this, however, there is no agreement regarding the kind of array, the distribution of probes across the genome, or the resolution that is most appropriate for clinical use. Although BAC arrays were initially used, most clinical laboratories now prefer oligonucleotide arrays because high-quality platforms that produce consistent results are reliably available from commercial sources. In addition, the use of larger numbers of smaller probes on oligonucleotide arrays permits more precise delineation of the breakpoints of CNVs that are detected, which facilitates genotype-phenotype correlation and clinical interpretation. AGH with a SNP array provides the additional advantage of generating genotypes that can be used to verify family relationships and find uniparental disomy as well as a second method (in addition to hybridization intensity) for identifying genomic imbalance [10,42-44].
We previously reported that 100 K Affymetrix GeneChip® AGH is a robust platform for the detection of pathogenic CNVs in patients with ID . Here we show that the detection rate of CNVs among such patients is higher with 500 K GeneChip® AGH than with 100 K GeneChip® AGH. We made about four times as many CNV calls overall with the 500 K platform as with the 100 K platform when using the same method of bioinformatic analysis in 54 trios studied with both technologies (Table (Table1).1). We also found 18 instances of pathogenic genomic imbalance in 16 of 100 children with ID and normal cytogenetic analysis studied by 500 K GeneChip® AGH (Table (Table2),2), compared to 11 instances of pathogenic genomic imbalance in 11 of 100 similarly-ascertained children tested by 100 K GeneChip® AGH in our previous study. Although the higher detection rate we observed with the 500 K platform may have occurred by chance because we just happened to include a few more patients with such genomic changes in the new cohort than in the 100 K cohort, we also found two additional de novo CNVs by 500 K GeneChip® AGH among the 44 children whose 100 K AGH was interpreted as normal in our earlier studies (Figure (Figure11 and Table Table33).
We detected three apparently de novo CNVs smaller than 200 kb among the 100 trios tested with 500 K GeneChip® AGH in the new cohort (a 107 kb deletion of chromosome 14 in Patient 818, a 186 kb deletion of chromosome 21 in Patient 8619, and an 89 kb deletion of chromosome 6 in Patient 216 that was actually a homozygous loss inherited from two heterozygous parents), but none of these CNVs was clearly pathogenic. The overall size distribution of pathogenic CNVs detected by 500 K GeneChip® AGH among 154 children with idiopathic ID in the present study is similar to that observed by 100 K GeneChip® AGH among 100 children with idiopathic ID whom we studied previously  (see Additional File 2: Supplemental Figure S1). The higher detection rate on the 500 K array therefore appears to be related more to better probe coverage in relevant genomic regions and an improved ability to distinguish CNVs from background noise, rather than to a capacity to identify much smaller pathogenic CNVs. This is illustrated in Families 3890 and 4840 (Figure (Figure1),1), in which a 1.6 Mb duplication of 8q23.2q23.3 and a 1.5 Mb deletion of 4p16.3, respectively, are obvious on the 500 K AGH but were not called on the 100 K analysis. In retrospect, the 4p16.3 deletion in Patient 4840 can be seen on the 100 K AGH copy number plot despite the noisy data, but it was not called by either the automated analysis or visual inspection of these plots when the initial study was done. Our failure to detect the de novo duplication of 8q23.2q23.3 in Patient 3890 was probably caused by the noisy data in the father's study.
Distinguishing benign CNVs from those that cause ID and other birth defects is a critical issue in routine clinical use of AGH. Benign CNVs occur in all people and are a major source of genetic variation in the normal population [21,27,29]. Most apparently benign CNVs over 2 kb in size occur as polymorphisms with minor allele frequencies of at least 5%  and are inherited from a parent [21,23,45,46].
Benign and pathogenic CNVs can usually be distinguished in patients with ID and other birth defects by inheritance and genotype-phenotype correlation [5,33,47]. In this study, we identified a mean of about 10 CNVs per subject in the 154 ID trios tested by 500 K AGH. The vast majority of these CNVs were characterized as benign because they were inherited from a normal parent. Genomic imbalance that occurs de novo in a patient with ID whose parents are normal is more likely to be pathogenic than genomic imbalance that was inherited unchanged from a normal parent. We performed AGH on both parents of every child with ID to determine the inheritance of the CNVs found in the child, but this is sometimes not possible in clinical practice. In such instances it is necessary to infer likely de novo occurrence by information obtained from populations that have previously been studied [19,22-24,48,49]. Great care must be taken to avoid misinterpretation when this is done, especially if the available data were obtained with lower resolution AGH, the phenotypic characteristics of the comparison population are uncertain, reported polymorphic CNVs do not have exactly the same breakpoints as the CNV of interest, or the population frequency of a previously-reported CNV is unknown.
Compelling evidence that a CNV in a person with ID is pathogenic exists if the genomic imbalance is known to cause the patient's phenotype in other individuals, e.g., if a child with del 9q34.3 has features of the 9q subtelomeric deletion syndrome (Patient 523), a child with del 1p36.32p36.33 has features of the 1p36 deletion syndrome (Patient 9133), or a child with del 17q21.31 has features of the syndrome associated with this deletion (Patient 2106). Pathogenicity is also supported when a CNV includes a gene that is known to cause the patient's phenotype when inactivated (if the CNV is a deletion) or over-expressed (if the CNV is a duplication). On the other hand, a CNV is unlikely to be pathogenic if it involves a highly polymorphic region in which genomic loss (or gain, whichever is present in the patient) of the entire involved segment is known to occur in normal people.
If a direct genotype-phenotype correlation of this kind cannot be made in a particular case, certain genetic features of the CNV may provide clues to its pathogenicity. CNVs that are larger and those that involve gene-rich regions are more likely to be pathogenic than CNVs that are smaller and involve only gene-poor regions [5,47]. In addition, clinical experience suggests that deletions are more likely to be pathogenic than duplications . The genetic content of a CNV may also make pathogenicity more or less likely. For example, involvement of a gene that lies within a pathway that is known to contain other dosage-sensitive genes associated with a similar phenotype strengthens the possibility of pathogenicity, while a CNV that does not contain any genes that are expressed in relevant tissues during embryogenesis is unlikely to be pathogenic for ID.
There are, of course, exceptions to each of these "rules". Some benign CNVs arise de novo [21,23,45,46,50], as appears to have occurred in the de novo deletion of chromosome 14q11.2 we found in Patient 818. The 107 kb region involved is highly polymorphic and contains several T-cell receptor variable region genes. On the other hand, some CNVs that are inherited from a normal parent are pathogenic for ID. Examples include maternal transmission of a UBE3A deletion to a child with Angelman syndrome , maternal transmission of a MECP2 duplication to a son , and CNVs such as dup 22q11.2 [53,54] or del 1q21.2  that can cause ID but exhibit incomplete penetrance.
Although large (> 250 kb) CNVs are often pathogenic, they may be benign [19,29]. Most benign CNVs are small (< 250 kb) [19,27,29], and it seems probable that the smaller a CNV, the more likely it is to be benign. Nevertheless, no clear size distinction exists between benign and pathogenic CNVs. We found pathogenic CNVs as small as the 298 kb deletion of 9q34.3 in Patient 523 in this study, and others have reported even smaller pathogenic CNVs [13-15,56-58]. Expression patterns, functional annotation and animal models can provide important clues to pathogenicity in some cases, but without knowledge of the phenotypic effects of a copy number alteration in humans, one can rarely, if ever, be certain whether a novel gain or loss of a particular genomic region can produce ID or other birth defects.
In our 500 K AGH study of 154 ID trios, 58 de novo CNVs were called by bioinformatic analysis, and 33 of these CNVs were confirmed and shown to be de novo by an independent method. Because we could assess the phenotypes of our patients in detail and correlate the findings with those obtained by AGH, we were able to determine with confidence whether the genomic imbalance we observed was pathogenic or not in every case studied except three (Tables (Tables22 and and3).3). Such genotype-phenotype correlation is critical to determining the effects of novel CNVs detected by AGH in patients with ID.
A CNV of uncertain clinical significance was encountered in three (1.9%) of the 154 trios analyzed in this series - a 362 kb duplication of 22q11.21 in Patient 9979 (Table (Table2),2), a 186 kb deletion of 21q22.11 in Patient 8619 (Table (Table2)2) and a 1.6 Mb duplication of chromosome 8q23.2q23.3 in Patient 3890 (Table (Table3).3). This rate of CNVs of uncertain clinical significance is similar to that reported in large series of patients with ID and other birth defects studied by AGH with "targeted" chips [31,35,59].
We were uncertain of the clinical significance of either case of UPD that we detected. Although only a few liveborn children with UPD 16 have been recognized, the reported experience does not suggest that UPD 16 can cause the abnormalities observed in Patient 1658 [60,61]. Paternal UPD 11p15 can produce Beckwith-Wiedemann syndrome , but this phenotype is very different from that observed in the affected child in Family 6904. However, as both of these cases involved isodisomy of a portion of the chromosome, we cannot rule out the possibility that the abnormal phenotype was produced by homozygosity for a recessive mutant allele [63,64]. Although the detection of UPD in addition to alterations of copy number is a theoretical benefit of using an array that includes probes for SNPs, the clinical utility of genome-wide screening for UPD in patients with idiopathic ID and other birth defects is uncertain.
We detected more pathogenic CNVs with 500 K AGH than with 100 K AGH, but some CNVs that were present among our patients were not detected using the 500 K assay. For example, our 500 K GeneChip® analysis failed to identify a pathogenic 83 kb deletion of chromosome 16p13.3 (3,862,993 bp to 3,945,522 bp) involving the CREBBP gene (Patient 5121). This de novo deletion was found by AGH on the Agilent® 244 K platform and was confirmed by MLPA. The patient is an 8 year-old boy whose clinical features are characteristic of the Rubinstein-Taybi syndrome, which has been associated with deletions and other mutations of CREBBP in other patients [65,66]. The 83 kb genomic region deleted in our patient is poorly represented on the 500 K GeneChip® arrays, with a total of only 15 SNPs. SNP arrays have uneven genomic coverage, and the addition of non-polymorphic oligonucleotide probes to the design of arrays like the one used in this study has been shown to provide substantially better detection of CNVs .
Affymetrix GeneChip® 500 K array genomic hybridization performed in individuals with idiopathic intellectual disability detected pathogenic genomic imbalance in 10 of 10 patients in whom 100 K GeneChip® array genomic hybridization found genomic imbalance, 1 of 44 patients in whom 100 K GeneChip® array genomic hybridization had found no abnormality, and 16 of 100 patients who had not previously been tested. Further improvements in array design, ongoing improvements in AGH software, and continuing enhancement of resources like DECIPHER https://decipher.sanger.ac.uk/ and the Toronto Database of Genomic Variants http://projects.tcag.ca/variation/ are helping to establish AGH as the primary clinical tool for recognition of genomic imbalance that causes ID and other birth defects [33,34,37,41,67-70]. It seems likely, however, that no perfect AGH platform for detection of pathogenic CNVs may ever exist and that effective clinical interpretation of these studies will continue to require considerable skill and experience [33,50,71].
We studied 100 children with idiopathic ID who had not been studied by AGH before ("the new cohort") and 54 of the idiopathic ID patients whom we had previously tested with 100 K Affymetrix GeneChip® AGH ("the 100 K cohort"). Ten of the 54 patients in the 100 K cohort had previously been found to have pathogenic genomic imbalance; the other 44 patients had previously been reported to have normal 100 K GeneChip® AGH (see Additional File 1: Supplemental Table S1) [10,16]. We also performed 500 K AGH on both unaffected parents of each child.
All of the children were assessed by a clinical geneticist who was unable to determine the cause of the ID despite thorough clinical evaluation and testing that included routine karyotyping with at least 450-band resolution. Subjects were selected for AGH testing because they had ID or developmental delay and at least one of the following additional clinical features: one major malformation, microcephaly, abnormal growth, or multiple minor anomalies. Informed consent was obtained from each family, and assent was also obtained from the child, if possible. The study was approved by the University of British Columbia Clinical Research Ethics Board.
DNA was extracted from whole blood with a Gentra Puregene DNA Purification Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. The DNA was precipitated in 70% alcohol, resuspended in hydration solution, and stored at 4°C.
DNA degradation, labeling, hybridization, and scanning were performed according to the manufacturer's protocols http://www.affymetrix.com using an Affymetrix Fluidics Station 450, Affymetrix Hybridization Oven 640 and Affymetrix GeneChip Scanner 3000 (Affymetrix, Inc., Santa Clara, CA, USA).
Chip-to-chip normalization, standardization to a reference, genotype detection, and copy number estimation on a single SNP basis were performed using the Affymetrix Power Tools (version 1.6.0) software suite http://www.affymetrix.com, as previously described . Estimation of CNV boundary positions was done using Significance of Mean Difference (SMD), a method that we developed. Briefly, the mean of SNP copy number estimates (or log2 ratios) within a CNV was compared to the mean of those on the rest of the chromosome, and the probability that the null hypothesis of Student's t-test was true, i.e., that the means were from the same distribution, was calculated. A search using this statistic was conducted over different CNV lengths to find the position and length (in number of SNPs) that yielded the lowest probability, defining the boundaries of a putative CNV. For each sample, a random data set was produced by shuffling the genomic positions of the data, and an identical search was conducted. The results of this search represent false discoveries due to the random variation of the individual SNP data.
The p-value distribution of apparent CNVs detected by SMD in all 462 samples (both normal parents and the affected child from 154 ID trios) is shown in Figure Figure3.3. We observed that a CNV call with a p-value of 1 × 10-8 usually had a false discovery rate of less than 5%, while a call with a p-value of 1 × 10-7 often had a false discovery rate of 30%, with some variation from sample to sample.
Since SMD compares a putative CNV to the rest of the chromosome, aberrations on the X chromosome are detected equivalently well for males and females. The pseudo-autosomal regions of the X-chromosomes are exceptions, but CNVs there can be detected as a function of the reference set used in the analysis. GeneChip® 500 K arrays do not contain Y-chromosome probes.
The SMD search was performed on each child with each parent as a reference. Since our goal was to find duplications and deletions that were de novo, a criterion for selection of CNVs for further analysis was that the child have the same aberration with each parent as reference. Correct parental relationships were confirmed in all trios by use of the SNP genotyping calls.
Every putative de novo CNV call was evaluated by another analysis conducted using a large reference set of individuals. Early in the study, the 48 sample reference set available from Affymetrix was used, but later a set of 50 mothers from our own data was used. The use of a local reference set significantly reduced noise. The use of the large reference set was required to detect the occurrence of aberrations in both parents that were not inherited by the child. Putative de novo CNVs were further evaluated by visual examination of the copy number plots in comparison to both parents as well as of the plots for the child and both parents in comparison to the large reference set.
Putative de novo CNVs identified by 500 K AGH were validated by fluorescence in situ hybridization (FISH) or multiplex ligation-dependant probe amplification (MLPA), Agilent® 244 K AGH, or repeat cytogenetic analysis. FISH was performed with BAC or fosmid probes selected using the University of California at Santa Cruz Genome Browser  and the May 2006 assembly of the human genome sequence. BAC or fosmid DNA was isolated by small-scale (miniprep) preparation and was labeled with Spectrum Red or Green (Vysis, Abbott Molecular, Abbott Park, IL, USA) by use of a Vysis nick translation reagent kit. The labeled product was mixed with 3 mg of human Cot-1 DNA (Invitrogen, Life Technologies Corporation, Carlsbad, CA, USA) and was isolated by means of a standard DNA precipitation method. Cytogenetic pellets were prepared according to standard clinical procedures, and chromosomes and nuclei were visualized by counterstaining with 4',6-diamidino-2-phenylindole. For deletions, at least 10 metaphase cells were analyzed, and interphase nuclei were examined but not counted. For duplications, at least 10 metaphase cells and at least 50 interphase nuclei were analyzed. All FISH probes were tested on metaphase spreads from unaffected individuals to assure proper hybridization.
MLPA was performed using the P070 subtelomeric kit (MRC Holland, Amsterdam, The Netherlands); all positive results were confirmed with another kit (P036B). The procedure was conducted according to the manufacturer's recommendations. Briefly, the patient's DNA was diluted in PCR-grade water and quality was assessed via spectrophotometry (Nanodrop®, Thermo Scientific, Wilmington, DE, USA). The hybridization solution (SALSA probe-mix and MLPA buffer) was added to a final DNA concentration of 60 ng/μl. DNA was denatured at 90°C, then hybridized for 16-20 hours at 60°C. Ligation was performed at 54°C for 15 minutes, and the ligated product was denatured at 98°C for 5 minutes and then amplified by PCR (SALSA PCR buffer, PCR-primers and polymerase). The PCR product was run for fragment analysis on an ABI 3130 sequencer (Applied Biosystems, Life Technologies, Carlsbad, CA, USA). Normalization of the data and analysis of the MLPA results were conducted using Coffalyser v3.5 software provided by MRC Holland (Amsterdam, The Netherlands).
AGH was performed with Agilent® 244 K oligonucleotide arrays (Agilent Technologies Inc., Santa Clara, CA, USA) according to the manufacturer's instructions. Two arrays were used for each trio, one in which the child's DNA was hybridized against the father's DNA, and another in which the child's DNA was hybridized against the mother's DNA. Captured images were analysed with Feature Extraction v 9.1 and CGH Analytics 3.5.14 (Agilent Technologies Inc., Santa Clara, CA, USA).
Cytogenetic analysis was performed on peripheral blood cultures after preparation and G-banding of metaphase chromosomes using standard clinical methods.
Given the genotypes for the child, mother, and father, errors in mendelian transmission were identified and their frequency compared to normal or technical error rates, which were very low. Essentially, the occurrences of an AA parent with a BB child, or vice versa, were counted. Graphical tools were developed to distinguish heterodisomy and isodisomy by visualization and to determine the parent of origin .
Confirmation of uniparental disomy was obtained by genotyping microsatellite markers chosen for high heterozygosity values. The following markers were genotyped and were informative for chromosome 11: D11S1363 (AFMA134WH5), D11S4046 (AFMB042YF5), D11S1318 (AFM218XE1), D11S4088 (AFMA155TE9) and D11S4146 (AFMB072WE5). The following markers were genotyped and were informative for chromosome 16: D16S3093 (AFMB308ZH9), D16S409 (AFM161XA1), D16S515 (AFM340YE5) and D16S402 (AFM031XA5). 25 ng of DNA diluted in 10:1 TE buffer were amplified for each PCR reaction. Primers fluorescently labeled with either 6-fam or HEX (Applied Biosystems, Life Technologies, Carlsbad, CA, USA) were used in conjunction with AmpliTaq Gold® PCR kit reagents (Applied Biosystems, Life Technologies, Carlsbad, CA, USA). PCR was performed using the following steps: 95°C for 10 minutes; 30 cycles of 95°C for 30 seconds, 55°C for 30 seconds, and 72°C for 90 seconds; and a final step at 72°C for 7 minutes. The resulting PCR product in a volume of 1 μl was combined with 9 μl formamide (Applied Biosystems, Life Technologies, Carlsbad, CA, USA) and 0.3 ul GeneScan ROX 500 (Applied Biosystems, Life Technologies, Carlsbad, CA, USA), denatured at 95°C for 5 minutes, quickly chilled in ice and loaded onto a Genetic Analyzer 3130XL (Applied Biosystems, Life Technologies, Carlsbad, CA, USA). Data were visualized using GeneMapper v4.0 (Applied Biosystems, Life Technologies, Carlsbad, CA, USA).
AGH: array genomic hybridization; BAC: bacterial artificial chromosome; bp: base pairs; CNV: copy number variant; del: deletion; FISH: fluorescence in situ hybridization; ID: intellectual disability; K: thousand; Mb: megabase; MLPA: multiplex ligation-dependant probe amplification; PCR: polymerase chain reaction; SMD: significance of mean difference; SNP: single nucleotide polymorphism; UPD: uniparental disomy.
JMF and MAM conceived, designed and coordinated this study. PB developed and implemented the consent process, and JMF, SA, LA, LA, PB, CB, WTG, SL, EL, PM, JLM, BCM, MSP, GAR, MIVA and S-LY recruited patients. Clinical evaluations were performed by LA, LA, CB, JMF, WTG, SL, EL, PM, BCM, MSP, MIVA and S-LY. AB and SC performed the array genomic hybridizations, and AB, ADD, PE, SF, SL, HIL, HQ, MAM and JMF interpreted the AGH data. Confirmatory FISH and MLPA studies were done by PE, DC, EL and JM. LA, LA, CV, PE, WTG, SL, EL, PM, JLM, BCM, MSP, MIVA, S-LY, FRZ and JMF performed the genotype-phenotype correlations. JMF drafted the manuscript, which was critically reviewed and edited by AB, ADD, PE, SL, JLM, MSP, FRZ and MAM. All authors read and approved the final manuscript.
Supplemental Table S1. 54 ID trios from the 100 K cohort on whom 500 K GeneChip® AGH was performed. The findings in families in whom de novo CNVs were found are summarized in Table Table33.
Supplemental Figure S1. Size distribution of pathogenic CNVs detected by 100 K and 500 K AGH in idiopathic ID trios. The number of deletions (red bars) or duplications (blue bars) in each size class is shown for 8 submicroscopic deletions and 2 submicroscopic duplications found in 100 idiopathic ID trios studied with 100 K AGH (empty and stippled bars) and for 22 submicroscopic deletions and 5 submicroscopic duplications found in 154 idiopathic ID trios studied with 500 K AGH (filled bars and stippled bars). 46 trios in the 100 K cohort were studied only with 100 K GeneChip® AGH (empty bars), 54 other trios from the 100 K cohort were studied with both 100 K and 500 K GeneChip® AGH (stippled bars), and 100 additional trios constituted a new cohort who were studied only with 500 K GeneChip® AGH (solid bars). Data for the 100 K trios are from . The 500 K data include two unbalanced reciprocal translocations, for each of which both a deletion and a duplication are shown in the figure.
This work was supported by grants from Genome Canada (JMF, MAM),Genome British Columbia (JMF, MAM), the Canada Foundation for Innovation (JMF), British Columbia Knowledge Development Fund (JMF), Réseau de Médecine Génétique Appliquée of the Fonds de la Recherche en Santé du Québec (JLM, EL) and by the Fonds d'encouragement à la recherche clinique du CHU Sainte-Justine (JLM, EL). The authors thank Jacquie Schein for providing the BACs used for FISH confirmation and Tracy Tucker for her helpful comments regarding the manuscript.