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.

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.

General olfactory sensitivity database (GOSdb): candidate genes and their genomic variations.

Genetic variations in olfactory receptors likely contribute to the diversity of odorant-specific sensitivity phenotypes. Our working hypothesis is that genetic variations in auxiliary olfactory genes, including those mediating transduction and sensory neuronal development, may constitute the genetic basis for general olfactory sensitivity (GOS) and congenital general anosmia (CGA). We thus performed a systematic exploration for auxiliary olfactory genes and their documented variation. This included a literature survey, seeking relevant functional in vitro studies, mouse gene knockouts and human disorders with olfactory phenotypes, as well as data mining in published transcriptome and proteome data for genes expressed in olfactory tissues. In addition, we performed next-generation transcriptome sequencing (RNA-seq) of human olfactory epithelium and mouse olfactory epithelium and bulb, so as to identify sensory-enriched transcripts. Employing a global score system based on attributes of the 11 data sources utilized, we identified a list of 1,680 candidate auxiliary olfactory genes, of which 450 are shortlisted as having higher probability of a functional role. For the top-scoring 136 genes, we identified genomic variants (probably damaging single nucleotide polymorphisms, indels, and copy number deletions) gleaned from public variation repositories. This database of genes and their variants should assist in rationalizing the great interindividual variation in human overall olfactory sensitivity (http://genome.weizmann.ac.il/GOSdb).

Personal receptor repertoires: olfaction as a model.

BACKGROUND: Information on nucleotide diversity along completely sequenced human genomes has increased tremendously over the last few years. This makes it possible to reassess the diversity status of distinct receptor proteins in different human individuals. To this end, we focused on the complete inventory of human olfactory receptor coding regions as a model for personal receptor repertoires. RESULTS: By performing data-mining from public and private sources we scored genetic variations in 413 intact OR loci, for which one or more individuals had an intact open reading frame. Using 1000 Genomes Project haplotypes, we identified a total of 4069 full-length polypeptide variants encoded by these OR loci, average of ~10 per locus, constituting a lower limit for the effective human OR repertoire. Each individual is found to harbor as many as 600 OR allelic variants, ~50% higher than the locus count. Because OR neuronal expression is allelically excluded, this has direct effect on smell perception diversity of the species. We further identified 244 OR segregating pseudogenes (SPGs), loci showing both intact and pseudogene forms in the population, twenty-six of which are annotatively "resurrected" from a pseudogene status in the reference genome. Using a custom SNP microarray we validated 150 SPGs in a cohort of 468 individuals, with every individual genome averaging 36 disrupted sequence variations, 15 in homozygote form. Finally, we generated a multi-source compendium of 63 OR loci harboring deletion Copy Number Variations (CNVs). Our combined data suggest that 271 of the 413 intact OR loci (66%) are affected by nonfunctional SNPs/indels and/or CNVs. CONCLUSIONS: These results portray a case of unusually high genetic diversity, and suggest that individual humans have a highly personalized inventory of functional olfactory receptors, a conclusion that might apply to other receptor multigene families.

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.

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.

Fine mapping of AHI1 as a schizophrenia susceptibility gene: from association to evolutionary evidence.

In previous studies, we identified a locus for schizophrenia on 6q23.3 and proposed the Abelson helper integration site 1 (AHI1) as the candidate gene. AHI1 is expressed in the brain and plays a key role in neurodevelopment, is involved in Joubert syndrome, and has been recently associated with autism. The neurodevelopmental role of AHI1 fits with etiological hypotheses of schizophrenia. To definitively confirm our hypothesis, we searched for associations using a dense map of the region. Our strongest findings lay within the AHI1 gene: single-nucleotide polymorphisms rs11154801 and rs7759971 showed significant associations (P=6.23E-06; P=0.84E-06) and haplotypes gave P values in the 10E-8 to 10E-10 range. The second highest significant region maps close to AHI1 and includes the intergenic region between BC040979 and PDE7B (rs2038549 at P=9.70E-06 and rs1475069 at P=6.97E-06), and PDE7B and MAP7. Using a sample of Palestinian Arab families to confirm these findings, we found isolated signals. While these results did not retain their significance after correction for multiple testing, the joint analysis across the 2 samples supports the role of AHI1, despite the presence of heterogeneity. Given the hypothesis of positive selection of schizophrenia genes, we resequenced a 11 kb region within AHI1 in ethnically defined populations and found evidence for a selective sweep. Network analysis indicates 2 haplotype clades, with schizophrenia-susceptibility haplotypes clustering within the major clade. In conclusion, our data support the role of AHI1 as a susceptibility gene for schizophrenia and confirm it has been subjected to positive selection, also shedding light on new possible candidate genes, MAP7 and PDE7B.

Single-nucleotide polymorphisms in the p53 pathway genes modify cancer risk in BRCA1 and BRCA2 carriers of Jewish-Ashkenazi descent.

Germline mutations in the BRCA1 and BRCA2 genes are associated with a significantly increased lifetime risk for developing breast and/or ovarian cancer. However, incomplete penetrance and substantial variability in age of disease onset among carriers of the same mutation suggests the involvement of additional modifier genes and/or environmental factors. Somatic inactivating mutations in the p53 gene and genes of the p53 pathway often accompany BRCA1/2-associated tumors. Therefore, we assessed whether these genes are modifiers of penetrance. We genotyped Jewish-Ashkenazi women for functional single-nucleotide polymorphisms (SNPs) in the AKT1 (C>T rs3730358) and the PERP (C>T rs2484067) genes that affect p53-mediated apoptosis, as well as two tag-SNPs in the CHEK2 (C>T rs743184) and the ZBRK1/ZNF350 (G>A rs2278414) genes that encode for proteins involved in growth arrest following DNA damage. The study population included 138 healthy women, 148 breast/ovarian cancer BRCA1/2 mutation carriers, 121 asymptomatic BRCA1/2 mutation carriers, and 210 sporadic noncarrier breast cancer patients. Utilizing lambda(2) and Kaplan-Meier analysis revealed a hazard ratio (HR) of 3.23 (95% CI: 1.44-54, P = 0.0184) for the TT genotype of AKT (rs3730358), HR = 2.105 (95% CI: 1.049-7.434, P = 0.039) for CHEK2 CC genotype (rs743184), and HR = 2.4743 (95% CI: 1.205-11.53, P = 0.022) for the AG genotype of ZBRK1/ZNF350 (rs2278414). No significant association between PERP variants and cancer was identified HR = 0.662 (95% CI: 0.289-1.324, P = 0.261). Our results suggest that genes that act upstream of p53, or participate in the DNA damage response, may modify the risk of cancer in women with mutant BRCA1/2 alleles.

Genome-wide linkage scan, fine mapping, and haplotype analysis in a large, inbred, Arab Israeli pedigree suggest a schizophrenia susceptibility locus on chromosome 20p13.

Linkage and association studies in schizophrenia have repeatedly drawn attention to several chromosomal regions and to genes within them. Conflicting patterns of association and the lack of a clear functional significance of the associated variants limit the interpretation of these results. The use of rare pedigrees, where genes with a major effect cause the disorder, has been proven beneficial in studies of other complex disorders. Our objective was to use this advantage by performing a genome wide linkage analysis for schizophrenia in a large, multiplex Israeli Arab pedigree. We genotyped 346 microsatellite markers in 24 pedigree members affected with schizophrenia spectrum disorders and 32 unaffected relatives. Two-point linkage analysis with SUPERLINK demonstrated a LOD score of 2.47 for D20S116 on chromosome 20p13 under an autosomal dominant mode of inheritance. Further fine mapping yielded a two-point LOD score of 2.56 for the adjacent marker D20S193 and narrowed down the linked region to 2-5 cM. A haplotype containing the markers D20S193, D20S889, and D20S116, 0.7 Mb in length, was found to be shared by most affected pedigree members. Genotyping of 43 SNPs in the interval supported these results with a multipoint LOD score of 2.7 around D20S193. We were also able to better define the boundaries of the shared haplotype which contains strong candidate genes for schizophrenia. Our study exemplifies the power of rare and unique pedigrees in drawing attention to novel regions for genetic studies of schizophrenia.

AHI1, a pivotal neurodevelopmental gene, and C6orf217 are associated with susceptibility to schizophrenia.

Schizophrenia, a severe neuropsychiatric disorder, is believed to involve multiple genetic factors. A significant body of evidence supports a pivotal role for abnormalities of brain development in the disorder. Linkage signals for schizophrenia map to human chromosome 6q. To obtain a finer localization, we genotyped 180 single nucleotide polymorphisms (SNPs) in a young, inbred Arab-Israeli family sample with a limited number of founders. The SNPs were mostly within a approximately 7 Mb region around the strong linkage peak at 136.2 Mb that we had previously mapped. The most significant genetic association with schizophrenia for single SNPs and haplotypes was within a 500 kb genomic region of high linkage disequilibrium (LD) at 135.85 Mb. In a different, outbred, nuclear family sample that was not appropriate for linkage analysis, under-transmitted haplotypes incorporating the same SNPs (but not the individual SNPs) were significantly associated with schizophrenia. The implicated genomic region harbors the Abelson Helper Integration Site 1 (AHI1) gene, which showed the strongest association signal, and an adjacent, primate-specific gene, C6orf217. Mutations in human AHI1 underlie the autosomal recessive Joubert Syndrome with brain malformation and mental retardation. Previous comparative genomic analysis has suggested accelerated evolution of AHI1 in the human lineage. C6orf217 has multiple splice isoforms and is expressed in brain but does not seem to encode a functional protein. The two genes appear in opposite orientations and their regulatory upstream regions overlap, which might affect their expression. Both, AHI1 and C6orf217 appear to be highly relevant candidate genes for schizophrenia.

Is the G72/G30 locus associated with schizophrenia? single nucleotide polymorphisms, haplotypes, and gene expression analysis.

BACKGROUND: The genes G72/G30 were recently implicated in schizophrenia in both Canadian and Russian populations. We hypothesized that 1) polymorphic changes in this gene region might be associated with schizophrenia in the Ashkenazi Jewish population and that 2) changes in G72/G30 gene expression might be expected in schizophrenic patients compared with control subjects. METHODS: Eleven single nucleotide polymorphisms (SNPs) encompassing the G72/G30 genes were typed in the genomic deoxyribonucleic acid (DNA) from 60 schizophrenic patients and 130 matched control subjects of Ashkenazi ethnic origin. Case-control comparisons were based on linkage disequilibrium (LD) and haplotype frequency estimations. Gene expression analysis of G72 and G30 was performed on 88 postmortem dorsolateral prefrontal cortex samples. RESULTS: Linkage disequilibrium analysis revealed two main SNP blocks. Haplotype analysis on block II, containing three SNPs external to the genes, demonstrated an association with schizophrenia. Gene expression analysis exhibited correlations between expression levels of the G72 and G30 genes, as well as a tendency toward overexpression of the G72 gene in schizophrenic brain samples of 44 schizophrenic patients compared with 44 control subjects. CONCLUSIONS: It is likely that the G72/G30 region is involved in susceptibility to schizophrenia in the Ashkenazi population. The elevation in expression of the G72 gene coincides with the glutamatergic theory of schizophrenia.

