Rare copy number variations (CNVs) are a recognized cause of common human disease. Predicting the genetic element(s) within a small CNV whose copy number loss or gain underlies a specific phenotype might be achieved reasonably rapidly for single patients. Identifying the biological processes that are commonly disrupted within a large patient cohort which possess larger CNVs, however, requires a more objective approach that exploits genomic resources. In this study, we first identified 98 large, rare CNVs within patients exhibiting multiple congenital anomalies. All patients presented with global developmental delay (DD), while other secondary symptoms such as cardiac defects, craniofacial features and seizures were varyingly presented. By applying a robust statistical procedure that matches patients’ clinical phenotypes to laboratory mouse gene knockouts, we were able to strongly implicate anomalies in brain morphology and, separately, in long-term potentiation as manifestations of these DD patients’ disorders. These and other significantly enriched model phenotypes provide insights into the pathoetiology of human DD and behavioral and anatomical secondary symptoms that are specific to DD patients. These enrichments set apart 103 genes, from among thousands overlapped by these CNVs, as strong candidates whose copy number change causally underlies approximately 46% of the cohort's DD syndromes and between 59 and 80% of the cohort's secondary symptoms. We also identified significantly enriched model phenotypes among genes overlapped by CNVs in both DD and learning disability cohorts, indicating a congruent etiology. These results demonstrate the high predictive potential of model organism phenotypes when implicating candidate genes for rare genomic disorders.

We report here on a normal-appearing male with pervasive developmental disorder who was found to have a de novo, apparently balanced complex rearrangement involving chromosomes 6, 10, and 21: 46,XY,ins(21;10)(q11.2;p11.2p13)t(6;21)(p23;q11.2). Further analysis by high-density oligonucleotide microarray was performed, showing an 8.8-Mb heterozygous deletion at 21q21.1-q21.3. Interestingly, the deletion is distal to the translocation breakpoint on chromosome 21. The deletion involves 19 genes, including NCAM2 and GRIK1, both of which are associated with normal brain development and function, and have been considered as possible candidate genes in autism and other neurobehavioral disorders. This case underscores the utility of genomewide microarray analysis for the detection of copy number alterations in patients with apparently balanced complex rearrangements and abnormal phenotypes.

In the present study, DNA from 28 pediatric low-grade astrocytomas was analyzed using Illumina HumanHap550K single-nucleotide polymorphism oligonucleotide arrays. A novel duplication in chromosome band 7q34 was identified in 17 of 22 juvenile pilocytic astrocytomas and three of six fibrillary astrocytomas. The 7q34 duplication spans 2.6 Mb of genomic sequence and contains approximately 20 genes, including two candidate tumor genes, HIPK2 and BRAF. There were no abnormalities in HIPK2, and analysis of two mutation hot-spots in BRAF revealed a V600E mutation in only one tumor without the duplication. Fluorescence in situ hybridization confirmed the 7q34 copy number change and was suggestive of a tandem duplication. Reverse transcription polymerase chain reaction-based sequencing revealed a fusion product between KIAA1549 and BRAF. The predicted fusion product includes the BRAF kinase domain and lacks the auto-inhibitory N-terminus. Western blot analysis revealed phosphorylated mitogen-activated protein kinase (MAPK) protein in tumors with the duplication, consistent with BRAF-induced activation of the pathway. Further studies are required to determine the role of this fusion gene in downstream MAPK signaling and its role in development of pediatric low-grade astrocytomas.

Translational Relevance

Previous reports suggested that abnormalities of INI1 could be detected in 70–75% of malignant rhabdoid tumors. The mechanism of inactivation in the other 25% remained unclear. The goal of this study was to perform a high-resolution genomic analysis of a large series of rhabdoid tumors with the expectation of identifying additional loci related to the initiation or progression of these malignancies. We also developed a comprehensive set of assays, including a new MLPA assay, to interrogate the INI1 locus in 22q11.2. Intragenic deletions could be detected using the Illumina 550K Beadchip, whereas single exon deletions could be detected using MLPA. The current study demonstrates that with a multi-platform approach, alterations at the INI1 locus can be detected in almost all cases. Thus, appropriate molecular genetic testing can be used as an aid in the diagnosis and for treatment planning for most patients.

