Microarrays in ecological research: A case study of a cDNA microarray for plant-herbivore interactions

BACKGROUND: Microarray technology allows researchers to simultaneously monitor changes in the expression ratios (ERs) of hundreds of genes and has thereby revolutionized most of biology. Although this technique has the potential of elucidating early stages in an organism's phenotypic response to complex ecological interactions, to date, it has not been fully incorporated into ecological research. This is partially due to a lack of simple procedures of handling and analyzing the expression ratio (ER) data produced from microarrays. RESULTS: We describe an analysis of the sources of variation in ERs from 73 hybridized cDNA microarrays, each with 234 herbivory-elicited genes from the model ecological expression system, Nicotiana attenuata, using procedures that are commonly used in ecologic research. Each gene is represented by two independently labeled PCR products and each product was arrayed in quadruplicate. We present a robust method of normalizing and analyzing ERs based on arbitrary thresholds and statistical criteria, and characterize a "norm of reaction" of ERs for 6 genes (4 of known function, 2 of unknown) with different ERs as determined across all analyzed arrays to provide a biologically-informed alternative to the use of arbitrary expression ratios in determining significance of expression. These gene-specific ERs and their variance (gene CV) were used to calculate array-based variances (array CV), which, in turn, were used to study the effects of array age, probe cDNA quantity and quality, and quality of spotted PCR products as estimates of technical variation. Cluster analysis and a Principal Component Analysis (PCA) were used to reveal associations among the transcriptional "imprints" of arrays hybridized with cDNA probes derived from mRNA from N. attenuata plants variously elicited and attacked by different herbivore species and from three congeners: N. quadrivalis, N. longiflora and N. clevelandii. Additionally, the PCA revealed the contribution of individual gene ERs to the associations among arrays. CONCLUSIONS: While the costs of 'boutique' array fabrication are rapidly declining, familiar methods for the analysis of the data they create are still missing. The case history illustrated here demonstrates the ease with which this powerful technology can be adapted to ecological research

Two-fold differences are the detection limit for determining transgene copy numbers in plants by real-time PCR

BACKGROUND: After transformation, plants that are homozygous and contain one copy of the transgene are typically selected for further study. If real-time PCR is to be used to determine copy number and zygosity, it must be able to distinguish hemizygous from homozygous and one-copy from two-copy plants. That is, it must be able to detect two-fold differences. RESULTS: When transgenic Nicotiana attenuata plants which had been previously determined by Southern analysis to contain one or two copies of the transgene, were analyzed by real-time PCR (2(-ΔΔCt )method), the method failed to confirm the results from the Southern analysis. In a second data set we analyzed offspring of a hemizygous one-copy plant, which were expected to segregate into three groups of offspring in a 1:2:1 ratio: no transgene, hemizygous, homozygous. Because it was not possible to distinguish homozygous from hemizygous plants with real-time PCR, we could not verify this segregation ratio. CONCLUSIONS: Detection of two-fold differences by real-time PCR is essential if this procedure is to be used for the characterization of transgenic plants. However, given the high variability between replicates, a detection of two-fold differences is in many cases not possible; in such cases Southern analysis is the more reliable procedure

A persistent major mutation in canonical jasmonate signaling is embedded in an herbivory-elicited gene network

When insect herbivores attack plants, elicitors from oral secretions and regurgitants (OS) enter wounds during feeding, eliciting defense responses. These generally require plant jasmonate (JA) signaling, specifically, a jasmonoyl-L-isoleucine (JA-Ile) burst, for their activation and are well studied in the native tobacco Nicotiana attenuata. We used intraspecific diversity captured in a 26-parent MAGIC population planted in nature and an updated genome assembly to impute natural variation in the OS-elicited JA-Ile burst linked to a mutation in the JA-Ile biosynthetic gene NaJAR4. Experiments revealed that NaJAR4 variants were associated with higher fitness in the absence of herbivores but compromised foliar defenses, with two NaJAR homologues (4 and 6) complementing each other spatially and temporally. From decade-long seed collections of natural populations, we uncovered enzymatically inactive variants occurring at variable frequencies, consistent with a balancing selection regime maintaining variants. Integrative analyses of OS-induced transcriptomes and metabolomes of natural accessions revealed that NaJAR4 is embedded in a nonlinear complex gene coexpression network orchestrating responses to OS, which we tested by silencing four hub genes in two connected coexpressed networks and examining their OS-elicited metabolic responses. Lines silenced in two hub genes (NaGLR and NaFB67) co-occurring in the NaJAR4/6 module showed responses proportional to JA-Ile accumulations; two from an adjacent module (NaERF and NaFB61) had constitutively expressed defenses with high resistance. We infer that mutations with large fitness consequences can persist in natural populations due to compensatory responses from gene networks, which allow for diversification in conserved signaling pathways and are generally consistent with predictions of an omnigene model

Molecular evolution and diversification of the Argonaute family of proteins in plants

