Evolutionary Transcriptomics and Proteomics: Insight into Plant Adaptation.

Comparative transcriptomics and proteomics (T&P) have brought biological insight into development, gene function, and physiological stress responses. However, RNA-seq and high-throughput proteomics remain underutilised in studies of plant adaptation. These methodologies have created discovery tools with the potential to significantly advance our understanding of adaptive diversification. We outline experimental recommendations for their application. We discuss analysis models and approaches that accelerate the identification of adaptive gene sets and integrate transcriptome, proteome, phenotypic, and environmental data. Finally, we encourage widespread uptake and future developments in T&P that will advance our understanding of evolution and adaptation.

De novo sequence assembly and characterisation of a partial transcriptome for an evolutionarily distinct reptile, the tuatara (Sphenodon punctatus).

BACKGROUND: The tuatara (Sphenodon punctatus) is a species of extraordinary zoological interest, being the only surviving member of an entire order of reptiles which diverged early in amniote evolution. In addition to their unique phylogenetic placement, many aspects of tuatara biology, including temperature-dependent sex determination, cold adaptation and extreme longevity have the potential to inform studies of genome evolution and development. Despite increasing interest in the tuatara genome, genomic resources for the species are still very limited. We aimed to address this by assembling a transcriptome for tuatara from an early-stage embryo, which will provide a resource for genome annotation, molecular marker development and studies of development and adaptation in tuatara. RESULTS: We obtained 30 million paired-end 50 bp reads from an Illumina Genome Analyzer and assembled them with Velvet and Oases using a range of kmers. After removing redundancy and filtering out low quality transcripts, our transcriptome dataset contained 32911 transcripts, with an N50 of 675 and a mean length of 451 bp. Almost 50% (15965) of these transcripts could be annotated by comparison with the NCBI non-redundant (NR) protein database or the chicken, green anole and zebrafish UniGene sets. A scan of candidate genes and repetitive elements revealed genes involved in immune function, sex differentiation and temperature-sensitivity, as well as over 200 microsatellite markers. CONCLUSIONS: This dataset represents a major increase in genomic resources for the tuatara, increasing the number of annotated gene sequences from just 60 to almost 16,000. This will facilitate future research in sex determination, genome evolution, local adaptation and population genetics of tuatara, as well as inform studies on amniote evolution.

Chips and tags suggest plant-environment interactions differ for two alpine Pachycladon species.

BACKGROUND: Expression profiling has been proposed as a means for screening non-model organisms in their natural environments to identify genes potentially important in adaptive diversification. Tag profiling using high throughput sequencing is a relatively low cost means of expression profiling with deep coverage. However the extent to which very short cDNA sequences can be effectively used in screening for candidate genes is unclear. Here we investigate this question using an evolutionarily distant as well as a closely related transcriptome for referencing tags. We do this by comparing differentially expressed genes and processes between two closely related allopolyploid species of Pachycladon which have distinct altitudinal preferences in the New Zealand Southern Alps. We validate biological inferences against earlier microarray analyses. RESULTS: Statistical and gene annotation enrichment analyses of tag profiles identified more differentially expressed genes of potential adaptive significance than previous analyses of array-based expression profiles. These include genes involved in glucosinolate metabolism, flowering time, and response to cold, desiccation, fungi and oxidation. In addition, despite the short length of 20mer tags, we were able to infer patterns of homeologous gene expression for 700 genes in our reference library of 7,128 full-length Pachycladon ESTs. We also demonstrate that there is significant information loss when mapping tags to the non-conspecific reference transcriptome of A. thaliana as opposed to P. fastigiatum ESTs but also describe mapping strategies by which the larger collection of A. thaliana ESTs can be used as a reference. CONCLUSION: When coupled with a reference transcriptome generated using RNA-seq, tag sequencing offers a promising approach for screening natural populations and identifying candidate genes of potential adaptive significance. We identify computational issues important for the successful application of tag profiling in a non-model allopolyploid plant species.

Cutoffs and k-mers: implications from a transcriptome study in allopolyploid plants.

