MLPA analysis in a cohort of patients with autism

Transcription

1 Peixoto et al. Molecular Cytogenetics (2017) 10:2 DOI /s z RESEARCH MLPA analysis in a cohort of patients with autism Sara Peixoto 1,2,3*, Joana B. Melo 1,4,5, José Ferrão 1, Luís M. Pires 1, Nuno Lavoura 1, Marta Pinto 1, Guiomar Oliveira 2,6 and Isabel M. Carreira 1,4,5* Open Access Abstract Background: is a global neurodevelopmental disorder which generally manifests during the first 2 years and continues throughout life, with a range of symptomatic variations. Epidemiological studies show an important role of genetic factors in autism and several susceptible regions and genes have been identified. The aim of our study was to validate a cost-effective set of commercial Multiplex Ligation dependent Probe Amplification (MLPA) and methylation specific multiplex ligation dependent probe amplification (MS-MLPA) test in autistic children refered by the neurodevelopmental center and autism unit of a Paediatric Hospital. Results: In this study 150 unrelated children with autism spectrum disorders were analysed for copy number variation in specific regions of chromosomes 15, 16 and 22, using MLPA. All the patients had been previously studied by conventional karyotype and fluorescence in situ hybridization (FISH) analysis for 15(q11.2q13) and, with these techniques, four alterations were identified. The MLPA technique confirmed these four and identified further six alterations by the combined application of the two different panels. Conclusions: Our data show that MLPA is a cost effective straightforward and rapid method for detection of imbalances in a clinically characterized population with autism. It contributes to strengthen the relationship between genotype and phenotype of children with autism, showing the clinical difference between deletions and duplications. Keywords:, spectrum disorders, Copy number variants, Genotype, Multiplex Ligation-dependent Probe Amplification (MLPA), Methylation Specific Multiplex Ligation-dependent Probe Amplification (MS-MLPA), Phenotype Background is a global neurodevelopment disorder which, in most cases, manifests during the first 2 years and prolongs throughout life, having clinical variances with aging. It belongs to a large family of disorders - autism spectrum disorder (ASD), clinically characterized by difficulties with communication and social interaction, verbal language deficiencies and by repetitive and stereotype behaviour [1]. Epidemiological studies revealed that the number of children with ASD has been increasing throughout the world. In Portugal, a published study in 2007, suggested a prevalence of 9,2/ cases, in the main land, and 15,6/ in the Azores Islands [2]. Several epidemiologic * Correspondence: saraccpeixoto@gmail.com; icarreira@fmed.uc.pt Equal contributors 1 Cytogenetics and Genomics Laboratory, Faculty of Medicine, University of Coimbra, Coimbra, Portugal Full list of author information is available at the end of the article studies identified the importance of genetic factors in autism, specifically the existence of a higher rate of similarities in monozygotic twins (60 to 90%) in contrast to only 3 to 10% in dizygotic twins [3 6]. Structural chromosomic imbalances - copy number variants (CNVs) seams to be a key player in the disorder and a risk factor for autism especially in the sporadic forms [7 11]. CNVs associated with autism may be inherited or de novo, affecting preferentially the chromosomes regions: 1q21, 2p16.3, 3p25-26, 7q36.2, 15q11-13, 15q24, 16p11.2, 16p13.11, 17q12 and 22q11.2 [12 14]. A recent genome wide copy number variation analysis projected between 156 and 280 genomic intervals contributing to autism. Exome sequencing of over 900 individuals provided an estimate of nearly contributing genes [15 18]. The most consistently reported submicroscopic chromosome abnormalities detected by chromosomal microarray (10 20%) are recurrent CNVs at 16p11.2, 15q11-13 and The Author(s) Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License ( which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

2 Peixoto et al. Molecular Cytogenetics (2017) 10:2 Page 2 of 10 22q11.2 [12, 19, 20]. Submicroscopic rearrangements as deletions or duplications in 15q11-13 region, and especially proximal 15q duplications containing the critical regions of Prader-Willi and Angelman Syndrome (PWS/AS), have been reported in various patients with autism representing the majority of chromosome alterations described in this population [12, 21]. Microdeletions of the short arm of chromosome 16p11.2 were also identified in up to 1% of patients with autism [22 25]. The microduplication of this region was also observed in a similar percentage of individuals, however the association with autism is less convincing due to the increased frequency observed in the control groups [25, 26]. Deletions and duplications in the region 22q13.3 also appear to be risk factors for ASD. Among the three genes (ACR, RABL2B, SHANK3) in 22q subtelomeric region, SHANK3 is the candidate gene for neurobehavioral symptoms observed in affected individuals with 22q13 deletions. This gene has a mutation frequency of 0.5 to 1% in individuals with ASD [27 29]. It is well accepted that array comparative genomic hybridization (CGH) should be the first genetic test to be offered for detection of genomic imbalances in patients with intellectual disability and ASD [30, 31], however access to these costly methods and techniques can be difficult especially given the current economic and financial context in many countries and consequently the pressure from the hospital administrations to contain expense. Clinicians often are confronted with this economical and financial issues and have difficulties to respond to families request to find the cause of these children pathology. One of the aims of this study was to evaluate the detection rate of MLPA technique that could allow a rapid and cost-saving response in a large consultation of children with ASD. We analysed the presence of CNVs in the most common described chromosomic regions associated with ASD (15q11-q13, 16p11.2 and 22q13) in 150 children using MLPA and MS-MLPA techniques. A clinical assessment of each case was also done to allow genotype-phenotype correlations in the cases with genetic alterations. Methods Clinical assessment For this study, a pilot population of 150 children was analysed. They were clinically diagnosed with autism by the Unit of Neurodevelopment and of the Paediatric Hospital, Coimbra Hospital and University Centre, Portugal. The diagnosis was based on a clinical observation by a multidisciplinary team coordinated by a neurodevelopmental paediatrician. ASD diagnosis was assigned on the basis of the gold standard instruments: parental or caregiver interview ( Diagnostic Interview-Revised, ADI-R [32]), direct structured proband assessment ( Diagnostic Observation Schedule, ADOS [33]), both history and observation for rating (Childhood Rating Scale, CARS [34]) and clinical examination performed by an experienced neurodevelopmental Paediatrician. The latter allows the classification of the degree of autism into mild to moderate (score between 30 and 37) and severe (38 60). The current diagnostic criteria for autism were revised according to the Diagnostic and Statistical Manual of Mental Disorders fifth edition, DSM-5 [1]. It was considered a diagnosis of autism (this term being used synonymously to ASD) any case where ADI-R and ADOS presented as positive scores and all patients met the criteria for ASD from the DSM-5. An EDTA blood sample from each child was collected for the genetic evaluation by MLPA and MS-MLPA for chromosomes 15q11-13, 16p11 and 22q13. A more detailed clinical evaluation was done for the ten children that presented alterations in the MLPA study, proceeding to the analysis of the relative clinical relevant data to the current clinical history of autism, accordingly to Table 1. Whenever possible, laboratory studies were done on the parents of the children with chromosome alterations. Conventional cytogenetic and fluorescence in situ hybridization All 150 children participating in this study, had been previously evaluated by high resolution conventional Table 1 Clinical data relevant to the clinical history of the cohort Personal history: pre and perinatal Parturition type Gestational age Apgar index Somatometry birth (weight/height/head circumference) Personal history: Acquisition of neurodevelopment Walking age First words First sentences Global developmental quotient (GDQ) - Griffiths scale (between 2 and 6 years of age) Global intelligence quotient (GIQ) - WISC-III (between 6 and 16 years of age) Pathological personal history Visual or auditory deficits Epilepsy (two or more critical episodes in apyrexia) Family history Physical exam Dysmorphisms and signs of neurocutaneous syndromes Collection of anthropometric measurements (actual growing) Classical neurologic exam

3 Peixoto et al. Molecular Cytogenetics (2017) 10:2 Page 3 of 10 cytogenetic using standard GTG-banding on prometaphases obtained from 72 h PHA stimulated peripheral blood lymphocyte cultures according to standard procedures [35]. FISH for the critical region 15q11.2 was performed using two commercially available DNA probes: LSI SNRPN/CEP15 (locus D15Z1)/LSI PML and LSI D15S10/ CEP15 (locus D15Z1)/LSI PML (Vysis, Chicago, IL) according to standard procedures and manufacturer s instructions [35]. Multiplex Ligation-dependent Probe Amplification (MLPA) MLPA (P343-B1) and MS-MLPA (ME028) (Fig. 1a) probe panels were applied as described in the protocols of the manufacturer (MRC Holland). The P343-B1 panel, applied to all cases, has 54 MLPA probes for the three regions 15q11-13, 16p11 and 22q13, implicated in autism. The 12 most relevant genes studied in chromosome 15 were: SNRPN-HB2-85, UBE3A, ATP10A, GABRB3, OCA2, APBA2, NDNL2, TJP1, TRPM1, KLF13, CHRNA7, SCG5). The MS-MLPA ME028 PWS/AS panel, was applied to the patients without alterations detected with the P343-B1 panel. This panel has 25 specific probes for the critical regions of Prader-Willi and Angelman syndromes, as well as five probes to assess the state of methylation, allowing the study of more proximal genes that are not included in P343-B1 panel, like NIPA1 and TUBGCP5. In the region of microdeletion 16p11 the nine studied genes were LAT, SPN, MAS, MVP, SEZ6L2, HIRIP3, DOC2A, MAPK3, CD2BP2. The gene SHANK3 in 22q13 was the one studied in chromosome 22. The products of amplification were identified and quantified by capillary electrophoresis in an ABI 3130 genetic analyser (Applied Biosystems, Japan) and the results were analysed using the GeneMapper v4.1 software (Applied Biosystems, Foster City, USA) and Coffalyser (MRC-Holland, Amsterdam, Holland). The linear ratios of deletion and duplication were fixed at 0.7 e 1.3 respectively. Results The study of the conventional karyotype, together with FISH analysis for 15(q11.2q13) detected duplication of this proximal region in three cases (I, IV, V) and triplication in one case (VI) (Table 2). The application of the MLPA Panel P343 confirmed two of the cases identified by cytogenetics (I, VI), redefined to a triplication a duplication previously diagnosed by FISH (case V) and identified alterations in three more cases - a duplication (15q11.2-q13.1) (case II) and a deletion (15q13.2-q13.3) (case IX) of the critical region of the chromosome 15 and a duplication in 22q13.33 (case X) (Table 2 and Fig. 1b). The use of Panel ME028 confirmed one of the alterations, diagnosed by FISH as a mosaic, to be a maternal duplication of the proximal region of 15q11.2 (case IV) and identified three more cases, with normal result after application of MLPA Panel P343: two microdeletions (cases VII and VIII) and one duplication (case III) of the most proximal region of 15q11.2 (NIPA1, TUBGCP5) (Table 2 and Fig. 1b). No alterations were observed in these 150 patients for the critical region 16p11.2. The main clinical data from these ten patients with structural chromosome alterations are described in Table 3. A B C Fig. 1 a Schematic representation of the studied genes by the two panels (MLPA panel ME028 and panel P343) of probes and of breakpoints on chromosome 15. Arrows indicate the genes identified by each panel (P343 dark arrows; MEO28 blue arrows), some of the genes (two arrows) are common to both panels; b Schematic representation of the extent of alterations detected in the region 15q11-15q13 by the two MLPA panels and the involved genes in the studied cases, with duplications in green boxes, triplications represented in yellow and deletions in red; c Schematic representation of the alteration detected in the region 22q13.33 by the MLPA P343 panel. BP-break point; C Case; CEN-centromere; Del-deletion; Dup-duplication; mos-mosaic; TEL-telomere; Trip-triplication

