Gene: a gene-centered information resource at NCBI

GEAR: A database of Genomic Elements Associated with drug Resistance

Drug resistance is becoming a serious problem that leads to the failure of standard treatments, which is generally developed because of genetic mutations of certain molecules. Here, we present GEAR (A database of Genomic Elements Associated with drug Resistance) that aims to provide comprehensive information about genomic elements (including genes, single-nucleotide polymorphisms and microRNAs) that are responsible for drug resistance. Right now, GEAR contains 1631 associations between 201 human drugs and 758 genes, 106 associations between 29 human drugs and 66 miRNAs, and 44 associations between 17 human drugs and 22 SNPs. These relationships are firstly extracted from primary literature with text mining and then manually curated. The drug resistome deposited in GEAR provides insights into the genetic factors underlying drug resistance. In addition, new indications and potential drug combinations can be identified based on the resistome. The GEAR database can be freely accessed through http://gear.comp-sysbio.org.

A Review of Bioinformatics Tools for Bio-Prospecting from Metagenomic Sequence Data

The microbiome can be defined as the community of microorganisms that live in a particular environment. Metagenomics is the practice of sequencing DNA from the genomes of all organisms present in a particular sample, and has become a common method for the study of microbiome population structure and function. Increasingly, researchers are finding novel genes encoded within metagenomes, many of which may be of interest to the biotechnology and pharmaceutical industries. However, such “bioprospecting” requires a suite of sophisticated bioinformatics tools to make sense of the data. This review summarizes the most commonly used bioinformatics tools for the assembly and annotation of metagenomic sequence data with the aim of discovering novel genes.

HRGFish: A database of hypoxia responsive genes in fishes

Several studies have highlighted the changes in the gene expression due to the hypoxia response in fishes, but the systematic organization of the information and the analytical platform for such genes are lacking. In the present study, an attempt was made to develop a database of hypoxia responsive genes in fishes (HRGFish), integrated with analytical tools, using LAMPP technology. Genes reported in hypoxia response for fishes were compiled through literature survey and the database presently covers 818 gene sequences and 35 gene types from 38 fishes. The upstream fragments (3,000 bp), covered in this database, enables to compute CG dinucleotides frequencies, motif finding of the hypoxia response element, identification of CpG island and mapping with the reference promoter of zebrafish. The database also includes functional annotation of genes and provides tools for analyzing sequences and designing primers for selected gene fragments. This may be the first database on the hypoxia response genes in fishes that provides a workbench to the scientific community involved in studying the evolution and ecological adaptation of the fish species in relation to hypoxia.

A New Chicken Genome Assembly Provides Insight into Avian Genome Structure

The importance of the Gallus gallus (chicken) as a model organism and agricultural animal merits a continuation of sequence assembly improvement efforts. We present a new version of the chicken genome assembly (Gallus_gallus-5.0; GCA_000002315.3), built from combined long single molecule sequencing technology, finished BACs, and improved physical maps. In overall assembled bases, we see a gain of 183 Mb, including 16.4 Mb in placed chromosomes with a corresponding gain in the percentage of intact repeat elements characterized. Of the 1.21 Gb genome, we include three previously missing autosomes, GGA30, 31, and 33, and improve sequence contig length 10-fold over the previous Gallus_gallus-4.0. Despite the significant base representation improvements made, 138 Mb of sequence is not yet located to chromosomes. When annotated for gene content, Gallus_gallus-5.0 shows an increase of 4679 annotated genes (2768 noncoding and 1911 protein-coding) over those in Gallus_gallus-4.0. We also revisited the question of what genes are missing in the avian lineage, as assessed by the highest quality avian genome assembly to date, and found that a large fraction of the original set of missing genes are still absent in sequenced bird species. Finally, our new data support a detailed map of MHC-B, encompassing two segments: one with a highly stable gene copy number and another in which the gene copy number is highly variable. The chicken model has been a critical resource for many other fields of study, and this new reference assembly will substantially further these efforts.

PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants

With the goal of providing a comprehensive, high-quality resource for both plant transcription factors (TFs) and their regulatory interactions with target genes, we upgraded plant TF database PlantTFDB to version 4.0 (http://planttfdb.cbi.pku.edu.cn/). In the new version, we identified 320 370 TFs from 165 species, presenting a more comprehensive genomic TF repertoires of green plants. Besides updating the pre-existing abundant functional and evolutionary annotation for identified TFs, we generated three new types of annotation which provide more directly clues to investigate functional mechanisms underlying: (i) a set of high-quality, non-redundant TF binding motifs derived from experiments; (ii) multiple types of regulatory elements identified from high-throughput sequencing data; (iii) regulatory interactions curated from literature and inferred by combining TF binding motifs and regulatory elements. In addition, we upgraded previous TF prediction server, and set up four novel tools for regulation prediction and functional enrichment analyses. Finally, we set up a novel companion portal PlantRegMap (http://plantregmap.cbi.pku.edu.cn) for users to access the regulation resource and analysis tools conveniently.

Target Safety Assessment: Strategies and Resources

An in-depth evaluation of target safety is an invaluable resource throughout drug discovery and development. The goal of a target safety evaluation is to identify potential unintended adverse consequences of target modulation, and to propose a risk evaluation and mitigation strategy to shepherd compounds through the discovery and development pipeline, to confirm and characterize unavoidable on-target toxicities in a timely manner to assist in early program advancement decisions, and to anticipate, monitor, and manage potential clinical adverse events. The role of an experienced discovery toxicologist in synthesizing the available information into an actionable set of recommendations for a safety evaluation strategy is critical to its successful application in early discovery programs. This chapter presents a summary of some of the information types and sources that should be investigated, and approaches that can be taken to generate an early assessment of potential safety liabilities.

A scored human protein-protein interaction network to catalyze genomic interpretation

Genome-scale human protein-protein interaction networks are critical to understanding cell biology and interpreting genomic data, but challenging to produce experimentally. Through data integration and quality control, we provide a scored human protein-protein interaction network (InWeb_InBioMap, or InWeb_IM) with severalfold more interactions (>500,000) and better functional biological relevance than comparable resources. We illustrate that InWeb_InBioMap enables functional interpretation of >4,700 cancer genomes and genes involved in autism.

Root-Knot and Cyst Nematodes Activate Procambium-Associated Genes in Arabidopsis Roots

Developmental plasticity is one of the most striking features of plant morphogenesis, as plants are able to vary their shapes in response to environmental cues. Biotic or abiotic stimuli often promote organogenesis events in plants not observed under normal growth conditions. Root-knot nematodes (RKNs) are known to parasitize multiple species of rooting plants and to induce characteristic tissue expansion called galls or root-knots on the roots of their hosts by perturbing the plant cellular machinery. Galls contain giant cells (GCs) and neighboring cells, and the GCs are a source of nutrients for the parasitizing nematode. Highly active cell proliferation was observed in galls. However, the underlying mechanisms that regulate the symptoms triggered by the plant-nematode interaction have not yet been elucidated. In this study, we deciphered the molecular mechanism of gall formation with an in vitro infection assay system using RKN Meloidogyne incognita, and the model plant Arabidopsis thaliana. By taking advantages of this system, we performed next-generation sequencing-based transcriptome profiling, and found that the expression of procambium identity-associated genes were enriched during gall formation. Clustering analyses with artificial xylogenic systems, together with the results of expression analyses of the candidate genes, showed a significant correlation between the induction of gall cells and procambium-associated cells. Furthermore, the promoters of several procambial marker genes such as ATHB8, TDR and WOX4 were activated not only in M. incognita-induced galls, but similarly in M. javanica induced-galls and Heterodera schachtii-induced syncytia. Our findings suggest that phytoparasitic nematodes modulate the host's developmental regulation of the vascular stem cells during gall formation.

On-the-fly selection of cell-specific enhancers, genes, miRNAs and proteins across the human body using SlideBase

Genomics consortia have produced large datasets profiling the expression of genes, micro-RNAs, enhancers and more across human tissues or cells. There is a need for intuitive tools to select subsets of such data that is the most relevant for specific studies. To this end, we present SlideBase, a web tool which offers a new way of selecting genes, promoters, enhancers and microRNAs that are preferentially expressed/used in a specified set of cells/tissues, based on the use of interactive sliders. With the help of sliders, SlideBase enables users to define custom expression thresholds for individual cell types/tissues, producing sets of genes, enhancers etc. which satisfy these constraints. Changes in slider settings result in simultaneous changes in the selected sets, updated in real time. SlideBase is linked to major databases from genomics consortia, including FANTOM, GTEx, The Human Protein Atlas and BioGPS.Database URL: http://slidebase.binf.ku.dk.

Integrative methods for analyzing big data in precision medicine

We provide an overview of recent developments in big data analyses in the context of precision medicine and health informatics. With the advance in technologies capturing molecular and medical data, we entered the area of "Big Data" in biology and medicine. These data offer many opportunities to advance precision medicine. We outline key challenges in precision medicine and present recent advances in data integration-based methods to uncover personalized information from big data produced by various omics studies. We survey recent integrative methods for disease subtyping, biomarkers discovery, and drug repurposing, and list the tools that are available to domain scientists. Given the ever-growing nature of these big data, we highlight key issues that big data integration methods will face.

OrthoDB v9.1: cataloging evolutionary and functional annotations for animal, fungal, plant, archaeal, bacterial and viral orthologs

OrthoDB is a comprehensive catalog of orthologs, genes inherited by extant species from a single gene in their last common ancestor. In 2016 OrthoDB reached its 9th release, growing to over 22 million genes from over 5000 species, now adding plants, archaea and viruses. In this update we focused on usability of this fast-growing wealth of data: updating the user and programmatic interfaces to browse and query the data, and further enhancing the already extensive integration of available gene functional annotations. Collating functional annotations from over 100 resources, and enabled us to propose descriptive titles for 87% of ortholog groups. Additionally, OrthoDB continues to provide computed evolutionary annotations and to allow user queries by sequence homology. The OrthoDB resource now enables users to generate publication-quality comparative genomics charts, as well as to upload, analyze and interactively explore their own private data. OrthoDB is available from http://orthodb.org.

Genenames.org: the HGNC and VGNC resources in 2017

The HUGO Gene Nomenclature Committee (HGNC) based at the European Bioinformatics Institute (EMBL-EBI) assigns unique symbols and names to human genes. Currently the HGNC database contains almost 40 000 approved gene symbols, over 19 000 of which represent protein-coding genes. In addition to naming genomic loci we manually curate genes into family sets based on shared characteristics such as homology, function or phenotype. We have recently updated our gene family resources and introduced new improved visualizations which can be seen alongside our gene symbol reports on our primary website http://www.genenames.org. In 2016 we expanded our remit and formed the Vertebrate Gene Nomenclature Committee (VGNC) which is responsible for assigning names to vertebrate species lacking a dedicated nomenclature group. Using the chimpanzee genome as a pilot project we have approved symbols and names for over 14 500 protein-coding genes in chimpanzee, and have developed a new website http://vertebrate.genenames.org to distribute these data. Here, we review our online data and resources, focusing particularly on the improvements and new developments made during the last two years.

PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants

With the goal of providing a comprehensive, high-quality resource for both plant transcription factors (TFs) and their regulatory interactions with target genes, we upgraded plant TF database PlantTFDB to version 4.0 (http://planttfdb.cbi.pku.edu.cn/). In the new version, we identified 320 370 TFs from 165 species, presenting a more comprehensive genomic TF repertoires of green plants. Besides updating the pre-existing abundant functional and evolutionary annotation for identified TFs, we generated three new types of annotation which provide more directly clues to investigate functional mechanisms underlying: (i) a set of high-quality, non-redundant TF binding motifs derived from experiments; (ii) multiple types of regulatory elements identified from high-throughput sequencing data; (iii) regulatory interactions curated from literature and inferred by combining TF binding motifs and regulatory elements. In addition, we upgraded previous TF prediction server, and set up four novel tools for regulation prediction and functional enrichment analyses. Finally, we set up a novel companion portal PlantRegMap (http://plantregmap.cbi.pku.edu.cn) for users to access the regulation resource and analysis tools conveniently.

Systems genetics analysis of pharmacogenomics variation during antidepressant treatment

Selective serotonin reuptake inhibitors (SSRIs) are the most widely used antidepressants, but the efficacy of the treatment varies significantly among individuals. It is believed that complex genetic mechanisms play a part in this variation. We have used a network based approach to unravel the involved genetic components. Moreover, we investigated the potential difference in the genetic interaction networks underlying SSRI treatment response over time. We found four hub genes (ASCC3, PPARGC1B, SCHIP1 and TMTC2) with different connectivity in the initial SSRI treatment period (baseline to week 4) compared with the subsequent period (4-8 weeks after initiation), suggesting that different genetic networks are important at different times during SSRI treatment. The strongest interactions in the initial SSRI treatment period involved genes encoding transcriptional factors, and in the subsequent period genes involved in calcium homeostasis. In conclusion, we suggest a difference in genetic interaction networks between initial and subsequent SSRI response.

AAgAtlas 1.0: a human autoantigen database

Autoantibodies refer to antibodies that target self-antigens, which can play pivotal roles in maintaining homeostasis, distinguishing normal from tumor tissue and trigger autoimmune diseases. In the last three decades, tremendous efforts have been devoted to elucidate the generation, evolution and functions of autoantibodies, as well as their target autoantigens. However, reports of these countless previously identified autoantigens are randomly dispersed in the literature. Here, we constructed an AAgAtlas database 1.0 using text-mining and manual curation. We extracted 45 830 autoantigen-related abstracts and 94 313 sentences from PubMed using the keywords of either ‘autoantigen’ or ‘autoantibody’ or their lexical variants, which were further refined to 25 520 abstracts, 43 253 sentences and 3984 candidates by our bio-entity recognizer based on the Protein Ontology. Finally, we identified 1126 genes as human autoantigens and 1071 related human diseases, with which we constructed a human autoantigen database (AAgAtlas database 1.0). The database provides a user-friendly interface to conveniently browse, retrieve and download human autoantigens as well as their associated diseases. The database is freely accessible at http://biokb.ncpsb.org/aagatlas/. We believe this database will be a valuable resource to track and understand human autoantigens as well as to investigate their functions in basic and translational research.

Network diffusion-based analysis of high-throughput data for the detection of differentially enriched modules

A relation exists between network proximity of molecular entities in interaction networks, functional similarity and association with diseases. The identification of network regions associated with biological functions and pathologies is a major goal in systems biology. We describe a network diffusion-based pipeline for the interpretation of different types of omics in the context of molecular interaction networks. We introduce the network smoothing index, a network-based quantity that allows to jointly quantify the amount of omics information in genes and in their network neighbourhood, using network diffusion to define network proximity. The approach is applicable to both descriptive and inferential statistics calculated on omics data. We also show that network resampling, applied to gene lists ranked by quantities derived from the network smoothing index, indicates the presence of significantly connected genes. As a proof of principle, we identified gene modules enriched in somatic mutations and transcriptional variations observed in samples of prostate adenocarcinoma (PRAD). In line with the local hypothesis, network smoothing index and network resampling underlined the existence of a connected component of genes harbouring molecular alterations in PRAD.

MINERVA-a platform for visualization and curation of molecular interaction networks

Our growing knowledge about various molecular mechanisms is becoming increasingly more structured and accessible. Different repositories of molecular interactions and available literature enable construction of focused and high-quality molecular interaction networks. Novel tools for curation and exploration of such networks are needed, in order to foster the development of a systems biology environment. In particular, solutions for visualization, annotation and data cross-linking will facilitate usage of network-encoded knowledge in biomedical research. To this end we developed the MINERVA (Molecular Interaction NEtwoRks VisuAlization) platform, a standalone webservice supporting curation, annotation and visualization of molecular interaction networks in Systems Biology Graphical Notation (SBGN)-compliant format. MINERVA provides automated content annotation and verification for improved quality control. The end users can explore and interact with hosted networks, and provide direct feedback to content curators. MINERVA enables mapping drug targets or overlaying experimental data on the visualized networks. Extensive export functions enable downloading areas of the visualized networks as SBGN-compliant models for efficient reuse of hosted networks. The software is available under Affero GPL 3.0 as a Virtual Machine snapshot, Debian package and Docker instance at http://r3lab.uni.lu/web/minerva-website/. We believe that MINERVA is an important contribution to systems biology community, as its architecture enables set-up of locally or globally accessible SBGN-oriented repositories of molecular interaction networks. Its functionalities allow overlay of multiple information layers, facilitating exploration of content and interpretation of data. Moreover, annotation and verification workflows of MINERVA improve the efficiency of curation of networks, allowing life-science researchers to better engage in development and use of biomedical knowledge repositories.

Meta-Analysis of Parkinson's Disease Transcriptome Data Using TRAM Software: Whole Substantia Nigra Tissue and Single Dopamine Neuron Differential Gene Expression

The understanding of the genetic basis of the Parkinson's disease (PD) and the correlation between genotype and phenotype has revolutionized our knowledge about the pathogenetic mechanisms of neurodegeneration, opening up exciting new therapeutic and neuroprotective perspectives. Genomic knowledge of PD is still in its early stages and can provide a good start for studies of the molecular mechanisms that underlie the gene expression variations and the epigenetic mechanisms that may contribute to the complex and characteristic phenotype of PD. In this study we used the software TRAM (Transcriptome Mapper) to analyse publicly available microarray data of a total of 151 PD patients and 130 healthy controls substantia nigra (SN) samples, to identify chromosomal segments and gene loci differential expression. In particular, we separately analyzed PD patients and controls data from post-mortem snap-frozen SN whole tissue and from laser microdissected midbrain dopamine (DA) neurons, to better characterize the specific DA neuronal expression profile associated with the late-stage Parkinson's condition. The default "Map" mode analysis resulted in 10 significantly over/under-expressed segments, mapping on 8 different chromosomes for SN whole tissue and in 4 segments mapping on 4 different chromosomes for DA neurons. In conclusion, TRAM software allowed us to confirm the deregulation of some genomic regions and loci involved in key molecular pathways related to neurodegeneration, as well as to provide new insights about genes and non-coding RNA transcripts not yet associated with the disease.

The Chromatin Remodelling Protein CHD1 Contains a Previously Unrecognised C-Terminal Helical Domain

The packaging of eukaryotic DNA into nucleosomes, and the organisation of these nucleosomes into chromatin, plays a critical role in regulating all DNA-associated processes. Chromodomain helicase DNA-binding protein 1 (CHD1) is an ATP-dependent chromatin remodelling protein that is conserved throughout eukaryotes and has an ability to assemble and organise nucleosomes both in vitro and in vivo. This activity is involved in the regulation of transcription and is implicated in mammalian development and stem cell biology. CHD1 is classically depicted as possessing a pair of tandem chromodomains that directly precede a core catalytic helicase-like domain that is then followed by a SANT-SLIDE DNA-binding domain. Here, we have identified an additional conserved domain C-terminal to the SANT-SLIDE domain and determined its structure by multidimensional heteronuclear NMR spectroscopy. We have termed this domain the CHD1 helical C-terminal (CHCT) domain as it is comprised of five α-helices arranged in a variant helical bundle topology. CHCT has a conserved, positively charged surface and is able to bind DNA and nucleosomes. In addition, we have identified another group of proteins, the as yet uncharacterised C17orf64 proteins, as also containing a conserved CHCT domain. Our data provide new structural insights into the CHD1 enzyme family.

Training and evaluation corpora for the extraction of causal relationships encoded in biological expression language (BEL)

Success in extracting biological relationships is mainly dependent on the complexity of the task as well as the availability of high-quality training data. Here, we describe the new corpora in the systems biology modeling language BEL for training and testing biological relationship extraction systems that we prepared for the BioCreative V BEL track. BEL was designed to capture relationships not only between proteins or chemicals, but also complex events such as biological processes or disease states. A BEL nanopub is the smallest unit of information and represents a biological relationship with its provenance. In BEL relationships (called BEL statements), the entities are normalized to defined namespaces mainly derived from public repositories, such as sequence databases, MeSH or publicly available ontologies. In the BEL nanopubs, the BEL statements are associated with citation information and supportive evidence such as a text excerpt. To enable the training of extraction tools, we prepared BEL resources and made them available to the community. We selected a subset of these resources focusing on a reduced set of namespaces, namely, human and mouse genes, ChEBI chemicals, MeSH diseases and GO biological processes, as well as relationship types ‘increases’ and ‘decreases’. The published training corpus contains 11 000 BEL statements from over 6000 supportive text excerpts. For method evaluation, we selected and re-annotated two smaller subcorpora containing 100 text excerpts. For this re-annotation, the inter-annotator agreement was measured by the BEL track evaluation environment and resulted in a maximal F-score of 91.18% for full statement agreement. In addition, for a set of 100 BEL statements, we do not only provide the gold standard expert annotations, but also text excerpts pre-selected by two automated systems. Those text excerpts were evaluated and manually annotated as true or false supportive in the course of the BioCreative V BEL track task.

Database URL:http://wiki.openbel.org/display/BIOC/Datasets

Family-based exome-wide assessment of maternal genetic effects on susceptibility to childhood B-cell acute lymphoblastic leukemia in hispanics

BACKGROUND

Children of Hispanic ancestry have a higher incidence of acute lymphoblastic leukemia (ALL) compared with other ethnic groups, but to the authors' knowledge, the genetic basis for these racial disparities remain incompletely understood. Genome-wide association studies of childhood ALL to date have focused on inherited genetic effects; however, maternal genetic effects (the role of the maternal genotype on phenotype development in the offspring) also may play a role in ALL susceptibility.

METHODS

The authors conducted a family-based exome-wide association study of maternal genetic effects among Hispanics with childhood B-cell ALL using the Illumina Infinium HumanExome BeadChip. A discovery cohort of 312 Guatemalan and Hispanic American families and an independent replication cohort of 152 Hispanic American families were used.

RESULTS

Three maternal single-nucleotide polymorphisms (SNPs) approached the study threshold for significance after correction for multiple testing (P<1.0 × 10^-6 ): MTL5 rs12365708 (testis expressed metallothionein-like protein [tesmin]) (relative risk [RR], 2.62; 95% confidence interval [95% CI], 1.61-4.27 [P = 1.8 × 10^-5 ]); ALKBH1 rs6494 (AlkB homolog 1, histone H2A dioxygenase) (RR, 3.77; 95% CI, 1.84-7.74 [P = 3.7 × 10^-5 ]); and NEUROG3 rs4536103 (neurogenin 3) (RR, 1.75; 95% CI, 1.30-2.37 [P = 1.2 × 10^-4 ]). Although effect sizes were similar, these SNPs were not nominally significant in the replication cohort in the current study. In a meta-analysis comprised of the discovery cohort and the replication cohort, these SNPs were still not found to be statistically significant after correction for multiple comparisons (rs12365708: pooled RR, 2.27 [95% CI, 1.48-3.50], P = 1.99 × 10^-4 ; rs6494: pooled RR, 2.31 [95% CI, 1.38-3.85], P = .001; and rs4536103: pooled RR, 1.67 [95% CI, 1.29-2.16] P = 9.23 × 10^-5 ).

CONCLUSIONS

In what to the authors' knowledge is the first family-based based exome-wide association study to investigate maternal genotype effects associated with childhood ALL, the results did not implicate a strong role of maternal genotype on disease risk among Hispanics; however, 3 maternal SNPs were identified that may play a modest role in susceptibility. Cancer 2016;122:3697-704. © 2016 American Cancer Society.

The Unique Evolutionary Signature of Genes Associated with Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a common heritable neurodevelopmental disorder, which is characterized by communication and social deficits that reduce the reproductive fitness of individuals with the disorder. Here, we studied the genomic characteristics of 651 ASD genes in a whole-exome sequencing dataset, to search for traces of the evolutionary forces that helped maintain ASD in the human population. We show that ASD genes are ~65 longer and ~20 % less variable than non-ASD genes. The mutational shortage in ASD genes was particularly eminent when considering only deleterious genetic variations, which is a hallmark of negative selection. We further show that these genomic characteristics are unique to ASD genes, as compared with brain-specific genes or with genes of other diseases. Our findings suggest that ASD genes have evolved under complex evolutionary forces, which have left a unique signature that can be used to identify new candidate ASD genes.

DTMiner: identification of potential disease targets through biomedical literature mining

Motivation: Biomedical researchers often search through massive catalogues of literature to look for potential relationships between genes and diseases. Given the rapid growth of biomedical literature, automatic relation extraction, a crucial technology in biomedical literature mining, has shown great potential to support research of gene-related diseases. Existing work in this field has produced datasets that are limited both in scale and accuracy.

Results: In this study, we propose a reliable and efficient framework that takes large biomedical literature repositories as inputs, identifies credible relationships between diseases and genes, and presents possible genes related to a given disease and possible diseases related to a given gene. The framework incorporates name entity recognition (NER), which identifies occurrences of genes and diseases in texts, association detection whereby we extract and evaluate features from gene–disease pairs, and ranking algorithms that estimate how closely the pairs are related. The F1-score of the NER phase is 0.87, which is higher than existing studies. The association detection phase takes drastically less time than previous work while maintaining a comparable F1-score of 0.86. The end-to-end result achieves a 0.259 F1-score for the top 50 genes associated with a disease, which performs better than previous work. In addition, we released a web service for public use of the dataset.

Availability and Implementation: The implementation of the proposed algorithms is publicly available at http://gdr-web.rwebox.com/public_html/index.php?page=download.php. The web service is available at http://gdr-web.rwebox.com/public_html/index.php.

Contact:jenny.wei@astrazeneca.com or kzhu@cs.sjtu.edu.cn

Supplementary information:Supplementary data are available at Bioinformatics online.

A gene browser of colorectal cancer with literature evidence and pre-computed regulatory information to identify key tumor suppressors and oncogenes

Colorectal cancer (CRC) is a cancer of growing incidence that associates with a high mortality rate worldwide. There is a poor understanding of the heterogeneity of CRC with regard to causative genetic mutations and gene regulatory mechanisms. Previous studies have identified several susceptibility genes in small-scale experiments. However, the information has not been comprehensively and systematically compiled and interpreted. In this study, we constructed the gbCRC, the first literature-based gene resource for investigating CRC-related human genes. The features of our database include: (i) manual curation of experimentally-verified genes reported in the literature; (ii) comprehensive integration of five reliable data sources; and (iii) pre-computed regulatory patterns involving transcription factors, microRNAs and long non-coding RNAs. In total, 2067 genes associating with 2819 PubMed abstracts were compiled. Comprehensive functional annotations associated with all the genes, including gene expression profiles, homologous genes in other model species, protein-protein interactions, somatic mutations, and potential methylation sites. These comprehensive annotations and this pre-computed regulatory information highlighted the importance of the gbCRC with regard to the unexplored regulatory network of CRC. This information is available in a plain text format that is free to download.

Subfunctionalization of peroxisome proliferator response elements accounts for retention of duplicated fabp1 genes in zebrafish

Background

In the duplication-degeneration-complementation (DDC) model, a duplicated gene has three possible fates: it may lose functionality through the accumulation of mutations (nonfunctionalization), acquire a new function (neofunctionalization), or each duplicate gene may retain a subset of functions of the ancestral gene (subfunctionalization). The role that promoter evolution plays in retention of duplicated genes in eukaryotic genomes is not well understood. Fatty acid-binding proteins (Fabp) belong to a multigene family that are highly conserved in sequence and function, but differ in their gene regulation, suggesting selective pressure is exerted via regulatory elements in the promoter.

Results

In this study, we describe the PPAR regulation of zebrafish fabp1a, fabp1b.1, and fabp1b.2 promoters and compare them to the PPAR regulation of the spotted gar fabp1 promoter, representative of the ancestral fabp1 gene. Evolution of the fabp1 promoter was inferred by sequence analysis, and differential PPAR-agonist activation of fabp1 promoter activity in zebrafish liver and intestine explant cells, and in HEK293A cells transiently transfected with wild-type and mutated fabp1promoter-reporter gene constructs. The promoter activity of spotted gar fabp1, representative of the ancestral fabp1, was induced by both PPARα- and PPARγ-specific agonists, but displayed a biphasic response to PPARα activation. Zebrafish fabp1a was PPARα-selective, fabp1b.1 was PPARγ-selective, and fabp1b.2 was not regulated by PPAR.

Conclusions

The zebrafish fabp1 promoters underwent two successive rounds of subfunctionalization with respect to PPAR regulation leading to retention of three zebrafish fabp1 genes with stimuli-specific regulation. Using a pharmacological approach, we demonstrated here the divergent regulation of the zebrafish fabp1a, fabp1b.1, and fabp1b.2 with regard to subfunctionalization of PPAR regulation following two rounds of gene duplication.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0717-x) contains supplementary material, which is available to authorized users.

Enrichr: a comprehensive gene set enrichment analysis web server 2016 update

Enrichment analysis is a popular method for analyzing gene sets generated by genome-wide experiments. Here we present a significant update to one of the tools in this domain called Enrichr. Enrichr currently contains a large collection of diverse gene set libraries available for analysis and download. In total, Enrichr currently contains 180 184 annotated gene sets from 102 gene set libraries. New features have been added to Enrichr including the ability to submit fuzzy sets, upload BED files, improved application programming interface and visualization of the results as clustergrams. Overall, Enrichr is a comprehensive resource for curated gene sets and a search engine that accumulates biological knowledge for further biological discoveries. Enrichr is freely available at: http://amp.pharm.mssm.edu/Enrichr.

Enrichr: a comprehensive gene set enrichment analysis web server 2016 update

Enrichment analysis is a popular method for analyzing gene sets generated by genome-wide experiments. Here we present a significant update to one of the tools in this domain called Enrichr. Enrichr currently contains a large collection of diverse gene set libraries available for analysis and download. In total, Enrichr currently contains 180 184 annotated gene sets from 102 gene set libraries. New features have been added to Enrichr including the ability to submit fuzzy sets, upload BED files, improved application programming interface and visualization of the results as clustergrams. Overall, Enrichr is a comprehensive resource for curated gene sets and a search engine that accumulates biological knowledge for further biological discoveries. Enrichr is freely available at: http://amp.pharm.mssm.edu/Enrichr.

The harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins

Genomics, epigenomics, transcriptomics, proteomics and metabolomics efforts rapidly generate a plethora of data on the activity and levels of biomolecules within mammalian cells. At the same time, curation projects that organize knowledge from the biomedical literature into online databases are expanding. Hence, there is a wealth of information about genes, proteins and their associations, with an urgent need for data integration to achieve better knowledge extraction and data reuse. For this purpose, we developed the Harmonizome: a collection of processed datasets gathered to serve and mine knowledge about genes and proteins from over 70 major online resources. We extracted, abstracted and organized data into ∼72 million functional associations between genes/proteins and their attributes. Such attributes could be physical relationships with other biomolecules, expression in cell lines and tissues, genetic associations with knockout mouse or human phenotypes, or changes in expression after drug treatment. We stored these associations in a relational database along with rich metadata for the genes/proteins, their attributes and the original resources. The freely available Harmonizome web portal provides a graphical user interface, a web service and a mobile app for querying, browsing and downloading all of the collected data. To demonstrate the utility of the Harmonizome, we computed and visualized gene-gene and attribute-attribute similarity networks, and through unsupervised clustering, identified many unexpected relationships by combining pairs of datasets such as the association between kinase perturbations and disease signatures. We also applied supervised machine learning methods to predict novel substrates for kinases, endogenous ligands for G-protein coupled receptors, mouse phenotypes for knockout genes, and classified unannotated transmembrane proteins for likelihood of being ion channels. The Harmonizome is a comprehensive resource of knowledge about genes and proteins, and as such, it enables researchers to discover novel relationships between biological entities, as well as form novel data-driven hypotheses for experimental validation.Database URL: http://amp.pharm.mssm.edu/Harmonizome.

Gene-Expression Profiling Suggests Impaired Signaling via the Interferon Pathway in Cstb-/- Microglia

Progressive myoclonus epilepsy of Unverricht-Lundborg type (EPM1, OMIM254800) is an autosomal recessive neurodegenerative disorder characterized by stimulus-sensitive and action-activated myoclonus, tonic-clonic epileptic seizures, and ataxia. Loss-of-function mutations in the gene encoding the cysteine protease inhibitor cystatin B (CSTB) underlie EPM1. The deficiency of CSTB in mice (Cstb^-/- mice) generates a phenotype resembling the symptoms of EPM1 patients and is accompanied by microglial activation at two weeks of age and an upregulation of immune system-associated genes in the cerebellum at one month of age. To shed light on molecular pathways and processes linked to CSTB deficiency in microglia we characterized the transcriptome of cultured Cstb^-/- mouse microglia using microarray hybridization and RNA sequencing (RNA-seq). The gene expression profiles obtained with these two techniques were in good accordance and not polarized to either pro- or anti-inflammatory status. In Cstb^-/- microglia, altogether 184 genes were differentially expressed. Of these, 33 genes were identified by both methods. Several interferon-regulated genes were weaker expressed in Cstb^-/- microglia compared to control. This was confirmed by quantitative real-time PCR of the transcripts Irf7 and Stat1. Subsequently, we explored the biological context of CSTB deficiency in microglia more deeply by functional enrichment and canonical pathway analysis. This uncovered a potential role for CSTB in chemotaxis, antigen-presentation, and in immune- and defense response-associated processes by altering JAK-STAT pathway signaling. These data support and expand the previously suggested involvement of inflammatory processes to the disease pathogenesis of EPM1 and connect CSTB deficiency in microglia to altered expression of interferon-regulated genes.

QuIN: A Web Server for Querying and Visualizing Chromatin Interaction Networks

Recent studies of the human genome have indicated that regulatory elements (e.g. promoters and enhancers) at distal genomic locations can interact with each other via chromatin folding and affect gene expression levels. Genomic technologies for mapping interactions between DNA regions, e.g., ChIA-PET and HiC, can generate genome-wide maps of interactions between regulatory elements. These interaction datasets are important resources to infer distal gene targets of non-coding regulatory elements and to facilitate prioritization of critical loci for important cellular functions. With the increasing diversity and complexity of genomic information and public ontologies, making sense of these datasets demands integrative and easy-to-use software tools. Moreover, network representation of chromatin interaction maps enables effective data visualization, integration, and mining. Currently, there is no software that can take full advantage of network theory approaches for the analysis of chromatin interaction datasets. To fill this gap, we developed a web-based application, QuIN, which enables: 1) building and visualizing chromatin interaction networks, 2) annotating networks with user-provided private and publicly available functional genomics and interaction datasets, 3) querying network components based on gene name or chromosome location, and 4) utilizing network based measures to identify and prioritize critical regulatory targets and their direct and indirect interactions. AVAILABILITY: QuIN’s web server is available at http://quin.jax.org QuIN is developed in Java and JavaScript, utilizing an Apache Tomcat web server and MySQL database and the source code is available under the GPLV3 license available on GitHub: https://github.com/UcarLab/QuIN/.

