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.

Further evidence for association of the RGS2 gene with antipsychotic-induced parkinsonism: protective role of a functional polymorphism in the 3'-untranslated region.

RGS2 (regulator of G-protein signaling 2) modulates dopamine receptor signal transduction. Functional variants in the gene may influence susceptibility to extrapyramidal symptoms (EPS) induced by antipsychotic drugs. To further investigate our previous report of association of the RGS2 gene with susceptibility to antipsychotic-induced EPS, we performed a replication study. EPS were rated in 184 US patients with schizophrenia (115 African Americans, 69 Caucasian) treated for at least a month with typical antipsychotic drugs (n=45), risperidone (n=46), olanzapine (n=50) or clozapine (n=43). Six single nucleotide polymorphisms (SNPs) within or flanking RGS2 were genotyped (rs1933695, rs2179652, rs2746073, rs4606, rs1819741 and rs1152746). Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated by logistic regression. Our results indicate association of SNP rs4606 with antipsychotic-induced parkinsonism (AIP), as measured by the Simpson Angus scale, in the overall sample and in the African-American subsample, the G (minor) allele having a protective effect. ORs for AIP among rs4606 G-allele carriers were 0.23 (95% CI 0.10-0.54, P=0.001) in the overall sample, and 0.20 (0.07-0.57, P=0.003) in the African-American subsample. In the previously studied Israeli sample the OR was 0.31 (0.11-0.84, P=0.02). We completely sequenced the RGS2 gene in nine patients with AIP and nine patients without, from the Israeli sample. No common coding polymorphisms or additional regulatory variants were revealed, suggesting that association of the rs4606 C/G polymorphism with AIP is biologically meaningful and not a consequence of linkage disequilibrium with another functional variant. Taken together, the findings of the current study support the association of RGS2 with AIP and focus on a possible protective effect of the minor G allele of SNP rs4606. This SNP is located in the 3'-regulatory region of the gene, and is known to influence RGS2 mRNA levels and protein expression.

Why do young women smoke? V. Role of direct and interactive effects of nicotinic cholinergic receptor gene variation on neurocognitive function.

Previous work suggests that young women who smoke cigarettes regularly, or did so in the past, manifest a neurocognitive profile that is characterized by small but significant impairments of response inhibition and attention. The present study sought to determine whether variation in nicotinic cholinergic receptor (nAchR) genes impacts upon cognitive function in these domains by overall or differential effects on the performance of current, former and non-smokers. The study sample consisted of 100 female college students, current or past smokers, and 144 who had never smoked. All performed a computerized neurocognitive test battery and were genotyped for 39 single nucleotide polymorphisms in 11 nAchR genes. The results, derived from linear or logistic regression, show significant direct and interactive relationships between single nucleotide polymorphisms and haplotypes in several nAchR genes and performance on the Matching Familiar Figures Test (MFFT) Stroop test, Continuous Performance Task (CPT) and Tower of London (TOL) test. Response inhibition (MFFT, Stroop, CPT Loading Phase, TOL) was associated with variants in CHRNA2, CHRNA4, CHRNA5, CHRNA7, CHRNA9, CHRNA10, CHRNB2 and CHRNB3. Selective attention (Stroop) was associated with CHRNA4, CHRNA5, CHRNA9 and CHRNB2. Sustained attention (CPT Boring Phase) was associated with CHRNA4, CHRNA5, CHRNA7, CHRNA10 and CHRNB3. Up to 37% of the variance among the smokers and up to 47% of the variance among the non-smokers on the test measures was explained. Differences between smokers and non-smokers in neurocognitive function, putatively implicated in susceptibility to nicotine dependence, may be modulated by variants in nAchR genes, with potential implications for prevention and treatment.

Search for hand osteoarthritis susceptibility locus on chromosome 6p12.3-p12.1.

The existence of osteoarthritis susceptibility loci on chromosome 6 for individuals suffering from hip and knee osteoarthritis has been suggested. We determined whether radiographic hand osteoarthritis in a demographically homogeneous population of European origin can be linked to loci on chromosome 6p12.3-p12.1. Nine single nucleotide polymorphisms (SNPs) were genotyped in 764 individuals (members of 189 nuclear and more complex two- or three-generation families). Radiographic hand osteoarthritis was characterized by two traits: (1) the total individual osteoarthritis score (PC1-OA) and (2) the osteophytes score (PC1-OS), obtained from the principal components analysis of sums of the Kellgren and Lawrence grade and of the osteophyte grades, respectively, for 14 joints on each hand. The contribution of genetic and environmental factors and of covariates such as age and body mass index to hand osteoarthritis was evaluated by variance components analysis. The association between the studied traits and selected DNA markers was evaluated by three types of transmission disequilibrium tests. The parent-offspring and sib-sib correlations were statistically significant for all studied traits. The additive genetic effects for PC1-OA and PC1-OS were estimated to be 43% and 37.9%, respectively. Transmission disequilibrium tests consistently revealed a statistically significant association (p values ranged from 0.017 to 0.030) between SNP rs1508632 and PC1-OS. In the tested cohort the putative genetic factors are influential enough to determine interindividual differences regarding the extent of hand osteoarthritis. SNP rs1508632 lies in immediate proximity to the TINAG gene, implicating it as a possible hand osteoarthritis susceptibility gene.

