Identification of a RAI1-associated disease network through integration of exome sequencing, transcriptomics, and 3D genomics
- Maria Nicla Loviglio†1,
- Christine R. Beck†2,
- Janson J. White2,
- Marion Leleu3, 4,
- Tamar Harel2,
- Nicolas Guex4,
- Anne Niknejad4,
- Weimin Bi2,
- Edward S. Chen2,
- Isaac Crespo4,
- Jiong Yan2, 8,
- Wu-Lin Charng2,
- Shen Gu2,
- Ping Fang2, 9,
- Zeynep Coban-Akdemir2,
- Chad A. Shaw2,
- Shalini N. Jhangiani5,
- Donna M. Muzny5,
- Richard A. Gibbs2, 5,
- Jacques Rougemont3, 4,
- Ioannis Xenarios1, 4,
- James R. Lupski2, 5, 6, 7 and
- Alexandre Reymond1Email author
© The Author(s). 2016
Received: 29 March 2016
Accepted: 16 September 2016
Published: 1 November 2016
Smith-Magenis syndrome (SMS) is a developmental disability/multiple congenital anomaly disorder resulting from haploinsufficiency of RAI1. It is characterized by distinctive facial features, brachydactyly, sleep disturbances, and stereotypic behaviors.
We investigated a cohort of 15 individuals with a clinical suspicion of SMS who showed neither deletion in the SMS critical region nor damaging variants in RAI1 using whole exome sequencing. A combination of network analysis (co-expression and biomedical text mining), transcriptomics, and circularized chromatin conformation capture (4C-seq) was applied to verify whether modified genes are part of the same disease network as known SMS-causing genes.
Potentially deleterious variants were identified in nine of these individuals using whole-exome sequencing. Eight of these changes affect KMT2D, ZEB2, MAP2K2, GLDC, CASK, MECP2, KDM5C, and POGZ, known to be associated with Kabuki syndrome 1, Mowat-Wilson syndrome, cardiofaciocutaneous syndrome, glycine encephalopathy, mental retardation and microcephaly with pontine and cerebellar hypoplasia, X-linked mental retardation 13, X-linked mental retardation Claes-Jensen type, and White-Sutton syndrome, respectively. The ninth individual carries a de novo variant in JAKMIP1, a regulator of neuronal translation that was recently found deleted in a patient with autism spectrum disorder. Analyses of co-expression and biomedical text mining suggest that these pathologies and SMS are part of the same disease network. Further support for this hypothesis was obtained from transcriptome profiling that showed that the expression levels of both Zeb2 and Map2k2 are perturbed in Rai1 –/– mice. As an orthogonal approach to potentially contributory disease gene variants, we used chromatin conformation capture to reveal chromatin contacts between RAI1 and the loci flanking ZEB2 and GLDC, as well as between RAI1 and human orthologs of the genes that show perturbed expression in our Rai1 –/– mouse model.
These holistic studies of RAI1 and its interactions allow insights into SMS and other disorders associated with intellectual disability and behavioral abnormalities. Our findings support a pan-genomic approach to the molecular diagnosis of a distinctive disorder.
Smith-Magenis syndrome (SMS; MIM #182290) is a rare genomic disorder with a prevalence of 1 in 15,000. It is associated with specific craniofacial dysmorphology, developmental delay (DD), moderate to profound intellectual disability (ID), and self-injurious and stereotypic behaviors [1, 2]. SMS individuals show sleep disturbance with frequent daytime napping and night-time awakenings. They display restricted interest, obsessive thinking, and social responsiveness scale scores consistent with autism spectrum disorder (ASD) . They repetitively mouth objects, rock, spin, or twirl their body, and grind their teeth . This distinctive profile is complemented by specific lick and flip and self-hug behaviors, as well as attachment to people [5–7]. Challenging behaviors such as self-injuries, physical aggression, and destructive behavior are significantly more prevalent in SMS than in ID with mixed etiologies . Self-injuries are present in 70–97 % of individuals and include polyembolokoilamania (insertion of foreign objects into bodily orifices) and onychotillomania (pulling out finger and toe nails). Unusual behaviors can comprise poking others’ eyes, forceful hugging, and punching fists through walls and windows.
Whereas SMS is classically associated with a deletion within cytogenetic G-band 17p11.2 that includes the RAI1 gene (about 90 % of individuals) or a nucleotide variant in that gene (about 5 %) [1, 9–12], some reports suggested genetic heterogeneity as SMS-like individuals were found to recurrently harbor deletions of the 2q37.3 or 2q23.1 cytobands encompassing HDAC4 and MBD5, respectively [13–15]. Similarly, PITX3 was proposed to be responsible for the SMS-like neurobehavioral abnormalities observed in an individual .
Here we use recent advances in genome sequencing technologies to further assess the genetic heterogeneity of SMS and the possible clinical overlap of this syndrome with other intellectual disability and cognitive dysfunction disorders, as some of the seemingly characteristic phenotypic features are non-discriminating among ID syndromes. We also evaluate the pertinence of network interactions and provide experimental data in support of potential molecular diagnoses.
Each of the 149 patients was clinically assessed by their respective physicians. Patients were diagnosed as potentially affected by SMS through clinical assessment. Briefly, all individuals presented intellectual disability and/or developmental delay, and the majority (>75 %) also had sleep disturbances, stereotypies, or other endophenotypes common to SMS (e.g. distinctive facial features, tantrums, self-injurious behaviors, onychotillomania). The clinical presentation of SMS is heterogeneous; therefore, the indication of SMS by a clinician can be either premature in the case of a young infant or possibly a misdiagnosis in an individual with behavioral issues and ID.
Detailed SMS patients’ phenotypes
The detailed phenotype descriptions of 13 of the 15 patients without RAI1 genetic alteration are described in Additional file 1: Supplementary text and Additional file 2: Table S1. The remaining two individuals had no clinical data available.
Array comparative genomic hybridization
Targeted chromosome 17p array comparative genomic hybridization (aCGH) analyses were carried out on each proband as previously reported . Additional genome-wide aCGH was conducted on each person using Baylor Miraca Genetics Laboratory design version 10.1, an Agilent 180 K oligo array. All array data were analyzed as previously described .
To uncover genetic variants associated with the abnormalities shown by the 15 patients without RAI1 genetic alteration, we performed whole-exome sequencing of DNA extracted from blood of the proband and both their parents whenever possible (eight trios) at the Baylor College of Medicine (BCM) Human Genome Sequencing Center (HGSC) via the Baylor-Hopkins Center for Mendelian Genetics. Exomes were captured and sequenced on an Illumina HiSeq platform using previously described methods . Sequence analysis was performed using the HGSC Mercury analysis pipeline (https://www.hgsc.bcm.edu/software/mercury) . Variants were filtered based on inheritance patterns including autosomal recessive, X-linked, and de novo/autosomal dominant. Variants with MAF < 0.05 in control cohorts (Atherosclerosis Risk in Communities (ARIC, https://www2.cscc.unc.edu/aric/), 1000 Genomes project (http://www.1000genomes.org/), the NHLBI Exome Sequencing Project (http://evs.gs.washington.edu/EVS/), and our internal BCM control database of > 5000 exomes generated as a member of the Centers for Mendelian Genomics)  and predicted to be deleterious by SIFT10 and/or PolyPhen were prioritized . Sanger sequencing confirmed putatively causative variants and their familial segregation.
The sequencing variants identified in this manuscript were deposited in ClinVar (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/clinvar/).
The primary sequence of each candidate protein was loaded in Swiss-PdbViewer aligned onto suitable modeling templates retrieved from SWISS-MODEL and superposed in three-dimensional (3D) space using Swiss-PdbViewer [22, 23]. Each variant was modeled in the context of the overall 3D structure to evaluate its potential impact with respect to protein folding, as well as to position of known disease-associated variants. We also assessed if missense variants perturbing the protein function clustered in 3D around key regions of the protein .
The ZEB2 Zinc finger residues 995–1078 were modeled using the pdb entry 1mey as template . MAP2K2 was modeled using both MAP2K2 (pdb entry 1s9i 3.2A resolution  and MAP2K1 (pdb entry 3eqi, 1.9A resolution) structures . The GLDC residues were aligned on the Synechocystis sp. glycine decarboxylase model PCC 6803 (pdb entry 4LHD) . To model the CASK variants, two partial CASK crystal structures (pdb entries 1kwa, chain A  and 1kgd, chain A (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/pubmed/11729206?dopt=Abstract)) covering residues 487–572 and 739–914, respectively, were superposed on the crystal structure of PALS1/Crb (pdb entry 4wsi ) that present 35 % identity with CASK.
Because literature resources do not use entity name in a consistent way, we first checked each gene identifier by using UniProtKB (http://www.uniprot.org) or HUGO Gene Nomenclature Committee (HGNC) database (http://www.genenames.org) in order to retrieve the recommended/approved name, short name(s), alternative and synonymous name(s) if any for each targeted gene, as well as the name(s) of the encoded protein. These were used as singleton and/or pairwise strings to extract information from various literature resources: PubMed (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/pubmed), Google Scholar, iHOP (http://www.ihop-net.org/UniPub/iHOP/), and EVEX (http://evexdb.org/), to cite here the original source of reference for this project. The obtained results were curated and the reported relationships were visualized using Cytoscape (3.2.1; http://www.cytoscape.org/). The connectivity was assessed using the Knet-function, which is based on the adaptation of spatial statistics concepts to network analysis proposed in . The statistical significance of the obtained Knet-function value was calculated with respect to a population of permuted networks (n = 106) derived from the original prior knowledge network. It is worth noting here that the connectivity is not only based on direct but also on indirect connections through shortest paths.
Identification of RAI1 interacting proteins
We identified ZBTB17/MIZ1 and BRD2 as likely interactors for RAI1 with a yeast two-hybrid assay. The yeast two-hybrid assays were performed in collaboration with the company Proteinlinks. Briefly, two fragments of the carboxyl-terminus of mouse Rai1 (a.a 1246–1841 and a.a. 1246–1890) were cloned into pCWX200 as baits. Around 10 million independent complementary DNA (cDNA) library clones (10× library coverage) were screened for protein–protein interactions with both baits. We cultured the Y304 yeast strain on galactose selective medium without leucine, histidine, trytophan, and uracil. Positive clones were replicated onto the four selective plates and examined with URA3 (or LEU2) and LacZ reporters. From this analysis, we identified ZBTB17/MIZ1, BRD2, and SOGA3 as reasonable candidates (at least two clones, supported by both baits) for RAI1 interaction candidates. These interactions were further assessed using co-immunoprecipitation (co-IP) analysis in HEK293 cells. Full-length Rai1 was cloned in pCMV-3xFLAG vector while the three candidates were cloned into pCMV-HA vectors to confirm the yeast two-hybrid results. Lysate from the co-transfected HEK293 cells (RAI1 and one of the candidates) was purified with EZview FLAG-M2 beads (Sigma) and analyzed with rat anti-HA (Abcam) on western blot. The interaction between RAI1 and ZBTB17/MIZ1 was confirmed by co-IP, however BRD2 did not express well enough on western blot, and SOGA3 was too sticky to conduct co-IP with, as it bound to the beads in the absence of FLAG-RAI1 (Additional file 3: Figure S1).
Embryo collection and RNA extraction
Mice were housed in standard specific pathogen-free conditions. All animal studies were conducted under protocols approved by the Baylor Institutional Animal Care and Use Committee and followed NIH guidelines. Timed matings between Rai1 heterozygous females and males in F2 generation in the C57BL/6 Tyr c-Brd and 129SvEv mixed genetic background were implemented to generate Rai1 –/– embryos. To harvest embryos, pregnant females were sacrificed by cervical dislocation and the embryos were dissected from the uterus in ice-cold phosphate buffered saline (PBS) solution. Similar sized embryos at 10.5 days post conception (dpc) were collected in 1.5 mL Eppendorf tubes, frozen immediately in liquid nitrogen, and stored in –80 °C. Portions of the yolk sac were saved for genotyping as described previously . For RNA extraction, the whole embryos were homogenized in Trizol and RNA was extracted according to the manufacturer’s instructions (Invitrogen) followed by purification on columns using an RNeasy mini kit (Qiagen Sciences, Germantown, MD, USA). The RNA integrity, concentration, and overall quality were tested with an Agilent Bioanalyzer 2100 and a NanoDrop ND-1000 spectrophotometer.
Microarray processing and analysis
A total of 5–10 μg of total RNA from each individual embryo of three Rai1 –/– at 10.5 dpc and three wild-type controls were used to produce complementary RNA (cRNA) target microarray transcriptome analyses. Embryos at 10.5 dpc were chosen because Rai1 functions during this stage as indicated by its strong expression and embryonic lethality of Rai1 –/– embryos from 7.5 to 18.5 dpc . In addition, the size of the Rai1 –/– embryos at 10.5 dpc is comparable to that of their wild-type littermates whereas the few surviving Rai1 –/– mice at birth are significantly smaller than the wild-type . The integrity and quality of the extracted RNAs were assessed on a 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). The target was generated using a reverse transcription reaction to produce cDNA (SuperScript Choice System, Gibco), which was subsequently subjected to in vitro transcription with biotinylated cytidine-5′-triphosphate and uridine-5′-triphosphate using ENZo BioArray High Yield RNA Transcript Labeling kit to produce biotinylated cRNA. The target was then fragmented and hybridized to Affymetrix Mouse Genome 430 2.0 Array GeneChips (Affymetrix, Santa Clara, CA, USA) in duplicates using an Affymetrix GeneChip Fluidics Station 400. The arrays were stained with phycoerythrin-coupled avidin and scanned using a GeneArray Scanner 3000. The resultant output was analyzed using Affymetrix Microarray Suite software and examined for excessive background or evidence for RNA degradation. The chips were assessed by scaling factor, average background, percent of probe sets that are present, number of probes present, and the 3′-end to 5′-end probe intensity ratio for housekeeping probe sets (β-actin and GAPDH), as well as the number of probes present for the “spiked in” probe sets (BioB, BioC, BioD, and Crex). All the chips were of good quality, which is further supported by the observations that they have similar RNA degradation patterns and the chips were well replicated within the same genotype group as shown by scatter plot analyses. The criteria for genes differentially expressed are that the log ratio of the normalized expression values in the Rai1 deficient embryos versus the controls is > 0.5 and the P value < 0.05, which empirically gives a very low false detection rate (FDR). The probe sets with very low expression values were filtered out. We analyzed the chromosomal position of all the regulated genes using the chromosomal coordinates within recent genome assemblies of the mouse. The array data were analyzed using the GC-RMA program to estimate the expression measures from the probe level data . The program corrects the background, normalizes the raw perfect match data using the quantile normalization method, and summarizes the probe values to probe set values (expression values, one per probe set per chip), in log2 scale. The fold change for each probe is the log ratio of average expression value in the mutant samples divided by that in the wild type controls. The fold change is considered to be significant if P ≤ 0.05.
Reverse transcription polymerase chain reaction (RT-PCR) validation
For RT-PCR validation of relevant expression targets, 1 μg of total RNA (intact by gel and measured by NanoDrop) was used for RT reactions using the Quanta qScript cDNA synthesis kit. Three separate RT reactions were performed using RNA from both a Rai1 –/– embryo and a wild-type control littermate. The RT reactions and a non-RT reaction using wild-type RNA as well as a water-only control were then run on a gel and all reactions containing both RNA and RT had similar patterns and intensity. A total of 1 μL of each RT reaction was used for subsequent PCR reactions. Primers for PCR were designed to transcript regions of Zeb2, Map2k2, and Rai1 using the UCSC browser version of GRCm38/mm10. The primers (from 5′ to 3′) are as follows:
Zeb2-7 F: CTTCAAGTACAAGCACCACCTGAA
Map2k2-2 F: TGAGAGGATCTCAGAGCTGGGT
Rai1-4 F: ATGTATCCACACCTACCACTACCCAT
4C-seq and 3C-PCR validation assays
Circularized chromosome conformation capture (4C) libraries were prepared from lymphoblastoid cell lines (LCLs) of two age-matched female control individuals. Briefly, LCLs were grown at 37 °C. 5 × 107 exponentially growing cells were harvested and crosslinked with 1 % formaldehyde, lysed, and cut with DpnII, a 4-cutter restriction enzyme that allows higher resolution [34, 35]. After ligation and reversal of the crosslinks, the DNA was purified to obtain the 3C library. This 3C library was further digested with NlaIII and circularized to obtain a 4C library. The inverse PCR primers to amplify 4C-seq (4C combined with multiplexed high-throughput sequencing) templates were designed to contain Illumina adaptor tails, sample barcodes, and viewpoint-specific sequences. The selected viewpoint maps within the 5′ portion of the first intron of the RAI1 gene (700 bp from the donor site of exon 1), a region enriched in DNaseI hypersensitive and transcription factor binding sites . It corresponds to the closest suitable DpnII fragment relative to the transcriptional start sites of the targeted gene. The sequence of the 4C-seq primers is reported in Additional file 2: Table S2. We amplified at least 1.6 μg of 4C template (using about 100 ng of 4C template per inverse PCR reaction, for a total number of 16 PCRs). We multiplexed the two 4C-seq templates in equimolar ratios and analyzed them on a 100-bp single-end Illumina HiSeq flow cell. The numbers of raw, excluded, and mapped reads for each LCL sample are detailed in Additional file 2: Table S3.
To validate selected physical interactions and loop formations between non-neighboring chromatin fragments, 5 × 107 exponentially growing cells were used in conjunction with our 3C protocol as described . We tested primers positioned on the chromosome 17p11.2 sense strand 5′ to 3 for the cis-interactions and primers designed at 9p24.1 compared to control 16p11.2 region for the trans-interactions (Additional file 3: Figure S2 with primers tables). The presence of physical interactions was determined by PCR amplimer production. Control PCRs included no input (“water”) as well as DNA from chromatin digested with DpnII but without the subsequent religation step (“- Ligase”) (Additional file 3: Figure S2).
4C-seq data analysis
4C-seq data were analyzed as previously described [34, 35, 37] through the 4C-seq pipeline available at http://htsstation.epfl.ch/)  and visualized with gFeatBrowser (http://www.gfeatbrowser.com). Briefly, the multiplexed samples were separated, undigested, and self-ligated reads removed. Remaining reads were aligned and translated to a virtual library of DpnII fragments. Read counts were then normalized to the total number of reads and replicates combined by averaging the resulting signal densities (Additional file 3: Figures S3 and S4). The local correlation between the profiles of the two samples per viewpoint was calculated (Spearman correlation: 0.83). The combined profiles were then smoothed with a window size of 29 fragments. The region directly surrounding the viewpoint is usually highly enriched and can show considerable experimental variation, thereby influencing overall fragment count. To minimize these effects, the viewpoint itself and the directly neighboring “undigested” fragment were excluded during the procedure. In addition to this filtering, we modeled the data to apply a profile correction similar to the one described in  using a fit with a slope -1 in a log-log scale . Significantly interacting regions were detected by applying a domainogram analysis as described . We selected BRICKS (Blocks of Regulators In Chromosomal Kontext) with a p value threshold < 0.01 for both “cis” and “trans” interactions, and annotated the BRICKs overlapping genes as well as the closest upstream and the closest downstream genes, in a window of +/– 500 kb. The 4C libraries used to perform the circular PCR with RAI1 viewpoint’s primers had been previously tested in , with seven additional viewpoints’ primer pairs. The BRICKs genes GTDC1 and KDM4C (and the flanking genes ZEB2 and GLDC) were not called as significantly interacting regions for any of these viewpoints (see ; Additional file 2: Tables S6–S12). The raw sequencing files are available at GEO under accession number GSE83420.
Gene annotation was obtained through BioScript (http://gdv.epfl.ch/bs). Protein interaction networks for BRICKs genes were determined using STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) v9.1 (http://string-db.org/) . We exploited GO with Enrichr (http://amp.pharm.mssm.edu/Enrichr/) to assess if the chromatin-contacted genes were enriched in specific pathways and genes associated with Mendelian diseases and GIANT (http://giant.princeton.edu/) and Genemania (http://www.genemania.org/) to test tissue-specific functional interactions and produce association networks, respectively [43–46]. The significance of the connectivity of the GIANT co-expression networks was assessed as described for the literature-mining network (see above). We used Enrichr Chromosome Location tool and BRICKS count in different window sizes (5 Mb, 1 Mb, and 500 kb) to determine whether any cytogenetic band other than 17p11.2 was enriched for BRICKS. Other than 17p11.2, we identified significant enrichments at cytobands 17p12, 17p13, and 2q22, where the gene ZEB2 is located.
Hi-C matrices from Rao et al.  were prepared by first applying a KR normalization to the 5 kb and 100 kb resolution observed matrices and then by dividing each normalized score by the expected one extracted from the KR expected file (as described previously in section II.c of the Extended Experimental Procedures of reference ). KR expected values less than 1 were set to 1 to avoid long-distance interaction biases. HiC matrices from Dixon et al.  were generated from the normalized datasets at a 40 kb resolution and transformed to a 400 kb resolution by summing the contacts observed in 10 × 10 sub-matrices. Expected vectors represent the mean number of contacts observed at a given distance and were used to calculate the observed/expected matrices.
Clinical and molecular findings
Through physicians from a large network of medical genetics centers, we enrolled a cohort of 149 individuals presenting with a constellation of SMS features. High-density 17p11.2 aCGH and Sanger sequencing of RAI1 showed that 134 out of 149 individuals presented a genetic or genomic alteration of the RAI1 gene [9, 11, 17, 49–52], 96/134 (72 %) individuals carried the classic recurrent 3.7 Mb SMS deletion, ten (7.5 %) contained an uncommon recurrent 1 (UR1) or UR2 rearrangement, 24 (18 %) a non-recurrent RAI1 deletion, and four (3 %) had a de novo variant in RAI1 [9, 11, 49, 52, 53] (Additional file 2 Table S1). Whereas these proportions are similar to published results [12, 54], it is likely that some clinicians did not refer individuals with SMS features who were negative for SMS molecular diagnosis (via aCGH or fluorescence in situ hybridization, FISH) or who were positive for another potentially causative CNV, for example 1p36 deletion syndrome [15, 55, 56] that shares multiple similarities with SMS. Indeed, many individuals were molecularly diagnosed prior to sample submission. Consistent with this hypothesis, a separate study identified mutations affecting RAI1 in only 30 % of participants with a suspected diagnosis of SMS .
The remaining 15 individuals (10 %) showed no discernable perturbation of the RAI1 gene. The 13 with available clinical data presented the following classical SMS features: ID (12/12), DD (13/13), sleep disturbances (8/10), and/or self-injurious behavior (10/11), in particular onychotillomania (6/7) (Additional file 1: Supplementary Text, Additional file 1: Table S1). To identify the underlying cause of the phenotypes of these 15 individuals, the probands and their parents when available (eight cases) were subjected to high-resolution genome-wide aCGH and whole-exome sequencing. We identified potentially causative variants in ten individuals (Additional file 2: Table S1). These were grouped into five categories: (1) a 47, XYY karyotype (subject BAB2492); (2) de novo variants in ZEB2 (BAB2386), CASK (BAB2540), KMT2D (BAB2319), and JAKMIP1 (BAB2451); (3) compound heterozygote variants in GLDC (BAB4947); (4) a MECP2 variant in a woman with random X-inactivation (BAB2552) inherited from the individual’s mother, who presented with a skewed X-inactivation pattern (away from this allele) in her blood (Additional file 3: Figure S5); and (5) variants in POGZ (BAB2330, variant not maternally inherited), MAP2K2 (BAB2474), and the X-linked KDM5C (BAB2293), the origins of which could not be assessed. We confirmed the segregation of sequence variants in available family members by Sanger sequencing.
Variant analysis and modeling
The variants identified in KMT2D (p.E3418X) and MECP2 (p.P389fsX) are predicted to be loss-of-function alleles, which are likely pathogenic alleles as KMT2D and MECP2 are “extremely intolerant” and “intolerant” to loss-of-function variation according to the Exome Aggregation Consortium database version 0.3 (http://exac.broadinstitute.org) [10.2015] (pLI = 1.0 and 0.7, respectively) and as analogous loss-of-function variants in KMT2D and MECP2 were identified in KABUK1  and MRXS13  individuals, respectively (Additional file 2: Table S5 and S6). Additionally, the de novo variant in the candidate gene JAKMIP1 (p.D586H) occurs in a highly conserved residue and is predicted to be deleterious to the protein structure. JAKMIP1 is “extremely intolerant” to loss-of-function variation according to ExAC (pLI = 0.99). When possible, we used X-ray structures and/or cryo-EM modeling to obtain a 3D representation of the remaining encoded proteins and compared the variants we identified with those previously reported in MOWS, MICPCH, GCE, CFC4, and MRXSC individuals (Additional file 2: Table S7–S11). By and large, these models suggest that the variants identified in the current study are detrimental to the encoded proteins: (1) the ZEB2 p.H1049P variant substitutes a residue that participates in the coordination of the Zn++ atom of one of the Zinc fingers, similar to the variant p.H1045R identified in a MOWS individual (Additional file 3: Figure S6A; Additional file 2: Table S11); (2) the MAP2K2 p.D69del variant removes one of the two aspartic acid residues involved in the binding of a Ca++ ion in the conserved GELKDD loop (Additional file 3: Figure S6B); (3) the GLDC p.L726Q and p.P647L variants likely affect the packing of the encoded protein in the neighborhood of the catalytic lysine K754 residue similar to the 61 missense variants identified in GCE individuals (Additional file 3: Figure S6C, Additional file 2: Table S8); and (4) the CASK p.R489W variant places a bulky tryptophan sidechain that cannot be accommodated in the structure without changing the molecular surface (Additional file 3: Figure S6D). The possible impact of the KDM5C p.K1023R variant on this conserved position (Additional file 3: Figure S7) could not be evaluated as no template is available for this region.
The identified rare variants affect Rai1-associated genes
To gain further insight about the genes regulated by Rai1 during mouse embryonic development, we performed microarray analysis on total RNA prepared from three 10.5 dpc Rai1 –/– embryos and from three of their wild-type littermates. The two Rai1 transcripts present on the array are significantly downregulated in Rai1 –/– embryos compared to wild-type littermates (e.g. the AK013909 transcript with a fold change of 6.2 shows the largest downregulation among the 45,037 assessed probe sets). In fact, the expression values for both transcripts are within background levels in the Rai1 –/– embryos, indicating that both transcripts are not expressed in the Rai1 –/– mutants and further corroborating the contention that the engineered Rai1 mutant allele is a complete null allele . In total, 142 and 157 probe sets showed an over twofold increase or decrease, respectively (Additional file 2: Table S12; see “Methods”) in the mutant mice when compared to wild-type littermates. Consistent with the hypothesis that genes potentially causative of the SMS-like phenotypes are functionally associated to or transcriptionally regulated by RAI1, the expression levels of both Zeb2 (ENSMUSG00000026872) and Map2k2 (ENSMUSG00000035027) were perturbed in Rai1 –/– mice (Additional file 2 Table S12). These expression arrays results were subsequently confirmed by RT-PCR (Additional file 3: Figure S9). We then assessed the chromosomal position of the dysregulated genes. The enrichment score using a Pearson Chi-square goodness of fit statistic indicated that they showed a biased chromosome distribution with 22 % of the genes downregulated and 26 % of the genes upregulated in the Rai1 mutants mapping to mouse chromosome 11 (MMU11) where the Rai1 gene resides. Less than 5 % of the differentially expressed genes are located on any chromosome other than MMU11. This enrichment on MMU11 for downregulated and upregulated genes in Rai1 –/– embryos is reminiscent of our previous finding that the engineered MMU11 deletion and reciprocal duplication that mimic SMS and Potocki-Lupski syndrome rearrangements were associated with a MMU11-wide transcriptome perturbation in the five assessed adult male tissues .
The genes mapping within the RAI1-chromatin contacted genomic loci (BRICKs genes) are enriched for genes that encode proteins that interact together (82 observed interactions versus 35 expected; P = 6.41e–12). BRICKs genes are also enriched for the GO term “detection of light stimulus involved in sensory perception” in Enrichr (P = 5.45e–3) (see “Methods,” Additional file 2: Table S14). Similarly, Enrichr showed that chromosome contacts were enriched in interchromosomal and intrachromosomal cytobands (17p11, adjusted P < 1e–09; 17p12, adjusted P = 9.7e–09; 17p13, adjusted P = 1.8e–03; and 2q22 adjusted P = 4.95e–02). ZEB2, one of the eight genes found mutated in the SMS individuals, maps to the latter 2q22.3 region and is flanked by BRICKs. To further assess possible functional relationships between RAI1 and chromatin-contacted genes, we retrieved the list of 322 genes flanking the BRICKs (BRICKs flanking genes, i.e. the closest genes to be found upstream and downstream of a BRICK within a 500 kb window). The 4C assays in particular identified interchromosomal contacts with restriction fragments mapping 200 kb away from the ZEB2 and GLDC gene loci. We then compared the lists of BRICKs genes and BRICKs flanking genes with the list of genes whose expression levels were perturbed in Rai1 –/– mouse embryos. Although our analysis is restricted by a small sample size, we found a consistent trend of over-representation (Fisher’s enrichment test, P = 0.22, OR = 1.5, and P = 0.2, OR = 1.4) with 10 and 18 chromatin-contacted BRICKs genes and BRICKs flanking genes, respectively, differentially expressed in the mouse knockdown model. Interestingly, 6/10 of these BRICKs genes mapping at cytobands 17p13, 17p11 (2 genes), 17q21 and 17q23 (2 genes) have mouse orthologs that map on mouse chromosome MMU11, thus possibly explaining the enrichment of MMU11-mapping genes within genes differentially expressed in Rai1 –/– mouse embryos (Additional file 3: Figure S12).
Within a cohort of 149 individuals presenting clinical features of SMS we identified 90 % (134/149) of individuals with either a heterozygous deletion of RAI1 or a predicted deleterious variant of the RAI1 gene. We used recent advances in genome sequencing technologies to possibly identify genetic alteration(s) associated with SMS in the remaining individuals. These strategies were successfully applied to discover loci associated with ID . They revealed a large genic overlap between ID and ASD, schizophrenia, and epileptic encephalopathy , suggesting that some developmental disorders have highly variable clinical presentations. They similarly uncovered limitations to the phenotype-driven strategy and conventional clinical paradigm of identifying individuals with very similar presentations as they revealed an unsuspected phenotypic variance of known disorders [73, 74].
It is important to the medical community to identify phenotypic overlap between diseases, which suggests common causes and alterations of the same pathways, as this knowledge could be exploited therapeutically. In this report, we identify previously unappreciated relationships between SMS and its major driver RAI1 and other diseases that include MOWS, MICPCH, KABUK1, GCE, MRXS13, CFC4, and MRXSC. Literature mining, co-expression data, transcriptome profiling of Rai1 –/– animal models, and chromosomal contacts support the existence of a comprehensive “biological module”  or “disease network”  underlying these diseases.
Although none of the 15 individuals described in this study have traditional molecular diagnoses involving RAI1 haploinsufficiency and thus should formally be considered misdiagnoses, many have phenotypes with considerable overlap with SMS (Fig. 3, Additional file 2: Table S1, Additional file 1: Supplementary text). BAB4947 presented facial dysmorphisms, SMS-like behavioral disturbances that include sleep problems, polyembolokoilamania, onychotillomania, brachycephaly, and brachydactyly, as well as known GLDC-variants associated features such as seizures. His clinical diagnosis could possibly be confounded by the likely presence of two molecular diagnoses: compound heterozygous variants in GLDC and an inherited frameshift variant in TCOF1, a gene associated with Treacher Collins syndrome-1 (OMIM #154500) and possibly responsible for the down-slanting eyes, everted lateral eyelids, and malar hypoplasia. The clinical scenario is similar with cases BAB2474 and BAB2540, who did not show CFC4- (e.g. ectodermal anomalies, craniofacial features) and MICPCH-distinctive features (e.g. microcephaly and pontocerebellar hypoplasia). Likewise individual BA2492 has a constellation of symptoms (sleep disturbance, DD, cognitive impairment, brachydactyly) compatible with only the most severe 47,XYY sex chromosome aneuploidy cases . Consistent with the hypothesis of expanded phenotypes, the phenotypic variability of White-Sutton syndrome associated with variants in POGZ keeps extending with clinical features including ASD, DD, ID, schizophrenia, and microcephaly [57, 71, 77–83]. We can also not formally rule out that we have not yet determined the true genetic cause(s) of the phenotypic spectrum of these individuals or they occur in presence of more complex, blended phenotypes as exemplified by BAB4947 above. BAB2451 harbors a “probably pathogenic” variant (SIFT HumDiv score = 1; HumVar score = 0.982) in the gene JAKMIP1. Recent findings have linked the loss of JAKMIP1 to neuronal translation dysregulation during synaptic development; mice knocked out for the JAKMIP1 paralog display social deficits, stereotyped activity, altered vocal communication, increased impulsivity, and other autistic-like behaviors .
The presented results support the notion that at least some of the identified variants in candidate SMS contributory genes CASK, GLDC, KDM5C, KMT2D, MAP2K2, MECP2, POGZ, and ZEB2 are causative of the observed phenotypes and thus that modification of the function of these genes is associated with a greater phenotypic variability than previously expected (Fig. 3). Conversely, one and two carriers of damaging RAI1 variants were identified within a total of 6381 ASD [79, 85] and 2426 ID [71, 78, 86–89] individuals, respectively. Whereas the phenotype of one of the ID individuals was retrospectively found to be consistent with SMS , we lack detailed phenotypic information regarding the other two cases. If we assume that these two individuals do not present with typical SMS features that would have excluded them from these cohorts, it suggests that the phenotype of carriers of RAI1 deleterious variants is similarly more variable than anticipated.
Structural variations, especially large rearrangements involving several genes, shape tissue transcriptomes and impact the expression of genes mapping to their flanks [64, 90]. We show that the homozygous deletion of Rai1 in mouse embryos  influences the expression of several genes and in particular MMU11 genes. Furthermore, the RAI1 viewpoint contacts the orthologous genes at the chromatin level. As some of these genes contribute to phenotypes associated with RAI1 variation (e.g. KRT17 with “hoarse voice” (HP:0001609), B9D1 with “low-set, posteriorly rotated ears” (HP:0000368), “hypertelorism” (HP:0000316), and “microcornea” (HP:0000482)), they could be involved in RAI1 pathways. The relevance of using 3C-based approaches as unbiased tools to discover clinically related genes is reinforced by their successful application in assessing connected regions involved in similar phenotypes  and genes interacting with risk loci identified in genome-wide association studies (GWAS) [91, 92]. The contacted regions encompass candidate genes involved in “detection of light stimulus” and related gene ontology terms. These processes all refer to photodetection, which controls circadian rhythm and melatonin production from the pineal gland. RAI1 is an important player in this mechanism, by controlling the transcriptional levels of CLOCK, a key component of the mammalian circadian oscillator that transcriptionally regulates many critical circadian genes . Another gene mapping within the SMS critical region on chromosome 17p11.2 and linked to these processes is the subunit 3 of the COP9 signal transduction complex (COPS3), essential for the light control of gene expression . It is thus possible that the disruption of the orthologous locus in the Rai1 –/– mice perturbs chromatin loops and affects expression levels of RAI1-contacted/functionally associated genes. We are well aware of the limitations of using LCLs in this type of study, and particularly to assess chromatin contacts between genes whose expression specificity resides in other cell lineages. These experiments are nevertheless worth pursuing simply because: (1) the primary human target tissues often remain beyond reach; (2) we cannot exclude a broad to ubiquitous expression pattern for the genes involved in these disease processes; and (3) long range chromatin contacts were shown to be stable across cell lines and tissues regardless of expression status . Similar limitations apply to the use of embryonic stem cell-derived material, while animal tissues have a different set of shortcomings.
Our results strongly support a disease network associated with RAI1 and illustrate the utility of a comprehensive multifaceted diagnostic approach even in the presence of a distinctive disorder.
We thank the patients and families for their contribution to this study. We are grateful to the members of the Lausanne Genomic Technologies Facility for technical help.
MNL was awarded an EMBO fellowship (ASTF 153-2015). CRB was an HHMI fellow of the Damon Runyon Cancer Research Foundation (DRG 2155-13) and is supported by a grant from the National Institutes of General Medical Sciences (K99GM120453). TH and WLC are supported by the NIH T32 GM07526 Medical Genetics Research Fellowship Program and the CPRIT RP140102 training Program, respectively. This work was supported by grants from the Swiss National Science Foundation (31003A_160203) and the Simons Foundation (SFARI274424) to AR; the US National Human Genome Research Institute (NHGRI)/National Heart Lung and Blood Institute (NHLBI) grant no. HG006542 to the Baylor-Hopkins Center for Mendelian Genomics; the Smith-Magenis Syndrome Research Foundation (SMSRF); and the National Institute of Neurological Disorders and Stroke (NINDS) NS058529 to JRL. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and material
The 4C-seq raw sequencing files and identified variants are available at GEO and ClinVar under accession number GSE83420 and SCV000299207-SCV000299218, respectively.
MNL performed the 3C and 4C experiments and analyzed the data with ML and JR. AN, MNL, IC, and IX conducted the network analysis; NG performed the variant analysis and modeling. CRB, JW, TH, ZCA, SNJ, DMM, RAG, and JRL performed the exome sequencing, analysis of the data, and subsequent validation of variants. PF, CRB, and JW conducted and analyzed X-inactivation data. WB, CRB, ESC, SG, and CAS performed and analyzed expression data on mouse embryos. JY, WLC, CRB, and WB performed and analyzed yeast two-hybrid data. MNL, CRB, JRL, and AR wrote the manuscript with contributions from IX, NG, and TH. JRL and AR designed the study and obtained the necessary funding. All authors commented on and approved the manuscript.
JRL has stock ownership in 23andMe, is a paid consultant for Regeneron Pharmaceuticals, has stock options in Lasergen, Inc., is a member of the Scientific Advisory Board of Baylor Miraca Genetics Laboratories, and is a co-inventor on multiple United States and European patents related to molecular diagnostics for inherited neuropathies, eye diseases, and bacterial genomic fingerprinting. Baylor College of Medicine (BCM) and Miraca Holdings Inc. have formed a joint venture with shared ownership and governance of the Baylor Miraca Genetics Laboratories (BMGL), which performs clinical exome sequencing. The Department of Molecular and Human Genetics at Baylor College of Medicine derives revenue from the chromosomal microarray analysis (CMA) and clinical exome sequencing offered in the Baylor Miraca Genetics Laboratory (BMGL; http://www.bmgl.com/BMGL/Default.aspx). The remaining authors declare that they have no competing interests.
Consent for publication
We obtained the authorization to publish participants’ data from their parents or legal guardians.
Ethics approval and consent to participate
The institutional review board of the Baylor College of Medicine approved this study. Participants were enrolled after written informed consent was obtained from parents or legal guardians. This study conforms to the Helsinki Declaration.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Smith ACM, Boyd KE, Elsea SH, Finucane BM, Haas-Givler B, Gropman A, et al. Smith-Magenis Syndrome. In: Pagon RA, Adam MP, Ardinger HH, Wallace SE, Amemiya A, Bean LJH, et al., editors. GeneReviews(R). Seattle: University of Washington; 1993.Google Scholar
- Elsea SH, Girirajan S. Smith-Magenis syndrome. Eur J Hum Genet. 2008;16(4):412–21.PubMedView ArticleGoogle Scholar
- Laje G, Morse R, Richter W, Ball J, Pao M, Smith AC. Autism spectrum features in Smith-Magenis syndrome. Am J Med Genet C Semin Med Genet. 2010;154C(4):456–62.PubMedPubMed CentralView ArticleGoogle Scholar
- Moss J, Oliver C, Arron K, Burbidge C, Berg K. The prevalence and phenomenology of repetitive behavior in genetic syndromes. J Autism Dev Disord. 2009;39(4):572–88.PubMedView ArticleGoogle Scholar
- Dykens EM, Finucane BM, Gayley C. Brief report: cognitive and behavioral profiles in persons with Smith-Magenis syndrome. J Autism Dev Disord. 1997;27(2):203–11.PubMedView ArticleGoogle Scholar
- Dykens EM, Smith AC. Distinctiveness and correlates of maladaptive behaviour in children and adolescents with Smith-Magenis syndrome. J Intellect Disabil Res. 1998;42(Pt 6):481–9.PubMedView ArticleGoogle Scholar
- Finucane BM, Konar D, Haas-Givler B, Kurtz MB, Scott Jr CI. The spasmodic upper-body squeeze: a characteristic behavior in Smith-Magenis syndrome. Dev Med Child Neurol. 1994;36(1):78–83.PubMedView ArticleGoogle Scholar
- Sloneem J, Oliver C, Udwin O, Woodcock KA. Prevalence, phenomenology, aetiology and predictors of challenging behaviour in Smith-Magenis syndrome. J Intellect Disabil Res. 2011;55(2):138–51.PubMedView ArticleGoogle Scholar
- Bi W, Saifi GM, Shaw CJ, Walz K, Fonseca P, Wilson M, et al. Mutations of RAI1, a PHD-containing protein, in nondeletion patients with Smith-Magenis syndrome. Hum Genet. 2004;115(6):515–24.PubMedView ArticleGoogle Scholar
- Girirajan S, Elsas 2nd LJ, Devriendt K, Elsea SH. RAI1 variations in Smith-Magenis syndrome patients without 17p11.2 deletions. J Med Genet. 2005;42(11):820–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Slager RE, Newton TL, Vlangos CN, Finucane B, Elsea SH. Mutations in RAI1 associated with Smith-Magenis syndrome. Nat Genet. 2003;33(4):466–8.PubMedView ArticleGoogle Scholar
- Vieira GH, Rodriguez JD, Carmona-Mora P, Cao L, Gamba BF, Carvalho DR, et al. Detection of classical 17p11.2 deletions, an atypical deletion and RAI1 alterations in patients with features suggestive of Smith-Magenis syndrome. Eur J Hum Genet. 2012;20(2):148–54.PubMedView ArticleGoogle Scholar
- Mullegama SV, Pugliesi L, Burns B, Shah Z, Tahir R, Gu Y, et al. MBD5 haploinsufficiency is associated with sleep disturbance and disrupts circadian pathways common to Smith-Magenis and fragile X syndromes. Eur J Hum Genet. 2015;23(6):781–9.PubMedView ArticleGoogle Scholar
- Williams SR, Aldred MA, Der Kaloustian VM, Halal F, Gowans G, McLeod DR, et al. Haploinsufficiency of HDAC4 causes brachydactyly mental retardation syndrome, with brachydactyly type E, developmental delays, and behavioral problems. Am J Hum Genet. 2010;87(2):219–28.PubMedPubMed CentralView ArticleGoogle Scholar
- Williams SR, Girirajan S, Tegay D, Nowak N, Hatchwell E, Elsea SH. Array comparative genomic hybridisation of 52 subjects with a Smith-Magenis-like phenotype: identification of dosage sensitive loci also associated with schizophrenia, autism, and developmental delay. J Med Genet. 2010;47(4):223–9.PubMedView ArticleGoogle Scholar
- Derwinska K, Mierzewska H, Goszczanska A, Szczepanik E, Xia Z, Kusmierska K, et al. Clinical improvement of the aggressive neurobehavioral phenotype in a patient with a deletion of PITX3 and the absence of L-DOPA in the cerebrospinal fluid. Am J Med Genet B Neuropsychiatr Genet. 2012;159B(2):236–42.PubMedView ArticleGoogle Scholar
- Liu P, Lacaria M, Zhang F, Withers M, Hastings PJ, Lupski JR. Frequency of nonallelic homologous recombination is correlated with length of homology: evidence that ectopic synapsis precedes ectopic crossing-over. Am J Hum Genet. 2011;89(4):580–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Wiszniewska J, Bi W, Shaw C, Stankiewicz P, Kang SH, Pursley AN, et al. Combined array CGH plus SNP genome analyses in a single assay for optimized clinical testing. Eur J Hum Genet. 2014;22(1):79–87.PubMedView ArticleGoogle Scholar
- Yuan B, Pehlivan D, Karaca E, Patel N, Charng WL, Gambin T, et al. Global transcriptional disturbances underlie Cornelia de Lange syndrome and related phenotypes. J Clin Invest. 2015;125(2):636–51.PubMedPubMed CentralView ArticleGoogle Scholar
- Reid JG, Carroll A, Veeraraghavan N, Dahdouli M, Sundquist A, English A, et al. Launching genomics into the cloud: deployment of Mercury, a next generation sequence analysis pipeline. BMC Bioinformatics. 2014;15:30.PubMedPubMed CentralView ArticleGoogle Scholar
- Chong JX, Buckingham KJ, Jhangiani SN, Boehm C, Sobreira N, Smith JD, et al. The genetic basis of Mendelian phenotypes: discoveries, challenges, and opportunities. Am J Hum Genet. 2015;97(2):199–215.PubMedPubMed CentralView ArticleGoogle Scholar
- Guex N, Peitsch MC. SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis. 1997;18(15):2714–23.PubMedView ArticleGoogle Scholar
- Schwede T, Kopp J, Guex N, Peitsch MC. SWISS-MODEL: An automated protein homology-modeling server. Nucleic Acids Res. 2003;31(13):3381–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Kamburov A, Lawrence MS, Polak P, Leshchiner I, Lage K, Golub TR, et al. Comprehensive assessment of cancer missense mutation clustering in protein structures. Proc Natl Acad Sci U S A. 2015;112(40):E5486–95.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim CA, Berg JM. A 2.2 A resolution crystal structure of a designed zinc finger protein bound to DNA. Nat Struct Biol. 1996;3(11):940–5.PubMedView ArticleGoogle Scholar
- Ohren JF, Chen H, Pavlovsky A, Whitehead C, Zhang E, Kuffa P, et al. Structures of human MAP kinase kinase 1 (MEK1) and MEK2 describe novel noncompetitive kinase inhibition. Nat Struct Mol Biol. 2004;11(12):1192–7.PubMedView ArticleGoogle Scholar
- Fischmann TO, Smith CK, Mayhood TW, Myers JE, Reichert P, Mannarino A, et al. Crystal structures of MEK1 binary and ternary complexes with nucleotides and inhibitors. Biochemistry. 2009;48(12):2661–74.PubMedView ArticleGoogle Scholar
- Hasse D, Andersson E, Carlsson G, Masloboy A, Hagemann M, Bauwe H, et al. Structure of the homodimeric glycine decarboxylase P-protein from Synechocystis sp. PCC 6803 suggests a mechanism for redox regulation. J Biol Chem. 2013;288(49):35333–45.PubMedPubMed CentralView ArticleGoogle Scholar
- Daniels DL, Cohen AR, Anderson JM, Brunger AT. Crystal structure of the hCASK PDZ domain reveals the structural basis of class II PDZ domain target recognition. Nat Struct Biol. 1998;5(4):317–25.PubMedView ArticleGoogle Scholar
- Li Y, Spangenberg O, Paarmann I, Konrad M, Lavie A. Structural basis for nucleotide-dependent regulation of membrane-associated guanylate kinase-like domains. J Biol Chem. 2002;277(6):4159–65.PubMedView ArticleGoogle Scholar
- Cornish AJ, Markowetz F. SANTA: quantifying the functional content of molecular networks. PLoS Comput Biol. 2014;10(9):e1003808.PubMedPubMed CentralView ArticleGoogle Scholar
- Bi W, Ohyama T, Nakamura H, Yan J, Visvanathan J, Justice MJ, et al. Inactivation of Rai1 in mice recapitulates phenotypes observed in chromosome engineered mouse models for Smith-Magenis syndrome. Hum Mol Genet. 2005;14(8):983–95.PubMedView ArticleGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4(2):249–64.PubMedView ArticleGoogle Scholar
- Noordermeer D, Leleu M, Splinter E, Rougemont J, De Laat W, Duboule D. The dynamic architecture of Hox gene clusters. Science. 2011;334(6053):222–5.PubMedView ArticleGoogle Scholar
- Loviglio MN, Leleu M, Mannik K, Passeggeri M, Giannuzzi G, van der Werf I, et al. Chromosomal contacts connect loci associated with autism, BMI and head circumference phenotypes. Mol Psychiatry. 2016. doi:10.1038/mp.2016.84.
- ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(7414):57–74.View ArticleGoogle Scholar
- Gheldof N, Leleu M, Noordermeer D, Rougemont J, Reymond A. Detecting long-range chromatin interactions using the chromosome conformation capture sequencing (4C-seq) method. Methods Mol Biol. 2012;786:211–25.PubMedView ArticleGoogle Scholar
- David FP, Delafontaine J, Carat S, Ross FJ, Lefebvre G, Jarosz Y, et al. HTSstation: a web application and open-access libraries for high-throughput sequencing data analysis. PLoS One. 2014;9(1):e85879.PubMedPubMed CentralView ArticleGoogle Scholar
- Tolhuis B, Blom M, Kerkhoven RM, Pagie L, Teunissen H, Nieuwland M, et al. Interactions among Polycomb domains are guided by chromosome architecture. PLoS Genet. 2011;7(3):e1001343.PubMedPubMed CentralView ArticleGoogle Scholar
- Lieberman-Aiden E, van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, et al. Comprehensive mapping of long-range interactions reveals folding principles of the human genome. Science. 2009;326(5950):289–93.PubMedPubMed CentralView ArticleGoogle Scholar
- de Wit E, Braunschweig U, Greil F, Bussemaker HJ, van Steensel B. Global chromatin domain organization of the Drosophila genome. PLoS Genet. 2008;4(3):e1000045.PubMedPubMed CentralView ArticleGoogle Scholar
- Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41(Database issue):D808–15.PubMedView ArticleGoogle Scholar
- Alexa A, Rahnenfuhrer J. topGO: Enrichment Analysis for Gene Ontology. R package version 2.24.0. 2016. https://bioconductor.org/packages/release/bioc/html/topGO.html.
- da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.View ArticleGoogle Scholar
- Reimand J, Arak T, Vilo J. g:Profiler--a web server for functional interpretation of gene lists (2011 update). Nucleic Acids Res. 2011;39(Web Server issue):W307–15.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.PubMedPubMed CentralView ArticleGoogle Scholar
- Rao SS, Huntley MH, Durand NC, Stamenova EK, Bochkov ID, Robinson JT, et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell. 2014;159(7):1665–80.PubMedView ArticleGoogle Scholar
- Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485(7398):376–80.PubMedPubMed CentralView ArticleGoogle Scholar
- Greenberg F, Guzzetta V, Montes de Oca-Luna R, Magenis RE, Smith AC, Richter SF, et al. Molecular analysis of the Smith-Magenis syndrome: a possible contiguous-gene syndrome associated with del(17)(p11.2). Am J Hum Genet. 1991;49(6):1207–18.PubMedPubMed CentralGoogle Scholar
- Juyal RC, Figuera LE, Hauge X, Elsea SH, Lupski JR, Greenberg F, et al. Molecular analyses of 17p11.2 deletions in 62 Smith-Magenis syndrome patients. Am J Hum Genet. 1996;58(5):998–1007.PubMedPubMed CentralGoogle Scholar
- Moncla A, Piras L, Arbex OF, Muscatelli F, Mattei MG, Mattei JF, et al. Physical mapping of microdeletions of the chromosome 17 short arm associated with Smith-Magenis syndrome. Hum Genet. 1993;90(6):657–60.PubMedView ArticleGoogle Scholar
- Shaw CJ, Withers MA, Lupski JR. Uncommon deletions of the Smith-Magenis syndrome region can be recurrent when alternate low-copy repeats act as homologous recombination substrates. Am J Hum Genet. 2004;75(1):75–81.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen KS, Manian P, Koeuth T, Potocki L, Zhao Q, Chinault AC, et al. Homologous recombination of a flanking repeat gene cluster is a mechanism for a common contiguous gene deletion syndrome. Nat Genet. 1997;17(2):154–63.PubMedView ArticleGoogle Scholar
- Vilboux T, Ciccone C, Blancato JK, Cox GF, Deshpande C, Introne WJ, et al. Molecular analysis of the Retinoic Acid Induced 1 gene (RAI1) in patients with suspected Smith-Magenis syndrome without the 17p11.2 deletion. PLoS One. 2011;6(8):e22861.PubMedPubMed CentralView ArticleGoogle Scholar
- Battaglia A. Commentary: Recognizing syndromes with overlapping features: How difficult is it? Considerations generated by the article on differential diagnosis of Smith-Magenis syndrome by Vieira and colleagues. Am J Med Genet A. 2011;155A(5):986–7.PubMedView ArticleGoogle Scholar
- Vieira GH, Rodriguez JD, Boy R, de Paiva IS, DuPont BR, Moretti-Ferreira D, et al. Differential diagnosis of Smith-Magenis syndrome: 1p36 deletion syndrome. Am J Med Genet A. 2011;155A(5):988–92.PubMedView ArticleGoogle Scholar
- White J, Beck CR, Harel T, Posey JE, Jhangiani SN, Tang S, et al. POGZ truncating alleles cause syndromic intellectual disability. Genome Med. 2016;8(1):3.PubMedPubMed CentralView ArticleGoogle Scholar
- Stessman HA, Willemsen MH, Fenckova M, Penn O, Hoischen A, Xiong B, et al. Disruption of POGZ is associated with intellectual disability and autism spectrum disorders. Am J Hum Genet. 2016;98(3):541–52.PubMedPubMed CentralView ArticleGoogle Scholar
- Edelman EA, Girirajan S, Finucane B, Patel PI, Lupski JR, Smith AC, et al. Gender, genotype, and phenotype differences in Smith-Magenis syndrome: a meta-analysis of 105 cases. Clin Genet. 2007;71(6):540–50.PubMedView ArticleGoogle Scholar
- Micale L, Augello B, Fusco C, Selicorni A, Loviglio MN, Silengo MC, et al. Mutation spectrum of MLL2 in a cohort of Kabuki syndrome patients. Orphanet J Rare Dis. 2011;6:38.PubMedPubMed CentralView ArticleGoogle Scholar
- Gomot M, Gendrot C, Verloes A, Raynaud M, David A, Yntema HG, et al. MECP2 gene mutations in non-syndromic X-linked mental retardation: phenotype-genotype correlation. Am J Med Genet A. 2003;123A(2):129–39.PubMedView ArticleGoogle Scholar
- Menche J, Sharma A, Kitsak M, Ghiassian SD, Vidal M, Loscalzo J, et al. Disease networks. Uncovering disease-disease relationships through the incomplete interactome. Science. 2015;347(6224):1257601.PubMedPubMed CentralView ArticleGoogle Scholar
- Greene CS, Krishnan A, Wong AK, Ricciotti E, Zelaya RA, Himmelstein DS, et al. Understanding multicellular function and disease with human tissue-specific networks. Nat Genet. 2015;47(6):569–76.PubMedPubMed CentralView ArticleGoogle Scholar
- Ricard G, Molina J, Chrast J, Gu W, Gheldof N, Pradervand S, et al. Phenotypic consequences of copy number variation: insights from Smith-Magenis and Potocki-Lupski syndrome mouse models. PLoS Biol. 2010;8(11):e1000543.PubMedPubMed CentralView ArticleGoogle Scholar
- de Laat W, Duboule D. Topology of mammalian developmental enhancers and their regulatory landscapes. Nature. 2013;502(7472):499–506.PubMedView ArticleGoogle Scholar
- Fanucchi S, Shibayama Y, Burd S, Weinberg MS, Mhlanga MM. Chromosomal contact permits transcription between coregulated genes. Cell. 2013;155(3):606–20.PubMedView ArticleGoogle Scholar
- Simonis M, Klous P, Splinter E, Moshkin Y, Willemsen R, de Wit E, et al. Nuclear organization of active and inactive chromatin domains uncovered by chromosome conformation capture-on-chip (4C). Nat Genet. 2006;38(11):1348–54.PubMedView ArticleGoogle Scholar
- Simonis M, Kooren J, de Laat W. An evaluation of 3C-based methods to capture DNA interactions. Nat Methods. 2007;4(11):895–901.PubMedView ArticleGoogle Scholar
- Gheldof N, Witwicki RM, Migliavacca E, Leleu M, Didelot G, Harewood L, et al. Structural variation-associated expression changes are paralleled by chromatin architecture modifications. PLoS One. 2013;8(11):e79973.PubMedPubMed CentralView ArticleGoogle Scholar
- de Wit E, de Laat W. A decade of 3C technologies: insights into nuclear organization. Genes Dev. 2012;26(1):11–24.PubMedPubMed CentralView ArticleGoogle Scholar
- Deciphering Developmental Disorders Study. Large-scale discovery of novel genetic causes of developmental disorders. Nature. 2015;519(7542):223–8.Google Scholar
- Vissers LE, Gilissen C, Veltman JA. Genetic studies in intellectual disability and related disorders. Nat Rev Genet. 2016;17(1):9–18.PubMedView ArticleGoogle Scholar
- D’Angelo D, Lebon S, Chen Q, Martin-Brevet S, Snyder LG, Hippolyte L, et al. Defining the effect of the 16p11.2 duplication on cognition, behavior, and medical comorbidities. JAMA Psychiatry. 2016;73(1):20–30.PubMedView ArticleGoogle Scholar
- Mannik K, Magi R, Mace A, Cole B, Guyatt AL, Shihab HA, et al. Copy number variations and cognitive phenotypes in unselected populations. JAMA. 2015;313(20):2044–54.PubMedPubMed CentralView ArticleGoogle Scholar
- Kochinke K, Zweier C, Nijhof B, Fenckova M, Cizek P, Honti F, et al. Systematic phenomics analysis deconvolutes genes mutated in intellectual disability into biologically coherent modules. Am J Hum Genet. 2016;98(1):149–64.PubMedPubMed CentralView ArticleGoogle Scholar
- Bardsley MZ, Kowal K, Levy C, Gosek A, Ayari N, Tartaglia N, et al. 47, XYY syndrome: clinical phenotype and timing of ascertainment. J Pediatr. 2013;163(4):1085–94.PubMedPubMed CentralView ArticleGoogle Scholar
- Fromer M, Pocklington AJ, Kavanagh DH, Williams HJ, Dwyer S, Gormley P, et al. De novo mutations in schizophrenia implicate synaptic networks. Nature. 2014;506(7487):179–84.PubMedPubMed CentralView ArticleGoogle Scholar
- Gilissen C, Hehir-Kwa JY, Thung DT, van de Vorst M, van Bon BW, Willemsen MH, et al. Genome sequencing identifies major causes of severe intellectual disability. Nature. 2014;511(7509):344–7.PubMedView ArticleGoogle Scholar
- Iossifov I, O’Roak BJ, Sanders SJ, Ronemus M, Krumm N, Levy D, et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature. 2014;515(7526):216–21.PubMedPubMed CentralView ArticleGoogle Scholar
- Iossifov I, Ronemus M, Levy D, Wang Z, Hakker I, Rosenbaum J, et al. De novo gene disruptions in children on the autistic spectrum. Neuron. 2012;74(2):285–99.PubMedPubMed CentralView ArticleGoogle Scholar
- Isidor B, Kury S, Rosenfeld JA, Besnard T, Schmitt S, Joss S, et al. De novo truncating mutations in the kinetochore-microtubules attachment gene CHAMP1 cause syndromic intellectual disability. Hum Mutat. 2016;37(4):354–8.PubMedView ArticleGoogle Scholar
- Neale BM, Kou Y, Liu L, Ma’ayan A, Samocha KE, Sabo A, et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012;485(7397):242–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Tan B, Zou Y, Zhang Y, Zhang R, Ou J, Shen Y, et al. A novel de novo POGZ mutation in a patient with intellectual disability. J Hum Genet. 2016;61(4):357–9.PubMedView ArticleGoogle Scholar
- Berg JM, Lee C, Chen L, Galvan L, Cepeda C, Chen JY, et al. JAKMIP1, a novel regulator of neuronal translation, modulates synaptic function and autistic-like behaviors in mouse. Neuron. 2015;88(6):1173–91.PubMedView ArticleGoogle Scholar
- De Rubeis S, He X, Goldberg AP, Poultney CS, Samocha K, Ercument Cicek A, et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature. 2014;515(7526):209–15.PubMedPubMed CentralView ArticleGoogle Scholar
- de Ligt J, Willemsen MH, van Bon BW, Kleefstra T, Yntema HG, Kroes T, et al. Diagnostic exome sequencing in persons with severe intellectual disability. N Engl J Med. 2012;367(20):1921–9.PubMedView ArticleGoogle Scholar
- Grozeva D, Carss K, Spasic-Boskovic O, Tejada MI, Gecz J, Shaw M, et al. Targeted next-generation sequencing analysis of 1,000 individuals with intellectual disability. Hum Mutat. 2015;36(12):1197–204.PubMedPubMed CentralView ArticleGoogle Scholar
- Rauch A, Wieczorek D, Graf E, Wieland T, Endele S, Schwarzmayr T, et al. Range of genetic mutations associated with severe non-syndromic sporadic intellectual disability: an exome sequencing study. Lancet. 2012;380(9854):1674–82.PubMedView ArticleGoogle Scholar
- Redin C, Gerard B, Lauer J, Herenger Y, Muller J, Quartier A, et al. Efficient strategy for the molecular diagnosis of intellectual disability using targeted high-throughput sequencing. J Med Genet. 2014;51(11):724–36.PubMedPubMed CentralView ArticleGoogle Scholar
- Henrichsen CN, Vinckenbosch N, Zollner S, Chaignat E, Pradervand S, Schutz F, et al. Segmental copy number variation shapes tissue transcriptomes. Nat Genet. 2009;41(4):424–9.PubMedView ArticleGoogle Scholar
- Du M, Yuan T, Schilter KF, Dittmar RL, Mackinnon A, Huang X, et al. Prostate cancer risk locus at 8q24 as a regulatory hub by physical interactions with multiple genomic loci across the genome. Hum Mol Genet. 2015;24(1):154–66.PubMedView ArticleGoogle Scholar
- Cai M, Kim S, Wang K, Farnham PJ, Coetzee GA, Lu W. 4C-seq revealed long-range interactions of a functional enhancer at the 8q24 prostate cancer risk locus. Sci Rep. 2016;6:22462.PubMedPubMed CentralView ArticleGoogle Scholar
- Williams SR, Zies D, Mullegama SV, Grotewiel MS, Elsea SH. Smith-Magenis syndrome results in disruption of CLOCK gene transcription and reveals an integral role for RAI1 in the maintenance of circadian rhythmicity. Am J Hum Genet. 2012;90(6):941–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Potocki L, Glaze D, Tan DX, Park SS, Kashork CD, Shaffer LG, et al. Circadian rhythm abnormalities of melatonin in Smith-Magenis syndrome. J Med Genet. 2000;37(6):428–33.PubMedPubMed CentralView ArticleGoogle Scholar
- Lupianez DG, Kraft K, Heinrich V, Krawitz P, Brancati F, Klopocki E, et al. Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions. Cell. 2015;161(5):1012–25.PubMedPubMed CentralView ArticleGoogle Scholar