The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses

GeneCards, the human gene compendium, enables researchers to effectively navigate and inter-relate the wide universe of human genes, diseases, variants, proteins, cells, and biological pathways. Our recently launched Version 4 has a revamped infrastructure facilitating faster data updates, better-targeted data queries, and friendlier user experience. It also provides a stronger foundation for the GeneCards suite of companion databases and analysis tools. Improved data unification includes gene-disease links via MalaCards and merged biological pathways via PathCards, as well as drug information and proteome expression. VarElect, another suite member, is a phenotype prioritizer for next-generation sequencing, leveraging the GeneCards and MalaCards knowledgebase. It automatically infers direct and indirect scored associations between hundreds or even thousands of variant-containing genes and disease phenotype terms. VarElect's capabilities, either independently or within TGex, our comprehensive variant analysis pipeline, help prepare for the challenge of clinical projects that involve thousands of exome/genome NGS analyses. © 2016 by John Wiley & Sons, Inc.

GCGene: a gene resource for gastric cancer with literature evidence

Gastric cancer (GC) is the fifth most common cancer and third leading cause of cancer-related deaths worldwide. Its lethality primarily stems from a lack of detection strategies for early stages of GC and a lack of noninvasive detection strategies for advanced stages. The development of early diagnostic biomarkers largely depends on understanding the biological pathways and regulatory mechanisms associated with putative GC genes. Unfortunately, the GC-implicated genes that have been identified thus far are scattered among thousands of published studies, and no systematic summary is available, which hinders the development of a large-scale genetic screen. To provide a publically accessible resource tool to meet this need, we constructed a literature-based database GCGene (Gastric Cancer Gene database) with comprehensive annotations supported by a user-friendly website. In the current release, we have collected 1,815 unique human genes including 1,678 protein-coding and 137 non-coding genes curated from extensive examination of 3,142 PubMed abstracts. The resulting database has a convenient web-based interface to facilitate both textual and sequence-based searches. All curated genes in GCGene are downloadable for advanced bioinformatics data mining. Gene prioritization was performed to rank the relative relevance of these genes in GC development. The 100 top-ranked genes are highly mutated according to the cohort of published studies we reviewed. By conducting a network analysis of these top-ranked GC-associated genes in the human interactome, we were able to identify strong links between 8 highly connected genes with low expression and patient survival time. GCGene is freely available to academic users at http://gcgene.bioinfo-minzhao.org/.

