Competing Endogenous RNA Networks Underlying Anatomical and Physiological Characteristics of Poplar Wood in Acclimation to Low Nitrogen Availability

How lignin biosynthesis responds to nitrogen in plants: a scoping review

Nitrogen (N) plays a critical role in the functioning of key amino acids and synthetic enzymes responsible for the various stages of lignin biosynthesis. However, the precise mechanisms through which N influences lignin biosynthesis have not been fully elucidated. This scoping review explores how lignin biosynthesis responds to N in plants. A systematic search of the literature in several databases was conducted using relevant keywords. Only 44 of the 1842 selected studies contained a range of plant species, experimental conditions, and research approaches. Lignin content, structure, and biosynthetic pathways in response to N are discussed, and possible response mechanisms of lignin under low N are proposed. Among the selected studies, 64.52% of the studies reter to lignin content found a negative correlation between N availability and lignin content. Usually, high N decreases the lignin content, delays cell lignification, increases p-hydroxyphenyl propane (H) monomer content, and regulates lignin synthesis through the expression of key genes (PAL, 4CL, CCR, CAD, COMT, LAC, and POD) encoding miRNAs and transcription factors (e.g., MYB, bHLH). N deficiency enhances lignin synthesis through the accumulation of phenylpropanoids, phenolics, and soluble carbohydrates, and indirect changes in phytohormones, secondary metabolites, etc. This review provides new insights and important references for future studies on the regulation of lignin biosynthesis.

Integrated small RNA, transcriptome and physiological approaches provide insight into Taxodium hybrid 'Zhongshanshan' roots in acclimation to prolonged flooding

Although Taxodium hybrid 'Zhongshanshan' 406 (Taxodium mucronatum Tenore × Taxodium distichum; Taxodium 406) is an extremely flooding-tolerant woody plant, the physiological and molecular mechanisms underlying acclimation of its roots to long-term flooding remain largely unknown. Thus, we exposed saplings of Taxodium 406 to either non-flooding (control) or flooding for 2 months. Flooding resulted in reduced root biomass, which is in line with lower concentrations of citrate, α-ketoglutaric acid, fumaric acid, malic acid and adenosine triphosphate (ATP) in Taxodium 406 roots. Flooding led to elevated activities of pyruvate decarboxylase, alcohol dehydrogenase and lactate dehydrogenase, which is consistent with higher lactate concentration in the roots of Taxodium 406. Flooding brought about stimulated activities of superoxide dismutase and catalase and elevated reduced glutathione (GSH) concentration and GSH/oxidized glutathione, which is in agreement with reduced concentrations of O2- and H2O2 in Taxodium 406 roots. The levels of starch, soluble protein, indole-3-acetic acid, gibberellin A4 and jasmonate were decreased, whereas the concentrations of glucose, total non-structural carbohydrates, most amino acids and 1-aminocyclopropane-1-carboxylate (ACC) were improved in the roots of flooding-treated Taxodium 406. Underlying these changes in growth and physiological characteristics, 12,420 mRNAs and 42 miRNAs were significantly differentially expressed, and 886 miRNA-mRNA pairs were identified in the roots of flooding-exposed Taxodium 406. For instance, 1-aminocyclopropane-1-carboxylate synthase 8 (ACS8) was a target of Th-miR162-3p and 1-aminocyclopropane-1-carboxylate oxidase 4 (ACO4) was a target of Th-miR166i, and the downregulation of Th-miR162-3p and Th-miR166i results in the upregulation of ACS8 and ACO4, probably bringing about higher ACC content in flooding-treated roots. Overall, these results indicate that differentially expressed mRNA and miRNAs are involved in regulating tricarboxylic acid cycle, ATP production, fermentation, and metabolism of carbohydrates, amino acids and phytohormones, as well as reactive oxygen species detoxification of Taxodium 406 roots. These processes play pivotal roles in acclimation to flooding stress. These results will improve our understanding of the molecular and physiological bases underlying woody plant flooding acclimation and provide valuable insights into breeding-flooding tolerant trees.

Network Analysis of Metabolome and Transcriptome Revealed Regulation of Different Nitrogen Concentrations on Hybrid Poplar Cambium Development