4 Peixoto et al. Molecular Cytogenetics (2017) 10:2 Page 4 of 10 Table 2 Detailed clinical features of all ten probands with chromosomal abnormalities Case/ Gender Age (Y) I/M 17 ADI-R/ADOS + CARS: 44,5 Severe Diagnosis Cognition Family History Karyotype and FISH (ISCN) Severe/Profound Intellectual Disability (QIG:29) a Brother has autism 46,XY.ish dup(15)(q11.2 q11.2)(snrpn++,ube3a++) P343 MLPA panel (ISCN) rsa15q11.2-q13.1(snrpn, UBE3A,ATP10A,GABRB3, OCA2)x3 ME028 MLPA panel (ISCN) Abnormality origin - Mother normal Father unavailable II/M 6 ADI-R/ADOS + Severe Mild Intellectual Adopted Biological Mother Disability (QDG:54) b has Depressive Disorder 46,XY rsa 15q11.2-q13.1(SNRPN, UBE3A,ATP10A,GABRB3, OCA2)x3 - unavailable III/M 6 ADI-R/ADOS + CARS: 35 Mild No Intellectual Disability (QDG: 102) b Father has drug addiction Brother has social disability 46,XY rsa(p343)x2 rsa 15q11.2(NIPA1, TUBGCP5)x3mat Mother IV/M 12 ADI-R/ADOS + CARS: 30 Mild Mild Intellectual Paternal grandfather has Disabilty (QDG:58) b schizophrenia Maternal grandmother is alcoholic Paternal Uncle has drug addiction mos 47,XY,+r.ish r(15)(d15z4 +,SNRPN-,UBE3A-)[69]/r(15) (D15Z4++,SNRPN-,UBE3A-) [2]/46,XY[105] rsa(p343)x2 rsa 15q11.2(NIPA1, TUBGCP5,MKRN3, MAGEL2,NDN)X3mat Mother V/F 9 ADI-R/ADOS + CARS: 37 Moderate Mild/Moderate Intellectual Disability (QDG:39) b Irrelevant 46,XX,dup(15)(q11.2q11.2) (SNRPN++,D15S10++) rsa 15q11.2-q13.1(SNRPN, UBE3A,ATP10A,GABRB3, OCA2)x4 - unavailable VI/F 10 ADI-R/ADOS + CARS: 32 Mild Severe/profound Intellectual Disability (QDG:29) b Irrelevant 46,XX,add(15)(q11.2).ish trp(15)(q11.2)(snrpn, UBE3A)x4 rsa 15q11.2-q13.1(SNRPN, UBE3A,ATP10A,GABRB3, OCA2,APBA2,NDNL2,TJP1) x4dn - de novo VII/M 7 ADI-R/ADOS + Mild No intellectual Irrelevant 46,XY rsa(p343)x2 rsa 15q11.2(NIPA1, Disability (QDG: 97) b TUBGCP5)x1 unavailable VIII/M 8 ADI-R/ADOS + CARS: 30 Mild Intellectual Disability Mother has learning (QDG:71) b difficulties 46,XY rsa(p343)x2 rsa 15q11.2(NIPA1, TUBGCP5)x1mat Mother IX/M 14 ADI-R/ADOS + CARS: 30 Mild Mild Intellectual Father has chromosomal Disability (QDG:61) b abnormalities but no clinical issues Paternal grandmother has Depressive Disorder 46,XY rsa 15q13.2-q13.3(MTMR15, TRPM, KLF13,CHRNA7)x1pat - Father X/M 5 ADI-R/ADOS + CARS: 31 Mild No intellectual Disability (QDG:99) Irrelevant 46,XY rsa 22q13.33(SHANK3)x3mat - Mother a WISC III; M-male, F-female b Griffiths ADI-R Diagnostic Interview Revised ADOS Diagnostic Observation Schedule ISCN International system for human cytogenetic nomenclature, 2016

