Genome-Wide Identification of the ERF Transcription Factor Family for Structure Analysis, Expression Pattern, and Response to Drought Stress in Populus alba × Populus glandulosa.

The Ethylene Responsive Factor (ERF) transcription factor family is important for regulating plant growth and stress responses. Although the expression patterns of ERF family members have been reported in many plant species, their role in Populus alba × Populus glandulosa, an important model plant for forest research, remains unclear. In this study, we identified 209 PagERF transcription factors by analyzing the P. alba × P. glandulosa genome. We analyzed their amino acid sequences, molecular weight, theoretical pI (Isoelectric point), instability index, aliphatic index, grand average of hydropathicity, and subcellular localization. Most PagERFs were predicted to localize in the nucleus, with only a few PagERFs localized in the cytoplasm and nucleus. Phylogenetic analysis divided the PagERF proteins into ten groups, Class I to X, with those belonging to the same group containing similar motifs. Cis-acting elements associated with plant hormones, abiotic stress responses, and MYB binding sites were analyzed in the promoters of PagERF genes. We used transcriptome data to analyze the expression patterns of PagERF genes in different tissues of P. alba × P. glandulosa, including axillary buds, young leaves, functional leaves, cambium, xylem, and roots, and the results indicated that PagERF genes are expressed in all tissues of P. alba × P. glandulosa, especially in roots. Quantitative verification results were consistent with transcriptome data. When P. alba × P. glandulosa seedlings were treated with 6% polyethylene glycol 6000 (PEG6000), the results of RT-qRCR showed that nine PagERF genes responded to drought stress in various tissues. This study provides a new perspective on the roles of PagERF family members in regulating plant growth and development, and responses to stress in P. alba × P. glandulosa. Our study provides a theoretical basis for ERF family research in the future.

Transcriptome Analysis of Populus euphratica under Salt Treatment and PeERF1 Gene Enhances Salt Tolerance in Transgenic Populus alba × Populus glandulosa.

Populus euphratica is mainly distributed in desert environments with dry and hot climate in summer and cold in winter. Compared with other poplars, P. euphratica is more resistant to salt stress. It is critical to investigate the transcriptome and molecular basis of salt tolerance in order to uncover stress-related genes. In this study, salt-tolerant treatment of P. euphratica resulted in an increase in osmo-regulatory substances and recovery of antioxidant enzymes. To improve the mining efficiency of candidate genes, the analysis combining both the transcriptome WGCNA and the former GWAS results was selected, and a range of key regulatory factors with salt resistance were found. The PeERF1 gene was highly connected in the turquoise modules with significant differences in salt stress traits, and the expression levels were significantly different in each treatment. For further functional verification of PeERF1, we obtained stable overexpression and dominant suppression transgenic lines by transforming into Populus alba × Populusglandulosa. The growth and physiological characteristics of the PeERF1 overexpressed plants were better than that of the wild type under salt stress. Transcriptome analysis of leaves of transgenic lines and WT revealed that highly enriched GO terms in DEGs were associated with stress responses, including abiotic stimuli responses, chemical responses, and oxidative stress responses. The result is helpful for in-depth analysis of the salt tolerance mechanism of poplar. This work provides important genes for poplar breeding with salt tolerance.

A genome wide transcriptional study of Populus alba x P. tremula var. glandulosa in response to nitrogen deficiency stress.

Poplar 84 K (Populus alba x P. tremula var. glandulosa) is a good resource for genetic engineering due to its rapid growth and wide adaptability, and it is also an excellent ornamental tree species. In this study, we used 84 K plantlets grown in the nitrogen-limited medium as experimental materials to explore the molecular mechanism in 84 K leaves under nitrogen deficiency. A total of 5,868 differentially expressed genes (DEGs) were identified using the transcriptional information from RNA-seq data. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment results revealed that the DEGs were mainly involved in energy metabolism and anthocyanin biosynthesis. We then identified differentially expressed transcription factors (TFs) and constructed TF centered gene co-expression networks for chlorophyll and anthocyanin biosynthesis pathway genes. Twenty potential regulators were finally identified. We speculated the transcription factors that control the pigmentation in leaves with the MYB-bHLH-WD40 (MBW) pigment regulatory model. Such identification will clarify the genetic basis of the secondary metabolism in 84 K, and being a source of candidate genes for future plant genetic engineering. Our work broadens the researchers' understanding of the regulation of anthocyanin synthesis in trees and provides new perspectives for ornamental 84 K poplar breeding. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12298-021-01012-3.