Secondary development is a key biological characteristic of woody plants and the basis of wood formation. Exogenous nitrogen can affect the secondary growth of poplar, and some regulatory mechanisms have been found in the secondary xylem. However, the effect of nitrogen on cambium has not been reported. Herein, we investigated the effects of different nitrogen concentrations on cambium development using combined transcriptome and metabolome analysis. The results show that, compared with 1 mM NH4NO3 (M), the layers of hybrid poplar cambium cells decreased under the 0.15 mM NH4NO3 (L) and 0.3 mM NH4NO3 (LM) treatments. However, there was no difference in the layers of hybrid poplar cambium cells under the 3 mM NH4NO3 (HM) and 5 mM NH4NO3 (H) treatments. Totals of 2365, 824, 649 and 398 DEGs were identified in the M versus (vs.) L, M vs. LM, M vs. HM and M vs. H groups, respectively. Expression profile analysis of the DEGs showed that exogenous nitrogen affected the gene expression involved in plant hormone signal transduction, phenylpropanoid biosynthesis, the starch and sucrose metabolism pathway and the ubiquitin-mediated proteolysis pathway. In M vs. L, M vs. LM, M vs. HM and M vs. H, differential metabolites were enriched in flavonoids, lignans, coumarins and saccharides. The combined analysis of the transcriptome and metabolome showed that some genes and metabolites in plant hormone signal transduction, phenylpropanoid biosynthesis and starch and sucrose metabolism pathways may be involved in nitrogen regulation in cambium development, whose functions need to be verified. In this study, from the point of view that nitrogen influences cambium development to regulate wood formation, the network analysis of the transcriptome and metabolomics of cambium under different nitrogen supply levels was studied for the first time, revealing the potential regulatory and metabolic mechanisms involved in this process and providing new insights into the effects of nitrogen on wood development.

Physiological and transcriptional regulation in Taxodium hybrid 'Zhongshanshan' leaves in acclimation to prolonged partial submergence

MAIN CONCLUSION

Taxodium 703 leaves activate fermentation, amino acids metabolism and ROS detoxification, and reduce TCA cycle and ABA biosynthesis in acclimation to prolonged partial submergence stress. Taxodium hybrid 'Zhongshanshan 703' (T. mucronatum × T. distichum; Taxodium 703) is a highly flooding-tolerant woody plant. To investigate the physiological and transcriptional regulatory mechanisms underlying its leaves in acclimation to long-term flooding, we exposed cuttings of Taxodium 703 to either non-flooding (control) or partial submergence for 2 months. The leaf tissues above (AL) and below (BL) flooding-water were separately harvested. Partial submergence decreased concentrations of chlorophyll (a + b) and dehydroascorbate (DHA) and lactate dehydrogenase (LDH) activity in AL, and reduced biomass, concentrations of succinic acid, fumaric acid and malic acid, and transcript levels of genes involved in tricarboxylic acid (TCA) cycle in BL. Under partial submergence, concentrations of starch, malondialdehyde and abscisic acid (ABA) decreased, and also mRNA levels of nine-cis-epoxycarotenoid dioxygenases that are involved in ABA biosynthesis in AL and BL of Taxodium 703. Partial submergence increased O2- content in AL, and improved concentrations of pyruvate and soluble sugars and activities of LDH and peroxidase in BL. In addition, partial submergence increased concentrations of ethanol, lactate, alanine, γ-aminobutyric acid, total amino acids and ascorbic acid (ASA), and ASA/DHA, activities of alcohol dehydrogenases (ADH) and ascorbate peroxidase, as well as transcript levels of ADH1A, ADH1B and genes involved in alanine biosynthesis and starch degradation in AL and BL of Taxodium 703. Overall, these results suggest that Taxodium 703 leaves activate fermentation, amino acids metabolism and reactive oxygen species detoxification, and maintain a steady supply of sugars, and reduce TCA cycle and ABA biosynthesis in acclimation to prolonged partial submergence stress.

Long non-coding RNA-mediated competing endogenous RNA regulatory network during flower development and color formation in Melastoma candidum

M. candidum, an evergreen shrubby flower known for its superior adaptation ability in South China, has gained increased attention in garden applications. However, scant attention has been paid to its flower development and color formation process at the non-coding RNA level. To fill this gap, we conducted a comprehensive analysis based on long non-coding RNA sequencing (lncRNA-seq), RNA-seq, small RNA sequencing (sRNA-seq), and widely targeted metabolome detection of three different flower developmental stages of M. candidum. After differentially expressed lncRNAs (DElncRNAs), differentially expressed mRNAs (DEmRNAs), differentially expressed microRNAs (DEmiRNAs), and differentially synthesized metabolites (DSmets) analyses between the different flower developmental stages, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were conducted to identify some key genes and metabolites in flavonoid, flavone, anthocyanin, carotenoid, and alkaloid-related GO terms and biosynthetic pathways. Three direct-acting models, including antisense-acting, cis-acting, and trans-acting between lncRNAs and mRNAs, were detected to illustrate the direct function of lncRNAs on target genes during flower development and color formation. Based on the competitive endogenous RNA (ceRNA) regulatory theory, we constructed a lncRNA-mediated regulatory network composed of DElncRNAs, DEmiRNAs, DEmRNAs, and DSmets to elucidate the indirect role of lncRNAs in the flower development and color formation of M. candidum. By utilizing correlation analyses between DERNAs and DSmets within the ceRNA regulatory network, alongside verification trials of the ceRNA regulatory mechanism, the study successfully illustrated the significance of lncRNAs in flower development and color formation process. This research provides a foundation for improving and regulating flower color at the lncRNA level in M. candidum, and sheds light on the potential applications of non-coding RNA in studies of flower development.