5 Peixoto et al. Molecular Cytogenetics (2017) 10:2 Page 5 of 10 Table 3 Prenatal and perinatal History and Neurodevelopment Case Prenatal History Actual Growing I Part/GA: Forceps/38w BW: 3350 g (P63); BL:52 cm (P90); BHC:35 cm (P72) AI :10 GPDD (1st year) March (months) 1st Words (months) 1st sentences (months) W: P75 H: P75 HC: >P50 No a Doesn t build sentences II Part/GA: Eutocic/? BW: 3380 g (P?); AI:10 W: P95 H: P50 HC: P95 Yes No III Part/GA: Ventouse/39w BW: 3450 g (P55); BL:50 cm W: P90 H: P75 HC: >P50 No a No (P47); BHC:35.5 cm (P69); AI :10 IV Part/GA: CS/ 40w BW: 3460 g (P41); BL:50.5 cm W: P75/90 H: P90 HC: P90 Yes No (P39); BHC:35 cm (P45); AI:10 V Part/GA: Ventouse/39w BW: 2485 g (P4); BL:47 cm W: P25 H: P25 HC: P50 Yes No (P8); BHC:32.5 cm (P11) AI :10 VI Part/GA: Eutocic/? BW: 2840 g (P?); BL:50 cm (P?); IA:10 W: >P75 H: P50 HC: P50 Yes Doesn t build Yes sentences VII Part/GA: Eutocic/35w BW: 2510 g(p46); BL:45 cm (P31); W: P50 H: P50 HC: P75 Yes No BHC:33 cm (P72) AI:10 VIII Part/GA: CS/ 38w High Risk in Pregnancy due to previous abortions (2) BW: 3470 g (P72); BL:49 cm (P47); BHC:36 cm (P88); AI:10 W: P95 H: P95 HC: > > P50 Yes No IX Unknown W: P50 H: P50 HC: >P95 Yes No X Part: CS/ 41w BW: 4610 g (P94); BL:51 cm (P35); BHC:37 cm (P76); AI:10 W: P95 H: P90 HC: >P95 No a No Epilepsy a only identified at the age of 2 years old - followed by regression Pre and perinatal history: AI Apgar Index at 5 min, BHC Birth head circumference, BL Birth length, BW Birth Weight, CS Caesarean Section, GA Gestational Age, Part Parturition type, P percentile according Fenton growth chart Actual Growing: HC head circumference, H Height, W Weight GPDD Global psychomotor developmental delay, w pregnancy weeks, g grams, cm centimeters No The study of the patient s parents was possible in only six of the ten individuals with alterations: in four of them, the alteration is maternally inherited (cases III, IV, VIII and X), one has a paternal origin (case IX) and in case VI the alteration is de novo. For the remaining four cases, it was not possible to ascertain the parental origin because either child was adopted or because some parents rejected to be studied. Cases I and II, are carriers of the same duplication on chromosome 15 [15q11.2q13.1(SNRPN,UBE3A,ATP10A, GABRB3,OCA2) x3] (Fig. 1b), both have severe autism but with different levels of cognition. Case II is an adopted boy with a global psychomotor developmental delay from the first year of age and currently with mild intellectual disability. Cases III and IV, evaluated by the second MLPA probe panel have more proximal duplications of chromosome 15 (Fig. 1b). Both children have mild autism. Case III which has the more proximal duplication shows a normal global intelligence quotient (GIQ) and a normal motor and language development. Case IV has a mild intellectual deficit. Cases V and VI (Fig. 1b), both presented mild to moderate autism but with distinct levels of intellectual disability, varying from mild in case V to severe in case VI, which is coincident with the extent of the triplicated region that includes a greater number of genes. In both patients, strabismus is present as well as motor development delay, not found in previously described patients, carriers of duplications. This strongly suggests agenedosageeffect. Cases VII and VIII (Fig. 1b) both have microdeletions of the proximal genes, both have similar levels of neurodevelopment features, presenting mild autism and mild (case VIII) or absent (case VII) intellectual disability. In case VIII, the alteration is of maternal origin and the progenitor is reported as having learning difficulties. In case IX (Fig. 1b), the deleted region is distal to the previous two cases reported and involves a greater number of genes, however the clinical presentation is similar between the three cases (VII, VIII and IX). In this case the deletion was inherited from the father, but was not possible to ascertain his phenotype. Lastly, from the 150 individual evaluated there was one case (X) identified with duplication of the SHANK3 gene on 22q13.33, associated with mild autism and average intellectual quotient and without other significant neurodevelopmental alterations except for clumsiness. The duplication was inherited from an asymptomatic mother. Concerning to the prenatal history, most cases did not reveal any anomalies, being born from gestations of term, without prenatal incidents. As for the early development markers, half of the children walked after 18 months, which suggests an early alteration in neurodevelopment. However, the growth with regards to

