Integrated Omic Approaches Reveal Molecular Mechanisms of Tolerance during Soybean and Meloidogyne incognita Interactions

The root-knot nematode (RKN), Meloidogyne incognita, is a devastating soybean pathogen worldwide. The use of resistant cultivars is the most effective method to prevent economic losses caused by RKNs. To elucidate the mechanisms involved in resistance to RKN, we determined the proteome and transcriptome profiles from roots of susceptible (BRS133) and highly tolerant (PI 595099) Glycine max genotypes 4, 12, and 30 days after RKN infestation. After in silico analysis, we described major defense molecules and mechanisms considered constitutive responses to nematode infestation, such as mTOR, PI3K-Akt, relaxin, and thermogenesis. The integrated data allowed us to identify protein families and metabolic pathways exclusively regulated in tolerant soybean genotypes. Among them, we highlighted the phenylpropanoid pathway as an early, robust, and systemic defense process capable of controlling M. incognita reproduction. Associated with this metabolic pathway, 29 differentially expressed genes encoding 11 different enzymes were identified, mainly from the flavonoid and derivative pathways. Based on differential expression in transcriptomic and proteomic data, as well as in the expression profile by RT-qPCR, and previous studies, we selected and overexpressed the GmPR10 gene in transgenic tobacco to assess its protective effect against M. incognita. Transgenic plants of the T₂ generation showed up to 58% reduction in the M. incognita reproduction factor. Finally, data suggest that GmPR10 overexpression can be effective against the plant parasitic nematode M. incognita, but its mechanism of action remains unclear. These findings will help develop new engineered soybean genotypes with higher performance in response to RKN infections.

Defining the combined stress response in wild Arachis

Nematodes and drought are major constraints in tropical agriculture and often occur simultaneously. Plant responses to these stresses are complex and require crosstalk between biotic and abiotic signaling pathways. In this study, we explored the transcriptome data of wild Arachis species subjected to drought (A-metaDEG) and the root-knot nematode Meloidogyne arenaria (B-metaDEG) via meta-analysis, to identify core-stress responsive genes to each individual and concurrent stresses in these species. Transcriptome analysis of a nematode/drought bioassay (cross-stress) showed that the set of stress responsive DEGs to concurrent stress is distinct from those resulting from overlapping A- and B-metaDEGs, indicating a specialized and unique response to combined stresses in wild Arachis. Whilst individual biotic and abiotic stresses elicit hormone-responsive genes, most notably in the jasmonic and abscisic acid pathways, combined stresses seem to trigger mainly the ethylene hormone pathway. The overexpression of a cross-stress tolerance candidate gene identified here, an endochitinase-encoding gene (AsECHI) from Arachis stenosperma, reduced up to 30% of M. incognita infection and increased post-drought recovery in Arabidopsis plants submitted to both stresses. The elucidation of the network of cross-stress responsive genes in Arachis contributes to better understanding the complex regulation of biotic and abiotic responses in plants facilitating more adequate crop breeding for combined stress tolerance.

Transcriptome Analysis and Knockdown of the Juvenile Hormone Esterase Gene Reveal Abnormal Feeding Behavior in the Sugarcane Giant Borer

The sugarcane giant borer (SGB), Telchin licus licus, is a pest that has strong economic relevance for sugarcane producers. Due to the endophytic behavior of the larva, current methods of management are inefficient. A promising biotechnological management option has been proposed based on RNA interference (RNAi), a process that uses molecules of double-stranded RNA (dsRNA) to specifically knock down essential genes and reduce insect survival. The selection of suitable target genes is often supported by omic sciences. Studies have shown that genes related to feeding adaptation processes are good candidates to be targeted by RNAi for pest management. Among those genes, esterases are highlighted because of their impact on insect development. In this study, the objective was to evaluate the transcriptome responses of the SGB’s gut in order to provide curated data of genes that could be used for pest management by RNAi in future studies. Further, we validated the function of an esterase-coding gene and its potential as a target for RNAi-based control. We sequenced the gut transcriptome of SGB larvae by Illumina HiSeq and evaluated its gene expression profiles in response to different diets (sugarcane stalk and artificial diet). We obtained differentially expressed genes (DEGs) involved in detoxification, digestion, and transport, which suggest a generalist mechanism of adaptation in SGB larvae. Among the DEGs, was identified and characterized a candidate juvenile hormone esterase gene (Tljhe). We knocked down the Tljhe gene by oral delivery of dsRNA molecules and evaluated gene expression in the gut. The survival and nutritional parameters of the larvae were measured along the developmental cycle of treated insects. We found that the gene Tljhe acts as a regulator of feeding behavior. The knockdown of Tljhe triggered a forced starvation state in late larval instars that significantly reduced the fitness of the larvae. However, the mechanism of action of this gene remains unclear, and the correlation between the expression of Tljhe and the levels of juvenile hormone (JH) metabolites in the hemolymph of the SGB must be assessed in future research.

SERK genes identification and expression analysis during somatic embryogenesis and sporogenesis of sexual and apomictic Brachiaria brizantha (Syn. Urochloa brizantha)

