Spatiotemporal Proteomic Analysis of Stress Granule Disassembly Using APEX Reveals Regulation by SUMOylation and Links to ALS Pathogenesis.

Stress granules (SGs) are cytoplasmic assemblies of proteins and non-translating mRNAs. Whereas much has been learned about SG formation, a major gap remains in understanding the compositional changes SGs undergo during normal disassembly and under disease conditions. Here, we address this gap by proteomic dissection of the SG temporal disassembly sequence using multi-bait APEX proximity proteomics. We discover 109 novel SG proteins and characterize distinct SG substructures. We reveal dozens of disassembly-engaged proteins (DEPs), some of which play functional roles in SG disassembly, including small ubiquitin-like modifier (SUMO) conjugating enzymes. We further demonstrate that SUMOylation regulates SG disassembly and SG formation. Parallel proteomics with amyotrophic lateral sclerosis (ALS)-associated C9ORF72 dipeptides uncovered attenuated DEP recruitment during SG disassembly and impaired SUMOylation. Accordingly, SUMO activity ameliorated C9ORF72-ALS-related neurodegeneration in Drosophila. By dissecting the SG spatiotemporal proteomic landscape, we provide an in-depth resource for future work on SG function and reveal basic and disease-relevant mechanisms of SG disassembly.

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

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

DAP5 drives translation of specific mRNA targets with upstream ORFs in human embryonic stem cells.

Death associated protein 5 (DAP5/eIF4G2/NAT1) is a member of the eIF4G translation initiation factors that has been shown to mediate noncanonical and/or cap-independent translation. It is essential for embryonic development and for differentiation of embryonic stem cells (ESCs), specifically its ability to drive translation of specific target mRNAs. In order to expand the repertoire of DAP5 target mRNAs, we compared ribosome profiles in control and DAP5 knockdown (KD) human ESCs (hESCs) to identify mRNAs with decreased ribosomal occupancy upon DAP5 silencing. A cohort of 68 genes showed decreased translation efficiency in DAP5 KD cells. Mass spectrometry confirmed decreased protein abundance of a significant portion of these targets. Among these was KMT2D, a histone methylase previously shown to be essential for ESC differentiation and embryonic development. We found that nearly half of the cohort of DAP5 target mRNAs displaying reduced translation efficiency of their main coding sequences upon DAP5 KD contained upstream open reading frames (uORFs) that are actively translated independently of DAP5. This is consistent with previously suggested mechanisms by which DAP5 mediates leaky scanning through uORFs and/or reinitiation at the main coding sequence. Crosslinking protein-RNA immunoprecipitation experiments indicated that a significant subset of DAP5 mRNA targets bound DAP5, indicating that direct binding between DAP5 protein and its target mRNAs is a frequent but not absolute requirement for DAP5-dependent translation of the main coding sequence. Thus, we have extended DAP5's function in translation of specific mRNAs in hESCs by a mechanism allowing translation of the main coding sequence following upstream translation of short ORFs.

rG4detector, a novel RNA G-quadruplex predictor, uncovers their impact on stress granule formation.

RNA G-quadruplexes (rG4s) are RNA secondary structures, which are formed by guanine-rich sequences and have important cellular functions. Existing computational tools for rG4 prediction rely on specific sequence features and/or were trained on small datasets, without considering rG4 stability information, and are therefore sub-optimal. Here, we developed rG4detector, a convolutional neural network to identify potential rG4s in transcriptomics data. rG4detector outperforms existing methods in both predicting rG4 stability and in detecting rG4-forming sequences. To demonstrate the biological-relevance of rG4detector, we employed it to study RNAs that are bound by the RNA-binding protein G3BP1. G3BP1 is central to the induction of stress granules (SGs), which are cytoplasmic biomolecular condensates that form in response to a variety of cellular stresses. Unexpectedly, rG4detector revealed a dynamic enrichment of rG4s bound by G3BP1 in response to cellular stress. In addition, we experimentally characterized G3BP1 cross-talk with rG4s, demonstrating that G3BP1 is a bona fide rG4-binding protein and that endogenous rG4s are enriched within SGs. Furthermore, we found that reduced rG4 availability impairs SG formation. Hence, we conclude that rG4s play a direct role in SG biology via their interactions with RNA-binding proteins and that rG4detector is a novel useful tool for rG4 transcriptomics data analyses.

