Hemizygous variants in protein phosphatase 1 regulatory subunit 3F (PPP1R3F) are associated with a neurodevelopmental disorder characterized by developmental delay, intellectual disability and autistic features.

Protein phosphatase 1 regulatory subunit 3F (PPP1R3F) is a member of the glycogen targeting subunits (GTSs), which belong to the large group of regulatory subunits of protein phosphatase 1 (PP1), a major eukaryotic serine/threonine protein phosphatase that regulates diverse cellular processes. Here, we describe the identification of hemizygous variants in PPP1R3F associated with a novel X-linked recessive neurodevelopmental disorder in 13 unrelated individuals. This disorder is characterized by developmental delay, mild intellectual disability, neurobehavioral issues such as autism spectrum disorder, seizures and other neurological findings including tone, gait and cerebellar abnormalities. PPP1R3F variants segregated with disease in affected hemizygous males that inherited the variants from their heterozygous carrier mothers. We show that PPP1R3F is predominantly expressed in brain astrocytes and localizes to the endoplasmic reticulum in cells. Glycogen content in PPP1R3F knockout astrocytoma cells appears to be more sensitive to fluxes in extracellular glucose levels than in wild-type cells, suggesting that PPP1R3F functions in maintaining steady brain glycogen levels under changing glucose conditions. We performed functional studies on nine of the identified variants and observed defects in PP1 binding, protein stability, subcellular localization and regulation of glycogen metabolism in most of them. Collectively, the genetic and molecular data indicate that deleterious variants in PPP1R3F are associated with a new X-linked disorder of glycogen metabolism, highlighting the critical role of GTSs in neurological development. This research expands our understanding of neurodevelopmental disorders and the role of PP1 in brain development and proper function.

Noncoding deletions reveal a gene that is critical for intestinal function.

Large-scale genome sequencing is poised to provide a substantial increase in the rate of discovery of disease-associated mutations, but the functional interpretation of such mutations remains challenging. Here we show that deletions of a sequence on human chromosome 16 that we term the intestine-critical region (ICR) cause intractable congenital diarrhoea in infants1,2. Reporter assays in transgenic mice show that the ICR contains a regulatory sequence that activates transcription during the development of the gastrointestinal system. Targeted deletion of the ICR in mice caused symptoms that recapitulated the human condition. Transcriptome analysis revealed that an unannotated open reading frame (Percc1) flanks the regulatory sequence, and the expression of this gene was lost in the developing gut of mice that lacked the ICR. Percc1-knockout mice displayed phenotypes similar to those observed upon ICR deletion in mice and patients, whereas an ICR-driven Percc1 transgene was sufficient to rescue the phenotypes found in mice that lacked the ICR. Together, our results identify a gene that is critical for intestinal function and underscore the need for targeted in vivo studies to interpret the growing number of clinical genetic findings that do not affect known protein-coding genes.

The best of both worlds: Blending cutting-edge research with clinical processes for a productive exome clinic.

Genomic medicine has been transformed by next-generation sequencing (NGS), inclusive of exome sequencing (ES) and genome sequencing (GS). Currently, ES is offered widely in clinical settings, with a less prevalent alternative model consisting of hybrid programs that incorporate research ES along with clinical patient workflows. We were among the earliest to implement a hybrid ES clinic, have provided diagnoses to 45% of probands, and have identified several novel candidate genes. Our program is enabled by a cost-effective investment by the health system and is unique in encompassing all the processes that have been variably included in other hybrid/clinical programs. These include careful patient selection, utilization of a phenotype-agnostic bioinformatics pipeline followed by manual curation of variants and phenotype integration by clinicians, close collaborations between the clinicians and the bioinformatician, pursuit of interesting variants, communication of results to patients in categories that are predicated upon the certainty of a diagnosis, and tracking changes in results over time and the underlying mechanisms for such changes. Due to its effectiveness, scalability to GS and its resource efficiency, specific elements of our paradigm can be incorporated into existing clinical settings, or the entire hybrid model can be implemented within health systems that have genomic medicine programs, to provide NGS in a scientifically rigorous, yet pragmatic setting.

