Next Article in Journal
Small Molecule Inhibitors of White Spot Syndrome Virus: Promise in Shrimp Seedling Culture
Previous Article in Journal
Adsorption of REEs from Aqueous Solution by EDTA-Chitosan Modified with Zeolite Imidazole Framework (ZIF-8)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 7
10.3390/ijms22073448
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Construction of a Full-Length cDNA Over-Expressing Library to Identify Valuable Genes from Populus tomentosa

by
Lingfei Kong
†,
Zeyu Li
†,
Qin Song
,
Xiaohong Li
and
Keming Luo
*
Chongqing Key Laboratory of Plant Resource Conservation and Germplasm Innovation, School of Life Sciences, Southwest University, Chongqing 400715, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2021, 22(7), 3448; https://doi.org/10.3390/ijms22073448
Submission received: 2 March 2021 / Revised: 17 March 2021 / Accepted: 24 March 2021 / Published: 26 March 2021
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Poplar wood is the main source of renewable biomass energy worldwide, and is also considered to be a model system for studying woody plants. The Full-length cDNA Over-eXpressing (FOX) gene hunting system is an effective method for generating gain-of-function mutants. Large numbers of novel genes have successfully been identified from many herbaceous plants according to the phenotype of gain-of-function mutants under normal or abiotic stress conditions using this system. However, the system has not been used for functional gene identification with high-throughput mutant screening in woody plants. In this study, we constructed a FOX library from the Chinese white poplar, Populus tomentosa. The poplar cDNA library was constructed into the plant expression vector pEarleyGate101 and further transformed into Arabidopsis thaliana (thale cress). We collected 1749 T1 transgenic plants identified by PCR. Of these, 593 single PCR bands from different transgenic lines were randomly selected for sequencing, and 402 diverse sequences of poplar genes were isolated. Most of these genes were involved in photosynthesis, environmental adaptation, and ribosome biogenesis based on Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation. We characterized in detail two mutant lines carrying PtoCPCa or PtoWRKY13 cDNA insertions. Phenotypic characterization showed that overexpression of these genes in A. thaliana affected trichome development or secondary cell wall (SCW) deposition, respectively. Together, the Populus-FOX-Arabidopsis library generated in our experiments will be helpful for efficient discovery of novel genes in poplar.

1. Introduction

Higher plants must resist the vagaries of different environments and their associated abiotic stresses, such as salinity, drought, and heavy metal toxicity, due to their sessile lifestyle and perennial growth habit [1,2]. Poplar spp. are considered to be excellent model plants for studying the mechanisms of secondary cell wall (SCW) deposition and wood development [3]. In addition to conventional breeding methods, large-scale screening of key functional genes in growth development and stress tolerance is also important for studying woody plants [4]. Unfortunately, there is a lack of genomic tools for identifying useful genes in different types of poplar.
In recent decades, loss-of-function mutant libraries constructed using T-DNA or transposon insertional mutagenesis approaches have been widely used in many plant species [5,6,7,8,9]. Numerous novel genes have been successfully cloned and characterized. However, due to the redundancy of genes, there may be only slight phenotypic differences between the wild-type (WT) and loss-of-function mutant plants [10,11,12]. Moreover, these methods usually generate recessive mutations, thus demanding a homozygote to display their phenotypes. Obviously, these methods cannot be applied to triploid tree species with long vegetative periods, such as hybrid Populus tomentosa (Chinese white poplar).
Gain-of-function mutation based on gene overexpression is an alternative and effective approach to evaluating gene function because it is unaffected by gene redundancy [13]. The gain-of-function mutant library has been developed using the activation tagging technique, which involves inserting tandem transcriptional enhancers into the plant genomes of many plant species, including Populus spp. [14,15,16,17]. By far the greatest number of valuable genes have been identified in poplar using this method [18,19,20,21]. The disadvantage of the activation tagging technique is that one enhancer can influence the transcription of genes on both sides of the insertion sites [22] and simultaneously increase the transcription of multiple genes, resulting in complex phenotypes. The Full-length cDNA Over-eXpressing (FOX) gene hunting system is one gain-of-function mutation approach, which can be used to avoid this problem [23]. Similar to the activation-tagging method, the FOX hunting system generates dominant mutations by overexpressing target genes [24]. Individual or limited numbers of normalized full-length cDNA segments from tissues or whole plants of interest under the control of the constitutive CaMV 35S promoter are directionally transformed into donor model plants, such as Arabidopsis spp. [24,25]. The ectopic expression of the full-length cDNAs may result in the expression of various phenotypes, which help to characterize the function of these candidate genes in plants [23]. To date, the FOX hunting system has been widely implemented for exploring gene function in Arabidopsis [23,26], rice [27,28], cotton [29], tomato [30], rape [31], and Lotus spp. [32], resulting in a high mutation rate, ranging from 10% to 20% [23,33].
To date, a large-scale FOX hunting system has not been generated in woody plants. However, recently, a xylem-specific mutant library of Populus was established based on this system [34]. According to the study, a total of only 113 transgenic lines were generated, whereas about 80% of these mutant lines were significantly affected in cellulose, lignin, and/or hemicellulose. Therefore, it is largely restricted to studying the function of most genes in poplars using this xylem-specific mutant library. In order to identify the helpful genes in all kinds of development processes and construct a genome-wide mutant library of Populus, in this study, we used the Gateway technology to generate a normalized genome-wide full-length cDNA library of P. tomentosa [35], and then the Populus FOX hunting library was transformed into Arabidopsis thaliana via the floral-dip method [36]. All of the transgenic plants were screened by phenotypic investigation. As a result, we collected 1749 positive T1 transgenic plants. Two mutant lines containing different poplar genes were selected for further functional exploration. Taken together, our results provide a rapid method for gene screening, which will provide insights into gene function in the complete genome of P. tomentosa.

2. Results

2.1. Construction of Populus-FOX-Arabidopsis Library

In order to eliminate the variance among individuals, we selected three independent WT plants of P. tomentosa to extract total RNA. After constructing the FOX library, 12 monoclonal colonies were randomly selected from the entry clone library and expression clone library and subjected to PCR determination for quality testing. As a result, the inserted fragment lengths ranged from 1.0 to 3.0 kilobases (kb) in both the entry and expression libraries (Figure 1A,B). Moreover, 11 colonies were positive for the insert, which is equivalent to a positive clone rate for the expression library of approximately 91.67% (Figure 1B). From the data, we can conclude that the expression library was effective.
To identify the transgenic plants, all T1 seedlings were sprayed with the herbicide Basta (phosphinothricin) (35% v/v) (Figure S1A,B). In total, we identified 2098 Basta-resistant plants, and 1749 of them were positive transgenic plants, which was confirmed by PCR determination (Figure S1C). Statistical analysis showed that the positive rate was up to 83.36%. We also compared the inserted fragment lengths randomly selected from transgenic plants with those from the expression clone library. The results showed that, in addition to the reduction in the enrichment of fragments longer than 2 kb in transgenic plants, the lengths of other fragments were proportional to those of the expression clone library (Figure 1C). To reveal the genes inserted into transgenic plants, a single PCR band was randomly selected for sequencing and isolated using a Basis Local Alignment Sequencing Tool (BLAST) search. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Orthology database was used to summarize the function of these genes (Figure 1D). A total of 312 gene sequences were mapped to 47 pathways, related mainly to photosynthesis, environmental adaptation, ribosome biogenesis, glutathione metabolism, protein catabolic, embryo development, and mitogen-activated protein kinase (MAPK) signaling. In general, the Populus-FOX-Arabidopsis library has been successfully established.

