# Detection of Chromosomal Breakpoints in Patients with Developmental Delay and Speech Disorders

Kagistia H. Utami¹,¹³, Axel M. Hillmer², Irene Aksoy³, Elaine G. Y. Chew², Audrey S. M. Teo², Zhenshui Zhang², Charlie W. H. Lee⁴, Pauline J. Chen⁴, Chan Chee Seng⁵, Pramila N. Ariyaratne⁴, Sigrid L. Rouam⁴, Lim Seong Soo¹, Saira Yousof⁶,⁷, Ivan Prokudin⁶,⁷, Gregory Peters⁸, Felicity Collins⁹, Meredith Wilson⁹, Alyson Kakakios¹⁰, Georges Haddad¹¹, Arnaud Menuet¹², Olivier Perche¹², Stacey Kiat Hong Tay¹³, Ken W. K. Sung⁴, Xiaoan Ruan¹⁴, Yijun Ruan¹⁴, Edison T. Liu², Sylvain Briault¹², Robyn V. Jamieson⁶, Sonia Davila¹, Valere Cacheux¹*

¹Human Genetics, Genome Institute of Singapore, Singapore, ²Cancer Therapeutics and Stratified Oncology, Genome Institute of Singapore, Singapore, ³Stem Cells and Developmental Biology, Genome Institute of Singapore, Singapore, ⁴Computational and Mathematical Biology, Genome Institute of Singapore, Singapore, ⁵Scientific & Research Computing, Genome Institute of Singapore, Singapore, ⁶Eye and Developmental Genetics Research, The Children's Hospital at Westmead, Children's Medical Research Institute and Save Sight Institute, Sydney, New South Wales, Australia, ⁷Disciplines of Paediatrics and Child Health and Genetic Medicine, Sydney Medical School, University of Sydney, Sydney, New South Wales, Australia, ⁸Department of Cytogenetics, The Children's Hospital at Westmead, Sydney, New South Wales, Australia, ⁹Department of Clinical Genetics, The Children's Hospital at Westmead, Sydney, New South Wales, Australia, ¹⁰Department of Immunology, The Children's Hospital at Westmead, Sydney, New South Wales, Australia, ¹¹Unité de Génétique, Centre Hospitalier, Blois, France, ¹²Service de Genetique INEM UMR7355 CNRS-University, Centre Hospitalier Régional d'Orléans, Orléans, France, ¹³Department of Paediatrics, Yong Loo Lin School of Medicine, National University of Singapore, Singapore, ¹⁴Genome Technology and Biology, Genome Institute of Singapore, Singapore

## Abstract

Delineating candidate genes at the chromosomal breakpoint regions in the apparently balanced chromosome rearrangements (ABCR) has been shown to be more effective with the emergence of next-generation sequencing (NGS) technologies. We employed a large-insert (7–11 kb) paired-end tag sequencing technology (DNA-PET) to systematically analyze genome of four patients harbouring cytogenetically defined ABCR with neurodevelopmental symptoms, including developmental delay (DD) and speech disorders. We characterized structural variants (SVs) specific to each individual, including those matching the chromosomal breakpoints. Refinement of these regions by Sanger sequencing resulted in the identification of five disrupted genes in three individuals: guanine nucleotide binding protein, q polypeptide (GNAQ), RNA-binding protein, fox-1 homolog (RBFOX3), unc-5 homolog D (C.elegans) (UNC5D), transmembrane protein 47 (TMEM47), and X-linked inhibitor of apoptosis (XIAP). Among them, XIAP is the causative gene for the immunodeficiency phenotype seen in the patient. The remaining genes displayed specific expression in the fetal brain and have known biologically relevant functions in brain development, suggesting putative candidate genes for neurodevelopmental phenotypes. This study demonstrates the application of NGS technologies in mapping individual gene disruptions in ABCR as a resource for deciphering candidate genes in human neurodevelopmental disorders (NDDs).

**Citation:** Utami KH, Hillmer AM, Aksoy I, Chew EGY, Teo ASM, et al. (2014) Detection of Chromosomal Breakpoints in Patients with Developmental Delay and Speech Disorders. PLoS ONE 9(3): e90852. doi:10.1371/journal.pone.0090852

**Editor:** Osman El-Maarri, University of Bonn, Institut of experimental hematology and transfusion medicine, Germany

**Received:** November 7, 2013; **Accepted:** February 4, 2014; **Published:** March 5, 2014

**Copyright:** © 2014 Utami et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This study is supported by grants from Biomedical Medical Research Council (BMRC) of the Agency for Science, Technology and Research (A\*STAR), Singapore. Support is also acknowledged from the National Health and Medical Research Australia (NHMRC). KHU is funded by Singapore International Graduate Award (SINGA) fellowship. Additional support was also provided by the Genome Institute of Singapore internal research funds from the BMRC. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: valerecacheux@gmail.com

## Introduction

Apparently balanced chromosomal rearrangements (ABCR) occur sporadically in the population or segregate within families, with a frequency of 1 in every 2000 live births. ABCR have been largely associated with infertility, ovarian failure and intellectual disability (ID) when detected during pre- and post-natal investigation or genetic counselling. [1,2,3,4,5,6] The risk of developing congenital anomalies or NDDs has been estimated to be 6.1% for *de novo* ABCR. [2] Although ABCR is also found in phenotypically normal individuals, there is an increased incidence of ABCR in NDDs patients, which may lead to novel candidate disease gene

identification through breakpoint cloning methods. Such method has been successfully applied in characterizing disease genes including *DMD* in Dystrophin for Duchenne Muscular Dystrophy [7], *DISC1* in Schizophrenia [8,9] and *ATP7A* in Menkes disease. [10,11] Karyotyping, Array-CGH (aCGH) and SNP arrays are currently first-tier diagnostic tools to investigate ABCR in pre- and post-natal settings. [6,12,13] Balanced chromosomal aberrations can only be detected through karyotype observation, although at a low resolution of approximately 5–10 Million base-pairs (Mb), [14] and subsequent Fluorescence *in situ* hybridization (FISH) mapping across the breakpoints for more detailed delineation is restricted to the investigated region. Microarrays technologies are able to
identify copy number changes at considerable resolution but fail to identify copy number neutral rearrangements. [15,16] Recent advances in NGS technologies render it possible to systematically identify genomic rearrangements including copy number neutral events at a base-pair resolution, thus facilitating candidate disease gene identification in the chromosomal breakpoint regions.

We have established a NGS pipeline, referred to as a DNA-PET sequencing, [17,18,19] in which we sequence the short ends (50 bp) of 5' and 3' tags of large-insert sizes between 7 to 11 kb genomic DNA fragments, followed by ligating the two paired-ends to form PET constructs and subjected to sequencing in a massive and highly parallel manner, (Figure S1, see Supporting Information S1) [17,20]. Hillmer and colleagues [17,18,19] have demonstrated that large-insert fragment sizes provided higher physical coverage with minimum sequencing efforts, and have the advantage over short-insert sizes in terms of large genomic SVs detection and covering complicated DNA sequence features such as repetitive regions. We implemented this technique to analyze the genome of four individuals with NDDs symptoms including Developmental Delay (DD), Speech Delay (SD), Language Delay (LD) and autistic disorder. These patients harbour cytogenetically defined ABCR (2 familial translocation, 1 familial inversion and 1 *de novo* inversion), and prior aCGH analysis did not reveal chromosomal imbalances. Our study shows that NGS technology enables rapid identification of individual gene disruptions and potential candidate genes in ABCR. We also demonstrated the correlation of disrupted gene *XLAP* in an inversion breakpoint to be causative for the patient's immunodeficiency phenotype.

