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.

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.

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.

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.

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.

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.

Risk Variants in the Exomes of Children With Critical Illness

IMPORTANCE

Diagnostic genetic testing can lead to changes in management in the pediatric intensive care unit. Genetic risk in children with critical illness but nondiagnostic exome sequencing (ES) has not been explored.

OBJECTIVE

To assess the association between loss-of-function (LOF) variants and pediatric critical illness.

DESIGN, SETTING, AND PARTICIPANTS

This genetic association study examined ES first screened for causative variants among 267 children at the Morgan Stanley Children's Hospital of NewYork-Presbyterian, of whom 22 were otherwise healthy with viral respiratory failure; 18 deceased children with bronchiolitis from the Office of the Chief Medical Examiner of New York City, of whom 14 were previously healthy; and 9990 controls from the Institute for Genomic Medicine at Columbia University Irving Medical Center. The ES data were generated between January 1, 2015, and December 31, 2020, and analyzed between January 1, 2017, and September 2, 2022.

EXPOSURE

Critical illness.

MAIN OUTCOMES AND MEASURES

Odds ratios and P values for genes and gene-sets enriched for rare LOF variants and the loss-of-function observed/expected upper bound fraction (LOEUF) score at which cases have a significant enrichment.

RESULTS

This study included 285 children with critical illness (median [range] age, 4.1 [0-18.9] years; 148 [52%] male) and 9990 controls. A total of 228 children (80%) did not receive a genetic diagnosis. After quality control (QC), 231 children harbored excess rare LOF variants in genes with a LOEUF score of 0.680 or less (intolerant genes) (P = 1.0 × 10-5). After QC, 176 children without a diagnosis harbored excess ultrarare LOF variants in intolerant genes but only in those without a known disease association (odds ratio, 1.8; 95% CI, 1.3-2.5). After QC, 25 children with viral respiratory failure harbored excess ultrarare LOF variants in intolerant genes but only in those without a known disease association (odds ratio, 2.8; 95% CI, 1.1-6.6). A total of 114 undiagnosed children were enriched for de novo LOF variants in genes without a known disease association (observed, 14; expected, 6.8; enrichment, 2.05).

CONCLUSIONS AND RELEVANCE

In this genetic association study, excess LOF variants were observed among critically ill children despite nondiagnostic ES. Variants lay in genes without a known disease association, suggesting future investigation may connect phenotypes to causative genes.

Alzheimer’s Disease Polygenic Risk Score Is Not Associated With Cognitive Decline Among Older Adults With Type 2 Diabetes

Objectives

Multiple risk loci for late-onset Alzheimer’s disease (LOAD) have been identified. Type 2 diabetes (T2D) is a risk factor for cognitive decline, dementia and Alzheimer’s disease (AD). We investigated the association of polygenic risk score (PRS) for LOAD with overall cognitive functioning and longitudinal decline, among older adults with T2D.

Methods

The study included 1046 Jewish participants from the Israel Diabetes and Cognitive Decline (IDCD) study, aged ≥ 65 years, diagnosed with T2D, and cognitively normal at baseline. The PRS included variants from 26 LOAD associated loci (at genome-wide significance level), and was calculated with and without APOE. Outcome measures, assessed in 18 months intervals, were global cognition and the specific domains of episodic memory, attention/working memory, executive functions, and language/semantic categorization. Random coefficient models were used for analysis, adjusting for demographic variables, T2D-related characteristics, and cardiovascular factors. Additionally, in a subsample of 202 individuals, we analyzed the association of PRS with the volumes of total gray matter, frontal lobe, hippocampus, amygdala, and white matter hyperintensities. Last, the association of PRS with amyloid beta (Aβ) burden was examined in 44 participants who underwent an 18F-flutemetamol PET scan.

Results

The PRS was not significantly associated with overall functioning or decline in global cognition or any of the specific cognitive domains. Similarly, following correction for multiple testing, there was no association with Aβ burden and other brain imaging phenotypes.

