Integrative Multi-omics Module Network Inference with Lemon-Tree

Module network inference is an established statistical method to reconstruct co-expression modules and their upstream regulatory programs from integrated multi-omics datasets measuring the activity levels of various cellular components across different individuals, experimental conditions or time points of a dynamic process. We have developed Lemon-Tree, an open-source, platform-independent, modular, extensible software package implementing state-of-the-art ensemble methods for module network inference. We benchmarked Lemon-Tree using large-scale tumor datasets and showed that Lemon-Tree algorithms compare favorably with state-of-the-art module network inference software. We also analyzed a large dataset of somatic copy-number alterations and gene expression levels measured in glioblastoma samples from The Cancer Genome Atlas and found that Lemon-Tree correctly identifies known glioblastoma oncogenes and tumor suppressors as master regulators in the inferred module network. Novel candidate driver genes predicted by Lemon-Tree were validated using tumor pathway and survival analyses. Lemon-Tree is available from http://lemon-tree.googlecode.com under the GNU General Public License version 2.0.

This is a PLOS Computational Biology Software Article

     

Inference of Population Splits and Mixtures from Genome-Wide Allele Frequency Data

Many aspects of the historical relationships between populations in a species are reflected in genetic data. Inferring these relationships from genetic data, however, remains a challenging task. In this paper, we present a statistical model for inferring the patterns of population splits and mixtures in multiple populations. In our model, the sampled populations in a species are related to their common ancestor through a graph of ancestral populations. Using genome-wide allele frequency data and a Gaussian approximation to genetic drift, we infer the structure of this graph. We applied this method to a set of 55 human populations and a set of 82 dog breeds and wild canids. In both species, we show that a simple bifurcating tree does not fully describe the data; in contrast, we infer many migration events. While some of the migration events that we find have been detected previously, many have not. For example, in the human data, we infer that Cambodians trace approximately 16% of their ancestry to a population ancestral to other extant East Asian populations. In the dog data, we infer that both the boxer and basenji trace a considerable fraction of their ancestry (9% and 25%, respectively) to wolves subsequent to domestication and that East Asian toy breeds (the Shih Tzu and the Pekingese) result from admixture between modern toy breeds and “ancient” Asian breeds. Software implementing the model described here, called TreeMix, is available at http://treemix.googlecode.com.     

Draft genome sequences and description of Lactobacillus rhamnosus strains L31, L34, and L35

Lactobacillus rhamnosus is a facultative, lactic acid bacterium in the phylum Firmicutes. Lactobacillus spp. are generally considered beneficial, and specific strains of L. rhamnosus are validated probiotics. We describe the draft genomes of three L. rhamnosus strains (L31, L34, and L35) isolated from the feces of Thai breastfed infants, which exhibit anti-inflammatory properties in vitro. The three genomes range between 2.8 – 2.9 Mb, and contain approximately 2,700 protein coding genes.     

Inference of the Properties of the Recombination Process from Whole Bacterial Genomes

Patterns of linkage disequilibrium, homoplasy, and incompatibility are difficult to interpret because they depend on several factors, including the recombination process and the population structure. Here we introduce a novel model-based framework to infer recombination properties from such summary statistics in bacterial genomes. The underlying model is sequentially Markovian so that data can be simulated very efficiently, and we use approximate Bayesian computation techniques to infer parameters. As this does not require us to calculate the likelihood function, the model can be easily extended to investigate less probed aspects of recombination. In particular, we extend our model to account for the bias in the recombination process whereby closely related bacteria recombine more often with one another. We show that this model provides a good fit to a data set of Bacillus cereus genomes and estimate several recombination properties, including the rate of bias in recombination. All the methods described in this article are implemented in a software package that is freely available for download at http://code.google.com/p/clonalorigin/.     

CREST – Classification Resources for Environmental Sequence Tags

