Genome-Wide Analysis of the miRNA-mRNAs Network Involved in Cold Tolerance in Populus simonii × P. nigra

Genome-Wide Identification of the ABC Gene Family and Its Expression in Response to the Wood Degradation of Poplar in Trametes gibbosa

Wood-rotting fungi's degradation of wood not only facilitates the eco-friendly treatment of organic materials, decreasing environmental pollution, but also supplies crucial components for producing biomass energy, thereby reducing dependence on fossil fuels. The ABC gene family, widely distributed in wood-rotting fungi, plays a crucial role in the metabolism of lignin, cellulose, and hemicellulose. Trametes gibbosa, as a representative species of wood-rotting fungi, exhibits robust capabilities in wood degradation. To investigate the function of the ABC gene family in wood degradation by T. gibbosa, we conducted a genome-wide analysis of T. gibbosa's ABC gene family. We identified a total of 12 Tg-ABCs classified into four subfamilies (ABCA, ABCB, ABCC, and ABCG). These subfamilies likely play significant roles in wood degradation. Scaffold localization and collinearity analysis results show that Tg-ABCs are dispersed on scaffolds and there is no duplication of gene sequences in the Tg-ABCs in the genome sequence of T. gibbosa. Phylogenetic and collinearity analyses of T. gibbosa along with four other wood-rotting fungi show that T. gibbosa shares a closer phylogenetic relationship with its same-genus fungus (Trametes versicolor), followed by Ganoderma leucocontextum, Laetiporus sulphureus, and Phlebia centrifuga in descending order of phylogenetic proximity. In addition, we conducted quantitative analyses of Tg-ABCs from T. gibbosa cultivated in both woody and non-woody environments for 10, 15, 20, 25, 30, and 35 days using an RT-qPCR analysis. The results reveal a significant difference in the expression levels of Tg-ABCs between woody and non-woody environments, suggesting an active involvement of the ABC gene family in wood degradation. During the wood degradation period of T. gibbosa, spanning from 10 to 35 days, the relative expression levels of most Tg-ABCs exhibited a trend of increasing, decreasing, and then increasing again. Additionally, at 20 and 35 days of wood degradation by T. gibbosa, the relative expression levels of Tg-ABCs peak, suggesting that at these time points, Tg-ABCs exert the most significant impact on the degradation of poplar wood by T. gibbosa. This study systematically reveals the biological characteristics of the ABC gene family in T. gibbosa and their response to woody environments. It establishes the foundation for a more profound comprehension of the wood-degradation mechanism of the ABC gene family and provides strong support for the development of more efficient wood-degradation strategies.

Role of bZIP transcription factors in the regulation of plant secondary metabolism

MAIN CONCLUSION

This study provides an overview of the structure, classification, regulatory mechanisms, and biological functions of the basic (region) leucine zipper transcription factors and their molecular mechanisms in flavonoid, terpenoid, alkaloid, phenolic acid, and lignin biosynthesis. Basic (region) leucine zippers (bZIPs) are evolutionarily conserved transcription factors (TFs) in eukaryotic organisms. The bZIP TFs are widely distributed in plants and play important roles in plant growth and development, photomorphogenesis, signal transduction, resistance to pathogenic microbes, biotic and abiotic stress, and secondary metabolism. Moreover, the expression of bZIP TFs not only promotes or inhibits the accumulation of secondary metabolites in medicinal plants, but also affects the stress response of plants to the external adverse environment. This paper describes the structure, classification, biological function, and regulatory mechanisms of bZIP TFs. In addition, the molecular mechanism of bZIP TFs regulating the biosynthesis of flavonoids, terpenoids, alkaloids, phenolic acids, and lignin are also elaborated. This review provides a summary for in-depth study of the molecular mechanism of bZIP TFs regulating the synthesis pathway of secondary metabolites and plant molecular breeding, which is of significance for the generation of beneficial secondary metabolites and the improvement of plant varieties.

High-throughput miRNA sequencing and identification of a novel ICE1-targeting miRNA in response to low temperature stress in Eucalyptus camaldulensis

