Conjoint Analysis of Genome-Wide lncRNA and mRNA Expression of Heteromorphic Leavesin Response to Environmental Heterogeneityin Populus euphratica

Unraveling the lncRNA-miRNA-mRNA Regulatory Network Involved in Poplar Coma Development through High-Throughput Sequencing

Poplar coma, the fluff-like appendages of seeds originating from the differentiated surface cells of the placenta and funicle, aids in the long-distance dispersal of seeds in the spring. However, it also poses hazards to human safety and causes pollution in the surrounding environment. Unraveling the regulatory mechanisms governing the initiation and development of coma is essential for addressing this issue comprehensively. In this study, strand-specific RNA-seq was conducted at three distinct stages of coma development, revealing 1888 lncRNAs and 52,810 mRNAs. The expression profiles of lncRNAs and mRNAs during coma development were analyzed. Subsequently, potential target genes of lncRNAs were predicted through co-localization and co-expression analyses. Integrating various types of sequencing data, lncRNA-miRNA-TF regulatory networks related to the initiation of coma were constructed. Utilizing identified differentially expressed genes encoding kinesin and actin, lncRNA-miRNA-mRNA regulatory networks associated with the construction and arrangement of the coma cytoskeleton were established. Additionally, relying on differentially expressed genes encoding cellulose synthase, sucrose synthase, and expansin, lncRNA-miRNA-mRNA regulatory networks related to coma cell wall synthesis and remodeling were developed. This study not only enhances the comprehension of lncRNA but also provides novel insights into the molecular mechanisms governing the initiation and development of poplar coma.

Morphological, Anatomical, and Physiological Characteristics of Heteroblastic Acacia melanoxylon Grown under Weak Light

Acacia melanoxylon is a fast-growing macrophanerophyte with strong adaptability whose leaf enables heteromorphic development. Light is one of the essential environmental factors that induces the development of the heteroblastic leaf of A. melanoxylon, but its mechanism is unclear. In this study, the seedlings of A. melanoxylon clones were treated with weak light (shading net with 40% of regular light transmittance) and normal light (control) conditions for 90 d and a follow-up observation. The results show that the seedlings' growth and biomass accumulation were inhibited under weak light. After 60 days of treatment, phyllodes were raised under the control condition while the remaining compound was raised under weak light. The balance of root, stem, and leaf biomass changed to 15:11:74 under weak light, while it was 40:15:45 under control conditions. After comparing the anatomical structures of the compound leaves and phyllode, they were shown to have their own strategies for staying hydrated, while phyllodes were more able to control water loss and adapt to intense light. The compound leaves exhibited elevated levels of K, Cu, Ca, and Mg, increased antioxidant enzyme activity and proline content, and higher concentrations of chlorophyll a, carotenoids, ABA, CTK, and GA. However, they displayed a relatively limited photosynthetic capacity. Phyllodes exhibited higher levels of Fe, cellulose, lignin, IAA content, and high photosynthetic capacity with a higher maximum net photosynthetic rate, light compensation point, dark respiration rate, and water use efficiency. The comparative analysis of compound leaves and phyllodes provides a basis for understanding the diverse survival strategies that heteroblastic plants employ to adapt to environmental changes.

Identification of Differentially Expressed lncRNAs in Response to Blue Light and Expression Pattern Analysis of Populus tomentosa Hybrid Poplar 741

Poplar is an important shelterbelt, timber stand, and city tree species that has been the focus of forestry research. The regulatory role of the long non-coding RNA molecule (lncRNA; length > 200 nt) has been a research hotspot in plants. In this study, seedlings of 741 poplar were irradiated with LED blue and white light, and the Illumina HiSeq 2000 sequencing platform was used to identify lncRNAs. |logFC| > 1 and p < 0.05 were considered to indicate differentially expressed lncRNAs, and nine differentially expressed lncRNAs were screened, the target genes of which were predicted, and three functionally annotated target genes were obtained. The differentially expressed lncRNAs were identified as miRNA targets. Six lncRNAs were determined to be target sites for twelve mRNAs in six miRNA families. LncRNAs and their target genes, including lncRNA MSTRG.20413.1-ptc-miR396e-5p-GRF9, were verified using quantitative real-time polymerase chain reaction analysis, and the expression patterns were analyzed. The analysis showed that the ptc-miR396e-5p expression was downregulated, while lncRNA MSTRG.20413.1 and GRF9 expression was upregulated, after blue light exposure. These results indicate that lncRNAs interact with miRNAs to regulate gene expression and affect plant growth and development.

Comprehensive analysis of lncRNA-mRNA regulatory network in Populus associated with xylem development

Long noncoding RNAs (lncRNAs) play essential roles in numerous biological processes in plants, such as regulating the gene expression. However, only a few studies have looked into their potential functions in xylem development. High-throughput sequencing of P. euramericana 'Zhonglin46' developing and mature xylem was performed in this study. Through sequencing analysis, 14,028 putative lncRNA transcripts were identified, including 4525 differentially expressed lncRNAs (DELs). Additional research revealed that in mature xylem, a total of 2320 DELs were upregulated and 2205 were downregulated compared to developing xylem. Meanwhile, there were a total of 8122 differentially expressed mRNAs (DEMs) that were upregulated and 16,424 that were downregulated in mature xylem compared with developing xylem. The cis- and trans-target genes of DELs were analyzed for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, which indicated that these DELs participate in controlling the phenylpropanoid and lignin biosynthesis pathway as well as the starch and sucrose metabolism pathway. Among the cis-regulated DELs, LNC_006291, LNC_006292, and LNC_006532 all participate in regulating multiple HCT gene family membranes. As targets, POPTR_001G045900v3 (CCR2) and POPTR_018G063500v3 (SUS) both have only one cis-regulatory lncRNA, referred to as LNC_000057 and LNC_006212, respectively. Moreover, LNC_004484 and two DELs named LNC_008014 and LNC_010781 were revealed to be important nodes in the co-expression network of trans-lncRNAs and mRNAs associated to the lignin biosynthesis pathway and cellulose and xylan biosynthetic pathways, respectively. Finally, quantitative real-time PCR (qRT-PCR) was used to confirme 34 pairs of lncRNA-mRNA. Taken together, these findings may help to clarify the regulatory role that lncRNAs play in xylem development and wood formation.

Drought stress-induced the formation of heteromorphic leaves of Populus euphratica Oliv: evidence from gene transcriptome

Populus euphratica Oliv., a dominant species of arid desert community, grows heteromorphic leaves at different crown positions. Whether heteromorphic leaves are a strategy of plant adaptation to drought stress is rarely reported. This study sequenced the transcriptome of three typical heteromorphic leaves (lanceolate, ovate and broad-ovate leaves) of P. euphratica, and measured their drought stress. We wanted to reveal the molecular mechanisms underlying the formation of heteromorphic leaves. Drought stress was increased significantly from lanceolate to ovate to broad-ovate leaves. Gene ontology (GO) and KEGG enrichment analysis showed that the MADs-box gene regulated the expression of peroxidase (POD) in the phenylpropane biosynthetic pathway. The up-regulated expression of the chalcone synthase (CHS) gene in broad-ovate leaves significantly activated the flavonoid biosynthetic pathway. In the process of leaf shape change, the different expressions of homeodomain leucine zipper (HD-ZIP) among the three heteromorphic leaves had potential interactions on the AUX and ABA pathways. The expression of Sucrose phosphate synthase (SPS) and sucrose synthase (SUS) increased from lanceolate to broad-ovate leaves, resulting in a consistent change in starch and sucrose content. We concluded that these resistance-related pathways are expressed in parallel with leaf formation genes, thereby inducing the formation of heteromorphic leaves. Our work provided a new insights for desert plants to adapt to drought stress.

Multi-omics analysis reveals spatiotemporal regulation and function of heteromorphic leaves in Populus