6 Peixoto et al. Molecular Cytogenetics (2017) 10:2 stature, weight and head circumference, did not reveal alterations (Table 3). Discussion Of the 150 children evaluated in this pilot study, only four (3%) had revealed alterations detected by conventional cytogenetics and FISH. This number increased after MLPA and MS-MLPA analysis with the identification of 6 more imbalances. Being a technique with higher resolution, MLPA also allowed the redefinition of a duplication, initially reported by conventional cytogenetics as a triplication (case V - Fig. 2) and allowed the identification of imbalances of six new cases (Table 2). Of these ten patients, only one presented duplication on the SHANK3 gene on 22q13.33, which is in accordance with the reported frequency in the population (1% of patients with ASD) [13, 28]. All the others showed imbalances on chromosome 15 (four duplications, two triplications and three deletions) [9, 36, 37]. No cases of deletion nor duplication on 16p11.2 were found in our sample, possibly due to the relative low frequency of up to 1% mentioned in other studies. Walsh and Page 6 of 10 Bracken in 2011 made a review of literature and a meta-analysis of 3613 ASD patients from seven studies redefining a prevalence of 16p11.2 microdeletion CNVs in 0.5% ( %, 95% CI) and 16p11.2 microdeletion CNVs in 0.28% ( %, 95% CI) not in disagreement with our results [25, 38]. CNVs in 15q11-15q13 region In cases I, II, III and IV we found duplications of specific regions of chromosome 15 that involve several break points (Fig. 1b). Two of the patients (case I and II) presented a duplication that affects BP2-BP3 (Fig. 1b). This alteration involves the critical region for Prader-Willi and Angelman Syndrome which is subjected to genomic imprinting and referred as the imbalance more frequently found in individuals with ASD [37, 39]. Numerous studies support that the majority of cases with this alteration are associated to maternal transmission or appear de novo, while the paternal inheritance leads to a normal phenotype [37, 40, 41]. It was not possible to ascertain the parental origin in these two patients (cases I and II), however the family history is Fig. 2 FISH and MLPA results for patient V. a Metaphase hybridization showing a gain in 15q11.2 (D15S10 probe in red) interpreted as a duplication. b Metaphase hybridization showing a gain in 15q11.2 (SNRPN probe in red) interpreted as a duplication. c Overview of MLPA P343 result using CoffalyserV7 software (MRC Holland, Amsterdam, Netherland) revealing the presence of 4 copies and not 3 of region 15q11.2q12 in patient V compatible with a triplication and not a duplication