ATM haplotypes and breast cancer risk in Jewish high-risk women.

While genetic factors clearly play a role in conferring breast cancer risk, the contribution of ATM gene mutations to breast cancer is still unsettled. To shed light on this issue, ATM haplotypes were constructed using eight SNPs spanning the ATM gene region (142 kb) in ethnically diverse non-Ashkenazi Jewish controls (n=118) and high-risk (n=142) women. Of the 28 haplotypes noted, four were encountered in frequencies of 5% or more and accounted for 85% of all haplotypes. Subsequently, ATM haplotyping of high-risk, non-Ashkenazi Jews was performed on 66 women with breast cancer and 76 asymptomatic. One SNP (rs228589) was significantly more prevalent among breast cancer cases compared with controls (P=4 x 10(-9)), and one discriminative ATM haplotype was significantly more prevalent among breast cancer cases (33.3%) compared with controls (3.8%), (P< or =10(-10)). There was no significant difference in the SNP and haplotype distribution between asymptomatic high-risk and symptomatic women as a function of disease status. We conclude that a specific ATM SNP and a specific haplotype are associated with increased breast cancer risk in high-risk non-Ashkenazi Jews.

Why do young women smoke? I. Direct and interactive effects of environment, psychological characteristics and nicotinic cholinergic receptor genes.

Despite the health hazards, cigarette smoking is disproportionately frequent among young women. A significant contribution of genetic factors to smoking phenotypes is well established. Efforts to identify susceptibility genes do not generally take into account possible interaction with environment, life experience and psychological characteristics. We recruited 501 female Israeli students aged 20-30 years, obtained comprehensive background data and details of cigarette smoking and administered a battery of psychological instruments. Smoking initiators (n=242) were divided into subgroups with high (n=127) and low (n=115) levels of nicotine dependence based on their scores on the Fagerstrom Tolerance Questionnaire and genotyped with noninitiators (n=142) for single nucleotide polymorphisms (SNPs) in 11 nicotinic cholinergic receptor genes. We found nominally significant (P<0.05) allelic and genotypic association with smoking initiation of SNP rs2072660 and multilocus haplotypes (P<0.007-0.05) in CHRNB2 and nominal (P<0.05) allelic or genotypic association of SNPs in CHRNA7 (rs1909884), CHRNA9 (rs4861065) and CHRNB3 (rs9298629) with nicotine dependence. Employing logistic regression and controlling for known risk factors, the best-fitting model for smoking initiation encompassed a 5 SNP haplotype in CHRNB2, neuroticism and novelty seeking (P=5.9 x 10(-14), Nagelkerke r(2)=0.30). For severity of nicotine dependence, two SNPs in CHRNA7 (rs1909884 and rs883473), one SNP in CHRNA5 (rs680244) and the interaction of a SNP in CHRNA7 (rs2337980) with neuroticism, were included in the model (P=2.24 x 10(-7), Nagelkerke r(2)=0.40). These findings indicate that background factors, psychological characteristics and genetic variation in nicotinic cholinergic receptors contribute independently or interactively to smoking initiation and to severity of nicotine dependence in young women.

The RUNX3 gene--sequence, structure and regulated expression.

The RUNX3 gene belongs to the runt domain family of transcription factors that act as master regulators of gene expression in major developmental pathways. In mammals the family includes three genes, RUNX1, RUNX2 and RUNX3. Here, we describe a comparative analysis of the human chromosome 1p36.1 encoded RUNX3 and mouse chromosome 4 encoded Runx3 genomic regions. The analysis revealed high similarities between the two genes in the overall size and organization and showed that RUNX3/Runx3 is the smallest in the family, but nevertheless exhibits all the structural elements characterizing the RUNX family. It also revealed that RUNX3/Runx3 bears a high content of the ancient mammalian repeat MIR. Together, these data delineate RUNX3/Runx3 as the evolutionary founder of the mammalian RUNX family. Detailed sequence analysis placed the two genes at a GC-rich H3 isochore with a sharp transition of GC content between the gene sequence and the downstream intergenic region. Two large conserved CpG islands were found within both genes, one around exon 2 and the other at the beginning of exon 6. RUNX1, RUNX2 and RUNX3 gene products bind to the same DNA motif, hence their temporal and spatial expression during development should be tightly regulated. Structure/function analysis showed that two promoter regions, designated P1 and P2, regulate RUNX3 expression in a cell type-specific manner. Transfection experiments demonstrated that both promoters were highly active in the GM1500 B-cell line, which endogenously expresses RUNX3, but were inactive in the K562 myeloid cell line, which does not express RUNX3.