2.2. Screening and Identification of Functional Genes in the FOX Library

In order to explore the genes affecting plant development via the FOX library, we observed and identified the phenotypes of the T2 generation from sequenced positive seedlings. The great majority of mutant plants showed distinct phenotypes, including lodging, dwarfing, increased branching, glabrous leaves, and so on (Figure 2A–L). Of the 402 FOX lines, 262 individuals showed obvious phenotypes. Interestingly, we found that many mutants exhibited multiple phenotypes, such as dwarfing associated with poor growth and small leaves associated with increased branching (Table 1), suggesting that these genes play important roles in plant growth and development.
P. tomentosa is a perennial tree species with broad ecological adaptations. For the purpose of screening genes that are conducive to stress tolerance, we planted the seeds of T2 transgenic plants and inspected the phenotypic modulation in traits such as root length as well as seed germination under different chemical treatments (Figure S2). Under normal conditions, several FOX lines exhibited longer root lengths (Figure 3A). In 150 mM NaCl, compared with the WT, the root growth of some FOX plants was greatly affected, while the other FOX lines were still able to grow well (Figure 3B). In 10 μM abscisic acid (ABA), some FOX plants exhibited an ABA-hypersensitive phenotype (Figure 3C). In 350 mM mannitol, some FOX plants produced fewer lateral roots, and the root lengths were shorter, which was a typical drought-escape response (Figure 3D). As a result, we screened out 89 stress-responsive plants after processing.

2.3. Characterization of Two Mutants with Different Phenotypes to Evaluate the FOX Library

To verify that the phenotypic changes were the result of the ectopic expression of the corresponding inserted genes, we randomly selected two genes for further functional characterization after comparing different phenotypes and identifying homozygous genes.
In all of the transgenic plants with altered variable leaves, K103 showed unique glabrous leaves. BLAST searching suggested that the inserted gene was a homolog of Potri.002G168900. Amino acid sequence alignment and phylogenetic tree analysis showed that the gene shares a high degree of homology with AtCPC (Figure S3) in A. thaliana, so we named it PtoCPCa. Compared to the WT plants, seedlings that overexpressed PtoCPCa presented smooth glabrous leaves (Figure 4A). Quantitative measurement of leaf trichomes by scanning electron microscopy (SEM) revealed a significant reduction in trichome density. The trichome density of PtoCPCa-OE seedlings was nearly 0/mm2, while the density of WT seedlings was 10.33/mm2 (Figure 4B). Trichome initiation and development are modulated by a complicated network that includes hormone signaling (gibberellin, jasmonic acid) [37,38], a MYB-bHLH-WD40 complex hub comprising GL1, TTG1, GL3, and EGL3 [39,40,41], as well as some key transcription factors, such as TRY, CPC, GL2, TTG2, etc. [42,43,44]. According to the qRT-PCR analysis we carried out, the relative expression level of the PtoCPCa gene was extremely high in the K103 lines (Figure 4C), while the levels of GL1, GL2, and TTG2 were reduced in the PtoCPCa-OE lines compared to those of the WT plants (Figure 4D). Suppression of these key factors may have been the cause of aberrant trichome development. This is consistent with the results of a previous study of Arabidopsis [42]. In addition, PtrTCL1 (homolog of PtoCPCa) ectopic expression in Arabidopsis also resulted in a similar phenotype [45]. This result suggested that our Populus-FOX-Arabidopsis library is an effective method for exploring functional genes.
We also wondered whether our FOX library could hunt out some novel genes. Indeed, according to the annotations, there are still some genes with unknown functions (Table 1). We found two individual mutant lines (K23 and K78) displaying the parallel phenotype of stiffer inflorescence stems, and both of them had the same gene Potri.005G086400 inserted, which encodes a WRKY transcription factor that had not previously been reported in poplar. Phylogenetic analysis indicated that Potri.005G086400 is the homologous gene of AtWRKY13 (Figure S4A), and we named it PtoWRKY13, based on a previous study [46]. Amino acid sequence alignment suggested that, except for a conserved WRKY domain, only a few amino acid residues showed no difference between PtoWRKY13 and AtWRKY13 (Figure S4B). It has been reported that AtWRKY13 affects SCW deposition by directed modulation of NST2 in Arabidopsis [47]. After staining with Toluidine Blue, observation of the basal cross-sections of the inflorescence stems from 6-week-old plants showed that, although PtoWRKY13-OE lines exhibit similar stem diameters and numbers of vascular bundles (Figure 5B,C), the SCWs were exaggeratedly thicker than those of WT plants (Figure 5A,D,E). We first detected the expression level of PtoWRKY13 in the two mutants by qRT-PCR analysis, and the results showed that both FOX lines had a high expression level of PtoWRKY13 (Figure 6A). In A. thaliana, the wrky13 mutant showed the phenotype of lignin reduction. Hence, we examined the expression of SCW biosynthetic genes in the PtoWRKY13-OE lines. The expression levels of four lignin pathway genes (HCT, PAL4, CAD6, and F5H) and four transcription factors (SND1, NST1, NST2, and C3H14) were all up-regulated in the PtoWRKY13-OE lines (Figure 6B,C). The results were fascinating because of the great differences between the amino acid sequences of PtoWRKY13 and AtWRKY13. It is interesting to further verify the phenotype and regulatory mechanism of PtoWRKY13 in poplar, as it may prove to be an important factor in improving the properties of poplar wood.

3. Discussion