In Brachiaria brizantha BbrizSERK1, BbrizSERK2 and BbrizSERK3 were identified. SERK expression marks somatic embryogenesis, sexual MMC, and sexual and apomictic PMC. BbrizSERK3 might have a regulatory role in reproductive development. Somatic embryogenesis receptor-like kinase (SERK) consists of plasma membrane receptor genes that have been characterized in various species, associated with several aspects of plant development, including reproduction. SERK genes are involved in anther development and in early embryo development in sexual and asexual seed formation. To comprehend the complexity of the SERK genes and their function in Brachiaria reproduction, we performed a homology-based search in a genomic database of a sexual B. brizantha and identified sequences of three SERK genes, BbrizSERK1, BbrizSERK2, and BbrizSERK3. RNASeq data showed equivalent abundance of BbrizSERK1 and BbrizSERK2 transcripts in ovaries at early megasporogenesis of sexuals and apomicts, while BbrizSERK3 transcripts were more abundant in ovaries of sexuals than in apomicts. BbrizSERK3 results in three coding sequences due to alternative splicing, among them Variant 1 results in a protein with all the predicted domains of a SERK. BbrizSERK transcripts were detected in male reproductive tissues of both sexual and apomictic plants, suggesting a role in controlling anther development. BbrizSERK transcripts were detected early in ovule development, in the integuments, and in the megaspore mother cell of the sexual plant, but not in the cells that give rise to apomictic embryo sacs, suggesting a role in female reproductive development of sexuals. This paper provides evidences that SERK genes plays a role in the onset and establishment of somatic embryogenesis and in the reproductive development of B. brizantha and suggests a distinct role of BbrizSERK in apomixis initiation.

Comparative gut transcriptome analysis of Diatraea saccharalis in response to the dietary source

The sugarcane borer (Diatraea saccharalis, Fabricius, 1794) is a devastating pest that causes millions of dollars of losses each year to sugarcane producers by reducing sugar and ethanol yields. The control of this pest is difficult due to its endophytic behavior and rapid development. Pest management through biotechnological approaches has emerged in recent years as an alternative to currently applied methods. Genetic information about the target pests is often required to perform biotechnology-based management. The genomic and transcriptomic data for D. saccharalis are very limited. Herein, we report a tissue-specific transcriptome of D. saccharalis larvae and a differential expression analysis highlighting the physiological characteristics of this pest in response to two different diets: sugarcane and an artificial diet. Sequencing was performed on the Illumina HiSeq 2000 platform, and a de novo assembly was generated. A total of 27,626 protein-coding unigenes were identified, among which 1,934 sequences were differentially expressed between treatments. Processes such as defence, digestion, detoxification, signaling, and transport were highly represented among the differentially expressed genes (DEGs). Furthermore, seven aminopeptidase genes were identified as candidates to encode receptors of Cry proteins, which are toxins of Bacillus thuringiensis used to control lepidopteran pests. Since plant-insect interactions have produced a considerable number of adaptive responses in hosts and herbivorous insects, the success of phytophagous insects relies on their ability to overcome challenges such as the response to plant defences and the intake of nutrients. In this study, we identified metabolic pathways and specific genes involved in these processes. Thus, our data strongly contribute to the knowledge advancement of insect transcripts, which can be a source of target genes for pest management.

Transcriptome and gene expression analysis of three developmental stages of the coffee berry borer, Hypothenemus hampei

Coffee production is a global industry valued at approximately 173 billion US dollars. One of the main challenges facing coffee production is the management of the coffee berry borer (CBB), Hypothenemus hampei, which is considered the primary arthropod pest of coffee worldwide. Current control strategies are inefficient for CBB management. Although biotechnological alternatives, including RNA interference (RNAi), have been proposed in recent years to control insect pests, characterizing the genetics of the target pest is essential for the successful application of these emerging technologies. In this study, we employed RNA-seq to obtain the transcriptome of three developmental stages of the CBB (larva, female and male) to increase our understanding of the CBB life cycle in relation to molecular features. The CBB transcriptome was sequenced using Illumina Hiseq and assembled de novo. Differential gene expression analysis was performed across the developmental stages. The final assembly produced 29,434 unigenes, of which 4,664 transcripts were differentially expressed. Genes linked to crucial physiological functions, such as digestion and detoxification, were determined to be tightly regulated between the reproductive and nonreproductive stages of CBB. The data obtained in this study help to elucidate the critical roles that several genes play as regulatory elements in CBB development.

Identification and expression profile of odorant-binding proteins in Halyomorpha halys (Hemiptera: Pentatomidae)

The brown marmorated stink bug, Halyomorpha halys, is a devastating invasive species in the USA. Similar to other insects, olfaction plays an important role in its survival and reproduction. As odorant-binding proteins (OBPs) are involved in the initial semiochemical recognition steps, we used RNA-Sequencing (RNA-Seq) to identify OBPs in its antennae, and studied their expression pattern in different body parts under semiochemical stimulation by either aggregation or alarm pheromone or food odorants. Thirty full-length putative HhalOBPs were identified, corresponding to 22 'classic' OBPs and eight 'Plus-C' OBPs. The similarity amongst them ranged from 4.95-70.92%, and with another 325 hemipteran OBPs similarity ranged from 1.94-91.51%, the highest levels being with other stink bug OBPs. Phylogenetic analysis confirmed the monophyly of seven groups of stink bug and other hemipteran OBPs. All 30 HhalOBPs were expressed and about 2/3 were expressed primarily in antennae. The expression of 21 HhalOBPs was higher in the antennae under alarm pheromone stimulus, indicating that multiple OBPs may be responding to this pheromone. Two were highest in antennae under aggregation pheromone stimulus. These findings should provide a basis for understanding the physiological functions of HhalOBPs and the chemosensory perception of this pest, which may help to uncover new control targets for behavioural interference.

