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.

     

The Year of the Wisent

Delving into European prehistory, two recent studies analyze ancient DNA from bison species depicted by our ancestors on the walls of their caves. The DNA tells a story of migrations driven by climate change but leaves some mystery clouding the genetic descent and climate preference of the still-extant wisent, otherwise known as the European bison.See research articles: https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-016-0317-7 http://www.nature.com/articles/ncomms13158     

The FOCUS,AFFINITY and EFFECTS trials studying the effect(s) of fluoxetine in patients with a recent stroke: statistical and health economic analysis plan for the trials and for the individual patient data meta-analysis

Background

Small trials have suggested that fluoxetine may improve neurological recovery from stroke. FOCUS, AFFINITY and EFFECTS are a family of investigator-led, multicentre, parallel group, randomised, placebo-controlled trials which aim to determine whether the routine administration of fluoxetine (20 mg daily) for six months after an acute stroke improves patients’ functional outcome.

Methods/Design

The core protocol for the three trials has been published (Mead et al., Trials 20:369, 2015). The trials include patients aged 18 years and older with a clinical diagnosis of stroke and persisting focal neurological deficits at randomisation 2–15 days after stroke onset. Patients are randomised centrally via each trials’ web-based randomisation system using a common minimisation algorithm. Patients are allocated fluoxetine 20 mg once daily or matching placebo capsules for six months. The primary outcome measure is the modified Rankin scale (mRS) at six months. Secondary outcomes include: living circumstances; the Stroke Impact Scale; EuroQol (EQ5D-5 L); the vitality subscale of the 36-Item Short Form Health Survey (SF36); diagnosis of depression; adherence to medication; serious adverse events including death and recurrent stroke; and resource use at six and 12 months and the mRS at 12 months.

Discussion

Minor variations have been tailored to the national setting in the UK (FOCUS), Australia, New Zealand and Vietnam (AFFINITY) and Sweden (EFFECTS). Each trial is run and funded independently and will report its own results. A prospectively planned individual patient data meta-analysis of all three trials will provide the most precise estimate of the overall effect and establish whether any effects differ between trials or subgroups. This statistical analysis plan describes the core analyses for all three trials and that for the individual patient data meta-analysis. Recruitment and follow-up in the FOCUS trial is expected to be completed by the end of 2018. AFFINITY and EFFECTS are likely to complete follow-up in 2020.

Trial registration

FOCUS: ISRCTN, ISRCTN83290762. Registered on 23 May 2012. EudraCT, 2011-005616-29. Registered on 3 February 2012.AFFINITY: Australian New Zealand Clinical Trials Registry, ACTRN12611000774921. Registered on 22 July 2011.EFFECTS: ISRCTN, ISRCTN13020412. Registered on 19 December 2014. Clinicaltrials.gov, NCT02683213. Registered on 2 February 2016. EudraCT, 2011-006130-16. Registered on 8 August 2014.

     

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/).     

Reliability of signal transmission in stochastic nerve axon equations

We introduce a method for computing probabilities for spontaneous activity and propagation failure of the action potential in spatially extended, conductance-based neuronal models subject to noise, based on statistical properties of the membrane potential. We compare different estimators with respect to the quality of detection, computational costs and robustness and propose the integral of the membrane potential along the axon as an appropriate estimator to detect both spontaneous activity and propagation failure. Performing a model reduction we achieve a simplified analytical expression based on the linearization at the resting potential (resp. the traveling action potential). This allows to approximate the probabilities for spontaneous activity and propagation failure in terms of (classical) hitting probabilities of one-dimensional linear stochastic differential equations. The quality of the approximation with respect to the noise amplitude is discussed and illustrated with numerical results for the spatially extended Hodgkin-Huxley equations. Python simulation code is supplied on GitHub under the link https://github.com/deristnochda/Hodgkin-Huxley-SPDE.     

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).

     

Paternal care behaviors of Japanese giant salamander <Emphasis Type="Italic">Andrias japonicus</Emphasis> in natural populations

