Antonio Muñoz-Mérida

Sma3s: A Three-Step Modular Annotator for Large Sequence Datasets

Automatic sequence annotation is an essential component of modern ‘omics’ studies, which aim to extract information from large collections of sequence data. Most existing tools use sequence homology to establish evolutionary relationships and assign putative functions to sequences. However, it can be difficult to define a similarity threshold that achieves sufficient coverage without sacrificing annotation quality. Defining the correct configuration is critical and can be challenging for non-specialist users. Thus, the development of robust automatic annotation techniques that generate high-quality annotations without needing expert knowledge would be very valuable for the research community. We present Sma3s, a tool for automatically annotating very large collections of biological sequences from any kind of gene library or genome. Sma3s is composed of three modules that progressively annotate query sequences using either: (i) very similar homologues, (ii) orthologous sequences or (iii) terms enriched in groups of homologous sequences. We trained the system using several random sets of known sequences, demonstrating average sensitivity and specificity values of ∼85%. In conclusion, Sma3s is a versatile tool for high-throughput annotation of a wide variety of sequence datasets that outperforms the accuracy of other well-established annotation algorithms, and it can enrich existing database annotations and uncover previously hidden features. Importantly, Sma3s has already been used in the functional annotation of two published transcriptomes.

Early and delayed long-term transcriptional changes and short-term transient responses during cold acclimation in olive leaves

Low temperature severely affects plant growth and development. To overcome this constraint, several plant species from regions having a cool season have evolved an adaptive response, called cold acclimation. We have studied this response in olive tree (Olea europaea L.) cv. Picual. Biochemical stress markers and cold-stress symptoms were detected after the first 24 h as sagging leaves. After 5 days, the plants were found to have completely recovered. Control and cold-stressed plants were sequenced by Illumina HiSeq 1000 paired-end technique. We also assembled a new olive transcriptome comprising 157,799 unigenes and found 6,309 unigenes differentially expressed in response to cold. Three types of response that led to cold acclimation were found: short-term transient response, early long-term response, and late long-term response. These subsets of unigenes were related to different biological processes. Early responses involved many cold-stress-responsive genes coding for, among many other things, C-repeat binding factor transcription factors, fatty acid desaturases, wax synthesis, and oligosaccharide metabolism. After long-term exposure to cold, a large proportion of gene down-regulation was found, including photosynthesis and plant growth genes. Up-regulated genes after long-term cold exposure were related to organelle fusion, nucleus organization, and DNA integration, including retrotransposons.

PreP+07: improvements of a user friendly tool to preprocess and analyse microarray data

BackgroundNowadays, microarray gene expression analysis is a widely used technology that scientists handle but whose final interpretation usually requires the participation of a specialist. The need for this participation is due to the requirement of some background in statistics that most users lack or have a very vague notion of. Moreover, programming skills could also be essential to analyse these data. An interactive, easy to use application seems therefore necessary to help researchers to extract full information from data and analyse them in a simple, powerful and confident way.ResultsPreP+07 is a standalone Windows XP application that presents a friendly interface for spot filtration, inter- and intra-slide normalization, duplicate resolution, dye-swapping, error removal and statistical analyses. Additionally, it contains two unique implementation of the procedures – double scan and Supervised Lowess-, a complete set of graphical representations – MA plot, RG plot, QQ plot, PP plot, PN plot – and can deal with many data formats, such as tabulated text, GenePix GPR and ArrayPRO. PreP+07 performance has been compared with the equivalent functions in Bioconductor using a tomato chip with 13056 spots. The number of differentially expressed genes considering p-values coming from the PreP+07 and Bioconductor Limma packages were statistically identical when the data set was only normalized; however, a slight variability was appreciated when the data was both normalized and scaled.ConclusionPreP+07 implementation provides a high degree of freedom in selecting and organizing a small set of widely used data processing protocols, and can handle many data formats. Its reliability has been proven so that a laboratory researcher can afford a statistical pre-processing of his/her microarray results and obtain a list of differentially expressed genes using PreP+07 without any programming skills. All of this gives support to scientists that have been using previous PreP releases since its first version in 2003.

Genetic and genome-wide transcriptomic analyses identify co-regulation of oxidative response and hormone transcript abundance with vitamin C content in tomato fruit