Sequencing of taxonomic or phylogenetic markers is becoming a fast and efficient method for studying environmental microbial communities. This has resulted in a steadily growing collection of marker sequences, most notably of the small-subunit (SSU) ribosomal RNA gene, and an increased understanding of microbial phylogeny, diversity and community composition patterns. However, to utilize these large datasets together with new sequencing technologies, a reliable and flexible system for taxonomic classification is critical. We developed CREST (Classification Resources for Environmental Sequence Tags), a set of resources and tools for generating and utilizing custom taxonomies and reference datasets for classification of environmental sequences. CREST uses an alignment-based classification method with the lowest common ancestor algorithm. It also uses explicit rank similarity criteria to reduce false positives and identify novel taxa. We implemented this method in a web server, a command line tool and the graphical user interfaced program MEGAN. Further, we provide the SSU rRNA reference database and taxonomy SilvaMod, derived from the publicly available SILVA SSURef, for classification of sequences from bacteria, archaea and eukaryotes. Using cross-validation and environmental datasets, we compared the performance of CREST and SilvaMod to the RDP Classifier. We also utilized Greengenes as a reference database, both with CREST and the RDP Classifier. These analyses indicate that CREST performs better than alignment-free methods with higher recall rate (sensitivity) as well as precision, and with the ability to accurately identify most sequences from novel taxa. Classification using SilvaMod performed better than with Greengenes, particularly when applied to environmental sequences. CREST is freely available under a GNU General Public License (v3) from http://apps.cbu.uib.no/crest and http://lcaclassifier.googlecode.com.     

The complex binding of PRDM9

A recent study investigates the in vitro DNA binding behavior of PRDM9, a zinc finger protein involved in the localization of recombination hotspots in mammals.Please see related research article: http://genomebiology.com/2013/14/4/R35     

Timber! Felling the loblolly pine genome

Conventional short read sequences derived from haploid DNA were extended into long super-reads enabling assembly of the massive 22 Gbp loblolly pine, Pinus taeda, genome.See related research http://genomebiology.com/2014/15/3/R59     

High quality draft genome sequence of the slightly halophilic bacterium Halomonas zhanjiangensis type strain JSM 078169T (DSM 21076T) from a sea urchin in southern China

Halomonas zhanjiangensis Chen et al. 2009 is a member of the genus Halomonas, family Halomonadaceae, class Gammaproteobacteria. Representatives of the genus Halomonas are a group of halophilic bacteria often isolated from salty environments. The type strain H. zhanjiangensis JSM 078169^T was isolated from a sea urchin (Hemicentrotus pulcherrimus) collected from the South China Sea. The genome of strain JSM 078169^T is the fourteenth sequenced genome in the genus Halomonas and the fifteenth in the family Halomonadaceae. The other thirteen genomes from the genus Halomonas are H. halocynthiae, H. venusta, H. alkaliphila, H. lutea, H. anticariensis, H. jeotgali, H. titanicae, H. desiderata, H. smyrnensis, H. salifodinae, H. boliviensis, H. elongata and H stevensii. Here, we describe the features of strain JSM 078169^T, together with the complete genome sequence and annotation from a culture of DSM 21076^T. The 4,060,520 bp long draft genome consists of 17 scaffolds with the 3,659 protein-coding and 80 RNA genes and is a part of Genomic Encyclopedia of Type Strains, Phase I: the one thousand microbial genomes (KMG) project.     

A Real-Time All-Atom Structural Search Engine for Proteins

Protein designers use a wide variety of software tools for de novo design, yet their repertoire still lacks a fast and interactive all-atom search engine. To solve this, we have built the Suns program: a real-time, atomic search engine integrated into the PyMOL molecular visualization system. Users build atomic-level structural search queries within PyMOL and receive a stream of search results aligned to their query within a few seconds. This instant feedback cycle enables a new “designability”-inspired approach to protein design where the designer searches for and interactively incorporates native-like fragments from proven protein structures. We demonstrate the use of Suns to interactively build protein motifs, tertiary interactions, and to identify scaffolds compatible with hot-spot residues. The official web site and installer are located at http://www.degradolab.org/suns/ and the source code is hosted at https://github.com/godotgildor/Suns (PyMOL plugin, BSD license), https://github.com/Gabriel439/suns-cmd (command line client, BSD license), and https://github.com/Gabriel439/suns-search (search engine server, GPLv2 license).

This is a PLOS Computational Biology Software Article

     

Birds do it,bees do it,worms and ciliates do it too: DNA methylation from unexpected corners of the tree of life