Long Non-Coding RNA lncWOX11a Suppresses Adventitious Root Formation of Poplar by Regulating the Expression of PeWOX11a

Long non-coding RNAs (lncRNAs), a class of poorly conserved transcripts without protein-encoding ability, are widely involved in plant organogenesis and stress responses by mediating the transmission and expression of genetic information at the transcriptional, posttranscriptional, and epigenetic levels. Here, we cloned and characterized a novel lncRNA molecule through sequence alignment, Sanger sequencing, transient expression in protoplasts, and genetic transformation in poplar. lncWOX11a is a 215 bp transcript located on poplar chromosome 13, ~50 kbp upstream of PeWOX11a on the reverse strand, and the lncRNA may fold into a series of complex stem–loop structures. Despite the small open reading frame (sORF) of 51 bp within lncWOX11a, bioinformatics analysis and protoplast transfection revealed that lncWOX11a has no protein-coding ability. The overexpression of lncWOX11a led to a decrease in the quantity of adventitious roots on the cuttings of transgenic poplars. Further, cis-regulatory module prediction and CRISPR/Cas9 knockout experiments with poplar protoplasts demonstrated that lncWOX11a acts as a negative regulator of adventitious rooting by downregulating the WUSCHEL-related homeobox gene WOX11, which is supposed to activate adventitious root development in plants. Collectively, our findings imply that lncWOX11a is essential for modulating the formation and development of adventitious roots.

Transcriptome profiling reveals the roles of pigment formation mechanisms in yellow Paeonia delavayi flowers