## Results

### Structural Variants Detection by DNA-PET Sequencing in Four Patients with Developmental Delay and Speech Disorders

By using patients' genomic DNA as starting material, libraries were generated and sequenced using a SOLiD platform, and we obtained an average of 35 million non-redundant paired-end reads (Table S1, see Supporting Information S1). The median physical coverage using our technique was 98x, with an average of 94x (Table S1, see Supporting Information S1). Majority of the Paired End Tags (PETs) were mapped to the reference genome NCBI Build 36 and referred to as concordant PETs (cPETS), which provided the copy number information based on sequencing read-depth. The remaining clustered PETs were referred to as discordantly mapped PETs (dPETs), which allowed the identification of structural variants (SVs). After filtering (see Methods section), approximately 96% of the SV calls generated for each library were shared with normal individuals published in the Database of Genomic Variants [21], 1000 Genome Project SVs Pilot release set [22,23] or previous paired-end sequencing studies of normal individuals [24,25,26].

Using our analysis pipeline for patient-specific SVs discovery, the extraction of normal SVs reduced this number to 7–19 events, with a mean of 14 SVs per patient (Tables S2–S6, see Supporting Information S1). We observed deletions as being the most frequent SVs, comprising 58% of the total patients-specific SVs (32 deletions out of a total of 55 SVs), while tandem duplication comprises 20%, and the remaining SVs (translocation, inversion and insertion) comprise 21% of total SVs (Table 1). We performed PCR analysis on 36 randomly chosen SVs with 2 sets of primers spanning predicted breakpoint junctions. Of these, 27 SVs were validated and produced a single, clear PCR band at expected size range, suggesting ~75% validation rate. In parallel to DNA-PET, we compared the copy number detection overlap to aCGH and

<table><caption>Table 1. List of SVs obtained after filtration of four patients.</caption><thead><tr><th rowspan="2">Family</th><th rowspan="2">Patient</th><th rowspan="2">Karyotype</th><th rowspan="2">dPET clusters</th><th rowspan="2">SVs after data curation</th><th rowspan="2">SVs overlap with normal</th><th colspan="7">Patient-specific SVs after filtration</th></tr><tr><th>Total</th><th>Del</th><th>Balance Transloc.</th><th>Isolated Transloc.</th><th>Unpaired Inversion</th><th>Paired Inversion</th><th>Tandem Duplication</th><th>Insertion</th></tr></thead><tbody><tr><td>1</td><td>CD5</td><td>46, XY,t(9;17) (q12;q24)</td><td>505</td><td>284</td><td>270 (95.07%)</td><td>14</td><td>12</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>2</td><td>CD10</td><td>46,XY, t(6;8) (q16.2,p11.2)</td><td>927</td><td>616</td><td>601 (97.50%)</td><td>15</td><td>9</td><td>2</td><td>1</td><td>0</td><td>1</td><td>2</td></tr><tr><td>3</td><td>CD8</td><td>46,XY, inv(X)(p22q26)</td><td>724</td><td>387</td><td>368 (95.09%)</td><td>19</td><td>9</td><td>0</td><td>0</td><td>1</td><td>2</td><td>7</td></tr><tr><td>4</td><td>CD9</td><td>46,XY, inv(5)(q22q35.1)</td><td>930</td><td>581</td><td>574 (98.79%)</td><td>7</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2</td><td>3</td></tr><tr><td colspan="2">Total</td><td></td><td>3,086</td><td>1,868</td><td>1,813</td><td>55</td><td>32 (58.1%)</td><td>4 (7.2%)</td><td>1 (1.8%)</td><td>4 (7.2%)</td><td>11 (20%)</td><td>2 (3.6%)</td></tr></tbody></table>

<footnote><sup>a</sup>Discordantly mapped PETs (dPETs) which connect the same two genomic regions are clustered together.</footnote>

<footnote><sup>b</sup>dPET clusters which has passed quality filters to remove sequencing artefacts provided the list of predicted SVs.</footnote>

<footnote><sup>c</sup>Approximately 95% of SVs overlapped with SVs found in normal population (see Methods).</footnote>

<footnote>doi:10.1371/journal.pone.0090852.t001</footnote>
observed an average of 93% overlap between both experiments (Table S7, see Supporting Information S1).

For each patient, we performed a case by case study to identify the potential genes which were disrupted by the cytogenetically visible rearrangements with the assumption that these were the most likely causative events. We validated DNA-PET breakpoints matching the cytogenetic rearrangements by Sanger sequencing (Table 2), FISH and evaluated gene expression by quantitative real-time PCR (RT-PCR). Apart from the specific chromosomal rearrangements, we checked other SVs obtained for each patient. Most of the SVs that coincide with coding regions have been reported in the Database of Genomic Variants (DGV), indicating that these SVs are likely benign (Tables S2–6, see Supporting Information S1). We checked for the presence of regulatory elements, such as transcription factor binding sites and epigenetic marks, using RegulomeDB and the ENCODE dataset, on SVs present in intra/intergenic regions and found no evidence of known regulatory elements in any of them.

From these extensive analyses, we identified five disrupted genes in three patients: GNAQ, RBFOX3, UNC5D, XIAP and TMEM47 (Table 2). Among them, XIAP has already been associated with X-linked lymphoproliferative disorder 2 (XLP-2: MIM [300635]). To evaluate the potential implication of the disrupted genes, we assessed the occurrence of pre-existing Copy Number Variations (CNVs) in DD cases from published studies and the DECIPHER Consortium [27,28,29] (Table 3). Enrichment of CNV counts in the cases versus controls was observed in RBFOX3 and UNC5D, suggesting that altered gene dosage in these genes might have a role in NDDs.

## Breakpoint Characterization through Detailed SVs Analysis

**Case 1 (Patient CD5).** In this family, the translocation segregates with a variable degree of phenotypic manifestation. Patient CD5 had language delay (LD) and DD during his childhood, and his two sons (CD21 and CD22) displayed absence of speech, DD and autistic disorder (Figure 1A). The translocation t(9;17) was represented by two abnormally oriented clusters corresponding to the balanced translocation between chromosome 9 and 17. Sanger sequencing refined the breakpoint coordinates in the translocation carriers, and revealed identical breakpoint patterns on chromosome 9 at position chr9:79,572,001–79,572,005 and on chromosome 17 at position chr17:74,764,927–74,767,964. There is a loss of 4 bp and 3,036 bp on chromosome 9 and 17, respectively, with a microhomology of 3 bp between paired breakpoints (Figure 1B).

The translocation disrupted GNAQ at intron 5 on chromosome 9 and RBFOX3 at intron 2 on chromosome 17. Both disrupted genes shared the same orientation, breakpoints lie within introns, and the resulting fusion is predicted to be in frame. However, RT-PCR analysis on the patient's lymphoblastoid cell line did not reveal any fusion transcript expression as neither gene is expressed in these cells, reflecting the importance of relevant cell type for validation of brain-specific genes. We checked the mRNA expression of both genes in human tissue panel and found that both genes are highly expressed in the fetal brain and the cerebellum (Figure 1C). GNAQ encodes a member of the Gαq heterotrimeric protein family, and null Gnaq mice exhibit cerebellar ataxia, and motor coordination deficit [30]. RBFOX3 is a neuronal specific splicing factor, which is exclusively expressed in the neuronal nuclei. There are 12 other SVs found in this patient, and nine of them were validated by PCR, in which 5 SVs were shared with his two affected sons (SV1, 2, 4, 6 and 13). These five SVs overlapped with known CNV regions (SV1, 2, 4) reported in healthy individuals or were in intergenic regions (SV6, 13) and were therefore excluded for further analysis (Table S2, see Supporting Information S1). Overall, the translocation segregates with LD and DD in this family with variable penetrance, hinting a potential functional impact.