Purpose

A high-resolution genomic profiling and comprehensive targeted analysis of INI1/SMARCB1 of a large series of pediatric rhabdoid tumors was performed. The aim was to identify regions of copy number change and loss of heterozygosity that might pinpoint additional loci involved in the development or progression of rhabdoid tumors, and define the spectrum of genomic alterations of INI1 in this malignancy.

Experimental Design

A multi-platform approach, utilizing Illumina single nucleotide polymorphism (SNP) based oligonucleotide arrays, multiplex ligation dependent probe amplification (MLPA), fluorescence in situ hybridization (FISH), and coding sequence analysis was used to characterize genome wide copy number changes, loss of heterozygosity, and genomic alterations of INI1/SMARCB1 in a series of pediatric rhabdoid tumors.

Results

The bi-allelic alterations of INI1 that led to inactivation were elucidated in 50 of 51 tumors. INI1 inactivation was demonstrated by a variety of mechanisms, including deletions, mutations, and loss of heterozygosity. The results from the array studies highlighted the complexity of rearrangements of chromosome 22, compared to the low frequency of alterations involving the other chromosomes.

Conclusions

The results from the genome wide SNP-array analysis suggest that INI1 is the primary tumor suppressor gene involved in the development of rhabdoid tumors with no second locus identified. In addition, we did not identify hot spots for the breakpoints in sporadic tumors with deletions of chromosome 22q11.2. By employing a multimodality approach, the wide spectrum of alterations of INI1 can be identified in the majority of patients, which increases the clinical utility of molecular diagnostic testing.

The use of array technology to define chromosome deletions and duplications is bringing us closer to establishing a genotype/phenotype map of genomic copy number alterations. We studied 21 patients and 5 relatives with deletions of the short arm of chromosome 20 using the Illumina HumanHap550 SNP array to 1) more accurately determine the deletion sizes, 2) identify and compare breakpoints, 3) establish genotype/phenotype correlations and 4) investigate the use of the HumanHap550 platform for analysis of chromosome deletions. Deletions ranged from 95kb to 14.62Mb, and all of the breakpoints were unique. Eleven patients had deletions between 95kb and 4Mb and these individuals had normal development, with no anomalies outside of those associated with Alagille syndrome. The proximal and distal boundaries of these eleven deletions constitute a 5.4MB region, and we propose that haploinsufficiency for only 1 of the 12 genes in this region causes phenotypic abnormalities. This defines the JAG1 associated critical region, in which deletions do not confer findings other than those associated with Alagille syndrome. The other 10 patients had deletions between 3.28Mb and 14.62Mb, which extended outside the critical region, and notably, all of these patients, had developmental delay. This group had other findings such as autism, scoliosis and bifid uvula. We identified 47 additional polymorphic genome-wide copy number variants (>20 SNPs), with 0–5 variants called per patient. Deletions of the short arm of chromosome 20 are associated with relatively mild and limited clinical anomalies. The use of SNP arrays provides accurate high-resolution definition of genomic abnormalities.