High-impact rare genetic variants in severe schizophrenia.

Extreme phenotype sequencing has led to the identification of high-impact rare genetic variants for many complex disorders but has not been applied to studies of severe schizophrenia. We sequenced 112 individuals with severe, extremely treatment-resistant schizophrenia, 218 individuals with typical schizophrenia, and 4,929 controls. We compared the burden of rare, damaging missense and loss-of-function variants between severe, extremely treatment-resistant schizophrenia, typical schizophrenia, and controls across mutation intolerant genes. Individuals with severe, extremely treatment-resistant schizophrenia had a high burden of rare loss-of-function (odds ratio, 1.91; 95% CI, 1.39 to 2.63; P = 7.8 × 10-5) and damaging missense variants in intolerant genes (odds ratio, 2.90; 95% CI, 2.02 to 4.15; P = 3.2 × 10-9). A total of 48.2% of individuals with severe, extremely treatment-resistant schizophrenia carried at least one rare, damaging missense or loss-of-function variant in intolerant genes compared to 29.8% of typical schizophrenia individuals (odds ratio, 2.18; 95% CI, 1.33 to 3.60; P = 1.6 × 10-3) and 25.4% of controls (odds ratio, 2.74; 95% CI, 1.85 to 4.06; P = 2.9 × 10-7). Restricting to genes previously associated with schizophrenia risk strengthened the enrichment with 8.9% of individuals with severe, extremely treatment-resistant schizophrenia carrying a damaging missense or loss-of-function variant compared to 2.3% of typical schizophrenia (odds ratio, 5.48; 95% CI, 1.52 to 19.74; P = 0.02) and 1.6% of controls (odds ratio, 5.82; 95% CI, 3.00 to 11.28; P = 2.6 × 10-8). These results demonstrate the power of extreme phenotype case selection in psychiatric genetics and an approach to augment schizophrenia gene discovery efforts.

Variants in SCAF4 Cause a Neurodevelopmental Disorder and Are Associated with Impaired mRNA Processing.

RNA polymerase II interacts with various other complexes and factors to ensure correct initiation, elongation, and termination of mRNA transcription. One of these proteins is SR-related CTD-associated factor 4 (SCAF4), which is important for correct usage of polyA sites for mRNA termination. Using exome sequencing and international matchmaking, we identified nine likely pathogenic germline variants in SCAF4 including two splice-site and seven truncating variants, all residing in the N-terminal two thirds of the protein. Eight of these variants occurred de novo, and one was inherited. Affected individuals demonstrated a variable neurodevelopmental disorder characterized by mild intellectual disability, seizures, behavioral abnormalities, and various skeletal and structural anomalies. Paired-end RNA sequencing on blood lymphocytes of SCAF4-deficient individuals revealed a broad deregulation of more than 9,000 genes and significant differential splicing of more than 2,900 genes, indicating an important role of SCAF4 in mRNA processing. Knockdown of the SCAF4 ortholog CG4266 in the model organism Drosophila melanogaster resulted in impaired locomotor function, learning, and short-term memory. Furthermore, we observed an increased number of active zones in larval neuromuscular junctions, representing large glutamatergic synapses. These observations indicate a role of CG4266 in nervous system development and function and support the implication of SCAF4 in neurodevelopmental phenotypes. In summary, our data show that heterozygous, likely gene-disrupting variants in SCAF4 are causative for a variable neurodevelopmental disorder associated with impaired mRNA processing.

A concurrent dual analysis of genomic data augments diagnoses: Experiences of 2 clinical sites in the Undiagnosed Diseases Network.

