Analysis and visualization of RNA-Seq expression data using RStudio, Bioconductor, and Integrated Genome Browser.

Sequencing costs are falling, but the cost of data analysis remains high, often because unforeseen problems arise, such as insufficient depth of sequencing or batch effects. Experimenting with data analysis methods during the planning phase of an experiment can reveal unanticipated problems and build valuable bioinformatics expertise in the organism or process being studied. This protocol describes using R Markdown and RStudio, user-friendly tools for statistical analysis and reproducible research in bioinformatics, to analyze and document the analysis of an example RNA-Seq data set from tomato pollen undergoing chronic heat stress.

Annotating genes of known and unknown function by large-scale coexpression analysis.

About 40% of the proteins encoded in eukaryotic genomes are proteins of unknown function (PUFs). Their functional characterization remains one of the main challenges in modern biology. In this study we identified the PUF encoding genes from Arabidopsis (Arabidopsis thaliana) using a combination of sequence similarity, domain-based, and empirical approaches.

Charting plant interactomes: possibilities and challenges.

Protein-protein interactions are essential for nearly all cellular processes. Therefore, an important goal of post-genomic research for defining gene function and understanding the function of macromolecular complexes involves creating 'interactome' maps from empirical or inferred datasets. Systematic efforts to conduct high-throughput surveys of protein-protein interactions in plants are needed to chart the complex and dynamic interaction networks that occur throughout plant development.

POFs: what we don't know can hurt us.

Over a quarter of all eukaryotic genes encode proteins with obscure features that lack currently defined motifs or domains (POFs). Interestingly, most of the differences in gene repertoire among species were recently found to be attributed to POFs. A comparison of the Arabidopsis, rice and poplar genomes reveals that Arabidopsis contains 5069 POFs, of which 2045 have no obvious homologs in rice or poplar and are likely to be involved in species- or phylogenetic-specific functions in Arabidopsis.

What makes species unique? The contribution of proteins with obscure features.

Background: Proteins with obscure features (POFs), which lack currently defined motifs or domains, represent between 18% and 38% of a typical eukaryotic proteome. To evaluate the contribution of this class of proteins to the diversity of eukaryotes, we performed a comparative analysis of the predicted proteomes derived from 10 different sequenced genomes, including budding and fission yeast, worm, fly, mosquito, Arabidopsis, rice, mouse, rat, and human.

Results: Only 1,650 protein groups were found to be conserved among these proteomes (BLAST E-value threshold of 10(-6)).

A phyloproteomic characterization of in vitro autophosphorylation in calcium-dependent protein kinases.

Calcium-dependent protein kinases (CDPKs) are a novel class of signaling molecules that have been broadly implicated in relaying specific calcium-mediated responses to biotic and abiotic stress as well as developmental cues in both plants and protists. Calcium-dependent autophosphorylation has been observed in almost all CDPKs examined, but a physiological role for autophosphorylation has not been demonstrated. To date, only a handful of autophosphorylation sites have been mapped to specific residues within CDPK amino acid sequences.

Transcriptome changes for Arabidopsis in response to salt, osmotic, and cold stress.

To identify genes of potential importance to cold, salt, and drought tolerance, global expression profiling was performed on Arabidopsis plants subjected to stress treatments of 4 degrees C, 100 mM NaCl, or 200 mM mannitol, respectively. RNA samples were collected separately from leaves and roots after 3- and 27-h stress treatments. Profiling was conducted with a GeneChip microarray with probe sets for approximately 8,100 genes.