BackgroundL-ascorbic acid (AsA; vitamin C) is essential for all living plants where it functions as the main hydrosoluble antioxidant. It has diverse roles in the regulation of plant cell growth and expansion, photosynthesis, and hormone-regulated processes. AsA is also an essential component of the human diet, being tomato fruit one of the main sources of this vitamin. To identify genes responsible for AsA content in tomato fruit, transcriptomic studies followed by clustering analysis were applied to two groups of fruits with contrasting AsA content. These fruits were identified after AsA profiling of an F8 Recombinant Inbred Line (RIL) population generated from a cross between the domesticated species Solanum lycopersicum and the wild relative Solanum pimpinellifollium.ResultsWe found large variability in AsA content within the RIL population with individual RILs with up to 4-fold difference in AsA content. Transcriptomic analysis identified genes whose expression correlated either positively (PVC genes) or negatively (NVC genes) with the AsA content of the fruits. Cluster analysis using SOTA allowed the identification of subsets of co-regulated genes mainly involved in hormones signaling, such as ethylene, ABA, gibberellin and auxin, rather than any of the known AsA biosynthetic genes. Data mining of the corresponding PVC and NVC orthologs in Arabidopis databases identified flagellin and other ROS-producing processes as cues resulting in differential regulation of a high percentage of the genes from both groups of co-regulated genes; more specifically, 26.6% of the orthologous PVC genes, and 15.5% of the orthologous NVC genes were induced and repressed, respectively, under flagellin22 treatment in Arabidopsis thaliana.ConclusionResults here reported indicate that the content of AsA in red tomato fruit from our selected RILs are not correlated with the expression of genes involved in its biosynthesis. On the contrary, the data presented here supports that AsA content in tomato fruit co-regulates with genes involved in hormone signaling and they are dependent on the oxidative status of the fruit.

Semirna: Searching for Plant miRNAs Using Target Sequences

Many plant genomes are already known, and new ones are being sequenced every year. The next step for researchers is to identify all of the functional elements in these genomes, including the important class of functional elements known as microRNAs (miRNAs), which are involved in posttranscriptional regulatory pathways. However, computational tools for predicting new plant miRNAs are limited, and there is a particular need for tools that can be used easily by laboratory researchers. We present semirna, a new tool for predicting miRNAs in plant genomes, available as a Web server. This tool takes a putative target sequence such as a messenger RNA (mRNA) as input, and allows users to search for miRNAs that target this sequence. It can also be used to determine whether small RNA sequences from massive sequencing analysis represent true miRNAs and to search for miRNAs in new genomes using homology. Semirna has shown a high level of accuracy using various test sets, and gives users the ability to search for miRNAs with several different adjustable parameters. Semirna, a user-friendly and intuitive Web server for predicting miRNA sequences, can be reached at http://www.bioinfocabd.upo.es/semirna/ . It is useful for researchers searching for miRNAs involved in particular pathways, as well as those searching for miRNAs in newly sequenced genomes.

Transcriptomic Analysis Using Olive Varieties and Breeding Progenies Identifies Candidate Genes Involved in Plant Architecture

Plant architecture is a critical trait in fruit crops that can significantly influence yield, pruning, planting density and harvesting. Little is known about how plant architecture is genetically determined in olive, were most of the existing varieties are traditional with an architecture poorly suited for modern growing and harvesting systems. In the present study, we have carried out microarray analysis of meristematic tissue to compare expression profiles of olive varieties displaying differences in architecture, as well as seedlings from their cross pooled on the basis of their sharing architecture-related phenotypes. The microarray used, previously developed by our group has already been applied to identify candidates genes involved in regulating juvenile to adult transition in the shoot apex of seedlings. Varieties with distinct architecture phenotypes and individuals from segregating progenies displaying opposite architecture features were used to link phenotype to expression. Here, we identify 2252 differentially expressed genes (DEGs) associated to differences in plant architecture. Microarray results were validated by quantitative RT-PCR carried out on genes with functional annotation likely related to plant architecture. Twelve of these genes were further analyzed in individual seedlings of the corresponding pool. We also examined Arabidopsis mutants in putative orthologs of these targeted candidate genes, finding altered architecture for most of them. This supports a functional conservation between species and potential biological relevance of the candidate genes identified. This study is the first to identify genes associated to plant architecture in olive, and the results obtained could be of great help in future programs aimed at selecting phenotypes adapted to modern cultivation practices in this species.

MOWServ: a web client for integration of bioinformatic resources.

The productivity of any scientist is affected by cumbersome, tedious and time-consuming tasks that try to make the heterogeneous web services compatible so that they can be useful in their research. MOWServ, the bioinformatic platform offered by the Spanish National Institute of Bioinformatics, was released to provide integrated access to databases and analytical tools. Since its release, the number of available services has grown dramatically, and it has become one of the main contributors of registered services in the EMBRACE Biocatalogue. The ontology that enables most of the web-service compatibility has been curated, improved and extended. The service discovery has been greatly enhanced by Magallanes software and biodataSF. User data are securely stored on the main server by an authentication protocol that enables the monitoring of current or already-finished user’s tasks, as well as the pipelining of successive data processing services. The BioMoby standard has been greatly extended with the new features included in the MOWServ, such as management of additional information (metadata such as extended descriptions, keywords and datafile examples), a qualified registry, error handling, asynchronous services and service replication. All of them have increased the MOWServ service quality, usability and robustness. MOWServ is available at http://www.inab.org/MOWServ/ and has a mirror at http://www.bitlab-es.com/MOWServ/.