PURPOSE: Next-generation sequencing (NGS) has revolutionized the diagnostic process for rare/ultrarare conditions. However, diagnosis rates differ between analytical pipelines. In the National Institutes of Health-Undiagnosed Diseases Network (UDN) study, each individual's NGS data are concurrently analyzed by the UDN sequencing core laboratory and the clinical sites. We examined the outcomes of this practice. METHODS: A retrospective review was performed at 2 UDN clinical sites to compare the variants and diagnoses/candidate genes identified with the dual analyses of the NGS data. RESULTS: In total, 95 individuals had 100 diagnoses/candidate genes. There was 59% concordance between the UDN sequencing core laboratories and the clinical sites in identifying diagnoses/candidate genes. The core laboratory provided more diagnoses, whereas the clinical sites prioritized more research variants/candidate genes (P < .001). The clinical sites solely identified 15% of the diagnoses/candidate genes. The differences between the 2 pipelines were more often because of variant prioritization disparities than variant detection. CONCLUSION: The unique dual analysis of NGS data in the UDN synergistically enhances outcomes. The core laboratory provided a clinical analysis with more diagnoses and the clinical sites prioritized more research variants/candidate genes. Implementing such concurrent dual analyses in other genomic research studies and clinical settings can improve both variant detection and prioritization.

The benefit of diagnostic whole genome sequencing in schizophrenia and other psychotic disorders.

Schizophrenia has a multifactorial etiology, involving a polygenic architecture. The potential benefit of whole genome sequencing (WGS) in schizophrenia and other psychotic disorders is not well studied. We investigated the yield of clinical WGS analysis in 251 families with a proband diagnosed with schizophrenia (N = 190), schizoaffective disorder (N = 49), or other conditions involving psychosis (N = 48). Participants were recruited in Israel and USA, mainly of Jewish, Arab, and other European ancestries. Trio (parents and proband) WGS was performed for 228 families (90.8%); in the other families, WGS included parents and at least two affected siblings. In the secondary analyses, we evaluated the contribution of rare variant enrichment in particular gene sets, and calculated polygenic risk score (PRS) for schizophrenia. For the primary outcome, diagnostic rate was 6.4%; we found clinically significant, single nucleotide variants (SNVs) or small insertions or deletions (indels) in 14 probands (5.6%), and copy number variants (CNVs) in 2 (0.8%). Significant enrichment of rare loss-of-function variants was observed in a gene set of top schizophrenia candidate genes in affected individuals, compared with population controls (N = 6,840). The PRS for schizophrenia was significantly increased in the affected individuals group, compared to their unaffected relatives. Last, we were also able to provide pharmacogenomics information based on CYP2D6 genotype data for most participants, and determine their antipsychotic metabolizer status. In conclusion, our findings suggest that WGS may have a role in the setting of both research and genetic counseling for individuals with schizophrenia and other psychotic disorders and their families.

Investigation into the genetics of fetal congenital lymphatic anomalies.

OBJECTIVE: Congenital lymphatic anomalies (LAs) arise due to defects in lymphatic development and often present in utero as pleural effusion, chylothorax, nuchal and soft tissue edema, ascites, or hydrops. Many LAs are caused by single nucleotide variants, which are not detected on routine prenatal testing. METHODS: Demographic data were compared between two subcohorts, those with clinically significant fetal edema (CSFE) and isolated fetal edema. A targeted variant analysis of LA genes was performed using American College of Medical Genetics criteria on whole exome sequencing (WES) data generated for 71 fetal edema cases who remained undiagnosed after standard workup. RESULTS: CSFE cases had poor outcomes, including preterm delivery, demise, and maternal preeclampsia. Pathogenic and likely pathogenic variants were identified in 7% (5/71) of cases, including variants in RASopathy genes, RASA1, SOS1, PTPN11, and a novel PIEZO1 variant. Variants of uncertain significance (VOUS) were identified in 45% (32/71) of cases. In CSFEs, VOUS were found in CELSR1, EPHB4, TIE1, PIEZO1, ITGA9, RASopathy genes, SOS1, SOS2, and RAF1. CONCLUSIONS: WES identified pathogenic and likely pathogenic variants and VOUS in LA genes in 51% of fetal edema cases, supporting WES and expanded hydrops panels in cases of idiopathic fetal hydrops and fluid collections.