Systematic analysis of the molecular mechanism underlying atherosclerosis using a text mining approach

Background

Atherosclerosis is one of the common health threats all over the world. It is a complex heritable disease that affects arterial blood vessels. Chronic inflammatory response plays an important role in atherogenesis. There has been little success in fully identifying functionally important genes in the pathogenesis of atherosclerosis.

Results

In the present study, we performed a systematic analysis of atherosclerosis-related genes using text mining. We identified a total of 1312 genes. Gene ontology (GO) analysis revealed that a total of 35 terms exhibited significance (p < 0.05) as overrepresented terms, indicating that atherosclerosis invokes many genes with a wide range of different functions. Pathway analysis demonstrated that the most highly enriched pathway is the Toll-like receptor signaling pathway. Finally, through gene network analysis, we prioritized 48 genes using the hub gene method.

Conclusions

Our study provides a valuable resource for the in-depth understanding of the mechanism underlying atherosclerosis.

Electronic supplementary material

The online version of this article (doi:10.1186/s40246-016-0075-1) contains supplementary material, which is available to authorized users.

Bioinformatics resources for pollen

Bioinformatics for Pollen. Pollen plays a key role in crop production, and its development is the most delicate phase in reproduction. Different metabolic pathways are involved in pollen development, and changes in the level of some metabolites, as well as responses to stress, are correlated with the reduction in pollen viability, leading consequently to a decrease in the fruit production. However, studies on pollen may be hard because gamete development and fertilization are complex processes that occur during a short window of time. The rise of the so-called -omics sciences provided key strategies to promote molecular research in pollen tissues, starting from model organisms and moving to increasing number of species. An integrated multi-level approach based on investigations from genomics, transcriptomics, proteomics and metabolomics appears now feasible to clarify key molecular processes in pollen development and viability. To this aim, bioinformatics has a fundamental role for data production and analysis, contributing varied and ad hoc methodologies, endowed with different sensitivity and specificity, necessary for extracting added-value information from the large amount of molecular data achievable. Bioinformatics is also essential for data management, organization, distribution and integration in suitable resources. This is necessary to catch the biological features of the pollen tissues and to design effective approaches to identifying structural or functional properties, enabling the modeling of the major involved processes in normal or in stress conditions. In this review, we provide an overview of the available bioinformatics resources for pollen, ranging from raw data collections to complete databases or platforms, when available, which include data and/or results from -omics efforts on the male gametophyte. Perspectives in the fields will also be described.