Gene function studies rely on the creation of mutant materials under existing experimental conditions. Up until now, the whole genome loss-of-function mutant libraries have successfully been developed in several plant species, including A. thaliana and rice, which greatly enable the discovery of functional genes [6,48]. With regard to polyploid plants, loss-of-function mutant generation via T-DNA, transposon insertion, and physical and chemical mutagenesis usually deprives the function of a limited number of genes in multigene families. The other members of the family may play similar roles in rescuing the function, so that the mutant appears to possess no observed phenotype. Furthermore, T-DNA insertion is a random event that occurs on the genome, which makes it difficult to screen out helpful mutants. Therefore, gain-of-function mutagenesis is a better choice for studying polyploid plants, of which the FOX hunting system is a representative technique. Compared with the loss-of-function method, it is convenient to discover the inserted genes from the FOX hunting mutant by PCR amplification and sequencing using universal primer pairs with no need for positional cloning or TAIL-PCR [23,27]. Moreover, the FOX system is a versatile tool for constructing a genome-wide library or specific development process-related library [26,34] due to the fact that the number of cDNAs in the library can be altered as needed. The screening of the FOX library can also be undertaken in conditions of interest, for example, high Zn concentrations [49], high salinity, or drought [29]. To date, the FOX hunting system has been widely used in polyploid plants such as cotton and rape [29,31]. These studies suggest that FOX is an effective method for exploring functional genes.
In this study, the genome-wide full-length cDNA of P. tomentosa was cloned to the pEarleyGate101 plant expression vector containing the Basta-resistant gene using the high-throughput Gateway technique. This FOX library was then transformed into A. thaliana to facilitate efficient screening of functional genes. As a result, we transformed more than 10,000 WT A. thaliana plants and obtained 2098 Basta-resistant plants, including 1749 positive transgenic plants by exploiting the Bar gene on the pEarleyGate100 binary expression vector as a selection marker and PCR identification, respectively (Supplementary Figure S1B,C). The positive rate was up to 83.36%, and the inserted fragment lengths ranged from 1to 3 kb. These results were consistent with those of previous studies using this technique. For example, the inserted fragments of a cotton-FOX-Arabidopsis overexpression library varied from 0.5 to 3.0 kb [29]; the average length of fragments from an Oryza-FOX-Arabidopsis plant was 1.5–2.0 kb [33]; in A. thaliana, the length ranged from 0.3 to 3.0 kb [23].
Rapidity and validity are two important metrics for a gene discovery technique. In our experiments, various phenotypes were observed, including lodging, wrinkled leaves, greater branching, and so on in the T1 generation (Figure 2A–L). Most of the phenotypes were inserted into a signal fragment, confirmed by PCR amplification. The results above indicated that we screened out a series of functional genes using the FOX system. Indeed, functional characterization of a randomly selected K103 line found that the phenotype was caused by overexpression of the corresponding inserted gene. K103 showed glabrous leaves, and blast searching and phylogenetic analysis demonstrated that PtoCPCa was inserted into K103. As shown in Figure 3D, the transcription levels of the trichome development key factors GL1, GL2, and TTG2 were reduced. This result is consistent with those of previous studies of A. thaliana [42], and were verified by experiments on P. trichocarpa [45]. Trichome development was affected by overexpression of eight poplar R3-MYB genes (including PtrCPCa, which was named PtrTCL1). Interestingly, one of the trichome development key transcription factors, AtGTL1, was first identified by FOX library screening [50], indicating that the FOX hunting system might be the unique method to identify various key factors from the specific development regulation network. In conclusion, we successfully screened out numerous functional genes via our FOX library.
Except for genes that had already been reported, we also hunted for a large number of novel genes. PtoWRKY13 can be considered a representative of genes with unknown functions. The formation of wood is an important and characteristic function of poplar, and the development of wood comes from the differentiation of the cambium cells into xylem cells, which then go on to cell expansion, SCW deposition, programmed cell death, and eventually wood formation [51]. The properties of wood are affected mainly by the deposition of SCWs. It is obvious that some genes of the NAC family play master switch roles in the regulation of SCW formation [52,53], while the existence of regulators that function upstream remains unknown. According to various studies, a few WRKY transcription factor family members play such a role. It has been reported that AtWRKY12 initiates pith SCW formation by directly regulating NST2 [54], and AtWRKY13 functions as an activator to directly regulate NST2 [47]. Additionally, AtWRKY15 suppresses tracheary element differentiation by regulating the upstream element VND7, possibly via auxin signaling [55]. The regulation of auxin in SCW formation has been widely reported [56,57] and might operate mainly through ASYMMETRIC LEAVES2 (AS2)/LATERAL ORGAN BOUNDARIES DOMAIN (LBD) family proteins. LBD29, which involves auxin signaling, represses the NAC master regulators and fiber wall biosynthesis [58], while LBD30/18 genes act downstream of ARF7 and are involved in a positive feedback loop for VND7 expression [59]. PtoWRKY13 overexpression lines K25/K78 showed stiffer stems and thickened SCWs (Figure 6A). Levels of lignin biosynthetic genes, including HCT, PAL4, CAD6, and F5H, were increased as well as those of SCW master switches SND1, NST1, and NST2 (Figure 6B,C). However, it is still not clear whether PtoWRKY13 directly regulates WNDs and the relationship with LBD proteins in P. tomentosa. In order to solve this problem, transgenic poplar will be generated in order to verify the phenotype and reveal the molecular mechanism in future studies.

4. Materials and Methods

4.1. Plant Materials and Growth Conditions

Populus tomentosa Carr. (clone 741) were grown in the greenhouse at 22 °C with a 16/8 h light/dark cycle with supplemental light (4000 lux), and the relative humidity should be around 60% for the optimum growth condition. Arabidopsis thaliana (Col-0) seedlings were grown in an incubator at 22 °C with 16 h 10,000 lux light and 8 h dark. The relative humidity was kept at 60% for P. tomentosa and 80% for Arabidopsis thaliana.
In order to observe the difference in sensitivity to salt stress during seed germination, seeds were evenly placed on MS medium supplemented with 150 mM NaCl. In order to observe the difference in sensitivity to stress during seedlings growth, seeds were evenly placed on MS medium supplemented with 3% (w/v) sucrose and 10 μM ABA, 150 mM NaCl, or 200 mM mannitol. Statistical germination rates were calculated from 1 to 9 days after the earliest germination. Root lengths were measured after 10 days of vertical culture.
The seeds of WT Col-0, homozygous PtoCPCa, and PtoWRKY13 overexpressing lines were kept at 4 °C for 3 days before being placed on MS medium supplemented with 3% (w/v) sucrose. After 10 days of germination on MS plates, the seedlings were transferred into soil and developed in a growth incubator at 22 °C under long-day conditions (16 h light/8 h dark).

4.2. RNA Extraction

Total RNA of poplar was extracted and purified using RNA RNeasy Plant Mini Kit (Qiagen, Hilden, Germany). The OD260/OD280 ratio measured via a NanoDrop 2000 Spectrophotometer (Thermo, West Palm Beach, FL, USA) and electrophoresis were used to judge the RNA quality. Total RNA (2 µg) treated with RNase-free DNase (Takara Dalian, China) was subjected to first-strand cDNA synthesis in a total volume of 20 µl using RT-AMV transcriptase (Takara, Dalian, China).

4.3. Construction of P. tomentosa Fl-cDNA Overexpression Library, Plant Transformation and Selection

Populus full-length cDNA was synthesized by reverse transcription with the isolated mRNA. DSN (Duplex-Specific Nuclease) enzyme was used to normalize the concentration of cDNA to generate a normalized FLcDNA library. The cDNA library was firstly cloned into an entry vector pDONR222, and then into an pEarleyGate101 plant expression vector using high-throughput Gateway technology
Populus-FOX-Arabidopsis library was produced by using the Gateway cloning technology [35] (Invitrogen, Carlsbad, CA, USA). Gateway technology contains two steps: BP recombination reaction (PCR fragment+ entry vector) and LR recombination reaction (entry vector after BP recombination reaction+ destination vector). After LR recombination reaction, the cDNAs were recombined into overexpression vector pEarleyGate101 (Invitrogen, Carlsbad, CA, USA), and the resulting construct resistance to glufosinate was introduced into Agrobacterium tumefaciens strain GV3101 by the freeze–thaw method.
The expression vector library was then transformed into Agrobacterium tumefaciens strain GV3101 and introduced into WT Arabidopsis seedlings using floral-dip method [36]. Subsequently, leaves selected from Basta-resistant T1 plants were used for DNA extraction and PCR analysis to confirm cDNA diversity. Transformed T1 seeds were selected with 350 mg/L of Basta solution. Phenotypes were scored based on morphological changes such as germination, leave size, shape, and color. All plants showing visible phenotypes were transferred to a new growing tray.

4.4. Genomic DNA Isolation, PCR and Sequencing

To identify integrated cDNAs, genomic DNA prepared from leaves of randomly selected 2100 T1 transgenic plants was the template for PCR amplification with primers complementary to vector sequences flanking the attB1 and attB2 sites. The PCR condition was 94 °C for 30 s for denaturation, 57 °C for 30 s for annealing and 72 °C for 120 s for elongation. The PCR products were gel purified and sequenced with the same primers. The identity of the transcript was revealed by sequence homology search using the Phytozome BLAST tool. To validate the phenotype conferred by the inducible expression of fl-cDNA, the cDNA was isolated and inserted into pEarleyGate101 expression vector driven by 35S promoter for Agrobacterium-mediated transformation of Arabidopsis wild-type (Col-0).

4.5. Phylogenetic Analysis and Sequence Alignment

