dna: An R package for differential network analysis

Differential network analysis provides a framework for examining if there is sufficient statistical evidence to conclude that the structure of a network differs under two experimental conditions or if the structures of two networks are different. The R package dna provides tools and procedures for differential network analysis of genomic data. The focus of this package is on gene-gene networks, but the methods are easily adaptable for more general biological processes. This package includes preprocessing tools for simultaneously preparing a pair of networks for analysis, procedures for computing connectivity scores between pairs of genes based on many available statistical techniques, and tools for handling modules of genes based on these scores. Also, procedures are provided for performing permutation tests based on these scores to determine if the connectivity of a gene differs between the two networks, to determine if the connectivity of a particular set of important genes differs between the two networks, and to determine if the overall module structure differs between the two networks. Several built-in options are available for the types of scores and distances used in the testing procedures, and additionally, the procedures provide flexible methods that allow the user to define custom scores and distances.

Availability

dna is freely available at The Comprehensive R Archive Network, http://CRAN.R-project.org/package=dna     

Profiling small RNA reveals multimodal substructural signals in a Boltzmann ensemble

As the biomedical impact of small RNAs grows, so does the need to understand competing structural alternatives for regions of functional interest. Suboptimal structure analysis provides significantly more RNA base pairing information than a single minimum free energy prediction. Yet computational enhancements like Boltzmann sampling have not been fully adopted by experimentalists since identifying meaningful patterns in this data can be challenging. Profiling is a novel approach to mining RNA suboptimal structure data which makes the power of ensemble-based analysis accessible in a stable and reliable way. Balancing abstraction and specificity, profiling identifies significant combinations of base pairs which dominate low-energy RNA secondary structures. By design, critical similarities and differences are highlighted, yielding crucial information for molecular biologists. The code is freely available via http://gtfold.sourceforge.net/profiling.html.     

ncRDeathDB: A comprehensive bioinformatics resource for deciphering network organization of the ncRNA-mediated cell death system

Programmed cell death (PCD) is a critical biological process involved in many important processes, and defects in PCD have been linked with numerous human diseases. In recent years, the protein architecture in different PCD subroutines has been explored, but our understanding of the global network organization of the noncoding RNA (ncRNA)-mediated cell death system is limited and ambiguous. Hence, we developed the comprehensive bioinformatics resource (ncRDeathDB, www.rna-society.org/ncrdeathdb) to archive ncRNA-associated cell death interactions. The current version of ncRDeathDB documents a total of more than 4600 ncRNA-mediated PCD entries in 12 species. ncRDeathDB provides a user-friendly interface to query, browse and manipulate these ncRNA-associated cell death interactions. Furthermore, this resource will help to visualize and navigate current knowledge of the noncoding RNA component of cell death and autophagy, to uncover the generic organizing principles of ncRNA-associated cell death systems, and to generate valuable biological hypotheses.     

Complete genome sequence of Granulicella mallensis type strain MP5ACTX8T,an acidobacterium from tundra soil

Granulicella mallensis MP5ACTX8^T is a novel species of the genus Granulicella in subdivision 1of Acidobacteria. G. mallensis is of ecological interest being a member of the dominant soil bacterial community active at low temperatures and nutrient limiting conditions in Arctic alpine tundra. G. mallensis is a cold-adapted acidophile and a versatile heterotroph that hydrolyzes a suite of sugars and complex polysaccharides. Genome analysis revealed metabolic versatility with genes involved in metabolism and transport of carbohydrates. These include gene modules encoding the carbohydrate-active enzyme (CAZyme) family involved in breakdown, utilization and biosynthesis of diverse structural and storage polysaccharides including plant based carbon polymers. The genome of Granulicella mallensis MP5ACTX8^T consists of a single replicon of 6,237,577 base pairs (bp) with 4,907 protein-coding genes and 53 RNA genes.     

RNAMotifScan: automatic identification of RNA structural motifs using secondary structural alignment