Two known recurrent constitutional translocations, t(11;22) and t(17;22), as well as a non-recurrent t(4;22), display derivative chromosomes that have joined to a common site within the low copy repeat B (LCR-B) region of 22q11.2. This breakpoint is located between two AT-rich inverted repeats that form a nearly perfect palindrome. Breakpoints within the 11q23, 17q11 and 4q35 partner chromosomes also fall near the center of palindromic sequences. In the present work the breakpoints of a fourth translocation involving LCR-B, a balanced ependymoma-associated t(1;22), were characterized not only to localize this junction relative to known genes, but also to further understand the mechanism underlying these rearrangements. FISH mapping was used to localize the 22q11.2 breakpoint to LCR-B and the 1p21 breakpoint to single BAC clones. STS mapping narrowed the 1p21.2 breakpoint to a 1990 bp AT-rich region, and junction fragments were amplified by nested PCR. Junction fragment-derived sequence indicates that the 1p21.2 breakpoint splits a 278 nt palindrome capable of forming stem–loop secondary structure. In contrast, the 1p21.2 reference genomic sequence from clones in the database does not exhibit this configuration, suggesting a predisposition for regional genomic instability perhaps etiologic for this rearrangement. Given its similarity to known chromosomal fragile site (FRA) sequences, this polymorphic 1p21.2 sequence may represent one of the FRA1 loci. Comparative analysis of the secondary structure of sequences surrounding translocation breakpoints that involve LCR-B with those not involving this region indicate a unique ability of the former to form stem–loop structures. The relative likelihood of forming these configurations appears to be related to the rate of translocation occurrence. Further analysis suggests that constitutional translocations in general occur between sequences of similar melting temperature and propensity for secondary structure.

Palindromic sequences are dispersed in the human genome and may cause chromosomal translocations in humans. They constitute unsequenced gaps in the human genome because of their resistance to PCR amplification, cloning into vectors, and sequencing. We have overcome these difficulties by using a combination of optimized PCR conditions, cloning in a recombination-deficient E. coli strain, and RNA polymerases in sequencing. Using these methods, we analyzed a palindromic AT-rich repeat (PATRR) in the neurofibromatosis type 1 (NF1) gene on chromosome 17 (17PATRR). The 17PATRR manifests a size polymorphism due to a highly variable length of (AT)n dinucleotide repeats within the PATRR. 17PATRRs can be categorized into two types: a longer one that comprises a nearly or completely perfect palindrome, and a shorter one that represents its deleted asymmetric derivative. In vitro analysis shows that the longer 17PATRR is more likely to form a cruciform structure than the shorter one. Two reported t(17;22)(q11;q11) patients with NF1, whose breakpoints were identified within the 17PATRR, have translocations that are derived from perfect or nearly perfect palindromic alleles. This implies that the symmetric structure of a PATRR can induce a translocation. We identified conserved PATRRs within the NF1 gene in great apes and similar inverted repeats in two Old World monkeys, but not in New World monkeys or other mammals. This indicates that the palindromic region appeared approximately 25 million years ago and elongated during primate evolution. Although such palindromic regions are usually unstable and disappear rapidly due to deletion, the 17PATRR in the NF1 gene was stably conserved during evolution for reasons that are still unknown.

Common copy number variations (CNVs) represent a significant source of genetic diversity, yet their influence on phenotypic variability, including disease susceptibility, remains poorly understood. To address this problem in cancer, we performed a genome-wide association study (GWAS) of CNVs in the childhood cancer neuroblastoma, a disease where SNP variations are known to influence susceptibility1,2. We first genotyped 846 Caucasian neuroblastoma patients and 803 healthy Caucasian controls at 550,000 single nucleotide polymorphisms, and performed a CNV-based test for association. We then replicated significant observations in two independent sample sets comprised of a total of 595 cases and 3,357 controls. We identified a common CNV at 1q21.1 associated with neuroblastoma in the discovery set, which was confirmed in both replication sets (Pcombined = 2.97 × 10−17; OR = 2.49, 95% CI: 2.02 to 3.05). This CNV was validated by quantitative PCR, fluorescent in situ hybridization, and analysis of matched tumor specimens, and was shown to be heritable in an independent set of 713 cancer-free trios. We identified a novel transcript within the CNV which showed high sequence similarity to several “Neuroblastoma breakpoint family” (NBPF) genes3,4 and represents a new member of this gene family (NBPFX). This transcript was preferentially expressed in fetal brain and fetal sympathetic nervous tissues, and expression level was strictly correlated with CNV state in neuroblastoma cells. These data demonstrate that inherited copy number variation at 1q21.1 is associated with neuroblastoma and implicate a novel NBPF gene in early tumorigenesis of this childhood cancer.