**Case 2 (Patient CD10).** In the second patient (CD10), the causative effect of the t(6;8) balanced translocation is unclear, due to variable expression of phenotype in the translocation carriers; his mother is asymptomatic and his younger sibling displayed schizencephaly and DD (CD11) (Figure 2A). DNA-PET identified abnormal paired reads matching the balanced translocation, and Sanger sequencing refined the breakpoint coordinates on chromosome 6 at position chr6:98,318,526–98,318,538 and chromosome 8 at position chr8:35,527,969–35,527,976. The translocation disrupted UNC5D at intron 5 on chromosome 8, while no coding genes were disrupted on chromosome 6 (Figure 2B). UNC5D encodes a member of human dependence receptor UNC5 family that is specifically expressed in the layer 4 of the developing neocortex in rats, which makes it a biologically plausible candidate [31,32]. Since the translocation does not fully segregate with the disease, all coding exons of UNC5D were sequenced in each translocation carriers, and we did not find additional point mutations in the other allele. UNC5D mRNA is exclusively expressed in the adult brain, fetal brain, and cerebellum (Figure 2C), thus expression changes could not be investigated in the patient's lymphoblastoid cell line. Notably, we observed enrichment of CNV counts in 15 cases described by DECIPHER (Figure 2D). Analysis of other SVs showed 5 deletions in intergenic

<table><caption>**Table 2.** Validated breakpoints and disrupted candidate genes from DNA-PET sequencing in four patients.</caption><thead><tr><th>Patient</th><th>Cytogenetic Analysis</th><th>Chr.</th><th>DNA-PET breakpoint predicted coordinate (SOLiD)</th><th>Validated breakpoint coordinate (Sanger)</th><th>Gene</th></tr></thead><tbody><tr><td rowspan="2">CD5</td><td rowspan="2">t(9;17)</td><td>9</td><td>79,571,716–79,573,787</td><td>79,572,001–79,572,005</td><td>GNAQ</td></tr><tr><td>17</td><td>74,764,804–74,768,401</td><td>74,764,927–74,767,964</td><td>RBFOX3</td></tr><tr><td rowspan="2">CD10</td><td rowspan="2">t(6;8)</td><td>6</td><td>98,318,059–98,318,840</td><td>98,318,526–98,318,538</td><td></td></tr><tr><td>8</td><td>35,527,808–35,528,282</td><td>35,527,964–35,527,976</td><td>UNC5D</td></tr><tr><td rowspan="2">CD8</td><td rowspan="2">inv(X)</td><td>X</td><td>119,984,597–119,984,844</td><td>119,984,613–119,984,615</td><td></td></tr><tr><td>X</td><td>122,839,027–122,845,538</td><td>122,839,398–122,844,862</td><td>XIAP</td></tr><tr><td rowspan="2">CD9</td><td rowspan="2">inv(S)</td><td>5</td><td>111,962,591–111,963,149</td><td>111,962,767</td><td></td></tr><tr><td>5</td><td>165,575,819–165,576,628</td><td>165,576,568–165,576,574</td><td></td></tr></tbody></table>

doi:10.1371/journal.pone.0090852.t002
<table><caption>**Table 3.** Screening of CNVs in cases/controls from published and public datasets.</caption><thead><tr><th rowspan="2">Patient</th><th colspan="2">Cases</th><th colspan="3">Control</th></tr><tr><th>Disrupted genes</th><th>Cooper et al [28]</th><th>DECIPHER Consortium</th><th>Total</th><th>Cooper et al [28]</th><th>1000 Genome</th><th>Total</th></tr></thead><tbody><tr><td>CD5</td><td>GNAQ</td><td>3</td><td>7</td><td>10</td><td>14</td><td>0</td><td>14</td></tr><tr><td>CD5</td><td>RBFOX3</td><td>68</td><td>11</td><td>79</td><td>1</td><td>4</td><td>5</td></tr><tr><td>CD10</td><td>UNC5D</td><td>2</td><td>15</td><td>17</td><td>0</td><td>0</td><td>0</td></tr><tr><td>CD8</td><td>TMEM47</td><td>0</td><td>10</td><td>10</td><td>0</td><td>4</td><td>4</td></tr></tbody></table>

The total number of cases in Cooper et al. [28] is 15,767 cases, and for DECIPHER is ~17,000 cases [29]. The total number of controls in Cooper et al. is 8329 controls and for 1000 Genome SV Release set is 185 controls [22].
doi:10.1371/journal.pone.0090852.t003

regions (SV1, 2, 11, 12, and 14) and 8 within genes (SV3, 4, 5, 6, 9, 10, 13 and 15), in which all of them overlapped with known CNVs regions (Table S3, see Supporting Information S1).

**Case 3 (Patient CD8).** The third patient (CD8), carried a maternal-derived pericentric inversion at chromosome X, between Xp22 and Xq26 (Figure 3A). His monozygotic twin brother harbors the same inversion. They both show signs of immunodeficiency combined with mild intellectual disability (ID) and LD. SVs analysis of this patient revealed a 1.2 Mb inversion located at chromosome Xq24–25 between paired breakpoints (Table S4, see Supporting Information S1), which did not correspond to the large inversion seen in the karyotype, inv(X)(p11.4;q24). Sanger sequencing allowed us to delineate the precise breakpoint on chromosome X at position chrX:119,984,613–119,984,615 and at position chrX:122,839,398–122,844,862 with a loss of 1 bp and 5,463 bp on each breakpoint junction, respectively. Due to the discrepancies between karyogram and NGS data, we traced other SVs in chromosome X by including lower confidence SVs (cluster size ≥4), and revealed a total of 10 SVs that clustered into multiple breakpoint hotspots on the p arm (SV14, 17, 20, 21), the centromeric region (SV22) and the q arm (SV13, 15, 16, 18, 19), suggesting a complex chromosomal rearrangement mechanism (Table S5, see Supporting Information S1). We verified the breakpoints in the inversion carriers (the mother and the twin boys) by PCR and FISH using probes RP1-315G1 (Xq25), RP1-296G17 (Xq24), RP11-330K13 (Xp21), W12-499N23 (Xq25) BAC and Fosmid clones (Figure 3B). These data suggested two large sequential breaks in each arm of chromosome X (Xp21 and Xq24–25) that resulted in tandem duplications, deletion and isolated breakpoints (Figure 3C). This double-inversion mechanism rearranged the landscape of chromosome X architecture, and appeared as a large inversion on G-banding karyotype. Combining DNA-PET sequencing and extensive FISH validations allowed us to characterize the complex chromosomal rearrangements in chromosome X, as a result of a for a sequential double inversion mechanism.

