Using The Arabidopsis Information Resource (TAIR) to Find Information About Arabidopsis Genes

Developing a bioinformatics pipeline for comparative protein classification analysis

BACKGROUND

Protein classification is a task of paramount importance in various fields of biology. Despite the great momentum of modern implementation of protein classification, machine learning techniques such as Random Forest and Neural Network could not always be used for several reasons: data collection, unbalanced classification or labelling of the data.As an alternative, I propose the use of a bioinformatics pipeline to search for and classify information from protein databases. Hence, to evaluate the efficiency and accuracy of the pipeline, I focused on the carotenoid biosynthetic genes and developed a filtering approach to retrieve orthologs clusters in two well-studied plants that belong to the Brassicaceae family: Arabidopsis thaliana and Brassica rapa Pekinensis group. The result obtained has been compared with previous studies on carotenoid biosynthetic genes in B. rapa where phylogenetic analysis was conducted.

RESULTS

The developed bioinformatics pipeline relies on commercial software and multiple databeses including the use of phylogeny, Gene Ontology terms (GOs) and Protein Families (Pfams) at a protein level. Furthermore, the phylogeny is coupled with "population analysis" to evaluate the potential orthologs. All the steps taken together give a final table of potential orthologs. The phylogenetic tree gives a result of 43 putative orthologs conserved in B. rapa Pekinensis group. Different A. thaliana proteins have more than one syntenic ortholog as also shown in a previous finding (Li et al., BMC Genomics 16(1):1-11, 2015).

CONCLUSIONS

This study demonstrates that, when the biological features of proteins of interest are not specific, I can rely on a computational approach in filtering steps for classification purposes. The comparison of the results obtained here for the carotenoid biosynthetic genes with previous research confirmed the accuracy of the developed pipeline which can therefore be applied for filtering different types of datasets.

Plant-necrotroph co-transcriptome networks illuminate a metabolic battlefield

A central goal of studying host-pathogen interaction is to understand how host and pathogen manipulate each other to promote their own fitness in a pathosystem. Co-transcriptomic approaches can simultaneously analyze dual transcriptomes during infection and provide a systematic map of the cross-kingdom communication between two species. Here we used the Arabidopsis-B. cinerea pathosystem to test how plant host and fungal pathogen interact at the transcriptomic level. We assessed the impact of genetic diversity in pathogen and host by utilization of a collection of 96 isolates infection on Arabidopsis wild-type and two mutants with jasmonate or salicylic acid compromised immunities. We identified ten B. cinereagene co-expression networks (GCNs) that encode known or novel virulence mechanisms. Construction of a dual interaction network by combining four host- and ten pathogen-GCNs revealed potential connections between the fungal and plant GCNs. These co-transcriptome data shed lights on the potential mechanisms underlying host-pathogen interaction.

Arabidopsis thaliana SPF1 and SPF2 are nuclear-located ULP2-like SUMO proteases that act downstream of SIZ1 in plant development

Post-translational modifiers such as the small ubiquitin-like modifier (SUMO) peptide act as fast and reversible protein regulators. Functional characterization of the sumoylation machinery has determined the key regulatory role that SUMO plays in plant development. Unlike components of the SUMO conjugation pathway, SUMO proteases (ULPs) are encoded by a relatively large gene family and are potential sources of specificity within the pathway. This study reports a thorough comparative genomics and phylogenetic characterization of plant ULPs, revealing the presence of one ULP1-like and three ULP2-like SUMO protease subgroups within plant genomes. As representatives of an under-studied subgroup, Arabidopsis SPF1 and SPF2 were subjected to functional characterization. Loss-of-function mutants implicated both proteins with vegetative growth, flowering time, and seed size and yield. Mutants constitutively accumulated SUMO conjugates, and yeast complementation assays associated these proteins with the function of ScUlp2 but not ScUlp1. Fluorescence imaging placed both proteins in the plant cell nucleoplasm. Transcriptomics analysis indicated strong regulatory involvement in secondary metabolism, cell wall remodelling, and nitrate assimilation. Furthermore, developmental defects of the spf1-1 spf2-2 (spf1/2) double-mutant opposed those of the major E3 ligase siz1 mutant and, most significantly, developmental and transcriptomic characterization of the siz1 spf1/2 triple-mutant placed SIZ1 as epistatic to SPF1 and SPF2.