BACKGROUND: Argonaute (AGO) proteins form the core of the RNA-induced silencing complex, a central component of the smRNA machinery. Although reported from several plant species, little is known about their evolution. Moreover, these genes have not yet been cloned from the ecological model plant, Nicotiana attenuata, in which the smRNA machinery is known to mediate important ecological traits. RESULTS: Here, we not only identify 11 AGOs in N. attenuata, we further annotate 133 genes in 17 plant species, previously not annotated in the Phytozome database, to increase the number of plant AGOs to 263 genes from 37 plant species. We report the phylogenetic classification, expansion, and diversification of AGOs in the plant kingdom, which resulted in the following hypothesis about their evolutionary history: an ancestral AGO underwent duplication events after the divergence of unicellular green algae, giving rise to four major classes with subsequent gains/losses during the radiation of higher plants, resulting in the large number of extant AGOs. Class-specific signatures in the RNA-binding and catalytic domains, which may contribute to the functional diversity of plant AGOs, as well as context-dependent changes in sequence and domain architecture that may have consequences for gene function were found. CONCLUSIONS: Together, the results demonstrate that the evolution of AGOs has been a dynamic process producing the signatures of functional diversification in the smRNA pathways of higher plants. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12870-014-0364-6) contains supplementary material, which is available to authorized users

Development of an inducible protein expression system based on the protozoan host Leishmania tarentolae

Production of functional eukaryotic proteins in recombinant form is a bottle-neck in various post-genomic applications and in life science in general. At least partially this is due to the problems associated with the use of endogenous RNA polymerase II for high-level transcription of heterologous genes in eukaryotic expression systems. To circumvent these problems we developed a new inducible protein expression system based on the protozoan host Leishmania tarentolae (Trypanosomatidae). We have created a strain of L. tarentolae constitutively co-expressing T7 RNA polymerase and tetracycline repressor. This strain could be stably transformed with the heterologous target gene under control of the T7 promoter/TET operator assembly, which can initiate transcription upon addition of tetracycline to the culture medium. Using this system, we demonstrated that enhanced green fluorescent protein (EGFP) could be overexpressed to a level of ca. 1% of total cellular protein. The developed system was tested for its ability to inducibly co-express multiple genes. Using two copies of the egfp gene integrated at two different genomic sites, we could obtain expression levels reaching 4% of total cellular protein. Further possible improvements and applications of the developed system are discussed

Supplemental Table S2

Results of diagnostic PCRs with 44 different transgenic N. attenuata lines

Ectopic Expression of AtJMT in Nicotiana attenuata: Creating a Metabolic Sink Has Tissue-Specific Consequences for the Jasmonate Metabolic Network and Silences Downstream Gene Expression1[W][OA]

To create a metabolic sink in the jasmonic acid (JA) pathway, we generated transgenic Nicotiana attenuata lines ectopically expressing Arabidopsis (Arabidopsis thaliana) jasmonic acid O-methyltransferase (35S-jmt) and additionally silenced in other lines the N. attenuata methyl jasmonate esterase (35S-jmt/ir-mje) to reduce the deesterification of methyl jasmonate (MeJA). Basal jasmonate levels did not differ between transgenic and wild-type plants; however, after wounding and elicitation with Manduca sexta oral secretions, the bursts of JA, jasmonoyl-isoleucine (JA-Ile), and their metabolites that are normally observed in the lamina, midvein, and petiole of elicited wild-type leaves were largely absent in both transformants but replaced by a burst of endogenous MeJA that accounted for almost half of the total elicited jasmonate pools. In these plants, MeJA became a metabolic sink that affected the jasmonate metabolic network and its spread to systemic leaves, with major effects on 12-oxo-phytodieonic acid, JA, and hydroxy-JA in petioles and on JA-Ile in laminas. Alterations in the size of jasmonate pools were most obvious in systemic tissues, especially petioles. Expression of threonine deaminase and trypsin proteinase inhibitor, two JA-inducible defense genes, was strongly decreased in both transgenic lines without influencing the expression of JA biosynthesis genes that were uncoupled from the wounding and elicitation with M. sexta oral secretions-elicited JA-Ile gradient in elicited leaves. Taken together, this study provides support for a central role of the vasculature in the propagation of jasmonates and new insights into the versatile spatiotemporal characteristics of the jasmonate metabolic network

RNA-Directed RNA Polymerase3 from Nicotiana attenuata Is Required for Competitive Growth in Natural Environments1[W][OA]

SDE1/SGS2/RdR6, a putative RNA-directed RNA polymerase, maintains plant defenses against viruses in Arabidopsis (Arabidopsis thaliana) and Nicotiana benthamiana, but its function has not been examined in natural habitats or with respect to other ecological stresses. We evaluated the organismic-level function of this gene (NaRdR3) in an ecological model species, Nicotiana attenuata, by transforming plants to stably silence RdR3 (irRdR3). Minor morphological changes (elongated leaves and reduced leaf number) and increased susceptibility to tobamoviruses typical of RdR6 silencing in other species were observed, but these changes did not alter the reproductive performance of singly grown plants (measured as seed and capsule production) or herbivore resistance in laboratory trials. 454-sequencing of irRdR3's small RNA (smRNA) transcriptome revealed that 21- and 24-nucleotide smRNAs were not affected, but the abundance of 22- to 23-nucleotide smRNAs was reduced. When planted in pairs with wild-type plants in N. attenuata's natural habitat in the Great Basin Desert, irRdR3 plants produced shorter stalks with significantly reduced flower and capsule numbers, but did not influence the ability of plants to resist the native herbivore community, indicating that silencing RdR3 reduced a plant's competitive ability. We tested this hypothesis in the glasshouse by planting irRdR3 and wild-type pairs in communal containers; again irRdR3 plants had severely reduced stalk elongation and reproductive measures. The reduced competitive ability of irRdR3 plants was associated with altered phytohormone homeostasis, especially as reflected in the distribution of auxin. We suggest that RdR3 helps to regulate hormone balance when plants compete with conspecifics in natural environments