7 Peixoto et al. Molecular Cytogenetics (2017) 10:2 Page 7 of 10 suggestive of neurodevelopmental delays, with a brother with autism (case I) and neuropsychiatric disorders, with mother with depressive disorder (case II), indicating a possible inheritance of the mentioned alterations (Table 2). These two cases (I and II) with the same genetic imbalances, both have severe autism although with different intellectual levels. According to Bolton PF et al. (2001) and Depienne C et al. (2009) the patients with duplication of the critical region (BP2-BP3), as occurs in these two cases, have frequently an abnormal phenotype that includes development delay, especially in speech and language, different degrees of intellectual disability, motor coordination difficulties and mild dysmorphisms (which are absent in our cases) [34, 40]. Case III, is a case of a more proximal duplication (between BP1-BP2) (Fig. 1b). The patient has mild autism, without intellectual impairment and with normal motor and language development. The common phenotype found in this alteration is highly variable, even between members of the same family [40, 42]. It can vary between autism, other neurodevelopmental disorders, mild socialization problems or learning difficulties, neuropsychiatric manifestations or even a normal phenotype, which is corroborated by other authors [37, 41]. The genetic alteration in this patient is also present in his mother who is phenotypically normal. One of the patient s brother, presents a serious socialization problem nevertheless he was not available for study. Case IV also presents a duplication of the proximal region involving more genes than the previous one (case III) (Fig. 1b). It corresponds to a mosaicism with three cell lines. This patient presents mild intellectual disability and mild autism, as in case III. In this case IV, a milder phenotype could be expected since this aberration is present in a mosaic state. In the blood lymphocytes, the normal cellular lineage is the most representative. The ring chromosome was inherited from the mother that is not a mosaic in the blood and does not present with any type of clinical alteration. Other tissues were not evaluated for the detection of the mosaic state. Regarding the two cases with triplication (V and VI), both have mild to moderate autism with distinct levels of intellectual disabilities which is, mild in case V but severe in case VI. The latter have a greater number of triplicated genes, being the distal breakpoint correspondent to BP4 (Fig. 1b). Genes APBA2, NDNL2 and TJP1 that are triplicated in case VI may be responsible for the more severe intellectual disability. Both patients had delay motor developmental skills (age of walking 24 months), which was not seen in the cases with duplication in our cohort, suggesting that a greater number of copies may be related to a more severe motor development delay. Both present vision deficit and one of them (case VI) also has epilepsy, which has been frequently associated with autism [43, 44]. Cases VII and VIII present microdeletion of proximal genes (NIPA1, TUBGCP5) and show similar levels of motor and cognitive development, both have mild autism with normal intelligence quotient and motor and language development delay, noted from the first year. These discrete phenotypes, as it happens in cases III and IV, that involve duplications of the same region, suggests that alterations involving this region, that not include the critical region of PWS and AS may not have any significant clinical effects. These alterations may also manifest as a milder form, like the mother in case VIII that, despite having the exact same deletion as her son, only had learning difficulties in her childhood which did not affect her adult independent life. Cases like this had been described repeatedly with significant phenotypic variability [39, 45]. Copy numbers losses and gains in this region have been identified with a frequency of 1% in the normal population, proving to be a challenge for genetic counselling [37, 39]. Case IX also corresponds to a microdeletion that is different from the previous ones since there are other breakpoints involved (BP4-BP5) (Fig. 1b) with a greater number of genes involved (MTMR15, TRPM1, KLF13, CHRNA7). The clinical condition is however similar to the two previous cases, contrary to what would be expected, based on the genetic differences. In this case, the deletion was inherited from an apparently normal father. This situation has already been reported in other cases, supporting the evidence of a higher frequency of deletions inherited from healthy parents, comparing to those that occur de novo [46, 47]. There are a great variety and heterogeneity of the phenotypic expression for this deletion of 1.6 Mb size. It ranges from normal phenotype (as it happens in the father), to intellectual disability with borderline intellectual quotient and autism (present in the patient), epilepsy, bipolar disorder, schizophrenia and other neuropsychiatric disorders [47, 48]. Depression that are frequently associated with this deletion, has been refered in the patient paternal grandmother but she was unavailable for study. In our sample, no cases of duplication of BP4-BP5 were found, which is in agreement with other reports [48 50]. If duplications involving this region cause phenotypical anomalies, its penetrance seems to be lower than it is in deletions which therefore may explain its minor frequency. The structural alterations presented in this study involving various regions of chromosome 15, have a significant role in the pathogenesis of other different conditions. They affect predominantly the brain function like in ASD and intellectual disability (that varies from mild to severe), epilepsy and neuropsychiatric disorders like schizophrenia and bipolar disease. This phenotypic variability suggests that other events may contribute for the manifestation of those conditions.