Overexpression Populus d-Type Cyclin Gene PsnCYCD1;1 Influences Cell Division and Produces Curved Leaf in Arabidopsis thaliana.

d-type cyclins (CYCDs) are a special class of cyclins and play extremely important roles in plant growth and development. In the plant kingdom, most of the existing studies on CYCDs have been done on herbaceous plants, with few on perennial woody plants. Here, we identified a Populus d-type cyclin gene, PsnCYCD1;1, which is mainly transcribed in leaf buds and stems. The promoter of PsnCYCD1;1 activated GUS gene expression and transgenic Arabidopsis lines were strongly GUS stained in whole seedlings and mature anthers. Moreover, subcellular localization analysis showed the fluorescence signal of PsnCYCD1;1-GFP fusion protein is present in the nucleus. Furthermore, overexpression of the PsnCYCD1;1 gene in Arabidopsis can promote cell division and lead to small cell generation and cytokinin response, resulting in curved leaves and twisted inflorescence stems. Moreover, the transcriptional levels of endogenous genes, such as ASs, KNATs, EXP10, and PHB, were upregulated by PsnCYCD1;1. Together, our results indicated that PsnCYCD1;1 participates in cell division by cytokinin response, providing new information on controlling plant architecture in woody plants.

Genome sequence and evolution of Betula platyphylla.

Betula L. (birch) is a pioneer hardwood tree species with ecological, economic, and evolutionary importance in the Northern Hemisphere. We sequenced the Betula platyphylla genome and assembled the sequences into 14 chromosomes. The Betula genome lacks evidence of recent whole-genome duplication and has the same paleoploidy level as Vitis vinifera and Prunus mume. Phylogenetic analysis of lignin pathway genes coupled with tissue-specific expression patterns provided clues for understanding the formation of higher ratios of syringyl to guaiacyl lignin observed in Betula species. Our transcriptome analysis of leaf tissues under a time-series cold stress experiment revealed the presence of the MEKK1-MKK2-MPK4 cascade and six additional mitogen-activated protein kinases that can be linked to a gene regulatory network involving many transcription factors and cold tolerance genes. Our genomic and transcriptome analyses provide insight into the structures, features, and evolution of the B. platyphylla genome. The chromosome-level genome and gene resources of B. platyphylla obtained in this study will facilitate the identification of important and essential genes governing important traits of trees and genetic improvement of B. platyphylla.

Transcriptomic analyses of Pinus koraiensis under different cold stresses.

BACKGROUND: Pinus koraiensis is an evergreen tree species with strong cold resistance. However, the transcriptomic patterns in response to cold stress are poorly understood for P. koraiensis. In this study, global transcriptome profiles were generated for P. koraiensis under cold stress (- 20 °C) over time by high-throughput sequencing. RESULTS: More than 763 million clean reads were produced, which assembled into a nonredundant data set of 123,445 unigenes. Among them, 38,905 unigenes had homology with known genes, 18,239 were assigned to 54 gene ontology (GO) categories and 18,909 were assigned to 25 clusters of orthologous groups (COG) categories. Comparison of transcriptomes of P. koraiensis seedlings grown at room temperature (20 °C) and low temperature (- 20 °C) revealed 9842 differential expressed genes (DEGs) in the 6 h sample, 9250 in the 24 h sample, and 9697 in the 48 h sample. The number of DEGs in the pairwise comparisons of 6 h, 24 h and 48 h was relatively small. The accuracy of the RNA-seq was validated by analyzing the expression patterns of 12 DEGs by quantitative real-time PCR (qRT-PCR). In this study, 34 DEGs (22 upregulated and 12 downregulated) were involved in the perception and transmission of cold signals, 96 DEGs (41 upregulated and 55 downregulated) encoding 8 transcription factors that regulated cold-related genes expression, and 27 DEGs (17 upregulated and 10 downregulated) were involved in antioxidant mechanisms in response to cold stress. Among them, the expression levels of c63631_g1 (annexin D1), c65620_g1 (alpha-amylase isozyme 3C), c61970_g1 (calcium-binding protein KIC), c51736_g1 (ABA), c58408_g1 (DREB3), c66599_g1 (DREB3), c67548_g2 (SOD), c55044_g1 (CAT), c71938_g2 (CAT) and c11358_g1 (GPX) first increased significantly and then decreased significantly with the extension of stress time. CONCLUSIONS: A large number of DEGs were identified in P. koraiensis under cold stress, especially the DEGs involved in the perception and transmission of cold signals, the DEGs encoding TFs related to cold regulation and the DEGs removing ROS in antioxidation mechanisms. The transcriptome and digital expression profiling of P. koraiensis could facilitate the understanding of the molecular control mechanism related to cold responses and provide the basis for the molecular breeding of conifers.