Sustainable funding for biocuration: The Arabidopsis Information Resource (TAIR) as a case study of a subscription-based funding model

Databases and data repositories provide essential functions for the research community by integrating, curating, archiving and otherwise packaging data to facilitate discovery and reuse. Despite their importance, funding for maintenance of these resources is increasingly hard to obtain. Fueled by a desire to find long term, sustainable solutions to database funding, staff from the Arabidopsis Information Resource (TAIR), founded the nonprofit organization, Phoenix Bioinformatics, using TAIR as a test case for user-based funding. Subscription-based funding has been proposed as an alternative to grant funding but its application has been very limited within the nonprofit sector. Our testing of this model indicates that it is a viable option, at least for some databases, and that it is possible to strike a balance that maximizes access while still incentivizing subscriptions. One year after transitioning to subscription support, TAIR is self-sustaining and Phoenix is poised to expand and support additional resources that wish to incorporate user-based funding strategies. Database URL: www.arabidopsis.org.

An Introduction to Genome Annotation

Genome projects have evolved from large international undertakings to tractable endeavors for a single lab. Accurate genome annotation is critical for successful genomic, genetic, and molecular biology experiments. These annotations can be generated using a number of approaches and available software tools. This unit describes methods for genome annotation and a number of software tools commonly used in gene annotation.

The Arabidopsis information resource: Making and mining the "gold standard" annotated reference plant genome

The Arabidopsis Information Resource (TAIR) is a continuously updated, online database of genetic and molecular biology data for the model plant Arabidopsis thaliana that provides a global research community with centralized access to data for over 30,000 Arabidopsis genes. TAIR's biocurators systematically extract, organize, and interconnect experimental data from the literature along with computational predictions, community submissions, and high throughput datasets to present a high quality and comprehensive picture of Arabidopsis gene function. TAIR provides tools for data visualization and analysis, and enables ordering of seed and DNA stocks, protein chips, and other experimental resources. TAIR actively engages with its users who contribute expertise and data that augments the work of the curatorial staff. TAIR's focus in an extensive and evolving ecosystem of online resources for plant biology is on the critically important role of extracting experimentally based research findings from the literature and making that information computationally accessible. In response to the loss of government grant funding, the TAIR team founded a nonprofit entity, Phoenix Bioinformatics, with the aim of developing sustainable funding models for biological databases, using TAIR as a test case. Phoenix has successfully transitioned TAIR to subscription-based funding while still keeping its data relatively open and accessible.

The Importance of Biological Databases in Biological Discovery

Biological databases play a central role in bioinformatics. They offer scientists the opportunity to access a wide variety of biologically relevant data, including the genomic sequences of an increasingly broad range of organisms. This unit provides a brief overview of major sequence databases and portals, such as GenBank, the UCSC Genome Browser, and Ensembl. Model organism databases, including WormBase, The Arabidopsis Information Resource (TAIR), and those made available through the Mouse Genome Informatics (MGI) resource, are also covered. Non-sequence-centric databases, such as Online Mendelian Inheritance in Man (OMIM), the Protein Data Bank (PDB), MetaCyc, and the Kyoto Encyclopedia of Genes and Genomes (KEGG), are also discussed.

Arabidopsis database and stock resources

The volume of Arabidopsis information has increased enormously in recent years as a result of the sequencing of the reference genome and other large-scale functional genomics projects. Much of the data is stored in public databases, where data are organized, analyzed, and made freely accessible to the research community. These databases are resources that researchers can utilize for making predictions and developing testable hypotheses. The methods in this chapter describe ways to access and utilize Arabidopsis data and genomic resources found in databases and stock centers.