Root Transcriptome Analysis of Wild Peanut Reveals Candidate Genes for Nematode Resistance

Wild peanut relatives (Arachis spp.) are genetically diverse and were adapted to a range of environments during the evolution course, constituting an important source of allele diversity for resistance to biotic and abiotic stresses. The wild diploid A. stenosperma harbors high levels of resistance to a variety of pathogens, including the root-knot nematode (RKN) Meloidogyne arenaria, through the onset of the Hypersensitive Response (HR). In order to identify genes and regulators triggering this defense response, a comprehensive root transcriptome analysis during the first stages of this incompatible interaction was conducted using Illumina Hi-Seq. Overall, eight cDNA libraries were produced generating 28.2 GB, which were de novo assembled into 44,132 contigs and 37,882 loci. Differentially expressed genes (DEGs) were identified and clustered according to their expression profile, with the majority being downregulated at 6 DAI, which coincides with the onset of the HR. Amongst these DEGs, 27 were selected for further qRT-PCR validation allowing the identification of nematode-responsive candidate genes that are putatively related to the resistance response. Those candidates are engaged in the salycilic (NBS-LRR, lipocalins, resveratrol synthase) and jasmonic (patatin, allene oxidase cyclase) acids pathways, and also related to hormonal balance (auxin responsive protein, GH3) and cellular plasticity and signaling (tetraspanin, integrin, expansin), with some of them showing contrasting expression behavior between Arachis RKN-resistant and susceptible genotypes. As these candidate genes activate different defensive signaling systems, the genetic (HR) and the induced resistance (IR), their pyramidding in one genotype via molecular breeding or transgenic strategy might contribute to a more durable resistance, thus improving the long-term control of RKN in peanut.

Transcriptome-Based Identification of Highly Similar Odorant-Binding Proteins among Neotropical Stink Bugs and Their Egg Parasitoid

Olfaction plays a fundamental role in insect survival through resource location and intra and interspecific communications. We used RNA-Seq to analyze transcriptomes for odorant-binding proteins (OBPs) from major stink bug pest species in Brazil, Euschistus heros, Chinavia ubica, and Dichelops melacanthus, and from their egg parasitoid, Telenomus podisi. We identified 23 OBPs in E. heros, 25 OBPs in C. ubica, 9 OBPs in D. melacanthus, and 7 OBPs in T. podisi. The deduced amino acid sequences of the full-length OBPs had low intraspecific similarity, but very high similarity between two pairs of OBPs from E. heros and C. ubica (76.4 and 84.0%) and between two pairs of OBPs from the parasitoid and its preferred host E. heros (82.4 and 88.5%), confirmed by a high similarity of their predicted tertiary structures. The similar pairs of OBPs from E. heros and C. ubica may suggest that they have derived from a common ancestor, and retain the same biological function to bind a ligand perceived or produced in both species. The T. podisi OBPs similar to E. heros were not orthologous to any known hymenopteran OBPs, and may have evolved independently and converged to the host OBPs, providing a possible basis for the host location of T. podisi using E. heros semiochemical cues.

Transcriptome Profiling of Wild Arachis from Water-Limited Environments Uncovers Drought Tolerance Candidate Genes

Peanut (Arachis hypogaea L.) is an important legume cultivated mostly in drought-prone areas where its productivity can be limited by water scarcity. The development of more drought-tolerant varieties is, therefore, a priority for peanut breeding programs worldwide. In contrast to cultivated peanut, wild relatives have a broader genetic diversity and constitute a rich source of resistance/tolerance alleles to biotic and abiotic stresses. The present study takes advantage of this diversity to identify drought-responsive genes by analyzing the expression profile of two wild species, Arachis duranensis and Arachis magna (AA and BB genomes, respectively), in response to progressive water deficit in soil. Data analysis from leaves and roots of A. duranensis (454 sequencing) and A. magna (suppression subtractive hybridization (SSH)) stressed and control complementary DNA (cDNA) libraries revealed several differentially expressed genes in silico, and 44 of them were selected for further validation by quantitative RT-PCR (qRT-PCR). This allowed the identification of drought-responsive candidate genes, such as Expansin, Nitrilase, NAC, and bZIP transcription factors, displaying significant levels of differential expression during stress imposition in both species. This is the first report on identification of differentially expressed genes under drought stress and recovery in wild Arachis species. The generated transcriptome data, besides being a valuable resource for gene discovery, will allow the characterization of new alleles and development of molecular markers associated with drought responses in peanut. These together constitute important tools for the peanut breeding program and also contribute to a better comprehension of gene modulation in response to water deficit and rehydration.

The genome sequence of Pseudoplusia includens single nucleopolyhedrovirus and an analysis of p26 gene evolution in the baculoviruses

BACKGROUND