This complex rearrangement disrupted two genes: *XLAP* and *TMEM47*. The latter encodes a member of PMP22/EMP/ Claudin protein family is ubiquitously expressed in human tissues including adult and fetal brain (Figure 3D). Two other genes were affected by tandem duplications as part of the complex rearrangement (SV18, 19): *SH2D1A*, encodes an SH2 domain containing 1A protein (SAP), which has been associated with X-linked lymphoproliferative syndrome type 1 (XLP-1 [MIM: 308240]). This gene was initially suspected as the primary cause of immunodeficiency phenotype seen in the patient, although it was rejected due to normal SAP expression in patient's blood lysates (see Patients and Methods). *ODZ1* is the second gene affected by this rearrangement. It encodes an Odd Oz/Ten-M homolog 1 of Drosophila pair-rule gene involved in post-segmentation processes, including embryonic development of the central nervous system (CNS), eyes, and limbs in Drosophila [33,34].

RT-PCR of the four genes revealed a total absence of expression for *XIAP1* and *TMEM47* and a moderate to normal expression for *SH2D1A* and *ODZ1* in patient's fibroblast cell line compared to sex-matched controls, supporting the absence of functional *XIAP1* and *TMEM47* genes in the male patient (Figure 3E). Considering the normal expression of SAP protein and the absence of functional gene expression of *XLAP*, these data clearly suggest that XLP-2 underlies the immunodeficiency symptoms of the affected twin boys.

**Case 4 (Patient CD9).** In patient 4 (CD9), sequencing analysis identified two paired-inversion clusters corresponding to
<figure><img src="image_1.png" /><figcaption>**Figure 1.** Patient CD5 with translocation t(9;17). A) The pedigree of patient CD5 is indicated. The translocation is transmitted to his two sons (CD21 and CD22). B) Translocation between chromosome 9 and 17 were validated by Sanger sequencing in three translocation carriers. The reference sequence is indicated, showing the fusion of two genes at the genomic level: the first five exons of GNAQ fused to exon 3–14 of *RBFOX3* and the first two exons of *RBFOX3* fused to exon 6–7 of GNAQ. C) mRNA expression of GNAQ and *RBFOX3* showed high expression in fetal brain, adult brain and cerebellum in human tissue panel.
doi:10.1371/journal.pone.0090852.g001</figcaption></figure>

the large inversion on chromosome 5. This inversion spans 53 Mb on chromosome 5q22.2 (chr5:111,966,591-111,963,149) and 5q34 (chr5:165,575,819-165,576,628), with 1 bp deletion, and micro-homologies of 3 bp between paired breakpoints. No gene was disrupted at the breakpoint junctions (Table S6, see Supporting Information S1), and no evidence of regulatory elements sitting at both breakpoint coordinates was found based on RegulomeDB and ENCODE databases. Additionally, there are 5 other patient-specific SVs; two were located in the intergenic regions (SV5 and 6); two others were tandem duplications overlapped with known CNVs (SV4 and 7); and one intronic deletion of 4 kb in *RNF19B* gene, which has not been listed in DGV (SV1) (Table S6, see Supporting Information S1). Based on this analysis, it seems unlikely that the large inversion on chromosome 5 causes the clinical phenotype. Other patient specific SVs or mutations that cannot be identified by DNA-PET such as point mutations or exposures of environmental factors might be the underlying cause of the phenotypic features seen in this patient.

Discussion

In the past decade, we have seen substantial progresses for identification of novel candidate genes in NDDs with the recent development in technologies. Two studies of large cohorts of patients with developmental delay described the enrichment of large CNVs in 15% of the cases [28], and highlighted the presence of additional large CNVs that co-exist with primary microdeletion/duplication syndrome in 10% of the cases as an additive contributing factor to more severe phenotype [35]. These CNVs have been useful to provide a better classification of microdeletion/duplication syndromes; however these regions often encompass multiple genes and thus make it challenging to identify plausible candidates. Recent exome sequencing study in individuals with ID identified potentially causative *de novo* single nucleotide variants (SNVs) with a diagnostic yield of 16%, comparable to the CNV burden obtained in copy number studies [36]. A more conventional approach in candidate gene identification involves delineating candidate genes in the chromosomal breakpoints of the ABCR [9,37,38,39]. In contrast to the downstream effects of SNVs or CNVs, genes that are disrupted
<figure><img src="image_2.png" /><figcaption>**Figure 2.** Patient CD10 with translocation t(6;8). A) The pedigree of patient CD10 is indicated. The familial translocation is inherited from asymptomatic carrier mother and shared with his affected sister (CD11). B) Sanger sequencing analysis refined the chromosomal breakpoint regions and revealed a loss of 11 bp on chromosome 6 and 8 bp on chromosome 8, with a microhomology of 3 bp between the paired breakpoints. C) *UNC5D* mRNA expression in human tissue panel showed high expression in the fetal brain, adult brain and cerebellum compared to other tissues. D) The translocation breakpoint is located at intron 1 of *UNC5D* indicated by the black arrow, encompasses 15 CNVs cases described in the DECIPHER. doi:10.1371/journal.pone.0090852.g002</figcaption></figure>

by translocations or inversion are presumably more severely affected, resulting generally in protein truncation or heterozygous inactivation of the affected allele.

Recent studies have described the feasibility of using NGS technologies to map the ABCR breakpoints in patients with neurodevelopmental abnormalities [37,38,40,41,42,43,44]. These technologies include (i) a shotgun sequencing approach by using

<figure><img src="image_3.png" /><figcaption>**Figure 3.** Patient CD8 with a complex chromosomal inversion. A) Karyogram of normal chromosome X compared to der(X) in patient CD8. B) FISH validation of 10 SVs shown in Table S5 (see Supporting Information S1) with the respective FISH probes: Hybridization of RP1-296G17-Biot (SV15) and RP1-315-Dig (SV16) were localized on the centromere of the patient's metaphase. Probes for SV17 and SV21 on Xp21 (RP11-330K13-Biot) and Xq25 (W12-499N23-Dig), respectively resulted in a split signal between Xp21 and the centromeric region in the patient's chromosome. Further FISH analysis was performed by using probe RP11-762M23-Biot on Xq11.1 (SV22) that was found to localize on the upper chromosomal arm. Probe RP11-655E22 on Xp11.2 was localized on the lower arm of derivative chromosome X. C) Reconstructed derivative chromosome X for patient CD8. Normal human chromosome X according to ISCN 2009 with the arrow orientation from a to d and the proposed mechanism of sequential double inversion in patient CD8. Based on our FISH analysis, an inversion occurred first between Xp21 and Xq25, changing the orientation of p and q arm with a shift of the centromere position towards the lower q-arm shown by inverted red arrow b and c. This was followed by the second inversion that occurred between Xq11.1 and Xq25, altering the orientation of the q-arm (inverted green arrow c). D) Expression of *TMEM47* in human tissue panel assessed by qRT-PCR. E) Expression analysis of four disrupted genes in patient CD8 assessed by qRT-PCR. doi:10.1371/journal.pone.0090852.g003</figcaption></figure>
the flow-sorted derivative chromosomes [40], (ii) custom jumping libraries coupled to targeted breakpoint capture [38], which are limited to the chromosomal breakpoint regions, (iii) standard paired-end sequencing or large-insert jumping libraries of 3–4 kb [37,38,41], and (iv) mate-pair sequencing of 2–3 kb insert sizes [42,43].

Our work described the use of a genome paired end tag (DNA-PET) sequencing with larger insert sizes between 7–11 kb [17,18,19,45,46]. The use of an approximately 8–15 kb insert sizes has been shown to be more advantageous in terms of SVs detection in more complicated DNA sequence features, such as repetitive regions or large genomic rearrangements, and also provides higher physical coverage with minimum sequencing efforts compared to smaller insert sizes [19,47]. We implemented this technique to map the breakpoints of four patients with DD harbouring cytogenetically defined ABCR, and identified specific breakpoints for all of them. Sanger sequencing was required to refine the breakpoints at the base pair level. The observed cryptic breakpoint anomalies were deletions ranging from 3 bp to 5,462 bp and/or microhomologies of 2–4 bp suggesting a microhomology-mediated end joining (MMEJ) [48] or non-homologous end joining mechanism (NHEJ) [49]. Both are major pathways for double strand break repair that often occurs in non-homologous regions of non-recurrent chromosomal breakpoints emphasizing the unique characteristics of these rearrangements.

