A mapping of an ensemble of mitochondrial sequences for various organisms into 3D space based on the word composition

To visualize a bird's-eye view of an ensemble of mitochondrial genome sequences for various species, we recently developed a novel method of mapping a biological sequence ensemble into Three-Dimensional (3D) vector space. First, we represented a biological sequence of a species s by a word-composition vector x(s), where its length ∣x(s)∣represents the sequence length, and its unit vector x(s)/∣x(s)∣represents the relative composition of the K-tuple words through the sequence and the size of the dimension, N=4(K), is the number of all possible words with the length of K. Second, we mapped the vector x(s) to the 3D position vector y(s), based on the two following simple principles: (1) ∣y(s)∣=∣x(s)∣and (2) the angle between y(s) and y(t) maximally correlates with the angle between x(s) and x(t). The mitochondrial genome sequences for 311 species, including 177 Animalia, 85 Fungi and 49 Green plants, were mapped into 3D space by using K=7. The mapping was successful because the angles between vectors before and after the mapping highly correlated with each other (correlation coefficients were 0.92-0.97). Interestingly, the Animalia kingdom is distributed along a single arc belt (just like the Milky Way on a Celestial Globe), and the Fungi and Green plant kingdoms are distributed in a similar arc belt. These two arc belts intersect at their respective middle regions and form a cross structure just like a jet aircraft fuselage and its wings. This new mapping method will allow researchers to intuitively interpret the visual information presented in the maps in a highly effective manner.     

Comparative protein analysis of the chitin metabolic pathway in extant organisms: A complex network approach

Chitin is a structural endogenous carbohydrate, which is a major component of fungal cell walls and arthropod exoskeletons. A renewable resource and the second most abundant polysaccharide in nature after cellulose, chitin is currently used for waste water clearing, cosmetics, medical, and veterinary applications. This work comprises data mining of protein sequences related to the chitin metabolic pathway of completely sequenced genomes of extant organisms pertaining to the three life domains, followed by meta-analysis using traditional sequence similarity comparison and complex network approaches. Complex networks involving proteins of the chitin metabolic pathway in extant organisms were constructed based on protein sequence similarity. Several usual network indices were estimated in order to obtain information on the topology of these networks, including those related to higher order neighborhood properties. Due to the assumed evolutionary character of the system, we also discuss issues related to modularity properties, with the concept of edge betweenness playing a particularly important role in our analysis. Complex network approach correctly identifies clusters of organisms that belong to phylogenetic groups without any a priori knowledge about the biological features of the investigated protein sequences. We envisage the prospect of using such a complex network approach as a high-throughput phylogenetic method.     

Dissecting of the FHB resistance QTL on the short arm of wheat chromosome 2D using a comparative genomic approach: from QTL to candidate gene

Colinearity in gene content and order between rice and closely related cereal crops has been a powerful tool for gene identification. Using a comparative genomic approach, we have identified the rice genomic region syntenous to the region of the short arm of wheat chromosome 2D, on which quantitative trait loci (QTLs) for Fusarium head blight (FHB) resistance and for controlling accumulation of the mycotoxin deoxynivalenol (DON) are closely located. Utilizing markers known to reside near the FHB resistance QTL and data from several wheat genetic maps, we have limited the syntenous region to 6.8 Mb of the short arm of rice chromosome 4. From the 6.8-Mb sequence of rice chromosome 4, we found three putative rice genes that could have a role in detoxification of mycotoxins. DNA sequences of these putative rice genes were used in BLAST searches to identify wheat expressed sequence tags (ESTs) exhibiting significant similarity. Combined data from expression analysis and gene mapping of wheat homologues and results of analysis of DON accumulation using doubled haploid populations revealed that a putative gene for multidrug resistance-associated protein (MRP) is a possible candidate for the FHB resistance and/or DON accumulation controlling QTLs on wheat chromosome 2DS and can be used as a molecular marker to eliminate the susceptible allele when the Chinese wheat variety Sumai 3 is used as a resistance source. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A novel giant gene CSMD3 encoding a protein with CUB and sushi multiple domains: a candidate gene for benign adult familial myoclonic epilepsy on human chromosome 8q23.3-q24.1