A secretion-enhancing cis regulatory targeting element (SECReTE) involved in mRNA localization and protein synthesis.

The localization of mRNAs encoding secreted/membrane proteins (mSMPs) to the endoplasmic reticulum (ER) likely facilitates the co-translational translocation of secreted proteins. However, studies have shown that mSMP recruitment to the ER in eukaryotes can occur in a manner that is independent of the ribosome, translational control, and the signal recognition particle, although the mechanism remains largely unknown. Here, we identify a cis-acting RNA sequence motif that enhances mSMP localization to the ER and appears to increase mRNA stability, and both the synthesis and secretion of secretome proteins. Termed SECReTE, for secretion-enhancing cis regulatory targeting element, this motif is enriched in mRNAs encoding secretome proteins translated on the ER in eukaryotes and on the inner membrane of prokaryotes. SECReTE consists of ≥10 nucleotide triplet repeats enriched with pyrimidine (C/U) every third base (i.e. NNY, where N = any nucleotide, Y = pyrimidine) and can be present in the untranslated as well as the coding regions of the mRNA. Synonymous mutations that elevate the SECReTE count in a given mRNA (e.g. SUC2, HSP150, and CCW12) lead to an increase in protein secretion in yeast, while a reduction in count led to less secretion and physiological defects. Moreover, the addition of SECReTE to the 3'UTR of an mRNA for an exogenously expressed protein (e.g. GFP) led to its increased secretion from yeast cells. Thus, SECReTE constitutes a novel RNA motif that facilitates ER-localized mRNA translation and protein secretion.

A unified nomenclature for vertebrate olfactory receptors.

BACKGROUND: Olfactory receptors (ORs) are G protein-coupled receptors with a crucial role in odor detection. A typical mammalian genome harbors ~ 1000 OR genes and pseudogenes; however, different gene duplication/deletion events have occurred in each species, resulting in complex orthology relationships. While the human OR nomenclature is widely accepted and based on phylogenetic classification into 18 families and further into subfamilies, for other mammals different and multiple nomenclature systems are currently in use, thus concealing important evolutionary and functional insights. RESULTS: Here, we describe the Mutual Maximum Similarity (MMS) algorithm, a systematic classifier for assigning a human-centric nomenclature to any OR gene based on inter-species hierarchical pairwise similarities. MMS was applied to the OR repertoires of seven mammals and zebrafish. Altogether, we assigned symbols to 10,249 ORs. This nomenclature is supported by both phylogenetic and synteny analyses. The availability of a unified nomenclature provides a framework for diverse studies, where textual symbol comparison allows immediate identification of potential ortholog groups as well as species-specific expansions/deletions; for example, Or52e5 and Or52e5b represent a rat-specific duplication of OR52E5. Another example is the complete absence of OR subfamily OR6Z among primate OR symbols. In other mammals, OR6Z members are located in one genomic cluster, suggesting a large deletion in the great ape lineage. An additional 14 mammalian OR subfamilies are missing from the primate genomes. While in chimpanzee 87% of the symbols were identical to human symbols, this number decreased to ~ 50% in dog and cow and to ~ 30% in rodents, reflecting the adaptive changes of the OR gene superfamily across diverse ecological niches. Application of the proposed nomenclature to zebrafish revealed similarity to mammalian ORs that could not be detected from the current zebrafish olfactory receptor gene nomenclature. CONCLUSIONS: We have consolidated a unified standard nomenclature system for the vertebrate OR superfamily. The new nomenclature system will be applied to cow, horse, dog and chimpanzee by the Vertebrate Gene Nomenclature Committee and its implementation is currently under consideration by other relevant species-specific nomenclature committees.

An active ß-lactamase is a part of an orchestrated cell wall stress resistance network of Bacillus subtilis and related rhizosphere species.