In this study, five genes were identified within the ABCR-breakpoint regions in three patients (*GNAQ*, *RBFOX3*, *UNC5D*, *XLAP*, and *TMEM47*), and we observed various complications in attempting to correlate these genes to the expressed phenotypes. For one case (Patient CD8), the complex inversion pattern seen in chromosome X resembles chromothripsis, phenomenon commonly found in cancer and recently, in congenital diseases, which resulted from localized shattering of one or few chromosomes and assembly of chromosomal pieces by non-homologous end-joining (NHEJ), and thus appeared as complex chromosome rearrangements, involving three or more breakpoints [41,42,43,50,51]. This complex events highlight the advantage but also the limitation of the DNA-PET technology as FISH experiments were necessary to better understand the high order structure of the rearrangement. The inversion breakpoint disrupted *XLAP* gene, which is the primary cause of a rare form of XLP2 [52,53] (XLP2 [MIM: 300635]) and likely to be causative for the immunodeficiency phenotype of the patient. More interestingly, these patients also presented mild ID and LD, features that are occasionally not associated with XLP syndromes, suggesting possibilities of additional contributing candidate genes within the complex rearrangements. *TMEM47* is among the potentially interesting candidate, with total absence of transcript expression in patient's cell line and expression in fetal to adult brain. We proposed the need to perform mutational screening of *TMEM47* in X-linked ID patients for future studies.

For one patient (Patient CD9), there were no disrupted genes or regulatory elements at the breakpoint regions, neither potential other interesting SVs, showing the limitation for candidate gene detection merely through ABCR breakpoint cloning methods. Alternatively, exome sequencing of parent-offspring trio can be considered as a method of choice to further investigate possible causal variants especially in cases where there is a lack of association between ABCR and disease.

For the two remaining familial translocations, we identified fusion gene at the genomic level, although unverifiable at the transcript level due to lack of expression in available cell lines between *GNAQ* and *RBFOX3* genes in a t(9;17) translocation in three affected members of a family (Patient CD5) with different

level of neurological symptoms (mild to severe) and *UNC5D* disruption in a family harbouring a t(6;8) translocation carried by two affected siblings (Patient CD10) and his mother. The translocation carriers in latter family displayed a broad range of clinical presentations; an asymptomatic mother, her first child with mild DD, and her second child presenting schizencephaly, polymicrogyria and LD, suggesting additional etiological factor underlying these features apart from the t(6;8) translocation. Besides the large phenotypic variability between translocation carriers seen in both families, we observed an enrichment of CNV counts for *UNC5D* and *RBFOX3* in NDDs-associated cases [28]. Therefore, we cannot totally exclude the possibility of high level of variability effects in these genes. Independent validation screening in isolated neurologically-affected patients, rather than collective multi-symptoms cohorts, would be significantly useful to provide significant association with specific neurodevelopmental symptoms in these patients.

Despite the potential functions of four identified genes (*GNAQ*, *RBFOX3*, *UNC5D* and *TMEM47*) in brain development, further functional validations are required to establish a correlation between these disrupted genes and patients' clinical phenotype. Furthermore, reporting such ABCR-disrupted genes is crucial to determine the clinical relevance of newly identified CNVs or SNVs encompassing these genes.

In this study, we showed the potential of a cost-effective NGS technology to rapidly pinpoint disrupted genes within ABCR breakpoint regions. High overlap (93%) between both NGS-and array-based technologies shown in this study emphasizes the importance of employing a single technique that can provide high coverage of genome information with consistent validation rate as a routine clinical diagnosis tool. The stringency of our filtering pipeline can be optimized and adjusted depending on the sequencing platforms, and this would revolutionize the characterization of individual genomes of patients without prior karyotypic observation. For the implementation in the clinical settings, unaffected parents or siblings should be included for genome investigation to reduce the number of familial, non-pathogenic SVs. Whilst the technology continues to improve, validation with longer read depth by Sanger sequencing is still necessary for better SVs annotation. Our study complements the existing application of NGS technology in unexplained NDDs patients for better characterization of chromosomal rearrangements and discovery of potential candidate genes.

## Patients and Methods

### Patients

All samples and information were collected after written informed consent from patient's parents was obtained and in accordance with local institutional review board approved protocols from National University of Singapore in Singapore, Children's Hospital Westmead in Sydney, Australia and Centre Hospitalier Regional in Orleans, France. DNA samples were obtained from peripheral blood lymphocytes and cultured skin fibroblasts obtained from patients seen at the participating institutes.

**Patient CD5.** This is a familial balanced translocation presenting variable degree of DD and autistic features. The first son displayed an autistic behavior and global DD at three years of age with an absence of speech, feeding and sleeping difficulties, habit disorders, and stereotypic movements. At 4 years, there is an improvement in communication and speech seen in the first son. Chromosome analysis revealed a translocation t(9;17) (see Table 1), which is shared with his father and sibling. During genetic
counselling, the father reported that he suffered from LD and DD during childhood that was not explored at that time, and this has resolved by adolescence. His second son was found to have DD and autistic features at the age of two years. DNA tested for the translocation was obtained from the father who was referred to as patient CD5. Chromosome analysis in the phenotypically normal mother revealed a normal karyotype.

**Patient CD10.** The patient was born at term by lower segment caesarean section, due to breech presentation. Delay in developmental milestones was noted at 18 months of age affecting both walking and speech. At 9 years old, comprehension and behavioural difficulties were noted at school. Karyotype analysis revealed a balanced translocation t(6;8) (see Table 1). Metabolic screening and FRAXA testing were normal. Karyotypic analysis in his parents revealed that his mother carried the same balanced translocation. She had no intellectual difficulties, but was reported to have had a 'hole in the heart' in childhood, which closed spontaneously. His younger sister had the same translocation detected by amniocentesis during pregnancy. She was noted to have plagiocephaly soon after birth. She had feeding difficulties in early infancy, which gradually resolved. At 2 years old, she had LD with low-average fine motor and gross motor skills. A right intermittent exotropia was noted. An MRI head scan showed a closed lip schizencephaly involving the right frontal lobe with polymicrogyria in the right Sylvian fissure. MRI brain scans in her brother and mother were normal.

**Patient CD8.** The patient was one of monozygotic twins. At the age of three years, he had an acute EBV infection with prolonged hepatosplenomegaly and abnormality of liver function tests. Immunological investigations showed reduced IgM, decreased CD4 T-helper cells and decreased natural killer (NK) cell function. Patient was suspected for XLP-1, but western blot analysis on patient's whole blood cell lysates showed normal expression of SAP. At the age of five years, mild ID was diagnosed, as well as a specific language disorder affecting his receptive and expressive skills. Karyotype analysis revealed a complex inversion inv(X) (see Table 1), derived from his phenotypically normal mother. His identical twin showed similar clinical features and carrying the same karyotype. The maternal uncle was reported to have recurrent infections and a maternal great uncle died at eight months of age, apparently due to liver abnormalities.