Pseudoplusia includens single nucleopolyhedrovirus (PsinSNPV-IE) is a baculovirus recently identified in our laboratory, with high pathogenicity to the soybean looper, Chrysodeixis includens (Lepidoptera: Noctuidae) (Walker, 1858). In Brazil, the C. includens caterpillar is an emerging pest and has caused significant losses in soybean and cotton crops. The PsinSNPV genome was determined and the phylogeny of the p26 gene within the family Baculoviridae was investigated.

RESULTS

The complete genome of PsinSNPV was sequenced (Roche 454 GS FLX - Titanium platform), annotated and compared with other Alphabaculoviruses, displaying a genome apparently different from other baculoviruses so far sequenced. The circular double-stranded DNA genome is 139,132 bp in length, with a GC content of 39.3 % and contains 141 open reading frames (ORFs). PsinSNPV possesses the 37 conserved baculovirus core genes, 102 genes found in other baculoviruses and 2 unique ORFs. Two baculovirus repeat ORFs (bro) homologs, bro-a (Psin33) and bro-b (Psin69), were identified and compared with Chrysodeixis chalcites nucleopolyhedrovirus (ChchNPV) and Trichoplusia ni single nucleopolyhedrovirus (TnSNPV) bro genes and showed high similarity, suggesting that these genes may be derived from an ancestor common to these viruses. The homologous repeats (hrs) are absent from the PsinSNPV genome, which is also the case in ChchNPV and TnSNPV. Two p26 gene homologs (p26a and p26b) were found in the PsinSNPV genome. P26 is thought to be required for optimal virion occlusion in the occlusion bodies (OBs), but its function is not well characterized. The P26 phylogenetic tree suggests that this gene was obtained from three independent acquisition events within the Baculoviridae family. The presence of a signal peptide only in the PsinSNPV p26a/ORF-20 homolog indicates distinct function between the two P26 proteins.

CONCLUSIONS

PsinSNPV has a genomic sequence apparently different from other baculoviruses sequenced so far. The complete genome sequence of PsinSNPV will provide a valuable resource, contributing to studies on its molecular biology and functional genomics, and will promote the development of this virus as an effective bioinsecticide.

Analysis of the leaf transcriptome of Musa acuminata during interaction with Mycosphaerella musicola: gene assembly, annotation and marker development

BACKGROUND

Although banana (Musa sp.) is an important edible crop, contributing towards poverty alleviation and food security, limited transcriptome datasets are available for use in accelerated molecular-based breeding in this genus. 454 GS-FLX Titanium technology was employed to determine the sequence of gene transcripts in genotypes of Musa acuminata ssp. burmannicoides Calcutta 4 and M. acuminata subgroup Cavendish cv. Grande Naine, contrasting in resistance to the fungal pathogen Mycosphaerella musicola, causal organism of Sigatoka leaf spot disease. To enrich for transcripts under biotic stress responses, full length-enriched cDNA libraries were prepared from whole plant leaf materials, both uninfected and artificially challenged with pathogen conidiospores.

RESULTS

The study generated 846,762 high quality sequence reads, with an average length of 334 bp and totalling 283 Mbp. De novo assembly generated 36,384 and 35,269 unigene sequences for M. acuminata Calcutta 4 and Cavendish Grande Naine, respectively. A total of 64.4% of the unigenes were annotated through Basic Local Alignment Search Tool (BLAST) similarity analyses against public databases.Assembled sequences were functionally mapped to Gene Ontology (GO) terms, with unigene functions covering a diverse range of molecular functions, biological processes and cellular components. Genes from a number of defense-related pathways were observed in transcripts from each cDNA library. Over 99% of contig unigenes mapped to exon regions in the reference M. acuminata DH Pahang whole genome sequence. A total of 4068 genic-SSR loci were identified in Calcutta 4 and 4095 in Cavendish Grande Naine. A subset of 95 potential defense-related gene-derived simple sequence repeat (SSR) loci were validated for specific amplification and polymorphism across M. acuminata accessions. Fourteen loci were polymorphic, with alleles per polymorphic locus ranging from 3 to 8 and polymorphism information content ranging from 0.34 to 0.82.

CONCLUSIONS

A large set of unigenes were characterized in this study for both M. acuminata Calcutta 4 and Cavendish Grande Naine, increasing the number of public domain Musa ESTs. This transcriptome is an invaluable resource for furthering our understanding of biological processes elicited during biotic stresses in Musa. Gene-based markers will facilitate molecular breeding strategies, forming the basis of genetic linkage mapping and analysis of quantitative trait loci.

Characterization of soil bacterial assemblies in Brazilian savanna-like vegetation reveals acidobacteria dominance

The Brazilian Cerrado is the second largest biome in Brazil and is considered a biodiversity hotspot. In this work, we compared the bacterial communities in Cerrado soil associated with four types of native vegetation (Cerrado Denso, Cerrado sensu stricto, Campo Sujo, and Mata de Galeria) by ribosomal RNA intergenic spacer analysis, terminal fragment restriction length polymorphism and pyrosequencing. The fingerprinting results were very similar. The bacterial communities of Cerrado Denso and Cerrado sensu stricto grouped together and were distinct from those in Campo Sujo and Mata de Galeria. Pyrosequencing generated approximately 40,000 16S rRNA gene sequences per sample and allowed the identification of 17 phyla in soil samples under Cerrado vegetation. Acidobacteria were dominant in all areas studied with a relative frequency of 40-47 %, followed closely by Proteobacteria accounting for 34-40 % of the sequences. Results from all molecular techniques used suggested that the bacterial communities of Cerrado sensu stricto and Cerrado Denso are very similar to each other, while Campo Sujo forms a separate group, and Mata de Galeria is the most distinct with higher species richness. This is the first extensive study of native Cerrado soil microbiota, an important but endangered biome.