BACKGROUND: Transcriptome analysis is increasingly being used to study the evolutionary origins and ecology of non-model plants. One issue for both transcriptome assembly and differential gene expression analyses is the common occurrence in plants of hybridisation and whole genome duplication (WGD) and hybridization resulting in allopolyploidy. The divergence of duplicated genes following WGD creates near identical homeologues that can be problematic for de novo assembly and also reference based assembly protocols that use short reads (35 - 100 bp). RESULTS: Here we report a successful strategy for the assembly of two transcriptomes made using 75 bp Illumina reads from Pachycladon fastigiatum and Pachycladon cheesemanii. Both are allopolyploid plant species (2n = 20) that originated in the New Zealand Alps about 0.8 million years ago. In a systematic analysis of 19 different coverage cutoffs and 20 different k-mer sizes we showed that i) none of the genes could be assembled across all of the parameter space ii) assembly of each gene required an optimal set of parameter values and iii) these parameter values could be explained in part by different gene expression levels and different degrees of similarity between genes. CONCLUSIONS: To obtain optimal transcriptome assemblies for allopolyploid plants, k-mer size and k-mer coverage need to be considered simultaneously across a broad parameter space. This is important for assembling a maximum number of full length ESTs and for avoiding chimeric assemblies of homeologous and paralogous gene copies.

Transcript and protein profiling identify candidate gene sets of potential adaptive significance in New Zealand Pachycladon.

BACKGROUND: Transcript profiling of closely related species provides a means for identifying genes potentially important in species diversification. However, the predictive value of transcript profiling for inferring downstream-physiological processes has been unclear. In the present study we use shotgun proteomics to validate inferences from microarray studies regarding physiological differences in three Pachycladon species. We compare transcript and protein profiling and evaluate their predictive value for inferring glucosinolate chemotypes characteristic of these species. RESULTS: Evidence from heterologous microarrays and shotgun proteomics revealed differential expression of genes involved in glucosinolate hydrolysis (myrosinase-associated proteins) and biosynthesis (methylthioalkylmalate isomerase and dehydrogenase), the interconversion of carbon dioxide and bicarbonate (carbonic anhydrases), water use efficiency (ascorbate peroxidase, 2 cys peroxiredoxin, 20 kDa chloroplastic chaperonin, mitochondrial succinyl CoA ligase) and others (glutathione-S-transferase, serine racemase, vegetative storage proteins, genes related to translation and photosynthesis). Differences in glucosinolate hydrolysis products were directly confirmed. Overall, prediction of protein abundances from transcript profiles was stronger than prediction of transcript abundance from protein profiles. Protein profiles also proved to be more accurate predictors of glucosinolate profiles than transcript profiles. The similarity of species profiles for both transcripts and proteins reflected previously inferred phylogenetic relationships while glucosinolate chemotypes did not. CONCLUSIONS: We have used transcript and protein profiling to predict physiological processes that evolved differently during diversification of three Pachycladon species. This approach has also identified candidate genes potentially important in adaptation, which are now the focus of ongoing study. Our results indicate that protein profiling provides a valuable tool for validating transcript profiles in studies of adaptive divergence.

Within and between whorls: comparative transcriptional profiling of Aquilegia and Arabidopsis.

BACKGROUND: The genus Aquilegia is an emerging model system in plant evolutionary biology predominantly because of its wide variation in floral traits and associated floral ecology. The anatomy of the Aquilegia flower is also very distinct. There are two whorls of petaloid organs, the outer whorl of sepals and the second whorl of petals that form nectar spurs, as well as a recently evolved fifth whorl of staminodia inserted between stamens and carpels. METHODOLOGY/PRINCIPAL FINDINGS: We designed an oligonucleotide microarray based on EST sequences from a mixed tissue, normalized cDNA library of an A. formosa x A. pubescens F2 population representing 17,246 unigenes. We then used this array to analyze floral gene expression in late pre-anthesis stage floral organs from a natural A. formosa population. In particular, we tested for gene expression patterns specific to each floral whorl and to combinations of whorls that correspond to traditional and modified ABC model groupings. Similar analyses were performed on gene expression data of Arabidopsis thaliana whorls previously obtained using the Ath1 gene chips (data available through The Arabidopsis Information Resource). CONCLUSIONS/SIGNIFICANCE: Our comparative gene expression analyses suggest that 1) petaloid sepals and petals of A. formosa share gene expression patterns more than either have organ-specific patterns, 2) petals of A. formosa and A. thaliana may be independently derived, 3) staminodia express B and C genes similar to stamens but the staminodium genetic program has also converged on aspects of the carpel program and 4) staminodia have unique up-regulation of regulatory genes and genes that have been implicated with defense against microbial infection and herbivory. Our study also highlights the value of comparative gene expression profiling and the Aquilegia microarray in particular for the study of floral evolution and ecology.