Despite the high economic and ecological importance of forests, our knowledge of the adaptive evolution of leaf traits remains very limited. Euphrates poplar (Populus euphratica), which has high tolerance to arid environment, has evolved four heteromorphic leaf forms, including narrow (linear and lanceolate) and broad (ovate and broad-ovate) leaves on different crowns. Here, we revealed the significant functional divergence of four P. euphratica heteromorphic leaves at physiological and cytological levels. Through global analysis of transcriptome and DNA methylation across tree and leaf developmental stages, we revealed that gene expression and DNA epigenetics differentially regulated key processes involving development and functional adaptation of heteromorphic leaves, such as hormone signaling pathways, cell division, and photosynthesis. Combined analysis of gene expression, methylation, ATAC-seq, and Hi-C-seq revealed longer interaction of 3D genome, hypomethylation, and open chromatin state upregulates IAA-related genes (such as PIN-FORMED1 and ANGUSTIFOLIA3) and promotes the occurrence of broad leaves while narrow leaves were associated with highly concentrated heterochromatin, hypermethylation, and upregulated abscisic acid pathway genes (such as Pyrabactin Resistance1-like10). Therefore, development of P. euphratica heteromorphic leaves along with functional divergence was regulated by differentially expressed genes, DNA methylation, chromatin accessibility, and 3D genome remodeling to adapt to the arid desert. This study advances our understanding of differential regulation on development and functional divergence of heteromorphic leaves in P. euphratica at the multi-omics level and provides a valuable resource for investigating the adaptive evolution of heteromorphic leaves in Populus.

Conjoint Analysis of Genome-Wide lncRNA and mRNA Expression during the Salicylic Acid Response in Populus × euramericana

Long noncoding RNAs (lncRNAs) participate in a wide range of biological processes, but lncRNAs in plants remain largely unknown; in particular, we lack a systematic identification of plant lncRNAs involved in hormone responses. To explore the molecular mechanism of the response of poplar to salicylic acid (SA), the changes in protective enzymes, which are closely related to plant resistance induced by exogenous SA, were studied, and the expression of mRNA and lncRNA were determined by high-throughput RNA sequencing. The results showed that the activities of phenylalanine ammonia lyase (PAL) and polyphenol oxidase (PPO), in the leaves of Populus × euramericana, were significantly increased by exogenous SA application. High-throughput RNA sequencing showed that 26,366 genes and 5690 lncRNAs were detected under the different treatment conditions: SA and H₂O application. Among these, 606 genes and 49 lncRNAs were differentially expressed. According to target prediction, lncRNAs and target genes involved in light response, stress response, plant disease resistance, and growth and development, were differentially expressed in SA-treated leaves. Interaction analysis showed that lncRNA-mRNA interactions, following exogenous SA, were involved in the response of poplar leaves to the external environment. Our study provides a comprehensive view of Populus × euramericana lncRNAs and offers insights into the potential functions and regulatory interactions of SA-responsive lncRNAs, thus forming the foundation for future functional analysis of SA-responsive lncRNAs in Populus × euramericana.

Genome-wide analysis of long non-coding RNAs in shoot apical meristem and vascular cambium in Populus tomentosa

Shoot apical and lateral meristems play essential roles in the formation and development of primary and secondary growth in plants. A delicate regulatory mechanism is needed to maintain homeostatic balance between the primary and secondary growth, as well as the self-renewal of meristems with the rate of cell division and differentiation of new meristems. However, little is known about the roles of long non-coding RNAs (lncRNAs) in the regulation of maintenance and differentiation of primary and secondary growth in Populus, especially in the cambium division and differentiation into secondary xylem. Here, 1298 lncRNAs were identified both in the apical meristem and vascular cambium, with 80 lncRNAs being expressed only in shoot apical meristem and 45 only in vascular cambium. There are 410 differentially expressed lncRNAs in shoot apical meristem and vascular cambium, among which 271 lncRNAs were up-regulated and 139 were down-regulated in cambium. The GO enrichment analysis revealed that differentially expressed lncRNAs mainly influenced the expression of lncRNAs related to the ribosome pathway, plant hormone signal pathway and photosynthesis pathway. The differentially expressed lncRNAs mainly target mRNA through cis-regulation in the vascular cambium. In addition, six key lncRNAs and also their significantly upregulated target genes were identified. Theses target genes are involved in plant secondary metabolites, cellulose and lignin synthesis, hormone and signal transduction. In addition, six key lncRNAs were identified, their significantly upregulated target genes are related to plant secondary metabolites, cellulose and lignin synthesis, hormone and signal transduction. Investigating lncRNA-mRNA interactions, we further found some genes that may be related to the development of vascular cambium, such as domain-containing transcription factors, cellulose synthesis genes, calcium dependent protein kinase 2, cytokinin receptor 1, glycosyl transferase and polyphenol oxidase. Our findings provide new insights into the lncRNA-mRNA networks in the development of vascular cambium of secondary growth in Populus.

Genetic Architecture of Heterophylly: Single and Multi-Leaf Genome-Wide Association Mapping in Populus euphratica