Conclusion

Our results suggest that the cumulative effect of LOAD susceptibility loci is not associated with a greater rate of cognitive decline in older adults with T2D, and other pathways may underlie this link.

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.

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.

Detection of a mosaic CDKL5 deletion and inversion by optical genome mapping ends an exhaustive diagnostic odyssey

BACKGROUND

Currently available structural variant (SV) detection methods do not span the complete spectrum of disease-causing SVs. Optical genome mapping (OGM), an emerging technology with the potential to resolve diagnostic dilemmas, was performed to investigate clinically-relevant SVs in a 4-year-old male with an epileptic encephalopathy of undiagnosed molecular origin.

METHODS

OGM was utilized to image long, megabase-size DNA molecules, fluorescently labeled at specific sequence motifs throughout the genome with high sensitivity for detection of SVs greater than 500 bp in size. OGM results were confirmed in a CLIA-certified laboratory via mate-pair sequencing.

RESULTS

OGM identified a mosaic, de novo 90 kb deletion and inversion on the X chromosome disrupting the CDKL5 gene. Detection of the mosaic deletion, which had been previously undetected by chromosomal microarray, an infantile epilepsy panel including exon-level microarray for CDKL5, exome sequencing as well as genome sequencing, resulted in a diagnosis of X-linked dominant early infantile epileptic encephalopathy-2.

CONCLUSION

OGM affords an effective technology for the detection of SVs, especially those that are mosaic, since these remain difficult to detect with current NGS technologies and with conventional chromosomal microarrays. Further research in undiagnosed populations with OGM is warranted.

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.

Association of the CD2AP locus with cognitive functioning among middle-aged individuals with a family history of Alzheimer's disease

First-degree family history is an established risk factor for Alzheimer's disease (AD). We investigated the association of late-onset AD risk loci with cognitive functioning among 315 offspring of AD patients. Participants were cognitively normal Jewish individuals, aged 40-65 years, from the Israel Registry for Alzheimer's Prevention (IRAP) study. Twenty-two single-nucleotide polymorphisms (SNPs) within these loci and the APOE E4 allele were included in the final analyses, and a polygenic risk score was also calculated. Using linear regression (assuming an additive genetic model), we found a significant association only for SNP rs9473117, located near the CD2-associated protein (CD2AP) gene, with global cognition. Controlling for demographic variables (age, sex, years of education, and ancestry), the late-onset AD risk allele C was associated with lower global cognitive functioning (p = 0.0005), and withstood correction for multiple testing. After adjusting for additional characteristics (APOE E4 status and then also for cardiovascular factors), the results remained essentially unchanged (p = 0.0003 and p = 0.0005, respectively). In secondary analyses examining specific cognitive domains, rs9473117 was similarly associated with episodic memory (p = 0.005), language (p = 0.009), and working memory/attention (p = 0.018) but not with executive functions (p = 0.27). Again, the results were similar after adjusting for APOE E4 status and cardiovascular factors. The polygenic risk score was not associated with global cognitive functioning or with any of the 4 domains. In conclusion, our findings suggest a contribution of the CD2AP locus to cognitive functioning in middle-aged individuals with a parental history of AD. Further validations, including in longitudinal studies, are required.

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.

TCF7L2 polymorphisms are associated with amygdalar volume in elderly individuals with Type 2 Diabetes