8 Peixoto et al. Molecular Cytogenetics (2017) 10:2 Page 8 of 10 Various factors have been proposed to explain this variability in phenotype. Besides from incomplete penetrance and the contribution of adjacent genes, the presence of a masked recessive mutation or functional polymorphism on one of the genes may be another factor, environmental and perinatal factors, maternal conditions [51, 52]. The possibility of different phenotypic expressions between diverse members of the family with a deletion, caused by defects of imprinting, would be another hypothesis to explain this variability, however it seems unlikely on the alterations that affect BP1-BP2 and BP4-BP5, as none of the genes in these regions have been described as imprinted (Genomic Imprinting Website: genes-by-species, September 2008). Another hypothesis is the possible capability that some individuals may have to overcome one or more deficiencies caused by haploinsufficiency of specific genes during pre or postnatal development [53]. This supports the fact that some carriers have normal development from birth, while others have learning difficulties during their childhood that cause no impact in the adult life. This is enforced by the fact that some parents who present the same alterations as their children, have no problem related to social integration. In addition to all these possible explanations, one cannot exclude the presence of other CNVs not detected with this approach, neither the environmental factors that may also play a vital role in the child s development. CNVs in SHANK3 at 22q13.33 region A single duplication was found on gene SHANK3 at 22q13.33 which is associated to mild autism with normal intellectual quotient. Microdeletions and point mutation involving SHANK3 region have been reported as cause of a spectrum of neuropsychiatric disorders including 22q13 deletion syndrome (also known as Phelan- Mcdermid syndrome) and ASD [28, 29]. In the same way, 22q13 duplications involving SHANK3 were reported in patients with Asperger syndrome, attention deficit-hyperactivity disorder (ADHD) or schizophrenia, suggesting that an overexpression of this gene could be also pathogenic [27]. In our sample, case X is a patient with a SHANK3 duplication that was inherited from a healthy mother which leads us to question of its pathogenicity. A situation of incomplete penetrance or expression variability cannot be ruled out in the same way as in 22q11.2 duplication syndrome [54], neither can be excluded the presence of other CNVs not detected with this approach. Conclusion In this study, it was also possible to establish some correlations between certain genotypes and corresponding phenotypes, despite of some variability. The study of autism is a challenge due to the great variety of behavioural and neurodevelopmental manifestations with variable severity and many co-morbidities. With the development of new technologies that allow an approach to the whole genome, certainly there will be a greater contribution for a better comprehension of its aetiology. It is vital that children with ASD are followedup by specialized professionals in the neurodevelopment area, so that clinical details can be well accurate to allow strong genotype-phenotype correlations. Our study indicates that MLPA can be a cost-effective method for detection of microdeletions and microduplications in ASD population. In patients with relevant phenotypic characteristics, this approach could replace other expensive and laborious techniques in clinical diagnosis. Additionally, in some cases karyotype and/or FISH analysis should be considered, mainly for the detection of supernumerary marker chromosomes. For example, the invdup(15) is a well documented cause of ASD with distinct phenotype and prognosis. It is well established that array-cgh should be the first-tier test for neurodevelopmental disorders, MLPA, although with some limitations, since it is a target method, can be considered as a test in large populations where the costs and economic pressures limits a more expensive method for diagnosis. Abbreviations ADHD: Attention deficit-hyperactivity disorder; ADIR: Diagnostic Interview Revised; ADOS: Diagnostic Observation Schedule; AS: Angelman Syndrome; ASD: spectrum disorder; BP: Break point; CARS: Childhood Rating Scale; CGH: Comparative genomic hybridization; CNV: Copy number variants; DSM-5: Diagnostic and Statistical Manual of Mental Disorders; EDTA: Ethylenediaminetetraacetic acid; FISH: Fluorescence in situ hybridization; GIQ: Global intelligence quotient; MLPA: Multiplex ligation dependent probe amplification; MS-MLPA: Methylation specific multiplex ligation dependent probe amplification; PWS: Prader Willi Syndrome Acknowledgements The authors wish to thank the patients and their families. Funding Not applicable for that section. Availability of data and materials The datasets generated and analysed during the current study are not publicly available due individual privacy could be compromised, but are available from the corresponding author on reasonable request and with permission of parents/guardians of all younger participants and Cytogenetics and Genomics Laboratory. Authors contributions IMC, GO have designed the study and drafted the manuscript; SP have clinically evaluated the patients and their families, provided clinical information for the interpretation of results, partial participation in MLPA analysis and interpretation and drafted the manuscript; JBM data acquisition and drafted the manuscript; LMP, JF have been involved in MLPA data acquisition and analysis and drafted the manuscript. MP, NL have performed conventional cytogenetics procedures, including karyotype and FISH analysis and helped to draft the manuscript. All authors read and approved the final manuscript.