Heterophylly is an adaptive strategy used by some plants in response to environmental changes. Due to the lack of representative plants with typical heteromorphic leaves, little is known about the genetic architecture of heterophylly in plants and the genes underlying its control. Here, we investigated the genetic characteristics underlying changes in leaf shape based on the model species, Populus euphratica, which exhibits typical heterophylly. A set of 401,571 single-nucleotide polymorphisms (SNPs) derived from whole-genome sequencing of 860 genotypes were associated with nine leaf traits, which were related to descriptive and shape data using single- and multi-leaf genome-wide association studies (GWAS). Multi-leaf GWAS allows for a more comprehensive understanding of the genetic architecture of heterophylly by considering multiple leaves simultaneously. The single-leaf GWAS detected 140 significant SNPs, whereas the multi-leaf GWAS detected 200 SNP-trait associations. Markers were found across 19 chromosomes, and 21 unique genes were implicated in traits and serve as potential targets for selection. Our results provide novel insights into the genomic architecture of heterophylly, and provide candidate genes for breeding or engineering P. euphratica. Our observations also improve understanding of the intrinsic mechanisms of plant growth, evolution, and adaptation in response to climate change.

iTRAQ-based proteomic analysis of heteromorphic leaves reveals eco-adaptability of Populus euphratica Oliv

BACKGROUND

Heterophylly is regard as adaptation to different environments in plant, and Populus euphratica is an important heterophyllous woody plant. However, information on its molecular mechanism in eco-adaptability remains obscure.

RESULTS

In this research, proteins were identified by isobaric tags for relative and absolute quantitation (iTRAQ) technology in lanceolate, ovate, and dentate broad-ovate leaves from adult P. euphratica trees, respectively. Besides, chlorophyll content, net photosynthetic rate, stomatal conductance, transpiration rate and peroxidase activity in these heteromorphic leaves were investigated. A total number of 2,689 proteins were detected in the heteromorphic leaves, of which 56, 73, and 222 differential abundance proteins (DAPs) were determined in ovate/lanceolate, dentate broad-ovate/lanceolate, and dentate broad-ovate/ovate comparison groups. Bioinformatics analysis suggested these altered proteins related to photosynthesis, stress tolerance, respiration and primary metabolism accumulated in dentate broad-ovate and ovate leaves, which were consistent with the results of physiological parameters and Real-time Quantitative PCR experiments.

CONCLUSION

This research demonstrated the mechanism of the differential abundance proteins in providing an optimal strategy of resource utilization and survival for P. euphratica, that could offer clues for further investigations into eco-adaptability of heterophyllous woody plants.

Non-Coding RNA Analyses of Seasonal Cambium Activity in Populus tomentosa

Non-coding RNA, known as long non-coding RNA (lncRNA), circular RNA (circRNA) and microRNA (miRNA), are taking part in the multiple developmental processes in plants. However, the roles of which played during the cambium activity periodicity of woody plants remain poorly understood. Here, lncRNA/circRNA-miRNA-mRNA regulatory networks of the cambium activity periodicity in Populus tomentosa was constructed, combined with morphologic observation and transcriptome profiling. Light microscopy and Periodic Acid Schiff (PAS) staining revealed that cell walls were much thicker and number of cell layers was increased during the active-dormant stage, accompanied by abundant change of polysaccharides. The novel lncRNAs and circRNAs were investigated, and we found that 2037 lncRNAs and 299 circRNAs were differentially expression during the vascular cambium period, respectively. Moreover, 1046 genes were identified as a target gene of 2037 novel lncRNAs, and 89 of which were the miRNA precursors or targets. By aligning miRNA precursors to the 7655 lncRNAs, 21 lncRNAs were identified as precursors tof 19 known miRNAs. Furthermore, the target mRNA of lncRNA/circRNA-miRNA network mainly participated in phytohormone, cell wall alteration and chlorophyll metabolism were analyzed by GO enrichment and KEGG pathway. Especially, circRNA33 and circRNA190 taking part in the phytohormone signal pathway were down-regulated during the active-dormant transition. Xyloglucan endotransglucosylase/hydrolase protein 24-like and UDP-glycosyltransferase 85A1 involved in the cell wall modification were the targets of lncRNA MSTRG.11198.1 and MSTRG.1050.1. Notably, circRNA103 and MSTRG.10851.1 regulate the cambium periodicity may interact with the miR482. These results give a new light into activity–dormancy regulation, associated with transcriptional dynamics and non-coding RNA networks of potential targets identification.