Analysis of binding properties and specificity through identification of the interface forming residues (IFR) for serine proteases in silico docked to different inhibitors

Background

Enzymes belonging to the same super family of proteins in general operate on variety of substrates and are inhibited by wide selection of inhibitors. In this work our main objective was to expand the scope of studies that consider only the catalytic and binding pocket amino acids while analyzing enzyme specificity and instead, include a wider category which we have named the Interface Forming Residues (IFR). We were motivated to identify those amino acids with decreased accessibility to solvent after docking of different types of inhibitors to sub classes of serine proteases and then create a table (matrix) of all amino acid positions at the interface as well as their respective occupancies. Our goal is to establish a platform for analysis of the relationship between IFR characteristics and binding properties/specificity for bi-molecular complexes.

Results

We propose a novel method for describing binding properties and delineating serine proteases specificity by compiling an exhaustive table of interface forming residues (IFR) for serine proteases and their inhibitors. Currently, the Protein Data Bank (PDB) does not contain all the data that our analysis would require. Therefore, an in silico approach was designed for building corresponding complexes

The IFRs are obtained by "rigid body docking" among 70 structurally aligned, sequence wise non-redundant, serine protease structures with 3 inhibitors: bovine pancreatic trypsin inhibitor (BPTI), ecotine and ovomucoid third domain inhibitor. The table (matrix) of all amino acid positions at the interface and their respective occupancy is created. We also developed a new computational protocol for predicting IFRs for those complexes which were not deciphered experimentally so far, achieving accuracy of at least 0.97.

Conclusions

The serine proteases interfaces prefer polar (including glycine) residues (with some exceptions). Charged residues were found to be uniquely prevalent at the interfaces between the "miscellaneous-virus" subfamily and the three inhibitors. This prompts speculation about how important this difference in IFR characteristics is for maintaining virulence of those organisms.

Our work here provides a unique tool for both structure/function relationship analysis as well as a compilation of indicators detailing how the specificity of various serine proteases may have been achieved and/or could be altered. It also indicates that the interface forming residues which also determine specificity of serine protease subfamily can not be presented in a canonical way but rather as a matrix of alternative populations of amino acids occupying variety of IFR positions.

Expression, purification and use of the soluble domain of Lactobacillus paracasei beta-fructosidase to optimise production of bioethanol from grass fructans

Microbial inulinases find application in food, pharmaceutical and biofuel industries. Here, a novel Lactobacillus paracasei beta-fructosidase was overexpressed as truncated cytosolic protein ((t)fosEp) in Escherichia coli. Purified (t)fosEp was thermostable (10-50 degrees C) with a pH optimum of 5; it showed highest affinity for bacterial levan (beta[2-6] linked fructose) followed by nystose, chicory inulin, 1-kestose (beta[2-1] linkages) and sucrose (K(m) values of 0.5, 15, 15.6, 49 and 398 mM, respectively). Hydrolysis of polyfructose moieties in agriculturally-sourced grass juice (GJ) with (t)fosEp resulted in the release of >13 mg/ml more bioavailable fructose than was measured in untreated GJ. Bioethanol yields from fermentation experiments with Brewer's yeast and GJ+(t)fosEp were >25% higher than those achieved using untreated GJ feedstock (36.5[+/-4.3] and 28.2[+/-2.7]mg ethanol/ml, respectively). This constitutes the first specific study of the potential to ferment ethanol from grass juice and the utility of a novel core domain of beta-fructosidase from L. paracasei.

Screening non-coding RNAs in transcriptomes from neglected species using PORTRAIT: case study of the pathogenic fungus Paracoccidioides brasiliensis

Background

Transcriptome sequences provide a complement to structural genomic information and provide snapshots of an organism's transcriptional profile. Such sequences also represent an alternative method for characterizing neglected species that are not expected to undergo whole-genome sequencing. One difficulty for transcriptome sequencing of these organisms is the low quality of reads and incomplete coverage of transcripts, both of which compromise further bioinformatics analyses. Another complicating factor is the lack of known protein homologs, which frustrates searches against established protein databases. This lack of homologs may be caused by divergence from well-characterized and over-represented model organisms. Another explanation is that non-coding RNAs (ncRNAs) may be caught during sequencing. NcRNAs are RNA sequences that, unlike messenger RNAs, do not code for protein products and instead perform unique functions by folding into higher order structural conformations. There is ncRNA screening software available that is specific for transcriptome sequences, but their analyses are optimized for those transcriptomes that are well represented in protein databases, and also assume that input ESTs are full-length and high quality.

Results