Recent studies have shown that RNA structural motifs play essential roles in RNA folding and interaction with other molecules. Computational identification and analysis of RNA structural motifs remains a challenging task. Existing motif identification methods based on 3D structure may not properly compare motifs with high structural variations. Other structural motif identification methods consider only nested canonical base-pairing structures and cannot be used to identify complex RNA structural motifs that often consist of various non-canonical base pairs due to uncommon hydrogen bond interactions. In this article, we present a novel RNA structural alignment method for RNA structural motif identification, RNAMotifScan, which takes into consideration the isosteric (both canonical and non-canonical) base pairs and multi-pairings in RNA structural motifs. The utility and accuracy of RNAMotifScan is demonstrated by searching for kink-turn, C-loop, sarcin-ricin, reverse kink-turn and E-loop motifs against a 23S rRNA (PDBid: 1S72), which is well characterized for the occurrences of these motifs. Finally, we search these motifs against the RNA structures in the entire Protein Data Bank and the abundances of them are estimated. RNAMotifScan is freely available at our supplementary website (http://genome.ucf.edu/RNAMotifScan).     

Evidence for a Nonendosomal Function of the Saccharomyces cerevisiae ESCRT-III-Like Protein Chm7

Endosomal sorting complex required for transport (ESCRT) proteins are involved in a number of cellular processes, such as endosomal protein sorting, HIV budding, cytokinesis, plasma membrane repair, and resealing of the nuclear envelope during mitosis. Here we explored the function of a noncanonical member of the ESCRT-III protein family, the Saccharomyces cerevisiae ortholog of human CHMP7. Very little is known about this protein. In silico analysis predicted that Chm7 (yeast ORF YJL049w) is a fusion of an ESCRT-II and ESCRT-III-like domain, which would suggest a role in endosomal protein sorting. However, our data argue against a role of Chm7 in endosomal protein sorting. The turnover of the endocytic cargo protein Ste6 and the vacuolar protein sorting of carboxypeptidase S (CPS) were not affected by CHM7 deletion, and Chm7 also responded very differently to a loss in Vps4 function compared to a canonical ESCRT-III protein. Our data indicate that the Chm7 function could be connected to the endoplasmic reticulum (ER). In line with a function at the ER, we observed a strong negative genetic interaction between the deletion of a gene function (APQ12) implicated in nuclear pore complex assembly and messenger RNA (mRNA) export and the CHM7 deletion. The patterns of genetic interactions between the APQ12 deletion and deletions of ESCRT-III genes, two-hybrid interactions, and the specific localization of mCherry fusion proteins are consistent with the notion that Chm7 performs a novel function at the ER as part of an alternative ESCRT-III complex.     

Using the Fast Fourier Transform to Accelerate the Computational Search for RNA Conformational Switches

Using complex roots of unity and the Fast Fourier Transform, we design a new thermodynamics-based algorithm, FFTbor, that computes the Boltzmann probability that secondary structures differ by base pairs from an arbitrary initial structure of a given RNA sequence. The algorithm, which runs in quartic time and quadratic space , is used to determine the correlation between kinetic folding speed and the ruggedness of the energy landscape, and to predict the location of riboswitch expression platform candidates. A web server is available at http://bioinformatics.bc.edu/clotelab/FFTbor/.     

Analysis of PHA-1 Reveals a Limited Role in Pharyngeal Development and Novel Functions in Other Tissues

PHA-1 encodes a cytoplasmic protein that is required for embryonic morphogenesis and attachment of the foregut (pharynx) to the mouth (buccal capsule). Previous reports have in some cases suggested that PHA-1 is essential for the differentiation of most or all pharyngeal cell types. By performing mosaic analysis with a recently acquired pha-1 null mutation (tm3671), we found that PHA-1 is not required within most or all pharyngeal cells for their proper specification, differentiation, or function. Rather, our evidence suggests that PHA-1 acts in the arcade or anterior epithelial cells of the pharynx to promote attachment of the pharynx to the future buccal capsule. In addition, PHA-1 appears to be required in the epidermis for embryonic morphogenesis, in the excretory system for osmoregulation, and in the somatic gonad for normal ovulation and fertility. PHA-1 activity is also required within at least a subset of intestinal cells for viability. To better understand the role of PHA-1 in the epidermis, we analyzed several apical junction markers in pha-1(tm3671) homozygous embryos. PHA-1 regulates the expression of several components of two apical junction complexes including AJM-1–DLG-1/discs large complex and the classical cadherin–catenin complex, which may account for the role of PHA-1 in embryonic morphogenesis.     

Complete genome sequence of Terriglobus saanensis type strain SP1PR4T,an Acidobacteria from tundra soil

Terriglobus saanensis SP1PR4^T is a novel species of the genus Terriglobus. T. saanensis is of ecological interest because it is a representative of the phylum Acidobacteria, which are dominant members of bacterial soil microbiota in Arctic ecosystems. T. saanensis is a cold-adapted acidophile and a versatile heterotroph utilizing a suite of simple sugars and complex polysaccharides. The genome contained an abundance of genes assigned to metabolism and transport of carbohydrates including gene modules encoding for carbohydrate-active enzyme (CAZyme) family involved in breakdown, utilization and biosynthesis of diverse structural and storage polysaccharides. T. saanensis SP1PR4^T represents the first member of genus Terriglobus with a completed genome sequence, consisting of a single replicon of 5,095,226 base pairs (bp), 54 RNA genes and 4,279 protein-coding genes. We infer that the physiology and metabolic potential of T. saanensis is adapted to allow for resilience to the nutrient-deficient conditions and fluctuating temperatures of Arctic tundra soils.     

Nucleoporin FG Domains Facilitate mRNP Remodeling at the Cytoplasmic Face of the Nuclear Pore Complex

Directional export of messenger RNA (mRNA) protein particles (mRNPs) through nuclear pore complexes (NPCs) requires multiple factors. In Saccharomyces cerevisiae, the NPC proteins Nup159 and Nup42 are asymmetrically localized to the cytoplasmic face and have distinct functional domains: a phenylalanine-glycine (FG) repeat domain that docks mRNP transport receptors and domains that bind the DEAD-box ATPase Dbp5 and its activating cofactor Gle1, respectively. We speculated that the Nup42 and Nup159 FG domains play a role in positioning mRNPs for the terminal mRNP-remodeling steps carried out by Dbp5. Here we find that deletion (Δ) of both the Nup42 and Nup159 FG domains results in a cold-sensitive poly(A)+ mRNA export defect. The nup42ΔFG nup159ΔFG mutant also has synthetic lethal genetic interactions with dbp5 and gle1 mutants. RNA cross-linking experiments further indicate that the nup42ΔFG nup159ΔFG mutant has a reduced capacity for mRNP remodeling during export. To further analyze the role of these FG domains, we replaced the Nup159 or Nup42 FG domains with FG domains from other Nups. These FG “swaps” demonstrate that only certain FG domains are functional at the NPC cytoplasmic face. Strikingly, fusing the Nup42 FG domain to the carboxy-terminus of Gle1 bypasses the need for the endogenous Nup42 FG domain, highlighting the importance of proximal positioning for these factors. We conclude that the Nup42 and Nup159 FG domains target the mRNP to Gle1 and Dbp5 for mRNP remodeling at the NPC. Moreover, these results provide key evidence that character and context play a direct role in FG domain function and mRNA export.     

Draft genome sequence of marine alphaproteobacterial strain HIMB11, the first cultivated representative of a unique lineage within the Roseobacter clade possessing an unusually small genome

Strain HIMB11 is a planktonic marine bacterium isolated from coastal seawater in Kaneohe Bay, Oahu, Hawaii belonging to the ubiquitous and versatile Roseobacter clade of the alphaproteobacterial family Rhodobacteraceae. Here we describe the preliminary characteristics of strain HIMB11, including annotation of the draft genome sequence and comparative genomic analysis with other members of the Roseobacter lineage. The 3,098,747 bp draft genome is arranged in 34 contigs and contains 3,183 protein-coding genes and 54 RNA genes. Phylogenomic and 16S rRNA gene analyses indicate that HIMB11 represents a unique sublineage within the Roseobacter clade. Comparison with other publicly available genome sequences from members of the Roseobacter lineage reveals that strain HIMB11 has the genomic potential to utilize a wide variety of energy sources (e.g. organic matter, reduced inorganic sulfur, light, carbon monoxide), while possessing a reduced number of substrate transporters.     

Non-contiguous finished genome sequence and description of Oceanobacillus massiliensis sp. nov.

Oceanobacillus massiliensis strain N’Diop^T sp. nov. is the type strain of O. massiliensis sp. nov., a new species within the genus Oceanobacillus. This strain, whose genome is described here, was isolated from the fecal flora of a healthy patient. O. massiliensis is an aerobic rod. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 3,532,675 bp long genome contains 3,519 protein-coding genes and 72 RNA genes, including between 6 and 8 rRNA operons.     

Non-contiguous finished genome sequence and description of Fenollaria massiliensis gen. nov., sp. nov., a new genus of anaerobic bacterium

Fenollaria massiliensis strain 9401234^T, is the type strain of Fenollaria massiliensis gen. nov., sp. nov., a new species within a new genus Fenollaria. This strain, whose genome is described here, was isolated from an osteoarticular sample. F. massiliensis strain 9401234^T is an obligate anaerobic Gram-negative bacillus. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 1.71 Mbp long genome exhibits a G+C content of 34.46% and contains 1,667 protein-coding and 30 RNA genes, including 3 rRNA genes.     

Draft genome sequence of Rhodococcus rhodochrous strain ATCC 17895

Rhodococcus rhodochrous ATCC 17895 possesses an array of mono- and dioxygenases, as well as hydratases, which makes it an interesting organism for biocatalysis. R. rhodochrous is a Gram-positive aerobic bacterium with a rod-like morphology. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 6,869,887 bp long genome contains 6,609 protein-coding genes and 53 RNA genes. Based on small subunit rRNA analysis, the strain is more likely to be a strain of Rhodococcus erythropolis rather than Rhodococcus rhodochrous.