Biologia Futura: progress and future perspectives of long non-coding RNAs in forest trees

Forest trees are affected by climate change, anthropogenic pressure, as well as abiotic and biotic stresses. Conventional tree breeding has so far been limited to enhance overall productivity, and our understanding of the genetic basis of quantitative traits is still inadequate. Quantum leaps in next-generation sequencing technologies and bioinformatics have permitted the exploration and identification of various non-coding regions of the genome other than protein coding genes. These genomic regions produce various types of non-coding RNAs and regulate myriads of biological functions at epigenetic, transcriptional and translational levels. Recently, long non-coding RNAs (lncRNAs) which act as molecular switch have been identified to be pivotal molecules in forest trees. This review focuses on progress made in regulatory mechanisms in various developmental phases like wood formation, adventitious rooting and flowering and stress responses. It was predicted that complex regulatory interactions among lncRNA, miRNA and gene exist. LncRNAs can function as a sponge for miRNAs, reducing the suppressive effect of miRNAs on target mRNAs and perhaps adding a new layer of regulatory interactions among non-coding RNA classes in trees. Furthermore, network analysis revealed the interactions of lncRNA and genes during the expression of several important genes. The insights generated about lncRNAs in forest trees would enable improvement of economically important traits including the devastating abiotic and biotic stresses. In addition, solid understanding on the wide range of regulatory functions of lncRNAs on traits influencing biomass productivity and adaptation would aid the applications of biotechnology in genetic improvement of forest trees.

Transcriptome Profile Analysis of Arabidopsis Reveals the Drought Stress-Induced Long Non-coding RNAs Associated With Photosynthesis, Chlorophyll Synthesis, Fatty Acid Synthesis and Degradation

Long non-coding RNAs (lncRNAs) play an important role in the response of plants to drought stress. The previous studies have reported that overexpression of LEA3 and VOC could enhance drought tolerance and improve the oil content in Brassica napus and Arabidopsis thaliana, and most of the efforts have been invested in the gene function analysis, there is little understanding of how genes that involved in these important pathways are regulated. In the present study, the transcriptomic results of LEA3 and VOC over-expressed (OE) lines were compared with the RNAi lines, mutant lines and control lines under long-term and short-term drought treatment, a series of differentially expressed lncRNAs were identified, and their regulation patterns in mRNA were also investigated in above mentioned materials. The regulation of the target genes of differentially expressed lncRNAs on plant biological functions was studied. It was revealed that the mutant lines had less drought-response related lncRNAs than that of the OE lines. Functional analysis demonstrated that multiple genes were involved in the carbon-fixing and chlorophyll metabolism, such as CDR1, CHLM, and CH1, were regulated by the upregulated lncRNA in OE lines. In LEA-OE, AT4G13180 that promotes the fatty acid synthesis was regulated by five lncRNAs that were upregulated under both long-term and short-term drought treatments. The key genes, including of SHM1, GOX2, and GS2, in the methylglyoxal synthesis pathway were all regulated by a number of down-regulated lncRNAs in OE lines, thereby reducing the content of such harmful compounds produced under stress in plants. This study identified a series of lncRNAs related to the pathways that affect photosynthesis, chlorophyll synthesis, fatty acid synthesis, degradation, and other important effects on drought resistance and oil content. The present study provided a series of lncRNAs for further improvement of crop varieties, especially drought resistant and oil content traits.

Progress in studying heteromorphic leaves in Populus euphratica: leaf morphology, anatomical structure, development regulation and their ecological adaptation to arid environments

Populus euphratica Oliv. is a tree that is strongly resistant to drought and salt stress, which is primarily distributed in arid and semiarid lands. The leaves of the species exhibit a special feature that causes them to be designated as heterophylly. In this brief review, we primarily discuss the heteromorphic leaf development and anatomical features, such as the differentiation of spongy and palisade tissues, in heteromorphic leaves of the species. Furthermore, we also discuss the different physiological characteristics in heteromorphic leaves related to the ecological adaptation of P. euphratica to drought environments. These traits in P. euphratica may help researchers evaluate its ecological value in arid areas and evaluate its scientific merit in understanding the mechanism of development of heteromorphic leaves in plants.