Partial Activation of SA- and JA-Defensive Pathways in Strawberry upon Colletotrichum acutatum Interaction

Understanding the nature of pathogen host interaction may help improve strawberry (Fragaria × ananassa) cultivars. Plant resistance to pathogenic agents usually operates through a complex network of defense mechanisms mediated by a diverse array of signaling molecules. In strawberry, resistance to a variety of pathogens has been reported to be mostly polygenic and quantitatively inherited, making it difficult to associate molecular markers with disease resistance genes. Colletotrichum acutatum spp. is a major strawberry pathogen, and completely resistant cultivars have not been reported. Moreover, strawberry defense network components and mechanisms remain largely unknown and poorly understood. Assessment of the strawberry response to C. acutatum included a global transcript analysis, and acidic hormones SA and JA measurements were analyzed after challenge with the pathogen. Induction of transcripts corresponding to the SA and JA signaling pathways and key genes controlling major steps within these defense pathways was detected. Accordingly, SA and JA accumulated in strawberry after infection. Contrastingly, induction of several important SA, JA, and oxidative stress-responsive defense genes, including FaPR1-1, FaLOX2, FaJAR1, FaPDF1, and FaGST1, was not detected, which suggests that specific branches in these defense pathways (those leading to FaPR1-2, FaPR2-1, FaPR2-2, FaAOS, FaPR5, and FaPR10) were activated. Our results reveal that specific aspects in SA and JA dependent signaling pathways are activated in strawberry upon interaction with C. acutatum. Certain described defense-associated transcripts related to these two known signaling pathways do not increase in abundance following infection. This finding suggests new insight into a specific putative molecular strategy for defense against this pathogen.

Usefulness of a New Large Set of High Throughput EST-SNP Markers as a Tool for Olive Germplasm Collection Management

Germplasm collections are basic tools for conservation, characterization, and efficient use of olive genetic resources. The identification of the olive cultivars maintained in the collections is an important ongoing task which has been performed by both, morphological and molecular markers. In the present study, based on the sequencing results of previous genomic projects, a new set of 1,043 EST-SNP markers has been identified. In order to evaluate its discrimination capacity and utility in diversity studies, this set of markers was used in a representative number of accessions from 20 different olive growing countries and maintained at the World Olive Germplasm Collection of IFAPA Centre ‘Alameda del Obispo’ (Córdoba, Spain), one of the world’s largest olive germplasm bank. Thus, the cultivated material included: cultivars belonging to previously defined core collections by means of SSR markers and agronomical traits, well known homonymy cases, possible redundancies previously identified in the collection, and recently introduced accessions. Marker stability was tested in repeated analyses of a selected number of accessions, as well as in different trees and accessions belonging to the same cultivar. In addition, 15 genotypes from a cross ‘Picual’ × ‘Arbequina’ cultivars from the IFAPA olive breeding program and a set of 89 wild genotypes were also included in the study. Our results indicate that, despite their relatively wide variability, the new set of EST-SNPs displayed lower levels of genetic diversity than SSRs in the set of olive core collections tested. However, the EST-SNP markers displayed consistent and reliable results from different plant material sources and plant propagation events. The EST-SNPs revealed a clear cut off between inter- and intra-cultivar variation in olive. Besides, they were able to reliably discriminate among different accessions, to detect possible homonymy cases as well as efficiently ascertain the presence of redundant germplasm in the collection. Additionally, these markers were highly transferable to the wild genotypes. These results, together with the low genotyping error rates and the easy and fully automated procedure used to get the genotyping data, validate the new set of EST-SNPs as possible markers of choice for olive cultivar identification.

Statistical significance for NGS reads similarities

In this work we present a significance curve to segregate random alignments from true matches in by identity sequence comparison, especially suitable for sequencing data produced by NGS-technologies. The experimental approach reproduces the random local ungapped similarities distribution by score and length from which it is possible to asses the statistical significance of any particular ungapped similarity. This work includes the study of the distribution behaviour as a function of the experimental technology used to produce the raw sequences, as well as the scoring system used in the comparison. Our approach reproduces the expected behaviour and completes the proposal of Rost and Sander for homology based sequence comparisons. Results can be exploited by computational applications to reduce the computational cost and memory usage.