Integration of latex protein sequence data provides comprehensive functional overview of latex proteins

The laticiferous system is one of the most important conduit systems in higher plants, which produces a milky-like sap known as latex. Latex contains diverse secondary metabolites with various ecological functions. To obtain a comprehensive overview of the latex proteome, we integrated available latex proteins sequences and constructed a comprehensive dataset composed of 1,208 non-redundant latex proteins from 20 various latex-bearing plants. The results of functional analyses revealed that latex proteins are involved in various biological processes, including transcription, translation, protein degradation and the plant response to environmental stimuli. The results of the comparative analysis showed that the functions of the latex proteins are similar to those of phloem, suggesting the functional conservation of plant vascular proteins. The presence of latex proteins in mitochondria and plastids suggests the production of diverse secondary metabolites. Furthermore, using a BLAST search, we identified 854 homologous latex proteins in eight plant species, including three latex-bearing plants, such as papaya, caster bean and cassava, suggesting that latex proteins were newly evolved in vascular plants. Taken together, this study is the largest and most comprehensive in silico analysis of the latex proteome. The results obtained here provide useful resources and information for characterizing the evolution of the latex proteome.

Identification of stable, high copy number, medium-sized RNA degradation intermediates that accumulate in plants under non-stress conditions

It is becoming increasingly evident that the RNA degradome is a crucial component of the total cellular RNA pool. Here, we present an analysis of the medium-sized RNAs (midi RNAs) that form in Arabidopsis thaliana. Our analyses revealed that the midi RNA fraction contained mostly 20-70-nt-long fragments derived from various RNA species, including tRNA, rRNA, mRNA and snRNA. The majority of these fragments could be classified as stable RNA degradation intermediates (RNA degradants). Using two dimensional polyacrylamide gel electrophoresis, we demonstrated that high copy number RNA (hcn RNA) degradants appear in plant cells not only during stress, as it was earlier suggested. They are continuously produced also under physiological conditions. The data collected indicated that the accumulation pattern of the hcn RNA degradants is organ-specific and can be affected by various endogenous and exogenous factors. In addition, we demonstrated that selected degradants efficiently inhibit translation in vitro. Thus, the results of our studies suggest that hcn RNA degradants are likely to be involved in the regulation of gene expression in plants.

Small RNA and transcriptome deep sequencing proffers insight into floral gene regulation in Rosa cultivars

Background

Roses (Rosa sp.), which belong to the family Rosaceae, are the most economically important ornamental plants—making up 30% of the floriculture market. However, given high demand for roses, rose breeding programs are limited in molecular resources which can greatly enhance and speed breeding efforts. A better understanding of important genes that contribute to important floral development and desired phenotypes will lead to improved rose cultivars. For this study, we analyzed rose miRNAs and the rose flower transcriptome in order to generate a database to expound upon current knowledge regarding regulation of important floral characteristics. A rose genetic database will enable comprehensive analysis of gene expression and regulation via miRNA among different Rosa cultivars.

Results

We produced more than 0.5 million reads from expressed sequences, totalling more than 110 million bp. From these, we generated 35,657, 31,434, 34,725, and 39,722 flower unigenes from Rosa hybrid: ‘Vital’, ‘Maroussia’, and ‘Sympathy’ and Rosa rugosa Thunb. , respectively. The unigenes were assigned functional annotations, domains, metabolic pathways, Gene Ontology (GO) terms, Plant Ontology (PO) terms, and MIPS Functional Catalogue (FunCat) terms. Rose flower transcripts were compared with genes from whole genome sequences of Rosaceae members (apple, strawberry, and peach) and grape. We also produced approximately 40 million small RNA reads from flower tissue for Rosa, representing 267 unique miRNA tags. Among identified miRNAs, 25 of them were novel and 242 of them were conserved miRNAs. Statistical analyses of miRNA profiles revealed both shared and species-specific miRNAs, which presumably effect flower development and phenotypes.