A phylogenetic tree was developed with MEGA6 (Lynnon Biosoft, Quebec, QC, Canada) using the neighbor-joining method (Tamura et al. 2013). Multiple alignments of protein sequences were conducted with the software DNAMAN 7 (Lynnon Biosoft, San Ramon, CA, USA).

4.6. Gene Expression Analysis

For quantitative real time PCR (qRT-PCR) analysis, total RNA from leaf tissues or stem tissues was extracted using RNA RNeasy Plant Mini Kit (Qiagen, Hilden, Germany) and treated with RNase-free DNase (TaKaRA, Dalian, China). Samples from at least three plants were pooled for analysis. The quality or integrity of RNA was checked by agarose gel electrophoresis and P100 spectrophotometer (Pultton, Ann Arbor, MI, USA). The criteria of high-quality total RNA include the following: (1) sharp distinct 28S and 18S rRNA bands, with the 28S band approximately twice as intense as the 18S band; (2) the value of D260/OD280 between 1.9 and 2.0; and (3) no detected genomic DNA band. The qualified RNA samples were reversely transcribed using RT-AMV (Avian Myeloblastosis Virus) transcriptase (TaKaRa, Dalian, China). Subsequently, qRT-PCR was performed in a volume of 20 μL containing 10 μL of SYBR Premix ExTaq TM (TaKaRa, Dalian, China). ACTIN7 and UBQ were used as an internal control in Arabidopsis and Populus, respectively. Three biological replicates of each sample and three technical replicates of each biological experiment were conducted. Primers used for qRT-PCR were listed in Table S1.

4.7. Microscopy and Histochemistry

The basal inflorescence stems of 6-week-old WT and transgenic Arabidopsis seedlings were cross-sectioned by using a Vibrating blade microtome (VT1000S, Leica, Weitzlar, Germany). To detect lignified cell walls, stem cross-sections (100 μm thickness) were stained with Toluidine Blue-O (TBO, 0.5% w/v). Observations were performed using an Olympus BX53 microscope.

4.8. Scanning Electron Microscopy (SEM) Observation

The abaxial epidermis of 3-week-old WT and PtoCPCa-OE transgenic Arabidopsis leaves were attached to sample stage using double-sided sticky tape and then directly viewed under a scanning electron microscope (Phenom™ Pure FEI, Eindhoven, The Netherland). The images were captured digitally. At least three biological replicates were measured for each test.

5. Conclusions

In this study, a genome-wide full-length Populus-FOX-Arabidopsis overexpression library was constructed using Gateway technology. A total of 1749 positive transgenic plants were obtained and various phenotypes were observed, such as lodging, poor growth, dwarfism, and so on in the T1 generation. Genes screened out from the library were related mainly to photosynthesis, environmental adaptation, and ribosome biogenesis based on KEGG pathway annotation. A total of 89 stress-responsive plants were screened out under NaCl, mannitol, or ABA conditions. This library provides a novel resource for exploring genes related to perennial growth and abiotic stress tolerance. To evaluate the FOX library, two mutant lines carrying PtoCPCa or PtoWRKY13 cDNA insertions were selected for further functional exploration in trichome development and secondary wall deposition respectively. The results demonstrated that our FOX library could be used to identify valuable genes from various development processes. In brief, our study showed a rapid gene discovery method that will provide important insights into the gene function of poplar in the whole genome. Future work will be focused mainly on novel gene characterization and molecular mechanism analysis.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22073448/s1. Figure S1: Screening and identification of positive FOX transgenic Arabidopsis thaliana seedlings. Figure S2: Screening of salt tolerance FOX mutants. Figure S3: Phylogenetic analysis and amino acid alignment of PtoCPCa and other R3 MYB transcription factors from different species. Figure S4: Phylogenetic analysis and amino acid alignment of PtoWRKY13 and other WRKY transcription factors from different species. Supplementary Data Set 1: Primers used in this article.

Author Contributions

K.L. conceived the research; L.K., Z.L., X.L. and Q.S. performed the experiments; L.K. drafted the manuscript; K.L. corrected the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key Project for Research on Transgenic Plant (2016ZX08010-003 and 2016ZX08010-003), the National Natural Science Foundation of China (31870657, 31870175, 32071791, and 31800505). We also obtained financial support from the Fundamental Research Funds for the Central Universities of China (2572018CL01).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Nannan Li and Shengting Li (Southwest University, China) for providing technical help with the Arabidopsis transformation.

Conflicts of Interest

The authors declare that the research was carried out without any potential conflict of interest.