We identified a novel giant gene encoding a transmembrane protein with CUB and sushi multiple domains on the human chromosome 8q23.3-q24.1 in which benign adult familial myoclonic epilepsy type 1 (BAFME1/FAME, OMIM:601068) has been mapped. This giant gene consists of 73 exons and spans over 1.2Mb on the genomic DNA region. It showed significant homology to two genes, CSMD1 gene on 8p23 and CSMD2 gene on 1p34, at reduced amino acid sequence level and hence we designated as CSMD3. The CSMD3 gene was expressed mainly in adult and fetal brains. We performed mutation analysis on the CSMD3 gene for seven patients with BAFME1/FAME, but no mutation was found in the coding sequence of the CSMD3 gene. Comparative genomic analysis revealed a conserved family of CSMD genes in the mouse and fugu genomes. Possible functions of the CSMD gene family are discussed.     

The complete genome of Burkholderia phenoliruptrix strain BR3459a,a symbiont of Mimosa flocculosa: highlighting the coexistence of symbiotic and pathogenic genes

Background

Burkholderia species play an important ecological role related to xenobiosis, the promotion of plant growth, the biocontrol of agricultural diseases, and symbiotic and non-symbiotic biological nitrogen fixation. Here, we highlight our study as providing the first complete genome of a symbiotic strain of B. phenoliruptrix, BR3459a (=CLA1), which was originally isolated in Brazil from nodules of Mimosa flocculosa and is effective in fixing nitrogen in association with this leguminous species.

Results

Genomic comparisons with other pathogenic and non-pathogenic Burkholderia strains grouped B. phenoliruptrix BR3459a with plant-associated beneficial and environmental species, although it shares a high percentage of its gene repertoire with species of the B. cepacia complex (Bcc) and "pseudomallei" group. The genomic analyses showed that the bce genes involved in exopolysaccharide production are clustered together in the same genomic region, constituting part of the Group III cluster of non-pathogenic bacteria. Regarding environmental stresses, we highlight genes that might be relevant in responses to osmotic, heat, cold and general stresses. Furthermore, a number of particularly interesting genes involved in the machinery of the T1SS, T2SS, T3SS, T4ASS and T6SS secretion systems were identified. The xenobiotic properties of strain BR3459a were also investigated, and some enzymes involved in the degradation of styrene, nitrotoluene, dioxin, chlorocyclohexane, chlorobenzene and caprolactam were identified. The genomic analyses also revealed a large number of antibiotic-related genes, the most important of which were correlated with streptomycin and novobiocin. The symbiotic plasmid showed high sequence identity with the symbiotic plasmid of B. phymatum. Additionally, comparative analysis of 545 housekeeping genes among pathogenic and non-pathogenic Burkholderia species strongly supports the definition of a new genus for the second branch, which would include BR3459a.

Conclusions

The analyses of B. phenoliruptrix BR3459a showed key property of fixing nitrogen that together with genes for high tolerance to environmental stresses might explain a successful strategy of symbiosis in the tropics. The strain also harbours interesting sets of genes with biotechnological potential. However, the resemblance of certain genes to those of pathogenic Burkholderia raise concerns about large-scale applications in agriculture or for bioremediation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-535) contains supplementary material, which is available to authorized users.     

Selectivity of bacterial proteome fractionation based on differential solubility: a mass spectrometry evaluation

We investigate the selectivity achieved after differential solubilization of bacterial proteomes following two procedures, both based on successive extraction of proteins in solutions of increasing solubilizing power. Recently, these procedures have gained notable popularity and several commercial kits are now available. A total of 225 proteins in one case and 227 proteins in the other were identified by LC MSMS analysis; 146 of them were identified in both procedures. The proportions of proteins identified as present in only one fraction were 64 and 57%, respectively. The distribution of cytosolic, membrane, and ribosomal proteins among the successive extracts was analyzed in detail. The effect of (1) replacement of low-speed with high-speed centrifugation, (2) omission of detergents in urea solutions, (3) successive washes of pellets, and (4) reproducibility was evaluated. Proteins with positive grand averages of hydropathicity values and membrane proteins were found in all fractions. This study highlights the benefits and limitations of differential solubilization methods, focusing on practical aspects that may strongly influence their selectivity.     

