Transcriptome sequencing of transgenic poplar (Populus × euramericana 'Guariento') expressing multiple resistance genes
- Weixi Zhang†1, 2,
- Yanguang Chu†1, 2,
- Changjun Ding1, 2,
- Bingyu Zhang1, 2,
- Qinjun Huang1, 2,
- Zanmin Hu3,
- Rongfeng Huang4,
- Yingchuan Tian5 and
- Xiaohua Su1, 2Email author
© Zhang et al.; licensee BioMed Central Ltd. 2014
Published: 20 June 2014
Transgenic poplar (Populus × euramericana 'Guariento') plants harboring five exogenous, stress-related genes exhibit increased tolerance to multiple stresses including drought, salt, waterlogging, and insect feeding, but the complex mechanisms underlying stress tolerance in these plants have not been elucidated. Here, we analyzed the differences in the transcriptomes of the transgenic poplar line D5-20 and the non-transgenic line D5-0 using high-throughput transcriptome sequencing techniques and elucidated the functions of the differentially expressed genes using various functional annotation methods.
We generated 11.80 Gb of sequencing data containing 63, 430, 901 sequences, with an average length of 200 bp. The processed sequences were mapped to reference genome sequences of Populus trichocarpa. An average of 62.30% and 61.48% sequences could be aligned with the reference genomes for D5-20 and D5-0, respectively. We detected 11,352 (D5-20) and 11,372 expressed genes (D5-0), 7,624 (56.61%; D5-20) and 7,453 (65.54%; D5-0) of which could be functionally annotated. A total of 782 differentially expressed genes in D5-20 were identified compared with D5-0, including 628 up-regulated and 154 down-regulated genes. In addition, 196 genes with putative functions related to stress responses were also annotated. Gene Ontology (GO) analysis revealed that 346 differentially expressed genes are mainly involved in 67 biological functions, such as DNA binding and nucleus. KEGG annotation revealed that 36 genes (21 up-regulated and 15 down-regulated) were enriched in 51 biological pathways, 9 of which are linked to glucose metabolism. KOG functional classification revealed that 475 genes were enriched in 23 types of KOG functions.
These results suggest that the transferred exogenous genes altered the expression of stress (biotic and abiotic) response genes, which were distributed in different metabolic pathways and were linked to some extent. Our results provide a theoretic basis for investigating the functional mechanisms of exogenous genes in transgenic plants.
Forest trees are renewable natural resources that are vital to the balance of the terrestrial ecosystem and have important commercial applications, including timber wood, paper and pulp, and biofuel production. The growth and development of forest trees are frequently challenged by biotic (such as pests and diseases) and abiotic (such as drought, soil salinity, and flooding) stresses, although natural forests have evolved a certain ability to cope with these adverse environmental factors. Genetically engineered plants have been developed in a wide range of tree species, and numerous transgenic clones with improved traits have been generated, some of which have undergone field trials to the environmental release stage; the genetic stability of genetically modified forests has been verified in field trials . Forest molecular breeding research has been greatly enhanced by the completion of genome-wide sequencing projects, the development of genome-wide chips, and the establishment of genetic maps in poplar.
Genes of interest that can be used in the genetic engineering of abiotic stress-resistant plants can be divided into two categories: 1) genes that directly respond to stress, which mainly include genes encoding functional proteins that protect the cell against stress damage, such as enzymes for the synthesis of osmolytes, and enzymes for the removal of active oxygen; 2) genes that regulate genetic expression and signal transduction under stress, which mainly include genes encoding transcription factors, protein kinases, and others. For example, some genes are involved in the expression of osmolytes, such as proline, glycine, betaine, sucrose, fructan, and so on, and they help plants accumulate osmolytes to maintain the osmotic balance and body moisture levels, which improves the drought tolerance of plants under drought and salt stress . Indeed, the exogenous expression of the fructan synthase gene (sacB) and trehalose-6-phosphate synthase genes (otsB, otsA, and TPS1) can improve the drought-resistance, salt-tolerance, and low-temperature-resistance of transgenic poplar , rice , spinach beet , and tobacco [6, 7]. Transporter genes can help alter the ionic and osmotic balance in plants by up-regulating the expression levels of proteins to control ion transport under abiotic stress ; for example, HAL1 expression improves the salt tolerance of transgenic tomatoes . Vitreoscilla hemoglobin (VHb), which has obligate aerobic properties, can be synthesized at substantial levels in plants under hypoxic conditions ; VHb overexpression can improve cellular oxygen ion levels and terminal oxidase activity under waterlogging stress. VHb gene expression can facilitate the growth of transgenic Nicotiana tabacum [11, 12], Datura innoxia , and Petunia hybrida Vilm .
Transcriptional regulation is a major stress response mechanism in plants. Under stress conditions, transcription factors involved in stress resistance can regulate the simultaneous expression of multiple stress-resistance genes and genes involved in the transport of stress signals . Therefore, regulating the expression of transcription factors has been proposed as a way to improve the stress resistance of plants. Numerous studies have shown that members of the MYB, MYC, ERF, bZIP, and WRKY transcription factor families are involved in stress response regulation . For example, overexpression of CaPF1 , AtDREB1A [17, 18], CBF3/DREB1A , SNAC1 , and others can increase the stress resistance of transgenic plants. Protein kinase plays an important role in responses to environmental changes and signal transduction in plants. For example, the expression of calcium-dependent and calmodulin-independent protein kinase (CDPK) genes, including NtCDPK4 , OsCPK6, OsCPK13, OsCPK17, and OsCPK25 , increases under drought, high salt, or low temperature stress, and tolerance to drought and high salt stress is significantly improved in transgenic Arabidopsis thaliana harboring AtCPK23 . Receptor-like kinases (RLKs) on cell membranes can sense external stimuli and are involved in intracellular signal transduction. Overexpression of RLK genes, e.g., CaMRP1  and StLRPK1 , can significantly improve the survival ability of plants under environmental stress conditions. Additionally, studies involving biotic stress resistance genes have shown that Bacillus thuringiensis (Bt) expresses proteins that are toxic to a variety of insects, and these proteins have been broadly applied in anti-insect plant research . Protease inhibitors, such as cysteine proteinase inhibitor (OCI)  and cowpea trypsin inhibitor (CPTI) , have been used independently or in conjunction with Bt genes to significantly increase the insect resistance of important agricultural plants . In addition, overexpression of the chitinase gene can significantly increase disease resistance in plants .
Numerous transgenic plants that are resistant to abiotic stress (drought, waterlogging, salinity) and biotic stress (disease, pests) have been produced using genetic transformation techniques, and extensive research on their exogenous gene expression has been performed [2, 8, 30, 31]. Nevertheless, many agronomic traits, such as insect and salt resistance, are usually jointly controlled by different genes. Therefore, some studies have focused on plant agronomic trait improvement and molecular breeding by means of gene connection or co-transformation to produce new, excellent varieties. However, most genetic transformation is limited to 1-3 major gene(s) in the same or related pathways. For example, Ye et al.  co-transformed rice with two major plant enzyme genes involved in regulating the synthesis of previtamin A or lycopene in endosperm chromatophore, producing transgenic plants that synthesize B-carotene in the endosperm, namely, golden rice. Most studies examining how transgenic plants respond to stress have been limited to examining to the expression of exogenous genes and the accumulation of metabolic substances that are directly related to stress responses. Few studies have examined the expression of stress resistance genes in plants on a genome-wide level. For example, Chan et al.  examined the transcriptome of Arabidopsis thaliana harboring the mannose-6-phosphate reductase gene M6PR an found that this gene activates the downstream abscisic acid (ABA) pathway by up-regulating the expression of ABA receptor genes (PYL4, PYL5, and PYL6) and down-regulating the expression of protein phosphatase 2C genes (ABI1 and ABI2) under salt stress. Maria et al. found that half of the differentially expressed genes in the transcriptome of transgenic rice (compared to non-transgenic rice) are linked to transferred, exogenous genes. Coll et al. evaluated the transcriptome differences between the commercial transgenic maize line MON810 and non-transgenic maize (in a 1/3 maize gene expression level) and found few differentially expressed genes between these lines. However, the different expression levels between transcriptomes of transgenic and non-transgenic plants have not been discovered.
Populus × euramericana is an interspecific hybrid produced from the cross of Populus nigra and Populus deltoides. Many P. × euramericana clones have been commercialized and used in forestry production and to promote ecosystem stability. We previously generated transgenic Populus × euramericana 'Guariento' harboring five effect genes . These five genes include the following: Vgb, encoding aerobic Vitreoscilla hemoglobin (VHb); SacB, encoding levansucrase, which is involved in Bacillus subtilis fructan biosynthesis; BtCry3A, encoding the Bt endotoxin from Coleoptera; OC-I, an anti-insect gene, encoding rice cystatin; and JERF36, a tomato gene, encoding the jasmonate (JA)/ethylene (ET) response factor protein. Greenhouse/laboratory experiments and field trials of two transgenic clones (D5-20 and D5-21), including stress (drought, soil salinity, and flooding) resistance and field experiments, have shown that the transgenic plants have increased tolerance to multiple stresses including drought, salt, waterlogging and insect feeding . In the present study, we intensively sequenced the transgenic clone D5-20 and the non-transgenic clone D5-0 using high-throughput transcriptome sequencing techniques. We then examined the expression of differentially expressed genes in transgenic vs. non-transgenic lines at the genome-wide level, and we explored the mechanisms used by the exogenous genes in the transgenic plants.
Illumina sequencing and alignment to the reference genome
Comparative statistics of sequencing data output and reference genomes
Global analysis of gene expression
To obtain robust information about the biological differences among samples, it is important to utilize biological replications. Here, we utilized a reproducible experimental design (Group A: Aa + A2 and so on; two groups in total) to analyze the expression profiles (a total of three alignment groups); variances were calculated at the levels of genes and gene isoforms separately to obtain typical samples for intensive sequencing analysis. The results presented below highlight the differences in gene expression levels between D5-20 and D5-0.
We calculated the Fragments Per Kilobase of transcript per Million fragments mapped (FPKM) values to indicate the gene expression levels within reference gene regions in standard read mode using Cufflink software. The FPKM approach is consistent with the RPKM calculation method, and it can be used to eliminate the influence of gene length and sequencing amount on the calculation of gene expression levels. Calculated gene expression levels can be directly applied when comparing the gene expression differences between species . The expression of mRNA has two notable features: heterogeneity and redundancy. A small portion of mRNA is highly expressed, while most mRNA is expressed at lower levels. Therefore, data concerning the gene expression density distribution within samples can be used to evaluate the normality of RNA-Seq data. Additional file 1 shows that the sequencing data in the present study are distributed normally and can be used for further analysis.
Distribution of gene sequences detected in P. × euramericana 'Guariento' via RNA-Seq
Gene length (bp)
Number of sequences
We examined the gene sequences using the KOG annotation system (http://genome.jgi.doe.gov/ Tutorial/tutorial/kog.html). Gene functional annotation information was provided for expressed P. × euramericana 'Guariento' genes. We found that 7,624 genes (D5-20) and 7,453 genes (D5-0) were annotated in the two specimens, respectively, accounting for 56.61% (11,352) and 65.54 % (11,372) of the total number of genes.
Analysis of differentially expressed genes in transgenic clone compared with nontransgenic poplar
Functional annotation of differentially expressed genes
Biological pathway analysis of differentially expressed genes
We performed pathway enrichment analysis to further investigate the biological functions of differentially expressed genes. We mapped all of these genes to reference pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.ad.jp/kegg/) to identify the biological pathways in which the genes may be involved. Of the 782 differentially expressed genes, 36 genes (21 up-regulated and 15 down-regulated) could be assigned to 51 KEGG pathways (Additional file 4). Notably, 18 genes were assigned to 9 glucose metabolic pathways. Among these genes, 13 (seven up- and six down-regulated), four (two up- and two down-regulated), and two (down-regulated) were predicted to be involved in "starch and sucrose metabolic (KO00500)", "Pentose and glucuronate interconversions (KO00040)", and "glycolysis/gluconeogenesis (KO00010)". Six other down-regulated genes encoding fructose-bisphosphate aldolases were sorted into glucose metabolic-related pathways (KO00010, KO00030, and KO00051). Additionally, 10 genes were linked to seven KEGG pathways related to amino acid metabolism, including two up-regulated genes involved in "cysteine and methionine metabolism" (KO00270). Moreover, we also identified one and two up-regulated genes that could be associated with "glutathione metabolic (KO00480)" and "Drug metabolism-other enzymes (KO00983)", respectively (Additional file 5). These KEGG annotations provide important clues for investigating specific biological processes that can be influenced by the expression of foreign genes in poplar transformed with multiple genes.
Stress response-related genes
We performed homologous protein alignment and annotation analysis for differentially expressed genes using SwissProt and the Uniport nonredundant protein database. The translated protein sequences encoded by the differentially expressed genes were aligned using Blastp analysis of protein databank sequences; the optimally aligned protein with e value < 1E-5 were used to identify the candidate names of differentially expressed genes/proteins. A total of 730 significantly differentially expressed genes could be annotated to homologous proteins, and the candidate names of 492 differentially expressed genes/proteins were identified. Among the 492 genes, 196 were considered to be putative stress response genes based on GO and KOG analysis. These stress-related genes could be sorted into seven major classes (Additional file 6); of these, serine/threonine protein kinase (39 genes), signaling molecule (27), and transcription factor (34) were the three classes with largest number of genes. Other classes including genes involved in oxygen metabolism, transporters, molecular chaperones, AAA+-type ATPases, and chitinases were also identified. Interestingly, members of several gene families were overrepresented among the stress-related genes, such as 17 up-regulated genes encoding receptor-like kinases, 11 encoding ubiquitin-protein ligases (nine up- and two down-regulated), 14 encoding AP2/ERF transcription factors (11 up- and three down-regulated), 11 encoding cytochrome P450 enzymes (four up- and seven down-regulated), and nine up-regulated genes encoding esterase/lipases.
Phylogenetic analysis of stress-related gene
Validation of RNA-Seq results by qRT-RCR
We performed gene expression analysis of transgenic poplar using Illumina HiSeq 2000 paired-end sequencing (RNA-Seq). We have generated 63, 430, 901 sequence reads (200 bp in length) corresponding to 11.80 Gb of sequencing data. Approximately 73.33% of the sequences could be mapped onto the reference genome sequence of P. trichocarpa. In grape, a total of 59,372,544 sequences reads corresponding to 2.2 Gb of sequencing data were generated by RNA-seq sequencing from three developmental stages, and approximately 82.4% of the sequences were aligned with the reference genome (Pinot Noir 40024) . This result suggests that significant genetic divergence exists in the genus Populus, and the degree of genomic difference in poplars may be higher than that in grape. Global analysis of gene expression revealed transcriptome reprogramming in transgenic poplar (D5-20) compared with the non-transgenic clone D5-0. A total of 782 genes were found to be differentially expressed, including 628 up-regulated (80.3%) and 154 down-regulated genes (19.7%), indicating that transcriptional activators play a major role of in the insertion/expression of the five transgenes in D5-20. Moreover, 197 putative biotic or abiotic stress-responsive genes were also identified. The data from RNA-Seq were validated by qRT-PCR analysis of 12 genes. The expression of nine out of 12 genes showed diverse trends under stress (drought, salt, or flooding) treatment, as revealed by qRT-PCR analysis (Figure 8), which may be associated with the introduction of exogenous genes and may involve different regulatory networks of these genes. Compared with non-transgenic poplar, the genes exhibited similar trends in expression under non-stressed conditions and under at least one stress condition, as revealed by qRT-PCR and RNA-Seq analysis. These genes may be involved in different regulatory networks under different stress conditions. These complex results suggested that there are complex interactions among multiple genes and the exogenously expressed genes. The magnitude of the number differentially expressed genes found in this study is similar to previously reported ranges for transgenic plants harboring single transformation constructs [34, 41, 42].
Fructans are involved in osmoregulation. Fructans accumulate under stress conditions as a result of the fructosyltransferase effect. The SacB in Bacillus subtilis encodes secreted levansucrase, which catalyzes the synthesis of fructans from substrate sucrose. The expression of sacB can increase plant cold and drought tolerance by significantly enhancing the accumulation of fructans, as observed in sugar beet, tobacco, corn, and other transgenic plants [2, 43]. In plants, sucrose synthase (EC 220.127.116.11, SuSy) catalyzes the decomposition of sucrose into fructose and glucose (sugar + UDP ↔ Fructose + UDPG); such a reaction is reversible [44–46]. Different combinations of fructosyltransferases catalyze the synthesis of different fructans . Fructans can be used as storage carbohydrates, they are capable of regulating the sucrose pool size in photosynthetic tissues and during sucrose metabolism, and so on; they release fructose and reduce the water freezing point in tissues via depolymerization of fructans . Therefore, increasing the fructan content is thought to influence the contents of other soluble sugars, total soluble sugars, and sugar metabolic products. For example, the contents of non-structural carbohydrates (glucose, fructose, sucrose, starch, fructan) are significantly higher in trans-sacB potatoes than in the wild type . We previously found that the fructan accumulation levels increased in transgenic poplars (D5-20) under drought stress . Transcriptome analysis of transgenic poplar (D5-20) revealed that the differentially expressed genes are enriched in 11 biological pathways linked to glucose metabolism, primarily including starch and sucrose metabolism pathway (KO00500), glucose metabolism pathway (KO00040), and glycolysis/gluconeogenesis pathway (KO00010). The enriched expression of these genes may result in elevated levels of fructan in D5-20. Furthermore, polysaccharide bases in fructan can become inserted between phospholipid bilayer molecules in the cell membrane and can protect lipids and interact with phospholipids to maintain the stability of the cell membrane [50, 51].
Plant glycosyltransferase (GT, EC 2.4.x.y) is a key enzyme for glycosylation that catalyzes the synthesis of specific secondary metabolic products and serves as an important component of secondary metabolism. UDP glycosyltransferases (UGTs) catalyzed the transfer of glycosyl groups from UDP saccharides to a variety of receptor molecules, such as carbohydrates (monosaccharides, oligosaccharides, polysaccharides), non-carbohydrates (proteins, lipids, antibiotics, plant hormones, plant toxins, and others), some exogenous substances (herbicides and pesticides), and so on . Glycosyltransferase are therefore thought to be involved in the tolerance of plants to biotic and abiotic stresses. For example, overexpression of UGT73C5 in transgenic Arabidopsis can improve the resistance against fungal toxins . Moreover, under drought stress, the expression of UGT74E2 (UDP-glucosyltransferase) in Arabidopsis thaliana can improve the rooting capacity and alter anthotaxy traits by regulating IBA and NAA activities, thereby improving resistance against drought and salt stress . In the present study, numerous glycosyltransferase-related genes were differentially expressed, including six UDP-glucuronosyl/UDP-glucosyl transferase genes (four up-regulated and two down-regulated), four galacturonosyl transferase genes (one up-regulated and down-regulated), one UDP-glucose 4-epimerase/UDP-sulfoquinovose synthase gene (up-regulated), four plant glycosyltransferase genes, two glycosyl hydrolase genes, two glucanase genes, two fructose bisphosphate aldolase genes, and so on (down-regulated). According to KOG functional classification, these genes mainly function in energy production and conversion, carbohydrate transport and metabolism, and cell wall/membrane/envelope biogenesis. The phenomena of up-regulation and down-regulation of similar glycosyltransferase genes may be due to the fact that different protein products are encoded by such genes, or that they have different catalytic substrates. For example, phylogenetic analysis of six UDP-glucuronosyl/UDP glycosyltransferase genes revealed that the amino acid sequences in a single group were identical in two up-regulated genes (Potri.017G052400 and Potri.017G052000), and the other two up-regulated genes were in a group as well. These differences may generate different catalytic activities or different catalytic substrates in the enzyme. Studies of the poplar transcriptome under salt stress have shown some glucose glucosyltransferase genes expressed, suggesting that such genes are involved in the salt stress resistance response in plants .
The ERF (ethylene-responsive factor) transcription factors belong to the plant AP2/EREBP (APETALA2/ethylene-responsive element banding protein) superfamily and play important roles in regulating growth, development, and stress response. For example, exogenous JERF expression activates the expression of a large number of downstream genes that function in stress resistance and helps increase resistance against salt, drought, and low temperature stress by activating multiple stress-related cis-acting elements in transgenic tobacco plants[57, 58]. We previously demonstrated that exogenous JERF36 significantly improves the salt tolerance of transgenic poplar (P. alba × Populus berolinensis) . Previous physiological greenhouse tests have shown that overexpression of JERF36 in D5-20 under salt stress can regulate instantaneous water use efficiency (iWUE) and root growth, thus improving salt resistance . In the present study, 14 AP2/ERF transcription factors in D5-20 were found to be differentially expressed, including 11 that were up-regulated and three that were down-regulated. Among these, seven genes were identified as members of the ERF subfamily (six in group B1 and one in group B3) according to genome-wide analysis of the poplar AP2/ERF superfamily . In addition, three genes were in the dehydration-responsive element-bonding protein (DREB) subfamily, and one up-regulated gene was in the AP2 subfamily (Group A2), one was in the DREB subfamily (Group A4), and two were down-regulated and in the AP2 subfamily. These findings suggest that the functions of ERF subfamily members are highly conserved, while the functions of genes in the DREB and AP2 subfamilies may be highly divergent. In addition, the differential expression of these AP2/ERF genes may have resulted from the expression of the transferred exogenous gene JERF36. Indeed, ERF subfamily genes regulate the activities of downstream genes as well as stress signaling molecules, such as ethylene (ET), jasmonate (JA), and salicylic acid (SA). For example, expression profiling has demonstrated that nine ERF subfamily members are involved in the abscisic acid (ABA), SA, JA, and ET signal transduction pathways, and they are also involved in biotic and abiotic stress responses . Other widely recognized transcription factors involved in plant responses to biotic and abiotic stresses , such as WRKY, basic leucine zipper (bZIP) protein, and MYC/MYB, were also identified in our dataset. These factors may participate in ABA and other signal transduction pathways, bind with DNA elements, and induce downstream gene expression in response to stress. Moreover, protein kinases, which act as regulatory factors in plant stress responses, such as receptor-like kinase (RLK), calcium-dependent and calmodulin-independent protein kinase (CDPK), and mitogen-activated protein kinase (MAPK), participate in hormonal signaling transduction processes under various stress conditions. Plant protein kinases can allow signals to be amplified in a step-wise manner through phosphorylation, and they then activate the downstream transcription factors and further induce the expression of resistance genes after they are induced by intracellular second messenger signaling molecules in the signal transduction process. For example, wheat TaMAPK1 can improve TaERF1 activity, promote binding with GCC-box and DRE/CRT elements, and further induce the expression of genes of interest [61, 62]. The JA-regulated MAPKK kinase JAM1 and Arabidopsis thaliana PKS33 kinase participate in the phosphorylation of the tobacco ERF protein ORC1, and Arabidopsis ERF7 further increases ORC1 and ERF7 protein activities . Moreover, the bZIP transcription factor SnRK pathway can be phosphorylated to regulate downstream gene expression . Moreover, recent evidence has shown that protein kinase genes are involved in various biotic and abiotic stress response pathways in an interdependent manner. For example, the CDPK pathway is thought to cross with the MAPK and SnRK pathways [65, 66]. Ca2+- or Ca2+ receptor-calmodulins (CaMs; such as CaM-like Ca-binding proteins with the EF-hand motif, and Ca2+-regulated protein kinases and Ca-binding proteins without the EF-hand motif) play important roles in the signal transduction process in these biotic and abiotic response pathways . In the present study, we detected 20 stress-resistance transcription factors in addition to ERF transcription factors, including five WRKY transcription factors, three zinc-finger transcription factor, three bZip transcription factors, and one MYB transcription factor that were significantly up-regulated. We also detected three MYB transcription factors, one WRKY transcription factor, and two bHLH transcription factor that were significantly down-regulated. We also found that 34 genes encoding serine/threonine protein kinases were up-regulated in the transgenic clones (D5-20), including 17 receptor-like kinases with leucine-rich replicated region, two MAPKs, and two lectin-like kinase; five other serine/threonine protein kinases in D5-20 were down-regulated. A number of genes encoding signaling molecules, including 11 calmodulin binding proteins (nine up- and two down-regulated), 11 ubiquitin ligases (nine up- and two down-regulated), and two up-regulated serine/threonine protein phosphatases, two down-regulated JA amino synthases, and one down-regulated ubiquitin protein were also identified in D5-20 (Additional file 6). These findings suggest that significant, multiple dimensions of the networks involving transcriptional regulation and signal transduction were altered in the D5-20 transcriptome under normal growth condition. Such a phenomenon could largely be attributed to the ectopic expression of transgenes (particularly the transcriptional regulator JERF36) followed by the modulation of the activities of numerous genes via transcriptional regulation, phosphorylation, or ubiquitination processes. Nevertheless, the specific mechanisms underpinning this altered gene expression remain unclear and need to be determined.
Vitreoscilla hemoglobin (VHb) is involved in oxygen-related metabolic pathways in the oxygenation state, such as the respiratory chain. This molecule transfers oxygen to the respiratory chain, regulates the activity of terminal oxidase, alters the efficiency of oxidative phosphorylation and metabolic pathways under hypoxic conditions, and promotes cell growth and the expression of some genes . The expression of the vgb gene can increase the effective oxygen concentration in a plant and enhance chlorophyll biosynthesis [69, 70] and can also promote plant growth [71, 72], improve waterlogging resistance , and increase the detoxification capacity of plants . We detected 21 differentially expressed genes encoding oxygen metabolic pathway enzymes, including nine (five glutaredoxin-related proteins, two thioredoxins, and two multicopper oxidases) that were up-regulated and 12 (three carbonic anhydrases, two peroxidases, two ferric reduction oxidases, three oxidases, and two iron/ascorbate family oxidoreductases) that were down-regulated. We previously reported that the Fv/Fm ratio and chlorophyll content in D5-20 were higher than those of D5-0 under flooding stress, suggesting that the expression of vgb can help maintain the photosynthetic system stability of a plant through the collection and transport of oxygen, thereby further improving flooding resistance. The large number of differentially expressed oxygen metabolic enzyme genes revealed in the present study may therefore be related to the expression of vgb in transgenic poplar.
RNA-Seq of D5-20 also revealed some genes with multiple functions, such as those encoding molecular chaperones, cytochrome P450 transporters, AAA +-ATPase, and chitinase. These genes play essential roles in plant growth and development. For example, DnaJ-like protein is involved in many metabolic processes, functioning as both a molecular chaperone and a regulatory protein . Cytochrome P450 is a large plant protein family that is involved in an array of catalytic reactions, including the synthesis and metabolism of macromolecules (lignin, keratin, and suberin), hormones, signaling molecules, natural pigments, and defensive substances . Transporter proteins are transmembrane proteins that are involved in the uptake, transport, and isolation of ions and small molecules. For example, the uptake and transport of phosphorus in plants relies on phosphorus transporters in the cell membranes of roots . The ATP-binding Cassette (ABC) family proteins play important roles in transmembrane transport of hormones, lipids, carbohydrates, inorganic acids, heavy metals, secondary metabolites, and xenobiotics  and are related to chlorophyll synthesis, Fe/S complex formation, and stomatal movement . An AAA+-ATPase gene (At1g64110) in Arabidopsis thaliana is involved in the abiotic stress response. Chitinase is one of main pathogenesis-related proteins in plants, and it can inhibit the growth of fungi and is involved in the interaction between the host plant and fungal pathogens. Resistance to insects also occurs in plants expressing chitinase genes .
BtCry3A encodes an endotoxin from Bt in Coleoptera and is a type of insecticidal crystal protein (ICP); ICPs are major, active insecticidal ingredients that cannot be directly synthesized in plant cells. Transcriptome analysis has revealed two major biological processes (stress/defense responses and amino acids metabolism) that are influenced by differentially expressed gene in rice expressing a synthetic Cry1b gene . However, such phenomena were not observed in the current study of poplar, which may be particularly associated with the potentially complex interactions among genes influenced by multiple transgenes in transgenic D5-20. Additionally, OC-I encodes a rice cysteine proteinase inhibitor, and its expression can be induced by JA, ABA, and wounding . The OC-I polypeptide may activate the lipid-based signal transduction pathway, by which plant cells release and convert linolenic acid into the oxylipin signaling molecule JA . Wounding can induce intracellular cascades and plasma membrane depolarization, which opens ion channels, increases intracellular Ca2+ concentrations, and activates MAP kinase and phospholipase A. These processes promote the release of linolenic acid in the cell membrane, transfer linolenic acid into the oxylipin signaling molecule JA, and activate numerous defensive genes . In the present study, we detected a number of esterase genes (up-regulated), JA amino synthase genes (down-regulated), and oligopeptide transporter genes (down-regulated), all of which are possibly involved in the abovementioned processes. We previously demonstrated that the expression of BtCry3A and OC-I leads to improved insect resistance (mainly to Plagiodera versicolora) under both laboratory and field conditions [37, 85, 86]. However, the effects of the insertion/expression of BtCry3A and OC-I on the transcriptome of D5-20 remain to be investigated.
Plants face a variety of adverse environmental stresses in the natural growth environment, including biotic stresses (such as fungi, bacteria, viruses, pests, and others) and abiotic stresses (such as high salt, drought, waterlogging, heat, cold, injury, and others). Plants have evolved a systematic defensive mechanism that includes multiple independent and crossed signaling pathways to overcome these adverse environmental conditions. These crossed signaling pathways in plants function during biotic , abiotic , and combined stress resistance responses . Therefore, the majority of plant traits, such as drought tolerance, salt tolerance, and insect resistance are controlled by multiple genes. These genes interact via signaling pathways in response to biotic and abiotic stresses. Such a phenomenon was also demonstrated in the present study; for example, the transferred SacB gene increased the content of fructan, an osmolyte, in transgenic clones, and further altered the expression of glucose metabolism genes. Glycoside metabolism is associated with disease resistance in plants. For example, UDP-glycosyltransferase is involved in resistance against biotic and abiotic stresses through the activities of glycosylated hormones and secondary metabolites in plants [90, 91]. Tomato JERFs play an important role in regulatory networks that integrate ET and ABA signaling pathways; JERFs improve salt, drought, and cold tolerance and disease resistance in transgenic plants by regulating the expression of stress-related genes . ERFs play important regulatory roles in molecular responses to hormones (ET), pathogens, cold, drought, high salt, and other stressors . In the present study, the transferred JERF36 gene caused the differential expression of a series of transcription factors, regulated gene expression, and led to hormonal changes. There is crosstalk between the ABA-, SA-, JA-, and ET-dependent signaling pathways. These signaling molecules achieve precise regulation of defensive responses in plants through synergistic or antagonistic effects . The transferred vgb gene can cause changes in oxygen metabolism in transgenic clones, induce the differential expression of antioxidant protective genes, increase the potential of oxidation stress resistance in transgenic clones, and further improve drought, salt, and reactive oxygen species (ROS) resistance in transgenic clones. A remarkable feature during the early interaction between plants and bacteria is the generation of ROS. ROS can directly help the plant resist pathogenic microorganisms and induce the activation of defensive genes. Therefore, changes in oxygen metabolism can also alter defense responses in plants.
The stress response in plants is a complicated process associated with multiple genes, multiple signal transduction pathways, and multiple gene expression products. Such processes involve the perception and transduction of stress signals, stress signal identification and transduction in appropriate receptors, and the expression of stress resistance genes. Therefore, comprehensively modify the traits of forest trees requires the transfer of genes in dissimilar metabolic pathways and the synergy of multiple genes. In the present study, transcriptome sequencing of transgenic poplar at the genome-wide level revealed that the transferred exogenous genes caused differential expression of stress (biotic and abiotic) response genes; these genes are distributed in different and associated metabolic pathways. The results of this study provide a theoretical basis for investigating the effects of exogenous genes in transgenic poplar lines. We also detected the differential expression of genes for some unexpected traits, such as genes involved in improved wood properties, nutrient utilization efficiency, and others. However, these results require further observation and verification.
The transcriptomes of transgenic (D5-20) and nontransgenic clones (D5-0) of P. × euramericana 'Guariento'  were compared in the current study. The transfected exogenous genes included the following: (1) the vgb gene, encoding VHb and related to the improvement of plant adaptive ability to lean oxygen environments; (2) the SacB gene, encoding levansucrase and related to plant cell osmotic regulation and the improvement of plant drought resistance; (3) the BtCry3A gene, encoding Bt insecticidal crystal protein, which confers highly specific resistance to Coleoptera (including poplar longicorn beetles); (4) the OC-I gene, encoding rice cysteine protease inhibitors that inhibit the growth of most Coleoptera pests; and (5) the JERF36 gene, encoding AP2/EREBP plant transcription factors, which are related to plant stress resistance . The BtCry3A and OC-I genes formed bivalent genes in the same vector. The neomycin phosphotransferase II gene (NPT II) derived from E. coli transposon Tn5 was used as a marker, which provided the plants with kanamycin resistance.
Preparation of cDNA library for RNA-Seq
Cuttings of transgenic (D5-20) and nontransgenic clones (D5-0) grown under the same conditions were selected and cultured in the greenhouse. When the plants grew to approximately 20 cm, with 7-9 leaves, the fourth fully expanded leaves were collected, rapidly frozen in liquid nitrogen, and used for extraction of RNA. Total RNA was extracted using Trizol Reagent (Invitrogen) according to the manufacturer's instructions, and contaminating DNA was digested using DNase I (Promega). The integrity, purity, and quality of total RNA were determined using 1.2% agarose gel electrophoresis and an RNA 6000 Kit (Agilent), and the mRNA was purified using a Dynabeads mRNA Purification Kit (Invitrogen) according to the manufacturer's instructions. The mRNA was randomly broken into fragments after the addition of Fragmentation Buffer (LC-Bio); the mRNA was incubated in 200 μL thin-walled PCR tubes at 94°C for 5 min. The RNA fragments were reverse transcribed into first-strand cDNA, then second-strand cDNA, using random primers and transcriptase. The double-stranded DNA fragments were purified using a QIAquick PCR Purification Kit (Qiagen, Hilden, Germany). The resultant double-stranded cDNA fragments were end-repaired, and an 'A' base was added to their at 3' ends; specific sequencing adaptors were connected to both ends of the DNA fragments. Each DNA fragment was sequenced at a specific joint; 300 bp DNA fragments were recovered through electrophoresis for PCR amplification to enrich the sequencing specimens. The PCR products were purified using a QIAquick Gel Purification Kit (Qiagen). The above procedures were performed according to the instructions for the mRNA-Seq 8-Sample Prep Kit (Illumina).
Sequencing of transcriptome and processing of sequence data
Libraries of both clones (D5-0 and D5-20) were sequenced using a HiSeq 2000 platform (Illumina) under double-end 100 bp sequencing mode. The quality-reliable sequencing fragments were selected, the poor-quality bases were dynamically removed from the 3' end, and the base fragments with high sequencing quality were ultimately reserved after the primer and adaptor sequences were removed from the raw data and the 3' end quality of sequencing fragments was tested. The remaining sequence fragments were used for subsequent analysis in the RNA-Seq module.
Mapping reads to the reference genome and annotation
The pretreated sequences were aligned with the poplar reference genomes (http://genome.jgi-psf.org/Poptr1/Poptr1.home.html/). RNA structural features need to be considered with multiple alignments statistical modes during sequence alignment. The alignment parameters and conditions were as follows: (1) Allowable to 2-base unalignment mode; (2) Allowable to maximally optimal 20-base alignment record for each read sequence; (3) Considerable to alternative spicing condition, set the segment length equal to 1/2 read length; (4) Allowable to at most 1 bp unaligned bases in a segment; (5) Set 3 bp insertion and deletion length to the maximum; (6) Allowable to 0-base unalignment in alternative spicing location, namely, complete alignment; (7) Set the min. isoform fraction to 0.15. The above setting indicates that exon A is presumably D higher than that of exon B in sequencing depth if one junction across two exons is overlapped by S sequencing fragments. Exon junctions are not believed to occur currently if the S/D ratio is less than the set minimum isoform fraction (here, 0.15); (8) Set the minimum intron length to 50 bp; (9) Set the maximum intron length to 50,000 bp; and (10) Allowable to 40-base alignment record and 2 bp unaligned bases under junction probing conditions.
Normalization of gene expression levels
RNA expression levels were calculated using an FPKM indicator. FPKM signifies the fragments per kilo base of transcript per million fragments mapped. The expression levels were calculated using multiple standard methods. In our study, the expression levels and differential gene expression were calculated using standard read mode within reference gene regions using Cufflink software .
Analysis of differentially expressed genes
The expression profiles between repeated experiments (Group A, including four groups, such as Aa + A2) were analyzed in comparison to their differences with three control groups. In difference analysis, the calculation was performed for gene and gene isoform levels. To test the significant variance in expression levels, the log ratios were compared between experimental specimens (conditions) with the log value of a conditional expression level. Fold-changed values between samples were calculated according to FPKM value. The levels of "| log2 Ratio | ≥ 1 and false discovery rate (FDR) ≤ 0.05" were used as the threshold to assess the significance of differential gene expression.
GO functional enrichment analysis was carried out by annotating differentially expressed genes using the Gene Ontology (GO) database and the poplar genome annotation database (http://genome.jgi-psf.org/programs/plants/index.jsf). KOG classifications of differentially expressed genes were performed according to polar genomics annotated information (http://www.phytozome.net/poplar). KEGG pathway analyses of differentially expressed genes were performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.ad.jp/kegg/) database. The enrichment of differentially expressed genes was calculated using Fisher's test and expressed as individual p values.
Annotation and analysis of homologous genes among differentially expressed genes
Homologous protein alignment and annotation analysis were performed for differentially expressed genes using SwissProt and the Uniport non-redundant protein database. The translated protein sequences appropriate to differentially expressed genes were aligned using a Blastp approach and protein databank sequences; the optimally aligned protein sequences with e values < 1E-5 were used to derive candidate names of differentially expressed genes/proteins.
Phylogenetic analysis was performed with the AP2/ERF genes with reference to the evolutionary relationships between poplar AP2/ERF genes; the amino acid sequences of UDP-glucuronosyl/UDP glycosyltransferase were completely aligned using default parameters in MEGA5.1 software. The phylogenetic tree was constructed using the neighbor-joining method.
Quantitative real-time (qRT)-PCR
Cuttings of the transgenic lines D5-20 and the control line D5-0 were cultured in a greenhouse. When the plants grew to approximately 20 cm, with seven to nine leaves, they were exposed to 200 mM NaCl, 20% PEG-6000, or water (water level at 1 cm above the soil surface) for 48 h, and untreated, well-watered plants were used as controls. After each treatment, a mixture of leaves of five plants was collected, frozen immediately in liquid nitrogen, and stored until use. Total RNA was extracted from leaves using an AmbionH Plant RNA Isolation Aid (Applied Biosystems, CA, USA) according to the manufacturer's instructions. Then, cDNA was synthesized using a PrimerScript RT Reagent Kit (TaKaRa, Dalian, China). The qRT-PCR was performed in a LightCycler® 480 II Real-time PCR Instrument (Roche, Swiss) with SYBR Green Realtime PCR Master Mix (TaKaRa). Gene-specific primers were designed to amplify 120-130 bp fragments of foreign genes, and parallel PCR was carried out using a gene-specific primer pair for poplar ACTIN1 (GenBank Accession XM_002298674), which was used as a reference gene. Primer sequences for the real-time PCR assay of the five genes and ACTIN1 are listed in Additional file 7. Five trees were tested per line, and four PCR replicates were performed for each RNA sample. The expression levels of mRNAs were normalized to ACTIN1 and were calculated using the 2-ΔΔCt method .
The publication charges this article were funded by the National High-Tech Research and Development Program of China, 863 Program (Grant No. 2011AA100201) and Twelfth Five National Key Technology R&D program (2012BAD01B03)
This article has been published as part of BMC Genetics Volume 15 Supplement 1, 2014: Selected articles from the International Symposium on Quantitative Genetics and Genomics of Woody Plants. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenet/supplements/15/S1.
- Walter C, Fladung M, Boerjan W: The 20-year environmental safety record of GM trees. Nat Biotechnol. 2010, 28: 656-658. 10.1038/nbt0710-656.View ArticlePubMedGoogle Scholar
- Bartels D, Sunkar R: Drought and salt tolerance in plants. Critical Reviews in Plant Sciences. 2005, 24: 23-58. 10.1080/07352680590910410.View ArticleGoogle Scholar
- Li Y, Su X, Zhang B, Zhang Z: Molecular detection and drought tolerance of SacB-transgenic Populus alba × P. glandulosa. Journal of Beijing Forestry University. 2007, 29: 1-6.Google Scholar
- Garg AK, Kim JK, Owens TG, Ranwala AP, Choi YD, Kochian LV, Wu RJ: Trehalose accumulation in rice plants confers high tolerance levels to different abiotic stresses. Proc Natl Acad Sci USA. 2002, 99: 15898-15903. 10.1073/pnas.252637799.PubMed CentralView ArticlePubMedGoogle Scholar
- Pilon-Smits EA, Hwang S, Mel LC, Zhu Y, Tai JC, Bravo RC, Chen Y, Leustek T, Terry N: Overexpression of ATP sulfurylase in indian mustard leads to increased selenate uptake, reduction, and tolerance. Plant Physiol. 1999, 119: 123-132. 10.1104/pp.119.1.123.PubMed CentralView ArticlePubMedGoogle Scholar
- Konstantinova T, Parvanova D, Atanassov A, Djilianov D: Freezing tolerant tobacco, transformed to accumulate osmoprotectants. Plant Science. 2002, 163: 157-164. 10.1016/S0168-9452(02)00090-0.View ArticleGoogle Scholar
- Almeida AEM, Silva AB, Araujo SS, Cardoso LISA, Santos DM, Torn E JEM, Silva JM, Paul MJ, Fevereiro PS: Responses to water withdrawal of tobacco plants genetically engineered with the AtTPS1 gene: a special reference to photosynthetic parameters. Euphytica. 2007, 154: 113-126. 10.1007/s10681-006-9277-2.View ArticleGoogle Scholar
- Bhatnagar-Mathur P, Vadez V, Sharma KK: Transgenic approaches for abiotic stress tolerance in plants: retrospect and prospects. Plant Cell Rep. 2008, 27: 411-424. 10.1007/s00299-007-0474-9.View ArticlePubMedGoogle Scholar
- Gisbert C, Rus AM, Bolarin MC, Lopez-Coronado JM, Arrillaga I, Montesinos C, Caro M, Serrano R, Moreno V: The yeast HAL1 gene improves salt tolerance of transgenic tomato. Plant Physiol. 2000, 123: 393-402. 10.1104/pp.123.1.393.PubMed CentralView ArticlePubMedGoogle Scholar
- Dikshit KL, Dikshit RP, Webster DA: Study of Vitreoscilla globin (vgb) gene expression and promoter activity in E. coli through transcriptional fusion. Nucleic Acids Res. 1990, 18: 4149-4155. 10.1093/nar/18.14.4149.PubMed CentralView ArticlePubMedGoogle Scholar
- Holmberg N, Lilius G, Bailey JE, Bulow L: Transgenic tobacco expressing Vitreoscilla hemoglobin exhibits enhanced growth and altered metabolite production. Nat Biotechnol. 1997, 15: 244-247. 10.1038/nbt0397-244.View ArticlePubMedGoogle Scholar
- Farrés J, Kallio PT: Improved Cell Growth in Tobacco Suspension Cultures Expressing Vitreoscilla Hemoglobin. Biotechnology Progress. 2002, 18: 229-233. 10.1021/bp010159v.View ArticlePubMedGoogle Scholar
- Bülow L, Holmberg N, Lilius G, Bailey JE: The metabolic effects of native and transgenic hemoglobins on plants. Trends in biotechnology. 1999, 17: 21-24. 10.1016/S0167-7799(98)01252-9.View ArticlePubMedGoogle Scholar
- Mao Z, Hu Y, Zhong J, Wang L, Guo J, Lin Z: Improvement of the hydroponic growth and waterlogging tolerance of petunias by the introduction of vhb gene. ACTA BOTANICA SINICA-CHINESE EDITION-. 2003, 45: 205-210.Google Scholar
- Singh K, Foley RC, Onate-Sanchez L: Transcription factors in plant defense and stress responses. Curr Opin Plant Biol. 2002, 5: 430-436. 10.1016/S1369-5266(02)00289-3.View ArticlePubMedGoogle Scholar
- Tang W, Newton RJ, Li C, Charles TM: Enhanced stress tolerance in transgenic pine expressing the pepper CaPF1 gene is associated with the polyamine biosynthesis. Plant Cell Rep. 2007, 26: 115-124.View ArticlePubMedGoogle Scholar
- Behnam B, Kikuchi A, Celebi-Toprak F, Yamanaka S, Kasuga M, Yamaguchi-Shinozaki K, Watanabe KN: The Arabidopsis DREB1A gene driven by the stress-inducible rd29A promoter increases salt-stress tolerance in proportion to its copy number in tetrasomic tetraploid potato (Solanum tuberosum). Plant biotechnology. 2006, 23: 169-177. 10.5511/plantbiotechnology.23.169.View ArticleGoogle Scholar
- Oh SJ, Song SI, Kim YS, Jang HJ, Kim SY, Kim M, Kim YK, Nahm BH, Kim JK: Arabidopsis CBF3/DREB1A and ABF3 in transgenic rice increased tolerance to abiotic stress without stunting growth. Plant Physiol. 2005, 138: 341-351. 10.1104/pp.104.059147.PubMed CentralView ArticlePubMedGoogle Scholar
- Hu H, Dai M, Yao J, Xiao B, Li X, Zhang Q, Xiong L: Overexpressing a NAM, ATAF, and CUC (NAC) transcription factor enhances drought resistance and salt tolerance in rice. Proc Natl Acad Sci USA. 2006, 103: 12987-12992. 10.1073/pnas.0604882103.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang M, Liang S, Lu YT: Cloning and functional characterization of NtCPK4, a new tobacco calcium-dependent protein kinase. Biochim Biophys Acta. 2005, 1729: 174-185. 10.1016/j.bbaexp.2005.04.006.View ArticlePubMedGoogle Scholar
- Wan B, Lin Y, Mou T: Expression of rice Ca2+-dependent protein kinases (CDPKs) genes under different environmental stresses. FEBS letters. 2007, 581: 1179-1189. 10.1016/j.febslet.2007.02.030.View ArticlePubMedGoogle Scholar
- Ma SY, Wu WH: AtCPK23 functions in Arabidopsis responses to drought and salt stresses. Plant Mol Biol. 2007, 65: 511-518. 10.1007/s11103-007-9187-2.View ArticlePubMedGoogle Scholar
- An SH, Choi HW, Hwang IS, Hong JK, Hwang BK: A novel pepper membrane-located receptor-like protein gene CaMRP1 is required for disease susceptibility, methyl jasmonate insensitivity and salt tolerance. Plant Mol Biol. 2008, 67: 519-533. 10.1007/s11103-008-9337-1.View ArticlePubMedGoogle Scholar
- Wu T, Tian Z, Liu J, Xie C: A novel leucine-rich repeat receptor-like kinase gene in potato, StLRPK1, is involved in response to diverse stresses. Mol Biol Rep. 2009, 36: 2365-2374. 10.1007/s11033-009-9459-9.View ArticlePubMedGoogle Scholar
- Filho OG, Denolf P, Peferoen M, Decazy B, Eskes AB, Frutos R: Susceptibility of the coffee leaf miner (Perileucoptera spp.) to bacillus thuringiensis delta-endotoxins: A model for transgenic perennial crops resistant to endocarpic insects. Curr Microbiol. 1998, 36: 175-179. 10.1007/PL00006763.View ArticlePubMedGoogle Scholar
- Leplé JC, Bonadé-Bottino M, Augustin S, Pilate G, Tan VDL, Delplanque A, Cornu D, Jouanin L: Toxicity toChrysomela tremulae (Coleoptera: Chrysomelidae) of transgenic poplars expressing a cysteine proteinase inhibitor. Molecular Breeding. 1995, 1: 319-328. 10.1007/BF01248409.View ArticleGoogle Scholar
- Zhang Q, Zhang ZY, Lin SZ, Lin YZ: Resistance of transgenic hybrid triploids in Populus tomentosa Carr. against 3 species of lepidopterans following two winter dormancies conferred by high level expression of cowpea trypsin inhibitor gene. Silvae genética. 2005, 54: 108-115.Google Scholar
- Hilder VA, Boulter D: Genetic engineering of crop plants for insect resistance - a critical review. Crop Protection. 1999, 18: 177-191. 10.1016/S0261-2194(99)00028-9.View ArticleGoogle Scholar
- Jia Z, Sun Y, Yuan L, Tian Q, Luo K: The chitinase gene (Bbchit1) from Beauveria bassiana enhances resistance to Cytospora chrysosperma in Populus tomentosa Carr. Biotechnol Lett. 2010, 32: 1325-1332. 10.1007/s10529-010-0297-6.View ArticlePubMedGoogle Scholar
- Sairam RK, Kumutha D, Ezhilmathi K, Deshmukh PS, Srivastava GC: Physiology and biochemistry of waterlogging tolerance in plants. Biologia Plantarum. 2008, 52: 401-412. 10.1007/s10535-008-0084-6.View ArticleGoogle Scholar
- Ferry N, Edwards MG, Gatehouse J, Capell T, Christou P, Gatehouse AM: Transgenic plants for insect pest control: a forward looking scientific perspective. Transgenic Res. 2006, 15: 13-19. 10.1007/s11248-005-4803-x.View ArticlePubMedGoogle Scholar
- Ye X, Al-Babili S, Kloti A, Zhang J, Lucca P, Beyer P, Potrykus I: Engineering the provitamin A (beta-carotene) biosynthetic pathway into (carotenoid-free) rice endosperm. Science. 2000, 287: 303-305. 10.1126/science.287.5451.303.View ArticlePubMedGoogle Scholar
- Chan Z, Grumet R, Loescher W: Global gene expression analysis of transgenic, mannitol-producing, and salt-tolerant Arabidopsis thaliana indicates widespread changes in abiotic and biotic stress-related genes. J Exp Bot. 2011, 62: 4787-803. 10.1093/jxb/err130.PubMed CentralView ArticlePubMedGoogle Scholar
- Montero M, Coll A, Nadal A, Messeguer J, Pla M: Only half the transcriptomic differences between resistant genetically modified and conventional rice are associated with the transgene. Plant biotechnology journal. 2011, 9: 693-702. 10.1111/j.1467-7652.2010.00572.x.View ArticlePubMedGoogle Scholar
- Coll A, Nadal A, Collado R, Capellades G, Kubista M, Messeguer J, Pla M: Natural variation explains most transcriptomic changes among maize plants of MON810 and comparable non-GM varieties subjected to two N-fertilization farming practices. Plant molecular biology. 2010, 73: 349-362. 10.1007/s11103-010-9624-5.View ArticlePubMedGoogle Scholar
- Wang J, Su X, Ji L, Zhang B, Hu Z, Huang R, Tian Y: Multiple transgenes Populus xeuramericana 'Guariento' plants obtained by biolistic bombardment. Chinese Science Bulletin. 2007, 52: 224-230.View ArticleGoogle Scholar
- Su X, Chu Y, Li H, Hou Y, Zhang B, Huang Q, Hu Z, Huang R, Tian Y: Expression of multiple resistance genes enhances tolerance to environmental stressors in transgenic poplar (Populus × euramericana 'Guariento'). PLoS One. 2011, 6: e24614-10.1371/journal.pone.0024614.PubMed CentralView ArticlePubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature methods. 2008, 5: 621-628. 10.1038/nmeth.1226.View ArticlePubMedGoogle Scholar
- Zhuang J, Cai B, Peng RH, Zhu B, Jin XF, Xue Y, Gao F, Fu XY, Tian YS, Zhao W, Qiao YS, Zhang Z, Xiong AS, Yao QH: Genome-wide analysis of the AP2/ERF gene family in Populus trichocarpa. Biochem Biophys Res Commun. 2008, 371: 468-474. 10.1016/j.bbrc.2008.04.087.View ArticlePubMedGoogle Scholar
- Zenoni S, Ferrarini A, Giacomelli E, Xumerle L, Fasoli M, Malerba G, Bellin D, Pezzotti M, Delledonne M: Characterization of transcriptional complexity during berry development in Vitis vinifera using RNA-Seq. Plant Physiology. 2010, 152: 1787-1795. 10.1104/pp.109.149716.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu X, Liu J: Transcriptional profiling of canker-resistant transgenic sweet orange (Citrus sinensis Osbeck) constitutively overexpressing a spermidine synthase gene. BioMed research international. 2013, 2013:918136: 10.1155/2013/918136.Google Scholar
- Kondrak M, Marincs F, Kalapos B, Juhasz Z, Banfalvi Z: Transcriptome analysis of potato leaves expressing the trehalose-6-phosphate synthase 1 gene of yeast. PloS one. 2011, 6: e23466-10.1371/journal.pone.0023466.PubMed CentralView ArticlePubMedGoogle Scholar
- Caimi PG, McCole LM, Klein TM, Hershey HP: Cytosolic expression of the Bacillus amyloliquefaciens SacB protein inhibits tissue development in transgenic tobacco and potato. New phytologist. 1997, 136: 19-28. 10.1046/j.1469-8137.1997.00719.x.View ArticleGoogle Scholar
- Chourey PS, Nelson OE: Interallelic complementation at the sh locus in maize at the enzyme level. Genetics. 1979, 91: 317-325.PubMed CentralPubMedGoogle Scholar
- Geigenberger P, Stitt M: Sucrose synthase catalyses a readily reversible reaction in vivo in developing potato tubers and other plant tissues. Planta. 1993, 189: 329-339. 10.1007/BF00194429.View ArticlePubMedGoogle Scholar
- Heim U, Weber H, Bäumlein H, Wobus U: A sucrose-synthase gene of Vicia faba L.: expression pattern in developing seeds in relation to starch synthesis and metabolic regulation. Planta. 1993, 191: 394-401.View ArticlePubMedGoogle Scholar
- Cairns AJ: Fructan biosynthesis in transgenic plants. Journal of Experimental Botany. 2003, 54: 549-567. 10.1093/jxb/erg056.View ArticlePubMedGoogle Scholar
- Nelson CJ, Spollen WG: Fructans. Physiologia Plantarum. 1987, 71: 512-516. 10.1111/j.1399-3054.1987.tb02892.x.View ArticleGoogle Scholar
- Van Der Meer IM, Ebskamp MJ, Visser RG, Weisbeek PJ, Smeekens SC: Fructan as a new carbohydrate sink in transgenic potato plants. The Plant Cell Online. 1994, 6: 561-570. 10.1105/tpc.6.4.561.View ArticleGoogle Scholar
- Vereyken IJ, Chupin V, Demel RA, Smeekens S, De Kruijff B: Fructans insert between the headgroups of phospholipids. Biochimica et Biophysica Acta (BBA)-Biomembranes. 2001, 1510: 307-320. 10.1016/S0005-2736(00)00363-1.View ArticleGoogle Scholar
- Vereyken IJ, Albert Van Kuik J, Evers TH, Rijken PJ, de Kruijff B: Structural requirements of the fructan-lipid interaction. Biophysical journal. 2003, 84: 3147-3154. 10.1016/S0006-3495(03)70039-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Vogt T, Jones P: Glycosyltransferases in plant natural product synthesis: characterization of a supergene family. Trends in plant science. 2000, 5: 380-386. 10.1016/S1360-1385(00)01720-9.View ArticlePubMedGoogle Scholar
- Poppenberger B, Berthiller F, Lucyshyn D, Sieberer T, Schuhmacher R, Krska R, Kuchler K, Glössl J, Luschnig C, Adam G: Detoxification of the Fusarium mycotoxin deoxynivalenol by a UDP-glucosyltransferase from Arabidopsis thaliana. Journal of Biological Chemistry. 2003, 278: 47905-47914. 10.1074/jbc.M307552200.View ArticlePubMedGoogle Scholar
- Tognetti VB, Van Aken O, Morreel K, Vandenbroucke K, Van De Cotte B, De Clercq I, Chiwocha S, Fenske R, Prinsen E, Boerjan W: Perturbation of indole-3-butyric acid homeostasis by the UDP-glucosyltransferase UGT74E2 modulates Arabidopsis architecture and water stress tolerance. The Plant Cell Online. 2010, 22: 2660-2679. 10.1105/tpc.109.071316.View ArticleGoogle Scholar
- Chen S, Jiang J, Li H, Liu G: The salt-responsive transcriptome of Populus simonii × Populus nigra via DGE. Gene. 2012, 504: 203-212. 10.1016/j.gene.2012.05.023.View ArticlePubMedGoogle Scholar
- Allen MD, Yamasaki K, Ohme-Takagi M, Tateno M, Suzuki M: A novel mode of DNA recognition by a β-sheet revealed by the solution structure of the GCC-box binding domain in complex with DNA. The EMBO Journal. 1998, 17: 5484-5496. 10.1093/emboj/17.18.5484.PubMed CentralView ArticlePubMedGoogle Scholar
- Li WZ, Zhang HW, Wang JY, Huang RF: Ethylene responsive factors and their application in improving tobacco tolerance to biotic and abiotic stresses. Biotechnol. Bull. 2006, 4: 30-34.Google Scholar
- Wu L, Zhang Z, Zhang H, Wang X, Huang R: Transcriptional modulation of ethylene response factor protein JERF3 in the oxidative stress response enhances tolerance of tobacco seedlings to salt, drought, and freezing. Plant physiology. 2008, 148: 1953-1963. 10.1104/pp.108.126813.PubMed CentralView ArticlePubMedGoogle Scholar
- Li Y, Su X, Zhang B, Huang Q, Zhang X, Huang R: Expression of jasmonic ethylene responsive factor gene in transgenic poplar tree leads to increased salt tolerance. Tree physiology. 2009, 29: 273-279.View ArticlePubMedGoogle Scholar
- Zhang G, Chen M, Chen X, Xu Z, Guan S, Li L, Li A, Guo J, Mao L, Ma Y: Phylogeny, gene structures, and expression patterns of the ERF gene family in soybean (Glycine max L.). Journal of experimental botany. 2008, 59: 4095-4107. 10.1093/jxb/ern248.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu Z, Xia L, Chen M, Cheng X, Zhang R, Li L, Zhao Y, Lu Y, Ni Z, Liu L: Isolation and molecular characterization of the Triticum aestivum L. ethylene-responsive factor 1 (TaERF1) that increases multiple stress tolerance. Plant molecular biology. 2007, 65: 719-732. 10.1007/s11103-007-9237-9.View ArticlePubMedGoogle Scholar
- Xu Z, Ni Z, Liu L, Nie L, Li L, Chen M, Ma Y: Characterization of the TaAIDFa gene encoding a CRT/DRE-binding factor responsive to drought, high-salt, and cold stress in wheat. Molecular Genetics and Genomics. 2008, 280: 497-508. 10.1007/s00438-008-0382-x.View ArticlePubMedGoogle Scholar
- De Boer K, Tilleman S, Pauwels L, Vanden Bossche R, De Sutter V, Vanderhaeghen R, Hilson P, Hamill JD, Goossens A: APETALA2/ETHYLENE RESPONSE FACTOR and basic helix-loop-helix tobacco transcription factors cooperatively mediate jasmonate-elicited nicotine biosynthesis. The Plant Journal. 2011, 66: 1053-1065. 10.1111/j.1365-313X.2011.04566.x.View ArticlePubMedGoogle Scholar
- Fujii H, Verslues PE, Zhu J: Arabidopsis decuple mutant reveals the importance of SnRK2 kinases in osmotic stress responses in vivo. Proceedings of the National Academy of Sciences. 2011, 108: 1717-1722. 10.1073/pnas.1018367108.View ArticleGoogle Scholar
- Wurzinger B, Mair A, Pfister B, Teige M: Cross-talk of calcium-dependent protein kinase and MAP kinase signaling. Plant signaling & behavior. 2011, 6: 8-12. 10.4161/psb.6.1.14012.View ArticleGoogle Scholar
- Luan S: The CBL-CIPK network in plant calcium signaling. Trends in plant science. 2009, 14: 37-42. 10.1016/j.tplants.2008.10.005.View ArticlePubMedGoogle Scholar
- Yang T, Poovaiah BW: Calcium/calmodulin-mediated signal network in plants. Trends in plant science. 2003, 8: 505-512. 10.1016/j.tplants.2003.09.004.View ArticlePubMedGoogle Scholar
- Chen W, Hughes DE, Bailey JE: Intracellular expression of Vitreoscilla hemoglobin alters the aerobic metabolism of Saccharomyces cerevisiae. Biotechnology progress. 1994, 10: 308-313. 10.1021/bp00027a011.View ArticlePubMedGoogle Scholar
- Bulow L, Holmberg N, Lilius G, Bailey JE: The metabolic effects of native and transgenic hemoglobins on plants. Trends Biotechnol. 1999, 17: 21-24. 10.1016/S0167-7799(98)01252-9.View ArticlePubMedGoogle Scholar
- Chereskin BM, Castelfranco PA: Effects of Iron and Oxygen on Chlorophyll Biosynthesis : II. OBSERVATIONS ON THE BIOSYNTHETIC PATHWAY IN ISOLATED ETIOCHLOROPLASTS. Plant Physiol. 1982, 69: 112-116. 10.1104/pp.69.1.112.PubMed CentralView ArticlePubMedGoogle Scholar
- Li X, Peng R, Fan H, Xiong A, Yao Q, Cheng Z, Li Y: Vitreoscilla hemoglobin overexpression increases submergence tolerance in cabbage. Plant cell reports. 2005, 23: 710-715. 10.1007/s00299-004-0872-1.View ArticlePubMedGoogle Scholar
- Häggman H, Frey AD, Ryynänen L, Aronen T, Julkunen Tiitto R, Tiimonen H, Pihakaski Maunsbach K, Jokipii S, Chen X, Kallio PT: Expression of Vitreoscilla haemoglobin in hybrid aspen (Populus tremula× tremuloides). Plant biotechnology journal. 2003, 1: 287-300. 10.1046/j.1467-7652.2003.00027.x.View ArticlePubMedGoogle Scholar
- Wang Z, Xiao Y, Chen W, Tang K, Zhang L: Functional expression of Vitreoscilla hemoglobin (VHb) in Arabidopsis relieves submergence, nitrosative, photo-oxidative stress and enhances antioxidants metabolism. Plant Science. 2009, 176: 66-77. 10.1016/j.plantsci.2008.09.011.View ArticleGoogle Scholar
- Frey AD, Farrés J, Bollinger CJ, Kallio PT: Bacterial hemoglobins and flavohemoglobins for alleviation of nitrosative stress in Escherichia coli. Applied and environmental microbiology. 2002, 68: 4835-4840. 10.1128/AEM.68.10.4835-4840.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Fliss AE, Rao J, Melville MW, Cheetham ME, Caplan AJ: Domain requirements of DnaJ-like (Hsp40) molecular chaperones in the activation of a steroid hormone receptor. Journal of Biological Chemistry. 1999, 274: 34045-34052. 10.1074/jbc.274.48.34045.View ArticlePubMedGoogle Scholar
- Werck-Reichhart D, Bak S, Paquette S: Cytochromes P450. The Arabidopsis book/American Society of Plant Biologists. 2002, 1 (e0028): 10.1199/tab.0028.Google Scholar
- Holford I: Soil phosphorus: its measurement, and its uptake by plants. Australian Journal of Soil Research. 1997, 35: 227-240. 10.1071/S96047.View ArticleGoogle Scholar
- Higgins CF: ABC transporters: from microorganisms to man. Annual review of cell biology. 1992, 8: 67-113. 10.1146/annurev.cb.08.110192.000435.View ArticlePubMedGoogle Scholar
- Martinoia E, Klein M, Geisler M, Bovet L, Forestier C, Kolukisaoglu Ü, Müller-Röber B, Schulz B: Multifunctionality of plant ABC transporters-more than just detoxifiers. Planta. 2002, 214: 345-355. 10.1007/s004250100661.View ArticlePubMedGoogle Scholar
- Ali MA, Plattner S, Radakovic Z, Wieczorek K, Elashry A, Grundler FM, Ammelburg M, Siddique S, Bohlmann H: An Arabidopsis ATPase gene involved in nematode-induced syncytium development and abiotic stress responses. The Plant Journal. 2013, 74: 852-866. 10.1111/tpj.12170.PubMed CentralView ArticlePubMedGoogle Scholar
- Ding X, Gopalakrishnan B, Johnson LB, White FF, Wang X, Morgan TD, Kramer KJ, Muthukrishnan S: Insect resistance of transgenic tobacco expressing an insect chitinase gene. Transgenic Research. 1998, 7: 77-84. 10.1023/A:1008820507262.View ArticlePubMedGoogle Scholar
- Lu Q, Zhang Y, Cao G, Zhang L, Liang G, Lu Y, Wu K, Gao X, Guo Y: A Fragment of Cadherin-Like Protein Enhances Bacillus thuringiensis Cry1B and Cry1C Toxicity to Spodoptera exigua (Lepidoptera: Noctuidae). Journal of Integrative Agriculture. 2012, 11: 628-638. 10.1016/S2095-3119(12)60050-1.View ArticleGoogle Scholar
- Koiwa H, Bressan RA, Hasegawa PM: Regulation of protease inhibitors and plant defense. Trends in Plant Science. 1997, 2: 379-384. 10.1016/S1360-1385(97)90052-2.View ArticleGoogle Scholar
- Ryan CA: The systemin signaling pathway: differential activation of plant defensive genes. Biochimica et Biophysica Acta (BBA)-Protein Structure and Molecular Enzymology. 2000, 1477: 112-121. 10.1016/S0167-4838(99)00269-1.View ArticleGoogle Scholar
- Zhang B, Chen M, Zhang X, Luan H, Diao S, Tian Y, Su X: Laboratory and field evaluation of the transgenic Populus alba × Populus glandulosa expressing double coleopteran-resistance genes. Tree Physiol. 2011, 31: 567-573. 10.1093/treephys/tpr032.View ArticlePubMedGoogle Scholar
- Zhang B, Chen M, Zhang X, Luan H, Tian Y, Su X: Expression of Bt-Cry3A in transgenic Populus alba × P. glandulosa and its effects on target and non-target pests and the arthropod community. Transgenic Res. 2011, 20: 523-532. 10.1007/s11248-010-9434-1.View ArticlePubMedGoogle Scholar
- Kunkel BN, Brooks DM: Cross talk between signaling pathways in pathogen defense. Current opinion in plant biology. 2002, 5: 325-331. 10.1016/S1369-5266(02)00275-3.View ArticlePubMedGoogle Scholar
- Chinnusamy V, Schumaker K, Zhu JK: Molecular genetic perspectives on cross-talk and specificity in abiotic stress signalling in plants. J Exp Bot. 2004, 55: 225-236.View ArticlePubMedGoogle Scholar
- Fujita M, Fujita Y, Noutoshi Y, Takahashi F, Narusaka Y, Yamaguchi-Shinozaki K, Shinozaki K: Crosstalk between abiotic and biotic stress responses: a current view from the points of convergence in the stress signaling networks. Current opinion in plant biology. 2006, 9: 436-442. 10.1016/j.pbi.2006.05.014.View ArticlePubMedGoogle Scholar
- O'Donnell PJ, Truesdale MR, Calvert CM, Dorans A, Roberts MR, Bowles DJ: A novel tomato gene that rapidly responds to wound-and pathogen-related signals. The Plant Journal. 1998, 14: 137-142. 10.1046/j.1365-313X.1998.00110.x.View ArticlePubMedGoogle Scholar
- Roberts MR, Warner SA, Darby R, Lim E, Draper J, Bowles DJ: Differential regulation of a glucosyl transferase gene homologue during defence responses in tobacco. Journal of experimental botany. 1999, 50: 407-410. 10.1093/jxb/50.332.407.View ArticleGoogle Scholar
- Yang H, Shen H, Chen L, Xing Y, Wang Z, Zhang J, Hong M: The OsEBP-89 gene of rice encodes a putative EREBP transcription factor and is temporally expressed in developing endosperm and intercalary meristem. Plant molecular biology. 2002, 50: 379-391. 10.1023/A:1019859612791.View ArticlePubMedGoogle Scholar
- Pieterse CC, Ton J, Van Loon LC: Cross-talk between plant defence signalling pathways: boost or burden?. AgBiotechNet. 2001, 3: 1-8.Google Scholar
- Livak KJ, Schmittgen TD: Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 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.