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

MAIN CONCLUSION: 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.

The genome of Eucalyptus grandis.

Eucalypts are the world's most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

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.

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.

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.

Quantitative prey species detection in predator guts across multiple trophic levels by mapping unassembled shotgun reads.

Quantifying species trophic interaction strengths is crucial for understanding community dynamics and has significant implications for pest management and species conservation. DNA-based methods to identify species interactions have revolutionized these efforts, but a significant limitation is the poor ability to quantify the strength of trophic interactions, that is the biomass or number of prey consumed. We present an improved pipeline, called Lazaro, to map unassembled shotgun reads to a comprehensive arthropod mitogenome database and show that the number of prey reads detected is quantitatively predicted from the prey biomass consumed, even for indirect predation. Two feeding bioassays were performed: starved coccinellid larvae consuming different numbers of aphids (Prey Quantity bioassay), and starved coccinellid larvae consuming a chrysopid larvae that had consumed aphids (Direct and Indirect Predation bioassay). Prey taxonomic assignment against a mitochondrial genome database had high accuracy (99.8% positive predictive value) and the number of prey reads was directly related to the number of prey consumed and inversely related to the elapsed time since consumption with high significance (r2  = .932, p = 4.92E-6). Aphids were detected up to 6 h after direct predation plus 3 h after indirect predation (9 h in total) and detection was related to the predator-specific decay rates. Lazaro enabled quantitative predictions of prey consumption across multiple trophic levels with high taxonomic resolution while eliminating all false positives, except for a few confirmed contaminants, and may be valuable for characterizing prey consumed by field-sampled predators. Moreover, Lazaro is readily applicable for species diversity determination from any degraded environmental DNA.

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.

Low levels of sulfur and cobalt during the pre- and periconceptional periods affect the oocyte yield of donors and the DNA methylome of preimplantation bovine embryos.

Maternal nutrition is critical in mammalian development, influencing the epigenetic reprogramming of gametes, embryos, and fetal programming. We evaluated the effects of different levels of sulfur (S) and cobalt (Co) in the maternal diet throughout the pre- and periconceptional periods on the biochemical and reproductive parameters of the donors and the DNA methylome of the progeny in Bos indicus cattle. The low-S/Co group differed from the control with respect to homocysteine, folic acid, B12, insulin growth factor 1, and glucose. The oocyte yield was lower in heifers from the low S/Co group than that in the control heifers. Embryos from the low-S/Co group exhibited 2320 differentially methylated regions (DMRs) across the genome compared with the control embryos. We also characterized candidate DMRs linked to the DNMT1 and DNMT3B genes in the blood and sperm cells of the adult progeny. A DMR located in DNMT1 that was identified in embryos remained differentially methylated in the sperm of the progeny from the low-S/Co group. Therefore, we associated changes in specific compounds in the maternal diet with DNA methylation modifications in the progeny. Our results help to elucidate the impact of maternal nutrition on epigenetic reprogramming in livestock, opening new avenues of research to study the effect of disturbed epigenetic patterns in early life on health and fertility in adulthood. Considering that cattle are physiologically similar to humans with respect to gestational length, our study may serve as a model for studies related to the developmental origin of health and disease in humans.

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 T2 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.

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.

Complete Genome Sequences of Seven New Chrysodeixis includens Nucleopolyhedrovirus Isolates from Minas Gerais and Mato Grosso States in Brazil.

We report the complete genomic sequences of seven viral isolates from the soybean looper (Chrysodeixis includens) from midwestern and southeastern Brazil. The genomes range from 138,760 to 139,637 bp in length with a G+C content of 39.2% and 140 open reading frames (ORFs).

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 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.

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.

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.

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 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.

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 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.