References

  1. Miller, G.; Schlauch, K.; Tam, R.; Cortes, D.; Torres, M.A.; Shulaev, V.; Dangl, J.L.; Mittler, R. The Plant NADPH Oxidase RBOHD Mediates Rapid Systemic Signaling in Response to Diverse Stimuli. Sci. Signal. 2009, 2, ra45. [Google Scholar] [CrossRef] [Green Version]
  2. Harfouche, A.; Meilan, R.; Altman, A. Molecular and physiological responses to abiotic stress in forest trees and their relevance to tree improvement. Tree Physiol. 2014, 34, 1181–1198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jansson, S.; Douglas, C.J. Populus: A model system for plant biology. Annu. Rev. Plant Biol. 2007, 58, 435–458. [Google Scholar] [CrossRef] [PubMed]
  4. Osakabe, Y.; Kajita, S.; Osakabe, K. Genetic engineering of woody plants: Current and future targets in a stressful environment. Physiol. Plant 2011, 142, 105–117. [Google Scholar] [CrossRef]
  5. Alonso, J.M.; Stepanova, A.N.; Leisse, T.J.; Kim, C.J.; Chen, H.; Shinn, P.; Stevenson, D.K.; Zimmerman, J.; Barajas, P.; Cheuk, R. Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science 2003, 301, 653–657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Piffanelli, P.; Droc, G.; Mieulet, D.; Lanau, N.; Bès, M.; Bourgeois, E.; Rouvière, C.; Gavory, F.; Cruaud, C.; Ghesquière, A.; et al. Large-scale characterization of Tos17 insertion sites in a rice T-DNA mutant library. Plant Mol. Biol. 2007, 65, 587–601. [Google Scholar] [CrossRef] [PubMed]
  7. Ram, H.; Soni, P.; Salvi, P.; Gandass, N.; Sharma, A.; Kaur, A.; Sharma, T.R. Insertional Mutagenesis Approaches and Their Use in Rice for Functional Genomics. Plants 2019, 8, 310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Nandety, R.S.; Serrani-Yarce, J.C.; Gill, U.S.; Oh, S.; Lee, H.K.; Zhang, X.; Dai, X.; Zhang, W.; Krom, N.; Wen, J.; et al. Insertional mutagenesis of Brachypodium distachyon using the Tnt1 retrotransposable element. Plant J. 2020, 103, 1924–1936. [Google Scholar] [CrossRef]
  9. Tadege, M.; Wen, J.; He, J.; Tu, H.; Kwak, Y.; Eschstruth, A.; Cayrel, A.; Endre, G.; Zhao, P.X.; Chabaud, M.; et al. Large-scale insertional mutagenesis using the Tnt1 retrotransposon in the model legume Medicago truncatula. Plant J. 2008, 54, 335–347. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, X.; Liu, S.; Tian, H.; Wang, S.; Chen, J.-G. The Small Ethylene Response Factor ERF96 is Involved in the Regulation of the Abscisic Acid Response in Arabidopsis. Front. Plant Sci. 2015, 6. [Google Scholar] [CrossRef] [Green Version]
  11. Jiao, B.; Zhao, X.; Lu, W.; Guo, L.; Luo, K. The R2R3 MYB transcription factor MYB189 negatively regulates secondary cell wall biosynthesis in Populus. Tree Physiol. 2019, 39, 1187–1200. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, X.; Yu, R.; Wang, J.; Lin, Z.; Han, X.; Deng, Z.; Fan, L.; He, H.; Deng, X.W.; Chen, H. The Asymmetric Expression of SAUR Genes Mediated by ARF7/19 Promotes the Gravitropism and Phototropism of Plant Hypocotyls. Cell Rep. 2020, 31, 107529. [Google Scholar] [CrossRef]
  13. Zhang, J.Z. Overexpression analysis of plant transcription factors. Curr. Opin. Plant Biol. 2003, 6, 430–440. [Google Scholar] [CrossRef]
  14. Waki, T.; Miyashima, S.; Nakanishi, M.; Ikeda, Y.; Hashimoto, T.; Nakajima, K. A GAL4-based targeted activation tagging system inArabidopsis thaliana. Plant J. 2013, 73, 357–367. [Google Scholar] [CrossRef] [PubMed]
  15. Manimaran, P.; Venkata Reddy, S.; Moin, M.; Raghurami Reddy, M.; Yugandhar, P.; Mohanraj, S.S.; Balachandran, S.M.; Kirti, P.B. Activation-tagging in indica rice identifies a novel transcription factor subunit, NF-YC13 associated with salt tolerance. Sci. Rep. 2017, 7, 9341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Davies, J.P.; Reddy, V.S.; Liu, X.L.; Reddy, A.S.; Ainley, W.M.; Folkerts, O.; Marri, P.; Jiang, K.; Wagner, D.R. Development of an activation tagging system for maize. Plant Direct 2019, 3, e00118. [Google Scholar] [CrossRef]
  17. Fladung, M.; Polak, O. Ac/Ds-transposon activation tagging in poplar: A powerful tool for gene discovery. BMC Genom. 2012, 13, 61. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Busov, V.B.; Meilan, R.; Pearce, D.W.; Ma, C.; Rood, S.B.; Strauss, S.H. Activation tagging of a dominant gibberellin catabolism gene (GA 2-oxidase) from poplar that regulates tree stature. Plant Physiol. 2003, 132, 1283–1291. [Google Scholar] [CrossRef] [Green Version]
  19. Trupiano, D.; Yordanov, Y.; Regan, S.; Meilan, R.; Tschaplinski, T.; Scippa, G.S.; Busov, V. Identification, characterization of an AP2/ERF transcription factor that promotes adventitious, lateral root formation in Populus. Planta 2013, 238, 271–282. [Google Scholar] [CrossRef]
  20. Yordanov, Y.S.; Ma, C.; Strauss, S.H.; Busov, V.B. EARLY BUD-BREAK 1 (EBB1) is a regulator of release from seasonal dormancy in poplar trees. Proc. Natl. Acad. Sci. USA 2014, 111, 10001–10006. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Dash, M.; Yordanov, Y.S.; Georgieva, T.; Tschaplinski, T.J.; Yordanova, E.; Busov, V. Poplar PtabZIP1-like enhances lateral root formation and biomass growth under drought stress. Plant J. 2017, 89, 692–705. [Google Scholar] [CrossRef] [Green Version]
  22. Ichikawa, T.; Nakazawa, M.; Kawashima, M.; Muto, S.; Gohda, K.; Suzuki, K.; Ishikawa, A.; Kobayashi, H.; Yoshizumi, T.; Tsumoto, Y. Sequence database of 1172 T-DNA insertion sites in Arabidopsis activation-tagging lines that showed phenotypes in T1 generation. Plant J. 2003, 36, 421–429. [Google Scholar] [CrossRef] [Green Version]
  23. Ichikawa, T.; Nakazawa, M.; Kawashima, M.; Iizumi, H.; Kuroda, H.; Kondou, Y.; Tsuhara, Y.; Suzuki, K.; Ishikawa, A.; Seki, M. The FOX hunting system: An alternative gain-of-function gene hunting technique. Plant J. 2006, 48, 974–985. [Google Scholar] [CrossRef] [PubMed]
  24. 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–393. [Google Scholar] [CrossRef]
  25. Curtis, M.D.; Grossniklaus, U. A gateway cloning vector set for high-throughput functional analysis of genes in planta. Plant Physiol. 2003, 133, 462–469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Fujita, M.; Mizukado, S.; Fujita, Y.; Ichikawa, T.; Nakazawa, M.; Seki, M.; Matsui, M.; Yamaguchi-Shinozaki, K.; Shinozaki, K. Identification of stress-tolerance-related transcription-factor genes via mini-scale Full-length cDNA Over-eXpressor (FOX) gene hunting system. Biochem. Biophys. Res. Commun. 2007, 364, 250–257. [Google Scholar] [CrossRef]
  27. Kondou, Y.; Higuchi, M.; Takahashi, S.; Sakurai, T.; Ichikawa, T.; Kuroda, H.; Yoshizumi, T.; Tsumoto, Y.; Horii, Y.; Kawashima, M. Systematic approaches to using the FOX hunting system to identify useful rice genes. Plant J. 2009, 57, 883–894. [Google Scholar] [CrossRef] [PubMed]
  28. Sakurai, T.; Kondou, Y.; Akiyama, K.; Kurotani, A.; Higuchi, M.; Ichikawa, T.; Kuroda, H.; Kusano, M.; Mori, M.; Saitou, T.; et al. RiceFOX: A Database of Arabidopsis Mutant Lines Overexpressing Rice Full-Length cDNA that Contains a Wide Range of Trait Information to Facilitate Analysis of Gene Function. Plant Cell Physiol. 2011, 52, 265–273. [Google Scholar] [CrossRef] [PubMed]
  29. Li, S.; Chen, H.; Hou, Z.; Li, Y.; Yang, C.; Wang, D.; Song, C.P. Screening of abiotic stress-responsive cotton genes using a cotton full-length cDNA overexpressing Arabidopsis library. J. Integr. Plant Biol. 2019, 62, 998–1016. [Google Scholar] [CrossRef]
  30. Wang, Y.-H.; Mosebach, C.M.; Kibbey, A.S.; Ryhal, M.K.; Jones, A.D.; Palmer, J.A.; Kochian, L.V. Generation of Arabidopsis Mutants by Heterologous Expression of a Full-Length cDNA Library from Tomato Fruits. Plant Mol. Biol. Rep. 2009, 27, 454–461. [Google Scholar] [CrossRef]
  31. Ling, J.; Li, R.; Nwafor, C.C.; Cheng, J.; Li, M.; Xu, Q.; Wu, J.; Gan, L.; Yang, Q.; Liu, C.; et al. Development of iFOX-hunting as a functional genomic tool and demonstration of its use to identify early senescence-related genes in the polyploid Brassica napus. Plant Biotechnol. J. 2018, 16, 591–602. [Google Scholar] [CrossRef] [Green Version]
  32. Himuro, Y.; Tanaka, H.; Hashiguchi, M.; Ichikawa, T.; Nakazawa, M.; Seki, M.; Fujita, M.; Shinozaki, K.; Matsui, M.; Akashi, R.; et al. FOX-superroots of Lotus corniculatus, overexpressing Arabidopsis full-length cDNA, show stable variations in morphological traits. J. Plant Physiol. 2011, 168, 181–187. [Google Scholar] [CrossRef]
  33. Nakamura, H.; Hakata, M.; Amano, K.; Miyao, A.; Toki, N.; Kajikawa, M.; Pang, J.; Higashi, N.; Ando, S.; Toki, S.; et al. A genome-wide gain-of function analysis of rice genes using the FOX-hunting system. Plant Mol. Biol. 2007, 65, 357–371. [Google Scholar] [CrossRef] [Green Version]
  34. Rauschendorfer, J.; Yordanov, Y.; Dobrev, P.; Vankova, R.; Sykes, R.; Kulheim, C.; Busov, V. Overexpression of a developing xylem cDNA library in transgenic poplar generates high mutation rate specific to wood formation. Plant Biotechnol. J. 2020, 18, 1434–1443. [Google Scholar] [CrossRef] [PubMed]
  35. Reece-Hoyes, J.S.; Walhout, A.J.M. Gateway Recombinational Cloning. Cold Spring Harb. Protoc. 2018, 2018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Zhang, X.; Henriques, R.; Lin, S.-S.; Niu, Q.-W.; Chua, N.-H. Agrobacterium-mediated transformation of Arabidopsis thaliana using the floral dip method. Nature Protoc. 2006, 1, 641–646. [Google Scholar] [CrossRef] [PubMed]
  37. An, L.; Zhou, Z.; Su, S.; Yan, A.; Gan, Y. GLABROUS INFLORESCENCE STEMS (GIS) is required for trichome branching through gibberellic acid signaling in Arabidopsis. Plant Cell Physiol. 2012, 53, 457–469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Qi, T.; Song, S.; Ren, Q.; Wu, D.; Huang, H.; Chen, Y.; Fan, M.; Peng, W.; Ren, C.; Xie, D. The Jasmonate-ZIM-domain proteins interact with the WD-Repeat/bHLH/MYB complexes to regulate Jasmonate-mediated anthocyanin accumulation and trichome initiation in Arabidopsis thaliana. Plant Cell 2011, 23, 1795–1814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Bloomer, R.H.; Juenger, T.E.; Symonds, V.V. Natural variation in GL1 and its effects on trichome density in Arabidopsis thaliana. Mol. Ecol. 2012, 21, 3501–3515. [Google Scholar] [CrossRef] [PubMed]
  40. Morohashi, K.; Zhao, M.; Yang, M.; Read, B.; Lloyd, A.; Lamb, R.; Grotewold, E. Participation of the Arabidopsis bHLH factor GL3 in trichome initiation regulatory events. Plant Physiol. 2007, 145, 736–746. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Shen, B.; Sinkevicius, K.W.; Selinger, D.A.; Tarczynski, M.C. The homeobox gene GLABRA2 affects seed oil content in Arabidopsis. Plant Mol. Biol. 2006, 60, 377–387. [Google Scholar] [CrossRef] [PubMed]
  42. Schellmann, S.; Schnittger, A.; Kirik, V.; Wada, T.; Okada, K.; Beermann, A.; Thumfahrt, J.; Jurgens, G.; Hulskamp, M. TRIPTYCHON and CAPRICE mediate lateral inhibition during trichome and root hair patterning in Arabidopsis. EMBO J. 2002, 21, 5036–5046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kirik, V.; Simon, M.; Huelskamp, M.; Schiefelbein, J. The ENHANCER OF TRY AND CPC1 gene acts redundantly with TRIPTYCHON and CAPRICE in trichome and root hair cell patterning in Arabidopsis. Dev. Biol. 2004, 268, 506–513. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Johnson, C.S.; Kolevski, B.; Smyth, D.R. TRANSPARENT TESTA GLABRA2, a Trichome and Seed Coat Development Gene of Arabidopsis, Encodes a WRKY Transcription Factor. Plant Cell 2002, 14, 1359–1375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Zhou, L.; Zheng, K.; Wang, X.; Tian, H.; Wang, X.; Wang, S. Control of trichome formation in Arabidopsis by poplar single-repeat R3 MYB transcription factors. Front. Plant Sci. 2014, 5, 262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Jiang, Y.; Duan, Y.; Yin, J.; Ye, S.; Zhu, J.; Zhang, F.; Lu, W.; Fan, D.; Luo, K. Genome-wide identification and characterization of the Populus WRKY transcription factor family and analysis of their expression in response to biotic and abiotic stresses. J. Exp. Bot. 2014, 65, 6629–6644. [Google Scholar] [CrossRef]
  47. Li, W.; Tian, Z.; Yu, D. WRKY13 acts in stem development in Arabidopsis thaliana. Plant Sci. 2015, 236, 205–213. [Google Scholar] [CrossRef] [PubMed]
  48. Ahn, J.H.; Kim, J.; Yoo, S.J.; Yoo, S.Y.; Roh, H.; Choi, J.H.; Choi, M.S.; Chung, K.S.; Han, E.J.; Hong, S.M. Isolation of 151 mutants that have developmental defects from T-DNA tagging. Plant Cell Physiol. 2007, 48, 169–178. [Google Scholar] [CrossRef] [Green Version]
  49. Ricachenevsky, F.K.; Punshon, T.; Lee, S.; Oliveira, B.H.N.; Trenz, T.S.; Maraschin, F.D.S.; Hindt, M.N.; Danku, J.; Salt, D.E.; Fett, J.P.; et al. Elemental Profiling of Rice FOX Lines Leads to Characterization of a New Zn Plasma Membrane Transporter, OsZIP7. Front. Plant Sci. 2018, 9, 865. [Google Scholar] [CrossRef] [Green Version]
  50. Breuer, C.; Kawamura, A.; Ichikawa, T.; Tominaga-Wada, R.; Wada, T.; Kondou, Y.; Muto, S.; Matsui, M.; Sugimoto, K. The Trihelix Transcription Factor GTL1 Regulates Ploidy-Dependent Cell Growth in the Arabidopsis Trichome. Plant Cell 2009, 21, 2307–2322. [Google Scholar] [CrossRef] [Green Version]
  51. Ye, Z.-H.; Zhong, R. Molecular control of wood formation in trees. J. Exp. Bot. 2015, 66, 4119–4131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Zhong, R.; Lee, C.; Ye, Z.-H. Functional Characterization of Poplar Wood-Associated NAC Domain Transcription Factors. Plant Physiol. 2010, 152, 1044–1055. [Google Scholar] [CrossRef] [Green Version]
  53. Lin, Y.C.; Li, W.; Sun, Y.H.; Kumari, S.; Wei, H.; Li, Q.; Tunlaya-Anukit, S.; Sederoff, R.R.; Chiang, V.L. SND1 Transcription Factor-Directed Quantitative Functional Hierarchical Genetic Regulatory Network in Wood Formation in Populus trichocarpa. Plant Cell 2013, 25, 4324–4341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wang, H.; Avci, U.; Nakashima, J.; Hahn, M.G.; Chen, F.; Dixon, R.A. Mutation of WRKY transcription factors initiates pith secondary wall formation and increases stem biomass in dicotyledonous plants. Proc. Natl. Acad. Sci. USA 2010, 107, 22338–22343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Ge, S.; Han, X.; Xu, X.; Shao, Y.; Zhu, Q.; Liu, Y.; Du, J.; Xu, J.; Zhang, S. WRKY15 Suppresses Tracheary Element Differentiation Upstream of VND7 During Xylem Formation. Plant Cell 2020, 32, 2307–2324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Didi, V.; Jackson, P.; Hejatko, J. Hormonal regulation of secondary cell wall formation. J. Exp. Bot. 2015, 66, 5015–5027. [Google Scholar] [CrossRef]
  57. Johnsson, C.; Jin, X.; Xue, W.; Dubreuil, C.; Lezhneva, L.; Fischer, U. The plant hormone auxin directs timing of xylem development by inhibition of secondary cell wall deposition through repression of secondary wall NAC-domain transcription factors. Physiol. Plant 2019, 165, 673–689. [Google Scholar] [CrossRef] [Green Version]
  58. Lee, K.H.; Du, Q.; Zhuo, C.; Qi, L.; Wang, H. LBD29-Involved Auxin Signaling Represses NAC Master Regulators and Fiber Wall Biosynthesis. Plant Physiol. 2019, 181, 595–608. [Google Scholar] [CrossRef] [Green Version]
  59. Soyano, T.; Thitamadee, S.; Machida, Y.; Chua, N.H. ASYMMETRIC LEAVES2-LIKE19/LATERAL ORGAN BOUNDARIES DOMAIN30 and ASL20/LBD18 regulate tracheary element differentiation in Arabidopsis. Plant Cell 2008, 20, 3359–3373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 22 03448 g001 550