Conclusions

In this study, we constructed a Rose miRNA and transcriptome database, and we analyzed the miRNAs and transcriptome generated from the flower tissues of four Rosa cultivars. The database provides a comprehensive genetic resource which can be used to better understand rose flower development and to identify candidate genes for important phenotypes.

TE-Locate: A Tool to Locate and Group Transposable Element Occurrences Using Paired-End Next-Generation Sequencing Data

Transposable elements (TEs) are common mobile DNA elements present in nearly all genomes. Since the movement of TEs within a genome can sometimes have phenotypic consequences, an accurate report of TE actions is desirable. To this end, we developed TE-Locate, a computational tool that uses paired-end reads to identify the novel locations of known TEs. TE-Locate can utilize either a database of TE sequences, or annotated TEs within the reference sequence of interest. This makes TE-Locate useful in the search for any mobile sequence, including retrotransposed gene copies. One major concern is to act on the correct hierarchy level, thereby avoiding an incorrect calling of a single insertion as multiple events of TEs with high sequence similarity. We used the (super)family level, but TE-Locate can also use any other level, right down to the individual transposable element. As an example of analysis with TE-Locate, we used the Swedish population in the 1,001 Arabidopsis genomes project, and presented the biological insights gained from the novel TEs, inducing the association between different TE superfamilies. The program is freely available, and the URL is provided in the end of the paper.

In vivo visualization of RNA in plants cells using the λN₂₂ system and a GATEWAY-compatible vector series for candidate RNAs

The past decade has seen a tremendous increase in RNA research, which has demonstrated that RNAs are involved in many more processes than were previously thought. The dynamics of RNA synthesis towards their regulated activity requires the interplay of RNAs with numerous RNA binding proteins (RBPs). The localization of RNA, a mechanism for controlling translation in a spatial and temporal fashion, requires processing and assembly of RNA into transport granules in the nucleus, transport towards cytoplasmic destinations and regulation of its activity. Compared with animal model systems little is known about RNA dynamics and motility in plants. Commonly used methods to study RNA transport and localization are time-consuming, and require expensive equipment and a high level of experimental skill. Here, we introduce the λN₂₂ RNA stem-loop binding system for the in vivo visualization of RNA in plant cells. The λN₂₂ system consists of two components: the λN₂₂ RNA binding peptide and the corresponding box-B stem loops. We generated fusions of λN₂₂ to different fluorophores and a GATEWAY vector series for the simple fusion of any target RNA 5' or 3' to box-B stem loops. We show that the λN₂₂ system can be used to detect RNAs in transient expression assays, and that it offers advantages compared with the previously described MS2 system. Furthermore, the λN₂₂ system can be used in combination with the MS2 system to visualize different RNAs simultaneously in the same cell. The toolbox of vectors generated for both systems is easy to use and promises significant progress in our understanding of RNA transport and localization in plant cells.

Mining Arabidopsis thaliana RNA-seq data with Integrated Genome Browser reveals stress-induced alternative splicing of the putative splicing regulator SR45a

PREMISE OF THE STUDY

High-throughput sequencing of cDNA libraries prepared from diverse samples (RNA-seq) can reveal genome-wide changes in alternative splicing. Using RNA-seq data to assess splicing at the level of individual genes requires the ability to visualize read alignments alongside genomic annotations. To meet this need, we added RNA-seq visualization capability to Integrated Genome Browser (IGB), a free desktop genome visualization tool. To illustrate this capability, we present an in-depth analysis of abiotic stresses and their effects on alternative splicing of SR45a (AT1G07350), a putative splicing regulator from Arabidopsis thaliana.

METHODS