MicroRNAs (miRNAs) play a crucial role in the growth, development, morphogenesis, signal transduction, and stress response in plants. The ICE (Inducer of CBF expression)-CBF (C-repeat binding factor)-COR (Cold-regulated gene) regulatory cascade is an important signalling pathway in plant response to low temperature stress, and it remains unknown whether this pathway is regulated by miRNAs. In this study, high-throughput sequencing was employed for predicting and identifying the miRNAs that were likely to target the ICE-CBF-COR pathway in Eucalyptus camaldulensis. A novel ICE1-targeting miRNA, eca-novel-miR-259-5p (nov-miR259), was further analysed. A total of 392 conserved miRNAs and 97 novel miRNAs were predicted, including 80 differentially expressed miRNAs. Of these, 30 miRNAs were predicted to be associated with the ICE-CBF-COR pathway. The full-length of mature nov-miR259 was 22 bp and its precursor gene was 60 bp in length, with a typical hairpin structure. The RNA ligase-mediated 5' amplification of cDNA ends (5'-RLM-RACE) and Agrobacterium-mediated tobacco transient expression assays demonstrated that nov-miR259 could cleave EcaICE1 in vivo. Moreover, qRT-PCR and Pearson's correlation analysis further revealed that the expression levels of nov-miR259 were almost significantly negatively correlated with those of its target gene, EcaICE1, and the other genes in the ICE-CBF-COR pathway. We first identified the nov-miR259 as a novel ICE1-targeting miRNA, and the nov-miR259-ICE1 module may be involved in regulating the cold stress response in E. camaldulensis.

Effects of differentially expressed microRNAs induced by rootstocks and silicon on improving chilling tolerance of cucumber seedlings (Cucumis sativus L.)

BACKGROUND

Rootstocks can improve the chilling tolerance of grafted cucumbers, but their effectiveness varies. Rootstocks with strong de-blooming capacity may result in lower chilling tolerance of grafted cucumbers compared to those with weak de-blooming capacity, while also reducing the silicon absorption. However, it remains unclear whether this reduction in chilling tolerance is due to differences in rootstock genotypes or the reduction in silicon absorption.

RESULTS

The chilling tolerance of cucumber seedlings was improved by using rootstocks and silicon nutrition. Rootstocks had a more significant effect than silicon nutrition, and the weak de-blooming rootstock 'Yunnan figleaf gourd' was superior to the strong de-blooming rootstock 'Huangchenggen No. 2'. Compared to self-rooted cucumber, twelve miRNAs were regulated by two rootstocks, including seven identical miRNAs (novel-mir23, novel-mir26, novel-mir30, novel-mir37, novel-mir46, miR395a and miR398a-3p) and five different miRNAs (novel-mir32, novel-mir38, novel-mir65, novel-mir78 and miR397a). Notably, four of these miRNAs (novel-mir38, novel-mir65, novel-mir78 and miR397a) were only identified in 'Yunnan figleaf gourd'-grafted cucumbers. Furthermore, six miRNAs (miR168a-5p, miR390a-5p, novel-mir26, novel-mir55, novel-mir67 and novel-mir70) were found to be responsive to exogenous silicon. Target gene prediction for 20 miRNAs resulted in 520 genes. Functional analysis of these target genes showed that 'Yunnan figleaf gourd' improves the chilling tolerance of cucumber by regulating laccase synthesis and sulfate metabolism, while 'Huangchenggen No. 2' and exogenous silicon reduced chilling stress damage to cucumber by regulating ROS scavenging and protein protection, respectively.

CONCLUSION

Among the identified miRNAs, novel-mir46 and miR398a-3p were found in cucumbers in response to chilling stress and two types of rootstocks. However, no identical miRNAs were identified in response to chilling stress and silicon. In addition, the differential expression of novel-mir38, novel-mir65, novel-mir78 and miR397a may be one of the important reasons for the differences in chilling tolerance of grafted cucumbers caused by two types of rootstocks.

Genome-Wide Analysis of the Expansin Gene Family in Populus and Characterization of Expression Changes in Response to Phytohormone (Abscisic Acid) and Abiotic (Low-Temperature) Stresses