Figure 1. Construction of a genome-wide Populus-Full-length cDNA Over-eXpressing (FOX)-Arabidopsis library. (A) Size distribution of full-length cDNA transgenes in entry vector. (B) Size distribution of full-length cDNA transgenes in expression vector. (C) Sequence lengths of 200 randomly selected clones from the expression clone library and 200 cDNAs from T1 transgenic plants. (D) Functional annotation and classification according to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway assignments of full-length cDNA transgenes in Arabidopsis thaliana.
Figure 1. Construction of a genome-wide Populus-Full-length cDNA Over-eXpressing (FOX)-Arabidopsis library. (A) Size distribution of full-length cDNA transgenes in entry vector. (B) Size distribution of full-length cDNA transgenes in expression vector. (C) Sequence lengths of 200 randomly selected clones from the expression clone library and 200 cDNAs from T1 transgenic plants. (D) Functional annotation and classification according to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway assignments of full-length cDNA transgenes in Arabidopsis thaliana.
Ijms 22 03448 g001
Ijms 22 03448 g002 550
Figure 2. Morphological phenotypes of transgenic plants from the FOX library under normal conditions. (A) Wild-type (WT). (B) Lodging. (C) Early senescence and shorter siliques. (D) Smaller florescence leaves and more branches. (E) Poor growth and smaller siliques. (F) Stiffer stems. (G) Bigger leaves and later bolting. (H) Fewer branches and enhanced anthocyanin accumulation. (I) Thick leaves and later bolting. (J) Glabrous leaves. (K) Curled leaves. (L) Wrinkled leaves. Scale bars: (A–L) 2 cm.
Figure 2. Morphological phenotypes of transgenic plants from the FOX library under normal conditions. (A) Wild-type (WT). (B) Lodging. (C) Early senescence and shorter siliques. (D) Smaller florescence leaves and more branches. (E) Poor growth and smaller siliques. (F) Stiffer stems. (G) Bigger leaves and later bolting. (H) Fewer branches and enhanced anthocyanin accumulation. (I) Thick leaves and later bolting. (J) Glabrous leaves. (K) Curled leaves. (L) Wrinkled leaves. Scale bars: (A–L) 2 cm.
Ijms 22 03448 g002
Ijms 22 03448 g003 550
Figure 3. Screening of abiotic stress-resistant FOX mutants. (A) Mock. (B) 150 mmol/L NaCl treatment for 6 days. (C) 10 μmol/L ABA treatment for 6 days. (D) 200 mmol/L mannitol treatment for 10 days.
Figure 3. Screening of abiotic stress-resistant FOX mutants. (A) Mock. (B) 150 mmol/L NaCl treatment for 6 days. (C) 10 μmol/L ABA treatment for 6 days. (D) 200 mmol/L mannitol treatment for 10 days.
Ijms 22 03448 g003
Ijms 22 03448 g004 550
Figure 4. PtoCPCa overexpression inhibits trichome development in A thaliana. (A) Morphological phenotypes of 3-week-old WT and PtoCPCa-OE lines. (B) Scanning electron microscopy (SEM) of trichomes of 3-week-old WT and PtoCPCa-OE lines. (C) Expression levels of PtoCPCa in WT and PtoCPCa-OE lines. (D) Expression levels of trichome development marker genes in WT and PtoCPCa-OE lines, including GL1 (AT3G27920), GL2 (AT1G79840), and TTG2 (AT2G37260). Actin7 was used as an internal control. Error bars represent the standard error (S.E.) of three independent biological replicates. Asterisks indicate a significant difference with respect to WT plants using the Student’s t-test (**, p < 0.01; ***, p < 0.001). Scale bars: (A) 1 cm; (B) 200 μm.
Figure 4. PtoCPCa overexpression inhibits trichome development in A thaliana. (A) Morphological phenotypes of 3-week-old WT and PtoCPCa-OE lines. (B) Scanning electron microscopy (SEM) of trichomes of 3-week-old WT and PtoCPCa-OE lines. (C) Expression levels of PtoCPCa in WT and PtoCPCa-OE lines. (D) Expression levels of trichome development marker genes in WT and PtoCPCa-OE lines, including GL1 (AT3G27920), GL2 (AT1G79840), and TTG2 (AT2G37260). Actin7 was used as an internal control. Error bars represent the standard error (S.E.) of three independent biological replicates. Asterisks indicate a significant difference with respect to WT plants using the Student’s t-test (**, p < 0.01; ***, p < 0.001). Scale bars: (A) 1 cm; (B) 200 μm.
Ijms 22 03448 g004
Ijms 22 03448 g005 550
Figure 5. PtoWRKY13 overexpression promotes secondary cell wall (SCW) deposition in transgenic A. thaliana inflorescence stems. (A) Cross-sectioning and staining with Toluidine Blue of basal stems of 6-week-old WT and PtoWRKY13-OE A. thaliana. (B) Stem diameter of WT and PtoWRKY13-OE lines. (C) Number of vessel bundles in WT and PtoWRKY13-OE lines. (D) Percentage of sclerenchyma cells in the stems of WT and PtoWRKY13-OE lines. (E) Sclerenchyma area in the stems of WT and PtoWRKY13-OE lines. (*, p < 0.05); ph, phloem; xf, xylem fiber; ve, vessel; if, interfascicular fiber. Scale bars: (A) 100 μm.
Figure 5. PtoWRKY13 overexpression promotes secondary cell wall (SCW) deposition in transgenic A. thaliana inflorescence stems. (A) Cross-sectioning and staining with Toluidine Blue of basal stems of 6-week-old WT and PtoWRKY13-OE A. thaliana. (B) Stem diameter of WT and PtoWRKY13-OE lines. (C) Number of vessel bundles in WT and PtoWRKY13-OE lines. (D) Percentage of sclerenchyma cells in the stems of WT and PtoWRKY13-OE lines. (E) Sclerenchyma area in the stems of WT and PtoWRKY13-OE lines. (*, p < 0.05); ph, phloem; xf, xylem fiber; ve, vessel; if, interfascicular fiber. Scale bars: (A) 100 μm.
Ijms 22 03448 g005
Ijms 22 03448 g006 550
Figure 6. Lignin biosynthetic genes and secondary wall-related genes are altered in PtoWRKY13-OE lines. (A) Expression levels of PtoWRKY13 in WT and PtoWRKY13-OE seedlings. (B) Expression levels of lignin pathway genes HCT (AT5G48930), PAL4 (AT3G10340), CAD 6(AT4G37070), and F5H (AT4G36220) in WT and transgenic A. thaliana. (C) Expression levels of secondary wall-related key transcription factors SND1 (AT1G32770), NST1 (AT2G46770), NST2 (AT3G61910), and C3H14 (AT1G66810) in WT and transgenic A. thaliana. Actin7 was used as an internal reference gene. Error bars represent the standard error (S.E.) of three independent biological replicates. Asterisks indicate a significant difference with respect to WT plants using the Student’s t-test (*, p < 0.05; **, p < 0.01; ***, p < 0.001).
Figure 6. Lignin biosynthetic genes and secondary wall-related genes are altered in PtoWRKY13-OE lines. (A) Expression levels of PtoWRKY13 in WT and PtoWRKY13-OE seedlings. (B) Expression levels of lignin pathway genes HCT (AT5G48930), PAL4 (AT3G10340), CAD 6(AT4G37070), and F5H (AT4G36220) in WT and transgenic A. thaliana. (C) Expression levels of secondary wall-related key transcription factors SND1 (AT1G32770), NST1 (AT2G46770), NST2 (AT3G61910), and C3H14 (AT1G66810) in WT and transgenic A. thaliana. Actin7 was used as an internal reference gene. Error bars represent the standard error (S.E.) of three independent biological replicates. Asterisks indicate a significant difference with respect to WT plants using the Student’s t-test (*, p < 0.05; **, p < 0.01; ***, p < 0.001).
Ijms 22 03448 g006
Table
Table 1. Phenotypes and functional annotation of mutant FOX cDNAs.
Table 1. Phenotypes and functional annotation of mutant FOX cDNAs.
LinePt Gene IDAt Gene IDQuery CoverPhenotypeAnnotation
N1Potri.002G212700.1AT1G65410.10.728Poor growth, thin stemNON-INTRINSIC ABC PROTEIN 11 (ATNAP11)
N3Potri.005G175900.1AT1G44835.20.715Poor growth, thin stemYbaK/aminoacyl-tRNA synthetase-associated domain-containing protein
N46Potri.008G075200.1AT3G10985.10.602Poor growth, senilismSENESCENCE ASSOCIATED GENE 20 (ATSAG20)
A33Potri.006G224000.1AT5G26180.10.7DwarfTRNA METHYLTRANSFERASE 4H (ATTRM4H)
N66Potri.006G130900.1AT5G60390.30.953DwarfGTP binding Elongation factor Tu family protein
N4Potri.003G126100.1AT4G11600.10.677Dwarf, poor growthGLUTATHIONE PEROXIDASE 6 (ATGPX6)
N9Potri.010G013400.1AT3G05540.10.792Dwarf, poor growthTRANSLATIONALLY CONTROLLED TUMOR PROTEIN 2 (ATTCTP2)
N15Potri.012G114200.1AT5G61670.20.772Dwarf, Poor growthProtein disulfide-isomerase/S-S rearrangase
N10Potri.001G123200.1AT4G16530.10.453Dwarf, poor growthDUF577
N19Potri.018G119600.1AT1G75400.10.354Dwarf, delayed floweringRING/U-box superfamily protein
N13Potri.011G045100.1AT2G20840.10.788Dwarf, delayed floweringSECRETORY CARRIER MEMBRANE PROTEIN 3 (ATSCAMP3)
N11Potri.008G145600.1AT1G17960.10.542Dwarf, fewer branchesThreonyl-tRNA synthetase
N6Potri.007G011200.1AT4G37870.10.873Dwarf, fewer branchesPHOSPHOENOLPYRUVATE CARBOXYKINASE 1 (ATPCK1)
N40Potri.004G233800.1AT5G61510.10.753Fewer branchesGroES-like zinc-binding alcohol dehydrogenase family protein
N65Potri.019G058000.1AT3G56360.10.542Fewer branchesHypothetical protein
A5Potri.008G178700.1AT1G24360.10.693Rapid growthNAD(P)-binding Rossmann-fold superfamily protein
A11Potri.016G091200.1AT5G03290.10.89Rapid growth, big leavesISOCITRATE DEHYDROGENASE V (ATIDHV)
A4Potri.012G047200.1AT1G73590.10.823Slender petiolesPIN-FORMED 1 (ATPIN1)
A1Potri.017G153500.1AT1G17860.10.665Slender petioles, dwarfTHALIANA KUNITZ TRYPSIN INHIBITOR 5 (ATKTI5)
A2Potri.001G162800.1AT1G17290.10.841Slender petioles, small leavesALANINE AMINOTRANSFERASE (ATAlaAT1)
A12Potri.019G054600.1AT1G24706.20.742Slender petiolesTHO/TREX complex
A7Potri.016G084800.1AT3G09640.10.86Small leaves, more branchesASCORBATE PEROXIDASE 2 (ATAPX2)
N21Potri.004G177500.1AT4G38510.40.98Small leaves, more branchesV-ATPASE B SUBUNIT 2 (ATVAB2)
A3Potri.009G042400.1AT3G54140.10.833Small leaves, more branchesPEPTIDE TRANSPORTER 1 (ATPTR1)
A6Potri.005G242500.1AT5G42150.10.775Small leavesGlutathione S-transferase family protein
A10Potri.009G133000.1AT4G38630.10.617Small leavesMULTIUBIQUITIN-CHAIN-BINDING PROTEIN 1 (ATMCB1)
N23Potri.008G056300.1AT3G55440.10.744Fewer cauline leavesCYTOSOLIC TRIOSE PHOSPHATE ISOMERASE (ATCTIMC)
N18Potri.015G029000.1AT1G48410.10.915Serrated leaves, delayed flowering, malformed siliquesARGONAUTE 1 (ATAGO1)
N94Potri.002G015100.1AT3G03190.10.642Small leaves, thin stemGLUTATHIONE S-TRANSFERASE F11 (ATGSTF11)
N96Potri.001G413500.1AT2G31670.10.649Thin stemStress responsive alpha-beta barrel domain protein.
N73Potri.004G015300.1AT4G21960.10.843Stiffer stemPeroxidase superfamily protein
A9Potri.008G161200.1AT4G14550.10.756Thicker stemINDOLE-3-ACETIC ACID INDUCIBLE 14 (ATIAA14)
N8Potri.005G146800.1AT5G67500.10.83LodgingVOLTAGE DEPENDENT ANION CHANNEL 2 (ATVDAC2)
A14Potri.004G081700.1AT1G27190.10.819Seed abortionBAK1-INTERACTING RECEPTOR-LIKE KINASE 3 (ATBIR3)
Query cover based on BLAST analysis against the Populus trichocarpa V3.0 database on phytozome. At—Arabidopsis thaliana, Pt—P. trichocarpa.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kong, L.; Li, Z.; Song, Q.; Li, X.; Luo, K. Construction of a Full-Length cDNA Over-Expressing Library to Identify Valuable Genes from Populus tomentosa. Int. J. Mol. Sci. 2021, 22, 3448. https://doi.org/10.3390/ijms22073448

AMA Style

Kong L, Li Z, Song Q, Li X, Luo K. Construction of a Full-Length cDNA Over-Expressing Library to Identify Valuable Genes from Populus tomentosa. International Journal of Molecular Sciences. 2021; 22(7):3448. https://doi.org/10.3390/ijms22073448

Chicago/Turabian Style

Kong, Lingfei, Zeyu Li, Qin Song, Xiaohong Li, and Keming Luo. 2021. "Construction of a Full-Length cDNA Over-Expressing Library to Identify Valuable Genes from Populus tomentosa" International Journal of Molecular Sciences 22, no. 7: 3448. https://doi.org/10.3390/ijms22073448

APA Style

Kong, L., Li, Z., Song, Q., Li, X., & Luo, K. (2021). Construction of a Full-Length cDNA Over-Expressing Library to Identify Valuable Genes from Populus tomentosa. International Journal of Molecular Sciences, 22(7), 3448. https://doi.org/10.3390/ijms22073448

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop