Publications

Protocells at life’s origin are often conceived as bilayer-enclosed precursors of life, whose self-reproduction rests on the early advent of replicating catalytic biopolymers. This Perspective describes an alternative scenario, wherein reproducing nanoscopic lipid micelles with catalytic capabilities were forerunners of biopolymer-containing protocells. This postulate gains considerable support from experiments describing micellar catalysis and autocatalytic proliferation, and, more recently, from reports on cross-catalysis in mixed micelles that lead to life-like steady-state dynamics. Such results, along with evidence for micellar prebiotic compatibility, synergize with predictions of our chemically stringent computer-simulated model, illustrating how mutually catalytic lipid networks may enable micellar compositional reproduction that could underlie primal selection and evolution. Finally, we highlight studies on how endogenously catalysed lipid modifications could guide further protocellular complexification, including micelle to vesicle transition and monomer to biopolymer progression. These portrayals substantiate the possibility that protocellular evolution could have been seeded by pre-RNA lipid assemblies.

Circulating inflammatory markers are essential to human health and disease, and they are often dysregulated or malfunctioning in cancers as well as in cardiovascular, metabolic, immunologic and neuropsychiatric disorders. However, the genetic contribution to the physiological variation of levels of circulating inflammatory markers is largely unknown. Here we report the results of a genome-wide genetic study of blood concentration of ten cytokines, including the hitherto unexplored calcium-binding protein (S100B). The study leverages a unique sample of neonatal blood spots from 9,459 Danish subjects from the iPSYCH initiative. We estimate the SNP-heritability of marker levels as ranging from essentially zero for Erythropoietin (EPO) up to 73% for S100B. We identify and replicate 16 associated genomic regions (p

"The Lipid World" was published in 2001, stemming from a highly effective collaboration with David Deamer during a sabbatical year 20 years ago at the Weizmann Institute of Science in Israel. The present review paper highlights the benefits of this scientific interaction and assesses the impact of the lipid world paper on the present understanding of the possible roles of amphiphiles and their assemblies in the origin of life. The lipid world is defined as a putative stage in the progression towards life's origin, during which diverse amphiphiles or other spontaneously aggregating small molecules could have concurrently played multiple key roles, including compartment formation, the appearance of mutually catalytic networks, molecular information processing, and the rise of collective self-reproduction and compositional inheritance. This review brings back into a broader perspective some key points originally made in the lipid world paper, stressing the distinction between the widely accepted role of lipids in forming compartments and their expanded capacities as delineated above. In the light of recent advancements, we discussed the topical relevance of the lipid worldview as an alternative to broadly accepted scenarios, and the need for further experimental and computer-based validation of the feasibility and implications of the individual attributes of this point of view. Finally, we point to possible avenues for exploring transition paths from small molecule-based noncovalent structures to more complex biopolymer-containing proto-cellular systems.

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

We studied the simulated replication and growth of prebiotic vesicles composed of 140 phospholipids and cholesterol using our R-GARD (Real Graded Autocatalysis Replication Domain) formalism that utilizes currently extant lipids that have known rate constants of lipid-vesicle interactions from published experimental data. R-GARD normally modifies kinetic parameters of lipid-vesicle interactions based on vesicle composition and properties. Our original R-GARD model tracked the growth and division of one vesicle at a time in an environment with unlimited lipids at a constant concentration. We explore here a modified model where vesicles compete for a finite supply of lipids. We observed that vesicles exhibit complex behavior including initial fast unrestricted growth, followed by intervesicle competition for diminishing resources, then a second growth burst driven by better-adapted vesicles, and ending with a final steady state. Furthermore, in simulations without kinetic parameter modifications (invariant kinetics), the initial replication was an order of magnitude slower, and vesicles' composition variability at the final steady state was much lower. The complex kinetic behavior was not observed either in the previously published R-GARD simulations or in additional simulations presented here with only one lipid component. This demonstrates that both a finite environment (inducing selection) and multiple components (providing variation for selection to act upon) are crucial for portraying evolution-like behavior. Such properties can improve survival in a changing environment by increasing the ability of early protocellular entities to respond to rapid environmental fluctuations likely present during abiogenesis both on Earth and possibly on other planets. This in silico simulation predicts that a relatively simple in vitro chemical system containing only lipid molecules might exhibit properties that are relevant to prebiotic processes.

Background: A key challenge in the realm of human disease research is next generation sequencing (NGS) interpretation, whereby identified filtered variant-harboring genes are associated with a patient's disease phenotypes. This necessitates bioinformatics tools linked to comprehensive knowledgebases. The GeneCards suite databases, which include GeneCards (human genes), MalaCards (human diseases) and PathCards (human pathways) together with additional tools, are presented with the focus on MalaCards utility for NGS interpretation as well as for large scale bioinformatic analyses. Results: VarElect, our NGS interpretation tool, leverages the broad information in the GeneCards suite databases. MalaCards algorithms unify disease-related terms and annotations from 69 sources. Further, MalaCards defines hierarchical relatedness- aliases, disease families, a related diseases network, categories and ontological classifications. GeneCards and MalaCards delineate and share a multi-tiered, scored gene-disease network, with stringency levels, including the definition of elite status- high quality gene-disease pairs, coming from manually curated trustworthy sources, that includes 4500 genes for 8000 diseases. This unique resource is key to NGS interpretation by VarElect. VarElect, a comprehensive search tool that helps infer both direct and indirect links between genes and user-supplied disease/phenotype terms, is robustly strengthened by the information found in MalaCards. The indirect mode benefits from GeneCards' diverse gene-to-gene relationships, including SuperPathsintegrated biological pathways from 12 information sources. We are currently adding an important information layer in the form of "disease SuperPaths", generated from the gene-disease matrix by an algorithm similar to that previously employed for biological pathway unification. This allows the discovery of novel gene-disease and disease-disease relationships. The advent of whole genome sequencing necessitates

2016 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.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.

Postgenomics data are produced in large volumes by life sciences and clinical applications of novel omics diagnostics and therapeutics for precision medicine. To move from "data-to-knowledge-to-innovation," a crucial missing step in the current era is, however, our limited understanding of biological and clinical contexts associated with data. Prominent among the emerging remedies to this challenge are the gene set enrichment tools. This study reports on GeneAnalytics (TM) (geneanalytics.genecards.org), a comprehensive and easy-to-apply gene set analysis tool for rapid contextualization of expression patterns and functional signatures embedded in the postgenomics Big Data domains, such as Next Generation Sequencing (NGS), RNAseq, and microarray experiments. GeneAnalytics' differentiating features include in-depth evidence-based scoring algorithms, an intuitive user interface and proprietary unified data. GeneAnalytics employs the LifeMap Science's GeneCards suite, including the GeneCards (R)-the human gene database; the MalaCards-the human diseases database; and the PathCards-the biological pathways database. Expression-based analysis in GeneAnalytics relies on the LifeMap Discovery (R)-the embryonic development and stem cells database, which includes manually curated expression data for normal and diseased tissues, enabling advanced matching algorithm for gene-tissue association. This assists in evaluating differentiation protocols and discovering biomarkers for tissues and cells. Results are directly linked to gene, disease, or cell "cards" in the GeneCards suite. Future developments aim to enhance the GeneAnalytics algorithm as well as visualizations, employing varied graphical display items. Such attributes make GeneAnalytics a broadly applicable postgenomics data analyses and interpretation tool for translation of data to knowledge-based innovation in various Big Data fields such as precision medicine, ecogenomics, nutrigenomics, pharmacogenomics, vaccinomic

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 phen

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 pseudoge

Autophagy dysfunction has been implicated in a group of progressive neurodegenerative diseases, and has been reported to play a major role in the pathogenesis of these disorders. We have recently reported a recessive mutation in TECPR2, an autophagy-implicated WD repeat-containing protein, in five individuals with a novel form of monogenic hereditary spastic paraparesis (HSP). We found that diseased skin fibroblasts had a decreased accumulation of the autophagy-initiation protein MAP1LC3B/LC3B, and an attenuated delivery of both LC3B and the cargo-recruiting protein SQSTM1/p62 to the lysosome where they are subject to degradation. The discovered TECPR2 mutation reveals for the first time a role for aberrant autophagy in a major class of Mendelian neurodegenerative diseases, and suggests mechanisms by which impaired autophagy may impinge on a broader scope of neurodegeneration.

Comprehensive disease classification, integration and annotation are crucial for biomedical discovery. At present, disease compilation is incomplete, heterogeneous and often lacking systematic inquiry mechanisms. We introduce MalaCards, an integrated database of human maladies and their annotations, modeled on the architecture and strategy of the GeneCards database of human genes. MalaCards mines and merges 44 data sources to generate a computerized card for each of 16 919 human diseases. Each MalaCard contains disease-specific prioritized annotations, as well as inter-disease connections, empowered by the GeneCards relational database, its searches and GeneDecks set analyses. First, we generate a disease list from 15 ranked sources, using disease-name unification heuristics. Next, we use four schemes to populate MalaCards sections: (i) directly interrogating disease resources, to establish integrated disease names, synonyms, summaries, drugs/therapeutics, clinical features, genetic tests and anatomical context; (ii) searching GeneCards for related publications, and for associated genes with corresponding relevance scores; (iii) analyzing disease-associated gene sets in GeneDecks to yield affiliated pathways, phenotypes, compounds and GO terms, sorted by a composite relevance score and presented with GeneCards links; and (iv) searching within MalaCards itself, e. g. for additional related diseases and anatomical context. The latter forms the basis for the construction of a disease network, based on shared MalaCards annotations, embodying associations based on etiology, clinical features and clinical conditions. This broadly disposed network has a power-law degree distribution, suggesting that this might be an inherent property of such networks. Work in progress includes hierarchical malady classification, ontological mapping and disease set analyses, striving to make MalaCards an even more effective tool for biomedical research.

We propose an automaton, a theoretical framework that demonstrates how to improve the yield of the synthesis of branched chemical polymer reactions. This is achieved by separating substeps of the path of synthesis into compartments. We use chemical containers (chemtainers) to carry the substances through a sequence of fixed successive compartments. We describe the automaton in mathematical terms and show how it can be configured automatically in order to synthesize a given branched polymer target. The algorithm we present finds an optimal path of synthesis in linear time. We discuss how the automaton models compartmentalized structures found in cells, such as the endoplasmic reticulum and the Golgi apparatus, and we show how this compartmentalization can be exploited for the synthesis of branched polymers such as oligosaccharides. Lastly, we show examples of artificial branched polymers and discuss how the automaton can be configured to synthesize them with maximal yield.

In this paper we explore the question of whether there is an optimal set up for a putative prebiotic system leading to open-ended evolution (OEE) of the events unfolding within this system. We do so by proposing two key innovations. First, we introduce a new index that measures OEE as a function of the likelihood of events unfolding within a universe given its initial conditions. Next, we apply this index to a variant of the graded autocatalysis replication domain (GARD) model, Segre et al. (P Natl Acad Sci USA 97(8):4112-4117, 2000; Markovitch and Lancet Artif Life 18(3), 2012), and use it to study - under a unified and concise prebiotic evolutionary framework - both a variety of initial conditions of the universe and the OEE of species that evolve from them.

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 similar to 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, similar to 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.

Paenibacillus dendritiformis is a Gram-positive, soil-dwelling, spore-forming social microorganism. An intriguing collective faculty of this strain is manifested by its ability to switch between different morphotypes, such as the branching ( T) and the chiral (C) morphotypes. Here we report the 6.3-Mb draft genome sequence of the P. dendritiformis C454 chiral morphotype.

It is widely accepted that autocatalysis constitutes a crucial facet of effective replication and evolution (e.g., in Eigen's hypercycle model). Other models for early evolution (e.g., by Dyson, Ganti, Varela, and Kauffman) invoke catalytic networks, where cross-catalysis is more apparent. A key question is how the balance between auto- (self-) and cross- (mutual) catalysis shapes the behavior of model evolving systems. This is investigated using the graded autocatalysis replication domain (GARD) model, previously shown to capture essential features of reproduction, mutation, and evolution in compositional molecular assemblies. We have performed numerical simulations of an ensemble of GARD networks, each with a different set of lognormally distributed catalytic values. We asked what is the influence of the catalytic content of such networks on beneficial evolution. Importantly, a clear trend was observed, wherein only networks with high mutual catalysis propensity (p(mc)) allowed for an augmented diversity of composomes, quasi-stationary compositions that exhibit high replication fidelity. We have reexamined a recent analysis that showed meager selection in a single GARD instance and for a few nonstationaty target compositions. In contrast, when we focused here on compotypes (clusters of composomes) as targets for selection in populations of compositional assemblies, appreciable selection response was observed for a large portion of the networks simulated. Further, stronger selection response was seen for high p(mc) values. Our simulations thus demonstrate that GARD can help analyze important facets of evolving systems, and indicate that excess mutual catalysis over self-catalysis is likely to be important for the emergence of molecular systems capable of evolutionlike behavior.

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.22x10-(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.-Alkelai, A., Lupoli, S., Greenbaum, L., Giegling, I., Kohn, Y., Sarner-Kanyas, K., Ben-Asher, E., Lancet, D., Rujescu, D., Macciardi, F., Lerer, B. Identification of new schizophrenia susceptibility

Copy-number variations (CNVs) are widespread in the human genome, but comprehensive assignments of integer locus copy-numbers (i.e., copy-number genotypes) that, for example, enable discrimination of homozygous from heterozygous CNVs, have remained challenging. Here we present CopySeq, a novel computational approach with an underlying statistical framework that analyzes the depth-of-coverage of high-throughput DNA sequencing reads, and can incorporate paired-end and breakpoint junction analysis based CNV-analysis approaches, to infer locus copy-number genotypes. We benchmarked CopySeq by genotyping 500 chromosome 1 CNV regions in 150 personal genomes sequenced at low-coverage. The assessed copy-number genotypes were highly concordant with our performed qPCR experiments (Pearson correlation coefficient 0.94), and with the published results of two microarray platforms (95-99% concordance). We further demonstrated the utility of CopySeq for analyzing gene regions enriched for segmental duplications by comprehensively inferring copy-number genotypes in the CNV-enriched >800 olfactory receptor (OR) human gene and pseudogene loci. CopySeq revealed that OR loci display an extensive range of locus copy-numbers across individuals, with zero to two copies in some OR loci, and two to nine copies in others. Among genetic variants affecting OR loci we identified deleterious variants including CNVs and SNPs affecting similar to 15% and similar to 20% of the human OR gene repertoire, respectively, implying that genetic variants with a possible impact on smell perception are widespread. Finally, we found that for several OR loci the reference genome appears to represent a minor-frequency variant, implying a necessary revision of the OR repertoire for future functional studies. CopySeq can ascertain genomic structural variation in specific gene families as well as at a genome-wide scale, where it may enable the quantitative evaluation of CNVs in genome-wide association studies in

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.-Torri, F., Akelai, A., Lupoli, S., Sironi, M., Amann-Zalcenstein, D., Fumagalli, M., Dal Fiume, C., BenAsher, E., Kanyas, K., Cagliani, R., Cozzi, P., Trombetti, G., Lievers, L. S., Salvi, E., Orro, A., Beckmann, J. S., Lancet, D., Kohn, Y., Milanesi, L., Ebste

Background: An important facet of early biological evolution is the selection of chiral enantiomers for molecules such as amino acids and sugars. The origin of this symmetry breaking is a long-standing question in molecular evolution. Previous models addressing this question include particular kinetic properties such as autocatalysis or negative cross catalysis. Results: We propose here a more general kinetic formalism for early enantioselection, based on our previously described Graded Autocatalysis Replication Domain (GARD) model for prebiotic evolution in molecular assemblies. This model is adapted here to the case of chiral molecules by applying symmetry constraints to mutual molecular recognition within the assembly. The ensuing dynamics shows spontaneous chiral symmetry breaking, with transitions towards stationary compositional states (composomes) enriched with one of the two enantiomers for some of the constituent molecule types. Furthermore, one or the other of the two antipodal compositional states of the assembly also shows time-dependent selection. Conclusion: It follows that chiral selection may be an emergent consequence of early catalytic molecular networks rather than a prerequisite for the initiation of primeval life processes. Elaborations of this model could help explain the prevalent chiral homogeneity in present-day living cells.

The Graded Autocatalysis Replication Domain (GARD) model describes an origin of life scenario which involves non-covalent compositional assemblies, made of monomeric mutually catalytic molecules. GARD constitutes an alternative to informational biopolymers as a mechanism of primordial inheritance. In the present work, we examined the effect of mutations, one of the most fundamental mechanisms for evolution, in the context of the networks of mutual interaction within GARD prebiotic assemblies. We performed a systematic analysis analogous to single and double gene deletions within GARD. While most deletions have only a small effect on both growth rate and molecular composition of the assemblies, similar to 10% of the deletions caused lethality, or sometimes showed enhanced fitness. Analysis of 14 different network properties on 2,000 different GARD networks indicated that lethality usually takes place when the deleted node has a high molecular count, or when it is a catalyst for such node. A correlation was also found between lethality and node degree centrality, similar to what is seen in real biological networks. Addressing double knockout mutations, our results demonstrate the occurrence of both synthetic lethality and extragenic suppression within GARD networks, and convey an attempt to correlate synthetic lethality to network node-pair properties. The analyses presented help establish GARD as a workable alternative prebiotic scenario, suggesting that life may have begun with large molecular networks of low fidelity, that later underwent evolutionary compaction and fidelity augmentation.

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. (C) 2009 Wiley-Liss, Inc.

The sense of smell is a complex molecular device, encompassing several hundred olfactory receptor proteins (ORs). These receptors, encoded by the largest human gene superfamily, integrate odorant signals into an accurate 'odor image' in the brain. Widespread phenotypic diversity in human olfaction is, in part, attributable to prevalent genetic variation in OR genes, owing to copy number variation, deletion alleles and deleterious single nucleotide polymorphisms. The development of new genomic tools, including next generation sequencing and CNV assays; provides opportunities to characterize the genetic variations of this system. The advent of large-scale functional screens of expressed ORs, combined with genetic association studies, has the potential to link variations in ORs to human chemosensory phenotypes. This promises to provide a genome-wide view of human olfaction, resulting in a deeper understanding of personalized odor coding, with the potential to decipher flavor and fragrance preferences.

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. (c) 2007 Wiley-Liss, Inc.

Background: Improvements in genome sequence annotation revealed discrepancies in the original probeset/gene assignment in Affymetrix microarray and the existence of differences between annotations and effective alignments of probes and transcription products. In the current generation of Affymetrix human GeneChips, most probesets include probes matching transcripts from more than one gene and probes which do not match any transcribed sequence. Results: We developed a novel set of custom Chip Definition Files (CDF) and the corresponding Bioconductor libraries for Affymetrix human GeneChips, based on the information contained in the GeneAnnot database. GeneAnnot-based CDFs are composed of unique custom-probesets, including only probes matching a single gene. Conclusion: GeneAnnot-based custom CDFs solve the problem of a reliable reconstruction of expression levels and eliminate the existence of more than one probeset per gene, which often leads to discordant expression signals for the same transcript when gene differential expression is the focus of the analysis. GeneAnnot CDFs are freely distributed and fully compliant with Affymetrix standards and all available software for gene expression analysis. The CDF libraries are available from http://www.xlab.unimo.it/GA CDF, along with supplementary information (CDF libraries, installation guidelines and R code, CDF statistics, and analysis results).

The coevolution of environment and living organisms is well known in nature. Here, it is suggested that similar processes can take place before the onset of life, where protocellular entities, rather than full-fledged living systems, coevolve along with their surroundings. Specifically, it is suggested that the chemical composition of the environment may have governed the chemical repertoire generated within molecular assemblies, compositional protocells, while compounds generated within these protocells altered the chemical composition of the environment. We present an extension of the graded autocatalysis replication domain (GARD) model-the environment exchange polymer GARD (EE-GARD) model. In the new model, molecules, which are formed in a protocellular assembly, may be exported to the environment that surrounds the protocell. Computer simulations of the model using an infinite-sized environment showed that EE-GARD assemblies may assume several distinct quasi-stationary compositions (composomes), similar to the observations in previous variants of the GARD model. A statistical analysis suggested that the repertoire of composomes manifested by the assemblies is independent of time. In simulations with a finite environment, this was not the case. Composomes, which were frequent in the early stages of the simulation disappeared, while others emerged. The change in the frequencies of composomes was found to be correlated with changes induced on the environment by the assembly. The EE-GARD model is the first GARD model to portray a possible time evolution of the composomes repertoire.

Objectives To investigate the role of genes encoding regulators of G protein signaling in early therapeutic response to antipsychotic drugs and in susceptibility to drug-induced extrapyramidal symptoms. As regulators of G protein signaling and regulators of G protein signaling-like proteins play a pivotal role in dopamine receptor signaling, genetically based, functional variation could contribute to interindividual variability in therapeutic and adverse effects. Methods Consecutively hospitalized, psychotic patients with Diagnostic and Statistical Manual of Mental Disorder-IV schizophrenia (n=121) were included in the study if they received treatment with typical antipsychotic medication (n=72) or typical antipsychotic drugs and risperidone (n=49) for at least 2 weeks. Clinical state and adverse effects were rated at baseline and after 2 weeks. Twenty-four single nucleotide polymorphisms were genotyped in five regulators of G protein signaling genes. Results None of the single nucleotide polymorphisms were related to clinical response to antipsychotic treatment at 2 weeks. Five out of six single nucleotide polymorphisms within or flanking the RGS2 gene were nominally associated with development or worsening of parkinsonian symptoms (PARK+) as measured by the Simpson Angus Scale, one of them after correction for multiple testing (rs4606, P=0.002). A GCCTG haplotype encompassing tagging single nucleotide polymorphisms within and flanking RGS2 was significantly overrepresented among PARK+ compared with PARK- patients (0.23 vs. 0.08, P=0.003). A second, 'protective', GTGCA haplotype was significantly overrepresented in PARK- patients (0.13 vs. 0.30, P=0.009). Both haplotype associations survive correction for multiple testing. Conclusions Subject to replication, these findings suggest that genetic variation in the RGS2 gene is associated with susceptibility to extrapyramidal symptoms induced by antipsychotic drugs.

The MDM2 protein is an ubiquitin ligase that plays a critical role in regulating the levels and activity of the p53 protein, which is a central tumor suppressor. A SNP in the human MDM2 gene (SNP309 T/G) occurs at frequencies dependent on demographic history and has been shown to have important differential effects on the activity of the MDM2 and p53 proteins and to associate with altered risk for the development of several cancers. In this report, the haplotype structure of the MDM2 gene is determined by using 14 different SNPs across the gene from three different population samples: Caucasians, African Americans, and the Ashkenazi Jewish ethnic group. The results presented in this report indicate that there is a substantially reduced variability of the deleterious SNP309 G allele haplotype in all three populations studied, whereas multiple common T allele haplotypes were found in all three populations. This observation, coupled with the relatively high frequency of the G allele haplotype in both and Caucasian and Ashkenazi Jewish population data sets, suggests that this haplotype could have undergone a recent positive selection sweep. An entropy-based selection test is presented that explicitly takes into account the correlations between different SNPs, and the analysis of MDM2 reveals a significant departure from the standard assumptions of selective neutrality.

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 similar to 7Mb region around the strong linkage peak at 136.2Mb 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.

Background: Olfactory receptors (ORs), the largest mammalian gene superfamily (900-1400 genes), has > 50% pseudogenes in humans. While most of these inactive genes are identified via coding frame (nonsense) disruptions, seemingly intact genes may also be inactive due to other deleterious (missense) mutations. An ultimate assessment of the actual size of the functional human OR repertoire thus requires an accurate distinction between genes and pseudogenes. Results: To characterize inactive ORs with intact open reading frame, we have developed a probabilistic Classifier for Olfactory Receptor Pseudogenes (CORP). This algorithm is based on deviations from a functionally crucial consensus, constituting sixty highly conserved positions identified by a comparison of two evolutionarily-constrained OR repertoires ( mouse and dog) with a small pseudogene fraction. We used a logistic regression analysis to assign appropriate coefficients to the conserved position and thus achieving maximal separation between active and inactive ORs. Consequently, the algorithms identified only 5% of the mouse functional ORs as pseudogenes, setting an upper limit of 0.05 to the false positive detection. Finally we used this algorithm to classify the 384 purportedly intact human OR genes. Of these, 135 were predicted as likely encoding non-functional proteins, and 38 were segregating between active and inactive forms due to missense polymorphisms. Conclusion: We demonstrated that the CORP algorithm is capable to distinguish between functional and non-functional OR genes with high precision even when the encoded protein would differ by a single amino acid. Using the CORP algorithm, we predict that similar to 70% of human OR genes are likely non-functional pseudogenes, a much higher number than hitherto suspected. The method we present may be employed for better annotation of inactive members in other gene families as well.

Background: Olfactory receptors (ORs) are the largest gene family in the human genome. Although they are expected to be expressed specifically in olfactory tissues, some ectopic expression has been reported, with special emphasis on sperm and testis. The present study systematically explores the expression patterns of OR genes in a large number of tissues and assesses the potential functional implication of such ectopic expression. Results: We analyzed the expression of hundreds of human and mouse OR transcripts, via EST and microarray data, in several dozens of human and mouse tissues. Different tissues had specific, relatively small OR gene subsets which had particularly high expression levels. In testis, average expression was not particularly high, and very few highly expressed genes were found, none corresponding to ORs previously implicated in sperm chemotaxis. Higher expression levels were more common for genes with a non-OR genomic neighbor. Importantly, no correlation in expression levels was detected for human-mouse orthologous pairs. Also, no significant difference in expression levels was seen between intact and pseudogenized ORs, except for the pseudogenes of subfamily 7E which has undergone a human-specific expansion. Conclusion: The OR superfamily as a whole, show widespread, locus-dependent and heterogeneous expression, in agreement with a neutral or near neutral evolutionary model for transcription control. These results cannot reject the possibility that small OR subsets might play functional roles in different tissues, however considerable care should be exerted when offering a functional interpretation for ectopic OR expression based only on transcription information.

We hypothesize that life began not with the first self-reproducing molecule or metabolic network, but as a prebiotic ecology of co-evolving populations of macromolecular aggregates (composomes). Each composome species had a particular molecular composition resulting from molecular complementarity among environmentally available prebiotic compounds. Natural selection acted on composomal species that varied in properties and functions such as stability, catalysis, fission, fusion and selective accumulation of molecules from solution. Fission permitted molecular replication based on composition rather than linear structure, while fusion created composomal variability. Catalytic functions provided additional chemical novelty resulting eventually in autocatalytic and mutually catalytic networks within composomal species. Composomal autocatalysis and interdependence allowed the Darwinian co-evolution of content and control (metabolism). The existence of chemical interfaces within complex composomes created linear templates upon which self-reproducing molecules (such as RNA) could be synthesized, permitting the evolution of informational replication by molecular templating. Mathematical and experimental tests are proposed.

We have systematically examined the domain composition across a comprehensive set of tissue-specific, midrange and housekeeping genes as defined by their mode of expression in 52 normal mouse tissues. We show a definite correlation between the number of domains and the degree of tissue specificity. This trend is further supported by a novel analysis involving the time of origin of each domain. Genes containing metazoan-specific domains are more prevalent in signal transduction and cell-communication pathways, and are depleted in primary metabolism. Our analyses suggest that highly modular gene products have been recruited for tissue-specific functions that are required in complex organisms.

The basic Graded Autocatalysis Replication Domain (GARD) model consists of a repertoire of small molecules, typically amphiphiles, which join and leave a non-covalent micelle-like assembly. Its replication behavior is due to occasional fission, followed by a homeostatic growth process governed by the assembly's composition. Limitations of the basic GARD model are its small finite molecular repertoire and the lack of a clear path from a 'monomer world' towards polymer-based living entities. We have now devised an extension of the model (polymer GARD or P-GARD), where a monomer-based GARD serves as a 'scaffold' for oligomer formation, as a result of internal chemical rules. We tested this concept with computer simulations of a simple case of monovalent monomers, whereby more complex molecules (dimers) are formed internally, in a manner resembling biosynthetic metabolism. We have observed events of dimer 'take-over' - the formation of compositionally stable, replication-prone quasi stationary states (composomes) that have appreciable dimer content. The appearance of novel metabolism-like networks obeys a time-dependent power law, reminiscent of evolution under punctuated equilibrium. A simulation under constant population conditions shows the dynamics of takeover and extinction of different composomes, leading to the generation of different population distributions. The P-GARD model offers a scenario whereby biopolymer formation may be a result of rather than a prerequisite for early life-like processes.

Quantitative variation in gene expression in humans is the outcome of various factors, including differences in genetic background, gender, age, and environment. However, the extent of the influence of these factors on gene expression is not clear. We attempted to address this issue by carrying out gene expression profiling in blood leukocytes with 13 individuals (including 5 pairs of monozygotic twins) on 10,000 genes using HG-U95Av2 oligonucleotide microarrays. The proportion of differentially expressed genes between monozygotic twins was low (up to 1.76%). Most of the variations belonged to the least variable category. These genes, exhibiting "random variations," did not show clear preference to any functional class, although " signaling and communication" and " immune and related functions" generally topped the list. The extent of variation in gene expression increased in comparisons between unrelated individuals (up to 14.13%). Most of the genes (89%) exhibiting random variations in twins also varied in expression in unrelated individuals. As with twins, signaling and communication topped the list, and substantial variations were observed in all three categories: least variable, moderately variable, and most variable. An important outcome of this study was that the housekeeping genes were nearly insensitive to random variations but appeared to be more susceptible to genetic differences. However, the highly expressed housekeeping genes exhibited low variation and appeared to be insensitive to all known factors. Gene expression profiling in monozygotic twins can provide useful data for the assessment of natural variation in gene expression in humans.

Genomic segments that do not code for proteins yet show high conservation among vertebrates have recently been identified by various computational methodologies. We refer to them as ANCORs (ancestral non-coding conserved regions). The frequency of individual ANCORs within the genome, along with their (correlated) inter-species identity scores, helps in assessing the probability that they function in transcription regulation or RNA coding.

Autoimmune diseases seem to have strong genetic attributes, and are affected to some extent by shared susceptibility loci. The latter potentially amount to hundreds of candidate genes (CG), creating the need for a prioritization strategy in genetic association studies. To form such a strategy, 26 autoimmune-related CG were genotyped for a total of 72 single nucleotide polymorphisms (SNPs) in three distinct Israeli ethnic populations: Ashkenazi Jews, Sephardic Jews and Arabs. Four quantitative criteria reflecting population stratification were analyzed: allele frequencies, haplotype frequencies, the F-st statistic for homozygotes distribution and linkage disequilibrium extents. According to the consequent interpopulation genomic diversity profiles, the genes were classified into conserved, intermediate and diversified gene groups. Our results demonstrate a correlation between the biological role of autoimmune-related CG and their interpopulation diversity profiles as classified by the different analyses. Annotation analysis suggests that genes more readily influenced by environmental conditions, such as immunological mediators, are 'population specific'. Conversely, genes showing genetic conservation across all populations are characterized by apoptotic and cleaving functions. We suggest a research strategy by which CG association studies should focus first on likely conserved gene categories, to increase the likelihood of attaining significant results and promote the development of gene-based therapies.

Molecular Chirality is of central interest in biological studies because enantiomeric compounds, while indistinguishable by most inanimate systems, show profoundly different properties in biochemical environments. Enantioselective separation methods, based on the differential recognition of two optical isomers by a chiral selector, have been amply documented. Also, great effort has been directed towards a theoretical understanding of the fundamental mechanisms underlying the chiral recognition process. Here we report a comprehensive data examination of enantioseparation measurements for over 72,000 chiral selector-selectand pairs from the chiral selection compendium CHIRBASE. The distribution of alpha = k'(D)/k'(L) values was found to follow a power law, equivalent to an exponential decay for chiral differential free energies. This observation is experimentally relevant in terms of the number of different individual or combinatorial selectors that need to be screened in order to observe a values higher than a preset minimum. A string model for enantiorecognition (SMED) formalism is proposed to account for this observation on the basis of an extended Ogston three-point interaction model. Partially overlapping molecular interaction domains are analyzed in terms of a string complementarity model for ligand-receptor complementarity. The results suggest that chiral selection statistics may be interpreted in terms of more general concepts related to biomolecular recognition. (C) 2004 Wiley-Liss, Inc.

While the last century brought an exquisite understanding of the molecular basis of life, very little is known about the detailed chemical mechanisms that afforded the emergence of life on early earth. There is a broad agreement that the problem lies in the realm of chemistry, and likely resides in the formation and mutual interactions of carbon-based molecules in aqueous medium. Yet, present-day experimental approaches can only capture the synthesis and behavior of a few molecule types at a time. On the other hand, experimental simulations of prebiotic syntheses, as well as chemical analyses of carbonaceous meteorites, suggest that the early prebiotic hydrosphere contained many thousands of different compounds. The present paper explores the idea that given the limitations of test-tube approaches with regards to such a 'random chemistry' scenario, an alternative mode of analysis should be pursued. It is argued that as computational tools for the reconstruction of molecular interactions improve rapidly, it may soon become possible to perform adequate computer-based simulations of prebiotic evolution. We thus propose to launch a computational origin of life endeavor (http://ool.weizmann.ac.il/CORE), involving computer simulations of realistic complex prebiotic chemical networks. In the present paper we provide specific examples, based on a novel algorithmic approach, which constitutes a hybrid of molecular dynamics and stochastic chemistry. As one potential solution for the immense hardware requirements dictated by this approach, we have begun to implement an idle CPU harvesting scheme, under the title ool@home.

The question of the origin of life is addressed by artificial life research, particularly in the realm of artificial chemistry. Such artificial chemistry is described by our Graded Autocatalysis Replication Domain (GARD) model. GARD depicts an unorthodox scenario suggested for emergence of life - the 'lipid world'. The model concerns molecular assemblies with mutual catalysis in an environment containing a plethora of molecular species. Many aspects of GARD were amply discussed. Here we concentrate on the importance of size constraints as depicted by the basic model and several of its variants. Occasional fission of a GARD assembly, which restricts the assembly size, is crucial for generating compositional quasi-stationary states ('composomes'). In a spatial version of GARD, bounded environments yield spontaneous emergence of different ecologies. Limiting the size of a population of GARD assemblies gives rise to a complex population dynamics. The last example, with possible wider impact to chemistry and nano-technology, suggests that size limit can give rise to spontaneous symmetry breaking. This latter result is compared to the classic Frank's model for homo-chirality, which requires explicit inhibition. We conclude that size restrictions are fundamental in the field of origin of life and artificial life, not only in order to facilitate evolutionary processes, as previously suggested, but also, for augmenting the dynamics portrayed by different scenarios and models.

The olfactory receptor ( OR) genes constitute the largest gene family in mammalian genomes. Humans have 11,000 OR genes, of which only similar to 40% have an intact coding region and are therefore putatively functional. In contrast, the fraction of intact OR genes in the genomes of the great apes is significantly greater (68% - 72%), suggesting that selective pressures on the OR repertoire vary among these species. We have examined the evolutionary forces that shaped the OR gene family in humans and chimpanzees by resequencing 20 OR genes in 16 humans, 16 chimpanzees, and one orangutan. We compared the variation at the OR genes with that at intergenic regions. In both humans and chimpanzees, OR pseudogenes seem to evolve neutrally. In chimpanzees, patterns of variability are consistent with purifying selection acting on intact OR genes, whereas, in humans, there is suggestive evidence for positive selection acting on intact OR genes. These observations are likely due to differences in lifestyle, between humans and great apes, that have led to distinct sensory needs.

We propose an algorithm for use with multisensor systems that is capable of the following: (a) identify an analyte independently of its concentration; (b) estimate the concentration of the analyte, even if the system was not previously exposed to this concentration; (c) tell when an analyte is of a chemical type not previously presented to the system. The algorithm, based upon recent work of Hopfield, uses the multiplicity of sensors explicitly, and is intuitive and easy to implement. We have tested it against real data, and it exhibits high quality performance. (C) 2003 Elsevier Science B.V. All rights reserved.

Of more than 1,000 human olfactory receptor genes, more than half seem to be pseudogenes. We investigated whether the most recent of these disruptions might still segregate with the intact form by genotyping 51 candidate genes in 189 ethnically diverse humans. The results show an unprecedented prevalence of segregating pseudogenes, identifying one of the most pronounced cases of functional population diversity in the human genome.

We propose a setup for an odor communication system. Its different parts are described, and ways to realize them are outlined. Our scheme enables an output device-the whiffer-to release an imitation of an odorant read in by an input device-the sniffer-upon command. The heart of the system is the novel algorithmic scheme that makes the scheme feasible. We are currently at work researching and developing some of the components that constitute the algorithm, and we hope that the description of the overall scheme in this paper will help to get other groups to join in this effort. (C) 2002 Elsevier Science Ltd. All rights reserved.

In addition to the visible complexity expressed in the morphogenesis of multicellular organisms, two levels of microscopic complexity may be discerned within every living cell. The first level is related to covalently bonded structures, namely molecules. The second level has to do with the generation of non-covalent molecular assemblies. Origin of life research has largely focused on the first complexity level, i.e. the appearance of covalent biopolymers. We present a life emergence scenario based mainly on the second complexity level. We argue that homeostatic molecular ensembles, for which we have coined the term "mesobiotic," have assumed a half-way position between prebiotic organic synthesis and full-fledged cellular (biotic) life.

Recent enhancements and current research in the GeneCards (GC) (http://bioinfo.weizmann.ac.il/cards/) project are described, including the addition of gene expression profiles and integrated gene locations. Also highlighted are the contributions of specialized associated human gene-centric databases developed at the Weizmann Institute. These include the Unified Database (UDB) (http://bioinfo.weizmann.ac.il/udb) for human genome mapping, the human Chromosome 21 database at the Weizmann Insti-tute (CroW 21) (http://bioinfo.weizmann.ac.il/crow21), and the Human Olfactory Receptor Data Explora-torium (HORDE) (http://bioinfo.weizmann.ac.il/HORDE). The synergistic relationships amongst these efforts have positively impacted the quality, quantity and usefulness of the GeneCards gene compendium.

The goal of pharmacogenetics is to identify "genetic fingerprints" that may predict a patient's response to pharmaceutical treatment. The use of pharmacogenetics replaces the trial-and-error strategy, which governs much of our clinical decision-making regarding treatment allocation in current medical practice, with individually tailored therapy. We review a pharmacogenetic research model, which implements high-throughput single nucleotide polymorphism technology to establish the correlation between drug-responsiveness and genetic polymorphisms of Copaxone(R)-treated multiple sclerosis patients. Implementation of similar pharmacogenetic approaches may promote the development of personalized medicine in multiple sclerosis as well as in other diseases. (C) 2002 Prous Science. All rights reserved.

We introduce and study an artificial neural network inspired by the probabilistic receptor affinity distribution model of olfaction. Our system consists of N sensory neurons whose outputs converge on a single processing linear threshold element. The system's aim is to model discrimination of a single target odorant from a large number p of background odorants within a range of odorant concentrations. We show that this is possible provided p does not exceed a critical value p(c) and calculate the critical capacity a(c) = p(c)/N. The critical capacity depends on the range of concentrations in which the discrimination is to be accomplished. If the olfactory bulb may be thought of as a collection of such processing elements, each responsible for the discrimination of a single odorant, our study provides a quantitative analysis of the potential computational properties of the olfactory bulb. The mathematical formulation of the problem we consider is one of determining the capacity for linear separability of continuous curves, embedded in a large-dimensional space. This is accomplished here by a numerical study, using a method that signals whether the discrimination task is realizable, together with a finite-size scaling analysis.

A chance encounter between members of a random repertoire and a molecular target is characteristic of different biological systems, including the immune and olfactory pathways as well as combinatorial libraries. In such systems, the affinity between the target and members of the repertoire is distributed with a probability function describing the propensity of. obtaining a particular affinity value. We have previously proposed a phenomenological receptor affinity distribution (RAD) formalism, which describes this probability function based on simple statistical considerations. In the present analysis, we use published data from diverse experimental systems, including phage display libraries, immunoglobulins and enzymes, to test the RAD model and to compare it to other affinity distribution formalisms. The RAD model is found to provide the best description for binding data for over eight orders of magnitude on the affinity scale, and to account for a relationship between repertoire size and the maximal obtainable affinity within different repertoires. This approach points to a potential universality of the rules that govern affinity distributions in biology. (C) 2002 Elsevier Science Ltd. All rights reserved.

We investigated the population differences in patterns of single nucleotide polymorphisms (SNPs) for a 400 kb olfactory receptor (OR) gene cluster on human chromosome 17p13.3. Samples were drawn from 35 individuals, of four different ethnogeographical origins: Pygmies, Bedouins, Yemenite Jews and Ashkenazi Jews. Of the 74 SNPs identified, two segregated between pseudogenized and intact ORs, while a third involved a change in a highly conserved motif proposed to mediate ligand-induced signal transduction. Linkage disequilibrium (LD) was computed based on phase inference across the cluster using Clark's haplotype subtraction algorithm. We also calculated LD directly from the genotypes using the expectation-maximization (EM) algorithm. Both methods yielded very similar results. Our analyses revealed substantial differences in nucleotide diversity, haplotype distribution and LD patterns among the different human populations. In particular, the two Jewish populations had low haplotype diversity and negligible decay of LD across the entire genomic region. Intriguingly, the three functional SNPs segregated at different frequencies in the different ethnogeographical groups, with the Pygmies having higher frequencies of the intact OR genes. Our data suggests that OR genes may have evolved to create different functional repertoires in distinct human populations.

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) Ca2+-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 Ca2+ 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 Ca2+ binding.

The concept of shape space, which has been successfully implemented in immunology, is used here to construct a model for the discrimination power of the olfactory system. Using reasonable assumptions on the behaviour of the biological system, we are able to estimate the number of distinct olfactory receptor types. Our estimated value of around 1000 receptor types is in good agreement with experimental data. (C) 2001 Society for Mathematical Biology.

The olfactory receptor (OR) subgenome harbors the largest known gene family in mammals, disposed in clusters on numerous chromosomes. One of the best characterized OR clusters, located at human chromosome 17p13.3, has previously been studied by us in human and in other primates, revealing a conserved set of 17 OR genes. Here, we report the identification of a syntenic OR cluster in the mouse and the partial DNA sequence of many of its OR genes. A probe for the mouse M5 gene, orthologous to one of the OR genes in the human cluster (OR17-25), was used to isolate six PAC clones, all mapping by in situ hybridization to mouse chromosome 11B3-11B5, a region of shared synteny with human chromosome 17p13.3. Thirteen mouse OR sequences amplified and sequenced from these PACs allowed us to construct a putative physical map of the OR gene cluster at the mouse Olfr1 locus. Several points of evidence, including a strong similarity in subfamily composition and at least four cases of gene orthology, suggest that the mouse Olfr1 and the human 17p13.3 clusters are orthologous. A detailed comparison of the OR sequences within the two clusters helps trace their independent evolutionary history in the two species. Two types of evolutionary scenarios are discerned: cases of "true orthologous genes" in which high sequence similarity suggests a shared conserved function, as opposed to instances in which orthologous genes may have undergone independent diversification in the realm of "free reign" repertoire expansion. (C) 2001 Academic Press.

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

Olfactory receptors (ORs) constitute the largest multigene family in multicellular organisms. Their evolutionary proliferation has been driven by the need to provide recognition capacity for millions of potential odorants with arbitrary chemical configurations. Human genome sequencing has provided a highly informative picture of die "olfactory subgenome", the repertoire of OR genes. We describe here an analysis of 224 human OR genes, a much larger number than hitherto systematically analyzed. These are derived by literature survey, data mining at 14 genomic clusters, and by an OR-targeted experimental sequencing strategy. The presented set contains at least 53% pseudogenes and is minimally divided into 11 gene families. One of these (no. 7) has undergone a particularly extensive expansion in primates. The analysis of this collection leads to insight into the origin of OR genes, suggesting a graded expansion through mammalian evolution. It also allows us to delineate a structural map of the respective proteins. A sequence database and analysis package is provided (http://bioinformatics.weizmann.ac.il/HORDE), which will be useful for analyzing human OR sequences genome-wide.

The genomic and cDNA structures were studied for eight human olfactory receptor (OR) genes within the chromosome 17p13.3 cluster. A common gene structure was revealed, which included an similar to1-kb intronless coding region terminated by a signal for polyadenylation and a variable number of upstream noncoding exons. The latter were found to be alternatively spliced, giving rise to different isoforms of OR mRNA. While the initial exons mostly agreed with previous computer predictions and were conserved within OR subfamilies, other upstream exons were novel and idiosyncratic. In some cases, repetitive sequences were involved in the generation of splice sites and putative transcription control elements. Such gene structure is consistent with early repertoire enhancement by retrogene generation, which was likely followed by extensive genomic duplication. Each OR gene had a unique signature of transcription factor elements, consistent with a combinatorial expression control mechanism. (C) 2000 Academic Press.

Substantial efforts are focused on identifying single-nucleotide polymorphisms (SNPs) throughout the human genome, particularly in coding regions (cSNPs), for both linkage disequilibrium and association studies(1,2). Less attention, however, has been directed to the clarification of evolutionary processes that are responsible for the variability in nucleotide diversity among different regions of the genome(3). We report here the population sequence diversity of genomic segments within a 450-kb cluster(4,5) of olfactory receptor (OR) genes(6,7) on human chromosome 17. We found a dichotomy in the pattern of nucleotide diversity between OR pseudogenes and introns on the one hand and the closely interspersed intact genes on the other. We suggest that weak positive selection is responsible for the observed patterns of genetic variation. This is inferred from a lower ratio of polymorphism to divergence in genes compared with pseudogenes or introns, high non-synonymous substitution rates in OR genes. and a small but significant overall reduction in variability in the entire OR gene cluster compared with other genomic regions. The dichotomy among functionally different segments within a short genomic distance requires high recombination rates within this OR cluster. Our work demonstrates the impact of weak positive selection on human nucleotide diversity, and has implications for the evolution of the olfactory repertoire.

The GESALT Workbench is a WWW-based tool for genomic sequence analysis, comparisn and annotation, with strong emphasis on visualization. GESALT integrates graphically the output of diverse sequence analysis algorithms producing an information-rich, interactive genomic map. Availability: The GESALT Workbench, as well as a more detailed description, are available at http://bioinfo.weizmann.ac.il/GESALT/. ContactL Gustavo @ bioinfo.weizmann.ac.il, Doron.Lancet @ weizmann.ac.il.

The primordial emergence of biopolymers, agents of the genetic machinery in modern cells, is not less enigmatic than the emergence of the genetic code itself. Here we discuss how potential early replicating protocellular systems based on a rudimentary form of inheritance, a "compositional genome", could evolve towards the emergence of "alphabetic" polymers, predating the genetic code. A computer simulated evolutionary process based on our previously proposed kinetic model may help understand the appearance of chemical combinatorics through early natural selection.

The olfactory receptor (OR) subgenome harbors the largest known gene family in mammals, disposed in clusters on numerous chromosomes. We have carried out a comparative evolutionary analysis of the best characterized genomic OR gene cluster, on human chromosome 17p13. Fifteen orthologs from chimpanzee (localized to chromosome 19p15), as well as key OR counterparts from other primates, have been identified and sequenced. Comparison among orthologs and paralogs revealed a multiplicity of gene conversion events, which occurred exclusively within OR subfamilies. These appear to lead to segment shuffling in the odorant binding site, an evolutionary process reminiscent of somatic combinatorial diversification in the immune system We also demonstrate that the functional mammalian OR repertoire has undergone a rapid decline in the past 10 million years: while for the common ancestor of all great apes an intact OR cluster is inferred, in present-day humans and great apes the cluster includes nearly 40% pseudogenes. (C) 1999 Academic Press.

Olfactory receptors (OR) are encoded by a large multigene family including hundreds of members dispersed throughout the human genome. Cloning and mapping studies have determined that a large proportion of the olfactory receptor genes are located on human chromosomes 6, 11, and 17, as well as distributed on other chromosomes. In this paper, we describe and characterize the organization of olfactory receptor genes on human chromosome 11 by using degenerate PCR-based probes to screen chromosome 11-specific and whole genome clone libraries for members of the OR gene family. OR genes were identified by DNA sequencing and then localized to regions of chromosome 11. Physical maps of several gene clusters were constructed to determine the chromosomal relationships between various members of the family. This work identified 25 new OR genes located on chromosome 11 in at least seven distinct regions. Three of these regions contain gene clusters that include additional members of this gene family not yet identified by sequencing. Phylogenetic analysis of the newly described OR genes suggests a mechanism for the generation of genetic diversity. (C) 1998 Academic Press.

A description of the emergence of life should delineate a chemically rigorous gradual transition from random collections of simple organic molecules to spatially confined assemblies displaying rudimentary self-reproduction capacity. It has been suggested that large sets of mutually catalytic molecules, and not self-replicating information-carrying biopolymers, could have been the precursors of life. We present here a stochastic model in which the mutually catalytic molecules are spontaneously aggregating amphiphiles. When such amphiphiles exert on each other random catalytic effects, biased molecular compositions emerge, that are endowed with replication-like properties. This approach may have important consequences to the understanding of very early chemical evolution. It could also guide a search for extraterrestrial forms of very primitive life.

This review cc represents an update of the nomenclature system for the UDP glucuronosyltransferase gene superfamily, which is based on divergent evolution. Since the previous review in 1991, sequences of many related UDP glycosyltransferases from lower organisms have appeared in the database, which expand our database considerably. At latest count, in animals, yeast, plants and bacteria there are 110 distinct cDNAs/genes whose protein products all contain a characteristic 'signature sequence' and, thus, are regarded as members of the same superfamily, Comparison of a relatedness tree of proteins leads to the definition of 33 families, It should be emphasized that at least six cloned UDP-GlcNAc N-acetylglucosaminyltransferases are not sufficiently homologous to be included as members of this superfamily and may represent an example of convergent evolution, For naming each gene, it is recommended that the root symbol UGT for human (Ugt for mouse and Drosophila), denoting 'UDP glycosyltransferase,' be followed by an Arabic number representing the family, a letter designating the subfamily, and an Arabic numeral denoting the individual gene within the family or subfamily, e.g. 'human UGT2B4' and 'mouse Ugt2b5'. We recommend the name 'UDP glycosyltransferase' because many of the proteins do not preferentially use UDP glucuronic acid, or their nucleotide sugar preference is unknown. Whereas the gene is italicized, the corresponding cDNA, transcript, protein and enzyme activity should be written with upper-case letters and without italics, e.g. 'human or mouse UGT1A1. 'The UGT1 gene (spanning > 500 kb) contains at least 12 promoters/first exons, which call be spliced and joined with common exons 2 through 5, leading to different N-terminal halves but identical C-terminal halves of the gene products; in this scheme each first exon is regarded as a distinct gene (e.g. UGT1A1, UGT1A2, ... UGT1A12). When an orthologous gene between species cannot be identified with certaint

The superfamily of olfactory receptor genes, whose products are thought to be activated by odorant ligands, is critical for odor recognition. Two-olfactory receptors, olp4 from rat and OR17-4 from human, were overexpressed in Sf9 insect cells, The presence of the proteins in cell membranes was monitored by immunoblotting with peptide-specific polyclonal antibodies directed against the C-terminal sequences of these receptors and with a mAb against an N-terminal octapeptide epitope tag. A DNA sequence that codes for a His(6) tag, which binds tightly to a Ni2+-chelale-affinity column, was incorporated into the N-termini of both genes. The expressed olfactory receptors were found mainly in the cell-membrane fraction. The proteins were difficult to solubilize by many detergents and only lysophosphatidylcholine was found to be both suitable for efficient solubilization of the overexpressed olfactory receptors and compatible with the purification system used. After solubilization, the olfactory receptors were purified to near homogeneity by affinity chromatography on nickel nitrilotriacetic acid resin and by cation-exchange chromatography. Electrophoresis of the purified proteins and visualization with Coomassie Blue staining or by immunoblotting with specific antibodies, revealed bands of 32, 69 and 94 kDa, which were identified as the monomeric, dimeric and trimeric forms of the receptor proteins. The oligomeric forms were resistant to reduction and alkylation, and are therefore thought to be held together by non-covalent hydrophobic interactions that are resistant to SDS. This finding is similar to previous observations for other guanine-nucleotide-binding-regulator-protein-coupled receptors. Reconstitution in phospholipid vesicles showed that the purified olfactory receptors insert specifically into the lipid bilayer. This provides a means to study functional reconstitution with putative transduction components such as olfactory guanine-nucleotide-binding-regulatory

A rat olfactory epithelium cDNA library was screened for olfactory receptor clones. One of the positively hybridizing cDNA clones was sequenced and found to encode a new member of the olfactory receptor superfamily. This cDNA, termed olp4, was used as a model of olfactory receptor for expression, both in vitro and in vivo. Expression of olp4, as well as of another previously cloned olfactory receptor (F5), was monitored by immunoprecipitation with a monoclonal antibody directed against a Flag peptide epitope tag, inserted at the N-terminus of the open reading frame, and a specific polyclonal antibody against a C-terminal peptide of olp4. Translation in vitro, followed by immunoprecipitation, showed a major olp4-specific band of 27-29 kDa. The olp4 and F5 polypeptides were found to be inserted into microsomal membranes as expected for integral membrane proteins. Expression in vivo of Flag-olp4 in Sf9 insect cells, using the baculovirus expression system, showed a specific polypeptide of the same size as the in vitro species, with an additional band of 34 kDa, which is most likely a glycosylated form. Fluorescence cytometry and immunohistochemical assays demonstrated the localization of the Flag-olp4 product on the cell surface of the infected host Sf9 cells, with the N-terminus and C-terminus in the proper orientation. Affinity chromatography was used for the partial purification of the olp4 polypeptide from infected Sf9 cells. The identification and purification of this expressed olfactory receptor polypeptide could open the way for further characterization and functional studies of the olfactory receptor superfamily members.

The olfactory epithelium is exposed to a variety of xenobiotic chemicals, including odorants and airborne toxic compounds. Recently, two novel, highly abundant, olfactory-specific biotransformation enzymes have been identified: cytochrome P-450olf1 and olfactory UDP-glucuronosyltransferase (UGT(olf)). The latter is a phase II biotransformation enzyme which catalyses the glucuronidation of alcohols, thiols, amines and carboxylic acids. Such covalent modification, which markedly affects lipid solubility and agonist potency, may be particularly important in the rapid termination of odorant signals. We report here the identification and characterization of a second olfactory phase II biotransformation enzyme, a glutathione S-transferase (GST). The olfactory epithelial cytosol shows the highest GST activity among the extrahepatic tissues examined. Significantly, olfactory epithelium had an activity 4-7 times higher than in other airway tissues, suggesting a role for this enzyme in chemoreception. The olfactory GST has been affinity-purified to homogeneity, and shown by h.p.l.c. and N-terminal amino acid sequencing to constitute mainly the Yb1 and Yb2 subunits, different from most other tissues that have mixtures of more enzyme classes. The identity of the olfactory enzymes was confirmed by PCR cloning and restriction enzyme analysis. Most importantly, the olfactory GSTs were found to catalyse glutathione.conjugation of several odorant classes, including many unsaturated aldehydes and ketones, as well as epoxides. Together with UGT(olf), olfactory GST provides the necessary broad coverage of covalent modification capacity, which may be crucial for the acuity of the olfactory process.

A generalized phenomenological model is presented for stereospecific recognition between biological receptors and their ligands. We ask what is the distribution of binding constants PSI(K) between an arbitrary ligand and members of a large receptor repertoire, such as immunoglobulins or olfactory receptors. For binding surfaces with B potential subsite and S different types of subsite configurations, the number of successful elementary interactions obeys a binomial distribution. The discrete probability function PSI(K) is then derived with assumptions on alpha, the free energy contribution per elementary interaction. The functional form of PSI(K) may be universal, although the parameter values could vary for different ligand types. An estimate of the parameter values of PSI(K) for iodovanillin, an analog of odorants and immunological haptens, is obtained by equilibrium dialysis experiments with nonimmune antibodies. Based on a simple relationship, predicted by the model, between the size of a receptor repertoire and its average maximal affinity toward an arbitrary ligand, the size of the olfactory receptor repertoire (N(olf)) is calculated as 300-1000, in very good agreement with recent molecular biological studies. A very similar estimate, N(olf) = 500, is independently derived by relating a theoretical distribution of maxima for PSI(K) with published human olfactory threshold variations. The present model also has implications to the question of olfactory coding and to the analysis of specific anosmias, genetic deficits in perceiving particular odorants. More generally, the proposed model provides a better understanding of ligand specificity in biological receptors and could help in understanding their evolution.

A nomenclature system for the UDP glucuronosyltransferase superfamily is proposed, based on divergent evolution of the genes. A total of 26 distinct cDNAs in five mammalian species have been sequenced to date. Comparison of the deduced amino acid sequences leads to the definition of two families and a total of three subfamilies. For naming each gene, we propose that the root symbol UGT for human (Ugt for mouse), representing "UDP glucuronosyltransferase," be followed by an Arabic number denoting the family, a letter designating the subfamily, and an Arabic numeral representing the individual gene within the family or subfamily (hyphen before the Arabic number for mouse), e.g., human UGT2B1 and murine Ugt2b-1. Whereas the gene and cDNA should be italicized, the corresponding transcript, protein, and enzyme activity should not be written with lowercase letters or in italics, e.g., human or murine UGT2B1. Recent experimental evidence suggests that several exons of the UGT1 gene might be shared, indicating that distinct UGT1 transcripts and proteins may arise via alternative splicing; the gene and gene product of alternative splicing will be designated with an asterisk, e.g., UGT1*6 and UGT1*6, respectively. When an orthologous gene between species cannot be identified with certainty, as occurs in the UGT2B subfamily, we recommend sequential naming of the genes chronologically as they become characterized. We suggest that the human nomenclature system be used for species other than the mouse. We anticipate that this UGT gene nomenclature system will require updating on a regular basis.

THE onset of olfactory transduction has been extensively studied 1-7, but considerably less is known about the molecular basis of olfactory signal termination 6,8,9. It has been suggested that the highly active cytochrome P450 monooxygenases of olfactory neuroepithelium 10-12 are termination enzymes 5,8,11,12, a notion supported by the identification and molecular cloning of olfactory-specific cytochrome P450s (refs. 13-16). But as reactions catalysed by cytochrome P450 (refs 17, 18) often do not significantly alter volatility, lipophilicity or odour properties 9,11, cytochrome P450 may not be solely responsible for olfactory signal termination. In liver and other tissues, drug hydroxylation by cytochrome P450 is frequently followed by phase II biotransformation, for example by UDP glucuronosyl transferase (UGT), resulting in a major change of solubility and chemical properties 19. We report here the molecular cloning and expression of an olfactory-specific UGT. The olfactory enzyme, but not the one in liver microsomes, shows preference for odorants over standard UGT substrates. Furthermore, glucuronic acid conjugation abolishes the ability of odorants 1,20 to stimulate olfactory adenylyl cyclase. This, together with the known broad spectrum of drug-detoxification enzymes 17,19, supports a role for olfactory UGT in terminating diverse odorant signals.