Commonalities across computational workflows for uncovering explanatory variants in undiagnosed cases.

PURPOSE: Genomic sequencing has become an increasingly powerful and relevant tool to be leveraged for the discovery of genetic aberrations underlying rare, Mendelian conditions. Although the computational tools incorporated into diagnostic workflows for this task are continually evolving and improving, we nevertheless sought to investigate commonalities across sequencing processing workflows to reveal consensus and standard practice tools and highlight exploratory analyses where technical and theoretical method improvements would be most impactful. METHODS: We collected details regarding the computational approaches used by a genetic testing laboratory and 11 clinical research sites in the United States participating in the Undiagnosed Diseases Network via meetings with bioinformaticians, online survey forms, and analyses of internal protocols. RESULTS: We found that tools for processing genomic sequencing data can be grouped into four distinct categories. Whereas well-established practices exist for initial variant calling and quality control steps, there is substantial divergence across sites in later stages for variant prioritization and multimodal data integration, demonstrating a diversity of approaches for solving the most mysterious undiagnosed cases. CONCLUSION: The largest differences across diagnostic workflows suggest that advances in structural variant detection, noncoding variant interpretation, and integration of additional biomedical data may be especially promising for solving chronically undiagnosed cases.

Phenotypic expansion of CACNA1C-associated disorders to include isolated neurological manifestations.

PURPOSE: CACNA1C encodes the alpha-1-subunit of a voltage-dependent L-type calcium channel expressed in human heart and brain. Heterozygous variants in CACNA1C have previously been reported in association with Timothy syndrome and long QT syndrome. Several case reports have suggested that CACNA1C variation may also be associated with a primarily neurological phenotype. METHODS: We describe 25 individuals from 22 families with heterozygous variants in CACNA1C, who present with predominantly neurological manifestations. RESULTS: Fourteen individuals have de novo, nontruncating variants and present variably with developmental delays, intellectual disability, autism, hypotonia, ataxia, and epilepsy. Functional studies of a subgroup of missense variants via patch clamp experiments demonstrated differential effects on channel function in vitro, including loss of function (p.Leu1408Val), neutral effect (p.Leu614Arg), and gain of function (p.Leu657Phe, p.Leu614Pro). The remaining 11 individuals from eight families have truncating variants in CACNA1C. The majority of these individuals have expressive language deficits, and half have autism. CONCLUSION: We expand the phenotype associated with CACNA1C variants to include neurodevelopmental abnormalities and epilepsy, in the absence of classic features of Timothy syndrome or long QT syndrome.

Expansion of the GRIA2 phenotypic representation: a novel de novo loss of function mutation in a case with childhood onset schizophrenia.

Childhood-onset schizophrenia (COS) is a rare form of schizophrenia with an onset before 13 years of age. There is rising evidence that genetic factors play a major role in COS etiology, yet, only a few single gene mutations have been discovered. Here we present a diagnostic whole-exome sequencing (WES) in an Israeli Jewish female with COS and additional neuropsychiatric conditions such as obsessive-compulsive disorder (OCD), anxiety, and aggressive behavior. Variant analysis revealed a de novo novel stop gained variant in GRIA2 gene (NM_000826.4: c.1522 G > T (p.Glu508Ter)). GRIA2 encodes for a subunit of the AMPA sensitive glutamate receptor (GluA2) that functions as ligand-gated ion channel in the central nervous system and plays an important role in excitatory synaptic transmission. GluA2 subunit mutations are known to cause variable neurodevelopmental phenotypes including intellectual disability, autism spectrum disorder, epilepsy, and OCD. Our findings support the potential diagnostic role of WES in COS, identify GRIA2 as possible cause to a broad psychiatric phenotype that includes COS as a major manifestation and expand the previously reported GRIA2 loss of function phenotypes.