cDNA libraries prepared from Arabidopsis plants that were subjected to heat and dehydration stresses were sequenced on an Illumina GAIIx sequencer, yielding more than 511 million high-quality 75-base, single-end sequence reads. Reads were aligned onto the reference genome and visualized in IGB.

KEY RESULTS

Using IGB, we confirmed exon-skipping alternative splicing in SR45a. Exon-skipped variant AT1G07350.1 encodes full-length SR45a protein with intact RS and RNA recognition motifs, while nonskipped variant AT1G07350.2 lacks the C-terminal RS region due to a frameshift in the alternative exon. Heat and drought stresses increased both transcript abundance and the proportion of exon-skipped transcripts encoding the full-length protein. We identified new splice sites and observed frequent intron retention flanking the alternative exon.

CONCLUSIONS

This study underlines the importance of visual inspection of RNA-seq alignments when investigating alternatively spliced genes. We showed that heat and dehydration stresses increase overall abundance of SR45a mRNA while also increasing production of transcripts encoding the full-length SR45a protein relative to other splice variants.

mirEX: a platform for comparative exploration of plant pri-miRNA expression data

mirEX is a comprehensive platform for comparative analysis of primary microRNA expression data. RT–qPCR-based gene expression profiles are stored in a universal and expandable database scheme and wrapped by an intuitive user-friendly interface. A new way of accessing gene expression data in mirEX includes a simple mouse operated querying system and dynamic graphs for data mining analyses. In contrast to other publicly available databases, the mirEX interface allows a simultaneous comparison of expression levels between various microRNA genes in diverse organs and developmental stages. Currently, mirEX integrates information about the expression profile of 190 Arabidopsis thaliana pri-miRNAs in seven different developmental stages: seeds, seedlings and various organs of mature plants. Additionally, by providing RNA structural models, publicly available deep sequencing results, experimental procedure details and careful selection of auxiliary data in the form of web links, mirEX can function as a one-stop solution for Arabidopsis microRNA information. A web-based mirEX interface can be accessed at http://bioinfo.amu.edu.pl/mirex.

High-density kinetic analysis of the metabolomic and transcriptomic response of Arabidopsis to eight environmental conditions

The time-resolved response of Arabidopsis thaliana towards changing light and/or temperature at the transcriptome and metabolome level is presented. Plants grown at 21°C with a light intensity of 150 μE m⁻² sec⁻¹ were either kept at this condition or transferred into seven different environments (4°C, darkness; 21°C, darkness; 32°C, darkness; 4°C, 85 μE m⁻² sec⁻¹; 21 °C, 75 μE m⁻² sec⁻¹; 21°C, 300 μE m⁻² sec⁻¹ ; 32°C, 150 μE m⁻² sec⁻¹). Samples were taken before (0 min) and at 22 time points after transfer resulting in (8×) 22 time points covering both a linear and a logarithmic time series totaling 177 states. Hierarchical cluster analysis shows that individual conditions (defined by temperature and light) diverge into distinct trajectories at condition-dependent times and that the metabolome follows different kinetics from the transcriptome. The metabolic responses are initially relatively faster when compared with the transcriptional responses. Gene Ontology over-representation analysis identifies a common response for all changed conditions at the transcriptome level during the early response phase (5-60 min). Metabolic networks reconstructed via metabolite-metabolite correlations reveal extensive environment-specific rewiring. Detailed analysis identifies conditional connections between amino acids and intermediates of the tricarboxylic acid cycle. Parallel analysis of transcriptional changes strongly support a model where in the absence of photosynthesis at normal/high temperatures protein degradation occurs rapidly and subsequent amino acid catabolism serves as the main cellular energy supply. These results thus demonstrate the engagement of the electron transfer flavoprotein system under short-term environmental perturbations.

Next generation functional proteomics in non-model plants: A survey on techniques and applications for the analysis of protein complexes and post-translational modifications