The association between several Single Nucleotide Polymorphisms (SNPs) within the transcription factor 7-like 2 (TCF7L2) gene and Type 2 Diabetes (T2D) as well as additional T2D-related traits is well established. Since alteration in total and regional brain volumes are consistent findings among T2D individuals, we studied the association of four T2D susceptibility SNPS within TCF7L2 (rs7901695, rs7903146, rs11196205, and rs12255372) with volumes of white matter hyperintensities (WMH), gray matter, and regional volumes of amygdala and hippocampus obtained from structural MRI among 191 T2D elderly Jewish individuals. Under recessive genetic model (controlling for age, sex and intracranial volume), we found that for all four SNPs, carriers of two copies of the T2D risk allele (homozygous genotype) had significantly smaller amygdalar volume: rs7901695- CC genotype vs. CT + TT genotypes, p = 0.002; rs7903146-TT vs. TC + CC, p = 0.003; rs11196205- CC vs. CG + GG, p = 0.0003; and rs12255372- TT vs. TG + GG, p = 0.003. Adjusting also for T2D-related covariates, body mass index (BMI), and ancestry did not change the results substantively (rs7901695, p = 0.003; rs7903146, p = 0.005; rs11196205, p = 0.001; and rs12255372, p = 0.005). Conditional analysis demonstrated that only rs11196205 was independently associated with amygdalar volume at a significant level. Separate analysis of left and right amygdala revealed stronger results for left amygdalar volume. Taken together, we report association of TCF7L2 SNPs with amygdalar volume among T2D elderly Jewish patients. Further studies in other populations are required to support these findings and reach more definitive conclusions.

The CADM2 gene is associated with processing speed performance - evidence among elderly with type 2 diabetes

Objectives: Recent large-scale meta-analysis of genome-wide association studies (GWAS) from multiple cohorts, demonstrated the association of the single nucleotide polymorphism (SNP) rs17518584, with processing speed (measured by the Digit Symbol Substitution Test (DSST) or the Letter Digit Substitution Test (LDST)), at GWAS significance level. This SNP is located within the cell adhesion molecule 2 (CADM2) gene. We aimed to validate this finding in our sample of 944 cognitively normal Jewish elderly individuals with type 2 diabetes (T2D), a population which is at risk for cognitive decline and dementia.Methods: Using linear regression, we studied the association of rs17518584 with DSST performance, adjusting for demographic, T2D-related characteristics and cardiovascular factors. In secondary analyses, associations with performance in four cognitive domains (episodic memory, language/semantic categorisation, attention/working memory and executive function) and overall cognition were examined.Results: Controlling for sex, age at cognitive assessment, years of education and ancestry, we found a significant association of rs17518584 with DSST performance (P = 0.013), consistent with the originally reported effect direction. Results remained significant even when the additional covariates (T2D-related and cardiovascular factors) were included in the analysis (P = 0.034). Moreover, this SNP was significantly associated with performance in the cognitive domains of language/semantic categorisation and executive function, as well as overall cognition.Conclusions: Taken together, irrespective of T2D-related characteristics and cardiovascular factors, our findings provide independent support for the association of CADM2 SNP rs17518584 with processing speed (and demonstrate association with additional cognitive phenotypes), among cognitively normal elderly individuals with T2D.

New insights into tardive dyskinesia genetics: Implementation of whole-exome sequencing approach