CATSPER2, a human autosomal nonsyndromic male infertility gene.

In the course of positional cloning of the Congenital Dyserythropoietic Anemia type I (CDAI) [MIM 224120] gene on 15q15.1-15.3, we examined a family of French origin, in which the propositus suffered from asthenoteratozoospermia and nonsyndromic deafness in addition to CDAI. Two of his brothers had a similar phenotype. All three siblings were homozygous carriers of the CDA1 mutation as well as of a distally located approximately 70 kb deletion of the proximal copy of a 106 kb tandem repeat on chromosome 15q15. These repeats encode four genes whose distal copies may be considered pseudogenes. Lack of functional stereocilin and CATSPER2 (a voltage-gate cation channel expressed specifically in spermatozoa) may explain the observed deafness and male infertility phenotypes. To the best of our knowledge, the involvement of CATSPER2 in asthenoteratozoospermia is the first description of a human autosomal gene defect associated with nonsyndromic male infertility.

USH3A transcripts encode clarin-1, a four-transmembrane-domain protein with a possible role in sensory synapses.

Usher syndrome type 3 (USH3) is an autosomal recessive disorder characterised by the association of post-lingual progressive hearing loss, progressive visual loss due to retinitis pigmentosa and variable presence of vestibular dysfunction. Because the previously defined transcripts do not account for all USH3 cases, we performed further analysis and revealed the presence of additional exons embedded in longer human and mouse USH3A transcripts and three novel USH3A mutations. Expression of Ush3a transcripts was localised by whole mount in situ hybridisation to cochlear hair cells and spiral ganglion cells. The full length USH3A transcript encodes clarin-1, a four-transmembrane-domain protein, which defines a novel vertebrate-specific family of three paralogues. Limited sequence homology to stargazin, a cerebellar synapse four-transmembrane-domain protein, suggests a role for clarin-1 in hair cell and photoreceptor cell synapses, as well as a common pathophysiological pathway for different Usher syndromes.