LB3D: a protein three-dimensional substructure search program based on the lower bound of a root mean square deviation value

Searching for protein structure-function relationships using three-dimensional (3D) structural coordinates represents a fundamental approach for determining the function of proteins with unknown functions. Since protein structure databases are rapidly growing in size, the development of a fast search method to find similar protein substructures by comparison of protein 3D structures is essential. In this article, we present a novel protein 3D structure search method to find all substructures with root mean square deviations (RMSDs) to the query structure that are lower than a given threshold value. Our new algorithm runs in O(m + N/m(0.5)) time, after O(N log N) preprocessing, where N is the database size and m is the query length. The new method is 1.8-41.6 times faster than the practically best known O(N) algorithm, according to computational experiments using a huge database (i.e., >20,000,000 C-alpha coordinates).     

Theoretical 3D model of histamine N-methyltransferase: insights into the effects of a genetic polymorphism on enzymatic activity and thermal stability

Histamine N-methyltransferase (HNMT) catalyzes the N-methylation of histamine in mammals. The experimentally determined HNMT three-dimensional (3D) structure is not available. However, there is a common genetic polymorphism for human HNMT (Thr105Ile) that reduces enzymatic activity and is a risk factor for asthma. To obtain insights into mechanisms responsible for the effects of that polymorphism on enzymatic activity and thermal stability, we predicted the 3D structure of HNMT using the threading method and molecular dynamics simulations in water. Herein, we report a theoretical 3D model of human HNMT which reveals that polymorphic residue Thr105Ile is located in the turn between a beta strand and an alpha helix on the protein surface away from the active site of HNMT. Ile105 energetically destabilizes folded HNMT because of its low Chou-Fasman score for forming a turn conformation and the exposure of its hydrophobic side chain to aqueous solution. It thus promotes the formation of misfolded proteins that are prone to the clearance by proteasomes. This information explains, for the first time, how genetic polymorphisms can cause enhanced protein degradation and why the thermal stability of allozyme Ile105 is lower than that of Thr105. It also supports the hypothesis that the experimental observation of a significantly lower level of HNMT enzymatic activity for allozyme Ile105 than that with Thr105 is due to a decreased concentration of allozyme Ile105, but not an alternation of the active-site topology of HNMT caused by the difference at residue 105.     

Allergome: the characterization of allergens based on a 2D gel electrophoresis approach

Type I hypersensitivity reactions are in constant progression in industrialized countries. The physiopathologic mechanism of these diseases implicates the production of specific immunoglobulin (Ig)E to allergenic molecules, their binding to the Fcε receptor on the surface of mast cells and basophils, and the release of inflammatory mediators when allergens are introduced into the body and crosslink with the IgE bound to the cell surface. An allergen is defined as a molecule that induces the production of, and binds to, IgE. The identification of the allergenic molecules is an important goal to improve diagnosis and treatment of allergy. This characterization aims to extract proteins from the allergenic source, to analyze IgE specificity by immunoblotting and to identify the proteins that bind IgE.     

A workflow for peptide-based proteomics in a poorly sequenced plant: a case study on the plasma membrane proteome of banana

Membrane proteins are an interesting class of proteins because of their functional importance. Unfortunately their analysis is hampered by low abundance and poor solubility in aqueous media. Since shotgun methods are high-throughput and partly overcome these problems, they are preferred for membrane proteomics. However, their application in non-model plants demands special precautions to prevent false positive identification of proteins. In the current paper, a workflow for membrane proteomics in banana, a poorly sequenced plant, is proposed. The main steps of this workflow are (i) optimization of the peptide separation, (ii) performing de novo sequencing to allow a sequence homology search and (iii) visualization of identified peptide-protein associations using Cytoscape to remove redundancy and wrongly assigned peptides, based on species-specific information. By applying this workflow, integral plasma membrane proteins from banana leaves were successfully identified.     