Studies describing intricate patterns of DNA methylation in nematode and ciliate are controversial due to the uncertainty of genomic evolutionary conservation of DNA methylation enzymes.See related research articles http://genomebiology.com/2012/13/10/R99 and http://genomebiology.com/2012/13/10/R100     

A database for human Y chromosome protein data

The human Y chromosome is the sex determining chromosome. The number of proteins associated with this chromosome is 196 and 107 of the 196 proteins have yet not been characterised. Here, we describe the analysis of these 107 proteins by computing various physicochemical properties using sequence and predicted structural data to elucidate molecular function. We present the derived data in the form a form a database made freely available for download, review, refinement and update.

Availability

http://puratham.googlepages.com/ or http://puratham.googlepages.com/ftpconnection     

PathVisio 3: An Extendable Pathway Analysis Toolbox

PathVisio is a commonly used pathway editor, visualization and analysis software. Biological pathways have been used by biologists for many years to describe the detailed steps in biological processes. Those powerful, visual representations help researchers to better understand, share and discuss knowledge. Since the first publication of PathVisio in 2008, the original paper was cited more than 170 times and PathVisio was used in many different biological studies. As an online editor PathVisio is also integrated in the community curated pathway database WikiPathways.Here we present the third version of PathVisio with the newest additions and improvements of the application. The core features of PathVisio are pathway drawing, advanced data visualization and pathway statistics. Additionally, PathVisio 3 introduces a new powerful extension systems that allows other developers to contribute additional functionality in form of plugins without changing the core application.PathVisio can be downloaded from http://www.pathvisio.org and in 2014 PathVisio 3 has been downloaded over 5,500 times. There are already more than 15 plugins available in the central plugin repository. PathVisio is a freely available, open-source tool published under the Apache 2.0 license (http://www.apache.org/licenses/LICENSE-2.0). It is implemented in Java and thus runs on all major operating systems. The code repository is available at http://svn.bigcat.unimaas.nl/pathvisio. The support mailing list for users is available on https://groups.google.com/forum/#!forum/wikipathways-discuss and for developers on https://groups.google.com/forum/#!forum/wikipathways-devel.

This is a PLOS Computational Biology software article.

     

SNES makes sense? Single-cell exome sequencing evolves

Technologies for single-cell sequencing are improving steadily. A recent study describes a new method for interrogating all coding sequences of the human genome at single-cell resolution.See related research by Leung et al., http://genomebiology.com/2015/16/1/55     

Genome analysis of Desulfotomaculum gibsoniae strain GrollT a highly versatile Gram-positive sulfate-reducing bacterium

Desulfotomaculum gibsoniae is a mesophilic member of the polyphyletic spore-forming genus Desulfotomaculum within the family Peptococcaceae. This bacterium was isolated from a freshwater ditch and is of interest because it can grow with a large variety of organic substrates, in particular several aromatic compounds, short-chain and medium-chain fatty acids, which are degraded completely to carbon dioxide coupled to the reduction of sulfate. It can grow autotrophically with H₂ + CO₂ and sulfate and slowly acetogenically with H₂ + CO_2, formate or methoxylated aromatic compounds in the absence of sulfate. It does not require any vitamins for growth. Here, we describe the features of D. gibsoniae strain Groll^T together with the genome sequence and annotation. The chromosome has 4,855,529 bp organized in one circular contig and is the largest genome of all sequenced Desulfotomaculum spp. to date. A total of 4,666 candidate protein-encoding genes and 96 RNA genes were identified. Genes of the acetyl-CoA pathway, possibly involved in heterotrophic growth and in CO₂ fixation during autotrophic growth, are present. The genome contains a large set of genes for the anaerobic transformation and degradation of aromatic compounds, which are lacking in the other sequenced Desulfotomaculum genomes.     

Genome sequence of the clover-nodulating Rhizobium leguminosarum bv. trifolii strain SRDI943.

Rhizobium leguminosarum bv. trifolii SRDI943 (strain syn. V2-2) is an aerobic, motile, Gram-negative, non-spore-forming rod that was isolated from a root nodule of Trifolium michelianum Savi cv. Paradana that had been grown in soil collected from a mixed pasture in Victoria, Australia. This isolate was found to have a broad clover host range but was sub-optimal for nitrogen fixation with T. subterraneum (fixing 20-54% of reference inoculant strain WSM1325) and was found to be totally ineffective with the clover species T. polymorphum and T. pratense. Here we describe the features of R. leguminosarum bv. trifolii strain SRDI943, together with genome sequence information and annotation. The 7,412,387 bp high-quality-draft genome is arranged into 5 scaffolds of 5 contigs, contains 7,317 protein-coding genes and 89 RNA-only encoding genes, and is one of 100 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.     