A novel synthetic-genetic-array-based yeast one-hybrid system for high discovery rate and short processing time.

Eukaryotic gene expression is often tightly regulated by interactions between transcription factors (TFs) and their DNA cis targets. Yeast one-hybrid (Y1H) is one of the most extensively used methods to discover these interactions. We developed a high-throughput meiosis-directed yeast one-hybrid system using the Magic Markers of the synthetic genetic array analysis. The system has a transcription factor-DNA interaction discovery rate twice as high as the conventional diploid-mating approach and a processing time nearly one-tenth of the haploid-transformation method. The system also offers the highest accuracy in identifying TF-DNA interactions that can be authenticated in vivo by chromatin immunoprecipitation. With these unique features, this meiosis-directed Y1H system is particularly suited for constructing novel and comprehensive genome-scale gene regulatory networks for various organisms.

Complete chloroplast genome sequence of Betula platyphylla: gene organization, RNA editing, and comparative and phylogenetic analyses.

BACKGROUND: Betula platyphylla is a common tree species in northern China that has high economic and medicinal value. Our laboratory has been devoted to genome research on B. platyphylla for approximately 10 years. As primary organelle genomes, the complete genome sequences of chloroplasts are important to study the divergence of species, RNA editing and phylogeny. In this study, we sequenced and analyzed the complete chloroplast (cp) genome sequence of B. platyphylla. RESULTS: The complete cp genome of B. platyphylla was 160,518 bp in length, which included a pair of inverted repeats (IRs) of 26,056 bp that separated a large single copy (LSC) region of 89,397 bp and a small single copy (SSC) region of 19,009 bp. The annotation contained a total of 129 genes, including 84 protein-coding genes, 37 tRNA genes and 8 rRNA genes. There were 3 genes using alternative initiation codons. Comparative genomics showed that the sequence of the Fagales species cp genome was relatively conserved, but there were still some high variation regions that could be used as molecular markers. The IR expansion event of B. platyphylla resulted in larger cp genomes and rps19 pseudogene formation. The simple sequence repeat (SSR) analysis showed that there were 105 SSRs in the cp genome of B. platyphylla. RNA editing sites recognition indicated that at least 80 RNA editing events occurred in the cp genome. Most of the substitutions were C to U, while a small proportion of them were not. In particular, three editing loci on the rRNA were converted to more than two other bases that had never been reported. For synonymous conversion, most of them increased the relative synonymous codon usage (RSCU) value of the codons. The phylogenetic analysis suggested that B. platyphylla had a closer evolutionary relationship with B. pendula than B. nana. CONCLUSIONS: In this study, we not only obtained and annotated the complete cp genome sequence of B. platyphylla, but we also identified new RNA editing sites and predicted the phylogenetic relationships among Fagales species. These findings will facilitate genomic, genetic engineering and phylogenetic studies of this important species.

Variation Analysis of Physiological Traits in Betula platyphylla Overexpressing TaLEA-ThbZIP Gene under Salt Stress.

