SnapShot: S-Phase Entry and Exit.

S-phase entry and exit are regulated by hundreds of protein complexes that assemble "just in time," orchestrated by a multitude of distinct events. To help understand their interplay, we have created a tailored visualization based on the Minardo layout, highlighting over 80 essential events. This complements our earlier visualization of M-phase, and both can be displayed together, giving a comprehensive overview of the events regulating the cell division cycle. To view this SnapShot, open or download the PDF.

A nondegenerate code of deleterious variants in Mendelian loci contributes to complex disease risk.

Although countless highly penetrant variants have been associated with Mendelian disorders, the genetic etiologies underlying complex diseases remain largely unresolved. By mining the medical records of over 110 million patients, we examine the extent to which Mendelian variation contributes to complex disease risk. We detect thousands of associations between Mendelian and complex diseases, revealing a nondegenerate, phenotypic code that links each complex disorder to a unique collection of Mendelian loci. Using genome-wide association results, we demonstrate that common variants associated with complex diseases are enriched in the genes indicated by this "Mendelian code." Finally, we detect hundreds of comorbidity associations among Mendelian disorders, and we use probabilistic genetic modeling to demonstrate that Mendelian variants likely contribute nonadditively to the risk for a subset of complex diseases. Overall, this study illustrates a complementary approach for mapping complex disease loci and provides unique predictions concerning the etiologies of specific diseases.

FAVA: high-quality functional association networks inferred from scRNA-seq and proteomics data.

MOTIVATION: Protein networks are commonly used for understanding how proteins interact. However, they are typically biased by data availability, favoring well-studied proteins with more interactions. To uncover functions of understudied proteins, we must use data that are not affected by this literature bias, such as single-cell RNA-seq and proteomics. Due to data sparseness and redundancy, functional association analysis becomes complex. RESULTS: To address this, we have developed FAVA (Functional Associations using Variational Autoencoders), which compresses high-dimensional data into a low-dimensional space. FAVA infers networks from high-dimensional omics data with much higher accuracy than existing methods, across a diverse collection of real as well as simulated datasets. FAVA can process large datasets with over 0.5 million conditions and has predicted 4210 interactions between 1039 understudied proteins. Our findings showcase FAVA's capability to offer novel perspectives on protein interactions. FAVA functions within the scverse ecosystem, employing AnnData as its input source. AVAILABILITY AND IMPLEMENTATION: Source code, documentation, and tutorials for FAVA are accessible on GitHub at https://github.com/mikelkou/fava. FAVA can also be installed and used via pip/PyPI as well as via the scverse ecosystem https://github.com/scverse/ecosystem-packages/tree/main/packages/favapy.

Pharos 2023: an integrated resource for the understudied human proteome.

