Infiltration-RNAseq: transcriptome profiling of Agrobacterium-mediated infiltration of transcription factors to discover gene function and expression networks in plants
© The Author(s) 2016
Received: 9 June 2016
Accepted: 4 October 2016
Published: 19 October 2016
Transcription factors (TFs) coordinate precise gene expression patterns that give rise to distinct phenotypic outputs. The identification of genes and transcriptional networks regulated by a TF often requires stable transformation and expression changes in plant cells. However, the production of stable transformants can be slow and laborious with no guarantee of success. Furthermore, transgenic plants overexpressing a TF of interest can present pleiotropic phenotypes and/or result in a high number of indirect gene expression changes. Therefore, fast, efficient, high-throughput methods for assaying TF function are needed.
Agroinfiltration is a simple plant biology method that allows transient gene expression. It is a rapid and powerful tool for the functional characterisation of TF genes in planta. High throughput RNA sequencing is now a widely used method for analysing gene expression profiles (transcriptomes). By coupling TF agroinfiltration with RNA sequencing (named here as Infiltration-RNAseq), gene expression networks and gene function can be identified within a few weeks rather than many months. As a proof of concept, we agroinfiltrated Medicago truncatula leaves with M. truncatula LEGUME ANTHOCYANIN PRODUCITION 1 (MtLAP1), a MYB transcription factor involved in the regulation of the anthocyanin pathway, and assessed the resulting transcriptome. Leaves infiltrated with MtLAP1 turned red indicating the production of anthocyanin pigment. Consistent with this, genes encoding enzymes in the anthocyanin biosynthetic pathway, and known transcriptional activators and repressors of the anthocyanin biosynthetic pathway, were upregulated. A novel observation was the induction of a R3-MYB transcriptional repressor that likely provides transcriptional feedback inhibition to prevent the deleterious effects of excess anthocyanins on photosynthesis.
Infiltration-RNAseq is a fast and convenient method for profiling TF-mediated gene expression changes. We utilised this method to identify TF-mediated transcriptional changes and TF target genes in M. truncatula and Nicotiana benthamiana. This included the identification of target genes of a TF not normally expressed in leaves, and targets of TFs from other plant species. Infiltration-RNAseq can be easily adapted to other plant species where agroinfiltration protocols have been optimised. The ability to identify downstream genes, including positive and negative transcriptional regulators, will result in a greater understanding of TF function.
KeywordsTranscription factor Agrobacterium tumefaciens Infiltration Medicago truncatula MtLAP1 Anthocyanin RNAseq Transcriptome Nicotiana benthamiana
Transcription factors (TFs) play a major role in plant development and their response to the environment. An understanding of the function(s) of TFs has been elucidated by a number of reverse genetic tools, including gene knockout and gene overexpression strategies. As some TFs display unique properties, the overexpression strategy has been particularly effective in revealing TF function [1–3]. However, these methods rely on stable transformation of plant cells, but obtaining stable transformants of crop plants is a time-consuming and labour-intensive process that is often inefficient.
In contrast to making stable transgenics, Agrobacterium-mediated transient expression of genes provides a rapid and simple alternative for analysing gene function . The most commonly used transient expression method involves ‘agroinfiltration’—injection of an Agrobacterium suspension (harbouring a gene of interest on a T-DNA vector) through the stomata and into the mesophyll of expanded leaves using a needleless syringe. Copies of the T-DNA are transferred from the Agrobacterium into the leaf parenchyma cells. Although only a low number of T-DNA copies integrate into the plant chromosomes, the non-integrated T-DNAs are transiently expressed for several days .
Leaves of Nicotiana benthamiana or Nicotiana tabacum have proven to be a reliable system for agroinfiltration, where a large fraction of cells are transformed, and in extreme cases can result in 50 % of the total soluble leaf protein encoded by the transferred gene . This has led to applications where pharmaceutically active proteins are produced via agroinfiltration at a commercially viable scale . In addition to bulk protein production, agroinfiltration has been used to address a range of biological questions. For example, it provides an in vivo system to identify protein-protein interactions , examine protein localization , study host-pathogen interactions  and understand TFs that regulate specific promoter sequences [11–13] and specific gene expression profiles . Detailed aspects of gene regulation, such as the role of upstream open reading frames (uORFs), introns within 5′-untranslated regions (5′UTRs) and non-canonical translation initiation have also been investigated via agroinfiltration [15–17]. The key advantage of agroinfiltration over stable transgenics is that experimental results can be generated after a few days, rather than many months, and compared to using stably transformed plant lines, transient expression assays eliminate variation due to different chromosomal positions and epigenetic states of the transformed constructs.
Although leaves of Nicotiana species are widely used for agroinfiltration, other plants have been used with variable success [18–24]. We have recently established a robust method for agroinfiltration of Medicago truncatula , an important model legume that is used extensively by researchers throughout the world. Adapting the agroinfiltration method for use in Medicago provides a more relevant tool for researchers of legume species, as the production of stably transformed leguminous plants is difficult and time consuming. A powerful use of the transient assay system in this species is the identification of direct targets of a specific TF . For example, gene expression analysis via quantitative real-time PCR demonstrated that candidate downstream targets of M. truncatula LEGUME ANTHOCYANIN PRODUCITON 1 (MtLAP1), a MYB TF involved in the regulation of the anthocyanin biosynthesis pathway , were upregulated in leaves that were agroinfiltrated with 35S:MtLAP1 .
Constitutive expression of MtLAP1 in transgenic or agroinfiltrated M. truncatula induces massive accumulation of anthocyanin pigments [25, 26]. Microarray analysis of transgenic Medicago plants over-expressing MtLAP1 found over 70 genes were up-regulated, many of which were involved in anthocyanin biosynthesis . However, in addition to the time consuming nature of producing transgenic plants, stable overexpression of TFs can result in a high number of indirect gene expression changes being called, and pleiotropic phenotypes that need to be interpreted with caution [1–3]. We propose a more simple experiment to understand TF-mediated gene expression changes that does not rely on the production of stable transformants—agroinfiltration of a TF followed by high-throughput RNA sequencing. Due to the short time frame of this type of experiment, we speculate that the gene expression analysis will be less noisy and allow greater discovery of specific downstream genes and thus gene expression networks.
In this article, we present the method of ‘Infiltration-RNAseq’—a simple and rapid analysis of gene expression profiles that allows the discovery of novel downstream targets of a chosen TF. In order to establish this technique, RNA sequencing analysis was performed on M. truncatula leaves agroinfiltrated with 35S:MtLAP1, the primary TF involved in the anthocyanin biosynthetic pathway [25, 26]. Infiltration-RNAseq was further validated using Nicotiana benthamaina leaves to determine the transcriptional changes resulting from infiltration of Agrobacterium either alone or expressing TFs from another plant or not normally expressed in leaves.
Infiltration-RNAseq: for discovery of downstream targets of transcription factors
Infiltration of Agrobacterium tumefaciens containing the 35S:MtLAP1 construct results in high MtLAP1 expression in M. truncatula . Medicago leaves overexpressing the MtLAP1 TF turn red (Fig. 1b; [25, 26]), indicating the anthocyanin biosynthetic pathway is upregulated. Leaves from plants agroinfiltrated with 35S:MtLAP1 or control constructs were harvested 4 days post-infiltration (Fig. 1a; “Methods”). While high MtLAP1 expression was observed 3 days post infiltration , we wanted to capture the highest gene expression changes of MtLAP1 target genes. Therefore, we reasoned that maximal translation of MtLAP1 mRNA and thus function of MtLAP1 protein would occur about a day later. As with infiltration of N. benthamiana, it is not recommended to sample more than 5 days post-infiltration as RNA silencing mechanisms will act on the agroinfiltrated construct and result in low reproducibility between samples and experiments [12, 27]. RNAseq libraries were made with RNA from the different leaf samples and the resulting RNAseq data was analysed to determine the gene expression changes (see “Methods”). We have named this method of gene expression analysis and discovery ‘Infiltration-RNAseq’.
Transcriptional regulation of anthocyanin biosynthetic genes in Medicago truncatula leaves via agroinfiltration of 35S:MtLAP1
Diverse and stable anthocyanin pigments are derived from anthocyanidins that are immediately modified by glycosylation, acylation, and/or methylation . Overexpression of MtLAP1 induced a number of glucosyltransferases and malonyltransferases (anthocyanin acyl transferases), and a glutathione-S-transferase (Additional file 1: Table S1; ) that are likely to be responsible for these modifications. Modification of the anthocyanidin compounds is essential for their transport into the vacuole where the water-soluble anthocyanin pigments are stored. Transport of diverse flavonoid compounds into the vacuole involves Multidrug and Toxic compound Extrusion (MATE) proteins [33, 34]. Agroinfiltration of MtLAP1 into Medicago leaves resulted in a significant 10-fold induction of a MATE transporter (Medtr1g021740) (Additional file 1: Table S1). This MATE transporter has not yet been characterised, but it is distinct from MtMATE1 and MtMATE2, which transport proanthocyanidin monomers and various flavonoid-/anthocyanin-glycosides, respectively [33, 34]. Differential expression analysis via DESeq2  identified MtMATE2 (Medtr1g100180) to be up-regulated 29-fold upon MtLAP1 infiltration (Additional file 5: Figure S1), however MATE2 did not meet the criteria for the stringent Intensity Difference Filter (P value >0.05), most probably due to high variation in expression between the biological replicates for the different conditions.
MtLAP1-regulated transcriptional activators and repressors controlling anthocyanin pigmentation in Medicago truncatula
Known targets of the MYB transcriptional activation complex for anthocyanin biosynthesis include bHLH2 (an AtTT8 clade bHLH). A core MBW transcriptional activation complex containing bHLH2 forms, which activates the expression of anthocyanin biosynthesis genes and reinforces expression of itself . Upon MtLAP1 agroinfiltration, the M. truncatula bHLH2 gene, Medtr1g072320/MtTT8 , displayed a tenfold induction (Fig. 4; Additional file 1: Table S1), suggesting the resulting bHLH2 protein could become part of the core MBW complex. This complex can also activate the expression of the R2R3-MYB repressors [37, 39, 40]. We found two R2R3-MYB repressor genes that were upregulated by MtLAP1 agroinfiltration: Medtr4g485530/MtMYB530, homologous to TrMYB133, displayed a 93-fold induction, and Medtr5g079670/MtMYB2, homologous to TrMYB134, displayed a 67-fold induction [40, 41] (Fig. 4; Additional file 1: Table S1). The inclusion of the R2R3-MYB repressor into the core MBW complex results in active transcriptional repression of target promoters of the core MBW complex . Similarly, the core MBW complex activates the expression of a smaller R3-MYB repressor, which provides feedback inhibition by inhibiting the formation of new core MBW complexes by titrating bHLH factors [37, 39]. The Medicago R3-MYB likely to be involved in this process is Medtr2g088730/MtMYB730 , which displayed a 207-fold induction in MtLAP1 agroinfiltrated leaves (Fig. 4; Additional file 1: Table S1). The deduced amino acid sequence of this gene contains the bHLH-interacting motif , which is highly predictive for binding anthocyanin-related bHLH partners, and is phylogenetically related to characterised R3-MYB proteins that regulate anthocyanin/proanthocyanidin synthesis (Additional file 6: Figure S2). This suggests feedback repression for anthocyanin synthesis by R3-MYBs occurs in leguminous species.
Our Infiltration-RNAseq experiment was successful in identifying potential positive and negative regulators of anthocyanin biosysnthesis in M. truncatula. Importantly, all of the genes encoding the positive and negative transcriptional regulators of the anthocyanin pathway were present in cluster 2 (Fig. 2b; Additional file 3: Table S3), which includes genes that had little or no expression in control samples, but were greatly induced upon MtLAP1 agroinfiltration. This pattern of gene expression is consistent with the induction of visible anthocyanin pigmentation upon infiltration with MtLAP1 (Fig. 1b), as expected for a probable transcriptional activator. The positive and negative transcriptional regulators of this system, especially the R3-MYB repressor, act to provide homeostasis for the cell—in leaves, anthocyanins can provide photo-protection from excess light, but excessive anthocyanin accumulation can compromise photosynthesis .
Other genes regulated by MtLAP1 overexpression
Gene ontology enrichments  for the remaining significantly differentially expressed genes (Additional file 1: Table S1) primarily identified the following processes: carbohydrate metabolism; hydrolase activity; acting on glycosyl bonds. Anthocyanin molecules are decorated with sugar molecules (such as glucose, rhamnose, galactose, arabinose and xylose), which act to stabilise the anthocyanins . We speculate that MtLAP1-induced genes encoding enzymes involved in hydrolysis of glycosidic bonds between two or more carbohydrates could result in the supply of sugar molecules for the glycosylation process.
Homologues of MtLAP1 (AtMYB75/AtPAP1) have been shown to act as a repressor of lignin biosynthetic genes and secondary cell-wall metabolism genes [45–47]. Altered expression of MYB75/PAP1 results in a major shift of the carbon partitioning to major sinks . The Arabidopsis myb75-1 mutant has an increase in carbon accumulation of monolignols, precursors of lignin in secondary cell walls, which results in increased cell wall thickness. In contrast, a gain-of-function mutation of MYB75 results in an increase in carbon accumulation in anthocyanins. A possible source of the carbon needed to produce the increased number of anthocyanin pigments in response to the overexpression of MtLAP1 is the plant cell wall. MtLAP1 agroinfiltration induced the expression of a number of genes encoding enzymes associated with, or the breakdown of, the cell wall (such as expansion and pectinesterase). We speculate that these enzymes might contribute to supplying the carbon to the cytoplasm for MtLAP1-induced anthocyanin production. However, further work is needed to understand the role of these genes, if any, in MtLAP1-induced anthocyanin production.
Infiltration-RNAseq in Nicotiana benthamiana
We have also utilised Infiltration-RNAseq via agroinfiltration of N. benthamiana. Here, we investigated three questions: (1) what are the transcriptional changes induced when Agrobacterium is infiltrated into a leaf, compared to infiltration of buffer only; (2) can Infiltration-RNAseq identify the targets of a TF from another plant species; (3) can the targets of a TF not normally expressed in leaves be identified using Infiltration-RNAseq.
Infiltration of Agrobacterium containing no construct into N. benthamiana leaves resulted in significant differential expression of a large number of genes (compared to buffer only; Additional file 7: Table S5). Upregulated genes included those involved in pathogenesis, and systemic acquired resistance (SAR) to overcome the perceived pathogen attack, such as pathogenesis-related (PR) proteins, WRKY transcription factors and disease resistance proteins. Agrobacterium infiltration resulted in the downregulation of chlorophyll binding proteins and photosynthetic metabolism genes. These transcriptional changes are likely the cause of the mild necrotic responses and leaf chlorosis symptoms elicited by disarmed Agrobacterium lab strains , which highlights the importance of using appropriate controls for assessing transcriptional changes via Infiltration-RNAseq (see “Methods”).
To determine if Infiltration-RNAseq can be used to identify the targets of a TF from another plant species, N. benthamiana leaves were agroinfiltrated with kiwifruit (Actinidia chinensis) AcMYB10 , an anthocyanin regulator homologous to MtLAP1 and AtPAP1, and compared to infiltration with Agrobacterium containing no construct. This experiment identified anthocyanin biosynthetic genes and transcriptional regulators were significantly upregulated (Additional file 8: Table S6), which supports the findings from the Infiltration-RNAseq analysis of Medicago leaves infiltrated with MtLAP1 (Fig. 3; Additional file 1: Table S1). In addition, it highlights that Infiltration-RNAseq can be employed in a range of plant species and targets of TFs from non-model plants can be identified in plants with optimised agroinfiltration protocols.
Our third experiment with N. Benthamiana was performed to ascertain whether Infiltration-RNAseq can identify the targets of a TF not normally expressed in leaves. To do this, we agroinfiltrated the Arabidopsis thaliana seed specific TF, LEAFY COTYLEDON2 (AtLEC2) that plays a critical role during both early and late embryo development ; it acts to repress leaf traits and premature germination, and it is involved in the activation of seed storage proteins . Compared to infiltration of Agrobacterium alone, agroinfiltration of AtLEC2 resulted in the increased expression of seed storage proteins, including numerous oleosin proteins and 2S albumin storage protein, and various auxin-responsive genes (Additional file 9: Table S7). This result suggests that ectopic-expression of a TF via agroinfiltration can be used in conjunction with RNAseq to identify downstream target genes.
Agroinfiltration has been extensively used to study the biochemical function of proteins encoded by various plant genes; for example, it has been used to verify protein-protein interactions , identify protein subcellular localization , and to examine the ability of TFs to regulate promoter-reporter gene constructs . However, because infiltration of Agrobacterium into a leaf will cause a large number of transcriptional changes, such as activation of genes involved in the plant’s immune response (Additional file 7: Table S5), it was unclear if this method could be used to reliably identify endogenous genes that might be regulated by the transient overexpression of specific TFs. Here, we show that by comparing the genes induced by Agrobacterium carrying a TF expressed under the control of the 35S promoter, with those induced by Agrobacterium carrying control constructs, we were able to identify TF-specific transcriptional changes.
MtLAP1, the MYB transcription factor that regulates the anthocyanin biosynthetic pathway in Medicago [25, 26], was used to investigate if the downstream targets of a TF could be identified using our Infiltration-RNAseq method. Our results showed that the complex regulation of the anthocyanin pathway was revealed in a single experiment where RNA sequencing data from Medicago leaves infiltrated with MtLAP1 was compared to data obtained from control infiltrations. While studies using model plants have already provided a detailed understanding of the regulation of the anthocyanin pathway , identifying the corresponding genes in other plant species can be challenging. For example, in Medicago, there are more than 16 genes that potentially encode the enzyme Chalcone Synthase (CHS) (Additional file 10: Figure S3), which catalyses the first step in the biosynthesis of flavonoids. However, our results show that only two of these genes were strongly (approximately 30-fold) significantly upregulated by MtLAP1 infiltration, while another three genes displayed weaker (approximately sixfold) significant upregulation (P value <0.05; Fig. 3; Additional file 1: Table S1). Similarly, MtLAP1 infiltration caused the significant (P value <0.05) upregulation of other genes encoding enzymes in the anthocyanin pathway to varying degrees (approximately 8–40-fold) (Fig. 3; Additional file 1: Table S1). As well as identifying biosynthetic genes, genes associated with having higher levels of anthocyanins in the leaf, such as those encoding enzymes that modify or transport anthocyanins into the vacuole, were found to be upregulated (Additional file 1: Table S1). The spatiotemporal regulation of the anthocyanin pathway is precisely controlled by transcriptional activators and repressors from the MYB and basic helix-loop-helix (bHLH) protein families [52–54]. Our Medicago infiltration-RNAseq experiment identified MYBs and bHLH genes that were significantly upregulated by MtLAP1 infiltration and are likely transcriptional activators and repressors of the anthocyanin pathway (Fig. 4; Additional file 1: Table S1). Of particular interest was the discovery of a Medicago R3-MYB, MtMYB730 (Medtr2g088730), which is likely to play a key role in the feedback repression of the MBW transcriptional activation complex. The observation that MtMYB730 expression was activated by MtLAP1 infiltration strongly suggests that feedback repression by R3-MYBs is also conserved in legumes; a feature of the anthocyanin gene regulation network that is proposed to occur throughout the eudicots , if not more widely .
The Infiltration-RNAseq method was further validated using N. benthamiana, a plant commonly used for agroinfiltration experiments. In this experiment, we compared the changes in gene expression caused by infiltrating Agrobacterium into N. benthamiana leaves with those caused by infiltrating buffer alone. As expected, the expression of a large number of genes, such as those encoding pathogenesis related (PR) proteins were significantly upregulated (Additional file 7: Table S5). As infiltration of Agrobacterium increases the expression of plant defence genes, we cannot rule out that this could mask the interpretation of results for TFs involved in defence. We envisage, for some defence-related TFs, the relevant changes in gene expression may not be seen, or the levels of gene expression change may not portray an accurate measure of how the TF effects expression of target genes. As such, we suggest caution should be taken when using Infiltration-RNAseq to understand the function of defence related TFs. This experiment also highlights the importance of using sufficient controls for an Infiltration-RNAseq experiment—we suggest using Agrobacterium either with or without a construct containing a transcriptional regulator that is independent of the TF being studied.
The use of an agroinfiltration control enabled TF-specific gene expression changes to be identified. Like the results obtained when 35S:MtLAP1 was infiltrated into Medicago leaves, the infiltration of the kiwifruit homologue of MtLAP1, AcMYB10, resulted in the significant upregulation of genes involved in anthocyanin biosynthesis and regulation in N. benthamiana (Additional file 8: Table S6). Similarly, infiltration of seed-specific AtLEC2 resulted in genes normally only found in seeds, such as seed storage proteins (legumin, globulin, glutelin) and lipid binding or storage proteins (oleosins, glycine-rich lipid binding proteins, and Lipid transfer protein) being expressed in the infiltrated leaves (Additional file 9: Table S7). Recently, Grimberg et al., used a variation of Infiltration-RNAseq to identify transcriptional transitions in N. benthamiana leaves upon infiltration of five homologs of the TF WRINKLED1, which is involved in oil synthesis . These experiments show that Infiltration-RNAseq has the power to determine the target genes of a TF of interest, even when they are not normally expressed in leaves.
Given the simplicity of Infiltration-RNAseq, this method could be used as a high-throughput way to analyse a large number of TFs. Currently, agroinfiltration is being used as a way to screen for TFs that might be involved in activating the expression of a particular gene or pathway . This involves co-infiltrating leaves with Agrobacterium containing a promoter-reporter construct together with Agrobacterium containing the candidate TF construct. However, often multiple members of a TF family can activate a promoter-reporter gene construct making it difficult to know if the correct TF has been identified. Infiltration-RNAseq could be used to rapidly examine the in planta targets of a number of different TF family members. Our method would also be particularly valuable for studying TFs where a loss- or gain-of-function mutation induces an embryonic lethal phenotype.
Infiltration-RNAseq should be able to be used in any plant species where agroinfiltration protocols have been optimised, such as tomato , lettuce  potato , petunia  grapevine , grapefruit  and the medicinal plant Maesa lanceolata . For many plants, such as Medicago and other legumes, making transgenic plants can be time-consuming and challenging. Infiltration-RNAseq provides another way to understand gene function and gene expression networks in this important family of plants that possess symbiotic nitrogen fixing capabilities.
Infiltration-RNAseq provides a new and rapid way of identifying the targets of TFs. We have shown that this method works in both Medicago and N. benthamiana, that it can be used to identify TF target genes that are not normally expressed in leaves, and that targets of TFs from other plant species can be uncovered. This method should be applicable for any plant species in which transient agroinfiltration expression can be performed. It will be particular valuable for research in species, such as Medicago, lacking simple and reliable transformation protocols.
For Medicago infiltrations, all three leaflets of the 4th trifoliate leaf of healthy 3 weeks old plants were infiltrated with Agrobacterium (strain AGL-1) transformed with the 35S:MtLAP1 construct  or various control constructs (Fig. 1a). The control constructs include the following genes that are independent of the anthocyanin biosynthetic pathway and were being investigated as potential flowering time regulators [56, 57]: Medtr7g018170 (MtCOla), Medtr5g069480 (MtCOlf), Medtr7g084970 (MtFTa1) and Medtr5g022780 (MtFD). The coding sequences were sub-cloned into pCR8/GW/TOPO (Life technologies), checked by sequencing and then cloned into the Gateway compatible binary vector pB2GW7 [56–58]. This different set of experimental controls removes the effects of agroinfiltration i.e. genes involved in the pathogen disease response will be masked as they will be similarly expressed in leaves infiltrated with the control or test constructs. Preparation and delivery of the Agrobacterium suspension was performed as previously described .
For N. benthamiana infiltrations, preparation and delivery of the Agrobacterium suspension was performed as previously described [12, 27]. Briefly, plants were grown until they had six leaves and the youngest leaves over 1 cm long were infiltrated with Buffer, Agrobacterium (strain GV3101:pMP90), or Agrobacterium (strain GV3101:pMP90) transformed with 35S:AcMYB10 or 35S:AtLEC2 . Each type of infiltration was performed in triplicate.
Optimisation of agroinfiltration will be required for different plant species and this will involve understanding when maximal expression of an infiltrated construct(s) occurs. Therefore, as with previous agroinfiltration experiments, we recommend sampling 3–4 days post infiltration and no later than about 5 days post infiltration as silencing of agroinfiltrated constructs is to be expected [12, 27].
Transcriptome analysis of Medicago truncatula agroinfiltrations
Infiltrated leaf tissue was harvested 4 days post infiltration with the intention to capture the highest gene expression changes of MtLAP1 target genes. Total RNA was extracted from the infiltrated leaf tissue using RNeasy Plant Mini Kit (Qiagen). RNAseq libraries were prepared from three control infiltrations (where one leaf was infiltrated with 35S:MtCOla, 35S:MtCOlf or 35S:MtFTa1plus35S:MtFD) and three MtLAP1 infiltrations (where three leaves were infiltrated with 35S:MtLAP1 and named: 35S:MtLAP1_1, 35S:MtLAP1_2 and 35S:MtLAP1_3) using 500 ng of total RNA, and indexed with the TruSeq Stranded mRNA Library Prep Kit according to the protocol (Illumina). The same amount of each RNAseq library was pooled and run on one lane of HiSeq2000 100PE (Illumina).
The quality of the sequencing reads was assessed with FASTQC . RNA sequencing reads were trimmed and filtered on length (≥25 bp) using Cleanadaptors . Sequences (7–13 million reads per library) were aligned to the M. truncatula genome (version Mt4.0) using TopHat2 v2.0.14 , with default parameters and –library-type fr-firststrand, –b2-sensitive (version 18.104.22.168 of Bowtie2) (Additional file 11: Table S8). Visualisation and quantitation of the RNAseq libraries was performed with Seqmonk (v0.32.0) . The RNAseq quantitation pipeline (within SeqMonk) was employed with uniquely mapped reads to generate a set of probes covering every mRNA in the genome, which were quantitated based on the number of reads falling within the exons of those mRNAs (ignoring any reads found in introns). The log-transformed counts produced were corrected for the total number of sequences in each dataset. The three control samples (35S:MtCOla, 35S:MtCOlf and 35S:FTa1plus35S:MtFD) and the three test samples (35S:LAP1_1, 35S:LAP1_2 and 35S:LAP1_3) were combined to make the ‘Control’ and ‘LAP1’ replicate sets, respectively (Fig. 2a; Additional file 12: Figure S4). Differential expression analysis between the two replicate sets was performed on raw counts with annotated mRNAs via the DESeq2 Filter (within SeqMonk) . This output was combined with the Intensity Difference Filter (within SeqMonk), a statistical based fold change filter that works by constructing a local distribution of differences for mRNAs with similar average intensity to the mRNA being examined. Genes were considered significantly differentially expressed when the adjusted P value was <0.05 for both filters. In other words, the P value represents the likelihood that a given gene is significantly differentially expressed by chance, taking into consideration the absolute expression level and how variable other genes are that are expressed at a similar level. In general, highly expressed genes have low variability and lowly expressed genes have high variability. Hierarchical correlation clustering of the genes statistically differentially expressed between ‘Control’ and ‘LAP1’ was performed with SeqMonk (R-value threshold of 0.7). Differentially expressed genes were subjected to enriched GO categorization using the GO Analysis Toolkit and Agricultural Community (agriGO) database with default settings .
Transcriptome analysis of Nicotiana benthamiana infiltrations
Infiltrated leaf tissue was harvested 3 days post infiltration. Total RNA was extracted from the infiltrated leaf tissue using RNeasy Plant Mini Kit (Qiagen). RNAseq libraries were prepared from the triplicate samples for each type of infiltration, using 500 ng of total RNA, and indexed with the TruSeq mRNA Library Prep Kit according to the protocol (Illumina). The same amount of each RNAseq library was pooled and run on one lane of HiSeq2000 100PE (Illumina).
The quality of the sequencing reads was assessed with FASTQC . RNA sequencing reads were trimmed using DynamicTrim (Phred score ≥20) and filtered on length (≥ 25 bp) using LengthSort . RNAseq reads were aligned against the N. benthamiana transcriptome  using Bowtie2 v2.1.0 , and RSEM v1.2.3  generated raw read counts for each transcript. DESeq  was used to run three differential expression analysis tests: (1) between Buffer infiltrated and Agrobacterium infiltrated; (2) between Agrobacterium and Agrobacterium containing AcMYB10; (3) between Agrobacterium and Agrobacterium containing AtLEC2. Lists of differentially expressed transcripts with a FDR adjusted P value <0.001 are shown in Additional file 7: Table S5, Additional file 8: Table S6, and Additional file 9: Table S7, respectively.
LEGUME ANTHOCYANIN PRODUCTION 1
multidrug and toxic compound extrusion
RCM and RPH designed the experiments; DMB, RHL, RPH, GBG, CMB performed the experiments; DMB, RPH, GBG, CMB, NWA and RCM analysed and interpreted the data; DMB and RCM wrote the manuscript. All authors read and approved the final manuscript.
The authors would like to thank Craig Wood (CSIRO Agriculture) for supplying the 35S:AtLEC2 construct.
The authors declare that they have no competing interests.
Availability of data and materials
The RNA sequencing datasets supporting the conclusions of this article are available in the Sequence Read Archive (SRA) database under study accession SRP091342.
Consent for publication
Ethics approval and consent to participate
This work was supported by a grant from the New Zealand Marsden Fund (UOO1208).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Zhang JZ. Overexpression analysis of plant transcription factors. Curr Opin Plant Biol. 2003;6:430–40.View ArticlePubMedGoogle Scholar
- Prelich G. Gene overexpression: uses, mechanisms, and interpretation. Genetics. 2012;190:841–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Kondou Y, Higuchi M, Matsui M. High-throughput characterization of plant gene functions by using gain-of-function technology. Annu Rev Plant Biol. 2010;61:373–93.View ArticlePubMedGoogle Scholar
- Krenek P, Samajova O, Luptovciak I, Doskocilova A, Komis G, Samaj J. Transient plant transformation mediated by Agrobacterium tumefaciens: principles, methods and applications. Biotechnol Adv. 2015;33:1024–42.View ArticlePubMedGoogle Scholar
- Jones DH, Doherty A, Sparks AC. Transient transformation of plants. Meth Mol Biol 2009;513:131–52.View ArticleGoogle Scholar
- Joensuu JJ, Conley AJ, Lienemann M, Brandle JE, Linder MB, Menassa R. Hydrophobin fusions for high-level transient protein expression and purification in Nicotiana benthamiana. Plant Physiol. 2010;152:622–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Vézina L-P, Faye L, Lerouge P, D’Aoust M-A, Marquet-Blouin E, Burel C, et al. Transient co-expression for fast and high-yield production of antibodies with human-like N-glycans in plants. Plant Biotechnol J. 2009;7:442–55.View ArticlePubMedGoogle Scholar
- Walter M, Chaban C, Schütze K, Batistic O, Weckermann K, Näke C, et al. Visualization of protein interactions in living plant cells using bimolecular fluorescence complementation. Plant J. 2004;40:428–38.View ArticlePubMedGoogle Scholar
- Goodin MM, Dietzgen RG, Schichnes D, Ruzin S, Jackson AO. pGD vectors: versatile tools for the expression of green and red fluorescent protein fusions in agroinfiltrated plant leaves. Plant J. 2002;31:375–83.View ArticlePubMedGoogle Scholar
- Ma L, Lukasik E, Gawehns F, Takken FLW. The use of agroinfiltration for transient expression of plant resistance and fungal effector proteins in Nicotiana benthamiana leaves. In: Bolton DM, Thomma PHJB, editors. Plant fungal pathogens: methods protocols. Totowa: Humana Press; 2012. p. 61–74.View ArticleGoogle Scholar
- Dare AP, Schaffer RJ, Lin-Wang K, Allan AC, Hellens RP. Identification of a cis-regulatory element by transient analysis of co-ordinately regulated genes. Plant Methods. 2008;4:1–10.View ArticleGoogle Scholar
- Hellens RP, Allan AC, Friel EN, Bolitho K, Grafton K, Templeton MD, et al. Transient expression vectors for functional genomics, quantification of promoter activity and RNA silencing in plants. Plant Methods. 2005;1:13.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Y, Li R, Qi M. In vivo analysis of plant promoters and transcription factors by agroinfiltration of tobacco leaves. Plant J. 2000;22:543–51.View ArticlePubMedGoogle Scholar
- Grimberg Å, Carlsson AS, Marttila S, Bhalerao R, Hofvander P. Transcriptional transitions in Nicotiana benthamiana leaves upon induction of oil synthesis by WRINKLED1 homologs from diverse species and tissues. BMC Plant Biol. 2015;15:1–17.View ArticleGoogle Scholar
- Chung BYW, Simons C, Firth AE, Brown CM, Hellens RP. Effect of 5’UTR introns on gene expression in Arabidopsis thaliana. BMC Genom. 2006;7:1–13.View ArticleGoogle Scholar
- Simpson GG, Laurie RE, Dijkwel PP, Quesada V, Stockwell PA, Dean C, et al. Noncanonical translation initiation of the Arabidopsis flowering time and alternative polyadenylation regulator FCA. Plant Cell. 2010;22:3764–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Laing WA, Martinez-Sanchez M, Wright MA, Bulley SM, Brewster D, Dare AP, et al. An upstream open reading frame is essential for feedback regulation of ascorbate biosynthesis in Arabidopsis. Plant Cell. 2015;27:772–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Wroblewski T, Tomczak A, Michelmore R. Optimization of Agrobacterium-mediated transient assays of gene expression in lettuce, tomato and Arabidopsis. Plant Biotechnol J. 2005;3:259–73.View ArticlePubMedGoogle Scholar
- Bhaskar PB, Venkateshwaran M, Wu L, Ane J-M, Jiang J. Agrobacterium-mediated transient gene expression and silencing: a rapid tool for functional gene assay in potato. PLoS ONE. 2009;4:e5812.View ArticlePubMedPubMed CentralGoogle Scholar
- Santos-Rosa M, Poutaraud A, Merdinoglu D, Mestre P. Development of a transient expression system in grapevine via agro-infiltration. Plant Cell Rep. 2008;27:1053–63.View ArticlePubMedGoogle Scholar
- Figueiredo JFL, Römer P, Lahaye T, Graham JH, White FF, Jones JB. Agrobacterium-mediated transient expression in citrus leaves: a rapid tool for gene expression and functional gene assay. Plant Cell Rep. 2011;30:1339–45.View ArticlePubMedGoogle Scholar
- Faizal A, Geelen D. Agroinfiltration of intact leaves as a method for the transient and stable transformation of saponin producing Maesa lanceolata. Plant Cell Rep. 2012;31:1517–26.View ArticlePubMedGoogle Scholar
- Chen Q, Dent M, Hurtado J, Stahnke J, McNulty A, Leuzinger K, et al. Transient protein expression by agroinfiltration in lettuce. In: MacDonald J, Kolotilin I, Menassa R, editors. Recombinant proteins from plants methods Protocols New York. NY: Springer; 2016. p. 55–67.View ArticleGoogle Scholar
- Conner AJ, Albert NW, Deroles SC. Transformation and regeneration of Petunia. In: Gerats T, Strommer J, editors. Petunia evolutionary, developmental and physiological genetics. New York: Springer; 2009. p. 395–409.Google Scholar
- Picard K, Lee R, Hellens R, Macknight R. Transient gene expression in Medicago truncatula leaves via agroinfiltration. In: Rose JR, editor. Legume genomics methods protocols. Totowa: Humana Press; 2013. p. 215–26.View ArticleGoogle Scholar
- Peel GJ, Pang Y, Modolo LV, Dixon RA. The LAP1 MYB transcription factor orchestrates anthocyanidin biosynthesis and glycosylation in Medicago. Plant J. 2009;59(1):136–49.View ArticlePubMedGoogle Scholar
- Voinnet O, Rivas S, Mestre P, Baulcombe D. An enhanced transient expression system in plants based on suppression of gene silencing by the p19 protein of tomato bushy stunt virus. Plant J. 2003;33:949–56.View ArticlePubMedGoogle Scholar
- Li P, Dong Q, Ge S, He X, Verdier J, Li D, et al. Metabolic engineering of proanthocyanidin production by repressing the isoflavone pathways and redirecting anthocyanidin precursor flux in legume. Plant Biotechnol J. 2016. doi:10.1111/pbi.12524.PubMed CentralGoogle Scholar
- Li Y-G, Tanner G, Larkin P. The DMACA–HCl protocol and the threshold proanthocyanidin content for bloat safety in forage legumes. J Sci Food Agric. 1996;70:89–101.View ArticleGoogle Scholar
- Albert NW, Griffiths AG, Cousins GR, Verry IM, Williams WM. Anthocyanin leaf markings are regulated by a family of R2R3-MYB genes in the genus Trifolium. New Phytol. 2015;205:882–93.View ArticlePubMedGoogle Scholar
- Hancock KR, Collette V, Fraser K, Greig M, Xue H, Richardson K, et al. Expression of the R2R3-MYB transcription factor TaMYB14 from Trifolium arvense activates proanthocyanidin biosynthesis in the legumes Trifolium repens and Medicago sativa. Plant Physiol. 2012;159:1204–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y, Butelli E, Martin C. Engineering anthocyanin biosynthesis in plants. Curr Opin Plant Biol. 2014;19:81–90.View ArticlePubMedGoogle Scholar
- Zhao J, Dixon RA. MATE transporters facilitate vacuolar uptake of epicatechin 3′–O–glucoside for proanthocyanidin biosynthesis in Medicago truncatula and Arabidopsis. Plant Cell. 2009;21:2323–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao J, Huhman D, Shadle G, He X-Z, Sumner LW, Tang Y, et al. MATE2 mediates vacuolar sequestration of flavonoid glycosides and glycoside malonates in Medicago truncatula. Plant Cell. 2011;23:1536–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):1.View ArticleGoogle Scholar
- Koes R, Verweij W, Quattrocchio F. Flavonoids: a colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci. 2005;10:236–42.View ArticlePubMedGoogle Scholar
- Albert NW, Davies KM, Lewis DH, Zhang H, Montefiori M, Brendolise C, et al. A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots. Plant Cell. 2014;26:962–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Li P, Chen B, Zhang G, Chen L, Dong Q, Wen J, et al. Regulation of anthocyanin and proanthocyanidin biosynthesis by Medicago truncatula bHLH transcription factor MtTT8. New Phytol. 2016;210:905–21.View ArticlePubMedGoogle Scholar
- Albert NW, Davies KM, Schwinn KE. Repression: the dark side of anthocyanin regulation? Acta Hortic. 2014;1048:129–36.View ArticleGoogle Scholar
- Albert NW. Subspecialization of R2R3-MYB repressors for anthocyanin and proanthocyanidin regulation in forage legumes. Front Plant Sci. 2015;6:1–13.View ArticleGoogle Scholar
- Jun JH, Liu C, Xiao X, Dixon RA. The transcriptional repressor MYB2 regulates both spatial and temporal patterns of proanthocyandin and anthocyanin pigmentation in Medicago truncatula. Plant Cell. 2015;27(10):2860–79.PubMedPubMed CentralGoogle Scholar
- Zimmermann IM, Heim MA, Weisshaar B, Uhrig JF. Comprehensive identification of Arabidopsis thaliana MYB transcription factors interacting with R/B-like BHLH proteins. Plant J. 2004;40:22–34.View ArticlePubMedGoogle Scholar
- Albert NW, Lewis DH, Zhang H, Irving LJ, Jameson PE, Davies KM. Light-induced vegetative anthocyanin pigmentation in Petunia. J Exp Bot. 2009;60:2191–202.View ArticlePubMedPubMed CentralGoogle Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z. agriGO: A GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38:64–70.View ArticleGoogle Scholar
- Bhargava A, Mansfield SD, Hall HC, Douglas CJ, Ellis BE. MYB75 functions in regulation of secondary cell wall formation in the Arabidopsis inflorescence stem. Plant Physiol. 2010;154:1428–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Bhargava A, Ahad A, Wang S, Mansfield SD, Haughn GW, Douglas CJ, et al. The interacting MYB75 and KNAT7 transcription factors modulate secondary cell wall deposition both in stems and seed coat in Arabidopsis. Planta. 2013;237:1199–211.View ArticlePubMedGoogle Scholar
- Fei H, Ellis BE, Vessey JK. Carbon partitioning in tissues of a gain-of-function mutant (MYB75/PAP1-D) and a loss-of-function mutant (myb75-1) in Arabidopsis thaliana. Botany. 2013;92:93–9.View ArticleGoogle Scholar
- Shamloul M, Trusa J, Mett V, Yusibov V. Optimization and utilization of Agrobacterium-mediated transient protein production in Nicotiana. J Vis Exp. 2014. doi:10.3791/51204.PubMedPubMed CentralGoogle Scholar
- Fraser LG, Seal AG, Montefiori M, McGhie TK, Tsang GK, Datson PM, et al. An R2R3 MYB transcription factor determines red petal colour in an Actinidia(kiwifruit) hybrid population. BMC Genom. 2013;14:1–19.View ArticleGoogle Scholar
- Stone SL, Kwong LW, Yee KM, Pelletier J, Lepiniec L, Fischer RL, et al. LEAFY COTYLEDON2 encodes a B3 domain transcription factor that induces embryo development. Proc Natl Acad Sci. 2001;98:11806–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Stone SL, Braybrook SA, Paula SL, Kwong LW, Meuser J, Pelletier J, et al. Arabidopsis LEAFY COTYLEDON2 induces maturation traits and auxin activity: implications for somatic embryogenesis. Proc Natl Acad Sci. 2008;105:3151–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Davies KM, Albert NW, Schwinn KE. From landing lights to mimicry: the molecular regulation of flower colouration and mechanisms for pigmentation patterning. Funct Plant Biol. 2012;39(8):619–38.View ArticleGoogle Scholar
- Gonzalez A, Zhao M, Leavitt JM, Lloyd AM. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J. 2008;53:814–27.View ArticlePubMedGoogle Scholar
- Schwinn K, Venail J, Shang Y, Mackay S, Alm V, Butelli E, et al. A small family of MYB-regulatory genes controls floral pigmentation intensity and patterning in the genus Antirrhinum. Plant Cell. 2006;18:831–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Albert NW, Davies KM, Schwinn KE. Gene regulation networks generate diverse pigmentation patterns in plants. Plant Signal Behav Biosci. 2014;9:e29526.View ArticleGoogle Scholar
- Wong ACS, Hecht VFG, Picard K, Diwadkar P, Laurie RE, Wen J, et al. Isolation and functional analysis of CONSTANS-LIKE genes suggests that a central role for CONSTANS in flowering time control is not evolutionarily conserved in Medicago truncatula. Front Plant Sci. 2014;5:486.View ArticlePubMedPubMed CentralGoogle Scholar
- Laurie RE, Diwadkar P, Jaudal M, Zhang L, Hecht V, Wen J, et al. The Medicago FLOWERING LOCUS T homolog, MtFTa1, is a key regulator of flowering time. Plant Physiol. 2011;156:2207–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Karimi M, Inzé D, Depicker A. GATEWAY™ vectors for Agrobacterium-mediated plant transformation. Trends Plant Sci. 2002;7:193–5.View ArticlePubMedGoogle Scholar
- Petrie JR, Shrestha P, Liu Q, Mansour MP, Wood CC, Zhou X-R, et al. Rapid expression of transgenes driven by seed-specific constructs in leaf tissue: DHA production. Plant Methods. 2010;6:1–6.View ArticleGoogle Scholar
- Andrews S. FastQC: a quality control tool for high throughput sequence data. [Internet]. [cited 2016 Apr 14]. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
- Stockwell PA, Chatterjee A, Rodger EJ, Morison IM. DMAP: differential methylation analysis package for RRBS and WGBS data. Bioinformatics. 2014;30:1814–22.View ArticlePubMedGoogle Scholar
- Andrews S. SeqMonk [Internet]. [cited 2016 Apr 14]. http://www.bioinformatics.babraham.ac.uk/projects/seqmonk/.
- Cox MP, Peterson DA, Biggs PJ. SolexaQA: at-a-glance quality assessment of Illumina second-generation sequencing data. BMC Bioinform. 2010;11:1–6.View ArticleGoogle Scholar
- Nakasugi K, Crowhurst RN, Bally J, Wood CC, Hellens RP, Waterhouse PM. De novo transcriptome sequence assembly and analysis of RNA silencing genes of Nicotiana benthamiana. PLoS ONE. 2013;8:e59534.View ArticlePubMedPubMed CentralGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011;12:1–16.Google Scholar
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:1–12.View ArticleGoogle Scholar