Parental care among salamanders is typically provided by females. A rare case of parental care by male salamanders appears to occur in Cryptobranchidae. Yet, paternal behaviors have rarely been reported from natural populations of any Cryptobranchid salamanders, and their adaptive significance is poorly understood. The present study aimed to examine paternal care behaviors in a fully aquatic Japanese giant salamander (Andrias japonicus) in situ. At the beginning of the summer breeding season, large males, called den-masters, occupy burrows along stream banks for breeding and nesting. We videotaped post-breeding behaviors of two den-masters that stayed with the eggs, one in a natural and the other in an artificial nest in natural streams. We identified three behaviors, tail fanning, agitating and egg eating, to be parental care. Tail fanning provides oxygenated water for the eggs. We found that the den-master in the artificial nest, where dissolved oxygen level was lower, displayed tail fanning more frequently. Agitating the eggs with its head and body likely prevents yolk adhesions. The den-masters selectively ate whiter eggs that appeared to be dead or infected with water mold. This behavior, which we termed hygienic filial cannibalism, likely prevents water mold from spreading over healthy eggs. Digital video images relating to this article are available at http://www.momo-p.com/showdetail-e.php?movieid=momo140906aj01a, http://www.momo-p.com/showdetail-e.php?movieid=momo140906aj02a, http://www.momo-p.com/showdetail-e.php?movieid=momo140906aj03a and http://www.momo-p.com/showdetail-e.php?movieid=momo140906aj04a.     

Analyses of density-dependent effects are needed to understand how and when <Emphasis Type="Italic">Wolbachia</Emphasis> can control dengue vectors

Releases of Wolbachia-infected mosquitoes have been shown to be an effective method of controlling Aedes aegypti, the main vector of dengue fever, in Australia. A study in BMC Biology from Penelope Hancock and others shows that incorporation of density-dependent effects into population models can provide major improvements in understanding how and when the infected populations can become established.See research article: https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-016-0319-5.     

MM2S: personalized diagnosis of medulloblastoma patients and model systems

Background

Medulloblastoma (MB) is a highly malignant and heterogeneous brain tumour that is the most common cause of cancer-related deaths in children. Increasing availability of genomic data over the last decade had resulted in improvement of human subtype classification methods, and the parallel development of MB mouse models towards identification of subtype-specific disease origins and signaling pathways. Despite these advances, MB classification schemes remained inadequate for personalized prediction of MB subtypes for individual patient samples and across model systems. To address this issue, we developed the Medullo-Model to Subtypes ( MM2S ) classifier, a new method enabling classification of individual gene expression profiles from MB samples (patient samples, mouse models, and cell lines) against well-established molecular subtypes [Genomics 106:96-106, 2015]. We demonstrated the accuracy and flexibility of MM2S in the largest meta-analysis of human patients and mouse models to date. Here, we present a new functional package that provides an easy-to-use and fully documented implementation of the MM2S method, with additional functionalities that allow users to obtain graphical and tabular summaries of MB subtype predictions for single samples and across sample replicates. The flexibility of the MM2S package promotes incorporation of MB predictions into large Medulloblastoma-driven analysis pipelines, making this tool suitable for use by researchers.

Results

The MM2S package is applied in two case studies involving human primary patient samples, as well as sample replicates of the GTML mouse model. We highlight functions that are of use for species-specific MB classification, across individual samples and sample replicates. We emphasize on the range of functions that can be used to derive both singular and meta-centric views of MB predictions, across samples and across MB subtypes.

Conclusions

Our MM2S package can be used to generate predictions without having to rely on an external web server or additional sources. Our open-source package facilitates and extends the MM2S algorithm in diverse computational and bioinformatics contexts. The package is available on CRAN, at the following URL: https://cran.r-project.org/web/packages/MM2S/, as well as on Github at the following URLs: https://github.com/DGendoo and https://github.com/bhklab.

     

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.     

The scent of your enemy is my friend? The acquisition of large carnivore scent by a smaller carnivore