The Illuminating the Druggable Genome (IDG) project aims to improve our understanding of understudied proteins and our ability to study them in the context of disease biology by perturbing them with small molecules, biologics, or other therapeutic modalities. Two main products from the IDG effort are the Target Central Resource Database (TCRD) (http://juniper.health.unm.edu/tcrd/), which curates and aggregates information, and Pharos (https://pharos.nih.gov/), a web interface for fusers to extract and visualize data from TCRD. Since the 2021 release, TCRD/Pharos has focused on developing visualization and analysis tools that help reveal higher-level patterns in the underlying data. The current iterations of TCRD and Pharos enable users to perform enrichment calculations based on subsets of targets, diseases, or ligands and to create interactive heat maps and UpSet charts of many types of annotations. Using several examples, we show how to address disease biology and drug discovery questions through enrichment calculations and UpSet charts.

Metastatic Infiltration of Nervous Tissue and Periosteal Nerve Sprouting in Multiple Myeloma-Induced Bone Pain in Mice and Human.

Multiple myeloma (MM) is a neoplasia of B plasma cells that often induces bone pain. However, the mechanisms underlying myeloma-induced bone pain (MIBP) are mostly unknown. Using a syngeneic MM mouse model, we show that periosteal nerve sprouting of calcitonin gene-related peptide (CGRP+) and growth associated protein 43 (GAP43+) fibers occurs concurrent to the onset of nociception and its blockade provides transient pain relief. MM patient samples also showed increased periosteal innervation. Mechanistically, we investigated MM induced gene expression changes in the dorsal root ganglia (DRG) innervating the MM-bearing bone of male mice and found alterations in pathways associated with cell cycle, immune response and neuronal signaling. The MM transcriptional signature was consistent with metastatic MM infiltration to the DRG, a never-before described feature of the disease that we further demonstrated histologically. In the DRG, MM cells caused loss of vascularization and neuronal injury, which may contribute to late-stage MIBP. Interestingly, the transcriptional signature of a MM patient was consistent with MM cell infiltration to the DRG. Overall, our results suggest that MM induces a plethora of peripheral nervous system alterations that may contribute to the failure of current analgesics and suggest neuroprotective drugs as appropriate strategies to treat early onset MIBP.SIGNIFICANCE STATEMENT Multiple myeloma (MM) is a painful bone marrow cancer that significantly impairs the quality of life of the patients. Analgesic therapies for myeloma-induced bone pain (MIBP) are limited and often ineffective, and the mechanisms of MIBP remain unknown. In this manuscript, we describe cancer-induced periosteal nerve sprouting in a mouse model of MIBP, where we also encounter metastasis to the dorsal root ganglia (DRG), a never-before described feature of the disease. Concomitant to myeloma infiltration, the lumbar DRGs presented blood vessel damage and transcriptional alterations, which may mediate MIBP. Explorative studies on human tissue support our preclinical findings. Understanding the mechanisms of MIBP is crucial to develop targeted analgesic with better efficacy and fewer side effects for this patient population.

Opportunities and barriers in omics-based biomarker discovery for steatotic liver diseases.

The rising prevalence of liver diseases related to obesity and excessive use of alcohol is fuelling an increasing demand for accurate biomarkers aimed at community screening, diagnosis of steatohepatitis and significant fibrosis, monitoring, prognostication and prediction of treatment efficacy. Breakthroughs in omics methodologies and the power of bioinformatics have created an excellent opportunity to apply technological advances to clinical needs, for instance in the development of precision biomarkers for personalised medicine. Via omics technologies, biological processes from the genes to circulating protein, as well as the microbiome - including bacteria, viruses and fungi, can be investigated on an axis. However, there are important barriers to omics-based biomarker discovery and validation, including the use of semi-quantitative measurements from untargeted platforms, which may exhibit high analytical, inter- and intra-individual variance. Standardising methods and the need to validate them across diverse populations presents a challenge, partly due to disease complexity and the dynamic nature of biomarker expression at different disease stages. Lack of validity causes lost opportunities when studies fail to provide the knowledge needed for regulatory approvals, all of which contributes to a delayed translation of these discoveries into clinical practice. While no omics-based biomarkers have matured to clinical implementation, the extent of data generated has enabled the hypothesis-free discovery of a plethora of candidate biomarkers that warrant further validation. To explore the many opportunities of omics technologies, hepatologists need detailed knowledge of commonalities and differences between the various omics layers, and both the barriers to and advantages of these approaches.

S1000: a better taxonomic name corpus for biomedical information extraction.

MOTIVATION: The recognition of mentions of species names in text is a critically important task for biomedical text mining. While deep learning-based methods have made great advances in many named entity recognition tasks, results for species name recognition remain poor. We hypothesize that this is primarily due to the lack of appropriate corpora. RESULTS: We introduce the S1000 corpus, a comprehensive manual re-annotation and extension of the S800 corpus. We demonstrate that S1000 makes highly accurate recognition of species names possible (F-score =93.1%), both for deep learning and dictionary-based methods. AVAILABILITY AND IMPLEMENTATION: All resources introduced in this study are available under open licenses from https://jensenlab.org/resources/s1000/. The webpage contains links to a Zenodo project and three GitHub repositories associated with the study.

Identifying the genes impacted by cell proliferation in proteomics and transcriptomics studies.

Hypothesis-free high-throughput profiling allows relative quantification of thousands of proteins or transcripts across samples and thereby identification of differentially expressed genes. It is used in many biological contexts to characterize differences between cell lines and tissues, identify drug mode of action or drivers of drug resistance, among others. Changes in gene expression can also be due to confounding factors that were not accounted for in the experimental plan, such as change in cell proliferation. We combined the analysis of 1,076 and 1,040 cell lines in five proteomics and three transcriptomics data sets to identify 157 genes that correlate with cell proliferation rates. These include actors in DNA replication and mitosis, and genes periodically expressed during the cell cycle. This signature of cell proliferation is a valuable resource when analyzing high-throughput data showing changes in proliferation across conditions. We show how to use this resource to help in interpretation of in vitro drug screens and tumor samples. It informs on differences of cell proliferation rates between conditions where such information is not directly available. The signature genes also highlight which hits in a screen may be due to proliferation changes; this can either contribute to biological interpretation or help focus on experiment-specific regulation events otherwise buried in the statistical analysis.

TCRD and Pharos 2021: mining the human proteome for disease biology.

In 2014, the National Institutes of Health (NIH) initiated the Illuminating the Druggable Genome (IDG) program to identify and improve our understanding of poorly characterized proteins that can potentially be modulated using small molecules or biologics. Two resources produced from these efforts are: The Target Central Resource Database (TCRD) (http://juniper.health.unm.edu/tcrd/) and Pharos (https://pharos.nih.gov/), a web interface to browse the TCRD. The ultimate goal of these resources is to highlight and facilitate research into currently understudied proteins, by aggregating a multitude of data sources, and ranking targets based on the amount of data available, and presenting data in machine learning ready format. Since the 2017 release, both TCRD and Pharos have produced two major releases, which have incorporated or expanded an additional 25 data sources. Recently incorporated data types include human and viral-human protein-protein interactions, protein-disease and protein-phenotype associations, and drug-induced gene signatures, among others. These aggregated data have enabled us to generate new visualizations and content sections in Pharos, in order to empower users to find new areas of study in the druggable genome.

TIGA: target illumination GWAS analytics.

MOTIVATION: Genome-wide association studies can reveal important genotype-phenotype associations; however, data quality and interpretability issues must be addressed. For drug discovery scientists seeking to prioritize targets based on the available evidence, these issues go beyond the single study. RESULTS: Here, we describe rational ranking, filtering and interpretation of inferred gene-trait associations and data aggregation across studies by leveraging existing curation and harmonization efforts. Each gene-trait association is evaluated for confidence, with scores derived solely from aggregated statistics, linking a protein-coding gene and phenotype. We propose a method for assessing confidence in gene-trait associations from evidence aggregated across studies, including a bibliometric assessment of scientific consensus based on the iCite relative citation ratio, and meanRank scores, to aggregate multivariate evidence.This method, intended for drug target hypothesis generation, scoring and ranking, has been implemented as an analytical pipeline, available as open source, with public datasets of results, and a web application designed for usability by drug discovery scientists. AVAILABILITY AND IMPLEMENTATION: Web application, datasets and source code via https://unmtid-shinyapps.net/tiga/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CoCoScore: context-aware co-occurrence scoring for text mining applications using distant supervision.

MOTIVATION: Information extraction by mining the scientific literature is key to uncovering relations between biomedical entities. Most existing approaches based on natural language processing extract relations from single sentence-level co-mentions, ignoring co-occurrence statistics over the whole corpus. Existing approaches counting entity co-occurrences ignore the textual context of each co-occurrence. RESULTS: We propose a novel corpus-wide co-occurrence scoring approach to relation extraction that takes the textual context of each co-mention into account. Our method, called CoCoScore, scores the certainty of stating an association for each sentence that co-mentions two entities. CoCoScore is trained using distant supervision based on a gold-standard set of associations between entities of interest. Instead of requiring a manually annotated training corpus, co-mentions are labeled as positives/negatives according to their presence/absence in the gold standard. We show that CoCoScore outperforms previous approaches in identifying human disease-gene and tissue-gene associations as well as in identifying physical and functional protein-protein associations in different species. CoCoScore is a versatile text mining tool to uncover pairwise associations via co-occurrence mining, within and beyond biomedical applications. AVAILABILITY AND IMPLEMENTATION: CoCoScore is available at: https://github.com/JungeAlexander/cocoscore. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Yield and Integrity of RNA from Brain Samples are Largely Unaffected by Pre-analytical Procedures.

Gene expression studies are reported to be influenced by pre-analytical factors that can compromise RNA yield and integrity, which in turn may confound the experimental findings. Here we investigate the impact of four pre-analytical factors on brain-derived RNA: time-before-collection, tissue specimen size, tissue collection method, and RNA isolation method. We report no significant differences in RNA yield or integrity between 20 mg and 60 mg tissue samples collected in either liquid nitrogen or the RNAlater stabilizing solution. Isolation of RNA employing the TRIzol reagent resulted in a higher yield compared to isolation via the QIAcube kit while the latter resulted in RNA of slightly better integrity. Keeping brain tissue samples at room temperature for up to 160 min prior to collection and isolation of RNA resulted in no significant difference in yield or integrity. Our findings have significant practical and financial consequences for clinical genomic departments and other laboratory settings performing large-scale routine RNA expression analysis of brain samples.

Inferring disease-associated long non-coding RNAs using genome-wide tissue expression profiles.

MOTIVATION: Long non-coding RNAs (lncRNAs) are important regulators in wide variety of biological processes, which are linked to many diseases. Compared to protein-coding genes (PCGs), the association between diseases and lncRNAs is still not well studied. Thus, inferring disease-associated lncRNAs on a genome-wide scale has become imperative. RESULTS: In this study, we propose a machine learning-based method, DislncRF, which infers disease-associated lncRNAs on a genome-wide scale based on tissue expression profiles. DislncRF uses random forest models trained on expression profiles of known disease-associated PCGs across human tissues to extract general patterns between expression profiles and diseases. These models are then applied to score associations between lncRNAs and diseases. DislncRF was benchmarked against a gold standard dataset and compared to other methods. The results show that DislncRF yields promising performance and outperforms the existing methods. The utility of DislncRF is further substantiated on two diseases in which we find that top scoring candidates are supported by literature or independent datasets. AVAILABILITY AND IMPLEMENTATION: https://github.com/xypan1232/DislncRF. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ProtFus: A Comprehensive Method Characterizing Protein-Protein Interactions of Fusion Proteins.

Tailored therapy aims to cure cancer patients effectively and safely, based on the complex interactions between patients' genomic features, disease pathology and drug metabolism. Thus, the continual increase in scientific literature drives the need for efficient methods of data mining to improve the extraction of useful information from texts based on patients' genomic features. An important application of text mining to tailored therapy in cancer encompasses the use of mutations and cancer fusion genes as moieties that change patients' cellular networks to develop cancer, and also affect drug metabolism. Fusion proteins, which are derived from the slippage of two parental genes, are produced in cancer by chromosomal aberrations and trans-splicing. Given that the two parental proteins for predicted fusion proteins are known, we used our previously developed method for identifying chimeric protein-protein interactions (ChiPPIs) associated with the fusion proteins. Here, we present a validation approach that receives fusion proteins of interest, predicts their cellular network alterations by ChiPPI and validates them by our new method, ProtFus, using an online literature search. This process resulted in a set of 358 fusion proteins and their corresponding protein interactions, as a training set for a Naïve Bayes classifier, to identify predicted fusion proteins that have reliable evidence in the literature and that were confirmed experimentally. Next, for a test group of 1817 fusion proteins, we were able to identify from the literature 2908 PPIs in total, across 18 cancer types. The described method, ProtFus, can be used for screening the literature to identify unique cases of fusion proteins and their PPIs, as means of studying alterations of protein networks in cancers. Availability: http://protfus.md.biu.ac.il/.

Standardized benchmarking in the quest for orthologs.

Achieving high accuracy in orthology inference is essential for many comparative, evolutionary and functional genomic analyses, yet the true evolutionary history of genes is generally unknown and orthologs are used for very different applications across phyla, requiring different precision-recall trade-offs. As a result, it is difficult to assess the performance of orthology inference methods. Here, we present a community effort to establish standards and an automated web-based service to facilitate orthology benchmarking. Using this service, we characterize 15 well-established inference methods and resources on a battery of 20 different benchmarks. Standardized benchmarking provides a way for users to identify the most effective methods for the problem at hand, sets a minimum requirement for new tools and resources, and guides the development of more accurate orthology inference methods.

A comprehensive and quantitative comparison of text-mining in 15 million full-text articles versus their corresponding abstracts.

Across academia and industry, text mining has become a popular strategy for keeping up with the rapid growth of the scientific literature. Text mining of the scientific literature has mostly been carried out on collections of abstracts, due to their availability. Here we present an analysis of 15 million English scientific full-text articles published during the period 1823-2016. We describe the development in article length and publication sub-topics during these nearly 250 years. We showcase the potential of text mining by extracting published protein-protein, disease-gene, and protein subcellular associations using a named entity recognition system, and quantitatively report on their accuracy using gold standard benchmark data sets. We subsequently compare the findings to corresponding results obtained on 16.5 million abstracts included in MEDLINE and show that text mining of full-text articles consistently outperforms using abstracts only.

Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse.

The rich fossil record of equids has made them a model for evolutionary processes. Here we present a 1.12-times coverage draft genome from a horse bone recovered from permafrost dated to approximately 560-780 thousand years before present (kyr BP). Our data represent the oldest full genome sequence determined so far by almost an order of magnitude. For comparison, we sequenced the genome of a Late Pleistocene horse (43 kyr BP), and modern genomes of five domestic horse breeds (Equus ferus caballus), a Przewalski's horse (E. f. przewalskii) and a donkey (E. asinus). Our analyses suggest that the Equus lineage giving rise to all contemporary horses, zebras and donkeys originated 4.0-4.5 million years before present (Myr BP), twice the conventionally accepted time to the most recent common ancestor of the genus Equus. We also find that horse population size fluctuated multiple times over the past 2 Myr, particularly during periods of severe climatic changes. We estimate that the Przewalski's and domestic horse populations diverged 38-72 kyr BP, and find no evidence of recent admixture between the domestic horse breeds and the Przewalski's horse investigated. This supports the contention that Przewalski's horses represent the last surviving wild horse population. We find similar levels of genetic variation among Przewalski's and domestic populations, indicating that the former are genetically viable and worthy of conservation efforts. We also find evidence for continuous selection on the immune system and olfaction throughout horse evolution. Finally, we identify 29 genomic regions among horse breeds that deviate from neutrality and show low levels of genetic variation compared to the Przewalski's horse. Such regions could correspond to loci selected early during domestication.

Accurate Quantification of Site-specific Acetylation Stoichiometry Reveals the Impact of Sirtuin Deacetylase CobB on the E. coli Acetylome.

Lysine acetylation is a protein posttranslational modification (PTM) that occurs on thousands of lysine residues in diverse organisms from bacteria to humans. Accurate measurement of acetylation stoichiometry on a proteome-wide scale remains challenging. Most methods employ a comparison of chemically acetylated peptides to native acetylated peptides, however, the potentially large differences in abundance between these peptides presents a challenge for accurate quantification. Stable isotope labeling by amino acids in cell culture (SILAC)-based mass spectrometry (MS) is one of the most widely used quantitative proteomic methods. Here we show that serial dilution of SILAC-labeled peptides (SD-SILAC) can be used to identify accurately quantified peptides and to estimate the quantification error rate. We applied SD-SILAC to determine absolute acetylation stoichiometry in exponentially-growing and stationary-phase wild-type and Sirtuin deacetylase CobB-deficient cells. To further analyze CobB-regulated sites under conditions of globally increased or decreased acetylation, we measured stoichiometry in phophotransacetylase (ptaΔ) and acetate kinase (ackAΔ) mutant strains in the presence and absence of the Sirtuin inhibitor nicotinamide. We measured acetylation stoichiometry at 3,669 unique sites and found that the vast majority of acetylation occurred at a low stoichiometry. Manipulations that cause increased nonenzymatic acetylation by acetyl-phosphate (AcP), such as stationary-phase arrest and deletion of ackA, resulted in globally increased acetylation stoichiometry. Comparison to relative quantification under the same conditions validated our stoichiometry estimates at hundreds of sites, demonstrating the accuracy of our method. Similar to Sirtuin deacetylase 3 (SIRT3) in mitochondria, CobB suppressed acetylation to lower than median stoichiometry in WT, ptaΔ, and ackAΔ cells. Together, our results provide a detailed view of acetylation stoichiometry in E. coli and suggest an evolutionarily conserved function of Sirtuin deacetylases in suppressing low stoichiometry acetylation.

Pharos: Collating protein information to shed light on the druggable genome.

The 'druggable genome' encompasses several protein families, but only a subset of targets within them have attracted significant research attention and thus have information about them publicly available. The Illuminating the Druggable Genome (IDG) program was initiated in 2014, has the goal of developing experimental techniques and a Knowledge Management Center (KMC) that would collect and organize information about protein targets from four families, representing the most common druggable targets with an emphasis on understudied proteins. Here, we describe two resources developed by the KMC: the Target Central Resource Database (TCRD) which collates many heterogeneous gene/protein datasets and Pharos (https://pharos.nih.gov), a multimodal web interface that presents the data from TCRD. We briefly describe the types and sources of data considered by the KMC and then highlight features of the Pharos interface designed to enable intuitive access to the IDG knowledgebase. The aim of Pharos is to encourage 'serendipitous browsing', whereby related, relevant information is made easily discoverable. We conclude by describing two use cases that highlight the utility of Pharos and TCRD.