Identification of a Functional Risk Variant for Pemphigus Vulgaris in the ST18 Gene.

Pemphigus vulgaris (PV) is a life-threatening autoimmune mucocutaneous blistering disease caused by disruption of intercellular adhesion due to auto-antibodies directed against epithelial components. Treatment is limited to immunosuppressive agents, which are associated with serious adverse effects. The propensity to develop the disease is in part genetically determined. We therefore reasoned that the delineation of PV genetic basis may point to novel therapeutic strategies. Using a genome-wide association approach, we recently found that genetic variants in the vicinity of the ST18 gene confer a significant risk for the disease. Here, using targeted deep sequencing, we identified a PV-associated variant residing within the ST18 promoter region (p<0.0002; odds ratio = 2.03). This variant was found to drive increased gene transcription in a p53/p63-dependent manner, which may explain the fact that ST18 is up-regulated in the skin of PV patients. We then discovered that when overexpressed, ST18 stimulates PV serum-induced secretion of key inflammatory molecules and contributes to PV serum-induced disruption of keratinocyte cell-cell adhesion, two processes previously implicated in the pathogenesis of PV. Thus, the present findings indicate that ST18 may play a direct role in PV and consequently represents a potential target for the treatment of this disease.

The human olfactory transcriptome.

BACKGROUND: Olfaction is a versatile sensory mechanism for detecting thousands of volatile odorants. Although molecular basis of odorant signaling is relatively well understood considerable gaps remain in the complete charting of all relevant gene products. To address this challenge, we applied RNAseq to four well-characterized human olfactory epithelial samples and compared the results to novel and published mouse olfactory epithelium as well as 16 human control tissues. RESULTS: We identified 194 non-olfactory receptor (OR) genes that are overexpressed in human olfactory tissues vs. CONTROLS: The highest overexpression is seen for lipocalins and bactericidal/permeability-increasing (BPI)-fold proteins, which in other species include secreted odorant carriers. Mouse-human discordance in orthologous lipocalin expression suggests different mammalian evolutionary paths in this family. Of the overexpressed genes 36 have documented olfactory function while for 158 there is little or no previous such functional evidence. The latter group includes GPCRs, neuropeptides, solute carriers, transcription factors and biotransformation enzymes. Many of them may be indirectly implicated in sensory function, and ~70 % are over expressed also in mouse olfactory epithelium, corroborating their olfactory role. Nearly 90 % of the intact OR repertoire, and ~60 % of the OR pseudogenes are expressed in the olfactory epithelium, with the latter showing a 3-fold lower expression. ORs transcription levels show a 1000-fold inter-paralog variation, as well as significant inter-individual differences. We assembled 160 transcripts representing 100 intact OR genes. These include 1-4 short 5' non-coding exons with considerable alternative splicing and long last exons that contain the coding region and 3' untranslated region of highly variable length. Notably, we identified 10 ORs with an intact open reading frame but with seemingly non-functional transcripts, suggesting a yet unreported OR pseudogenization mechanism. Analysis of the OR upstream regions indicated an enrichment of the homeobox family transcription factor binding sites and a consensus localization of a specific transcription factor binding site subfamily (Olf/EBF). CONCLUSIONS: We provide an overview of expression levels of ORs and auxiliary genes in human olfactory epithelium. This forms a transcriptomic view of the entire OR repertoire, and reveals a large number of over-expressed uncharacterized human non-receptor genes, providing a platform for future discovery.

VarElect: the phenotype-based variation prioritizer of the GeneCards Suite.