**Patient CD9.** The patient was born by normal delivery after an uneventful pregnancy. Mild delay in early motor and speech developmental milestones was reported. At age twelve, he was found to have average intellectual ability, an expressive language disorder, and specific learning difficulties in maths, spellings and reading; as well as attention deficit disorder. Chromosome analysis revealed a paracentric inversion of chromosome 5 (see Table 1). Both parents have normal karyotypes and do not show any phenotypic abnormality.

## Cytogenetic and FISH Analysis

Karyotypes were determined from G-banding analysis using standard protocol according to the ISCN nomenclature. FISH analysis was carried out using protocols as described elsewhere. [54,55] Metaphase chromosomes were prepared from Epstein-Barr virus-transformed lymphoblastoid cell lines (EBV-LCL) or cultured skin fibroblast cells obtained from patients, parents or siblings carriers by standard techniques. BAC and fosmid probes were obtained from BACPAC Resources (Oakland, CA). Probes were labelled by nick-translation kit (Enzo) with biotin-16-dUTP or digoxigenin-11-dUTP (Roche). The probes were blocked with 1 µg/µl Cot-1 DNA (Life Technologies), and resuspended at a concentration of 5 ng/µl in hybridization buffer (2xSSC, 10%

Dextran Sulfate, 1x PBS, 50% Formamide). Fluorescent signals were visualized by avidin-conjugated fluorescein isothiocyanate (FITC) (Vector Laboratories, CA) or anti-Digoxigenin-Rhodamine (Roche). Chromosomes were counter-stained with DAPI and the signal was analysed using a Nikon Epifluorescence Microscope equipped with ISIS Metasystems for imaging analysis.

## Genomic DNA Preparation

Genomic DNA from lymphoblastoid cell lines, fibroblast, or blood from the patient was extracted by Qiagen Blood and Cell Culture DNA Kits (Qiagen), according to manufacturer's instruction. Lymphoblastoid and fibroblast cell lines were maintained to a minimum number of passages prior to DNA extraction. Quality and quantity of the extracted DNA were measured using NanoDrop 1000 Spectrophotometer and agarose gel electrophoresis.

### aCGH

aCGH was performed in all patient samples using the SurePrint G3 Human 2 x 400 k aCGH Microarray (Agilent Technologies Inc. Santa Clara) according to the manufacturer's instruction. The microarray slides were scanned on an Agilent Microarray Scanner. Data were processed by Genomic Workbench software, standard edition 5.0.14 (Agilent). We used the Aberration Detection Method (ADM-2) algorithm to identify DNA copy number variations (CNV). The ADM-2 algorithm identifies all aberrant intervals in a given sample with consistently high or low log ratio based on the statistical score. We applied a filtering option of minimum of three probes in region and centralization threshold of 6. We used NCBI Build 36 as a reference genome. For smaller CNVs that were identified by DNA-PET, we compared with the aCGH raw data and interpreted the copy number change by looking at the individual probe ratio, using a cutoff of 0.1 (<0.1 indicated deletion, >0.1 indicated duplication).

## DNA-PET Library Construction

We constructed the DNA-PET libraries for four patients according to Method described in Hillmer et al. [17]. Briefly, genomic DNA was hydrosheared to 7–11 kb DNA fragments. Long Mate Paired (LMP) cap adaptors were ligated to the hydrosheared and end-repaired DNA fragments. The cap adaptor-ligated DNA fragments were separated by agarose gel electrophoresis, recovered using the QIAEXII Gel Extraction Kit (QIAGEN) and circularized with a biotinylated adaptor that connects the cap adaptors at both ends of the DNA fragments. Missing 5' phosphate groups of cap adaptors created a nick on each strand after circularization of the DNA. Both nicks were translated outwards by >50 bp into the circularized genomic DNA fragment by DNA polymerase I (NEB). The nick-translated constructs were then digested with T7 exonuclease and S1 nuclease (NEB), to release paired-end tag (PETs) library constructs. These constructs were ligated with SOLiD sequencing adaptors P1 and P2 (Life Technologies), and amplified using 2x HF Phusion Master Mix (Finnzymes OY) for sequencing. High throughput sequencing of the 50 bp libraries was performed on SOLiD sequencers (v3plus and v4, respectively) according to the manufacturer's recommendation (Life Technologies). Sequence tags were mapped to the human reference sequence (NCBI Build 36) and paired using SOLiD System Analysis Pipeline Tool Bioscope, allowing up to 12 color code mismatches per 50 bp tag. For sample CD5 and CD8, two DNA size fractions were merged for library construction which resulted in a reduced sensitivity to identify small deletions.
## DNA-PET Data Curation

The majority of the PET sequences mapped accordingly to the reference genome (concordant PETs or cPETs) with expected mapping orientation (5' tag to 3' tag) and expected mapping distance (according to the selected fragment size) The distribution of cPET in the genome was used to retrieve copy number information.⁵ The remaining portion of the PETs mapped discordantly to the reference genome ( discordant PET or dPETs), classified as those with incorrect paired-tag orientation and incorrect genomic distances. These dPETs provided information to search for genomic rearrangements; with specific criteria for different types of SVs according to PET mapping orientation and genomic region as described by Yao et al. [19]. The overlapping dPETs representing similar SVs were clustered together as dPET clusters and counted as the cluster size, according to the procedure as described in Hillmer et al. [17] with refined data curation as described in Ng et al. [20].

## Filtering of Normal Structural Variations (SVs)

Comparison of clusters across different genomes was performed as described by Ng et al. [20]. We included DNA-PET data of 23 normal individuals (25 DNA-PET data sets) and the pilot release set of 1000 Genome Project [22] from Mills et al. in the cross-genome comparison, and identified SVs that were present in the normal libraries. In addition, we used the breakpoint locations to compare the identified SVs with published SVs based on paired-end sequencing studies of 18 additional normal individuals [24,25] and Database of Genomic Variants. The fraction of predicted SV which overlapped with a published SV was calculated by the percentage of overlap relative to the larger event. Thus, we categorized SVs that overlapped by 80% or more with those identified by these studies as normal SVs. Hence, SVs classified as normal have been excluded to identify rearrangements which underlie the diseases.

## CNVs Screening in Public Datasets

We used the CNVs map from published data by Cooper et al. with 15,767 cases of DD and 8,329 adult controls to screen for deletions or duplications in our candidate genes [28]. We also screened from DECIPHER databases, with >17,000 cases carrying CNVs disrupting individual genes. For additional control dataset, we screened for normal CNVs in the first release SVs set of 1000 Genome Project Consortium of 185 individuals [22]. CNVs were counted in both cases and controls spanning candidate genes in these datasets.

## In silico Analysis of Regulatory Regions

The hg18 coordinates of genomic regions from each patient's SVs list were converted to hg19 in LiftOver from UCSC genome browser. SVs coincide within introns or intergenic regions were assessed for the probability score of functional regulatory regions in RegulomeDB [56] and ENCODE data [57,58] in UCSC genome browser.

## References

1. Fantes JA, Boland E, Ramsay J, Donnai D, Splitt M, et al. (2008) FISH mapping of de novo apparently balanced chromosome rearrangements identifies characteristics associated with phenotypic abnormality. Am J Hum Genet 82: 916–926.

2. Warburton D (1991) De novo balanced chromosome rearrangements and extra marker chromosomes identified at prenatal diagnosis: clinical significance and distribution of breakpoints. Am J Hum Genet 49: 995–1013.

## Validations of Expected Breakpoints by PCR