A hallmark of the Gram-positive bacteria, such as the soil-dwelling bacterium Bacillus subtilis, is their cell wall. Here, we report that d-leucine and flavomycin, biofilm inhibitors targeting the cell wall, activate the ß-lactamase PenP. This ß-lactamase contributes to ampicillin resistance in B. subtilis under all conditions tested. In contrast, both Spo0A, a master regulator of nutritional stress, and the general cell wall stress response, differentially contribute to ß-lactam resistance under different conditions. To test whether ß-lactam resistance and ß-lactamase genes are widespread in other Bacilli, we isolated Bacillus species from undisturbed soils, and found that their genomes can encode up to five ß-lactamases with differentiated activity spectra. Surprisingly, the activity of environmental ß-lactamases and PenP, as well as the general stress response, resulted in a similarly reduced lag phase of the culture in the presence of ß-lactam antibiotics, with little or no impact on the logarithmic growth rate. The length of the lag phase may determine the outcome of the competition between ß-lactams and ß-lactamases producers. Overall, our work suggests that antibiotic resistance genes in B. subtilis and related species are ancient and widespread, and could be selected by interspecies competition in undisturbed soils.

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

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

The Disordered Landscape of the 20S Proteasome Substrates Reveals Tight Association with Phase Separated Granules.

Proteasomal degradation is the main route of regulated proteostasis. The 20S proteasome is the core particle (CP) responsible for the catalytic activity of all proteasome complexes. Structural constraints mean that only unfolded, extended polypeptide chains may enter the catalytic core of the 20S proteasome. It has been previously shown that the 20S CP is active in degradation of certain intrinsically disordered proteins (IDP) lacking structural constrains. Here, a comprehensive analysis of the 20S CP substrates in vitro is conducted. It is revealed that the 20S CP substrates are highly disordered. However, not all the IDPs are 20S CP substrates. The group of the IDPs that are 20S CP substrates, termed 20S-IDPome are characterized by having significantly more protein binding partners, more posttranslational modification sites, and are highly enriched for RNA binding proteins. The vast majority of them are involved in splicing, mRNA processing, and translation. Remarkably, it is found that low complexity proteins with prion-like domain (PrLD), which interact with GR or PR di-peptide repeats, are the most preferential 20S CP substrates. The finding suggests roles of the 20S CP in gene transcription and formation of phase-separated granules.

Characterization of proteome dynamics during growth in oleate reveals a new peroxisome-targeting receptor.

To optimally perform the diversity of metabolic functions that occur within peroxisomes, cells must dynamically regulate peroxisome size, number and content in response to the cell state and the environment. Except for transcriptional regulation little is known about the mechanisms used to perform this complicated feat. Focusing on the yeast Saccharomyces cerevisiae, we used complementary high-content screens to follow changes in localization of most proteins during growth in oleate. We found extensive changes in cellular architecture and identified several proteins that colocalized with peroxisomes that had not previously been considered peroxisomal proteins. One of the newly identified peroxisomal proteins, Ymr018w, is a protein with an unknown function that is similar to the yeast and human peroxisomal targeting receptor Pex5. We demonstrate that Ymr018w is a new peroxisomal-targeting receptor that targets a subset of matrix proteins to peroxisomes. We, therefore, renamed Ymr018w, Pex9, and suggest that Pex9 is a condition-specific targeting receptor that enables the dynamic rewiring of peroxisomes in response to metabolic needs. Moreover, we suggest that Pex5-like receptors might also exist in vertebrates.

The human olfactory transcriptome.

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

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

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

Regulation of receptor-type protein tyrosine phosphatases by their C-terminal tail domains.

Protein tyrosine phosphatases (PTPs) perform specific functions in vivo, despite being vastly outnumbered by their substrates. Because of this and due to the central roles PTPs play in regulating cellular function, PTP activity is regulated by a large variety of molecular mechanisms. We review evidence that indicates that the divergent C-terminal tail sequences (C-terminal domains, CTDs) of receptor-type PTPs (RPTPs) help regulate RPTP function by controlling intermolecular associations in a way that is itself subject to physiological regulation. We propose that the CTD of each RPTP defines an 'interaction code' that helps determine molecules it will interact with under various physiological conditions, thus helping to regulate and diversify PTP function.

The genome of a songbird.

The zebra finch is an important model organism in several fields with unique relevance to human neuroscience. Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chicken-the only bird with a sequenced genome until now. Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes. We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat-based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour.

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

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

Genome analysis of the platypus reveals unique signatures of evolution.