We propose an algorithm called PORTRAIT, which is suitable for ncRNA analysis of transcriptomes from poorly characterized species. Sequences are translated by software that is resistant to sequencing errors, and the predicted putative proteins, along with their source transcripts, are evaluated for coding potential by a support vector machine (SVM). Either of two SVM models may be employed: if a putative protein is found, a protein-dependent SVM model is used; if it is not found, a protein-independent SVM model is used instead. Only ab initio features are extracted, so that no homology information is needed. We illustrate the use of PORTRAIT by predicting ncRNAs from the transcriptome of the pathogenic fungus Paracoccidoides brasiliensis and five other related fungi.

Conclusion

PORTRAIT can be integrated into pipelines, and provides a low computational cost solution for ncRNA detection in transcriptome sequencing projects.

Analysis of non-TIR NBS-LRR resistance gene analogs in Musa acuminata Colla: isolation, RFLP marker development, and physical mapping

BACKGROUND

Many commercial banana varieties lack sources of resistance to pests and diseases, as a consequence of sterility and narrow genetic background. Fertile wild relatives, by contrast, possess greater variability and represent potential sources of disease resistance genes (R-genes). The largest known family of plant R-genes encode proteins with nucleotide-binding site (NBS) and C-terminal leucine-rich repeat (LRR) domains. Conserved motifs in such genes in diverse plant species offer a means for isolation of candidate genes in banana which may be involved in plant defence.

RESULTS

A computational strategy was developed for unbiased conserved motif discovery in NBS and LRR domains in R-genes and homologues in monocotyledonous plant species. Degenerate PCR primers targeting conserved motifs were tested on the wild cultivar Musa acuminata subsp. burmannicoides, var. Calcutta 4, which is resistant to a number of fungal pathogens and nematodes. One hundred and seventy four resistance gene analogs (RGAs) were amplified and assembled into 52 contiguous sequences. Motifs present were typical of the non-TIR NBS-LRR RGA subfamily. A phylogenetic analysis of deduced amino-acid sequences for 33 RGAs with contiguous open reading frames (ORFs), together with RGAs from Arabidopsis thaliana and Oryza sativa, grouped most Musa RGAs within monocotyledon-specific clades. RFLP-RGA markers were developed, with 12 displaying distinct polymorphisms in parentals and F1 progeny of a diploid M. acuminata mapping population. Eighty eight BAC clones were identified in M. acuminata Calcutta 4, M. acuminata Grande Naine, and M. balbisiana Pisang Klutuk Wulung BAC libraries when hybridized to two RGA probes. Multiple copy RGAs were common within BAC clones, potentially representing variation reservoirs for evolution of new R-gene specificities.

CONCLUSION

This is the first large scale analysis of NBS-LRR RGAs in M. acuminata Calcutta 4. Contig sequences were deposited in GenBank and assigned numbers ER935972 - ER936023. RGA sequences and isolated BACs are a valuable resource for R-gene discovery, and in future applications will provide insight into the organization and evolution of NBS-LRR R-genes in the Musa A and B genome. The developed RFLP-RGA markers are applicable for genetic map development and marker assisted selection for defined traits such as pest and disease resistance.

The Star STING server: a multiplatform environment for protein structure analysis

Star STING is the latest version of the STING suite of programs and corresponding database. We report on five important aspects of this package that have acquired some new characteristics, designed to add key advantages to the whole suite: 1) availability for most popular platforms and browsers, 2) introduction of the STING_DB quality assessment, 3) improvement in algorithms for calculation of three STING parameters, 4) introduction of five new STING modules, and 5) expansion of the existing modules. Star STING is freely accessible at: http://sms.cbi.cnptia.embrapa.br/SMS/, http://trantor.bioc.columbia.edu/SMS, http://www.es.embnet.org/SMS/, http://gibk26.bse.kyutech.ac.jp/SMS/ and http://www.ar.embnet.org/SMS.

The Diamond STING server

Diamond STING is a new version of the STING suite of programs for a comprehensive analysis of a relationship between protein sequence, structure, function and stability. We have added a number of new functionalities by both providing more structure parameters to the STING Database and by improving/expanding the interface for enhanced data handling. The integration among the STING components has also been improved. A new key feature is the ability of the STING server to handle local files containing protein structures (either modeled or not yet deposited to the Protein Data Bank) so that they can be used by the principal STING components: (Java)Protein Dossier ((J)PD) and STING Report. The current capabilities of the new STING version and a couple of biologically relevant applications are described here. We have provided an example where Diamond STING identifies the active site amino acids and folding essential amino acids (both previously determined by experiments) by filtering out all but those residues by selecting the numerical values/ranges for a set of corresponding parameters. This is the fundamental step toward a more interesting endeavor-the prediction of such residues. Diamond STING is freely accessible at http://sms.cbi.cnptia.embrapa.br and http://trantor.bioc.columbia.edu/SMS.

Analysis of expressed sequence tags from Musa acuminata ssp. burmannicoides, var. Calcutta 4 (AA) leaves submitted to temperature stresses