A missense mutation in a highly conserved region of CASQ2 is associated with autosomal recessive catecholamine-induced polymorphic ventricular tachycardia in Bedouin families from Israel.

Catecholamine-induced polymorphic ventricular tachycardia (PVT) is characterized by episodes of syncope, seizures, or sudden death, in response to physical activity or emotional stress. Two modes of inheritance have been described: autosomal dominant and autosomal recessive. Mutations in the ryanodine receptor 2 gene (RYR2), which encodes a cardiac sarcoplasmic reticulum (SR) Ca(2+)-release channel, were recently shown to cause the autosomal dominant form of the disease. In the present report, we describe a missense mutation in a highly conserved region of the calsequestrin 2 gene (CASQ2) as the potential cause of the autosomal recessive form. The CASQ2 protein serves as the major Ca(2+) reservoir within the SR of cardiac myocytes and is part of a protein complex that contains the ryanodine receptor. The mutation, which is in full segregation in seven Bedouin families affected by the disorder, converts a negatively charged aspartic acid into a positively charged histidine, in a highly negatively charged domain, and is likely to exert its deleterious effect by disrupting Ca(2+) binding.

Mucolipidosis type IV: novel MCOLN1 mutations in Jewish and non-Jewish patients and the frequency of the disease in the Ashkenazi Jewish population.

The gene MCOLN1 is mutated in Mucolipidosis type IV (MLIV), a neurodegenerative, recessive, lysosomal storage disorder. The disease is found in relatively high frequency among Ashkenazi Jews due to two founder mutations that comprise 95% of the MLIV alleles in this population [Bargal et al., 2000]. In this report we complete the mutation analysis of Jewish and non-Jewish MLIV patients whose DNA were available to us. Four novel mutations were identified in the MCOLN1 gene of severely affected patients: two missense, T232P and F465L; a nonsense, R322X; and an 11-bp insertion in exon 12. The nonsense mutation (R322X) was identified in two unrelated patients with different haplotypes in the MCOLN1 chromosomal region, indicating a mutation hotspot in this CpG site. An in-frame deletion (F408del) was identified in a patient with unusual mild psychomotor retardation. The frequency of MLIV in the general Jewish Ashkenazi population was estimated in a sample of 2,000 anonymous, unrelated individuals assayed for the two founder mutations. This analysis indicated a heterozygotes frequency of about 1/100. A preferred nucleotide numbering system for MCOLN1 mutations is presented and the issue of a screening program for the detection of high-risk families in the Jewish Ashkenazi population is discussed.

Architecture and anatomy of the genomic locus encoding the human leukemia-associated transcription factor RUNX1/AML1.

The RUNX1 gene on human chromosome 21q22.12 belongs to the 'runt domain' gene family of transcription factors (also known as AML/CBFA/PEBP2alpha). RUNX1 is a key regulator of hematopoiesis and a frequent target of leukemia associated chromosomal translocations. Here we present a detailed analysis of the RUNX1 locus based on its complete genomic sequence. RUNX1 spans 260 kb and its expression is regulated through two distinct promoter regions, that are 160 kb apart. A very large CpG island complex marks the proximal promoter (promoter-2), and an additional CpG island is located at the 3' end of the gene. Hitherto, 12 different alternatively spliced RUNX1 cDNAs have been identified. Genomic sequence analysis of intron/exon boundaries of these cDNAs has shown that all consist of properly spliced authentic coding regions. This indicates that the large repertoire of RUNX1 proteins, ranging in size between 20-52 kDa, are generated through usage of alternatively spliced exons some of which contain in frame stop codons. The gene's introns are largely depleted of repetitive sequences, especially of the LINE1 family. The RUNX1 locus marks the transition from a ~1 Mb of gene-poor region containing only pseudogenes, to a gene-rich region containing several functional genes. A search for RUNX1 sequences that may be involved in the high frequency of chromosomal translocations revealed that a 555 bp long segment originating in chromosome 11 FLI1 gene was transposed into RUNX1 intron 4.1. This intron harbors the t(8;21) and t(3;21) chromosomal breakpoints involved in acute myeloid leukemia. Interestingly, the FLI1 homologous sequence contains a breakpoint of the t(11;22) translocation associated with Ewing's tumors, and may have a similar function in RUNX1.