9 Peixoto et al. Molecular Cytogenetics (2017) 10:2 Page 9 of 10 Competing interests The authors declare that they have no competing interests. Consent for publication Informed consent was obtained from the parents/guardians of all younger participants. Children and adolescents also gave oral informed consent. Ethics approval and consent to participate This study and all the procedures were reviewed and approved by the Ethics Commission of our Paediatric Hospital and was conducted in accordance with the declaration of Helsinki. Author details 1 Cytogenetics and Genomics Laboratory, Faculty of Medicine, University of Coimbra, Coimbra, Portugal. 2 Neurodevelopmental and Unit from Child Developmental Center and Centro de Investigação e Formação Clinica, Hospital Pediátrico, Centro Hospitalar e Universitário de Coimbra, Coimbra, Portugal. 3 Department of Paediatrics of the Centro Hospitalar de Trás-os-Montes e Alto Douro, EPE, Vila Real, Portugal. 4 CNC.IBILI, University of Coimbra, Coimbra, Portugal. 5 CIMAGO - Centro Investigação em Meio Ambiente, Genética e Oncobiologia, Faculty of Medicine, University of Coimbra, Coimbra, Portugal. 6 University Clinic of Pediatrics and Institute for Biomedical Imaging and Life Science, Faculty of Medicine, University of Coimbra, Coimbra, Portugal. Received: 16 September 2016 Accepted: 4 January 2017 References 1. Association, A.P. spectrum disorder. In: Diagnostic and Statistical Manual of Mental Disorders (DSM-5) p Oliveira G, et al. Epidemiology of autism spectrum disorder in Portugal: prevalence, clinical characterization, and medical conditions. Dev Med Child Neurol. 2007;49(10): Colvert E, et al. Heritability of Spectrum Disorder in a UK Population- Based Twin Sample. JAMA Psychiatry. 2015;72(5): Faras H, Al Ateeqi N, Tidmarsh L. spectrum disorders. Ann Saudi Med. 2010;30(4): Kumar RA, Christian SL. Genetics of autism spectrum disorders. Curr Neurol Neurosci Rep. 2009;9(3): Sandin S, et al. The familial risk of autism. JAMA. 2014;311(17): Abrahams BS, Geschwind DH. Advances in autism genetics: on the threshold of a new neurobiology. Nat Rev Genet. 2008;9(5): Genome Project, C, et al. Mapping autism risk loci using genetic linkage and chromosomal rearrangements. Nat Genet. 2007;39(3): Pinto D, et al. Convergence of genes and cellular pathways dysregulated in autism spectrum disorders. Am J Hum Genet. 2014;94(5): Prasad A, et al. A discovery resource of rare copy number variations in individuals with autism spectrum disorder. G3 (Bethesda). 2012;2(12): Sebat J, et al. Strong association of de novo copy number mutations with autism. Science. 2007;316(5823): Carter MT, Scherer SW. spectrum disorder in the genetics clinic: a review. Clin Genet. 2013;83(5): Freitag CM, et al. Genetics of autistic disorders: review and clinical implications. Eur Child Adolesc Psychiatry. 2010;19(3): Mefford HC, Batshaw ML, Hoffman EP. Genomics, intellectual disability, and autism. N Engl J Med. 2012;366(8): Iossifov I, et al. De novo gene disruptions in children on the autistic spectrum. Neuron. 2012;74(2): Neale BM, et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012;485(7397): O Roak BJ, et al. Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations. Nat Genet. 2011;43(6): Sanders SJ, et al. De novo mutations revealed by whole-exome sequencing are strongly associated with autism. Nature. 2012;485(7397): Carreira IM, et al. Copy number variants prioritization after array-cgh analysis - a cohort of 1000 patients. Mol Cytogenet. 2015;8: Moreira DP, et al. Investigation of 15q11-q13, 16p11.2 and 22q13 CNVs in autism spectrum disorder Brazilian individuals with and without epilepsy. PLoS One. 2014;9(9):e Schaaf CP, et al. Truncating mutations of MAGEL2 cause Prader-Willi phenotypes and autism. Nat Genet. 2013;45(11): Marshall CR, et al. Structural variation of chromosomes in autism spectrum disorder. Am J Hum Genet. 2008;82(2): Tabet AC, et al. multiplex family with 16p11.2p12.2 microduplication syndrome in monozygotic twins and distal 16p11.2 deletion in their brother. Eur J Hum Genet. 2012;20(5): Weiss LA, et al. Association between microdeletion and microduplication at 16p11.2 and autism. N Engl J Med. 2008;358(7): Walsh KM, Bracken MB. Copy number variation in the dosage-sensitive 16p11.2 interval accounts for only a small proportion of autism incidence: a systematic review and meta-analysis. Genet Med. 2011;13(5): Lord C, Rutter M, Le Couteur A. Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Dev Disord. 1994;24(5): Han K, et al. SHANK3 overexpression causes manic-like behaviour with unique pharmacogenetic properties. Nature. 2013;503(7474): Leblond CS, et al. Meta-analysis of SHANK Mutations in Spectrum Disorders: a gradient of severity in cognitive impairments. PLoS Genet. 2014;10(9):e Uchino S, Waga C. SHANK3 as an autism spectrum disorder-associated gene. Brain Dev. 2013;35(2): Friedman JM, et al. Oligonucleotide microarray analysis of genomic imbalance in children with mental retardation. Am J Hum Genet. 2006;79(3): Miller DT, et al. Consensus statement: chromosomal microarray is a first-tier clinical diagnostic test for individuals with developmental disabilities or congenital anomalies. Am J Hum Genet. 2010;86(5): Lord C, et al. diagnostic observation schedule: a standardized observation of communicative and social behavior. J Dev Disord. 1989;19(2): Rellini E, et al. Childhood Rating Scale (CARS) and Behavior Checklist (ABC) correspondence and conflicts with DSM-IV criteria in diagnosis of autism. J Dev Disord. 2004;34(6): Depienne C, et al. Screening for genomic rearrangements and methylation abnormalities of the 15q11-q13 region in autism spectrum disorders. Biol Psychiatry. 2009;66(4): Melo JB, et al. Chromosome 5 derived small supernumerary marker: towards a genotype/phenotype correlation of proximal chromosome 5 imbalances. J Appl Genet. 2011;52(2): Bijlsma EK, et al. Extending the phenotype of recurrent rearrangements of 16p11.2: deletions in mentally retarded patients without autism and in normal individuals. Eur J Med Genet. 2009;52(2 3): Pujana MA, et al. Human chromosome 15q11-q14 regions of rearrangements contain clusters of LCR15 duplicons. Eur J Hum Genet. 2002;10(1): Browne CE, et al. Inherited interstitial duplications of proximal 15q: genotype-phenotype correlations. Am J Hum Genet. 1997;61(6): Burnside RD, et al. Microdeletion/microduplication of proximal 15q11. 2 between BP1 and BP2: a susceptibility region for neurological dysfunction including developmental and language delay. Hum Genet. 2011;130(4): Bolton PF, et al. The phenotypic manifestations of interstitial duplications of proximal 15q with special reference to the autistic spectrum disorders. Am J Med Genet. 2001;105(8): Piard J, et al. Clinical and molecular characterization of a large family with an interstitial 15q11q13 duplication. Am J Med Genet A. 2010;152A(8): Combi R, et al. Clinical and genetic evaluation of a family showing both autism and epilepsy. Brain Res Bull. 2010;82(1 2): De Wolf V, et al. Genetic counseling for susceptibility loci and neurodevelopmental disorders: the del15q11.2 as an example. Am J Med Genet A. 2013;161A(11): Tuchman R, Hirtz D, Mamounas LA. NINDS epilepsy and autism spectrum disorders workshop report. Neurology. 2013;81(18): Cox DM, Butler MG. The 15q11.2 BP1-BP2 microdeletion syndrome: a review. Int J Mol Sci. 2015;16(2): Sharp AJ, et al. A recurrent 15q13.3 microdeletion syndrome associated with mental retardation and seizures. Nat Genet. 2008;40(3): van Bon BW, et al. Further delineation of the 15q13 microdeletion and duplication syndromes: a clinical spectrum varying from non-pathogenic to a severe outcome. J Med Genet. 2009;46(8):