Scent marking is critical to intraspecific communication in many mammal species, but little is known regarding its role in communication among different species. We used 4 years of motion-triggered video to document the use of scent marking areas—termed “community scrapes”—by pumas (Puma concolor) (http://www.momo-p.com/showdetail-e.php?movieid=momo160812pc01a) and other carnivore species. We found that gray foxes (Urocyon cinereoargenteus) routinely rubbed their cheeks on puma scrapes (http://www.momo-p.com/showdetail-e.php?movieid=momo160812uc01a), and tested a series of hypotheses to determine its function. We found that gray foxes selected puma scrapes over other objects, and cheek rubbing by foxes was also correlated with how recently a puma had visited the scrape, suggesting that foxes were intent upon accumulating fresh puma scent. Cheek rubbing by foxes was not correlated with their breeding season or with how recently another fox had visited the site. Finally we found a cascading pattern in the occurrence of pumas, coyotes (Canis latrans) and gray foxes at community scrapes, suggesting that gray foxes may use puma scent to deter predation. This is the first published study to find evidence of a subordinate species using the scent of a dominant species to communicate with heterospecifics. The behavioral cascade we found in scent marking patterns also suggests that scent marking could be a mechanism that impacts the distribution and abundance of species. Additional videos pertaining to this article include http://www.momo-p.com/showdetail-e.php?movieid=momo160812uc02a, and http://www.momo-p.com/showdetail-e.php?movieid=momo160812uc03a.     

Interspecific responses of wild African carnivores to odour of 3-mercapto-3-methylbutanol,a component of wildcat and leopard urine

The scent of 3-mercapto-3-methylbutanol (3-M-3-MB), a volatile component of leopard (Panthera pardus) and domestic cat (Felis silvestris catus) urine, released at about 10 ng/s from slow-release dispensers, elicited scent-marking from African civet (Civettictis civetta), small-spotted genet (Genetta genetta) and slender mongoose (Galerella sanguinea), as well as African wildcat (F. s. cafra). A female leopard was apparently repelled by the scent. The scent-marking and scent-rubbing by species other than African wildcats and leopards were unexpected and have important implications for the design of studies to investigate chemical communication between wild mammals and the use of camera traps to estimate animal numbers. Videos showing the behaviours referred to in this article are available at; http://www.momo-p.com/showdetail-e.php?movieid=momo161223fs01a; http://www.momo-p.com/showdetail-e.php?movieid=momo161223gs01a; http://www.momo-p.com/showdetail-e.php?movieid=momo161223gg01a.     

PureCLIP: capturing target-specific protein–RNA interaction footprints from single-nucleotide CLIP-seq data

The iCLIP and eCLIP techniques facilitate the detection of protein–RNA interaction sites at high resolution, based on diagnostic events at crosslink sites. However, previous methods do not explicitly model the specifics of iCLIP and eCLIP truncation patterns and possible biases. We developed PureCLIP (https://github.com/skrakau/PureCLIP), a hidden Markov model based approach, which simultaneously performs peak-calling and individual crosslink site detection. It explicitly incorporates a non-specific background signal and, for the first time, non-specific sequence biases. On both simulated and real data, PureCLIP is more accurate in calling crosslink sites than other state-of-the-art methods and has a higher agreement across replicates.     

Speed and accuracy improvement of higher-order epistasis detection on CUDA-enabled GPUs

The discovery of higher-order epistatic interactions is an important task in the field of genome wide association studies which allows for the identification of complex interaction patterns between multiple genetic markers. Some existing bruteforce approaches explore the whole space of k-interactions in an exhaustive manner resulting in almost intractable execution times. Computational cost can be reduced drastically by restricting the search space with suitable preprocessing filters which prune unpromising candidates. Other approaches mitigate the execution time by employing massively parallel accelerators in order to benefit from the vast computational resources of these architectures. In this paper, we combine a novel preprocessing filter, namely SingleMI, with massively parallel computation on modern GPUs to further accelerate epistasis discovery. Our implementation improves both the runtime and accuracy when compared to a previous GPU counterpart that employs mutual information clustering for prefiltering. SingleMI is open source software and publicly available at: https://github.com/sleeepyjack/singlemi/.     

A Cladist is a systematist who seeks a natural classification: some comments on Quinn (2017)

In response to Quinn (Biol Philos, 2017.  https://doi.org/10.1007/s10539-017-9577-z) we identify cladistics to be about natural classifications and their discovery and thereby propose to add an eighth cladistic definition to Quinn’s list, namely the systematist who seeks to discover natural classifications, regardless of their affiliation, theoretical or methodological justifications.     

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.     

miRTrace reveals the organismal origins of microRNA sequencing data

We present here miRTrace, the first algorithm to trace microRNA sequencing data back to their taxonomic origins. This is a challenge with profound implications for forensics, parasitology, food control, and research settings where cross-contamination can compromise results. miRTrace accurately (>?99%) assigns real and simulated data to 14 important animal and plant groups, sensitively detects parasitic infection in mammals, and discovers the primate origin of single cells. Applying our algorithm to over 700 public datasets, we find evidence that over 7% are cross-contaminated and present a novel solution to clean these computationally, even after sequencing has occurred. miRTrace is freely available at https://github.com/friedlanderlab/mirtrace.     

Metastasis as supra-cellular selection? A reply to Lean and Plutynski

In response to Germain (Biol Philos 27:785–810, 2012. doi: 10.1007/s10539-012-9334-2) argument that evolution by natural selection has a limited explanatory power in cancer, Lean and Plutynski (Biol Philos 31:39–57, 2016. doi: 10.1007/s10539-015-9511-1) have recently argued that many adaptations in cancer only make sense at the tumor level, and that cancer progression mirrors the major evolutionary transitions. While we agree that selection could potentially act at various levels of organization in cancers, we argue that tumor-level selection (MLS2) is unlikely to actually play a relevant role in our understanding of the somatic evolution of human cancers.