We present a draft genome sequence of the platypus, Ornithorhynchus anatinus. This monotreme exhibits a fascinating combination of reptilian and mammalian characters. For example, platypuses have a coat of fur adapted to an aquatic lifestyle; platypus females lactate, yet lay eggs; and males are equipped with venom similar to that of reptiles. Analysis of the first monotreme genome aligned these features with genetic innovations. We find that reptile and platypus venom proteins have been co-opted independently from the same gene families; milk protein genes are conserved despite platypuses laying eggs; and immune gene family expansions are directly related to platypus biology. Expansions of protein, non-protein-coding RNA and microRNA families, as well as repeat elements, are identified. Sequencing of this genome now provides a valuable resource for deep mammalian comparative analyses, as well as for monotreme biology and conservation.

Multigenerational inheritance drives symbiotic interactions of the bacterium Bacillus subtilis with its plant host.

Bacillus subtilis is a beneficial bacterium that supports plant growth and protects plants from bacterial, fungal, and viral infections. Using a simplified system of B. subtilis and Arabidopsis thaliana interactions, we studied the fitness and transcriptome of bacteria detached from the root over generations of growth in LB medium. We found that bacteria previously associated with the root or exposed to its secretions had greater stress tolerance and were more competitive in root colonization than bacteria not previously exposed to the root. Furthermore, our transcriptome results provide evidence that plant secretions induce a microbial stress response and fundamentally alter signaling by the cyclic nucleotide c-di-AMP, a signature maintained by their descendants. The changes in cellular physiology due to exposure to plant exudates were multigenerational, as they allowed not only the bacterial cells that colonized a new plant but also their descendants to have an advance over naive competitors of the same species, while the overall plasticity of gene expression and rapid adaptation were maintained. These changes were hereditary but not permanent. Our work demonstrates a bacterial memory manifested by multigenerational reversible adaptation to plant hosts in the form of activation of the stressosome, which confers an advantage to symbiotic bacteria during competition.

PRDM16 co-operates with LHX2 to shape the human brain.

PRDM16 is a dynamic transcriptional regulator of various stem cell niches, including adipocytic, hematopoietic, cardiac progenitors, and neural stem cells. PRDM16 has been suggested to contribute to 1p36 deletion syndrome, one of the most prevalent subtelomeric microdeletion syndromes. We report a patient with a de novo nonsense mutation in the PRDM16 coding sequence, accompanied by lissencephaly and microcephaly features. Human stem cells were genetically modified to mimic this mutation, generating cortical organoids that exhibited altered cell cycle dynamics. RNA sequencing of cortical organoids at day 32 unveiled changes in cell adhesion and WNT-signaling pathways. ChIP-seq of PRDM16 identified binding sites in postmortem human fetal cortex, indicating the conservation of PRDM16 binding to developmental genes in mice and humans, potentially at enhancer sites. A shared motif between PRDM16 and LHX2 was identified and further examined through comparison with LHX2 ChIP-seq data from mice. These results suggested a collaborative partnership between PRDM16 and LHX2 in regulating a common set of genes and pathways in cortical radial glia cells, possibly via their synergistic involvement in cortical development.

Loss of EIF4G2 mediates aggressiveness in distinct human endometrial cancer subpopulations with poor survival outcome in patients.

The non-canonical translation initiation factor EIF4G2 plays essential roles in cellular stress responses via translation of selective mRNA cohorts. Currently there is limited and conflicting information regarding its involvement in cancer development and progression. Here we assessed its role in endometrial cancer (EC), in a cohort of 280 EC patients across different types, grades, and stages, and found that low EIF4G2 expression highly correlated with poor overall- and recurrence-free survival in Grade 2 EC patients, monitored over a period of up to 12 years. To establish a causative connection between low EIF4G2 expression and cancer progression, we stably knocked-down EIF4G2 in two human EC cell lines in parallel. EIF4G2 depletion resulted in increased resistance to conventional therapies and increased the prevalence of molecular markers for aggressive cell subsets, altering their transcriptional and proteomic landscapes. Prominent among the proteins with decreased abundance were Kinesin-1 motor proteins, KIF5B and KLC1, 2, 3. Multiplexed imaging of the EC patient tumor cohort showed a correlation between decreased expression of the kinesin proteins, and poor survival in patients with tumors of certain grades and stages. These findings reveal potential novel biomarkers for Grade 2 EC with ramifications for patient stratification and therapeutic interventions.

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

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