Deficiency of asparagine synthetase causes congenital microcephaly and a progressive form of encephalopathy.

We analyzed four families that presented with a similar condition characterized by congenital microcephaly, intellectual disability, progressive cerebral atrophy, and intractable seizures. We show that recessive mutations in the ASNS gene are responsible for this syndrome. Two of the identified missense mutations dramatically reduce ASNS protein abundance, suggesting that the mutations cause loss of function. Hypomorphic Asns mutant mice have structural brain abnormalities, including enlarged ventricles and reduced cortical thickness, and show deficits in learning and memory mimicking aspects of the patient phenotype. ASNS encodes asparagine synthetase, which catalyzes the synthesis of asparagine from glutamine and aspartate. The neurological impairment resulting from ASNS deficiency may be explained by asparagine depletion in the brain or by accumulation of aspartate/glutamate leading to enhanced excitability and neuronal damage. Our study thus indicates that asparagine synthesis is essential for the development and function of the brain but not for that of other organs.

Association of the type 2 diabetes mellitus susceptibility gene, TCF7L2, with schizophrenia in an Arab-Israeli family sample.

Many reports in different populations have demonstrated linkage of the 10q24-q26 region to schizophrenia, thus encouraging further analysis of this locus for detection of specific schizophrenia genes. Our group previously reported linkage of the 10q24-q26 region to schizophrenia in a unique, homogeneous sample of Arab-Israeli families with multiple schizophrenia-affected individuals, under a dominant model of inheritance. To further explore this candidate region and identify specific susceptibility variants within it, we performed re-analysis of the 10q24-26 genotype data, taken from our previous genome-wide association study (GWAS) (Alkelai et al, 2011). We analyzed 2089 SNPs in an extended sample of 57 Arab Israeli families (189 genotyped individuals), under the dominant model of inheritance, which best fits this locus according to previously performed MOD score analysis. We found significant association with schizophrenia of the TCF7L2 gene intronic SNP, rs12573128, (p = 7.01×10⁻6) and of the nearby intergenic SNP, rs1033772, (p = 6.59×10⁻6) which is positioned between TCF7L2 and HABP2. TCF7L2 is one of the best confirmed susceptibility genes for type 2 diabetes (T2D) among different ethnic groups, has a role in pancreatic beta cell function and may contribute to the comorbidity of schizophrenia and T2D. These preliminary results independently support previous findings regarding a possible role of TCF7L2 in susceptibility to schizophrenia, and strengthen the importance of integrating linkage analysis models of inheritance while performing association analyses in regions of interest. Further validation studies in additional populations are required.