The human short interspersed repeated element (SINE), Alu, amplifies through a poorly understood RNA-mediated mechanism, termed retroposition. There are over one million copies of Alu per haploid human genome. The copies show some internal variations in sequence and are very heterogeneous in chromosomal environment. However, very few Alu elements actively amplify. The amplification rate has decreased greatly in the last 40 million years. Factors influencing Alu transcription would directly affect an element’s retroposition capability. Therefore, we evaluated several features that might influence expression from individual Alu elements. The influence of various internal sequence variations and 3′ unique flanks on full-length Alu RNA steady-state levels was determined. Alu subfamily diagnostic mutations do not significantly alter the amount of Alu RNA observed. However, sequences containing random mutations throughout the right half of selected genomic Alu elements altered Alu RNA steady-state levels in cultured cells. In addition, sequence variations at the 3′ unique end of the transcript also significantly altered the Alu RNA levels. In general, sequence mutations and 3′ end sequences contribute to Alu RNA levels, suggesting that the master Alu element(s) have a multitude of individual differences that collectively gives them a selective advantage over other Alu elements.

Kabuki syndrome is a rare, multiple malformation disorder characterized by a distinctive facial appearance, cardiac anomalies, skeletal abnormalities, and mild to moderate intellectual disability. Simplex cases make up the vast majority of the reported cases with Kabuki syndrome, but parent-to-child transmission in more than a half-dozen instances indicates that it is an autosomal dominant disorder. We recently reported that Kabuki syndrome is caused by mutations in MLL2, a gene that encodes a Trithorax-group histone methyltransferase, a protein important in the epigenetic control of active chromatin states. Here, we report on the screening of 110 families with Kabuki syndrome. MLL2 mutations were found in 81/110 (74%) of families. In simplex cases for which DNA was available from both parents, 25 mutations were confirmed to be de novo, while a transmitted MLL2 mutation was found in two of three familial cases. The majority of variants found to cause Kabuki syndrome were novel nonsense or frameshift mutations that are predicted to result in haploinsufficiency. The clinical characteristics of MLL2 mutation-positive cases did not differ significantly from MLL2 mutation-negative cases with the exception that renal anomalies were more common in MLL2 mutation-positive cases. These results are important for understanding the phenotypic consequences of MLL2 mutations for individuals and their families as well as for providing a basis for the identification of additional genes for Kabuki syndrome.

Deletions of the 22q11.2 region distal to the 22q11.21 microdeletion syndrome region have recently been described in individuals with mental retardation and congenital anomalies. Because these deletions are mediated by low-copy repeats (LCRs), located distal to the 22q11.21 DiGeorge/velocardiofacial microdeletion region, duplications are predicted to occur with a frequency equal to the deletion. However, few microduplications of this region have been reported. We report the identification of 18 individuals with microduplications of 22q11.21–q11.23. The duplication boundaries for all individuals are within LCRs distal to the DiGeorge/velocardiofacial microdeletion region. Clinical records for nine subjects reveal shared characteristics, but also several examples of contradicting clinical features (e.g. macrocephaly versus microcephaly and upslanting versus downslanting palpebral fissures). Of 12 cases for whom parental DNA samples were available for testing, one is de novo and 11 inherited the microduplication from a parent, three of whom reportedly have learning problems or developmental delay. The variable phenotypes and preponderance of familial cases obfuscate the clinical relevance of the molecular data and emphasize the need for careful parental assessments and clinical correlations.

Background

Recent studies have shown that copy number variations (CNVs) are frequent in higher eukaryotes and associated with a substantial portion of inherited and acquired risk for various human diseases. The increasing availability of high-resolution genome surveillance platforms provides opportunity for rapidly assessing research and clinical samples for CNV content, as well as for determining the potential pathogenicity of identified variants. However, few informatics tools for accurate and efficient CNV detection and assessment currently exist.

Results