Primers were designed by Primer3 program, and the amplicons spanning the breakpoint were predicted by dPET clusters according to human genome assembly NCBI Build 36. PCR was carried out with JumpStart REDAccuTaq LA polymerase (Sigma Aldrich Inc., St. Louis, MO) in a 50 µl reaction volume and with 500 ng of genomic DNA as a template. The following program was used: 1) Initial denaturation at 96°C for 30 sec, 2) 40 cycles of 15 sec at 94°, 30 sec at 58°C, 10 min at 68°, 3) 68°C for 10 min. PCR products showing single bands were purified by Gel Extraction Kit (Qiagen) and used as templates for sequencing in both directions by Sanger sequencing. The sequences of junction fragments were aligned to the human genome reference sequence using Blat [59].

## Expression Analysis

Total RNA was extracted from patient's EBV-LCL or fibroblasts using RNAeasy kit (Qiagen). Reverse transcription of 2 µg RNA derived from patients, 2 lymphoblast controls, 2 fibroblast controls and commercially available human tissue panel RNA (Clontech) was performed in 20 µl of SuperScript III Reverse Transcriptase reaction buffer (Life Technologies) using random hexamer primers. Quantitative real-time PCR (qRT-PCR) reactions were performed in the ABI PRISM 7500 HT system (Life Technologies) with five-fold dilution of cDNA, 200 nM of each primer using the SybrGreen PCR Master Mix (Life Technologies). Data were analyzed using 2^ΔΔCt method, and normalized against control sample with human ACTB. Each measurement was performed in triplicate. The controls used in this study are derived from EBV-LCL or skin fibroblast cells of 4 normal individuals.

## Supporting Information

**Information S1** The file *SupportingInformation_S1.pdf* contains additional information to the manuscript. It consists of 9 pages, 1 Figure and 7 tables.
(DOCX)

## Acknowledgments

We thank the patients and family members for their participation. We thank Dr. Stuart Tangye for Western Blot analysis in patient CD8. This study makes use of data generated by the DECIPHER Consortium. A full list of centres who contributed to the generation of the data is available from http://decipher.sanger.ac.uk and via email from decipher@sanger.ac.uk. The sequencing data from this study have been submitted to NCBI Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo) under accession no. GSE51430.

## Author Contributions

Conceived and designed the experiments: AMH VC SD. Performed the experiments: KHU EGYC LSS ASMT ZZ LSS SY IP. Analyzed the data: KHU AMH CWHL PJC CCS PNA SLR. Wrote the paper: KHU AMH IA VC SD SKHT. Provided sequencing platform and facility: XR YR AMH ETL. Cpatient samples: GP FC MW AK GH AM OP SB RVJ. Provided analysis tools: KWKS.

3. MacGregor DJ, Imrie S, Tolmie JL (1989) Outcome of de novo balanced translocations ascertained prenatally. J Med Genet 26: 590–591.

4. Vandeweyer G, Kooy RF (2009) Balanced translocations in mental retardation. Hum Genet 126: 133–147.

5. Hochstenbach R, van Binsbergen E, Engelen J, Nieuwint A, Polstra A, et al. (2009) Array analysis and karyotyping: workflow consequences based on a retrospective study of 36,325 patients with idiopathic developmental delay in the Netherlands. Eur J Med Genet 52: 161–169.
6. Rauch A, Hoyer J, Guth S, Zweier C, Kraus C, et al. (2006) Diagnostic yield of various genetic approaches in patients with unexplained developmental delay or mental retardation. Am J Med Genet A 140: 2063-2074.

7. Ray PN, Belfall B, Duff C, Logan C, Kean V, et al. (1985) Cloning of the breakpoint of an X:21 translocation associated with Duchenne muscular dystrophy. Nature 318: 672-675.

8. Muir WJ, Gosden CM, Brookes AJ, Fantes J, Evans KL, et al. (1995) Direct microdissection and microcloning of a translocation breakpoint region, t(1;11)(q42.2;q21), associated with schizophrenia. Cytogenet Cell Genet 70: 35-40.

9. Millar JK, Wilson-Annan JC, Anderson S, Christie S, Taylor MS, et al. (2000) Disruption of two novel genes by a translocation co-segregating with schizophrenia. Hum Mol Genet 9: 1415-1423.

10. Chelly J, Tumer Z, Tonnesen T, Petterson A, Ishikawa-Brush Y, et al. (1993) Isolation of a candidate gene for Menkes disease that encodes a potential heavy metal binding protein. Nat Genet 3: 14-19.

11. Verga V, Hall BK, Wang SR, Johnson S, Higgins JV, et al. (1991) Localization of the translocation breakpoint in a female with Menkes syndrome to Xq13.2-q13.3 proximal to PGK-1. Am J Hum Genet 48: 1133-1138.

12. Miller DT, Adam MP, Aradhya S, Biesecker LG, Brothman AR, et al. (2010) Consensus statement: chromosomal microarray is a first-tier clinical diagnostic test for individuals with developmental disabilities or congenital anomalies. Am J Hum Genet 86: 749-764.

13. Rauch A, Ruschendorf F, Huang J, Trautmann U, Becker C, et al. (2004) Molecular karyotyping using an SNP array for genomewide genotyping. J Med Genet 41: 916-922.

14. Shaffer LG (2005) American College of Medical Genetics guideline on the cytogenetic evaluation of the individual with developmental delay or mental retardation. Genet Med 7: 650-654.

15. Le Scouarnec S, Gribble SM (2012) Characterising chromosome rearrangements: recent technical advances in molecular cytogenetics. Heredity (Edinb) 108: 75-85.

16. Savage MS, Mourad MJ, Wapner RJ (2011) Evolving applications of microarray analysis in prenatal diagnosis. Curr Opin Obstet Gynecol 23: 103-108.

17. Hillmer AM, Yao F, Inaki K, Lee WH, Ariyaratne PN, et al. (2011) Comprehensive long-span paired-end-tag mapping reveals characteristic patterns of structural variations in epithelial cancer genomes. Genome Res 21: 665-675.

18. Inaki K, Hillmer AM, Ukil L, Yao F, Woo XY, et al. (2011) Transcriptional consequences of genomic structural aberrations in breast cancer. Genome Res 21: 676-687.

19. Yao F, Ariyaratne PN, Hillmer AM, Lee WH, Li G, et al. (2012) Long span DNA paired-end-tag (DNA-PET) sequencing strategy for the interrogation of genomic structural mutations and fusion-point-guided reconstruction of amplicons. PLoS One 7: e46152.

20. Alam MM, Chattopadhyaya M, Chakrabarti S (2012) On the origin of large two-photon activity of DANS molecule. J Phys Chem A 116: 11034-11040.

21. Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, et al. (2004) Detection of large-scale variation in the human genome. Nat Genet 36: 949-951.

22. Mills RE, Walter K, Stewart C, Handsaker RE, Chen K, et al. (2011) Mapping copy number variation by population-scale genome sequencing. Nature 470: 59-65.

23. Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, et al. (2012) An integrated map of genetic variation from 1,092 human genomes. Nature 491: 56-65.

24. Kidd JM, Cooper GM, Donahue WF, Hayden HS, Sampas N, et al. (2008) Mapping and sequencing of structural variation from eight human genomes. Nature 453: 56-64.

25. Korbel JO, Urban AE, Affourtit JP, Godwin B, Grubert F, et al. (2007) Paired-end mapping reveals extensive structural variation in the human genome. Science 318: 420-426.

26. Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, et al. (2010) Origins and functional impact of copy number variation in the human genome. Nature 464: 704-712.

