Cohen Syndrome (COH1) is a rare autosomal recessive disorder, principally identified by ocular, neural and muscular deficits. We identified three large consanguineous Pakistani families with intellectual disability and in some cases with autistic traits.

Autism spectrum disorder (ASD) is a highly heritable disorder of complex and heterogeneous aetiology. It is primarily characterized by altered cognitive ability including impaired language and communication skills and fundamental deficits in social reciprocity. Despite some notable successes in neuropsychiatric genetics, overall, the high heritability of ASD (~90%) remains poorly explained by common genetic risk variants. However, recent studies suggest that rare genomic variation, in particular copy number variation, may account for a significant proportion of the genetic basis of ASD. We present a large scale analysis to identify candidate genes which may contain low-frequency recessive variation contributing to ASD while taking into account the potential contribution of population differences to the genetic heterogeneity of ASD. Our strategy, homozygous haplotype (HH) mapping, aims to detect homozygous segments of identical haplotype structure that are shared at a higher frequency amongst ASD patients compared to parental controls. The analysis was performed on 1,402 Autism Genome Project trios genotyped for 1 million single nucleotide polymorphisms (SNPs). We identified 25 known and 1,218 novel ASD candidate genes in the discovery analysis including CADM2, ABHD14A, CHRFAM7A, GRIK2, GRM3, EPHA3, FGF10, KCND2, PDZK1, IMMP2L and FOXP2. Furthermore, 10 of the previously reported ASD genes and 300 of the novel candidates identified in the discovery analysis were replicated in an independent sample of 1,182 trios. Our results demonstrate that regions of HH are significantly enriched for previously reported ASD candidate genes and the observed association is independent of gene size (odds ratio 2.10). Our findings highlight the applicability of HH mapping in complex disorders such as ASD and offer an alternative approach to the analysis of genome-wide association data.

Previously, we have described the clinical and molecular characterization of a de novo 14q13.1-q21.1 microdeletion, less than 3.5 Mb in size, in a patient with severe microcephaly, psychomotor retardation, and other clinical anomalies. Here we report the characterization of the genomic structure of the human tuberin-like protein gene 1 (TULIP1; approved gene symbol GARNL1), a CpGisland-associated, brain-expressed candidate gene for the neurological findings in our patient, and its murine homologue. The human TULIP1 gene was mapped to chromosome band 14q13.2 by fluorescence in situ hybridization of BAC clone RP11-355C3 (GenBank Accession No. AL160231), containing the 3' region of the gene. TULIP1 spans about 271 kb of human genomic DNA and is divided into 41 exons. An untranscribed, processed pseudogene of TULIP1 was found on human chromosome band 9q31.1. The active locus TULIP1, encoding a predicted protein of 2036 amino acids, is expressed ubiquitously in pre- and postnatal human tissues. The murine homologue Tulip1 spans about 220 kb of mouse genomic DNA and is also divided into 41 exons, encoding a predicted protein of 2035 amino acids. No pseudogene could be found in the available mouse sequence data. Several splicing variants were found. Considering the location, expression profile, and predicted function, TULIP1 is a strong candidate for several neurological features seen in 14q deletion patients. Additionally we searched for mutations in the coding region of TULIP1 in subjects from a family with idiopathic basal ganglia calcification (IBGC; Fahr disease), previously linked to chromosome 14q. We identified two novel SNPs in the intron-exon boundaries; however, they did not segregate only with affected subjects in the predicted model of an autosomal dominant disease such as IBGC.

Joubert syndrome and related disorders (JSRD) are primarily autosomal-recessive conditions characterized by hypotonia, ataxia, abnormal eye movements, and intellectual disability with a distinctive mid-hindbrain malformation. Variable features include retinal dystrophy, cystic kidney disease, and liver fibrosis. JSRD are included in the rapidly expanding group of disorders called ciliopathies, because all six gene products implicated in JSRD (NPHP1, AHI1, CEP290, RPGRIP1L, TMEM67, and ARL13B) function in the primary cilium/basal body organelle. By using homozygosity mapping in consanguineous families, we identify loss-of-function mutations in CC2D2A in JSRD patients with and without retinal, kidney, and liver disease. CC2D2A is expressed in all fetal and adult tissues tested. In ciliated cells, we observe localization of recombinant CC2D2A at the basal body and colocalization with CEP290, whose cognate gene is mutated in multiple hereditary ciliopathies. In addition, the proteins can physically interact in vitro, as shown by yeast two-hybrid and GST pull-down experiments. A nonsense mutation in the zebrafish CC2D2A ortholog (sentinel) results in pronephric cysts, a hallmark of ciliary dysfunction analogous to human cystic kidney disease. Knockdown of cep290 function in sentinel fish results in a synergistic pronephric cyst phenotype, revealing a genetic interaction between CC2D2A and CEP290 and implicating CC2D2A in cilium/basal body function. These observations extend the genetic spectrum of JSRD and provide a model system for studying extragenic modifiers in JSRD and other ciliopathies.

We have identified a large multigenerational Austrian family displaying a novel form of X-linked recessive myopathy. Affected individuals develop an adult-onset scapulo-axio-peroneal myopathy with bent-spine syndrome characterized by specific atrophy of postural muscles along with pseudoathleticism or hypertrophy and cardiac involvement. Known X-linked myopathies were excluded by simple-tandem-repeat polymorphism (STRP) and single-nucleotide polymorphism (SNP) analysis, direct gene sequencing, and immunohistochemical analysis. STRP analysis revealed significant linkage at Xq25-q27.1. Haplotype analysis based on SNP microarray data from selected family members confirmed this linkage region on the distal arm of the X chromosome, thereby narrowing down the critical interval to 12 Mb. Sequencing of functional candidate genes led to the identification of a missense mutation within the four and a half LIM domain 1 gene (FHL1), which putatively disrupts the fourth LIM domain of the protein. Mutation screening of FHL1 in a myopathy family from the UK exhibiting an almost identical phenotype revealed a 3 bp insertion mutation within the second LIM domain. FHL1 on Xq26.3 is highly expressed in skeletal and cardiac muscles. Western-blot analysis of muscle biopsies showed a marked decrease in protein expression of FHL1 in patients, in concordance with the genetic data. In summary, we have to our knowledge characterized a new disorder, X-linked myopathy with postural muscle atrophy (XMPMA), and identified FHL1 as the causative gene. This is the first FHL protein to be identified in conjunction with a human genetic disorder and further supports the role of FHL proteins in the development and maintenance of muscle tissue. Mutation screening of FHL1 should be considered for patients with uncharacterized myopathies and cardiomyopathies.

Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders.

Nucleolin (NCL/C23; OMIM: 164035) is a major nucleolar protein that plays a critical role in multiple processes, including ribosome assembly and maturation, chromatin decondensation, and pre-rRNA transcription. Due to its diverse functions, nucleolin has frequently been implicated in pathological processes, including cancer and viral infection. We recently identified a de novo frameshifting indel mutation of NCL, p.Gly664Glufs*70, through whole-exome sequencing of autism spectrum disorder trios. Through the transfection of constructs encoding either a wild-type human nucleolin or a mutant nucleolin with the same C-terminal sequence predicted for the autism proband, and by using co-localization with the nucleophosmin (NPM; B23) protein, we have shown that the nucleolin mutation leads to mislocalization of the NCL protein from the nucleolus to the nucleoplasm. Moreover, a construct with a nonsense mutation at the same residue, p.Gly664*, shows a very similar effect on the location of the NCL protein, thus confirming the presence of a predicted nucleolar location signal in this region of the NCL protein. Real-time fluorescence recovery experiments show significant changes in the kinetics and mobility of mutant NCL protein in the nucleoplasm of HEK293Tcells. Several other studies also report de novoNCL mutations in ASD or neurodevelopmental disorders. The altered mislocalization and dynamics of mutant NCL (p.G664Glufs*70/p.G664*) may have relevance to the etiopathlogy of NCL-related ASD and other neurodevelopmental phenotypes.

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

Patched domain-containing 1 (PTCHD1) is a well-established susceptibility gene for autism spectrum disorder (ASD) and intellectual disability (ID). Previous studies have suggested that alterations in the dosage of PTCHD1 may contribute to the etiology of both ASD and ID. However, there has not yet been a thorough investigation regarding mechanisms that regulate PTCHD1 expression. We sought to characterize the Ptchd1 promoter in a mouse neuronal model, as well as to identify and validate cis regulatory elements. We defined specific regions of the Ptchd1 promoter essential for robust expression in P19-induced neurons. Evolutionarily-conserved putative transcription factor binding sites within these regions were subsequently identified. Using a pairwise comparison of chromatin accessibility between mouse forebrain and liver tissues, a candidate regulatory region, ~ 9.1 kbp downstream of the Ptchd1 stop codon was defined. This region harbours two ENCODE-predicted enhancer cis-regulatory elements. Further, using DNase footprint analysis, a putative YY1-binding motif was also identified. Genomic deletion of the entire 8 kbp downstream open chromatin region attenuated Ptchd1 transcription by over 60% in our neuronal model, corroborating its predicted regulatory function. This study provides mechanistic insights related to the expression of PTCHD1, and provides important context to interpret genetic and genomic variation at this locus which may influence neurodevelopment.

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

We have used genome-wide genotyping to identify an overlapping homozygosity-by-descent locus on chromosome 9q34.3 (MRT15) in four consanguineous families affected by nonsyndromic autosomal-recessive intellectual disability (NS-ARID) and one in which the patients show additional clinical features. Four of the families are from Pakistan, and one is from Iran. Using a combination of next-generation sequencing and Sanger sequencing, we have identified mutations in the gene MAN1B1, encoding a mannosyl oligosaccharide, alpha 1,2-mannosidase. In one Pakistani family, MR43, a homozygous nonsense mutation (RefSeq number NM_016219.3: c.1418G>A [p.Trp473*]), segregated with intellectual disability and additional dysmorphic features. We also identified the missense mutation c. 1189G>A (p.Glu397Lys; RefSeq number NM_016219.3), which segregates with NS-ARID in three families who come from the same village and probably have shared inheritance. In the Iranian family, the missense mutation c.1000C>T (p.Arg334Cys; RefSeq number NM_016219.3) also segregates with NS-ARID. Both missense mutations are at amino acid residues that are conserved across the animal kingdom, and they either reduce k(cat) by ∼1300-fold or disrupt stable protein expression in mammalian cells. MAN1B1 is one of the few NS-ARID genes with an elevated mutation frequency in patients with NS-ARID from different populations.

Genome control is operated by transcription factors (TFs) controlling their target genes by binding to promoters and enhancers. Conceptually, the interactions between TFs, their binding sites, and their functional targets are represented by gene regulatory networks (GRNs). Deciphering in vivo GRNs underlying organ development in an unbiased genome-wide setting involves identifying both functional TF-gene interactions and physical TF-DNA interactions. To reverse engineer the GRNs of eye development in Drosophila, we performed RNA-seq across 72 genetic perturbations and sorted cell types and inferred a coexpression network. Next, we derived direct TF-DNA interactions using computational motif inference, ultimately connecting 241 TFs to 5,632 direct target genes through 24,926 enhancers. Using this network, we found network motifs, cis-regulatory codes, and regulators of eye development. We validate the predicted target regions of Grainyhead by ChIP-seq and identify this factor as a general cofactor in the eye network, being bound to thousands of nucleosome-free regions.

The functionality of sensory neurons is defined by the expression of specific sensory receptor genes. During the development of the Drosophila larval eye, photoreceptor neurons (PRs) make a binary choice to express either the blue-sensitive Rhodopsin 5 (Rh5) or the green-sensitive Rhodopsin 6 (Rh6). Later during metamorphosis, ecdysone signaling induces a cell fate and sensory receptor switch: Rh5-PRs are re-programmed to express Rh6 and become the eyelet, a small group of extraretinal PRs involved in circadian entrainment. However, the genetic and molecular mechanisms of how the binary cell fate decisions are made and switched remain poorly understood. We show that interplay of two transcription factors Senseless (Sens) and Hazy control cell fate decisions, terminal differentiation of the larval eye and its transformation into eyelet. During initial differentiation, a pulse of Sens expression in primary precursors regulates their differentiation into Rh5-PRs and repression of an alternative Rh6-cell fate. Later, during the transformation of the larval eye into the adult eyelet, Sens serves as an anti-apoptotic factor in Rh5-PRs, which helps in promoting survival of Rh5-PRs during metamorphosis and is subsequently required for Rh6 expression. Comparably, during PR differentiation Hazy functions in initiation and maintenance of rhodopsin expression. Hazy represses Sens specifically in the Rh6-PRs, allowing them to die during metamorphosis. Our findings show that the same transcription factors regulate diverse aspects of larval and adult PR development at different stages and in a context-dependent manner.

The major mood disorders, which include bipolar disorder and major depressive disorder (MDD), are considered heritable traits, although previous genetic association studies have had limited success in robustly identifying risk loci. We performed a meta-analysis of five case-control cohorts for major mood disorder, including over 13,600 individuals genotyped on high-density SNP arrays. We identified SNPs at 3p21.1 associated with major mood disorders (rs2251219, P = 3.63 x 10(-8); odds ratio = 0.87; 95% confidence interval, 0.83-0.92), with supportive evidence for association observed in two out of three independent replication cohorts. These results provide an example of a shared genetic susceptibility locus for bipolar disorder and MDD.

Structural variation (copy number variation [CNV] including deletion and duplication, translocation, inversion) of chromosomes has been identified in some individuals with autism spectrum disorder (ASD), but the full etiologic role is unknown. We performed genome-wide assessment for structural abnormalities in 427 unrelated ASD cases via single-nucleotide polymorphism microarrays and karyotyping. With microarrays, we discovered 277 unbalanced CNVs in 44% of ASD families not present in 500 controls (and re-examined in another 1152 controls). Karyotyping detected additional balanced changes. Although most variants were inherited, we found a total of 27 cases with de novo alterations, and in three (11%) of these individuals, two or more new variants were observed. De novo CNVs were found in approximately 7% and approximately 2% of idiopathic families having one child, or two or more ASD siblings, respectively. We also detected 13 loci with recurrent/overlapping CNV in unrelated cases, and at these sites, deletions and duplications affecting the same gene(s) in different individuals and sometimes in asymptomatic carriers were also found. Notwithstanding complexities, our results further implicate the SHANK3-NLGN4-NRXN1 postsynaptic density genes and also identify novel loci at DPP6-DPP10-PCDH9 (synapse complex), ANKRD11, DPYD, PTCHD1, 15q24, among others, for a role in ASD susceptibility. Our most compelling result discovered CNV at 16p11.2 (p = 0.002) (with characteristics of a genomic disorder) at approximately 1% frequency. Some of the ASD regions were also common to mental retardation loci. Structural variants were found in sufficiently high frequency influencing ASD to suggest that cytogenetic and microarray analyses be considered in routine clinical workup.

Autosomal-recessive inheritance is believed to be relatively common in mental retardation (MR), although only four genes for nonsyndromic autosomal-recessive mental retardation (ARMR) have been reported. In this study, we ascertained a consanguineous Pakistani family with ARMR in four living individuals from three branches of the family, plus an additional affected individual later identified as a phenocopy. Retinitis pigmentosa was present in affected individuals, but no other features suggestive of a syndromic form of MR were found. We used Affymetrix 500K microarrays to perform homozygosity mapping and identified a homozygous and haploidentical region of 11.2 Mb on chromosome 4p15.33-p15.2. Linkage analysis across this region produced a maximum two-point LOD score of 3.59. We sequenced genes within the critical region and identified a homozygous splice-site mutation segregating in the family, within a coiled-coil and C2 domain-containing gene, CC2D2A. This mutation leads to the skipping of exon 19, resulting in a frameshift and a truncated protein lacking the C2 domain. Conservation analysis for CC2D2A suggests a functional domain near the C terminus as well as the C2 domain. Preliminary functional studies of CC2D2A suggest a possible role in Ca(2+)-dependent signal transduction. Identifying the function of CC2D2A, and a possible common pathway with CC2D1A, in correct neuronal development and functioning may help identify possible therapeutic targets for MR.

Trivalent chromium (Cr) and bitter melon (Momordica charantia L., BM) have been shown to independently interact with the insulin signaling pathway leading to improvements in the symptoms of insulin resistance and diabetes in some animal models and human subjects. The aim of this study was to examine whether the combination of the two nutritional supplements could potentially have additive effects on treating these conditions in high-fat-fed streptozotocin (STZ)-induced diabetic rats. The experiment was conducted with 110 male Wistar rats divided into eleven groups and fed either a control or high-fat diet for 7 weeks. Half of the rats on the high-fat diet were injected with STZ (30 mg/kg body mass) to induce diabetes. The high-fat (HF) diets were then supplemented with a combination of Cr (as chromium(III) propionate complex, Cr3: either 10 or 50 mg Cr/kg diet) and bitter melon (lyophilized whole fruit: either 10 or 50 g/kg diet) for 6 weeks. After termination of the experiment, blood and internal organs were harvested for blood biochemical, hematological, and mineral (Cr) analyses using appropriate analytical methods. It was found that neither Cr(III) nor BM was able to significantly affect blood indices in HF and diabetic rats, but BM tended to improve body mass gain, blood glucose, and LDL cholesterol values, but decreased Cr content in the liver and kidneys of the Cr-co-supplemented type 2 diabetic model of rats. Supplementary Cr(III) had no appreciable effect on glucose and lipid metabolism in high-fat-fed STZ-induced diabetic rats. Supplementary BM fruit powder had some observable effects on body mass of high-fat-fed rats; these effects seem to be dampened when BM was co-administered with Cr. Cr(III) and BM appear to act as nutritional antagonists when both administered in food, probably due to binding of Cr by the polyphenol-type compounds present in the plant material. Graphical Abstract.

The risk of epilepsy among individuals with intellectual disability (ID) is approximately ten times that of the general population. From a cohort of >5,000 families affected by neurodevelopmental disorders, we identified six consanguineous families harboring homozygous inactivating variants in MBOAT7, encoding lysophosphatidylinositol acyltransferase (LPIAT1). Subjects presented with ID frequently accompanied by epilepsy and autistic features. LPIAT1 is a membrane-bound phospholipid-remodeling enzyme that transfers arachidonic acid (AA) to lysophosphatidylinositol to produce AA-containing phosphatidylinositol. This study suggests a role for AA-containing phosphatidylinositols in the development of ID accompanied by epilepsy and autistic features.

PTCHD1 has been implicated in Autism Spectrum Disorders (ASDs) and/or intellectual disability, where copy-number-variant losses or loss-of-function coding mutations segregate with disease in an X-linked recessive fashion. Missense variants of PTCHD1 have also been reported in patients. However, the significance of these mutations remains undetermined since the activities, subcellular localization, and regulation of the PTCHD1 protein are currently unknown. This paucity of data concerning PTCHD1 prevents the effective evaluation of sequence variants identified during diagnostic screening. Here, we characterize PTCHD1 protein binding partners, extending previously reported interactions with postsynaptic scaffolding protein, SAP102. Six rare missense variants of PTCHD1 were also identified from patients with neurodevelopmental disorders. After modelling these variants on a hypothetical three-dimensional structure of PTCHD1, based on the solved structure of NPC1, PTCHD1 variants harboring these mutations were assessed for protein stability, post-translational processing, and protein trafficking. We show here that the wild-type PTCHD1 post-translational modification includes complex N-glycosylation and that specific mutant proteins disrupt normal N-link glycosylation processing. However, regardless of their processing, these mutants still localized to PSD95-containing dendritic processes and remained competent for complexing SAP102.

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 22 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, YWHAE, DPH1, GSDMB, MED24, THRA, EEF1A2, and KCNQ2 in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance and transferability of BD polygenic risk scores across ancestrally diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).