The congruent development of computational technology, bioinformatics and analytical instrumentation makes proteomics ready for the next leap. Present-day state of the art proteomics grew from a descriptive method towards a full stake holder in systems biology. High throughput and genome wide studies are now made at the functional level. These include quantitative aspects, functional aspects with respect to protein interactions as well as post translational modifications and advanced computational methods that aid in predicting protein function and mapping these functionalities across the species border. In this review an overview is given of the current status of these aspects in plant studies with special attention to non-genomic model plants.

Arabidopsis plasmodesmal proteome

The multicellular nature of plants requires that cells should communicate in order to coordinate essential functions. This is achieved in part by molecular flux through pores in the cell wall, called plasmodesmata. We describe the proteomic analysis of plasmodesmata purified from the walls of Arabidopsis suspension cells. Isolated plasmodesmata were seen as membrane-rich structures largely devoid of immunoreactive markers for the plasma membrane, endoplasmic reticulum and cytoplasmic components. Using nano-liquid chromatography and an Orbitrap ion-trap tandem mass spectrometer, 1341 proteins were identified. We refer to this list as the plasmodesmata- or PD-proteome. Relative to other cell wall proteomes, the PD-proteome is depleted in wall proteins and enriched for membrane proteins, but still has a significant number (35%) of putative cytoplasmic contaminants, probably reflecting the sensitivity of the proteomic detection system. To validate the PD-proteome we searched for known plasmodesmal proteins and used molecular and cell biological techniques to identify novel putative plasmodesmal proteins from a small subset of candidates. The PD-proteome contained known plasmodesmal proteins and some inferred plasmodesmal proteins, based upon sequence or functional homology with examples identified in different plant systems. Many of these had a membrane association reflecting the membranous nature of isolated structures. Exploiting this connection we analysed a sample of the abundant receptor-like class of membrane proteins and a small random selection of other membrane proteins for their ability to target plasmodesmata as fluorescently-tagged fusion proteins. From 15 candidates we identified three receptor-like kinases, a tetraspanin and a protein of unknown function as novel potential plasmodesmal proteins. Together with published work, these data suggest that the membranous elements in plasmodesmata may be rich in receptor-like functions, and they validate the content of the PD-proteome as a valuable resource for the further uncovering of the structure and function of plasmodesmata as key components in cell-to-cell communication in plants.

Identification of drought-responsive universal stress proteins in viridiplantae

Genes encoding proteins that contain the universal stress protein (USP) domain are known to provide bacteria, archaea, fungi, protozoa, and plants with the ability to respond to a plethora of environmental stresses. Specifically in plants, drought tolerance is a desirable phenotype. However, limited focused and organized functional genomic datasets exist on drought-responsive plant USP genes to facilitate their characterization. The overall objective of the investigation was to identify diverse plant universal stress proteins and Expressed Sequence Tags (ESTs) responsive to water-deficit stress. We hypothesize that cross-database mining of functional annotations in protein and gene transcript bioinformatics resources would help identify candidate drought-responsive universal stress proteins and transcripts from multiple plant species. Our bioinformatics approach retrieved, mined and integrated comprehensive functional annotation data on 511 protein and 1561 ESTs sequences from 161 viridiplantae taxa. A total of 32 drought-responsive ESTs from 7 plant genera Glycine, Hordeum, Manihot, Medicago, Oryza, Pinus and Triticum were identified. Two Arabidopsis USP genes At3g62550 and At3g53990 that encode ATP-binding motif were up-regulated in a drought microarray dataset. Further, a dataset of 80 simple sequence repeats (SSRs) linked to 20 singletons and 47 transcript assembles was constructed. Integrating the datasets on SSRs and drought-responsive ESTs identified three drought-responsive ESTs from bread wheat (BE604157), soybean (BM887317) and maritime pine (BX682209). The SSR sequence types were CAG, ATA and AT respectively. The datasets from cross-database mining provide organized resources for the characterization of USP genes as useful targets for engineering plant varieties tolerant to unfavorable environmental conditions.