Expanding and Enriching the LncRNA Gene-Disease Landscape Using the GeneCaRNA Database.

The GeneCaRNA human gene database is a member of the GeneCards Suite. It presents ~280,000 human non-coding RNA genes, identified algorithmically from ~690,000 RNAcentral transcripts. This expands by ~tenfold the ncRNA gene count relative to other sources. GeneCaRNA thus contains ~120,000 long non-coding RNAs (LncRNAs, >200 bases long), including ~100,000 novel genes. The latter have sparse functional information, a vast terra incognita for future research. LncRNA genes are uniformly represented on all nuclear chromosomes, with 10 genes on mitochondrial DNA. Data obtained from MalaCards, another GeneCards Suite member, finds 1547 genes associated with 1 to 50 diseases. About 15% of the associations portray experimental evidence, with cancers tending to be multigenic. Preliminary text mining within GeneCaRNA discovers interactions of lncRNA transcripts with target gene products, with 25% being ncRNAs and 75% proteins. GeneCaRNA has a biological pathways section, which at present shows 131 pathways for 38 lncRNA genes, a basis for future expansion. Finally, our GeneHancer database provides regulatory elements for ~110,000 lncRNA genes, offering pointers for co-regulated genes and genetic linkages from enhancers to diseases. We anticipate that the broad vista provided by GeneCaRNA will serve as an essential guide for further lncRNA research in disease decipherment.

Genome-wide association study identifies 16 genomic regions associated with circulating cytokines at birth.

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 < 5 x 10-9), of which four are novel. We show that the associated variants map to enhancer elements, suggesting a possible transcriptional effect of genomic variants on the cytokine levels. The identification of the genetic architecture underlying the basic levels of cytokines is likely to prompt studies investigating the relationship between cytokines and complex disease. Our results also suggest that the genetic architecture of cytokines is stable from neonatal to adult life.

Author Correction: Introducing PIONEER: a project to harness big data in prostate cancer research.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Introducing PIONEER: a project to harness big data in prostate cancer research.

Prostate Cancer Diagnosis and Treatment Enhancement Through the Power of Big Data in Europe (PIONEER) is a European network of excellence for big data in prostate cancer, consisting of 32 private and public stakeholders from 9 countries across Europe. Launched by the Innovative Medicines Initiative 2 and part of the Big Data for Better Outcomes Programme (BD4BO), the overarching goal of PIONEER is to provide high-quality evidence on prostate cancer management by unlocking the potential of big data. The project has identified critical evidence gaps in prostate cancer care, via a detailed prioritization exercise including all key stakeholders. By standardizing and integrating existing high-quality and multidisciplinary data sources from patients with prostate cancer across different stages of the disease, the resulting big data will be assembled into a single innovative data platform for research. Based on a unique set of methodologies, PIONEER aims to advance the field of prostate cancer care with a particular focus on improving prostate-cancer-related outcomes, health system efficiency by streamlining patient management, and the quality of health and social care delivered to all men with prostate cancer and their families worldwide.

Genome analysis and knowledge-driven variant interpretation with TGex.

BACKGROUND: The clinical genetics revolution ushers in great opportunities, accompanied by significant challenges. The fundamental mission in clinical genetics is to analyze genomes, and to identify the most relevant genetic variations underlying a patient's phenotypes and symptoms. The adoption of Whole Genome Sequencing requires novel capacities for interpretation of non-coding variants. RESULTS: We present TGex, the Translational Genomics expert, a novel genome variation analysis and interpretation platform, with remarkable exome analysis capacities and a pioneering approach of non-coding variants interpretation. TGex's main strength is combining state-of-the-art variant filtering with knowledge-driven analysis made possible by VarElect, our highly effective gene-phenotype interpretation tool. VarElect leverages the widely used GeneCards knowledgebase, which integrates information from > 150 automatically-mined data sources. Access to such a comprehensive data compendium also facilitates TGex's broad variant annotation, supporting evidence exploration, and decision making. TGex has an interactive, user-friendly, and easy adaptive interface, ACMG compliance, and an automated reporting system. Beyond comprehensive whole exome sequence capabilities, TGex encompasses innovative non-coding variants interpretation, towards the goal of maximal exploitation of whole genome sequence analyses in the clinical genetics practice. This is enabled by GeneCards' recently developed GeneHancer, a novel integrative and fully annotated database of human enhancers and promoters. Examining use-cases from a variety of TGex users world-wide, we demonstrate its high diagnostic yields (42% for single exome and 50% for trios in 1500 rare genetic disease cases) and critical actionable genetic findings. The platform's support for integration with EHR and LIMS through dedicated APIs facilitates automated retrieval of patient data for TGex's customizable reporting engine, establishing a rapid and cost-effective workflow for an entire range of clinical genetic testing, including rare disorders, cancer predisposition, tumor biopsies and health screening. CONCLUSIONS: TGex is an innovative tool for the annotation, analysis and prioritization of coding and non-coding genomic variants. It provides access to an extensive knowledgebase of genomic annotations, with intuitive and flexible configuration options, allows quick adaptation, and addresses various workflow requirements. It thus simplifies and accelerates variant interpretation in clinical genetics workflows, with remarkable diagnostic yield, as exemplified in the described use cases. TGex is available at http://tgex.genecards.org/.