Lymphoblast and brain expression of AHI1 and the novel primate-specific gene, C6orf217, in schizophrenia and bipolar disorder.

Association with schizophrenia of the Abelson Helper Integration Site 1 (AHI1) gene on chromosome 6q23 and the adjacent primate-specific gene, C6orf217, was demonstrated in an inbred, Arab Israeli family sample and replicated in an Icelandic case control sample. Further support was provided by a second replication in a large European sample and a meta-analysis that supported association with schizophrenia of all seven alleles overtransmitted to affected subjects in the original study. We examined constitutive expression of AHI1 and C6orf217 in immortalized lymphoblasts of patients from the Arab Israeli family sample in which the association with schizophrenia was originally discovered and population-matched normal controls, and in post-mortem brain of patients with schizophrenia and bipolar (BP) disorder and control subjects from the Stanley Medical Research Institute Collection. We found a significant effect of diagnostic group in the lymphoblast sample (F=5.72; df=2,39; p=0.006). Patients with early age of onset had higher AHI1 expression than controls and later onset patients (p=0.002; 0.03 respectively). C6orf217 expression in lymphoblasts was too low to measure. We found no difference in brain expression of AHI1 in schizophrenia or BP patients compared to controls. However, there was a genotypic difference in AHI1 expression for SNP rs9321501, which was strongly associated with schizophrenia in the original study. Genotypes that included the undertransmitted C allele (CC/AC) showed lower expression than the homozygous AA genotype (F=4.73, df=2,83; p=0.028). There was no significant difference in brain expression of C6orf217 between patients and controls and no genotypic effect. This study provides further evidence for involvement of AHI1 in susceptibility to schizophrenia.