An effective sequence-alignment-free superpositioning of pairwise or multiple structures with missing data

Background

Superpositioning is an important problem in structural biology. Determining an optimal superposition requires a one-to-one correspondence between the atoms of two proteins structures. However, in practice, some atoms are missing from their original structures. Current superposition implementations address the missing data crudely by ignoring such atoms from their structures.

Results

In this paper, we propose an effective method for superpositioning pairwise and multiple structures without sequence alignment. It is a two-stage procedure including data reduction and data registration.

Conclusions

Numerical experiments demonstrated that our method is effective and efficient. The code package of protein structure superposition method for addressing the cases with missing data is implemented by MATLAB, and it is freely available from: http://sourceforge.net/projects/pssm123/files/?source=navbar

     

Nasal gene expression differentiates COPD from controls and overlaps bronchial gene expression

Background

Nasal gene expression profiling is a promising method to characterize COPD non-invasively. We aimed to identify a nasal gene expression profile to distinguish COPD patients from healthy controls. We investigated whether this COPD-associated gene expression profile in nasal epithelium is comparable with the profile observed in bronchial epithelium.

Methods

Genome wide gene expression analysis was performed on nasal epithelial brushes of 31 severe COPD patients and 22 controls, all current smokers, using Affymetrix Human Gene 1.0 ST Arrays. We repeated the gene expression analysis on bronchial epithelial brushes in 2 independent cohorts of mild-to-moderate COPD patients and controls.

Results

In nasal epithelium, 135 genes were significantly differentially expressed between severe COPD patients and controls, 21 being up- and 114 downregulated in COPD (false discovery rate?<?0.01). Gene Set Enrichment Analysis (GSEA) showed significant concordant enrichment of COPD-associated nasal and bronchial gene expression in both independent cohorts (FDR_GSEA <?0.001).

Conclusion

We identified a nasal gene expression profile that differentiates severe COPD patients from controls. Of interest, part of the nasal gene expression changes in COPD mimics differentially expressed genes in the bronchus. These findings indicate that nasal gene expression profiling is potentially useful as a non-invasive biomarker in COPD.

Trial registration

ClinicalTrials.gov registration number NCT01351792 (registration date May 10, 2011), ClinicalTrials.gov registration number NCT00848406 (registration date February 19, 2009), ClinicalTrials.gov registration number NCT00807469 (registration date December 11, 2008).

     

PHYSICO2: an UNIX based standalone procedure for computation of physicochemical,window-dependent and substitution based evolutionary properties of protein sequences along with automated block preparation tool,version 2

AvailabilityPHYSICO2: is freely available at http://sourceforge.net/projects/physico2/ along with its documentation at https://sourceforge.net/projects/physico2/files/Documentation.pdf/download for all users.     

Concept recognition for extracting protein interaction relations from biomedical text

Background:

Reliable information extraction applications have been a long sought goal of the biomedical text mining community, a goal that if reached would provide valuable tools to benchside biologists in their increasingly difficult task of assimilating the knowledge contained in the biomedical literature. We present an integrated approach to concept recognition in biomedical text. Concept recognition provides key information that has been largely missing from previous biomedical information extraction efforts, namely direct links to well defined knowledge resources that explicitly cement the concept's semantics. The BioCreative II tasks discussed in this special issue have provided a unique opportunity to demonstrate the effectiveness of concept recognition in the field of biomedical language processing.

Results:

Through the modular construction of a protein interaction relation extraction system, we present several use cases of concept recognition in biomedical text, and relate these use cases to potential uses by the benchside biologist.

Conclusion:

Current information extraction technologies are approaching performance standards at which concept recognition can begin to deliver high quality data to the benchside biologist. Our system is available as part of the BioCreative Meta-Server project and on the internet http://bionlp.sourceforge.net.

     

Larval host records of butterflies in Japan