Tardive dyskinesia (TD) is an adverse movement disorder induced by chronic treatment with antipsychotics drugs. The contribution of common genetic variants to TD susceptibility has been investigated in recent years, but with limited success. The aim of the current study was to investigate the potential contribution of rare variants to TD vulnerability. In order to identify TD risk genes, we performed whole-exome sequencing (WES) and gene-based collapsing analysis focusing on rare (allele frequency < 1%) and putatively deleterious variants (qualifying variants). 82 Jewish schizophrenia patients chronically treated with antipsychotics were included and classified as having severe TD or lack of any abnormal movements based on a rigorous definition of the TD phenotype. First, we performed a case-control, exome-wide collapsing analysis comparing 39 schizophrenia patients with severe TD to 3118 unrelated population controls. Then, we checked the potential top candidate genes among 43 patients without any TD manifestations. All the genes that were found to harbor one or more qualifying variants in patients without any TD features were excluded from the final list of candidate genes. Only one gene, regulating synaptic membrane exocytosis 2 (RIMS2), showed significant enrichment of qualifying variants in TD patients compared with unrelated population controls after correcting for multiple testing (Fisher's exact test p = 5.32E-08, logistic regression p = 2.50E-08). Enrichment was caused by a single variant (rs567070433) due to a frameshift in an alternative transcript of RIMS2. None of the TD negative patients had qualifying variants in this gene. In a validation cohort of 140 schizophrenia patients assessed for TD, the variant was also not detected in any individual. Some potentially suggestive TD genes were detected in the TD cohort and warrant follow-up in future studies. No significant enrichment in previously reported TD candidate genes was identified. To the best of our knowledge, this is the first WES study of TD, demonstrating the potential role of rare loss-of-function variant enrichment in this pharmacogenetic phenotype.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2722-2) contains supplementary material, which is available to authorized users.

A role for TENM1 mutations in congenital general anosmia

Congenital general anosmia (CGA) is a neurological disorder entailing a complete innate inability to sense odors. While the mechanisms underlying vertebrate olfaction have been studied in detail, there are still gaps in our understanding of the molecular genetic basis of innate olfactory disorders. Applying whole-exome sequencing to a family multiply affected with CGA, we identified three members with a rare X-linked missense mutation in the TENM1 (teneurin 1) gene (ENST00000422452:c.C4829T). In Drosophila melanogaster, TENM1 functions in synaptic-partner-matching between axons of olfactory sensory neurons and target projection neurons and is involved in synapse organization in the olfactory system. We used CRISPR-Cas9 system to generate a Tenm1 disrupted mouse model. Tenm1(-/-) and point-mutated Tenm1(A) (/A) adult mice were shown to have an altered ability to locate a buried food pellet. Tenm1(A) (/A) mice also displayed an altered ability to sense aversive odors. Results of our study, that describes a new Tenm1 mouse, agree with the hypothesis that TENM1 has a role in olfaction. However, additional studies should be done in larger CGA cohorts, to provide statistical evidence that loss-of-function mutations in TENM1 can solely cause the disease in our and other CGA cases.

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.

Common variants of IRF3 conferring risk of schizophrenia

Schizophrenia is a brain disorder with high heritability. Recent studies have implicated genes involved in the immune response pathway in the pathogenesis of schizophrenia. Interferon regulatory factor 3 (IRF3), a virus-immune-related gene, activates the transcription of several interferon-induced genes, and functionally interacts with several schizophrenia susceptibility genes. To test whether IRF3 is a schizophrenia susceptibility gene, we analyzed the associations of its SNPs with schizophrenia in independent population samples as well as reported data from expression quantitative trait loci (eQTL) in healthy individuals. We observed multiple independent SNPs in IRF3 showing nominally significant associations with schizophrenia (P < 0.05); more intriguingly, a SNP (rs11880923), which is significantly correlated with IRF3 expression in independent samples (P < 0.05), is also consistently associated with schizophrenia across different cohorts and in combined samples (odds ratio = 1.075, Pmeta = 2.08 × 10(-5)), especially in Caucasians (odds ratio = 1.078, Pmeta = 2.46 × 10(-5)). These results suggested that IRF3 is likely a risk gene for schizophrenia, at least in Caucasians. Although the clinical associations of IRF3 with diagnosis did not achieve genome-wide level of statistical significance, the observed odds ratio is comparable with other susceptibility loci identified through large-scale genetic association studies on schizophrenia, which could be regarded simply as small but detectable effects.

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.

Identification of new schizophrenia susceptibility loci in an ethnically homogeneous, family-based, Arab-Israeli sample

While the use of population-based samples is a common strategy in genome-wide association studies (GWASs), family-based samples have considerable advantages, such as robustness against population stratification and false-positive associations, better quality control, and the possibility to check for both linkage and association. In a genome-wide linkage study of schizophrenia in Arab-Israeli families with multiple affected individuals, we previously reported significant evidence for a susceptibility locus at chromosome 6q23.2-q24.1 and suggestive evidence at chromosomes 10q22.3-26.3, 2q36.1-37.3 and 7p21.1-22.3. To identify schizophrenia susceptibility genes, we applied a family-based GWAS strategy in an enlarged, ethnically homogeneous, Arab-Israeli family sample. We performed genome-wide single nucleotide polymorphism (SNP) genotyping and single SNP transmission disequilibrium test association analysis and found genome-wide significant association (best value of P=1.22×10(-11)) for 8 SNPs within or near highly reasonable functional candidate genes for schizophrenia. Of particular interest are a group of SNPs within and flanking the transcriptional factor LRRFIP1 gene. To determine replicability of the significant associations beyond the Arab-Israeli population, we studied the association of the significant SNPs in a German case-control validation sample and found replication of associations near the UGT1 subfamily and EFHD1 genes. Applying an exploratory homozygosity mapping approach as a complementary strategy to identify schizophrenia susceptibility genes in our Arab Israeli sample, we identified 8 putative disease loci. Overall, this GWAS, which emphasizes the important contribution of family based studies, identifies promising candidate genes for schizophrenia.

Why do young women smoke? VII COMT as a risk modifying gene for Nicotine dependence - role of gene-gene interaction, personality, and environmental factors

OBJECTIVES

Catechol-O-methyltransferase (COMT) may be a risk modifying gene for Nicotine dependence (ND) rather than a direct susceptibility gene for this phenotype. Brain nicotinic cholinergic receptors modulate dopaminergic transmission, and several variants within the neighboring CHRNA5-CHRNA3 genes have been associated with ND. Therefore, it is biologically reasonable to study the interactive contribution of COMT and the CHRNA5 and CHRNA3 genes to ND.

METHODS

Using a case-control sample of 90 young, Israeli, Jewish female smokers (FTND ≥ 4) and 108 controls (FTND = 0 during heaviest period of smoking), we studied association with ND of 8 COMT tagging SNPs, their interaction with tagging CHRNA5-A3 SNPs and the role of background, personality, and environmental factors.

RESULTS

None of the COMT SNPs were associated directly with ND. In pairwise interaction analysis of SNPs from the two loci (COMT SNP-CHRNA5-CHRNA3 SNP), the interaction of intronic COMT SNP, rs9332377, with CHRNA3 3'UTR SNP rs660652 was significantly associated with ND (p = 0.0005), withstanding correction for multiple testing.

CONCLUSION

Addition of the genetic interaction variable into a model of non-genetic ND predictors [parental smoking, novelty seeking (NS), and lifetime history of trauma], substantially increases the percentage of ND variance explained by the model, as well as the percentage of cases correctly identified by it.

Evidence for an interaction of schizophrenia susceptibility loci on chromosome 6q23.3 and 10q24.33-q26.13 in Arab Israeli families

A genome scan for schizophrenia related loci in Arab Israeli families by Lerer et al. [Lerer et al. (2003); Mol Psychiatry 8:488-498] detected significant evidence for linkage at chromosome 6q23. Subsequent fine mapping [Levi et al. (2005); Eur J Hum Genet 13:763-771], association [Amann-Zalcenstein et al. (2006); Eur J Hum Genet 14:1111-1119] and replication studies [Ingason et al. (2007); Eur J Hum Genet 15:988-991] identified AHI1 as a putative susceptibility gene. The same genome scan revealed suggestive evidence for a schizophrenia susceptibility locus in the 10q23-26 region. Genes at these two loci may act independently in the pathogenesis of the disease in our homogeneous sample of Arab Israeli families or may interact with each other and with other factors in a common biological pathway. The purpose of our current study was to test the hypothesis of genetic interaction between these two loci and to identify the type of interaction between them. The initial stage of our study focused on the 10q23-q26 region which has not been explored further in our sample. The second stage of the study included a test for possible genetic interaction between the 6q23.3 locus and the refined 10q24.33-q26.13 locus. A final candidate region of 19.9 Mb between markers D10S222 (105.3 Mb) and D10S587 (125.2 Mb) was found on chromosome 10 by non-parametric and parametric linkage analyses. These linkage findings are consistent with previous reports in the same chromosomal region. Two-locus multipoint linkage analysis under three complex disease inheritance models (heterogeneity, multiplicative, and additive models) yielded a best maximum LOD score of 7.45 under the multiplicative model suggesting overlapping function of the 6q23.3 and 10q24.33-q26.13 loci.

The XVth World Congress of Psychiatric Genetics, October 7-11, 2007: Rapporteur summaries of oral presentations

The World Congress of Psychiatric Genetics (WCPG) has become an annual event since the early 1990's sponsored by the International Society of Psychiatric Genetics (ISPG). Each year the latest published and unpublished findings are aired for discussion by representatives of the majority of research programs on this topic world-wide. The 2007 congress was held in New York City and attracted over 1000 researchers. The topics emphasized included results from whole genome association studies, the significance of copy number variation and the important contributions of epigenetic events to psychiatric disorders. There were over 20 oral sessions devoted to these and other topics of interest. Young investigator recipients of travel awards served as rapporteurs to summarize sessions and these summaries follow.

The CC2D1A, a member of a new gene family with C2 domains, is involved in autosomal recessive non-syndromic mental retardation

BACKGROUND

The molecular basis of autosomal recessive non-syndromic mental retardation (NSMR) is poorly understood, mostly owing to heterogeneity and absence of clinical criteria for grouping families for linkage analysis. Only two autosomal genes, the PRSS12 gene on chromosome 4q26 and the CRBN on chromosome 3p26, have been shown to cause autosomal recessive NSMR, each gene in only one family.

OBJECTIVE

To identify the gene causing autosomal recessive NSMR on chromosome 19p13.12.

RESULTS

The candidate region established by homozygosity mapping was narrowed down from 2.4 Mb to 0.9 Mb on chromosome 19p13.12. A protein truncating mutation was identified in the gene CC2D1A in nine consanguineous families with severe autosomal recessive NSMR. The absence of the wild type protein in the lymphoblastoid cells of the patients was confirmed. CC2D1A is a member of a previously uncharacterised gene family that carries two conserved motifs, a C2 domain and a DM14 domain. The C2 domain is found in proteins which function in calcium dependent phospholipid binding; the DM14 domain is unique to the CC2D1A protein family and its role is unknown. CC2D1A is a putative signal transducer participating in positive regulation of I-kappaB kinase/NFkappaB cascade. Expression of CC2D1A mRNA was shown in the embryonic ventricular zone and developing cortical plate in staged mouse embryos, persisting into adulthood, with highest expression in the cerebral cortex and hippocampus.

CONCLUSIONS

A previously unknown signal transduction pathway is important in human cognitive development.

Mapping of a new locus for autosomal recessive non-syndromic mental retardation in the chromosomal region 19p13.12-p13.2: further genetic heterogeneity

OBJECTIVE

To identify and clinically evaluate four consanguineous families of Israeli Arab origin with non-syndromic mental retardation (NSMR), comprising a total of 10 affected and 24 unaffected individuals.

PARTICIPANTS AND METHODS

All the families originated from the same small village and had the same family name. Association of the condition in these families with the two known autosomal recessive NSMR loci on chromosomes 3p25-pter and 4q24 (neurotrypsin gene) was excluded.

RESULTS

Linkage of the disease gene to chromosome 19p13.12-p13.2(Zmax = 7.06 at theta = 0.00) for the marker D19S840 was established. All the affected individuals were found to be homozygous for a common haplotype for the markers cen-RFX1-D19S840-D19S558-D19S221-tel.

CONCLUSIONS

The results suggest that the disease is caused by a single mutation derived from a single ancestral founder in all the families. Recombination events and a common disease bearing haplotype defined a critical region of 2.4 Mb, between the loci D19S547 proximally and D19S1165 distally.