Mechanisms for evolving hypervariability: the case of conopeptides.

Hypervariability is a prominent feature of large gene families that mediate interactions between organisms, such as venom-derived toxins or immunoglobulins. In order to study mechanisms for evolution of hypervariability, we examined an EST-generated assemblage of 170 distinct conopeptide sequences from the venoms of five species of marine Conus snails. These sequences were assigned to eight gene families, defined by conserved elements in the signal domain and untranslated regions. Order-of-magnitude differences were observed in the expression levels of individual conopeptides, with five to seven transcripts typically comprising over 50% of the sequenced clones in a given species. The conopeptide precursor alignments revealed four striking features peculiar to the mature peptide domain: (1) an accelerated rate of nucleotide substitution, (2) a bias for transversions over transitions in nucleotide substitutions, (3) a position-specific conservation of cysteine codons within the hypervariable region, and (4) a preponderance of nonsynonymous substitutions over synonymous substitutions. We propose that the first three observations argue for a mutator mechanism targeted to mature domains in conopeptide genes, combining a protective activity specific for cysteine codons and a mutagenic polymerase that exhibits transversion bias, such as DNA polymerase V. The high D:(n)/D:(s) ratio is consistent with positive or diversifying selection, and further analyses by intraspecific/interspecific gene tree contingency tests weakly support recent diversifying selection in the evolution of conopeptides. Since only the most highly expressed transcripts segregate in gene trees according to the feeding specificity of the species, diversifying selection might be acting primarily on these sequences. The combination of a targeted mutator mechanism to generate high variability with the subsequent action of diversifying selection on highly expressed variants might explain both the hypervariability of conopeptides and the large number of unique sequences per species.

Harvesting the human genome: the Israeli perspective.

The post-genome era is at our door, and soon the complete human genome sequence will be available for the next set of goals. Israel is well equipped and skilled to join the worldwide harvest of the human genome, but additional massive government investment is required. This will affect various domains of activity, including the fields of diagnostics and therapeutics. The technologies and know-how described above constitute the basis for future human genome applications in Israel.

Identification of the gene causing mucolipidosis type IV.

Mucolipidosis type IV (MLIV) is an autosomal recessive, neurodegenerative, lysosomal storage disorder characterized by psychomotor retardation and ophthalmological abnormalities including corneal opacities, retinal degeneration and strabismus. Most patients reach a maximal developmental level of 12?15 months. The disease was classified as a mucolipidosis following observations by electron microscopy indicating the lysosomal storage of lipids together with water-soluble, granulated substances. Over 80% of the MLIV patients diagnosed are Ashkenazi Jews, including severely affected and mildly affected patients. The gene causing MLIV was previously mapped to human chromosome 19p13.2-13.3 in a region of approximately 1 cM (ref. 7). Haplotype analysis in the MLIV gene region of over 70 MLIV Ashkenazi chromosomes indicated the existence of two founder chromosomes among 95% of the Ashkenazi MLIV families: a major haplotype in 72% and a minor haplotype in 23% of the MLIV chromosomes (ref. 7, and G.B., unpublished data). The remaining 5% are distinct haplotypes found only in single patients. The basic metabolic defect causing the lysosomal storage in MLIV has not yet been identified. Thus, positional cloning was an alternative to identify the MLIV gene. We report here the identification of a new gene in this human chromosomal region in which MLIV-specific mutations were identified.

Sequence, structure, and evolution of a complete human olfactory receptor gene cluster.

The olfactory receptor (OR) gene cluster on human chromosome 17p13.3 was subjected to mixed shotgun automated DNA sequencing. The resulting 412 kb of genomic sequence include 17 OR coding regions, 6 of which are pseudogenes. Six of the coding regions were discovered only upon genomic sequencing, while the others were previously reported as partial sequences. A comparison of DNA sequences in the vicinity of the OR coding regions revealed a common gene structure with an intronless coding region and at least one upstream noncoding exon. Potential gene control regions including specific pyrimidine:purine tracts and Olf-1 sites have been identified. One of the pseudogenes apparently has evolved into a CpG island. Four extensive CpG islands can be discerned within the cluster, not coupled to specific OR genes. The cluster is flanked at its telomeric end by an unidentified open reading frame (C17orf2) with no significant similarity to any known protein. A high proportion of the cluster sequence (about 60%) belongs to various families of interspersed repetitive elements, with a clear predominance of LINE repeats. The OR genes in the cluster belong to two families and seven subfamilies, which show a relatively high degree of intermixing along the cluster, in seemingly random orientations. This genomic organization may be best accounted for by a complex series of evolutionary events.