10 Peixoto et al. Molecular Cytogenetics (2017) 10:2 Page 10 of Miller DT, et al. Microdeletion/duplication at 15q13.2q13.3 among individuals with features of autism and other neuropsychiatric disorders. J Med Genet. 2009;46(4): Helbig I, et al. 15q13.3 microdeletions increase risk of idiopathic generalized epilepsy. Nat Genet. 2009;41(2): Kolevzon A, et al. Analysis of a purported SHANK3 mutation in a boy with autism: clinical impact of rare variant research in neurodevelopmental disabilities. Brain Res. 2011;1380: Lopez-Rangel E, Lewis ME. Loud and clear evidence for gene silencing by epigenetic mechanisms in autism spectrum and related neurodevelopmental disorders. Clin Genet. 2006;69(1): Krakowiak P, et al. Maternal metabolic conditions and risk for autism and other neurodevelopmental disorders. Pediatrics. 2012;129(5):e Seidman JG, Seidman C. Transcription factor haploinsufficiency: when half a loaf is not enough. J Clin Invest. 2002;109(4): Courtens W, Schramme I, Laridon A. Microduplication 22q11.2: a benign polymorphism or a syndrome with a very large clinical variability and reduced penetrance? Report of two families. Am J Med Genet A. 2008; 146A(6): Submit your next manuscript to BioMed Central and we will help you at every step: We accept pre-submission inquiries Our selector tool helps you to find the most relevant journal We provide round the clock customer support Convenient online submission Thorough peer review Inclusion in PubMed and all major indexing services Maximum visibility for your research Submit your manuscript at