27. Corpas M, Bragin E, Clayton S, Bevan P, Firth HV (2012) Interpretation of genomic copy number variants using DECIPHER. Curr Protoc Hum Genet Chapter 8: Unit 8 14.

28. Cooper GM, Coe BP, Girirajan S, Rosenfeld JA, Vu TH, et al. (2011) A copy number variation morbidity map of developmental delay. Nat Genet 43: 838-846.

29. Swaminathan GJ, Bragin E, Chatzimichali EA, Corpas M, Bevan AP, et al. (2012) DECIPHER: web-based, community resource for clinical interpretation of rare variants in developmental disorders. Hum Mol Genet 21: R37-44.

30. Offermanns S, Hashimoto K, Watanabe M, Sun W, Kurihara H, et al. (1997) Impaired motor coordination and persistent multiple climbing fiber innervation of cerebellar Purkinje cells in mice lacking Galphaq. Proc Natl Acad Sci U S A 94: 14089-14094.

31. Takemoto M, Hattori Y, Zhao H, Sato H, Tamada A, et al. (2011) Laminar and areal expression of unc5d and its role in cortical cell survival. Cereb Cortex 21: 1925-1934.

32. Zhu Y, Li Y, Haraguchi S, Yu M, Ohira M, et al. (2013) Dependence receptor UNC5D mediates nerve growth factor depletion-induced neuroblastoma regression. J Clin Invest 123: 2935-2947.

33. Ben-Zur T, Feige E, Motro B, Wides R (2000) The mammalian Odz gene family: homologs of a Drosophila pair-rule gene with expression implying distinct yet overlapping developmental roles. Dev Biol 217: 107-120.

34. Ben-Zur T, Wides R (1999) Mapping homologs of Drosophila odd Oz (odz): Doc4/Odz4 to mouse chromosome 7, Odz1 to mouse chromosome 11; and ODZ3 to human chromosome Xq25. Genomics 58: 102-103.

35. Girirajan S, Rosenfeld JA, Coe BP, Parikh S, Friedman N, et al. (2012) Phenotypic heterogeneity of genomic disorders and rare copy-number variants. N Engl J Med 367: 1321-1331.

36. de Ligt J, Willemse MH, van Bon BW, Kleefstra T, Yntema HG, et al. (2012) Diagnostic exome sequencing in persons with severe intellectual disability. N Engl J Med 367: 1921-1929.

37. Talkowski ME, Rosenfeld JA, Blumenthal I, Pillalamarri V, Chiang C, et al. (2012) Sequencing chromosomal abnormalities reveals neurodevelopmental loci that confer risk across diagnostic boundaries. Cell 149: 525-537.

38. Talkowski ME, Ernst C, Heilbut A, Chiang C, Hanscom C, et al. (2011) Next-generation sequencing strategies enable routine detection of balanced chromosome rearrangements for clinical diagnostics and genetic research. Am J Hum Genet 88: 469-481.

39. Millar JK, Christie S, Anderson S, Lawson D, Hsiao-Wei Loh D, et al. (2001) Genomic structure and localisation within a linkage hotspot of Disrupted In Schizophrenia 1, a gene disrupted by a translocation segregating with schizophrenia. Mol Psychiatry 6: 173-178.

40. Chen W, Ullmann R, Langnick C, Menzel C, Wotschofsky Z, et al. (2010) Breakpoint analysis of balanced chromosome rearrangements by next-generation paired-end sequencing. Eur J Hum Genet 18: 539-543.

41. Chiang C, Jacobsen JC, Ernst C, Hanscom C, Heilbut A, et al. (2012) Complex reorganization and predominant non-homologous repair following chromosomal breakage in karyotypically balanced germline rearrangements and transgenic integration. Nat Genet 44: 390-397, S391.

42. Kloosterman WP, Guryev V, van Roosmalen M, Duran KJ, de Bruijn E, et al. (2011) Chromothripsis as a mechanism driving complex de novo structural rearrangements in the germline. Hum Mol Genet 20: 1916-1924.

43. Kloosterman WP, Tavakoli-Yaraki M, van Roosmalen MJ, van Binsbergen E, Renkens I, et al. (2012) Constitutional chromothripsis rearrangements involve clustered double-stranded DNA breaks and nonhomologous repair mechanisms. Cell Rep 1: 648-655.

44. Schluth-Boland C, Labalme A, Cordier MP, Till M, Nadeau G, et al. (2013) Breakpoint mapping by next generation sequencing reveals causative gene disruption in patients carrying apparently balanced chromosome rearrangements with intellectual deficiency and/or congenital malformations. J Med Genet 50: 144-150.

45. Ng KP, Hillmer AM, Chuah CT, Juan WC, Ko TK, et al. (2012) A common BIM deletion polymorphism mediates intrinsic resistance and inferior responses to tyrosine kinase inhibitors in cancer. Nat Med.

46. Nagarajan N, Bertrand D, Hillmer AM, Zang ZJ, Yao F, et al. (2012) Whole-genome reconstruction and mutational signatures in gastric cancer. Genome Biol 13: R115.

47. van Heesch S, Kloosterman WP, Lansu N, Ruzius FP, Levandowsky E, et al. (2013) Improving mammalian genome scaffolding using large insert mate-pair next-generation sequencing. BMC Genomics 14: 257.

48. McVey M, Lee SE (2008) MMEJ repair of double-strand breaks (director's cut): deleted sequences and alternative endings. Trends Genet 24: 529-538.

49. Wang JC, Dang L, Fisker T (2010) Chromosome 6 between-arm intrachromosomal insertion with intrasegmental double inversion: a four-break model. Am J Med Genet A 152A: 209-211.

50. Kloosterman WP, Cuppen E (2013) Chromothripsis in congenital disorders and cancer: similarities and differences. Curr Opin Cell Biol.

51. Liu P, Erez A, Nagamani SC, Dhar SU, Kolodziejska KE, et al. (2011) Chromosome catastrophes involve replication mechanisms generating complex genomic rearrangements. Cell 146: 889-903.

52. Rigaud S, Fondaneche MC, Lambert N, Pasquier B, Mateo V, et al. (2006) XIAP deficiency in humans causes an X-linked lymphoproliferative syndrome. Nature 444: 110-114.

53. Rigaud S, Latour S (2007) [An X-linked lymphoproliferative syndrome (XLP) caused by mutations in the inhibitor-of-apoptosis gene XIAP]. Med Sci (Paris) 23: 235-237.

54. Lichter P, Tang CJ, Call K, Hermanson G, Evans GA, et al. (1990) High-resolution mapping of human chromosome 11 by in situ hybridization with cosmid clones. Science 247: 64-69.

55. Cherif D, Derre J, Berger R (1990) Fluorescence in situ hybridization on metaphase chromosomes with biotinylated probes. In situ hybridization, biotin labeling, cosmids, gene mapping, oncogene amplification. Nouv Rev Fr Hematol 32: 459-460.

56. Boyle AP, Hong EL, Hariharan M, Cheng Y, Schaub MA, et al. (2012) Annotation of functional variation in personal genomes using RegulomeDB. Genome Res 22: 1790-1797.

57. Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, et al. (2012) The accessible chromatin landscape of the human genome. Nature 489: 75-82.

58. Bernstein BE, Birney E, Dunham I, Green ED, Gunter C, et al. (2012) An integrated encyclopedia of DNA elements in the human genome. Nature 489: 57-74.

59. Kent WJ (2002) BLAT-the BLAST-like alignment tool. Genome Res 12: 656-664.