The aim of this study was to determine whether transgenic birch (Betula platyphylla) ectopic overexpressing a late embryogenesis abundant (LEA) gene and a basic leucine zipper (bZIP) gene from the salt-tolerant genus Tamarix (salt cedar) show increased tolerance to salt (NaCl) stress. Co-transfer of TaLEA and ThbZIP in birch under the control of two independent CaMV 35S promoters significantly enhanced salt stress. PCR and northern blot analyses indicated that the two genes were ectopically overexpressed in several dual-gene transgenic birch lines. We compared the effects of salt stress among three transgenic birch lines (L-4, L-5, and L-8) and wild type (WT). In all lines, the net photosynthesis values were higher before salt stress treatment than afterwards. After the salt stress treatment, the transgenic lines L-4 and L-8 showed higher values for photosynthetic traits, chlorophyll fluorescence, peroxidase and superoxide dismutase activities, and lower malondialdehyde and Na+ contents, compared with those in WT and L-5. These different responses to salt stress suggested that the transcriptional level of the TaLEA and ThbZIP genes differed among the transgenic lines, resulting in a variety of genetic and phenotypic effects. The results of this research can provide a theoretical basis for the genetic engineering of salt-tolerant trees.

Phosphorylation is an on/off switch for 5-hydroxyconiferaldehyde O-methyltransferase activity in poplar monolignol biosynthesis.

Although phosphorylation has long been known to be an important regulatory modification of proteins, no unequivocal evidence has been presented to show functional control by phosphorylation for the plant monolignol biosynthetic pathway. Here, we present the discovery of phosphorylation-mediated on/off regulation of enzyme activity for 5-hydroxyconiferaldehyde O-methyltransferase 2 (PtrAldOMT2), an enzyme central to monolignol biosynthesis for lignification in stem-differentiating xylem (SDX) of Populus trichocarpa. Phosphorylation turned off the PtrAldOMT2 activity, as demonstrated in vitro by using purified phosphorylated and unphosphorylated recombinant PtrAldOMT2. Protein extracts of P. trichocarpa SDX, which contains endogenous kinases, also phosphorylated recombinant PtrAldOMT2 and turned off the recombinant protein activity. Similarly, ATP/Mn(2+)-activated phosphorylation of SDX protein extracts reduced the endogenous SDX PtrAldOMT2 activity by ∼ 60%, and dephosphorylation fully restored the activity. Global shotgun proteomic analysis of phosphopeptide-enriched P. trichocarpa SDX protein fractions identified PtrAldOMT2 monophosphorylation at Ser(123) or Ser(125) in vivo. Phosphorylation-site mutagenesis verified the PtrAldOMT2 phosphorylation at Ser(123) or Ser(125) and confirmed the functional importance of these phosphorylation sites for O-methyltransferase activity. The PtrAldOMT2 Ser(123) phosphorylation site is conserved across 93% of AldOMTs from 46 diverse plant species, and 98% of the AldOMTs have either Ser(123) or Ser(125). PtrAldOMT2 is a homodimeric cytosolic enzyme expressed more abundantly in syringyl lignin-rich fiber cells than in guaiacyl lignin-rich vessel cells. The reversible phosphorylation of PtrAldOMT2 is likely to have an important role in regulating syringyl monolignol biosynthesis of P. trichocarpa.

Overexpression of two PsnAP1 genes from Populus simonii × P. nigra causes early flowering in transgenic tobacco and Arabidopsis.

In Arabidopsis, AP1 is a floral meristem identity gene and plays an important role in floral organ development. In this study, PsnAP1-1 and PsnAP1-2 were isolated from the male reproductive buds of poplar (Populus simonii × P. nigra), which are the orthologs of AP1 in Arabidopsis, by sequence analysis. Northern blot and qRT-PCR analysis showed that PsnAP1-1 and PsnAP1-2 exhibited high expression level in early inflorescence development of poplar. Subcellular localization showed the PsnAP1-1 and PsnAP1-2 proteins are localized in the nucleus. Overexpression of PsnAP1-1 and PsnAP1-2 in tobacco under the control of a CaMV 35S promoter significantly enhanced early flowering. These transgenic plants also showed much earlier stem initiation and higher rates of photosynthesis than did wild-type tobacco. qRT-PCR analysis further indicated that overexpression of PsnAP1-1 and PsnAP1-2 resulted in up-regulation of genes related to flowering, such as NtMADS4, NtMADS5 and NtMADS11. Overexpression of PsnAP1-1 and PsnAP1-2 in Arabidopsis also induced early flowering, but did not complement the ap1-10 floral morphology to any noticeable extent. This study indicates that PsnAP1-1 and PsnAP1-2 play a role in floral transition of poplar.