An approach to transcriptome analysis of non-model organisms using short-read sequences.

Transcriptome analysis using high-throughput short-read sequencing technology is straightforward when the sequenced genome is the same species or extremely similar to the reference genome. We present an analysis approach for when the sequenced organism does not have an already sequenced genome that can be used for a reference, as will be the case of many non-model organisms. As proof of concept, data from Solexa sequencing of the polyploid plant Pachycladon enysii was analysed using our approach with its nearest model reference genome being the diploid plant Arabidopsis thaliana. By using a combination of mapping and de novo assembly tools we could determine duplicate genes belonging to one or other of the genome copies. Our approach demonstrates that transcriptome analysis using high-throughput short-read sequencing need not be restricted to the genomes of model organisms.

Convergence, constraint and the role of gene expression during adaptive radiation: floral anthocyanins in Aquilegia.

Convergent phenotypes are testament to the role of natural selection in evolution. However, little is known about whether convergence in phenotype extends to convergence at the molecular level. We use the independent losses of floral anthocyanins in columbines (Aquilegia) to determine the degree of molecular convergence in gene expression across the anthocyanin biosynthetic pathway (ABP). Using a phylogeny of the North American Aquilegia clade, we inferred six independent losses of floral anthocyanins. Via semiquantitative reverse transcriptase-polymerase chain reaction (RT-PCR), we monitored developmental and tissue-specific variation in expression of the six major structural ABP loci in three Aquilegia species, two that produce anthocyanins (A+) and one that does not (A-). We then compared ABP expression in petals of old-bud and pre-anthesis flowers of 13 Aquilegia species, eight wild species and two horticultural lines representing seven independent A- lineages as well as three wild A+ species. We only found evidence of down-regulation of ABP loci in A- lineages and losses of expression were significantly more prevalent for genes late in the pathway. Independent contrast analysis indicates that changes in expression of dihydroflavonol reductase (DFR) and anthocyanidin synthase (ANS) are strongly phylogenetically correlated consistent with the multilocus targets of trans-regulatory elements in the ABP of other systems. Our findings strongly suggest that pleiotropy constrains the evolution of loss of floral anthocyanins to mutations affecting genes late in the ABP mostly through convergent changes in regulatory genes. These patterns support the hypothesis that rapid evolutionary change occurs largely through regulatory rather than structural mutations.

Specificity in ecological interactions: attack from the same lepidopteran herbivore results in species-specific transcriptional responses in two solanaceous host plants.

Model systems have proven enormously useful in elucidating the biochemical function of plant genes. However their ecological function, having been sculpted by evolutionary forces specific to a species, may be less conserved across taxa. Responses to wounding and herbivore attack differ among plant families and are known to be mediated by oxylipin, ethylene, and systemin-signaling networks. We analyzed transcriptional responses of two native Solanaceous species to the attack of an herbivore whose elicitors are known not to be influenced by diet. With The Institute for Genomic Research 10k-cDNA potato (Solanum tuberosum) microarray, we compared the transcriptional responses of Nicotiana attenuata with those of black nightshade (Solanum nigrum) when both were attacked by the Solanaceous generalist herbivore, Manduca sexta. Based on an NADH dehydrogenase subunit F phylogeny, S. nigrum is more closely related to potato than N. attenuata but responded significantly less to M. sexta attack. Apart from transcriptional differences anticipated from their differences in secondary metabolism, both species showed distinct transcriptional patterns (with only 10% overlap in significantly regulated genes), which point to fundamental differences in the signaling cascades and downstream genes mediating herbivore resistance. The lackluster transcriptional response of S. nigrum could not be attributed to its inability to respond to elicitation, because methyl jasmonate elicitation of S. nigrum resulted in a strong transcriptional response. Given that attack from the same herbivore elicits profoundly different responses in two Solanaceaous taxa, we conclude that blueprints for commonly regulated responses to plant-herbivore interactions appear unlikely.