BACKGROUND: Next generation sequencing (NGS) provides a key technology for deciphering the genetic underpinnings of human diseases. Typical NGS analyses of a patient depict tens of thousands non-reference coding variants, but only one or very few are expected to be significant for the relevant disorder. In a filtering stage, one employs family segregation, rarity in the population, predicted protein impact and evolutionary conservation as a means for shortening the variation list. However, narrowing down further towards culprit disease genes usually entails laborious seeking of gene-phenotype relationships, consulting numerous separate databases. Thus, a major challenge is to transition from the few hundred shortlisted genes to the most viable disease-causing candidates. RESULTS: We describe a novel tool, VarElect ( http://ve.genecards.org ), a comprehensive phenotype-dependent variant/gene prioritizer, based on the widely-used GeneCards, which helps rapidly identify causal mutations with extensive evidence. The GeneCards suite offers an effective and speedy alternative, whereby >120 gene-centric automatically-mined data sources are jointly available for the task. VarElect cashes on this wealth of information, as well as on GeneCards' powerful free-text Boolean search and scoring capabilities, proficiently matching variant-containing genes to submitted disease/symptom keywords. The tool also leverages the rich disease and pathway information of MalaCards, the human disease database, and PathCards, the unified pathway (SuperPaths) database, both within the GeneCards Suite. The VarElect algorithm infers direct as well as indirect links between genes and phenotypes, the latter benefitting from GeneCards' diverse gene-to-gene data links in GenesLikeMe. Finally, our tool offers an extensive gene-phenotype evidence portrayal ("MiniCards") and hyperlinks to the parent databases. CONCLUSIONS: We demonstrate that VarElect compares favorably with several often-used NGS phenotyping tools, thus providing a robust facility for ranking genes, pointing out their likelihood to be related to a patient's disease. VarElect's capacity to automatically process numerous NGS cases, either in stand-alone format or in VCF-analyzer mode (TGex and VarAnnot), is indispensable for emerging clinical projects that involve thousands of whole exome/genome NGS analyses.

Mutation in TECPR2 reveals a role for autophagy in hereditary spastic paraparesis.

We studied five individuals from three Jewish Bukharian families affected by an apparently autosomal-recessive form of hereditary spastic paraparesis accompanied by severe intellectual disability, fluctuating central hypoventilation, gastresophageal reflux disease, wake apnea, areflexia, and unique dysmorphic features. Exome sequencing identified one homozygous variant shared among all affected individuals and absent in controls: a 1 bp frameshift TECPR2 deletion leading to a premature stop codon and predicting significant degradation of the protein. TECPR2 has been reported as a positive regulator of autophagy. We thus examined the autophagy-related fate of two key autophagic proteins, SQSTM1 (p62) and MAP1LC3B (LC3), in skin fibroblasts of an affected individual, as compared to a healthy control, and found that both protein levels were decreased and that there was a more pronounced decrease in the lipidated form of LC3 (LC3II). siRNA knockdown of TECPR2 showed similar changes, consistent with aberrant autophagy. Our results are strengthened by the fact that autophagy dysfunction has been implicated in a number of other neurodegenerative diseases. The discovered TECPR2 mutation implicates autophagy, a central intracellular mechanism, in spastic paraparesis.

Massive Parallel DNA Sequencing of Patients with Inherited Cardiomyopathies in Cyprus and Suggestion of Digenic or Oligogenic Inheritance.

Inherited cardiomyopathies represent a highly heterogeneous group of cardiac diseases. DNA variants in genes expressed in cardiomyocytes cause a diverse spectrum of cardiomyopathies, ultimately leading to heart failure, arrythmias, and sudden cardiac death. We applied massive parallel DNA sequencing using a 72-gene panel for studying inherited cardiomyopathies. We report on variants in 25 families, where pathogenicity was predicted by different computational approaches, databases, and an in-house filtering analysis. All variants were validated using Sanger sequencing. Familial segregation was tested when possible. We identified 41 different variants in 26 genes. Analytically, we identified fifteen variants previously reported in the Human Gene Mutation Database: twelve mentioned as disease-causing mutations (DM) and three as probable disease-causing mutations (DM?). Additionally, we identified 26 novel variants. We classified the forty-one variants as follows: twenty-eight (68.3%) as variants of uncertain significance, eight (19.5%) as likely pathogenic, and five (12.2%) as pathogenic. We genetically characterized families with a cardiac phenotype. The genetic heterogeneity and the multiplicity of candidate variants are making a definite molecular diagnosis challenging, especially when there is a suspicion of incomplete penetrance or digenic-oligogenic inheritance. This is the first systematic study of inherited cardiac conditions in Cyprus, enabling us to develop a genetic baseline and precision cardiology.

Genetic insights into childhood-onset schizophrenia: The yield of clinical exome sequencing.

Childhood-onset schizophrenia (COS) is a rare form of schizophrenia with an onset prior to 13 years of age. Although genetic factors play a role in COS etiology, only a few causal variants have been reported to date. This study presents a diagnostic exome sequencing (ES) in 37 Israeli Jewish families with a proband diagnosed with COS. By implementing a trio/duo ES approach and applying a well-established diagnostic pipeline, we detected clinically significant variants in 7 probands (19 %). These single nucleotide variants and indels were mostly inherited. The implicated genes were ANKRD11, GRIA2, CHD2, CLCN3, CLTC, IGF1R and MICU1. In a secondary analysis that compared COS patients to 4721 healthy controls, we observed that patients had a significant enrichment of rare loss of function (LoF) variants in LoF intolerant genes associated with developmental diseases. Taken together, ES could be considered as a valuable tool in the genetic workup for COS patients.

Comprehensive Genetic Analysis of Druze Provides Insights into Carrier Screening.

BACKGROUND: Druze individuals, like many genetically homogeneous and isolated populations, harbor recurring pathogenic variants (PV) in autosomal recessive (AR) disorders. METHODS: Variant calling of whole-genome sequencing (WGS) of 40 Druze from the Human Genome Diversity Project (HGDP) was performed (HGDP-cohort). Additionally, we performed whole exome sequencing (WES) of 118 Druze individuals: 38 trios and 2 couples, representing geographically distinct clans (WES-cohort). Rates of validated PV were compared with rates in worldwide and Middle Eastern populations, from the gnomAD and dbSNP datasets. RESULTS: Overall, 34 PVs were identified: 30 PVs in genes underlying AR disorders, 3 additional PVs were associated with autosomal dominant (AD) disorders, and 1 PV with X-linked-dominant inherited disorder in the WES cohort. CONCLUSIONS: The newly identified PVs associated with AR conditions should be considered for incorporation into prenatal-screening options offered to Druze individuals after an extension and validation of the results in a larger study.

DOCK4 and CEACAM21 as novel schizophrenia candidate genes in the Jewish population.

It is well accepted that schizophrenia has a strong genetic component. Several genome-wide association studies (GWASs) of schizophrenia have been published in recent years; most of them population based with a case-control design. Nevertheless, identifying the specific genetic variants which contribute to susceptibility to the disorder remains a challenging task. A family-based GWAS strategy may be helpful in the identification of schizophrenia susceptibility genes since it is protected against population stratification, enables better accounting for genotyping errors and is more sensitive for identification of rare variants which have a very low frequency in the general population. In this project we implemented a family-based GWAS of schizophrenia in a sample of 107 Jewish-Israeli families. We found one genome-wide significant association in the intron of the DOCK4 gene (rs2074127, p value=1.134×10⁻7) and six additional nominally significant association signals with p<1×10⁻5. One of the top single nucleotide polymorphisms (p<1×10⁻5) which is located in the predicted intron of the CEACAM21 gene was significantly replicated in independent family-based sample of Arab-Israeli origin (rs4803480: p value=0.002; combined p value=9.61×10⁻8), surviving correction for multiple testing. Both DOCK4 and CEACAM21 are biologically reasonable candidate genes for schizophrenia although generalizability of the association of DOCK4 with schizophrenia should be investigated in further studies. In addition, gene-wide significant associations were found within three schizophrenia candidate genes: PGBD1, RELN and PRODH, replicating previously reported associations. By application of a family-based strategy to GWAS, our study revealed new schizophrenia susceptibility loci in the Jewish-Israeli population.