Generating Gene Ontology-Disease Inferences to Explore Mechanisms of Human Disease at the Comparative Toxicogenomics Database

Strategies for discovering common molecular events among disparate diseases hold promise for improving understanding of disease etiology and expanding treatment options. One technique is to leverage curated datasets found in the public domain. The Comparative Toxicogenomics Database (CTD; http://ctdbase.org/) manually curates chemical-gene, chemical-disease, and gene-disease interactions from the scientific literature. The use of official gene symbols in CTD interactions enables this information to be combined with the Gene Ontology (GO) file from NCBI Gene. By integrating these GO-gene annotations with CTD’s gene-disease dataset, we produce 753,000 inferences between 15,700 GO terms and 4,200 diseases, providing opportunities to explore presumptive molecular underpinnings of diseases and identify biological similarities. Through a variety of applications, we demonstrate the utility of this novel resource. As a proof-of-concept, we first analyze known repositioned drugs (e.g., raloxifene and sildenafil) and see that their target diseases have a greater degree of similarity when comparing GO terms vs. genes. Next, a computational analysis predicts seemingly non-intuitive diseases (e.g., stomach ulcers and atherosclerosis) as being similar to bipolar disorder, and these are validated in the literature as reported co-diseases. Additionally, we leverage other CTD content to develop testable hypotheses about thalidomide-gene networks to treat seemingly disparate diseases. Finally, we illustrate how CTD tools can rank a series of drugs as potential candidates for repositioning against B-cell chronic lymphocytic leukemia and predict cisplatin and the small molecule inhibitor JQ1 as lead compounds. The CTD dataset is freely available for users to navigate pathologies within the context of extensive biological processes, molecular functions, and cellular components conferred by GO. This inference set should aid researchers, bioinformaticists, and pharmaceutical drug makers in finding commonalities in disease mechanisms, which in turn could help identify new therapeutics, new indications for existing pharmaceuticals, potential disease comorbidities, and alerts for side effects.

ToxEvaluator: an integrated computational platform to aid the interpretation of toxicology study-related findings

Attempts are frequently made to investigate adverse findings from preclinical toxicology studies in order to better understand underlying toxicity mechanisms. These efforts often begin with limited information, including a description of the adverse finding, knowledge of the structure of the chemical associated with its cause and the intended pharmacological target. ToxEvaluator was developed jointly by Pfizer and the Comparative Toxicogenomics Database (http://ctdbase.org) team at North Carolina State University as an in silico platform to facilitate interpretation of toxicity findings in light of prior knowledge. Through the integration of a diverse set of in silico tools that leverage a number of public and proprietary databases, ToxEvaluator streamlines the process of aggregating and interrogating diverse sources of information. The user enters compound and target identifiers, and selects adverse event descriptors from a safety lexicon and mapped MeSH disease terms. ToxEvaluator provides a summary report with multiple distinct areas organized according to what target or structural aspects have been linked to the adverse finding, including primary pharmacology, structurally similar proprietary compounds, structurally similar public domain compounds, predicted secondary (i.e. off-target) pharmacology and known secondary pharmacology. Similar proprietary compounds and their associated in vivo toxicity findings are reported, along with a link to relevant supporting documents. For similar public domain compounds and interacting targets, ToxEvaluator integrates relationships curated in Comparative Toxicogenomics Database, returning all direct and inferred linkages between them. As an example of its utility, we demonstrate how ToxEvaluator rapidly identified direct (primary pharmacology) and indirect (secondary pharmacology) linkages between cerivastatin and myopathy.

Synteny Portal: a web-based application portal for synteny block analysis

Recent advances in next-generation sequencing technologies and genome assembly algorithms have enabled the accumulation of a huge volume of genome sequences from various species. This has provided new opportunities for large-scale comparative genomics studies. Identifying and utilizing synteny blocks, which are genomic regions conserved among multiple species, is key to understanding genomic architecture and the evolutionary history of genomes. However, the construction and visualization of such synteny blocks from multiple species are very challenging, especially for biologists with a lack of computational skills. Here, we present Synteny Portal, a versatile web-based application portal for constructing, visualizing and browsing synteny blocks. With Synteny Portal, users can easily (i) construct synteny blocks among multiple species by using prebuilt alignments in the UCSC genome browser database, (ii) visualize and download syntenic relationships as high-quality images, (iii) browse synteny blocks with genetic information and (iv) download the details of synteny blocks to be used as input for downstream synteny-based analyses, all in an intuitive and easy-to-use web-based interface. We believe that Synteny Portal will serve as a highly valuable tool that will enable biologists to easily perform comparative genomics studies by compensating limitations of existing tools. Synteny Portal is freely available at http://bioinfo.konkuk.ac.kr/synteny_portal.

High-performance web services for querying gene and variant annotation

Efficient tools for data management and integration are essential for many aspects of high-throughput biology. In particular, annotations of genes and human genetic variants are commonly used but highly fragmented across many resources. Here, we describe MyGene.info and MyVariant.info, high-performance web services for querying gene and variant annotation information. These web services are currently accessed more than three million times permonth. They also demonstrate a generalizable cloud-based model for organizing and querying biological annotation information. MyGene.info and MyVariant.info are provided as high-performance web services, accessible at http://mygene.info and http://myvariant.info . Both are offered free of charge to the research community.

Enrichr: a comprehensive gene set enrichment analysis web server 2016 update

Enrichment analysis is a popular method for analyzing gene sets generated by genome-wide experiments. Here we present a significant update to one of the tools in this domain called Enrichr. Enrichr currently contains a large collection of diverse gene set libraries available for analysis and download. In total, Enrichr currently contains 180 184 annotated gene sets from 102 gene set libraries. New features have been added to Enrichr including the ability to submit fuzzy sets, upload BED files, improved application programming interface and visualization of the results as clustergrams. Overall, Enrichr is a comprehensive resource for curated gene sets and a search engine that accumulates biological knowledge for further biological discoveries. Enrichr is freely available at: http://amp.pharm.mssm.edu/Enrichr.

Divergent evolution of cis-acting peroxisome proliferator-activated receptor elements that differentially control the tandemly duplicated fatty acid-binding protein genes, fabp1b.1 and fabp1b.2, in zebrafish

Gene duplication is thought to facilitate increasing complexity in the evolution of life. The fate of most duplicated genes is nonfunctionalization: functional decay resulting from the accumulation of mutations. According to the duplication-degeneration-complementation (DDC) model, duplicated genes are retained by subfunctionalization, where the functions of the ancestral gene are sub-divided between duplicate genes, or by neofunctionalization, where one of the duplicates acquires a new function. Here, we report the differential regulation of the zebrafish tandemly duplicated fatty acid-binding protein genes, fabp1b.1 and fabp1b.2, by peroxisome proliferator-activated receptors (PPAR). fabp1b.1 mRNA levels were induced in tissue explants of liver, but not intestine, by PPAR agonists. fabp1b.1 promoter activity was induced to a greater extent by rosiglitazone (PPARγ-selective agonist) compared to WY 14,643 (PPARα-selective agonist) in HEK293A cells. Mutation of a peroxisome proliferator response element (PPRE) at -1232 bp in the fabp1b.1 promoter reduced PPAR-dependent activation. fabp1b.2 promoter activity was not affected by PPAR agonists. Differential regulation of the duplicated fabp1b promoters may be the result of PPRE loss in fabp1b.2 during a meiotic crossing-over event. Retention of PPAR inducibility in fabp1b.1 and not fabp1b.2 suggests unique regulation and function of the fabp1b duplicates.

g:Profiler-a web server for functional interpretation of gene lists (2016 update)

Functional enrichment analysis is a key step in interpreting gene lists discovered in diverse high-throughput experiments. g:Profiler studies flat and ranked gene lists and finds statistically significant Gene Ontology terms, pathways and other gene function related terms. Translation of hundreds of gene identifiers is another core feature of g:Profiler. Since its first publication in 2007, our web server has become a popular tool of choice among basic and translational researchers. Timeliness is a major advantage of g:Profiler as genome and pathway information is synchronized with the Ensembl database in quarterly updates. g:Profiler supports 213 species including mammals and other vertebrates, plants, insects and fungi. The 2016 update of g:Profiler introduces several novel features. We have added further functional datasets to interpret gene lists, including transcription factor binding site predictions, Mendelian disease annotations, information about protein expression and complexes and gene mappings of human genetic polymorphisms. Besides the interactive web interface, g:Profiler can be accessed in computational pipelines using our R package, Python interface and BioJS component. g:Profiler is freely available at http://biit.cs.ut.ee/gprofiler/.

REGene: a literature-based knowledgebase of animal regeneration that bridge tissue regeneration and cancer

Regeneration is a common phenomenon across multiple animal phyla. Regeneration-related genes (REGs) are critical for fundamental cellular processes such as proliferation and differentiation. Identification of REGs and elucidating their functions may help to further develop effective treatment strategies in regenerative medicine. So far, REGs have been largely identified by small-scale experimental studies and a comprehensive characterization of the diverse biological processes regulated by REGs is lacking. Therefore, there is an ever-growing need to integrate REGs at the genomics, epigenetics, and transcriptome level to provide a reference list of REGs for regeneration and regenerative medicine research. Towards achieving this, we developed the first literature-based database called REGene (REgeneration Gene database). In the current release, REGene contains 948 human (929 protein-coding and 19 non-coding genes) and 8445 homologous genes curated from gene ontology and extensive literature examination. Additionally, the REGene database provides detailed annotations for each REG, including: gene expression, methylation sites, upstream transcription factors, and protein-protein interactions. An analysis of the collected REGs reveals strong links to a variety of cancers in terms of genetic mutation, protein domains, and cellular pathways. We have prepared a web interface to share these regeneration genes, supported by refined browsing and searching functions at http://REGene.bioinfo-minzhao.org/.

Copy Number Variants Associated with 14 Cases of Self-Injurious Behavior

Copy number variants (CNVs) were detected and analyzed in 14 probands with autism and intellectual disability with self-injurious behavior (SIB) resulting in tissue damage. For each proband we obtained a clinical history and detailed behavioral descriptions. Genetic anomalies were observed in all probands, and likely clinical significance could be established in four cases. This included two cases having novel, de novo copy number variants and two cases having variants likely to have functional significance. These cases included segmental trisomy 14, segmental monosomy 21, and variants predicted to disrupt the function of ZEB2 (encoding a transcription factor) and HTR2C (encoding a serotonin receptor). Our results identify variants in regions previously implicated in intellectual disability and suggest candidate genes that could contribute to the etiology of SIB.

Whole-exome sequencing in obsessive-compulsive disorder identifies rare mutations in immunological and neurodevelopmental pathways

Studies of rare genetic variation have identified molecular pathways conferring risk for developmental neuropsychiatric disorders. To date, no published whole-exome sequencing studies have been reported in obsessive-compulsive disorder (OCD). We sequenced all the genome coding regions in 20 sporadic OCD cases and their unaffected parents to identify rare de novo (DN) single-nucleotide variants (SNVs). The primary aim of this pilot study was to determine whether DN variation contributes to OCD risk. To this aim, we evaluated whether there is an elevated rate of DN mutations in OCD, which would justify this approach toward gene discovery in larger studies of the disorder. Furthermore, to explore functional molecular correlations among genes with nonsynonymous DN SNVs in OCD probands, a protein-protein interaction (PPI) network was generated based on databases of direct molecular interactions. We applied Degree-Aware Disease Gene Prioritization (DADA) to rank the PPI network genes based on their relatedness to a set of OCD candidate genes from two OCD genome-wide association studies (Stewart et al., 2013; Mattheisen et al., 2014). In addition, we performed a pathway analysis with genes from the PPI network. The rate of DN SNVs in OCD was 2.51 × 10(-8) per base per generation, significantly higher than a previous estimated rate in unaffected subjects using the same sequencing platform and analytic pipeline. Several genes harboring DN SNVs in OCD were highly interconnected in the PPI network and ranked high in the DADA analysis. Nearly all the DN SNVs in this study are in genes expressed in the human brain, and a pathway analysis revealed enrichment in immunological and central nervous system functioning and development. The results of this pilot study indicate that further investigation of DN variation in larger OCD cohorts is warranted to identify specific risk genes and to confirm our preliminary finding with regard to PPI network enrichment for particular biological pathways and functions.

Cell type-specific properties and environment shape tissue specificity of cancer genes

One of the biggest mysteries in cancer research remains why mutations in certain genes cause cancer only at specific sites in the human body. The poor correlation between the expression level of a cancer gene and the tissues in which it causes malignant transformations raises the question of which factors determine the tissue-specific effects of a mutation. Here, we explore why some cancer genes are associated only with few different cancer types (i.e., are specific), while others are found mutated in a large number of different types of cancer (i.e., are general). We do so by contrasting cellular functions of specific-cancer genes with those of general ones to identify properties that determine where in the body a gene mutation is causing malignant transformations. We identified different groups of cancer genes that did not behave as expected (i.e., DNA repair genes being tissue specific, immune response genes showing a bimodal specificity function or strong association of generally expressed genes to particular cancers). Analysis of these three groups demonstrates the importance of environmental impact for understanding why certain cancer genes are only involved in the development of some cancer types but are rarely found mutated in other types of cancer.

The clinical trial landscape in oncology and connectivity of somatic mutational profiles to targeted therapies

Background

Precision medicine in oncology relies on rapid associations between patient-specific variations and targeted therapeutic efficacy. Due to the advancement of genomic analysis, a vast literature characterizing cancer-associated molecular aberrations and relative therapeutic relevance has been published. However, data are not uniformly reported or readily available, and accessing relevant information in a clinically acceptable time-frame is a daunting proposition, hampering connections between patients and appropriate therapeutic options. One important therapeutic avenue for oncology patients is through clinical trials. Accordingly, a global view into the availability of targeted clinical trials would provide insight into strengths and weaknesses and potentially enable research focus. However, data regarding the landscape of clinical trials in oncology is not readily available, and as a result, a comprehensive understanding of clinical trial availability is difficult.

Results

To support clinical decision-making, we have developed a data loader and mapper that connects sequence information from oncology patients to data stored in an in-house database, the JAX Clinical Knowledgebase (JAX-CKB), which can be queried readily to access comprehensive data for clinical reporting via customized reporting queries. JAX-CKB functions as a repository to house expertly curated clinically relevant data surrounding our 358-gene panel, the JAX Cancer Treatment Profile (JAX CTP), and supports annotation of functional significance of molecular variants. Through queries of data housed in JAX-CKB, we have analyzed the landscape of clinical trials relevant to our 358-gene targeted sequencing panel to evaluate strengths and weaknesses in current molecular targeting in oncology. Through this analysis, we have identified patient indications, molecular aberrations, and targeted therapy classes that have strong or weak representation in clinical trials.

Conclusions

Here, we describe the development and disseminate system methods for associating patient genomic sequence data with clinically relevant information, facilitating interpretation and providing a mechanism for informing therapeutic decision-making. Additionally, through customized queries, we have the capability to rapidly analyze the landscape of targeted therapies in clinical trials, enabling a unique view into current therapeutic availability in oncology.

ECGene: A Literature-Based Knowledgebase of Endometrial Cancer Genes

Endometrial cancer (EC) ranks as the sixth common cancer for women worldwide. To better distinguish cancer subtypes and identify effective early diagnostic biomarkers, we need improved understanding of the biological mechanisms associated with EC dysregulated genes. Although there is a wealth of clinical and molecular information relevant to EC in the literature, there has been no systematic summary of EC‐implicated genes. In this study, we developed a literature‐based database ECGene (Endometrial Cancer Gene database) with comprehensive annotations. ECGene features manual curation of 414 genes from thousands of publications, results from eight EC gene expression datasets, precomputation of coexpressed long noncoding RNAs, and an EC‐implicated gene interactome. In the current release, we generated and comprehensively annotated a list of 458 EC‐implicated genes. We found the top‐ranked EC‐implicated genes are frequently mutated in The Cancer Genome Atlas (TCGA) tumor samples. Furthermore, systematic analysis of coexpressed lncRNAs provided insight into the important roles of lncRNA in EC development. ECGene has a user‐friendly Web interface and is freely available at http://ecgene.bioinfo‐minzhao.org/. As the first literature‐based online resource for EC, ECGene serves as a useful gateway for researchers to explore EC genetics.

Integrating 400 million variants from 80,000 human samples with extensive annotations: towards a knowledge base to analyze disease cohorts

BACKGROUND

Data from a plethora of high-throughput sequencing studies is readily available to researchers, providing genetic variants detected in a variety of healthy and disease populations. While each individual cohort helps gain insights into polymorphic and disease-associated variants, a joint perspective can be more powerful in identifying polymorphisms, rare variants, disease-associations, genetic burden, somatic variants, and disease mechanisms.

DESCRIPTION

We have set up a Reference Variant Store (RVS) containing variants observed in a number of large-scale sequencing efforts, such as 1000 Genomes, ExAC, Scripps Wellderly, UK10K; various genotyping studies; and disease association databases. RVS holds extensive annotations pertaining to affected genes, functional impacts, disease associations, and population frequencies. RVS currently stores 400 million distinct variants observed in more than 80,000 human samples.

CONCLUSIONS

RVS facilitates cross-study analysis to discover novel genetic risk factors, gene-disease associations, potential disease mechanisms, and actionable variants. Due to its large reference populations, RVS can also be employed for variant filtration and gene prioritization.

AVAILABILITY

A web interface to public datasets and annotations in RVS is available at https://rvs.u.hpc.mssm.edu/.

GEneSTATION 1.0: a synthetic resource of diverse evolutionary and functional genomic data for studying the evolution of pregnancy-associated tissues and phenotypes

Mammalian gestation and pregnancy are fast evolving processes that involve the interaction of the fetal, maternal and paternal genomes. Version 1.0 of the GEneSTATION database (http://genestation.org) integrates diverse types of omics data across mammals to advance understanding of the genetic basis of gestation and pregnancy-associated phenotypes and to accelerate the translation of discoveries from model organisms to humans. GEneSTATION is built using tools from the Generic Model Organism Database project, including the biology-aware database CHADO, new tools for rapid data integration, and algorithms that streamline synthesis and user access. GEneSTATION contains curated life history information on pregnancy and reproduction from 23 high-quality mammalian genomes. For every human gene, GEneSTATION contains diverse evolutionary (e.g. gene age, population genetic and molecular evolutionary statistics), organismal (e.g. tissue-specific gene and protein expression, differential gene expression, disease phenotype), and molecular data types (e.g. Gene Ontology Annotation, protein interactions), as well as links to many general (e.g. Entrez, PubMed) and pregnancy disease-specific (e.g. PTBgene, dbPTB) databases. By facilitating the synthesis of diverse functional and evolutionary data in pregnancy-associated tissues and phenotypes and enabling their quick, intuitive, accurate and customized meta-analysis, GEneSTATION provides a novel platform for comprehensive investigation of the function and evolution of mammalian pregnancy.