The yellow colour of ornamental varieties of tree peony originated from Paeonia delavayi. However, but P. delavayi and Paeonia suffruticosa belong to different subgroups, so hybridization is difficult and results in a long breeding cycle. However, no comprehensive transcriptomic profiling has focused on the colour formation mechanisms of yellow tree peony petals. Analysing the colour formation mechanism of yellow petals in P. delavayi is very important for directional molecular breeding. In this study, the transcriptional map of yellow pigment development in petals was used to analyse the mechanism of petal colour formation. We analysed the genes related to the metabolism of flavonoids and carotenoids and the transcription factors (TFs) involved in P. delavayi var. lutea (pure yellow individual) yellow pigment development using transcriptome sequence profiling. Transcriptome sequence profiles revealed three and four differentially expressed transcripts (DETs) involved in flavonoid biosynthesis and carotenoid biosynthesis, respectively. An analysis of DETs in the flavonoid pathway showed that chalcone synthase (CHS) and chalcone 2´-glucosyltransferases (THC2'GT) act in synergy to synthesize isosalipurposide (ISP). CHS and flavonol synthase (FLS) synergistically synthesize quercetin and kaempferol. DEG analysis of the carotenoid pathway revealed that phytoene synthase (PSY), carotenoid isomerase (CRTISO) and β-carotene hydroxylases (CHYB) play a key role in regulating lutein formation, and carotenoid cleavage dioxygenase (CCD) plays an important role in the degradation of carotenoids. These two pathways may be regulated by TF families such as bHLH, ARF, and MYB. The results of the transient overexpression of genes showed that CHS and CHI are regulated by PdMYB2. In this study, the molecular mechanism of ISP synthesis was analysed in depth, and the complete metabolic pathway of carotenoids in Paeonia L. was reported for the first time. By studying the formation mechanism of yellow pigment in P. delavayi petals, a breeding strategy for improving flavonol and carotenoid contents and reducing anthocyanin synthesis by genetic engineering was suggested.

Identification and Functional Prediction of Poplar Root circRNAs Involved in Treatment With Different Forms of Nitrogen

Circular RNAs (circRNAs) are a class of noncoding RNA molecules with ring structures formed by covalent bonds and are commonly present in organisms, playing an important regulatory role in plant growth and development. However, the mechanism of circRNAs in poplar root responses to different forms of nitrogen (N) is still unclear. In this study, high-throughput sequencing was used to identify and predict the function of circRNAs in the roots of poplar exposed to three N forms [1 mM NO₃ ^- (T1), 0.5 mM NH₄NO₃ (T2, control) and 1 mM NH₄ ⁺ (T3)]. A total of 2,193 circRNAs were identified, and 37, 24 and 45 differentially expressed circRNAs (DECs) were screened in the T1-T2, T3-T2 and T1-T3 comparisons, respectively. In addition, 30 DECs could act as miRNA sponges, and several of them could bind miRNA family members that play key roles in response to different N forms, indicating their important functions in response to N and plant growth and development. Furthermore, we generated a competing endogenous RNA (ceRNA) regulatory network in poplar roots treated with three N forms. DECs could participate in responses to N in poplar roots through the ceRNA regulatory network, which mainly included N metabolism, amino acid metabolism and synthesis, response to NO₃ ^- or NH₄ ⁺ and remobilization of N. Together, these results provide new insights into the potential role of circRNAs in poplar root responses to different N forms.

Genome-wide identification of lncRNAs, miRNAs, mRNAs and their regulatory networks involved in tapping panel dryness in rubber tree (Hevea brasiliensis)

Noncoding RNAs (ncRNAs) play pivotal roles in various biological processes in plants. However, the role of ncRNAs in tapping panel dryness (TPD) of rubber tree (Hevea brasiliensis Muell. Arg.) is largely unknown. Here, the whole transcriptome analyses of bark tissues from healthy and TPD trees were performed to identify differentially expressed long ncRNAs (DELs), microRNAs/miRNAs (DEMs), genes (DEGs) and their regulatory networks involved in TPD. A total of 263 DELs, 174 DEMs and 1574 DEGs were identified in the bark of TPD tree compared with that of healthy tree. Kyoto Encyclopedia of Genes and Genomes analysis revealed that most of the DEGs and targets of DELs and DEMs were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites and plant hormone signal transduction. Additionally, the majority of DEGs and DELs related to rubber biosynthesis were downregulated in TPD trees. Furthermore, 98 DEGs and 44 DELs were targeted by 54 DEMs, 190 DEGs were identified as putative targets of 56 DELs, and 2 and 44 DELs were predicted as precursors and endogenous target mimics of 2 and 6 DEMs, respectively. Based on these, the DEL-DEM-DEG regulatory network involved in TPD was constructed, and 13 hub DELs, 3 hub DEMs and 2 hub DEGs were identified. The results provide novel insights into the regulatory roles of ncRNAs underlying TPD and lay a foundation for future functional characterization of long ncRNAs, miRNAs and genes involved in TPD in rubber tree.

Nitrate/ammonium-responsive microRNA-mRNA regulatory networks affect root system architecture in Populus × canescens

BACKGROUND

Nitrate (NO₃^-) and ammonium (NH₄⁺) are the primary forms of inorganic nitrogen (N) taken up by plant roots, and a lack of these N sources commonly limits plant growth. To better understand how NO₃^- and NH₄⁺ differentially affect root system architecture, we analyzed the expression profiles of microRNAs and their targets in poplar roots treated with three forms of nitrogen S1 (NO₃^-), S2 (NH₄NO₃, normal), and S3 (NH₄⁺) via RNA sequencing.

RESULTS

The results revealed a total of 709 miRNAs. Among them, 57 significantly differentially expressed miRNAs and 28 differentially expressed miRNA-target pairs showed correlated expression profiles in S1 vs. S2. Thirty-six significantly differentially expressed miRNAs and 12 differentially expressed miRNA-target pairs showed correlated expression profiles in S3 vs. S2. In particular, NFYA3, a target of upregulated ptc-miR169i and ptc-miR169b, was downregulated in S1 vs. S2, while NFYA1, a target of upregulated ptc-miR169b, was downregulated in S3 vs. S2 and probably played an important role in the changes in root morphology observed when the poplar plants were treated with different N forms. Furthermore, the miRNA-target pairs ptc-miR169i/b-D6PKL2, ptc-miR393a-5p-AFB2, ptc-miR6445a-NAC14, ptc-miR172d-AP2, csi-miR396a-5p_R + 1_1ss21GA-EBP1, ath-miR396b-5p_R + 1-TPR4, and ptc-miR166a/b/c-ATHB-8 probably contributed to the changes in root morphology observed when poplar plants were treated with different N forms.

CONCLUSIONS

These results demonstrate that differentially expressed miRNAs and their targets play an important role in the regulation of the poplar root system architecture by different N forms.

Identification and expression analysis of the PtGATL genes under different nitrogen and carbon dioxide treatments in Populus trichocarpa

Pectin is one of the most important components of the plant cell wall. Galacturonosyltransferase-like (GATL) is an important enzyme involved in forming pectin in Arabidopsis thaliana. In this study, 12 PtGATL genes were identified and characterized based on the Populus trichocarpa genome using bioinformatics methods. The results showed that the PtGATLs contained four typical motifs, including DXD, LPPF, GLG, and HXXGXXKPW. According to phylogenetic analysis, PtGATLs were divided into six groups. Chromosome distribution and genome synteny analysis showed that there were 11 segmental-duplicated gene pairs with repeated fragments on chromosomes 2, 5, 7, 8, 10, and 14. Tissue-specific expression profiles indicated that these PtGATLs had different expression patterns. The transcription level of PtGATLs was regulated by different carbon dioxide and nitrogen concentrations. In conclusion, the identification and analysis of PtGATL genes in poplar provide important information on the gene function.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-022-03129-y.

Physiological Characteristics and Transcriptomic Dissection in Two Root Segments with Contrasting Net Fluxes of Ammonium and Nitrate of Poplar Under Low Nitrogen Availability

To investigate physiological and transcriptomic regulation mechanisms underlying the distinct net fluxes of NH4+ and NO3- in different root segments of Populus species under low nitrogen (N) conditions, we used saplings of Populus × canescens supplied with either 500 (normal N) or 50 (low N) μM NH4NO3. The net fluxes of NH4+ and NO3-, the concentrations of NH4+, amino acids and organic acids and the enzymatic activities of nitrite reductase (NiR) and glutamine synthetase (GS) in root segment II (SII, 35-70 mm to the apex) were lower than those in root segment I (SI, 0-35 mm to the apex). The net NH4+ influxes and the concentrations of organic acids were elevated, whereas the concentrations of NH4+ and NO3- and the activities of NiR and GS were reduced in SI and SII in response to low N. A number of genes were significantly differentially expressed in SII vs SI and in both segments grown under low vs normal N conditions, and these genes were mainly involved in the transport of NH4+ and NO3-, N metabolism and adenosine triphosphate synthesis. Moreover, the hub gene coexpression networks were dissected and correlated with N physiological processes in SI and SII under normal and low N conditions. These results suggest that the hub gene coexpression networks play pivotal roles in regulating N uptake and assimilation, amino acid metabolism and the levels of organic acids from the tricarboxylic acid cycle in the two root segments of poplars in acclimation to low N availability.

Genome-Wide Identification and Characterization of Long Noncoding RNAs in Populus × canescens Roots Treated With Different Nitrogen Fertilizers

Nitrate (NO₃ ^-) and ammonium (NH₄ ⁺) are the primary forms of inorganic nitrogen acquired by plant roots. LncRNAs, as key regulators of gene expression, are a class of non-coding RNAs larger than 200 bp. However, knowledge about the regulatory role of lncRNAs in response to different nitrogen forms remains limited, particularly in woody plants. Here, we performed strand-specific RNA-sequencing of P. × canescens roots under three different nitrogen fertilization treatments. In total, 324 lncRNAs and 6,112 mRNAs were identified as showing significantly differential expression between the NO₃ ^- and NH₄NO₃ treatments. Moreover, 333 lncRNAs and 6,007 mRNAs showed significantly differential expression between the NH₄ ⁺ and NH₄NO₃ treatments. Further analysis suggested that these lncRNAs and mRNAs have different response mechanisms for different nitrogen forms. In addition, functional annotation of cis and trans target mRNAs of differentially expressed lncRNAs indicated that 60 lncRNAs corresponding to 49 differentially expressed cis and trans target mRNAs were involved in plant nitrogen metabolism and amino acid biosynthesis and metabolism. Furthermore, 42 lncRNAs were identified as putative precursors of 63 miRNAs, and 28 differentially expressed lncRNAs were potential endogenous target mimics targeted by 96 miRNAs. Moreover, ceRNA regulation networks were constructed. MSTRG.6097.1, MSTRG.13550.1, MSTRG.2693.1, and MSTRG.12899.1, as hub lncRNAs in the ceRNA networks, are potential candidate lncRNAs for studying the regulatory mechanism in poplar roots under different nitrogen fertilization treatments. The results provide a basis for obtaining insight into the molecular mechanisms of lncRNA responses to different nitrogen forms in woody plants.

Transcriptomic Analysis Reveals Regulatory Networks for Osmotic Water Stress and Rewatering Response in the Leaves of Ginkgo biloba

To elucidate the transcriptomic regulation mechanisms that underlie the response of Ginkgo biloba to dehydration and rehydration, we used ginkgo saplings exposed to osmotically driven water stress and subsequent rewatering. When compared with a control group, 137, 1453, 1148, and 679 genes were differentially expressed in ginkgo leaves responding to 2, 6, 12, and 24 h of water deficit, and 796 and 1530 genes were differentially expressed responding to 24 and 48 h of rewatering. Upregulated genes participated in the biosynthesis of abscisic acid, eliminating reactive oxygen species (ROS), and biosynthesis of flavonoids and bilobalide, and downregulated genes were involved in water transport and cell wall enlargement in water stress-treated ginkgo leaves. Under rehydration conditions, the genes associated with water transport and cell wall enlargement were upregulated, and the genes that participated in eliminating ROS and the biosynthesis of flavonoids and bilobalide were downregulated in the leaves of G. biloba. Furthermore, the weighted gene coexpression networks were established and correlated with distinct water stress and rewatering time-point samples. Hub genes that act as key players in the networks were identified. Overall, these results indicate that the gene coexpression networks play essential roles in the transcriptional reconfiguration of ginkgo leaves in response to water stress and rewatering.

Morphological, physiological, and transcriptional responses to low nitrogen stress in Populus deltoides Marsh. clones with contrasting nitrogen use efficiency

Background

Nitrogen (N) is one of the main factors limiting the wood yield in poplar cultivation. Understanding the molecular mechanism of N utilization could play a guiding role in improving the nitrogen use efficiency (NUE) of poplar.

Results

In this study, three N-efficient genotypes (A1-A3) and three N-inefficient genotypes (C1-C3) of Populus deltoides were cultured under low N stress (5 μM NH₄NO₃) and normal N supply (750 μM NH₄NO₃). The dry matter mass, leaf morphology, and chlorophyll content of both genotypes decreased under N starvation. The low nitrogen adaptation coefficients of the leaves and stems biomass of group A were significantly higher than those of group C (p < 0.05). Interestingly, N starvation induced fine root growth in group A, but not in group C. Next, a detailed time-course analysis of enzyme activities and gene expression in leaves identified 2062 specifically differentially expressed genes (DEGs) in group A and 1118 in group C. Moreover, the sensitivity to N starvation of group A was weak, and DEGs related to hormone signal transduction and stimulus response played an important role in the low N response this group. Weighted gene co-expression network analysis identified genes related to membranes, catalytic activity, enzymatic activity, and response to stresses that might be critical for poplar’s adaption to N starvation and these genes participated in the negative regulation of various biological processes. Finally, ten influential hub genes and twelve transcription factors were identified in the response to N starvation. Among them, four hub genes were related to programmed cell death and the defense response, and PodelWRKY18, with high connectivity, was involved in plant signal transduction. The expression of hub genes increased gradually with the extension of low N stress time, and the expression changes in group A were more obvious than those in group C.

Conclusions

Under N starvation, group A showed stronger adaptability and better NUE than group C in terms of morphology and physiology. The discovery of hub genes and transcription factors might provide new information for the analysis of the molecular mechanism of NUE and its improvement in poplar.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07991-7.

Bioinformatics analysis of PAE family in Populus trichocarpa and responsiveness to carbon and nitrogen treatment

Plant Pectin acetylesterase (PAE) belongs to family CE13 of carbohydrate esterases in the CAZy database. The ability of PAE to regulate the degree of acetylation of pectin, an important polysaccharide in the cell wall, affects the structure of plant cell wall. In this study, ten PtPAE genes were identified and characterized in Populus trichocarpa genome using bioinformatics methods, and the physiochemical properties such as molecular weight, isoelectric points, and hydrophilicity, as well as the secondary and tertiary structure of the protein were predicted. According to phylogenetic analysis, ten PtPAEs can be divided into three evolutionary clades, each of which had similar gene structure and motifs. Tissue-specific expression profiles indicated that the PtPAEs had different expression patterns. Real-time quantitative PCR (RT-qPCR) analysis showed that transcription level of PtPAEs was regulated by different CO₂ and nitrogen concentrations. These results provide important information for the study of the phylogenetic relationship and function of PtPAEs in Populus trichocarpa.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02918-1.

Wood properties and transcriptional responses of poplar hybrids in mixed cropping with the nitrogen-fixing species Robinia pseudoacacia

Cultivation of fast-growing tree species is often confined to marginal land. Mixed cropping with Robinia pseudoacacia, a legume tree species that forms a symbiosis with N2-fixing bacteria, has been proposed to be a measure to improve soil fertility and to achieve beneficial effects on the cocultivated tree species. The goal of our study was to examine the influence of a Robinia mixture on wood chemistry, anatomy and gene expression in poplar. We hypothesized that annual wood growth is stimulated in species mixtures due to the positive effects of Robinia on nitrogen availability and complementary resource use. Alternatively, we hypothesized that competition, especially for water, has negative effects on the wood growth of poplar. We used two commercial biomass clones, Hybride 275 (H275, Populus trichocarpa × Populus maximowiczii) and Max1 (Populus nigra × P. maximowiczii), which were planted at two locations with contrasting soil fertility in monoculture or mixed plots with Robinia to investigate the annual wood increment, wood nitrogen and δ13C, wood anatomy (length, cell wall thickness, lumina and frequencies of fibers and vessels) and transcriptional profiles in the developing xylem of 4-year-old stems. In a mixture with Robinia, the annual stem increment was reduced, nitrogen in wood was enhanced, δ13C in wood was decreased, vessel and fiber frequencies were increased and fiber lengths and fiber lumina were decreased. Transcriptional profiles showed stronger differences between the genotypes and sites than between mono and mixed cultivation. The transcriptional abundances of only one gene (the putative nitrate transporter, NRT1.2) and one gene ontology term ('immune system process') were significantly enriched in wood-forming tissues in response to the mixture, irrespective of the poplar genotype and growth location. Weighted gene coexpression network analyses extracted gene modules that linked wood nitrogen mainly to vessel traits and wood δ13C with fiber traits. Collectively, molecular and anatomical changes in poplar wood suggest beneficial effects on the water and N supply in response to the mixture with Robinia. These alterations may render poplars less drought-susceptible. However, these benefits are accompanied by a reduced wood increment, emphasizing that other critical factors, presumably light competition or allelopathic effects, overrule a potential growth stimulation.

Synergies and Entanglement in Secondary Cell Wall Development and Abiotic Stress Response in Trees

A major challenge for sustainable food, fuel, and fiber production is simultaneous genetic improvement of yield, biomass quality, and resilience to episodic environmental stress and climate change. For Populus and other forest trees, quality traits involve alterations in the secondary cell wall (SCW) of wood for traditional uses, as well as for a growing diversity of biofuels and bioproducts. Alterations in wood properties that are desirable for specific end uses can have negative effects on growth and stress tolerance. Understanding of the diverse roles of SCW genes is necessary for the genetic improvement of fast-growing, short-rotation trees that face perennial challenges in their growth and development. Here, we review recent progress into the synergies and antagonisms of SCW development and abiotic stress responses, particularly, the roles of transcription factors, SCW biogenesis genes, and paralog evolution.

Genome-wide analysis of UGDH genes in Populus trichocarpa and responsiveness to nitrogen treatment

Plant UDP-glucose 6-dehydrogenase (UGDH) is an important enzyme for the formation of hemicellulose and pectin. Previous studies on UGDH have primarily focused on the biosynthesis of cell wall polysaccharides, while few studies have focused on their regulation by exogenous nitrogen. In the present study, four genes encoding PtUGDH proteins were analyzed by bioinformatics methods. And, the expression profiles of PtUGDH genes under different nitrogen treatments were evaluated with qRT-PCR. The results showed that PtUGDHs have conserved NAD coenzyme binding motif GAGYVGG and the catalytic motif GFGGSCFQKDIL. According to the phylogenetic analysis, PtUGDHs were divided into two subgroups. PtUGDH3 and PtUGDH4 were closely related to AtUGDH1 (important for normal development of Arabidopsis cell wall structure). Chromosomal distribution and genome synteny analysis revealed four segmental-duplicated gene pairs on chr4, 8, 10 and 17. Tissue-specific data from PlantGenIE showed that PtUGDH3 and PtUGDH4 were highly expressed in stems. The qRT-PCR detection showed that the expression of PtUGDH3 in the lower stem and PtUGDH2 of upper leaves were significantly increased induced by low ammonium or nitrate condition. This comprehensive analysis of the UGDH family in poplar provides new insights into their regulation by nitrogen, and would increase our understanding of the roles of UGDHs in hemicellulose and pectin biosynthesis in the cell wall and during poplar development.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02697-9.

Whole-Transcriptome RNA Sequencing Reveals the Global Molecular Responses and CeRNA Regulatory Network of mRNAs, lncRNAs, miRNAs and circRNAs in Response to Salt Stress in Sugar Beet (Beta vulgaris)

Sugar beet is an important sugar-yielding crop with some tolerance to salt, but the mechanistic basis of this tolerance is not known. In the present study, we have used whole-transcriptome RNA-seq and degradome sequencing in response to salt stress to uncover differentially expressed (DE) mRNAs, microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in both leaves and roots. A competitive endogenous RNA (ceRNA) network was constructed with the predicted DE pairs, which revealed regulatory roles under salt stress. A functional analysis suggests that ceRNAs are implicated in copper redistribution, plasma membrane permeability, glycometabolism and energy metabolism, NAC transcription factor and the phosphoinositol signaling system. Overall, we conducted for the first time a full transcriptomic analysis of sugar beet under salt stress that involves a potential ceRNA network, thus providing a basis to study the potential functions of lncRNAs/circRNAs.

Dissecting MicroRNA-mRNA Regulatory Networks Underlying Sulfur Assimilation and Cadmium Accumulation in Poplar Leaves

The process of cadmium (Cd) accumulation and detoxification under different sulfur levels remains largely unknown in woody plants. To investigate the physiological and transcriptomic regulation mechanisms of poplars in response to different sulfate (S) supply levels and Cd exposure, we exposed Populus deltoides saplings to one of the low, moderate and high S levels together with either 0 or 50 µM Cd. Cd accumulation was decreased in low S-treated poplar leaves, and it tended to be increased in high S-supplied leaves under the Cd exposure condition. Sulfur nutrition was deficient in low S-supplied poplars, and it was improved in high S-treated leaves. Cd exposure resulted in lower sulfur level in the leaves supplied with moderate S, it exacerbated a Cd-induced sulfur decrease in low S-treated leaves and it caused a higher sulfur concentration in high S-supplied leaves. In line with the physiological changes, a number of mRNAs and microRNAs (miRNAs) involved in Cd accumulation and sulfur assimilation were identified and the miRNA-mRNA networks were dissected. In the networks, miR395 and miR399 members were identified as hub miRNAs and their targets were ATP sulfurylase 3 (ATPS3) and phosphate 2 (PHO2), respectively. These results suggest that Cd accumulation and sulfur assimilation are constrained by low and enhanced by high S supply, and Cd toxicity is aggravated by low and relieved by high S in poplar leaves, and that miRNA-mRNA regulatory networks play pivotal roles in sulfur-mediated Cd accumulation and detoxification in Cd-exposed poplars.

MicroRNA comparison between poplar and larch provides insight into the different mechanism of wood formation

MiRNA transcriptome analysis of different tissues in poplar and larch suggests variant roles of miRNAs in regulating wood formation between two kinds of phyla. Poplar and larch belong to two different phyla. Both are ecological woody species and major resources for wood-related industrial applications. However, wood properties are different between these two species and the molecular basis is largely unknown. In this study, we performed high-throughput sequencing of microRNAs (miRNAs) in the three tissues, xylem, phloem and leaf of Populus alba × Populus glandulosa and Larix kaempferi. Differentially expressed miRNA (DEmiRNA) analysis identified 85 xylem-specific miRNAs in P. alba × P. glandulosa and 158 xylem-specific miRNAs in L. kaempferi. Among 36 common miRNAs, 12 were conserved between the two species. GO and KEGG analyses of the miRNA target genes showed similar metabolism in two species. Through KEGG and BLASTN, we predicted target genes of xylem differentially expressed (DEmiRNA) in the wood formation-related pathways and located DEmiRNAs in these pathways. A network was built for wood formation-related DEmiRNAs, their target genes and orthologous genes in Arabidopsis thaliana. Comparison of DEmiRNA and target gene annotation between P. alba × P. glandulosa and L. kaempferi suggested the different functions of DEmiRNAs and divergent mechanism in wood formation between two species, providing knowledge to understand wood formation mechanism in gymnosperm and angiosperm woody plants.

Genome-wide identification of BXL genes in Populus trichocarpa and their expression under different nitrogen treatments

β-d-xylosidase (BXL) hydrolyzes xylobiose and xylo-oligosaccharides into xylose monomers, and is a rate-limiting enzyme in the degradation of hemicellulose in the cell wall. In this study, ten genes encoding putative BXL proteins were identified in the Populus trichocarpa genome by bioinformatics methods. In the phylogenetic analysis, the PtBXLs formed two subfamilies. PtBXL8 and PtBXL9 were closely related to AtBXL1, an important enzyme in the normal development of the Arabidopsis cell wall structure. Chromosomal distribution and genome synteny analyses revealed two tandem-duplicated gene pairs PtBXL3/4 and PtBXL6/7 on chromosomes II and V, respectively, and six segmental-duplicated gene pairs on chromosomes II, V, VIII, X, and XIV among the PtBXL gene family. Tissue-specific expression data from PlantGenIE indicated that PtBXL2, 4, 5, and 10 were highly expressed in stems. Quantitative real-time RT-PCR analyses revealed that PtBXL4, 5, and 9 were up-regulated in the upper stem in response to the low and high ammonium and nitrate treatments. The influence of nitrogen on the expression of PtBXL4, 5, and 9 genes may affect the formation of the plant secondary cell wall. This comprehensive analysis of the BXL family in poplar provides new insights into their regulation by nitrogen and increases our understanding of the roles of BXLs in hemicellulose metabolism in the secondary cell wall and during plant development.

Identification and characterization of circular RNAs during wood formation of poplars in acclimation to low nitrogen availability

Circular RNA (circRNA) identification and expression profiles, and construction of circRNAs-miRNAs-mRNAs networks indicates that circRNAs are involved in wood formation of poplars in acclimation to low nitrogen availability. Circular RNAs (circRNAs) are covalently closed non-coding RNAs that play pivotal roles in various biological processes. However, circRNAs' roles in wood formation of poplars in acclimation to low nitrogen (N) availability are currently unknown. Here, we undertook a systematic identification and characterization of circRNAs in the wood of Populus × canescens exposed to either 50 (low N) or 500 (normal N) µM NH₄NO₃ using rRNA-depleted RNA-sequencing. A total of 2,509 unique circRNAs were identified, and 163 (ca. 6.5%) circRNAs were significantly differentially expressed (DE) under low N condition. We observed a positive correlation between the expression patterns of DE circRNAs and their hosting protein-coding genes. Moreover, circRNAs-miRNAs-mRNAs' networks were identified in the wood of poplars under low N availability. For instance, upregulated several circRNAs, such as circRNA1226, circRNA 1732, and circRNA392 induced increases in nuclear factor Y, subunit A1-A (NFYA1-A), NFYA1-B, and NFYA10 transcript levels via the mediation of miR169b members, which is in line with reduced xylem width and cell layers of the xylem in the wood of low N-supplied poplars. Upregulation of circRNA1006, circRNA1344, circRNA1941, circRNA901, and circRNA146 caused increased transcript level of MYB61 via the mediation of a miR5021 member, corresponding well to the higher lignin concentration in the wood of low N-treated poplars. Overall, these results indicated that DE circRNAs play an essential role in regulating gene expression via circRNAs-miRNAs-mRNAs' networks to modulate wood anatomical and chemical properties of poplars in acclimation to low N availability.