Next-generation sequencing of patients with congenital anosmia.

We performed whole exome or genome sequencing in eight multiply affected families with ostensibly isolated congenital anosmia. Hypothesis-free analyses based on the assumption of fully penetrant recessive/dominant/X-linked models obtained no strong single candidate variant in any of these families. In total, these eight families showed 548 rare segregating variants that were predicted to be damaging, in 510 genes. Three Kallmann syndrome genes (FGFR1, SEMA3A, and CHD7) were identified. We performed permutation-based analysis to test for overall enrichment of these 510 genes carrying these 548 variants with genes mutated in Kallmann syndrome and with a control set of genes mutated in hypogonadotrophic hypogonadism without anosmia. The variants were found to be enriched for Kallmann syndrome genes (3 observed vs. 0.398 expected, p = 0.007), but not for the second set of genes. Among these three variants, two have been already reported in genes related to syndromic anosmia (FGFR1 (p.(R250W)), CHD7 (p.(L2806V))) and one was novel (SEMA3A (p.(T717I))). To replicate these findings, we performed targeted sequencing of 16 genes involved in Kallmann syndrome and hypogonadotrophic hypogonadism in 29 additional families, mostly singletons. This yielded an additional 6 variants in 5 Kallmann syndrome genes (PROKR2, SEMA3A, CHD7, PROK2, ANOS1), two of them already reported to cause Kallmann syndrome. In all, our study suggests involvement of 6 syndromic Kallmann genes in isolated anosmia. Further, we report a yet unreported appearance of di-genic inheritance in a family with congenital isolated anosmia. These results are consistent with a complex molecular basis of congenital anosmia.

Personal receptor repertoires: olfaction as a model.

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

In-silico human genomics with GeneCards.