Using Japanese literature, we created a consolidated list of host records of butterflies in Japan. The list used the host records described in eight major illustrated reference books, two checklists, and 14 other pieces of literature. The presence of larvae on plants, the observation of larvae eating plants or insects in the field were considered as host records. We collected all species recorded in Japan. Scientific, family, and Japanese names of butterflies were consolidated using the BINRAN database (http://binran.lepimages.jp/). Scientific and Japanese names of host plants were based on the YList database (http://ylist.info/). If scientific names of host plants were not found in YList, we used scientific names based on The Plant List (http://www.theplantlist.org/). Family names of host plants were based on the Catalogue of Life database (http://www.catalogueoflife.org/). Scientific, family, and Japanese names of host insects were based on the MOKUROKU database (http://konchudb.agr.agr.kyushu-u.ac.jp/mokuroku/) for Hymenoptera and the catalogue of the Paraneoptera of Japan published by the Entomological Society of Japan for Hemiptera. We also provided the references of each host record and the original names described in the referred literature. Two datasets, HostDB and ReferenceDB, were created to include 3600 records of butterfly larval hosts in Japan, along with scientific and Japanese names of each species and a literature list. These datasets will be useful for basic and applied biological studies of butterflies. Data files are stored in the Ecological Research Data Archives (http://db.cger.nies.go.jp/JaLTER/ER_DataPapers/) and available from http://hostbj.lepumus.net/. These datasets are published under the Creative Commons License Attribution-ShareAlike 4.0 (CC BY-SA, https://creativecommons.org/licenses/by-sa/4.0/).     

UBCG: Up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction

Genome-based phylogeny plays a central role in the future taxonomy and phylogenetics of Bacteria and Archaea by replacing 16S rRNA gene phylogeny. The concatenated core gene alignments are frequently used for such a purpose. The bacterial core genes are defined as single-copy, homologous genes that are present in most of the known bacterial species. There have been several studies describing such a gene set, but the number of species considered was rather small. Here we present the up-to-date bacterial core gene set, named UBCG, and software suites to accommodate necessary steps to generate and evaluate phylogenetic trees. The method was successfully used to infer phylogenomic relationship of Escherichia and related taxa and can be used for the set of genomes at any taxonomic ranks of Bacteria. The UBCG pipeline and file viewer are freely available at https://www.ezbiocloud.net/tools/ubcg and https://www.ezbiocloud.net/tools/ubcg_viewer, respectively.     

PyPedia: using the wiki paradigm as crowd sourcing environment for bioinformatics protocols

Background

Today researchers can choose from many bioinformatics protocols for all types of life sciences research, computational environments and coding languages. Although the majority of these are open source, few of them possess all virtues to maximize reuse and promote reproducible science. Wikipedia has proven a great tool to disseminate information and enhance collaboration between users with varying expertise and background to author qualitative content via crowdsourcing. However, it remains an open question whether the wiki paradigm can be applied to bioinformatics protocols.

Results

We piloted PyPedia, a wiki where each article is both implementation and documentation of a bioinformatics computational protocol in the python language. Hyperlinks within the wiki can be used to compose complex workflows and induce reuse. A RESTful API enables code execution outside the wiki. Initial content of PyPedia contains articles for population statistics, bioinformatics format conversions and genotype imputation. Use of the easy to learn wiki syntax effectively lowers the barriers to bring expert programmers and less computer savvy researchers on the same page.

Conclusions

PyPedia demonstrates how wiki can provide a collaborative development, sharing and even execution environment for biologists and bioinformaticians that complement existing resources, useful for local and multi-center research teams.

Availability

PyPedia is available online at: http://www.pypedia.com. The source code and installation instructions are available at: https://github.com/kantale/PyPedia_server. The PyPedia python library is available at: https://github.com/kantale/pypedia. PyPedia is open-source, available under the BSD 2-Clause License.

     

Advances in transgenic animal models and techniques

On May 11th and 12th 2017 was held in Nantes, France, the international meeting “Advances in transgenic animal models and techniques” (http://www.trm.univ-nantes.fr/). This biennial meeting is the fifth one of its kind to be organized by the Transgenic Rats ImmunoPhenomic (TRIP) Nantes facility (http://www.tgr.nantes.inserm.fr/). The meeting was supported by private companies (SONIDEL, Scionics computer innovation, New England Biolabs, MERCK, genOway, Journal Disease Models and Mechanisms) and by public institutions (International Society for Transgenic Technology, University of Nantes, INSERM UMR 1064, SFR François Bonamy, CNRS, Région Pays de la Loire, Biogenouest, TEFOR infrastructure, ITUN, IHU-CESTI and DHU-Oncogeffe and Labex IGO). Around 100 participants, from France but also from different European countries, Japan and USA, attended the meeting.     

Mating and aggressive behaviour of <Emphasis Type="Italic">Brachyrhaphis olomina</Emphasis> (Cyprinodontiformes: Poeciliidae)

Despite the increasing interest in the use of intromittent male genitalia and coercive mating behaviour in poeciliids, detailed studies of the mating behaviour of most species in this family are lacking. We describe here the mating and aggressive behaviours of Brachyrhaphis olomina, and correlate them with the condition of the female’s ovum and embryos (immature, mature and pregnant). B. olomina performed a wide range of aggressive (sidle spread, tail beating, coordinate) and mating behaviours (approximation, touch, lateral display, touch-lateral display). Some behaviours (e.g. tail beating) are shared with other poeciliids, but two sexual behaviours (touch and lateral display) and one aggressive (coordinate) behaviour may be unique to B. olomina and were not reported in a previous study. Differences in male behaviour when paired with a female with mature ovum (more mating displays, no agonistic movements) suggest that males detect the female’s reproductive condition from some distance. The distinctive nature of mating behaviour in B. olomina highlights the importance of studying different species to have a better understanding of the evolution of mating and aggressive behaviours in poeciliids. Digital video images related to the article are available at http://www.momo-p.com/showdetail-e.php?movieid=momo170720bo01a, http://www.momo-p.com/showdetail-e.php?movieid=momo170720bo02a and http://www.momo-p.com/showdetail-e.php?movieid=momo170720bo03a.     

Long-term survival in patients treated with ruxolitinib for myelofibrosis: COMFORT-I and -II pooled analyses

Background

Myelofibrosis (MF) is associated with a variety of burdensome symptoms and reduced survival compared with age-/sex-matched controls. This analysis evaluated the long-term survival benefit with ruxolitinib, a Janus kinase (JAK)1/JAK2 inhibitor, in patients with intermediate-2 (int-2) or high-risk MF.

Methods

This was an exploratory analysis of 5-year data pooled from the phase 3 COMFORT-I and -II trials. In both trials, patients could cross over to ruxolitinib from the control group (COMFORT-I, placebo; COMFORT-II, best available therapy). All continuing patients in the control groups crossed over to ruxolitinib by the 3-year follow-up. Overall survival (OS; a secondary endpoint in both trials) was evaluated using pooled intent-to-treat data from patients randomized to ruxolitinib or the control groups. OS was also evaluated in subgroups stratified by baseline anemia and transfusion status at week 24.

Results

A total of 528 patients were included in this analysis; 301 were originally randomized to ruxolitinib (COMFORT-I, n?=?155; COMFORT-II, n?=?146) and 227 to control (n?=?154 and n?=?73, respectively). The risk of death was reduced by 30% among patients randomized to ruxolitinib compared with patients in the control group (median OS, 5.3 vs 3.8 years, respectively; hazard ratio [HR], 0.70 [95% CI, 0.54–0.91]; P?=?0.0065). After correcting for crossover using a rank-preserving structural failure time (RPSFT) method, the OS advantage was more pronounced for patients who were originally randomized to ruxolitinib compared with patients who crossed over from control to ruxolitinib (median OS, 5.3 vs 2.3 years; HR [ruxolitinib vs RPSFT], 0.35 [95% CI, 0.23–0.59]). An analysis of OS censoring patients at the time of crossover also demonstrated that ruxolitinib prolonged OS compared with control (median OS, 5.3 vs 2.4 years; HR [ruxolitinib vs censored at crossover], 0.53 [95% CI, 0.36–0.78]; P?=?0.0013). The survival benefit with ruxolitinib was observed irrespective of baseline anemia status or transfusion requirements at week 24.

Conclusions

These findings support ruxolitinib treatment for patients with int-2 or high-risk MF, regardless of anemia or transfusion status. Further analyses will be important for exploring ruxolitinib earlier in the disease course to assess the effect on the natural history of MF.

Trial registration

ClinicalTrials.gov identifiers, NCT00952289 and NCT00934544.

     

<Emphasis Type="Italic">CellFateScout</Emphasis> – a bioinformatics tool for elucidating small molecule signaling pathways that drive cells in a specific direction

Background

Small molecule effects can be represented by active signaling pathways within functional networks. Identifying these can help to design new strategies to utilize known small molecules, e.g. to trigger specific cellular transformations or to reposition known drugs.

Results

We developed CellFateScout that uses the method of Latent Variables to turn differential high-throughput expression data and a functional network into a list of active signaling pathways. Applying it to Connectivity Map data, i.e., differential expression data describing small molecule effects, we then generated a Human Small Molecule Mechanisms Database. Finally, using a list of active signaling pathways as query, a similarity search can identify small molecules from the database that may trigger these pathways. We validated our approach systematically, using expression data of small molecule perturbations, yielding better predictions than popular bioinformatics tools.

Conclusions

CellFateScout can be used to select small molecules for their desired effects. The CellFateScout Cytoscape plugin, a tutorial and the Human Small Molecule Mechanisms Database are available at https://sourceforge.net/projects/cellfatescout/ under LGPLv2 license.

     

Web-based design and analysis tools for CRISPR base editing

Background

As a result of its simplicity and high efficiency, the CRISPR-Cas system has been widely used as a genome editing tool. Recently, CRISPR base editors, which consist of deactivated Cas9 (dCas9) or Cas9 nickase (nCas9) linked with a cytidine or a guanine deaminase, have been developed. Base editing tools will be very useful for gene correction because they can produce highly specific DNA substitutions without the introduction of any donor DNA, but dedicated web-based tools to facilitate the use of such tools have not yet been developed.

Results

We present two web tools for base editors, named BE-Designer and BE-Analyzer. BE-Designer provides all possible base editor target sequences in a given input DNA sequence with useful information including potential off-target sites. BE-Analyzer, a tool for assessing base editing outcomes from next generation sequencing (NGS) data, provides information about mutations in a table and interactive graphs. Furthermore, because the tool runs client-side, large amounts of targeted deep sequencing data (<?1?GB) do not need to be uploaded to a server, substantially reducing running time and increasing data security. BE-Designer and BE-Analyzer can be freely accessed at http://www.rgenome.net/be-designer/ and http://www.rgenome.net/be-analyzer/, respectively.

Conclusion

We develop two useful web tools to design target sequence (BE-Designer) and to analyze NGS data from experimental results (BE-Analyzer) for CRISPR base editors.

     

Prediction of linear B-cell epitopes

The use of antigenicity scales based on physicochemical properties and the sliding window method in combination with an averaging algorithm and subsequent search for the maximum value is the classical method for B-cell epitope prediction. However, recent studies have demonstrated that the best classical methods provide a poor correlation with experimental data. We review both classical and novel algorithms and present our own implementation of the algorithms. The AAPPred software is available at http://www.bioinf.ru/aappred/.     

HeurAA: Accurate and Fast Detection of Genetic Variations with a Novel Heuristic Amplicon Aligner Program for Next Generation Sequencing

Next generation sequencing (NGS) of PCR amplicons is a standard approach to detect genetic variations in personalized medicine such as cancer diagnostics. Computer programs used in the NGS community often miss insertions and deletions (indels) that constitute a large part of known human mutations. We have developed HeurAA, an open source, heuristic amplicon aligner program. We tested the program on simulated datasets as well as experimental data from multiplex sequencing of 40 amplicons in 12 oncogenes collected on a 454 Genome Sequencer from lung cancer cell lines. We found that HeurAA can accurately detect all indels, and is more than an order of magnitude faster than previous programs. HeurAA can compare reads and reference sequences up to several thousand base pairs in length, and it can evaluate data from complex mixtures containing reads of different gene-segments from different samples. HeurAA is written in C and Perl for Linux operating systems, the code and the documentation are available for research applications at http://sourceforge.net/projects/heuraa/