Tetrahydropyrimidine derivatives inhibit binding of a Tat-like, arginine-containing peptide, to HIV TAR RNA in vitro.

The ability of a small molecule, 2-methyl,4-carboxy,5-hydroxy-3,4,5,6-tetrahydropyrimidine (THP(A)), which accumulates intracellularly in various streptomyces, to inhibit the interaction of Tat peptide (R52) with TAR RNA is presented. Using gel-shift assay, we found that the inhibition constant Ki of THP(A) is 50-100 nM, which is in the range of the binding constants of Tat peptide and protein. THP(A) is approximately 10(6) times more tightly bound than the free L-arginine. The high binding affinity may be attributed to the special delocalized positive charge on the NCN group and the hydroxyl group at the 5 position of this molecule. A model for THP(A)-TAR interaction, analogous to the arginine guanidinum group-TAR interaction, is presented. The relatively high uptake of THP(A) by mammalian cells warrants in vivo Tat/TAR inhibition studies.

Transcription elongation of the murine ornithine decarboxylase (ODC) gene is regulated in vitro at two downstream elements by different attenuation mechanisms.

Ornithine decarboxylase (ODC) plays an important role in cell proliferation. Its expression is tightly regulated at the mRNA and protein levels and is found to be deregulated in various malignancies. The rapid and dramatic induction of cellular ODC mRNA upon serum addition raised the possibility that a transcriptional attenuation mechanism may be involved in the regulation of ODC gene expression. Using transcription in HeLa nuclear extract and isolated transcription complexes, we have identified two sites of transcription arrest downstream to the transcription start site: Attenuator 1 (Att.1) located at +220, near two repeats of a USF/Myc-Max binding consensus sequence and attenuator 2 (Att.2) located at +1590 near a long stretch of T-residues. The two attenuators exhibit distinct properties as revealed by elongation of briefly initiated and partially purified transcription complexes: Att.1 serves as a transient pause site while arrest at Att.2 is more prolonged. The arrest at both attenuators is modulated by the general elongation factor TFIIS. In a promoter independent transcription system, using partially purified RNA polymerase II, only Att.2 was recognized efficiently. This suggests that the recognition of Att.2 is an intrinsic property of the polymerase while Att.1 recognition has to be facilitated by an auxiliary factor/s.

Minute virus of mice infection modifies cellular transcription elongation.

Our previous observations indicated that upon infection with minute virus of mice (MVM), Ehrlich ascites cells lose a transcription elongation activity which is essential for the readthrough of the MVM attenuator. This was monitored by the ability of extracts from uninfected but not from infected cells to support readthrough of the P4 attenuator when added to partially purified transcription elongation complexes. We have investigated the nature of this change in transcription elongation following MVM infection. In this communication, we show that infection of Ehrlich ascites cells with MVM leads to a general shift in the length of nascent mRNA synthesized in isolated nuclei and separated by sucrose gradients. Furthermore, infection leads to attenuation of transcription of the cellular gene c-fos but not c-myc. We show biochemical evidence to support a model by which, following MVM infection, there is a functional reduction in the activity of a TFIIS-like general transcriptional elongation activity.

Histone H1-mediated inhibition of transcription initiation of methylated templates in vitro.

The detailed mechanism underlying the inhibition of transcription by DNA methylation is still obscure. Chromatin structure has frequently been proposed as a role player in this mechanism. Histone H1 is a known key element in the formation and stabilization of chromatin fibers. We describe here experiments designed to examine the effect of DNA methylation on the binding of histone H1 to DNA and the consequent inhibitory effect of the bound histone H1 on in vitro transcription. The results of these experiments showed a clear preferential binding of histone H1 to methylated DNA as compared with unmethylated DNA. The in vitro transcription assay indicated that transcription of methylated templates was inhibited at a lower histone H1/DNA ratio than of unmethylated templates, and that the extent of inhibition depends on the density of methyl groups in the promoter region. This inhibition of in vitro transcription was alleviated efficiently by methylated competitor DNA, whereas, under similar conditions, almost no effect was observed with unmethylated competitor. Experiments designed to pinpoint the stage in the transcription process that was suppressed by the preferred binding of histone H1 to methylated template revealed that inhibition occurred at the initiation and not at the elongation level.