In order to discover genes expressed in leaves of Musa acuminata ssp. burmannicoides var. Calcutta 4 (AA), from plants submitted to temperature stress, we produced and characterized two full-length enriched cDNA libraries. Total RNA from plants subjected to temperatures ranging from 5 degrees C to 25 degrees C and from 25 degrees C to 45 degrees C was used to produce a COLD and a HOT cDNA library, respectively. We sequenced 1,440 clones from each library. Following quality analysis and vector trimming, we assembled 2,286 sequences from both libraries into 1,019 putative transcripts, consisting of 217 clusters and 802 singletons, which we denoted Musa acuminata assembled expressed sequence tagged (EST) sequences (MaAES). Of these MaAES, 22.87% showed no matches with existing sequences in public databases. A global analysis of the MaAES data set indicated that 10% of the sequenced cDNAs are present in both cDNA libraries, while 42% and 48% are present only in the COLD or in the HOT libraries, respectively. Annotation of the MaAES data set categorized them into 22 functional classes. Of the 2,286 high-quality sequences, 715 (31.28%) originated from full-length cDNA clones and resulted in a set of 149 genes.

STING Millennium Suite: integrated software for extensive analyses of 3d structures of proteins and their complexes

BACKGROUND

The integration of many aspects of protein/DNA structure analysis is an important requirement for software products in general area of structural bioinformatics. In fact, there are too few software packages on the internet which can be described as successful in this respect. We might say that what is still missing is publicly available, web based software for interactive analysis of the sequence/structure/function of proteins and their complexes with DNA and ligands. Some of existing software packages do have certain level of integration and do offer analysis of several structure related parameters, however not to the extent generally demanded by a user.

RESULTS

We are reporting here about new Sting Millennium Suite (SMS) version which is fully accessible (including for local files at client end), web based software for molecular structure and sequence/structure/function analysis. The new SMS client version is now operational also on Linux boxes and it works with non-public pdb formatted files (structures not deposited at the RCSB/PDB), eliminating earlier requirement for the registration if SMS components were to be used with user's local files. At the same time the new SMS offers some important additions and improvements such as link to ProTherm as well as significant re-engineering of SMS component ConSSeq. Also, we have added 3 new SMS mirror sites to existing network of global SMS servers: Argentina, Japan and Spain.

CONCLUSION

SMS is already established software package and many key data base and software servers worldwide, do offer either a link to, or host the SMS. SMS (Sting Millennium Suite) is web-based publicly available software developed to aid researches in their quest for translating information about the structures of macromolecules into knowledge. SMS allows to a user to interactively analyze molecular structures, cross-referencing visualized information with a correlated one, available across the internet. SMS is already used as a didactic tool by some universities. SMS analysis is now possible on Linux OS boxes and with no requirement for registration when using local files.

JavaProtein Dossier: a novel web-based data visualization tool for comprehensive analysis of protein structure

JavaProtein Dossier ((J)PD) is a new concept, database and visualization tool providing one of the largest collections of the physicochemical parameters describing proteins' structure, stability, function and interaction with other macromolecules. By collecting as many descriptors/parameters as possible within a single database, we can achieve a better use of the available data and information. Furthermore, data grouping allows us to generate different parameters with the potential to provide new insights into the sequence-structure-function relationship. In (J)PD, residue selection can be performed according to multiple criteria. (J)PD can simultaneously display and analyze all the physicochemical parameters of any pair of structures, using precalculated structural alignments, allowing direct parameter comparison at corresponding amino acid positions among homologous structures. In order to focus on the physicochemical (and consequently pharmacological) profile of proteins, visualization tools (showing the structure and structural parameters) also had to be optimized. Our response to this challenge was the use of Java technology with its exceptional level of interactivity. (J)PD is freely accessible (within the Gold Sting Suite) at http://sms.cbi.cnptia.embrapa.br, http://mirrors.rcsb.org/SMS, http://trantor.bioc.columbia.edu/SMS and http://www.es.embnet.org/SMS/ (Option: (Java)Protein Dossier).

STING Contacts: a web-based application for identification and analysis of amino acid contacts within protein structure and across protein interfaces

Amino acid contacts in terms of atomic interactions are essential factors to be considered in the analysis of the structure of a protein and its complexes. Consequently, molecular biologists do require specific tools for the identification and visualization of all such contacts. Graphical contacts (GC) and interface forming residue graphical contacts (IFRgc) presented here, calculate atomic contacts among amino acids based on a table of predefined pairs of the atom types and their distances, and then display them using number of different forms. The inventory of currently listed contact types by GC and IFRgc include hydrogen bonds (in nine different flavors), hydrophobic interactions, charge-charge interactions, aromatic stacking and disulfide bonds. Such extensive catalog of the interactions, representing the forces that govern protein folding, stability and binding, is the key feature of these two applications. GC and IFRgc are part of STING Millennium Suite.

AVAILABILITY

http://sms.cbi.cnptia.embrapa.br/SMS, http://trantor.bioc.columbia.edu/SMS, http://mirrors.rcsb.org//SMS, http://www.es.embnet.org/SMS and http://www.ar.embnet.org/SMS (Options: Graphical Contacts and IFR Graphical Contacts).

ConSSeq: a web-based application for analysis of amino acid conservation based on HSSP database and within context of structure

SUMMARY

A web-based application to analyze protein amino acids conservation-Consensus Sequence (ConSSeq) is presented. ConSSeq graphically represents information about amino acid conservation based on sequence alignments reported in homology-derived structures of proteins. Beyond the relative entropy for each position in the alignment, ConSSeq also presents the consensus sequence and information about the amino acids, which are predominant at each position of the alignment. ConSSeq is part of the STING Millennium Suite and is implemented as a Java Applet.