Expansins are a group of cell wall enzyme proteins that help to loosen cell walls by breaking hydrogen bonds between cellulose microfibrils and hemicellulose. Expansins are essential plant proteins that are involved in several key processes, including seed germination, the growth of pollen tubes and root hairs, fruit ripening and abscission processes. Currently, there is a lack of knowledge concerning the role of expansins in woody plants. In this study, we analyzed expansin genes using Populus genome as the study target. Thirty-six members of the expansin gene family were identified in Populus that were divided into four subfamilies (EXPA, EXPB, EXLA and EXLB). We analyzed the molecular structure, chromosome localization, evolutionary relationships and tissue specificity of these genes and investigated expression changes in responses to phytohormone and abiotic stresses of the expansin genes of Populus tremula L. (PtEXs). Molecular structure analysis revealed that each PtEX protein had several conserved motifs and all of the PtEXs genes had multiple exons. Chromosome structure analysis showed that the expansin gene family is distributed on 14 chromosomes. The PtEXs gene family expansion patterns showed segmental duplication. Transcriptome data of Populus revealed that 36 PtEXs genes were differently expressed in different tissues. Cis-element analysis showed that the PtEXs were closely associated with plant development and responses to phytohormone and abiotic stress. Quantitative real-time PCR showed that abscisic acid (ABA) and low-temperature treatment affected the expression of some PtEXs genes, suggesting that these genes are involved in responses to phytohormone and abiotic stress. This study provides a further understanding of the expansin gene family in Populus and forms a basis for future functional research studies.

Symbiotic Modulation as a Driver of Niche Expansion of Coastal Plants in the San Juan Archipelago of Washington State

Modern evolutionary theory and population genetics posit that adaptation and habitat expansion of plants result from processes exclusive to their genomes. Here, we present studies showing that plants can grow across complex habitat gradients by modulating symbiotic associations with Class 2 fungal endophytes. Endophyte analysis of three native (Leymus mollis, Distichlis spicata, and Salicornia pacifica) and one invasive (Spartina anglica) plant growing across adjacent microhabitats in the San Juan Archipelago altered associations with Class 2 fungal endophytes in response to soil salinity levels. At the microhabitat interfaces where the gradation of salinity varied, the plants were colonized by endophytes from both microhabitats. A reciprocal transplant study along a salt gradient demonstrated that Leymus mollis (dunegrass) required endophytes indigenous to each microhabitat for optimal fitness and/or survival. In contrast, when dunegrass and Grindelia integrifolia (gumweed) were found growing in low salinity, but high drought habitats, these plant species had their own unique dominant endophyte association regardless of geographic proximity and conferred drought but not high salt stress tolerance. Modulation of endophyte abundance occurred in planta based on the ability of the symbiont to confer tolerance to the stress imposed on plants. The ability of an endophyte to confer appropriate stress tolerance resulted in a significant increase of in planta fungal abundance. Conversely, the inability of an endophyte to confer stress tolerance resulted in a decrease of in planta fungal abundance. Our studies indicate that Class 2 fungal endophytes can provide a symbiotic mechanism for niche expansion and phenotypic plasticity across environmental gradients.

Small RNA and Degradome Deep Sequencing Reveals the Roles of microRNAs in Peanut (Arachis hypogaea L.) Cold Response

Cold stress is a major environmental factor that affects plant growth and development, as well as fruit postharvest life and quality. MicroRNAs (miRNAs) are a class of non-coding small RNAs that play crucial roles in various abiotic stresses. Peanuts (Arachis hypogaea L.), one of the most important grain legumes and source of edible oils and proteins, are cultivated in the semi-arid tropical and subtropical regions of the world. To date, there has been no report on the role of miRNAs in the response to cold stress in cultivated peanuts. In this study, we profiled cold-responsive miRNAs in peanuts using deep sequencing in cold-sensitive (WQL20) alongside a tolerant variety (WQL30). A total of 407 known miRNAs and 143 novel peanut-specific miRNAs were identified. The expression of selected known and novel miRNAs was validated by northern blotting and six known cold-responsive miRNAs were revealed. Degradome sequencing identified six cold-responsive miRNAs that regulate 12 target genes. The correlative expression patterns of several miRNAs and their target genes were further validated using qRT-PCR. Our data showed that miR160-ARF, miR482-WDRL, miR2118-DR, miR396-GRF, miR162-DCL, miR1511-SRF, and miR1511-SPIRAL1 modules may mediate cold stress responses. Transient expression analysis in Nicotiana benthamiana found that miR160, miR482, and miR2118 may play positive roles, and miR396, miR162, and miR1511 play negative roles in the regulation of peanut cold tolerance. Our results provide a foundation for understanding miRNA-dependent cold stress response in peanuts. The characterized correlations between miRNAs and their response to cold stress could serve as markers in breeding programs or tools for improving cold tolerance of peanuts.

Study on the differences of phyllosphere microorganisms between poplar hybrid offspring and their parents

The females and males of dioecious plants have evolved sex-specific characteristics in terms of their morphological and physiological properties. However, the differentiation of phyllosphere microorganism of dioecious plants between parents and hybrid offspring remain largely unexplored. Here, the phyllosphere bacterial and fungal community diversity and composition of female (Populus nigra 'DH5' (PNDH5)), male (P. simonii 'DH4' (PSDH4)), and the hybrid offspring (P. simonii × P. nigra 'DH1' (PSPNDH1), P. simonii × P. nigra 'DH2' (PSPNDH2), P. simonii × P. nigra 'DH3' (PSPNDH3)) were investigated using 16S rDNA/ITS rDNA gene-based Illumina NovaSeq 6000 sequencing. There was considerable variation of plant height, diameter at breast height, leaf area, length of petioles, leaf moisture content, and starch among different samples, and PSDH2 owned the highest plant height, diameter at breast height, and length of petioles. No distinct differences of phyllosphere bacterial community diversity were observed among PSDH4, PNDH5, PSPNDH1, PSPNDH2, and PSPNDH3; while, PSPNDH2 owned the highest fungal Pielou_e index, Shannon index, and Simpson index. Firmicutes and Ascomycota were the predominant phyllosphere bacterial and fungal community at the phylum level, respectively. Bacilli and Gammaproteobacteria were the two most dominant bacterial classes regardless of parent and the hybrid offspring. The predominant phyllosphere fungal community was Dothideomycetes at the class level. The NMDS demonstrated that phyllosphere microbial community obviously differed between parents and offspring, while the phyllosphere microbial community presented some similarities under different hybrid progeny. Also, leaf characteristics contributed to the differentiation of phyllosphere bacterial and fungal communities between parents and hybrid offspring. These results highlighted the discrimination of phyllosphere microorganisms on parent and hybrid offspring, which provided clues to potential host-related species in the phyllosphere environment.

Fine-Tuning Cold Stress Response Through Regulated Cellular Abundance and Mechanistic Actions of Transcription Factors

Inflictions caused by cold stress can result in disastrous effects on the productivity and survival of plants. Cold stress response in plants requires crosstalk between multiple signaling pathways including cold, heat, and reactive oxygen species (ROS) signaling networks. CBF, MYB, bHLH, and WRKY families are among the TFs that function as key players in the regulation of cold stress response at the molecular level. This review discusses some of the latest understanding on the regulation of expression and the mechanistic actions of plant TFs to address cold stress response. It was shown that the plant response consists of early and late responses as well as memory reprogramming for long-term protection against cold stress. The regulatory network can be differentiated into CBF-dependent and independent pathways involving different sets of TFs. Post-transcriptional regulation by miRNAs, control during ribosomal translation process, and post-translational regulation involving 26S proteosomic degradation are processes that affect the cellular abundance of key regulatory TFs, which is an important aspect of the regulation for cold acclimation. Therefore, fine-tuning of the regulation by TFs for adjusting to the cold stress condition involving the dynamic action of protein kinases, membrane ion channels, adapters, and modifiers is emphasized in this review.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Identification of MicroRNAs in Taxillus chinensis (DC.) Danser Seeds under Cold Stress

Taxillus chinensis (DC.) Danser, a parasitic plant that belongs to the Loranthaceae family, has a long history of being used in the Chinese medicine. We observed that the loranthus seeds were sensitive to temperature and could lose viability below 0°C quickly. Thus, we performed small RNA sequencing to study the microRNA (miRNA) regulation in the loranthus seeds under cold stress. In total, we identified 600 miRNAs, for the first time, in the loranthus seeds under cold stress. Then, we detected 224, 229, and 223 miRNAs (TPM > 1) in A0 (control), A1 (cold treatment for 12 h at 0°C), and A2 (cold treatment for 36 h at 0°C), respectively. We next identified 103 differentially expressed miRNAs (DEmiRs) in the loranthus seeds in response to cold. Notably, miR408 was upregulated during the cold treatment, which can regulate genes encoding phytocyanin family proteins and phytophenol oxidases. Some DEmiRs were specific to A1 and may function in early response to cold, such as gma-miR393b-3p, miR946, ath-miR779.2-3p, miR398, and miR9662. It is interesting that ICE3, IAA13, and multiple transcription factors (e.g., WRKY and CRF4 and TCP4) regulated by the DEmiRs have been reported to respond cold in other plants. We further identified 4, 3, and 4 DEmiRs involved in the pathways “responding to cold,” “responding to abiotic stimulus,” and “seed development/germination,” respectively. qRT-PCR was used to confirm the expression changes of DEmiRs and their targets in the loranthus seeds during the cold treatment. This is the first time to study cold-responsive miRNAs in loranthus, and our findings provide a valuable resource for future studies.