Herbivore-induced plant vaccination. Part II. Array-studies reveal the transience of herbivore-specific transcriptional imprints and a distinct imprint from stress combinations.

Summary Microarray technology has given plant biologists the ability to simultaneously monitor changes in the expression of hundreds of genes, and yet, to date, this technology has not been applied to ecological phenomena. In native tobacco (Nicotiana attenuata), prior attack of sap-feeding mirids (Tupiocoris notatus) results in vaccination of the plant against subsequent attacks by chewing hornworms (Manduca sexta). This vaccination is mediated by a combination of direct and indirect defenses and tolerance responses, which act in concert with the attack preferences of a generalist predator. Here, we use microarrays enriched in herbivore-elicited genes with a principal components analysis (PCA) to characterize transcriptional 'imprints' of single, sequential, or simultaneous attacks by these two main herbivores of N. attenuata. The PCA identified distinctly different imprints left by individual attack from the two species after 24 h, but not after 5 days. Moreover, imprints of sequential or simultaneous attacks differed significantly from those of single attack, suggesting the existence of a distinct gene expression program responsive to the combination of biological stressors. A dissection of the transcriptional imprints revealed responses in direct and indirect defense genes that were well correlated with observed increases in defense metabolites. Attack from both herbivores elicits a switch from growth- to defense-related transcriptional processes, and herbivore-specific changes occur largely in primary metabolism and signaling cascades. PCA of these polygenic transcriptional imprints characterizes the ephemeral changes in the transcriptome that occur during the maturation of ecologically relevant phenotypic responses.

Herbivore-induced ethylene burst reduces fitness costs of jasmonate- and oral secretion-induced defenses in Nicotiana attenuata.

Specialist herbivores are known to alter their host's wound-induced responses but the beneficiaries of these alterations are unknown. Nicotiana attenuata plants release a burst of ethylene specifically in response to feeding by Manduca sexta larvae, which is known to suppress wound- and methyl jasmonate (MeJA)-inducible nicotine accumulation. The ethylene burst may be a mechanism by which M. sexta larvae feed "stealthily" on their host plants or, alternatively, it may allow the plant to optimize its defense response against this specialist herbivore by reducing costs of induction. We examined the impact of the ethylene burst on defense-related fitness costs that are readily observed when plants are treated with MeJA and grown in competition with untreated plants. We elicited nicotine induction (with MeJA), the ethylene burst (with the ethylene releasing compound, ethephon) and inhibited the plant's ability to perceive ethylene (with applications of an antagonist of ethylene receptors, 1-methylcyclopropene, 1-MCP). By simultaneously applying MeJA and ethephon we mimicked the plant's hormonal response to larval attack. We hypothesized that if the ethylene burst benefited the plant, the fitness costs of MeJA induction should be reduced by ethephon and restored if the plants were additionally treated with 1-MCP. In a second experiment, we applied larval oral secretion (OS) to elicit endogenous hormone production and predicted that the 1-MCP treatment should reduce the fitness of OS-treated plants. Our measures of plant fitness, namely the rate of stalk elongation and lifetime capsule production, supported these predictions. We conclude that the ethylene burst elicited by this specialist herbivore can reduce MeJA-induced fitness costs and increase the competitive strength of OS-treated plants. Suppressed nicotine production is likely to contribute to, but is not sufficient to explain, the observed fitness outcomes. The intensity of intra-specific competition and herbivore attack will likely determine the adaptive value of the M. sexta-elicited ethylene response.