We developed a suite of software tools and resources (CNV Workshop) for automated, genome-wide CNV detection from a variety of SNP array platforms. CNV Workshop includes three major components: detection, annotation, and presentation of structural variants from genome array data. CNV detection utilizes a robust and genotype-specific extension of the Circular Binary Segmentation algorithm, and the use of additional detection algorithms is supported. Predicted CNVs are captured in a MySQL database that supports cohort-based projects and incorporates a secure user authentication layer and user/admin roles. To assist with determination of pathogenicity, detected CNVs are also annotated automatically for gene content, known disease loci, and gene-based literature references. Results are easily queried, sorted, filtered, and visualized via a web-based presentation layer that includes a GBrowse-based graphical representation of CNV content and relevant public data, integration with the UCSC Genome Browser, and tabular displays of genomic attributes for each CNV.

Conclusions

To our knowledge, CNV Workshop represents the first cohesive and convenient platform for detection, annotation, and assessment of the biological and clinical significance of structural variants. CNV Workshop has been successfully utilized for assessment of genomic variation in healthy individuals and disease cohorts and is an ideal platform for coordinating multiple associated projects.

Availability and Implementation

Available on the web at: http://sourceforge.net/projects/cnv

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACD/MPV) is a rare, neonatally lethal developmental disorder of the lung with defining histologic abnormalities typically associated with multiple congenital anomalies (MCA). Using array CGH analysis, we have identified six overlapping microdeletions encompassing the FOX transcription factor gene cluster in chromosome 16q24.1q24.2 in patients with ACD/MPV and MCA. Subsequently, we have identified four different heterozygous mutations (frameshift, nonsense, and no-stop) in the candidate FOXF1 gene in unrelated patients with sporadic ACD/MPV and MCA. Custom-designed, high-resolution microarray analysis of additional ACD/MPV samples revealed one microdeletion harboring FOXF1 and two distinct microdeletions upstream of FOXF1, implicating a position effect. DNA sequence analysis revealed that in six of nine deletions, both breakpoints occurred in the portions of Alu elements showing eight to 43 base pairs of perfect microhomology, suggesting replication error Microhomology-Mediated Break-Induced Replication (MMBIR)/Fork Stalling and Template Switching (FoSTeS) as a mechanism of their formation. In contrast to the association of point mutations in FOXF1 with bowel malrotation, microdeletions of FOXF1 were associated with hypoplastic left heart syndrome and gastrointestinal atresias, probably due to haploinsufficiency for the neighboring FOXC2 and FOXL1 genes. These differences reveal the phenotypic consequences of gene alterations in cis.

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACD/MPV) is a rare, neonatally lethal developmental disorder of the lung with defining histologic abnormalities typically associated with multiple congenital anomalies (MCA). Using array CGH analysis, we have identified six overlapping microdeletions encompassing the FOX transcription factor gene cluster in chromosome 16q24.1q24.2 in patients with ACD/MPV and MCA. Subsequently, we have identified four different heterozygous mutations (frameshift, nonsense, and no-stop) in the candidate FOXF1 gene in unrelated patients with sporadic ACD/MPV and MCA. Custom-designed, high-resolution microarray analysis of additional ACD/MPV samples revealed one microdeletion harboring FOXF1 and two distinct microdeletions upstream of FOXF1, implicating a position effect. DNA sequence analysis revealed that in six of nine deletions, both breakpoints occurred in the portions of Alu elements showing eight to 43 base pairs of perfect microhomology, suggesting replication error Microhomology-Mediated Break-Induced Replication (MMBIR)/Fork Stalling and Template Switching (FoSTeS) as a mechanism of their formation. In contrast to the association of point mutations in FOXF1 with bowel malrotation, microdeletions of FOXF1 were associated with hypoplastic left heart syndrome and gastrointestinal atresias, probably due to haploinsufficiency for the neighboring FOXC2 and FOXL1 genes. These differences reveal the phenotypic consequences of gene alterations in cis.