Genome sequence of the clover-nodulating Rhizobium leguminosarum bv. trifolii strain SRDI565.

Rhizobium leguminosarum bv. trifolii SRDI565 (syn. N8-J) is an aerobic, motile, Gram-negative, non-spore-forming rod. SRDI565 was isolated from a nodule recovered from the roots of the annual clover Trifolium subterraneum subsp. subterraneum grown in the greenhouse and inoculated with soil collected from New South Wales, Australia. SRDI565 has a broad host range for nodulation within the clover genus, however N₂-fixation is sub-optimal with some Trifolium species and ineffective with others. Here we describe the features of R. leguminosarum bv. trifolii strain SRDI565, together with genome sequence information and annotation. The 6,905,599 bp high-quality-draft genome is arranged into 7 scaffolds of 7 contigs, contains 6,750 protein-coding genes and 86 RNA-only encoding genes, and is one of 100 rhizobial genomes sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.     

Ascochyta rabiei: A threat to global chickpea production

The necrotrophic fungus Ascochyta rabiei causes Ascochyta blight (AB) disease in chickpea. A. rabiei infects all aerial parts of the plant, which results in severe yield loss. At present, AB disease occurs in most chickpea‐growing countries. Globally increased incidences of A. rabiei infection and the emergence of new aggressive isolates directed the interest of researchers toward understanding the evolution of pathogenic determinants in this fungus. In this review, we summarize the molecular and genetic studies of the pathogen along with approaches that are helping in combating the disease. Possible areas of future research are also suggested.Taxonomykingdom Mycota, phylum Ascomycota, class Dothideomycetes, subclass Coelomycetes, order Pleosporales, family Didymellaceae, genus Ascochyta, species rabiei. Primary host A. rabiei survives primarily on Cicer species.Disease symptoms A. rabiei infects aboveground parts of the plant including leaves, petioles, stems, pods, and seeds. The disease symptoms first appear as watersoaked lesions on the leaves and stems, which turn brown or dark brown. Early symptoms include small circular necrotic lesions visible on the leaves and oval brown lesions on the stem. At later stages of infection, the lesions may girdle the stem and the region above the girdle falls off. The disease severity increases at the reproductive stage and rounded lesions with concentric rings, due to asexual structures called pycnidia, appear on leaves, stems, and pods. The infected pod becomes blighted and often results in shrivelled and infected seeds.Disease management strategiesCrop failures may be avoided by judicious practices of integrated disease management based on the use of resistant or tolerant cultivars and growing chickpea in areas where conditions are least favourable for AB disease development. Use of healthy seeds free of A. rabiei, seed treatments with fungicides, and proper destruction of diseased stubbles can also reduce the fungal inoculum load. Crop rotation with nonhost crops is critical for controlling the disease. Planting moderately resistant cultivars and prudent application of fungicides is also a way to combat AB disease. However, the scarcity of AB‐resistant accessions and the continuous evolution of the pathogen challenges the disease management process.Useful websites https://www.ndsu.edu/pubweb/pulse‐info/resourcespdf/Ascochyta%20blight%20of%20chickpea.pdf https://saskpulse.com/files/newsletters/180531_ascochyta_in_chickpeas‐compressed.pdf http://www.pulseaus.com.au/growing‐pulses/bmp/chickpea/ascochyta‐blight http://agriculture.vic.gov.au/agriculture/pests‐diseases‐and‐weeds/plant‐diseases/grains‐pulses‐and‐cereals/ascochyta‐blight‐of‐chickpea http://www.croppro.com.au/crop_disease_manual/ch05s02.php https://www.northernpulse.com/uploads/resources/722/handout‐chickpeaascochyta‐nov13‐2011.pdf http://oar.icrisat.org/184/1/24_2010_IB_no_82_Host_Plant https://www.crop.bayer.com.au/find‐crop‐solutions/by‐pest/diseases/ascochyta‐blight