Plant biotechnology for lignocellulosic biofuel production.

Lignocelluloses from plant cell walls are attractive resources for sustainable biofuel production. However, conversion of lignocellulose to biofuel is more expensive than other current technologies, due to the costs of chemical pretreatment and enzyme hydrolysis for cell wall deconstruction. Recalcitrance of cell walls to deconstruction has been reduced in many plant species by modifying plant cell walls through biotechnology. These results have been achieved by reducing lignin content and altering its composition and structure. Reduction of recalcitrance has also been achieved by manipulating hemicellulose biosynthesis and by overexpression of bacterial enzymes in plants to disrupt linkages in the lignin-carbohydrate complexes. These modified plants often have improved saccharification yield and higher ethanol production. Cell wall-degrading (CWD) enzymes from bacteria and fungi have been expressed at high levels in plants to increase the efficiency of saccharification compared with exogenous addition of cellulolytic enzymes. In planta expression of heat-stable CWD enzymes from bacterial thermophiles has made autohydrolysis possible. Transgenic plants can be engineered to reduce recalcitrance without any yield penalty, indicating that successful cell wall modification can be achieved without impacting cell wall integrity or plant development. A more complete understanding of cell wall formation and structure should greatly improve lignocellulosic feedstocks and reduce the cost of biofuel production.

A robust chromatin immunoprecipitation protocol for studying transcription factor-DNA interactions and histone modifications in wood-forming tissue.

Woody cells and tissues are recalcitrant to standard chromatin immunoprecipitation (ChIP) procedures. However, we recently successfully implemented ChIP in wood-forming tissue of the model woody plant Populus trichocarpa. Here we provide the detailed ChIP protocol optimized for wood-forming tissue that we used in those studies. By using stem-differentiating xylem (SDX; a wood-forming tissue), we identified all steps that were ineffective in standard ChIP protocols and systematically modified them to develop and optimize a robust ChIP protocol. The protocol includes tissue collection, cross-linking, nuclear isolation, chromatin extraction, DNA fragmentation, immunoprecipitation, DNA purification and sequence analysis. The protocol takes 2.5 d to complete and allows a robust 8-10-fold enrichment of transcription factor (TF)-bound genomic fragments (~150 ng/g of SDX) over nonspecific DNAs. The enriched DNAs are of high quality and can be used for subsequent PCR and DNA-seq analyses. We used this protocol to identify genome-wide specific TF-DNA interactions during wood formation and histone modifications associated with regulation of wood formation. Our protocol, which may be suitable for many tissue types, is so far the only working ChIP system for wood-forming tissue.

A simple improved-throughput xylem protoplast system for studying wood formation.

Isolated protoplasts serve as a transient expression system that is highly representative of stable transgenics in terms of transcriptome responses. They can also be used as a cellular system to study gene transactivation and nucleocytoplasmic protein trafficking. They are particularly useful for systems studies in which stable transgenics and mutants are unavailable. We present a protocol for the isolation and transfection of protoplasts from wood-forming tissue, the stem-differentiating xylem (SDX), in the model woody plant Populus trichocarpa. The method involves tissue preparation, digestion of SDX cell walls, protoplast isolation and DNA transfection. Our approach is markedly faster and provides better yields than previous protocols; small (milligrams)- to large (20 g)-scale SDX preparations can be achieved in ~60 s, with isolation of protoplasts and their subsequent transfection taking ~50 min. Up to ten different samples can be processed simultaneously in this time scale. Our protocol gives a high yield (~2.5 × 10(7) protoplasts per g of SDX) of protoplasts sharing 96% transcriptome identity with intact SDX.

Overexpression of a MADS-box gene from birch (Betula platyphylla) promotes flowering and enhances chloroplast development in transgenic tobacco.

