smyRNA: A Novel Ab Initio ncRNA Gene Finder

Background

Non-coding RNAs (ncRNAs) have important functional roles in the cell: for example, they regulate gene expression by means of establishing stable joint structures with target mRNAs via complementary sequence motifs. Sequence motifs are also important determinants of the structure of ncRNAs. Although ncRNAs are abundant, discovering novel ncRNAs on genome sequences has proven to be a hard task; in particular past attempts for ab initio ncRNA search mostly failed with the exception of tools that can identify micro RNAs.

Methodology/Principal Findings

We present a very general ab initio ncRNA gene finder that exploits differential distributions of sequence motifs between ncRNAs and background genome sequences.

Conclusions/Significance

Our method, once trained on a set of ncRNAs from a given species, can be applied to a genome sequences of other organisms to find not only ncRNAs homologous to those in the training set but also others that potentially belong to novel (and perhaps unknown) ncRNA families. Availability: http://compbio.cs.sfu.ca/taverna/smyrna     

Homology modeling of major intrinsic proteins in rice,maize and <Emphasis Type="Italic">Arabidopsis</Emphasis>: comparative analysis of transmembrane helix association and aromatic/arginine selectivity filters

Background  

The major intrinsic proteins (MIPs) facilitate the transport of water and neutral solutes across the lipid bilayers. Plant MIPs are believed to be important in cell division and expansion and in water transport properties in response to environmental conditions. More than 30 MIP sequences have been identified in Arabidopsis thaliana, maize and rice. Plasma membrane intrinsic proteins (PIPs), tonoplast intrinsic proteins (TIPs), Nod26-like intrinsic protein (NIPs) and small and basic intrinsic proteins (SIPs) are subfamilies of plant MIPs. Despite sequence diversity, all the experimentally determined structures belonging to the MIP superfamily have the same "hour-glass" fold.     

MIPDB: a relational database dedicated to MIP family proteins

BACKGROUND INFORMATION: The MIPs (major intrinsic proteins) constitute a large family of membrane proteins that facilitate the passive transport of water and small neutral solutes across cell membranes. Since water is the most abundant molecule in all living organisms, the discovery of selective water-transporting channels called AQPs (aquaporins) has led to new knowledge on both the physiological and molecular mechanisms of membrane permeability. The MIPs are identified in Archaea, Bacteria and Eukaryota, and the rapid accumulation of new sequences in the database provides an opportunity for large-scale analysis, to identify functional and/or structural signatures or to infer evolutionary relationships. To help perform such an analysis, we have developed MIPDB (database for MIP proteins), a relational database dedicated to members of the MIP family. RESULTS: MIPDB is a motif-oriented database that integrates data on 785 MIP proteins from more than 200 organisms and contains 230 distinct sequence motifs. MIPDB proposes the classification of MIP proteins into three functional subgroups: AQPs, glycerol-uptake facilitators and aquaglyceroporins. Plant MIPs are classified into three specific subgroups according to their subcellular distribution in the plasma membrane, tonoplast or the symbiosome membrane. Some motifs of the database are highly selective and can be used to predict the transport function or subcellular localization of unknown MIP proteins. CONCLUSIONS: MIPDB offers a user-friendly and intuitive interface for a rapid and easy access to MIP resources and to sequence analysis tools. MIPDB is a web application, publicly accessible at http://idefix.univ-rennes1.fr:8080/Prot/index.html.     

Protein Aggregation Profile of the Bacterial Cytosol

Background

Protein misfolding is usually deleterious for the cell, either as a consequence of the loss of protein function or the buildup of insoluble and toxic aggregates. The aggregation behavior of a given polypeptide is strongly influenced by the intrinsic properties encoded in its sequence. This has allowed the development of effective computational methods to predict protein aggregation propensity.

Methodology/Principal Findings

Here, we use the AGGRESCAN algorithm to approximate the aggregation profile of an experimental cytosolic Escherichia coli proteome. The analysis indicates that the aggregation propensity of bacterial proteins is associated with their length, conformation, location, function, and abundance. The data are consistent with the predictions of other algorithms on different theoretical proteomes.

Conclusions/Significance

Overall, the study suggests that the avoidance of protein aggregation in functional environments acts as a strong evolutionary constraint on polypeptide sequences in both prokaryotic and eukaryotic organisms.     

Hidden localization motifs: naturally occurring peroxisomal targeting signals in non-peroxisomal proteins

Background

Can sequence segments coding for subcellular targeting or for posttranslational modifications occur in proteins that are not substrates in either of these processes? Although considerable effort has been invested in achieving low false-positive prediction rates, even accurate sequence-analysis tools for the recognition of these motifs generate a small but noticeable number of protein hits that lack the appropriate biological context but cannot be rationalized as false positives.

Results

We show that the carboxyl termini of a set of definitely non-peroxisomal proteins with predicted peroxisomal targeting signals interact with the peroxisomal matrix protein receptor peroxin 5 (PEX5) in a yeast two-hybrid test. Moreover, we show that examples of these proteins - chicken lysozyme, human tyrosinase and the yeast mitochondrial ribosomal protein L2 (encoded by MRP7) - are imported into peroxisomes in vivo if their original sorting signals are disguised. We also show that even prokaryotic proteins can contain peroxisomal targeting sequences.

Conclusions

Thus, functional localization signals can evolve in unrelated protein sequences as a result of neutral mutations, and subcellular targeting is hierarchically organized, with signal accessibility playing a decisive role. The occurrence of silent functional motifs in unrelated proteins is important for the development of sequence-based function prediction tools and the interpretation of their results. Silent functional signals have the potential to acquire importance in future evolutionary scenarios and in pathological conditions.     

Chemical and enzymatic fractionation of cell walls from Fucales: insights into the structure of the extracellular matrix of brown algae

Background and Aims

Brown algae are photosynthetic multicellular marine organisms evolutionarily distant from land plants, with a distinctive cell wall. They feature carbohydrates shared with plants (cellulose), animals (fucose-containing sulfated polysaccharides, FCSPs) or bacteria (alginates). How these components are organized into a three-dimensional extracellular matrix (ECM) still remains unclear. Recent molecular analysis of the corresponding biosynthetic routes points toward a complex evolutionary history that shaped the ECM structure in brown algae.

Methods

Exhaustive sequential extractions and composition analyses of cell wall material from various brown algae of the order Fucales were performed. Dedicated enzymatic degradations were used to release and identify cell wall partners. This approach was complemented by systematic chromatographic analysis to study polymer interlinks further. An additional structural assessment of the sulfated fucan extracted from Himanthalia elongata was made.

Key Results

The data indicate that FCSPs are tightly associated with proteins and cellulose within the walls. Alginates are associated with most phenolic compounds. The sulfated fucans from H. elongata were shown to have a regular α-(1→3) backbone structure, while an alternating α-(1→3), (1→4) structure has been described in some brown algae from the order Fucales.

Conclusions

The data provide a global snapshot of the cell wall architecture in brown algae, and contribute to the understanding of the structure–function relationships of the main cell wall components. Enzymatic cross-linking of alginates by phenols may regulate the strengthening of the wall, and sulfated polysaccharides may play a key role in the adaptation to osmotic stress. The emergence and evolution of ECM components is further discussed in relation to the evolution of multicellularity in brown algae.     

Trends in genome dynamics among major orders of insects revealed through variations in protein families

Background

Insects belong to a class that accounts for the majority of animals on earth. With over one million identified species, insects display a huge diversity and occupy extreme environments. At present, there are dozens of fully sequenced insect genomes that cover a range of habitats, social behavior and morphologies. In view of such diverse collection of genomes, revealing evolutionary trends and charting functional relationships of proteins remain challenging.

Results

We analyzed the relatedness of 17 complete proteomes representative of proteomes from insects including louse, bee, beetle, ants, flies and mosquitoes, as well as an out-group from the crustaceans. The analyzed proteomes mostly represented the orders of Hymenoptera and Diptera. The 287,405 protein sequences from the 18 proteomes were automatically clustered into 20,933 families, including 799 singletons. A comprehensive analysis based on statistical considerations identified the families that were significantly expanded or reduced in any of the studied organisms. Among all the tested species, ants are characterized by an exceptionally high rate of family gain and loss. By assigning annotations to hundreds of species-specific families, the functional diversity among species and between the major clades (Diptera and Hymenoptera) is revealed. We found that many species-specific families are associated with receptor signaling, stress-related functions and proteases. The highest variability among insects associates with the function of transposition and nucleic acids processes (collectively coined TNAP). Specifically, the wasp and ants have an order of magnitude more TNAP families and proteins relative to species that belong to Diptera (mosquitoes and flies).

Conclusions

An unsupervised clustering methodology combined with a comparative functional analysis unveiled proteomic signatures in the major clades of winged insects. We propose that the expansion of TNAP families in Hymenoptera potentially contributes to the accelerated genome dynamics that characterize the wasp and ants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1771-2) contains supplementary material, which is available to authorized users.     

A deep insight into the sialotranscriptome of the gulf coast tick, Amblyomma maculatum

Background

Saliva of blood sucking arthropods contains compounds that antagonize their hosts'' hemostasis, which include platelet aggregation, vasoconstriction and blood clotting; saliva of these organisms also has anti-inflammatory and immunomodullatory properties. Perhaps because hosts mount an active immune response against these compounds, the diversity of these compounds is large even among related blood sucking species. Because of these properties, saliva helps blood feeding as well as help the establishment of pathogens that can be transmitted during blood feeding.

Methodology/Principal Findings

We have obtained 1,626,969 reads by pyrosequencing a salivary gland cDNA library from adult females Amblyomma maculatum ticks at different times of feeding. Assembly of this data produced 72,441 sequences larger than 149 nucleotides from which 15,914 coding sequences were extracted. Of these, 5,353 had >75% coverage to their best match in the non-redundant database from the National Center for Biotechnology information, allowing for the deposition of 4,850 sequences to GenBank. The annotated data sets are available as hyperlinked spreadsheets. Putative secreted proteins were classified in 133 families, most of which have no known function.

Conclusions/Significance

This data set of proteins constitutes a mining platform for novel pharmacologically active proteins and for uncovering vaccine targets against A. maculatum and the diseases they carry.     

Analyses of phylogeny, evolution, conserved sequences and genome-wide expression of the ICK/KRP family of plant CDK inhibitors

Background and Aims

The cell cycle is controlled by cyclin-dependent kinases (CDKs), and CDK inhibitors are major regulators of their activities. The ICK/KRP family of CDK inhibitors has been reported in several plants, with seven members in arabidopsis; however, the phylogenetic relationship among members in different species is unknown. Also, there is a need to understand how these genes and proteins are regulated. Furthermore, little information is available on the functional differences among ICK/KRP family members.

Methods

We searched publicly available databases and identified over 120 unique ICK/KRP protein sequences from more than 60 plant species. Phylogenetic analysis was performed using 101 full-length sequences from 40 species and intron–exon organization of ICK/KRP genes in model species. Conserved sequences and motifs were analysed using ICK/KRP protein sequences from arabidopsis (Arabidopsis thaliana), rice (Orysa sativa) and poplar (Populus trichocarpa). In addition, gene expression was examined using microarray data from arabidopsis, rice and poplar, and further analysed by RT-PCR for arabidopsis.

Key Results and Conclusions

Phylogenetic analysis showed that plant ICK/KRP proteins can be grouped into three major classes. Whereas the C-class contains sequences from dicotyledons, monocotyledons and gymnosperms, the A- and B-classes contain only sequences from dicotyledons or monocotyledons, respectively, suggesting that the A- and B-classes might have evolved from the C-class. This classification is also supported by exon–intron organization. Genes in the A- and B- classes have four exons, whereas genes in the C-class have only three exons. Analysis of sequences from arabidopsis, rice and poplar identified conserved sequence motifs, some of which had not been described previously, and putative functional sites. The presence of conserved motifs in different family members is consistent with the classification. In addition, gene expression analysis showed preferential expression of ICK/KRP genes in certain tissues. A model has been proposed for the evolution of this gene family in plants.     

Structural relationships in the lysozyme superfamily: significant evidence for glycoside hydrolase signature motifs

Background

Chitin is a polysaccharide that forms the hard, outer shell of arthropods and the cell walls of fungi and some algae. Peptidoglycan is a polymer of sugars and amino acids constituting the cell walls of most bacteria. Enzymes that are able to hydrolyze these cell membrane polymers generally play important roles for protecting plants and animals against infection with insects and pathogens. A particular group of such glycoside hydrolase enzymes share some common features in their three-dimensional structure and in their molecular mechanism, forming the lysozyme superfamily.

Results

Besides having a similar fold, all known catalytic domains of glycoside hydrolase proteins of lysozyme superfamily (families and subfamilies GH19, GH22, GH23, GH24 and GH46) share in common two structural elements: the central helix of the all-α domain, which invariably contains the catalytic glutamate residue acting as general-acid catalyst, and a β-hairpin pointed towards the substrate binding cleft. The invariant β-hairpin structure is interestingly found to display the highest amino acid conservation in aligned sequences of a given family, thereby allowing to define signature motifs for each GH family. Most of such signature motifs are found to have promising performances for searching sequence databases. Our structural analysis further indicates that the GH motifs participate in enzymatic catalysis essentially by containing the catalytic water positioning residue of inverting mechanism.

Conclusions

The seven families and subfamilies of the lysozyme superfamily all have in common a β-hairpin structure which displays a family-specific sequence motif. These GH β-hairpin motifs contain potentially important residues for the catalytic activity, thereby suggesting the participation of the GH motif to catalysis and also revealing a common catalytic scheme utilized by enzymes of the lysozyme superfamily.     

Differences in the number of intrinsically disordered regions between yeast duplicated proteins, and their relationship with functional divergence

Background

Intrinsically disordered regions are enriched in short interaction motifs that play a critical role in many protein-protein interactions. Since new short interaction motifs may easily evolve, they have the potential to rapidly change protein interactions and cellular signaling. In this work we examined the dynamics of gain and loss of intrinsically disordered regions in duplicated proteins to inspect if changes after genome duplication can create functional divergence. For this purpose we used Saccharomyces cerevisiae and the outgroup species Lachancea kluyveri.

Principal Findings

We find that genes duplicated as part of a genome duplication (ohnologs) are significantly more intrinsically disordered than singletons (p<2.2^e-16, Wilcoxon), reflecting a preference for retaining intrinsically disordered proteins in duplicate. In addition, there have been marked changes in the extent of intrinsic disorder following duplication. A large number of duplicated genes have more intrinsic disorder than their L. kluyveri ortholog (29% for duplicates versus 25% for singletons) and an even greater number have less intrinsic disorder than the L. kluyveri ortholog (37% for duplicates versus 25% for singletons). Finally, we show that the number of physical interactions is significantly greater in the more intrinsically disordered ohnolog of a pair (p = 0.003, Wilcoxon).

Conclusion

This work shows that intrinsic disorder gain and loss in a protein is a mechanism by which a genome can also diverge and innovate. The higher number of interactors for proteins that have gained intrinsic disorder compared with their duplicates may reflect the acquisition of new interaction partners or new functional roles.     

Identifying Cis-Regulatory Sequences by Word Profile Similarity

Background

Recognizing regulatory sequences in genomes is a continuing challenge, despite a wealth of available genomic data and a growing number of experimentally validated examples.

Methodology/Principal Findings

We discuss here a simple approach to search for regulatory sequences based on the compositional similarity of genomic regions and known cis-regulatory sequences. This method, which is not limited to searching for predefined motifs, recovers sequences known to be under similar regulatory control. The words shared by the recovered sequences often correspond to known binding sites. Furthermore, we show that although local word profile clustering is predictive for the regulatory sequences involved in blastoderm segmentation, local dissimilarity is a more universal feature of known regulatory sequences in Drosophila.

Conclusions/Significance

Our method leverages sequence motifs within a known regulatory sequence to identify co-regulated sequences without explicitly defining binding sites. We also show that regulatory sequences can be distinguished from surrounding sequences by local sequence dissimilarity, a novel feature in identifying regulatory sequences across a genome. Source code for WPH-finder is available for download at http://rana.lbl.gov/downloads/wph.tar.gz.     

Efficacy of Mycobacterium indicus pranii Immunotherapy as an Adjunct to Chemotherapy for Tuberculosis and Underlying Immune Responses in the Lung

Background

The 9-month-long chemotherapy of tuberculosis often results in poor compliance and emergence of drug-resistant strains. So, improved therapeutic strategy is urgently needed. Immunotherapy could be beneficial for the effective management of the disease. Previously we showed the protective efficacy of Mycobacterium indicus pranii (MIP) when given as prophylactic vaccine in animal models of tuberculosis.

Methods

We sought to investigate whether MIP can be used as an adjunct to the chemotherapy in guinea pig models of tuberculosis. Efficacy of MIP was evaluated when given subcutaneously or by aerosol.

Results

MIP-therapy as an adjunct to the chemotherapy was found to be effective in accelerating bacterial killing and improving organ pathology. MIP-immunotherapy resulted in higher numbers of activated antigen-presenting cells and lymphocytes in the infected lungs and also modulated the granulomatous response. Early increase in protective Th1 immune response was observed in the immunotherapy group. Following subsequent doses of MIP, decrease in the inflammatory response and increase in the immunosuppressive response was observed, which resulted in the improvement of lung pathology.

Conclusion

MIP immunotherapy is a valuable adjunct to chemotherapy for tuberculosis. Aerosol route of immunotherapy can play a crucial role for inducing immediate local immune response in the lung.     

Analysis of the key elements of FFAT-like motifs identifies new proteins that potentially bind VAP on the ER, including two AKAPs and FAPP2

Background

Two phenylalanines (FF) in an acidic tract (FFAT)-motifs were originally described as having seven elements: an acidic flanking region followed by 6 residues (EFFDA–E). Such motifs are found in several lipid transfer protein (LTP) families, and they interact with a protein on the cytosolic face of the ER called vesicle-associated membrane protein-associated protein (VAP). Mutation of which causes ER stress and motor neuron disease, making it important to determine which proteins bind VAP. Among other proteins that bind VAP, some contain FFAT-like motifs that are missing one or more of the seven elements. Defining how much variation is tolerated in FFAT-like motifs is a preliminary step prior to the identification of the full range of VAP interactors.

Results

We used a quantifiable in vivo system that measured ER targeting in a reporter yeast strain that over-expressed VAP to study the effect of substituting different elements of FFAT-like motifs in turn. By defining FFAT-like motifs more widely than before, we found them in novel proteins the functions of which had not previously been directly linked to the ER, including: two PKA anchoring proteins, AKAP220 and AKAP110; a family of plant LTPs; and the glycolipid LTP phosphatidylinositol-four-phosphate adaptor-protein-2 (FAPP-2).

Conclusion

All of the seven essential elements of a FFAT motif tolerate variation, and weak targeting to the ER via VAP is still detected if two elements are substituted. In addition to the strong FFAT motifs already known, there are additional proteins with weaker FFAT-like motifs, which might be functionally important VAP interactors.     

Efficient prediction of human protein-protein interactions at a global scale

Background

Our knowledge of global protein-protein interaction (PPI) networks in complex organisms such as humans is hindered by technical limitations of current methods.

Results

On the basis of short co-occurring polypeptide regions, we developed a tool called MP-PIPE capable of predicting a global human PPI network within 3 months. With a recall of 23% at a precision of 82.1%, we predicted 172,132 putative PPIs. We demonstrate the usefulness of these predictions through a range of experiments.

Conclusions

The speed and accuracy associated with MP-PIPE can make this a potential tool to study individual human PPI networks (from genomic sequences alone) for personalized medicine.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0383-1) contains supplementary material, which is available to authorized users.     

Comparison of Different Magnetic Resonance Cholangiography Techniques in Living Liver Donors Including Gd-EOB-DTPA Enhanced T1-Weighted Sequences

Objectives

Preoperative evaluation of potential living liver donors (PLLDs) includes the assessment of the biliary anatomy to avoid postoperative complications. Aim of this study was to compare T2-weighted (T2w) and Gd-EOB-DTPA enhanced T1-weighted (T1w) magnetic resonance cholangiography (MRC) techniques in the evaluation of PLLDs.

Materials and Methods

30 PLLDs underwent MRC on a 1.5 T Magnetom Avanto (Siemens, Erlangen, Germany) using (A) 2D T2w HASTE (Half Fourier Acquisition Single Shot Turbo Spin Echo) fat saturated (fs) in axial plane, (B) 2D T2w HASTE fs thick slices in coronal plane, (C) free breathing 3D T2w TSE (turbo spin echo) RESTORE (high-resolution navigator corrected) plus (D) maximum intensity projections (MIPs), (E) T2w SPACE (sampling perfection with application optimized contrasts using different flip angle evolutions) plus (F) MIPs and (G) T2w TSE BLADE as well as Gd-EOB-DTPA T1w images without (G) and with (H) inversion recovery. Contrast enhanced CT cholangiography served as reference imaging modality. Two independent reviewers evaluated the biliary tract anatomy on a 5-point scale subjectively and objectively. Data sets were compared using a Mann-Whitney-U-test. Kappa values were also calculated.

Results

Source images and maximum intensity projections of 3D T2w TSE sequences (RESTORE and SPACE) proved to be best for subjective and objective evaluation directly followed by 2D HASTE sequences. Interobserver variabilities were good to excellent (k = 0.622–0.804).

Conclusions

3D T2w sequences are essential for preoperative biliary tract evaluation in potential living liver donors. Furthermore, our results underline the value of different MRCP sequence types for the evaluation of the biliary anatomy in PLLDs including Gd-EOB-DTPA enhanced T1w MRC.