AVAILABILITY

http://sms.cbi.cnptia.embrapa.br/SMS/STINGm/consseq/, http://trantor.bioc.columbia.edu/SMS/STINGm/consseq/, http://mirrors.rcsb.org//SMS/STINGm/consseq/, http://www.es.embnet.org/SMS/STINGm/consseq/ and http://www.ar.embnet.org/SMS/STINGm/consseq/

Defining 3D residue environment in protein structures using SCORPION and FORMIGA

SUMMARY

Two web-based applications to analyze amino acids three-dimensional (3D) local environment within protein structures-SCORPION and FORMIGA-are presented. SCORPION and FORMIGA produce a graphical presentation for simple statistical data showing the frequency of residue occurrence within a given sphere (defined here as the 3D contacts). The center of that sphere is placed at the Calpha and at the last heavy atom in the side chain of the selected amino acid. Further depth of detail is given in terms of a secondary structure to which the profiled amino acid belongs. Results obtained with those two applications are relevant for estimating the importance of the amino acid 3D local environment for protein folding and stability. Effectively, SCORPION and FORMIGA construct knowledge-based force fields. The difference between SCORPION and FORMIGA is in that the latter operates on protein interfaces, while the former only functions for a single protein chain. Both applications are implemented as stand-alone components of STING Millennium Suite.

AVAILABILITY

http://sms.cbi.cnptia.embrapa.br/SMS, http://trantor.bioc.columbia.edu/SMS, http://mirrors.rcsb.org/SMS, http://www.es.embnet.org/SMS and http://www.ar.embnet.org/SMS. [options: Scorpion, Formiga]

STING Millennium: A web-based suite of programs for comprehensive and simultaneous analysis of protein structure and sequence

STING Millennium Suite (SMS) is a new web-based suite of programs and databases providing visualization and a complex analysis of molecular sequence and structure for the data deposited at the Protein Data Bank (PDB). SMS operates with a collection of both publicly available data (PDB, HSSP, Prosite) and its own data (contacts, interface contacts, surface accessibility). Biologists find SMS useful because it provides a variety of algorithms and validated data, wrapped-up in a user friendly web interface. Using SMS it is now possible to analyze sequence to structure relationships, the quality of the structure, nature and volume of atomic contacts of intra and inter chain type, relative conservation of amino acids at the specific sequence position based on multiple sequence alignment, indications of folding essential residue (FER) based on the relationship of the residue conservation to the intra-chain contacts and Calpha-Calpha and Cbeta-Cbeta distance geometry. Specific emphasis in SMS is given to interface forming residues (IFR)-amino acids that define the interactive portion of the protein surfaces. SMS may simultaneously display and analyze previously superimposed structures. PDB updates trigger SMS updates in a synchronized fashion. SMS is freely accessible for public data at http://www.cbi.cnptia.embrapa.br, http://mirrors.rcsb.org/SMS and http://trantor.bioc.columbia.edu/SMS.

TMCompare: transmembrane region sequence and structure

TMCompare is an alignment and visualization tool for comparison of sequence information for membrane proteins contained in SWISS-PROT entries, with structural information contained in PDB files. The program can be used for: detection of breaks in alpha helical structure of transmembrane regions; examination of differences in coverage between PDB and SWISS-PROT files; examination of annotation differences between PDB files and associated SWISS-PROT files; examination and comparison of assigned PDB alpha helix regions and assigned SWISS-PROT transmembrane regions in linear sequence (one letter code) format; examination of these differences in 3D using the CHIME plugin, allowing; analysis of the alpha and non-alpha content of transmembrane regions.

AVAILABILITY

TMCompare is available for use through selection of a query protein via the internet (http://www.membraneproteins.org/TMCompare)

CONTACT

tmcompare@membraneproteins.org

Analysis of structural and physico-chemical parameters involved in the specificity of binding between alpha-amylases and their inhibitors

Enzyme-inhibitor specificity was studied for alpha-amylases and their inhibitors. We purified and cloned the cDNAs of two different alpha-amylase inhibitors from the common bean (Phaseolus vulgaris) and have recently cloned the cDNA of an alpha-amylase of the Mexican bean weevil (Zabrotes subfasciatus), which is inhibited by alpha-amylase inhibitor 2 but not by alpha-amylase inhibitor 1. The crystal structure of AI-1 complexed with pancreatic porcine alpha-amylase allowed us to model the structure of AI-2. The structure of Zabrotes subfasciatus alpha-amylase was modeled based on the crystal structure of Tenebrio molitor alpha-amylase. Pairwise AI-1 and AI-2 with PPA and ZSA complexes were modeled. For these complexes we first identified the interface forming residues. In addition, we identified the hydrogen bonds, ionic interactions and loss of hydrophobic surface area resulting from complex formation. The parameters we studied provide insight into the general scheme of binding, but fall short of explaining the specificity of the inhibition. We also introduce three new tools-software packages STING, HORNET and STINGPaint-which efficiently determine the interface forming residues and the ionic interaction data, the hydrogen bond net as well as aid in interpretation of multiple sequence alignment, respectively.