Since 1998, the bioinformatics, systems biology, genomics and medical communities have enjoyed a synergistic relationship with the GeneCards database of human genes (http://www.genecards.org). This human gene compendium was created to help to introduce order into the increasing chaos of information flow. As a consequence of viewing details and deep links related to specific genes, users have often requested enhanced capabilities, such that, over time, GeneCards has blossomed into a suite of tools (including GeneDecks, GeneALaCart, GeneLoc, GeneNote and GeneAnnot) for a variety of analyses of both single human genes and sets thereof. In this paper, we focus on inhouse and external research activities which have been enabled, enhanced, complemented and, in some cases, motivated by GeneCards. In turn, such interactions have often inspired and propelled improvements in GeneCards. We describe here the evolution and architecture of this project, including examples of synergistic applications in diverse areas such as synthetic lethality in cancer, the annotation of genetic variations in disease, omics integration in a systems biology approach to kidney disease, and bioinformatics tools.

Systems medicine and integrated care to combat chronic noncommunicable diseases.

We propose an innovative, integrated, cost-effective health system to combat major non-communicable diseases (NCDs), including cardiovascular, chronic respiratory, metabolic, rheumatologic and neurologic disorders and cancers, which together are the predominant health problem of the 21st century. This proposed holistic strategy involves comprehensive patient-centered integrated care and multi-scale, multi-modal and multi-level systems approaches to tackle NCDs as a common group of diseases. Rather than studying each disease individually, it will take into account their intertwined gene-environment, socio-economic interactions and co-morbidities that lead to individual-specific complex phenotypes. It will implement a road map for predictive, preventive, personalized and participatory (P4) medicine based on a robust and extensive knowledge management infrastructure that contains individual patient information. It will be supported by strategic partnerships involving all stakeholders, including general practitioners associated with patient-centered care. This systems medicine strategy, which will take a holistic approach to disease, is designed to allow the results to be used globally, taking into account the needs and specificities of local economies and health systems.

Synthetic lethal hubs associated with vincristine resistant neuroblastoma.

Chemotherapy of cancer experiences a number of shortcomings including development of drug resistance. This fact also holds true for neuroblastoma utilizing chemotherapeutics as vincristine. We performed a comparative analysis of molecular and cellular mechanisms associated with vincristine resistance utilizing cell line as well as human tissue data. Differential gene expression analysis revealed molecular features, processes and pathways afflicted with drug resistance mechanisms in general, and specifically with vincristine significantly involving actin associated features. However, specific mode of resistance as well as underlying genotype of parental, vincristine sensitive cells apparently exhibited significant heterogeneity. No consensus profile for vincristine resistance could be derived, but resistance-associated changes on the level of individual neuroblastoma cell lines as well as individual patient profiles became clearly evident. Based on these prerequisites we utilized the concept of synthetic lethality aimed at identifying hub proteins which when inhibited promise to induce cell death due to a synthetic lethal interaction with down-regulated, chemoresistance associated features. Our screening procedure identified synthetic lethal hub proteins afflicted with actin associated processes holding synthetic lethal interactions to down-regulated features individually found in all chemoresistant cell lines tested, therefore promising an improved therapeutic window. Verification of such synthetic lethal hub candidates in human neuroblastoma tissue expression profiles indicated the feasibility of this screening approach for addressing vincristine resistance in neuroblastoma.

GeneCards Version 3: the human gene integrator.

GeneCards (www.genecards.org) is a comprehensive, authoritative compendium of annotative information about human genes, widely used for nearly 15 years. Its gene-centric content is automatically mined and integrated from over 80 digital sources, resulting in a web-based deep-linked card for each of >73,000 human gene entries, encompassing the following categories: protein coding, pseudogene, RNA gene, genetic locus, cluster and uncategorized. We now introduce GeneCards Version 3, featuring a speedy and sophisticated search engine and a revamped, technologically enabling infrastructure, catering to the expanding needs of biomedical researchers. A key focus is on gene-set analyses, which leverage GeneCards' unique wealth of combinatorial annotations. These include the GeneALaCart batch query facility, which tabulates user-selected annotations for multiple genes and GeneDecks, which identifies similar genes with shared annotations, and finds set-shared annotations by descriptor enrichment analysis. Such set-centric features address a host of applications, including microarray data analysis, cross-database annotation mapping and gene-disorder associations for drug targeting. We highlight the new Version 3 database architecture, its multi-faceted search engine, and its semi-automated quality assurance system. Data enhancements include an expanded visualization of gene expression patterns in normal and cancer tissues, an integrated alternative splicing pattern display, and augmented multi-source SNPs and pathways sections. GeneCards now provides direct links to gene-related research reagents such as antibodies, recombinant proteins, DNA clones and inhibitory RNAs and features gene-related drugs and compounds lists. We also portray the GeneCards Inferred Functionality Score annotation landscape tool for scoring a gene's functional information status. Finally, we delineate examples of applications and collaborations that have benefited from the GeneCards suite. Database URL: www.genecards.org.

Common peptides shed light on evolution of Olfactory Receptors.

BACKGROUND: Olfactory Receptors (ORs) form the largest multigene family in vertebrates. Their evolution and their expansion in the vertebrate genomes was the subject of many studies. In this paper we apply a motif-based approach to this problem in order to uncover evolutionary characteristics. RESULTS: We extract deterministic motifs from ORs belonging to ten species using the MEX (Motif Extraction) algorithm, thus defining Common Peptides (CPs) characteristic to ORs. We identify species-specific CPs and show that their relative abundance is high only in fish and frog, suggesting relevance to water-soluble odorants. We estimate the origins of CPs according to the tree of life and track the gains and losses of CPs through evolution. We identify major CP gain in tetrapods and major losses in reptiles. Although the number of human ORs is less than half of the number of ORs in other mammals, the fraction of lost CPs is only 11%.By examining the positions of CPs along the OR sequence, we find two regions that expanded only in tetrapods. Using CPs we are able to establish remote homology relations between ORs and non-OR GPCRs.Selecting CPs according to their evolutionary age, we bicluster ORs and CPs for each species. Clean biclustering emerges when using relatively novel CPs. Evolutionary age is used to track the history of CP acquisition in the collection of mammalian OR families within HORDE (Human Olfactory Receptor Data Explorer). CONCLUSION: The CP method provides a novel perspective that reveals interesting traits in the evolution of olfactory receptors. It is consistent with previous knowledge, and provides finer details. Using available phylogenetic trees, evolution can be rephrased in terms of CP origins.

Pharmacogenetics of glatiramer acetate therapy for multiple sclerosis reveals drug-response markers.

Genetic-based optimization of treatment prescription is becoming a central research focus in the management of chronic diseases, such as multiple sclerosis, which incur a prolonged drug-regimen adjustment. This study was aimed to identify genetic markers that can predict response to glatiramer acetate (Copaxone) immunotherapy for relapsing multiple sclerosis. For this purpose, we genotyped fractional cohorts of two glatiramer acetate clinical trials for HLA-DRB1*1501 and 61 single nucleotide polymorphisms within a total of 27 candidate genes. Statistical analyses included single nucleotide polymorphism-by-single nucleotide polymorphism and haplotype tests of drug-by-genotype effects in drug-treated versus placebo-treated groups. We report the detection of a statistically significant association between glatiramer acetate response and a single nucleotide polymorphism in a T-cell receptor beta (TRB@) variant replicated in the two independent cohorts (odds ratio=6.85). Findings in the Cathepsin S (CTSS) gene (P=0.049 corrected for all single nucleotide polymorphisms and definitions tested, odds ratio=11.59) in one of the cohorts indicate a possible association that needs to be further investigated. Additionally, we recorded nominally significant associations of response with five other genes, MBP, CD86, FAS, IL1R1 and IL12RB2, which are likely to be involved in glatiramer acetate's mode-of-action, both directly and indirectly. Each of these association signals in and of itself is consistent with the no-association null-hypothesis, but the number of detected associations is surprising vis-à-vis chance expectation. Moreover, the restriction of these associations to the glatiramer acetate-treated group, rather than the placebo group, clearly demonstrates drug-specific genetic effects. These findings provide additional progress toward development of pharmacogenetics-based personalized treatment for multiple sclerosis.

Haplotype structure and selection of the MDM2 oncogene in humans.

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.

Mutations in olfactory signal transduction genes are not a major cause of human congenital general anosmia.

Anosmia affects the western world population, mostly the elderly, reaching to 5% in subjects over the age of 45 years and strongly lowering their quality of life. A smaller minority (about 0.01%) is born without a sense of smell, afflicted with congenital general anosmia (CGA). No causative genes for human CGA have been identified yet, except for some syndromic cases such as Kallman syndrome. In mice, however, deletion of any of the 3 main olfactory transduction components (guanidine triphosphate binding protein, adenylyl cyclase, and the cyclic adenosine monophosphate-gated channel) causes profound reduction of physiological responses to odorants. In an attempt to identify human CGA-related mutations, we performed whole-genome linkage analysis in affected families, but no significant linkage signals were observed, probably due to the small size of families analyzed. We further carried out direct mutation screening in the 3 main olfactory transduction genes in 64 unrelated anosmic individuals. No potentially causative mutations were identified, indicating that transduction gene variations underlie human CGA rarely and that mutations in other genes have to be identified. The screened genes were found to be under purifying selection, suggesting that they play a crucial functional role not only in olfaction but also potentially in additional pathways.

Genotype phenotype correlations in Israeli colorectal cancer patients.

While genetic factors clearly play a key role in colorectal cancer (CRC) pathogenesis and in determining its phenotypic features, the precise genes that involved are largely unknown. To gain insight into these genes, consecutive Israeli CRC patients were genotyped using SNPs from within candidate genes: APC, beta-Catenin, K-RAS, DCC, P16, PTEN, RB1, P15, APOE, ERCC2, P53, MTHFR and hMSH2. Genotyping of consecutive, unselected colorectal cancer patients was done mostly by utilizing the MassARRAY technology (Sequenom) and to a lesser extent DGGE, ARMS and direct DNA sequencing. Correlation of genotypes with specific phenotypic features was carried out for all patients and separately for the Ashkenazim. Overall, 456 patients were analyzed, the majority (64.25%) being of Ashkenazi origin; mean age at diagnosis was 65.6 +/- 14 (range 25-90 years), and the mean follow-up was 4.7 +/- 0.28 (range 0-30 years). Statistically significant associations were noted between SNPs in beta-catenin and APOE and a positive family history of cancer (beta-catenin: p=0.034, APOE: p=0.033); tumor location and a DCC SNP (p=0.038) and the P53 R72P mutation and survival (p=0.0336). In Ashkenazi patients, ERCC2 and MTHFR genes' SNPs were associated with age at diagnosis (ERCC2: p=0.025, MTHFR: p=0.0005); a P53 polymorphism, APOE and Rb SNPs with a family history of cancer (P53 p=0.034;APOE p=0.04, Rb p= 0.022); DCC SNP with tumor location (p=0.014); and p15 SNP with tumor grade (p=0.032). This preliminary study shows that genetic factors play a role in determining CRC phenotypic features and that a larger cohort with longer follow-up is clearly needed.

The olfactory receptor universe - from whole genome analysis to structure and evolution

Olfactory receptors (ORs) constitute the largest gene-family in the vertebrate genome. We have attempted to provide a comprehensive view of the OR universe through diverse tools of bioinformatics and computational biology. Among others, we have constructed the Human Olfactory Receptor Data Exploratorium (HORDE, http://bioportal.weizmann.ac.il/HORDE/) as a free online resource, which integrates information on ORs from different species. We studied the genomic organization of 853 human ORs and divided the repertoire into 135 clusters, accessible through our new cluster viewer feature. An analysis of intact and pseudogenized ORs in different clusters, as well as of OR expression patterns, is provided, relevant to OR transcription control. Coding single nucleotide polymorphisms were integrated; these are to be used for genotype-phenotype correlation studies. HORDE allows a unique opportunity for discerning protein structural and functional information of the individual OR proteins. By applying novel data analysis strategies to the >3000 OR genes of mouse, dog and human within HORDE, we have generated a set of refined rhodopsin-based homology models for ORs. For model improvement, we employed a novel analysis of specific positions along the seven transmembrane helices at which prolines generate helix-breaking kinks. The model shows family-specif ic structural features, including idiosyncratic kink patterns, which lead to significant differences in the inferred odorant binding site structure. Such analyses form a basis for a comprehensive sequence-based classification of OR proteins in terms of potential odorant binding specificities.

Congenital dyserythropoietic anemia type I is caused by mutations in codanin-1.

Congenital dyserythropoietic anemias (CDAs) constitute a rare group of inherited red-blood-cell disorders associated with dysplastic changes in late erythroid precursors. CDA type I (CDAI [MIM 224120], gene symbol CDAN1) is characterized by erythroid pathological features such as internuclear chromatin bridges, spongy heterochromatin, and invagination of the nuclear membrane, carrying cytoplasmic organelles into the nucleus. A cluster of 45 highly inbred Israeli Bedouin with CDAI enabled the mapping of the CDAN1 disease gene to a 2-Mb interval, now refined to 1.2 Mb, containing 15 candidate genes on human chromosome 15q15 (Tamary et al. 1998). After the characterization and exclusion of 13 of these genes, we identified the CDAN1 gene through 12 different mutations in 9 families with CDAI. This 28-exon gene, which is transcribed ubiquitously into 4738 nt mRNA, was reconstructed on the basis of gene prediction and homology searches. It encodes codanin-1, a putative o-glycosylated protein of 1,226 amino acids, with no obvious transmembrane domains. Codanin-1 has a 150-residue amino-terminal domain with sequence similarity to collagens and two shorter segments that show weak similarities to the microtubule-associated proteins, MAP1B (neuraxin) and synapsin. These findings, and the cellular phenotype, suggest that codanin-1 may be involved in nuclear envelope integrity, conceivably related to microtubule attachments. The specific mechanisms by which codanin-1 underlies normal erythropoiesis remain to be elucidated.