In this study, a MADS-box gene (BpMADS), which is an ortholog of AP1 from Arabidopsis, was isolated from birch (Betula platyphylla). Transgenic Arabidopsis containing a BpMADS promoter::GUS construct was produced, which exhibited strong GUS staining in sepal tissues. Ectopic expression of BpMADS significantly enhanced the flowering of tobacco (35S::BpMADS). In addition, the chloroplasts of transgenic tobacco exhibited much higher growth and division rates, as well rates of photosynthesis, than wild-type. A grafting experiment demonstrated that the flowering time of the scion was not affected by stock that overexpressed BpMADS. In addition, the overexpression of BpMADS resulted in the upregulation of some flowering-related genes in tobacco.

A novel zinc-finger-like gene from Tamarix hispida is involved in salt and osmotic tolerance.

In the present study, a zinc-finger-like cDNA (ThZFL) was cloned from the Tamarix hispida. Northern blot analysis showed that the expression of ThZFL can be induced by salt, osmotic stress and ABA treatment. Overexpression of the ThZFL confers salt and osmotic stress tolerance in both yeast Saccharomyces cerevisiae and tobacco. Furthermore, MDA levels in ThZFL transformed tobacco were significantly decreased compared with control plants under salt and osmotic stress, suggesting ThZFL may confer stress tolerance by decreasing membrane lipid peroxidation. Subcellular localization analysis showed the ThZFL protein is localized in the cell wall. Our results indicated the ThZFL gene is an excellent candidate for genetic engineering to improve salt and osmotic tolerance in agricultural plants.

Enhanced salt tolerance of transgenic poplar plants expressing a manganese superoxide dismutase from Tamarix androssowii.

Superoxide dismutases (SODs) play important role in stress tolerance of plants. In this study, an MnSOD gene (TaMnSOD) from Tamarix androssowii, under the control of the CaMV35S promoter, was introduced into poplar (Populus davidiana x P. bolleana). The physiological parameters, including SOD activity, malondialdehyde (MDA) content, relative electrical conductivity (REC) and relative weight gain, of transgenic lines and wild type (WT) plants, were measured and compared. The results showed that SOD activity was enhanced in transgenic plants, and the MDA content and REC were significantly decreased compared to WT plants when exposed to NaCl stress. In addition, the relative weight gains of the transgenic plants were 8- to 23-fold of those observed for WT plants after NaCl stress for 30 days. The data showed that the SOD activities that increased in transgenic lines are 1.3-4-folds of that increased in the WT plant when exposed to NaCl stress. Our analysis showed that increases in SOD activities as low as 0.15-fold can also significantly enhance salt tolerance in transgenic plants, suggesting an important role of increased SOD activity in plant salt tolerance

Cloning and morphological properties of Nicgl;CYCD3;1 gene in genetic tumors from interspecific hybrid of N. langsdorffii and N. glauca.

Plant genetic tumors represent neoplastic growths, which arise spontaneously in hybrid plants without apparent external induction. To understand the molecular nature of unregulated cell proliferation, a cyclin D cDNA clone encoding a cyclin D of 1104bp was isolated from a genetic tumor and designated Nicgl;CYCD3;1 gene. DNA gel blot analysis suggested that there are two copies of Nicgl;CYCD3;1 in the genetic tumors. Northern analysis showed that this gene had the highest expression level in genetic tumor compared to Nicotiana glauca, N. langsdorffii and hybrid plants. Plant morphology of hybrid plant was an intermediate between N. glauca and N. langsdorffii and was altered in the genetic tumors. The cell cycle distribution in N. glauca was G0/G1, 90.59; S, 0.60; G2/M, 8.81; in N. langsdorffii it was G 0/G1, 86.22; S, 6.90; G2/M, 6.88; in hybrid plants it was G 0/G1, 96.40; S, 1.79; G2/M, 1.81; and in genetic tumors G 0/G1, 74.70; S, 2.35; G2/M, 22.94. These data provide new insights into the molecular mechanisms underlying genetic tumor formation from interspecific hybrid between N. langsdorffii and N. glauca.