Variation and Genetic Parameters of Leaf Morphological Traits of Eight Families from Populus simonii × P. nigra

Leaf morphology in Populus L. varies extensively among sections, species and clones under strong genetic control. P. nigra L. (section Aigeiros), with large and triangular leaves, is a commercial forest tree of economic importance for fast growth and high yield in Europe. P. simonii Carr. (section Tacamahaca) with small land rhomboid ovate leaves performs cold and dry resistance/tolerance in the semi-arid region of Northern China. Leaf morphological traits could be used as early indicators to improve the efficiency of selection. In order to investigate the genetic variation pattern of leaf morphology traits, estimate breeding values (combining ability), as well as evaluate crossing combinations of parents, 1872 intersectional progenies from eight families (P. simonii × P. nigra) and their parents were planted with cuttings for the clonal replicate field trial in Northern China. Four leaf size traits (area, perimeter, length, width) and roundness were measured with leaf samples from the 1-year-old clonal plantation. Significant differences regarding leaf traits were found between and among three female clones of P. simonii from Inner Mongolia, China and six male clones of P. nigra from Casale Monferrato, Italy. The genetic variation coefficient, heritability and genetic variance component of most traits in male parents were greater than these of female parents. Heritability estimates of male and female parents were above 0.56 and 0.17, respectively. Plentiful leaf variations with normal and continuous distributions exited in the hybrid progenies among and within families with the genetic variation coefficient and heritability above 28.49 and 0.24, respectively. Heritability estimates showed that leaf area was the most heritable trait, followed by leaf width. The breeding value ranking of parents allowed us to select the parental clones for new crosses and extend the mating design. Two male parental clones (N430 and N429) had greater breeding values (general combining ability, GCA) of leaf size traits than other clones. The special combining ability (SCA) of the crossing combination between P. simonii cl. ZL-3 and P. nigra cl. N430 was greater than that of others. Eight putatively superior genotypes, most combined with the female parental clone ZL-3, can be selected for future testing under near-commercial conditions. Significant genetic and phenotypic correlations were found between five leaf morphology traits with the coefficients above 0.9, except for leaf roundness. The results showed that leaf morphology traits were under strong genetic control and the parental clones with high GCA and SCA effects could be utilized in heterosis breeding, which will provide a starting point for devising a new selection strategy of parents and progenies.

Underpinning miRNA-miRNA co-functional interaction patterns in the metabolism of Oryza sativa by genome-scale network analysis

MicroRNA (miRNA) is a class of non-coding small RNAs, which post-transcriptionally regulate a large number of genes and are now known to be important regulators in a wide variety of biological processes including metabolism. Thus, for better understanding these complex biological networks, and to derive their significance and inter-dependency, a systems biology approach enables us to explore and draw vital insights into these molecular network architectures. In this study, we aimed to understand the significance of synergistic miRNA-miRNA interactions in rice by constructing and analysing metabolic networks. The construction of the network involves target gene prediction of experimentally verified miRNAs of rice and then appending associated metabolic pathways to the network. A genome-scale miRNA-miRNA co-functional network (MFSN) is constructed based on co-regulatory interactions among the miRNAs and common target genes by applying transformational procedures. The analysis of the extracted MFSN modules identifies co-regulated target genes that are associated with corresponding interconnected metabolic pathways such as VALDEG-PWY (L-valine degradation I pathway was found to be targeted by multiple miRNA families, such as osa-miR812, osa-miR818, osa-miR821, and osa-miR5799 families while another pathway that was found to be associated with multiple miRNA families was PWY-6952 (glycerophosphodiester degradation pathway), PWY-6952 was found to be targeted by osa-miR812, osa-miR11344 and osa-miR5801 families. Such extensive study will help in systematically elucidating the regulatory networks in metabolism of rice, which in turn can be utilised to devise strategies for crop improvement and novel cultivar development.