The carbon and hydrogen stable isotope composition of bacteriogenic methane: A laboratory study using a landfill inoculum

Anaerobic bacterial degradation of landfill waste produces a globally significant source of the greenhouse gas methane. Stable isotopic measurements of methane [δ^I3C(CH₄) and δD(CH₄)] can often fingerprint different sources of methane (natural vs. anthro‐pogenic) and help identify the bacterial processes involved in methane production. Landfill microbial communities are complex and diverse, and hence so too is the biogeochem‐istry of methane formation. To investigate the influence of (l) the methane formation pathway (acetoclastic methanogenesis and CO₂ reduction), and (2) SD of water on the stable isotopic composition of landfill methane, two model butyrate‐degrading landfill systems were established. The systems were inoculated with domestic refuse from a landfill and incubated in the laboratory for 92 days. Both systems were identical except δD of water initially added to system 2 was 118% heavier than system 1. Between days 39 and 72 the systems were resupplemented with butyrate. Production of CH₄ and CO₂ and changes in volatile fatty acid concentration confirmed that active methanogenic populations had been established. CH₄ became ¹³C enriched in both incubations with time. Interpreting changes in acetate, butyrate, and propionate concentration during incubation is complicated, but these observations and other information suggest that the dominant methanogenic substrate changed front CO₂/H₂ to acetate as the experiment progressed. This is also consistent with the observed ¹³C enrichment of CH₄, as ¹³C discrimination during methane production from acetate is less than from CO₂. In contrast, δD(CH₄) remained relatively constant, suggesting that this measurement may not provide a reliable basis for distinguishing between CH₄ from CO₂ reduction and acetoclastic methanogenesis, as has previously been suggested.     

A requiem to the nuclear matrix: from a controversial concept to 3D organization of the nucleus

The first papers coining the term “nuclear matrix” were published 40 years ago. Here, we review the data obtained during the nuclear matrix studies and discuss the contribution of this controversial concept to our current understanding of nuclear architecture and three-dimensional organization of genome.     

Open questions on the 3D structures of collagen containing vertebrate mineralized tissues: A perspective

Our current understanding of the structures of vertebrate mineralized tissues is largely based on light microscopy/histology and projections of 3D structures onto 2D planes using electron microscopy. We know little about the fine details of these structures in 3D at the length scales of their basic building blocks, the inherent variations of structure within a tissue and the cell-extracellular tissue interfaces. This limits progress in understanding tissue formation, relating structure to mechanical and metabolic functions, and obtaining deeper insights into pathologies and the evolution of these tissues. In this perspective we identify and discuss a series of open questions pertaining to collagen containing vertebrate mineralized tissues that can be addressed using appropriate 3D structural determination methods. By so doing we hope to encourage more research into the 3D structures of mineralized vertebrate tissues.     

A framework for querying a database for structural information on 3D images of macromolecules: A web-based query-by-content prototype on the BioImage macromolecular server

Nowadays we are experiencing a remarkable growth in the number of databases that have become accessible over the Web. However, in a certain number of cases, for example, in the case of BioImage, this information is not of a textual nature, thus posing new challenges in the design of tools to handle these data. In this work, we concentrate on the development of new mechanisms aimed at "querying" these databases of complex data sets by their intrinsic content, rather than by their textual annotations only. We concentrate our efforts on a subset of BioImage containing 3D images (volumes) of biological macromolecules, implementing a first prototype of a "query-by-content" system. In the context of databases of complex data types the term query-by-content makes reference to those data modeling techniques in which user-defined functions aim at "understanding" (to some extent) the informational content of the data sets. In these systems the matching criteria introduced by the user are related to intrinsic features concerning the 3D images themselves, hence, complementing traditional queries by textual key words only. Efficient computational algorithms are required in order to "extract" structural information of the 3D images prior to storing them in the database. Also, easy-to-use interfaces should be implemented in order to obtain feedback from the expert. Our query-by-content prototype is used to construct a concrete query, making use of basic structural features, which are then evaluated over a set of three-dimensional images of biological macromolecules. This experimental implementation can be accessed via the Web at the BioImage server in Madrid, at http://www.bioimage.org/qbc/index.html.     

A method of focused classification, based on the bootstrap 3D variance analysis, and its application to EF-G-dependent translocation

The bootstrap-based method for calculation of the 3D variance in cryo-EM maps reconstructed from sets of their projections was applied to a dataset of functional ribosomal complexes containing the Escherichia coli 70S ribosome, tRNAs, and elongation factor G (EF-G). The variance map revealed regions of high variability in the intersubunit space of the ribosome: in the locations of tRNAs, in the putative location of EF-G, and in the vicinity of the L1 protein. This result indicated heterogeneity of the dataset. A method of focused classification was put forward in order to sort out the projection data into approximately homogenous subsets. The method is based on the identification and localization of a region of high variance that a subsequent classification step can be focused on by the use of a 3D spherical mask. After initial classification, template volumes are created and are subsequently refined using a multireference 3D projection alignment procedure. In the application to the ribosome dataset, the two resulting structures were interpreted as resulting from ribosomal complexes with bound EF-G and an empty A site, or, alternatively, from complexes that had no EF-G bound but had both A and P sites occupied by tRNA. The proposed method of focused classification proved to be a successful tool in the analysis of the heterogeneous cryo-EM dataset. The associated calculation of the correlations within the density map confirmed the conformational variability of the complex, which could be interpreted in terms of the ribosomal elongation cycle.     

Determination of n-3 HUFA content in Atlantic salmon flesh based on the lipid content, morphometric measurements and blood fatty acid composition: A modeling approach

Fish are the most important source of n-3 highly unsaturated fatty acids (HUFA) in the human diet and, with declining wild stocks, an increasing proportion is being provided by aquaculture. Paradoxically, aquaculture fish diets have traditionally incorporated fish oil and meal derived from wild fisheries. Continued aquaculture development requires fish oil to be replaced with vegetable oils, the only sustainable alternative. However, vegetable oils lack n-3 HUFA and so flesh from fish reared on these diets can also have reduced n-3 HUFA and thus reduced nutritional quality. This accepted, the flesh n-3 HUFA content should be an economically important trait, however to be included in the breeding goal the trait must be measurable. In the present study, we investigated whether the flesh n-3 HUFA content of salmon can be estimated in a non-fatal way. We showed that a general regression model based on flesh lipid content, morphometric and blood fatty acid measurements could estimate and predict flesh n-3 HUFA content. This would allow a choice from a range of selection methods, including mass selection or within-family selection, if this important flesh quality trait is to be included in future salmon breeding programmes.     

A molecular modeling study on the enantioselectivity of aryl alkyl ketone reductions by a NADPH-dependent carbonyl reductase

Automated structural analysis of Sporobolomyces salmonicolor carbonyl reductase (SSCR) indicates that the two largest potential receptor sites are in the vicinity of the nicotinamide reductant. The largest receptor site is a scalene triangle with sides of approximately 8 A by 9 A by 13 A, which is narrow in width; one corner is surrounded by hydrophilic residues that can favorably bond with the ketone oxygen. Docking aryl alkyl ketones shows a distinct preference for binding to the largest receptor site, and for conformations that place the carbonyl oxygen of the substrate in the hydrophilic corner of the largest receptor site. Favorable docking conformations for aryl alkyl ketones fall into two low-energy ensembles. These conformational ensembles are distinguished by the positions of the substituents, presenting either the Si- or Re-face of the ketone to the nicotinamide reductant. For the ketones investigated here, there is a correspondence between the major enantiomer of the alcohol obtained from the reduction of the ketone and the conformer found to have the most stable interaction energy with the receptor site in all cases. The receptor site modeling, docking simulations, molecular dynamics, and enzyme-substrate geometry optimizations lead to a model for understanding the enantioselectivity of this NADPH-dependent carbonyl reductase.