MotivationAccurate detection, genotyping and downstream analysis of genomic variants from high-throughput sequencing data are fundamental features in modern production pipelines for genetic-based diagnosis in medicine or genomic selection in plant and animal breeding. Our research group maintains the Next-Generation Sequencing Experience Platform (NGSEP) as a precise, efficient and easy-to-use software solution for these features.ResultsUnderstanding that incorrect alignments around short tandem repeats are an important source of genotyping errors, we implemented in NGSEP new algorithms for realignment and haplotype clustering of reads spanning indels and short tandem repeats. We performed extensive benchmark experiments comparing NGSEP to state-of-the-art software using real data from three sequencing protocols and four species with different distributions of repetitive elements. NGSEP consistently shows comparative accuracy and better efficiency compared to the existing solutions. We expect that this work will contribute to the continuous improvement of quality in variant calling needed for modern applications in medicine and agriculture.Availability and implementationNGSEP is available as open source software at http://ngsep.sf.net.Supplementary information Supplementary data are available at Bioinformatics online.

The subcellular localization of proteins is very important for characterizing its function in a cell. Accurate prediction of the subcellular locations in computational paradigm has been an active area of interest. Most of the work has been focused on single localization prediction. Only few studies have discussed the multi-target localization, but have not achieved good accuracy so far; in plant sciences, very limited work has been done. Here we report the development of a novel tool Plant-mSubP, which is based on integrated machine learning approaches to efficiently predict the subcellular localizations in plant proteomes. The proposed approach predicts with high accuracy 11 single localizations and three dual locations of plant cell. Several hybrid features based on composition and physicochemical properties of a protein such as amino acid composition, pseudo amino acid composition, auto-correlation descriptors, quasi-sequence-order descriptors and hybrid features are used to represent the protein. The performance of the proposed method has been assessed through a training set as well as an independent test set. Using the hybrid feature of the pseudo amino acid composition, N-Center-C terminal amino acid composition and the dipeptide composition (PseAAC-NCC-DIPEP), an overall accuracy of 81.97 %, 84.75 % and 87.88 % is achieved on the training data set of proteins containing the single-label, single- and dual-label combined, and dual-label proteins, respectively. When tested on the independent data, an accuracy of 64.36 %, 64.84 % and 81.08 % is achieved on the single-label, single- and dual-label, and dual-label proteins, respectively. The prediction models have been implemented on a web server available at http://bioinfo.usu.edu/Plant-mSubP/. The results indicate that the proposed approach is comparable to the existing methods in single localization prediction and outperforms all other existing tools when compared for dual-label proteins. The prediction tool will be a useful resource for better annotation of various plant proteomes.

The aerobic, Gram-negative motile bacillus, Burkholderia pseudomallei is a facultative intracellular bacterium causing melioidosis, a critical disease of public health importance, which is widely endemic in the tropics and subtropical regions of the world. Melioidosis is associated with high case fatality rates in animals and humans; even with treatment, its mortality is 20–50%. It also infects plants and is designated as a biothreat agent. B. pseudomallei is pathogenic due to its ability to invade, resist factors in serum and survive intracellularly. Despite its importance, to date only a few effector proteins have been functionally characterized, and there is not much information regarding the host–pathogen protein–protein interactions (PPI) of this system, which are important to studying infection mechanisms and thereby develop prevention measures. We explored two computational approaches, the homology-based interolog and the domain-based method, to predict genome-scale host–pathogen interactions (HPIs) between two different strains of B. pseudomallei (prototypical, and highly virulent) and human. In total, 76 335 common HPIs (between the two strains) were predicted involving 8264 human and 1753 B. pseudomallei proteins. Among the unique PPIs, 14 131 non-redundant HPIs were found to be unique between the prototypical strain and human, compared to 3043 non-redundant HPIs between the highly virulent strain and human. The protein hubs analysis showed that most B. pseudomallei proteins formed a hub with human dnaK complex proteins associated with tuberculosis, a disease similar in symptoms to melioidosis. In addition, drug-binding and carbohydrate-binding mechanisms were found overrepresented within the host–pathogen network, and metabolic pathways were frequently activated according to the pathway enrichment. Subcellular localization analysis showed that most of the pathogen proteins are targeting human proteins inside cytoplasm and nucleus. We also discovered the host targets of the drug-related pathogen proteins and proteins that form T3SS and T6SS in B. pseudomallei. Additionally, a comparison between the unique PPI patterns present in the prototypical and highly virulent strains was performed. The current study is the first report on developing a genome-scale host–pathogen protein interaction networks between the human and B. pseudomallei, a critical biothreat agent. We have identified novel virulence factors and their interacting partners in the human proteome. These PPIs can be further validated by high-throughput experiments and may give new insights on how B. pseudomallei interacts with its host, which will help medical researchers in developing better prevention measures.

Motivation Understanding the mechanisms underlying infectious diseases is fundamental to develop prevention strategies. Host-Pathogen Interactions (HPI) are actively studied worldwide to find potential genomic targets for the development of novel drugs, vaccines, and other therapeutics. Determining which proteins are involved in the interaction system behind an infectious process is the first step to develop an efficient disease control strategy. Very few computational methods have been implemented as web services to infer novel HPIs, and there is not a single framework which combines several of those approaches to produce and visualize a comprehensive analysis of host-pathogen interactions. Results Here, we introduce PredHPI, a powerful framework that integrates both the detection and visualization of interaction networks in a single web service, facilitating the apprehension of model and non-model host-pathogen systems to aid the biologists in building hypotheses and designing appropriate experiments. PredHPI is built on high-performance computing resources on the backend capable of handling proteome-scale sequence data from both the host as well as pathogen. Data are displayed in an information-rich and interactive visualization, which can be further customized with user-defined layouts. We believe PredHPI will serve as an invaluable resource to diverse experimental biologists and will help advance the research in the understanding of complex infectious diseases. Availability and implementation PredHPI tool is freely available at http://bioinfo.usu.edu/PredHPI/. Supplementary information All the supplementary data, figures S1–S3, and excel files S1–S3 are available at Bioinformatics online.

Recent developments in High Throughput Sequencing (HTS) technologies and bioinformatics, including improved read lengths and genome assemblers allow the reconstruction of complex genomes with unprecedented quality and contiguity. Sugarcane has one of the most complicated genomes among grassess with a haploid length of 1Gbp and a ploidies between 8 and 12. In this work, we present a genome assembly of the Colombian sugarcane hybrid CC 01-1940. Three types of sequencing technologies were combined for this assembly: PacBio long reads, Illumina paired short reads, and Hi-C reads. We achieved a median contig length of 34.94 Mbp and a total genome assembly of 903.2 Mbp. We annotated a total of 63,724 protein coding genes and performed a reconstruction and comparative analysis of the sucrose metabolism pathway. Nucleotide evolution measurements between orthologs with close species suggest that divergence between Saccharum officinarum and Saccharum spontaneum occurred <2 million years ago. Synteny analysis between CC 01-1940 and the S. spontaneum genome confirms the presence of translocation events between the species and a random contribution throughout the entire genome in current sugarcane hybrids. Analysis of RNA-Seq data from leaf and root tissue of contrasting sugarcane genotypes subjected to water stress treatments revealed 17,490 differentially expressed genes, from which 3,633 correspond to genes expressed exclusively in tolerant genotypes. We expect the resources presented here to serve as a source of information to improve the selection processes of new varieties of the breeding programs of sugarcane.