Morphological, phenological, and transcriptional analyses provide insight into the diverse flowering traits of a mutant of the relic woody plant Liriodendron chinense

Genomic survey and expression analysis of LcARFs reveal multiple functions to somatic embryogenesis in Liriodendron

BACKGROUND

Auxin response factors (ARFs) are critical transcription factors that mediate the auxin signaling pathway and are essential for regulating plant growth. However, there is a lack of understanding regarding the ARF gene family in Liriodendron chinense, a vital species in landscaping and economics. Thus, further research is needed to explore the roles of ARFs in L. chinense and their potential applications in plant development.

RESULT

In this study, we have identified 20 LcARF genes that belong to three subfamilies in the genome of L. chinense. The analysis of their conserved domains, gene structure, and phylogeny suggests that LcARFs may be evolutionarily conserved and functionally similar to other plant ARFs. The expression of LcARFs varies in different tissues. Additionally, they are also involved in different developmental stages of somatic embryogenesis. Overexpression of LcARF1, LcARF2a, and LcARF5 led to increased activity within callus. Additionally, our promoter-GFP fusion study indicated that LcARF1 may play a role in embryogenesis. Overall, this study provides insights into the functions of LcARFs in plant development and embryogenesis, which could facilitate the improvement of somatic embryogenesis in L. chinense.

CONCLUSION

The research findings presented in this study shed light on the regulatory roles of LcARFs in somatic embryogenesis in L. chinense and may aid in accelerating the breeding process of this tree species. By identifying the specific LcARFs involved in different stages of somatic embryogenesis, this study provides a basis for developing targeted breeding strategies aimed at optimizing somatic embryogenesis in L. chinense, which holds great potential for improving the growth and productivity of this economically important species.

Hydrangea arborescens 'Annabelle' Flower Formation and Flowering in the Current Year

The perennial woody plant Hydrangea arborescens 'Annabelle' is of great research value due to its unique mechanism of flower development that occurs in the current year, resulting in decorative flowers that can be enjoyed for a relatively long period of time. However, the mechanisms underlying the regulation of current-year flower development in H. arborescens 'Annabelle' are still not fully understood. In this study, we conducted an associated analysis to explore the core regulating network in H. arborescens 'Annabelle' by combining phenological observations, physiological assays, and transcriptome comparisons across seven flower developmental stages. Through this analysis, we constructed a gene co-expression network (GCN) based on the highest reciprocal rank (HRR), using 509 differentially expressed genes (DEGs) identified from seven flowering-related pathways, as well as the biosynthesis of eight flowering-related phytohormones and signal transduction in the transcriptomic analysis. According to the analysis of the GCN, we identified 14 key genes with the highest functional connectivity that played critical roles in specific development stages. We confirmed that 135 transcription factors (AP2/ERF, bHLH, CO-like, GRAS, MIKC, SBP, WRKY) were highly co-expressed with the 14 key genes, indicating their close associations with the development of current-year flowers. We further proposed a hypothetical model of a gene regulatory network for the development of the whole flower. This model suggested that the photoperiod, aging, and gibberellin pathways, along with the phytohormones abscisic acid (ABA), gibberellin (GA), brassinosteroid (BR), and jasmonic acid (JA), work synergistically to promote the floral transition. Additionally, auxin, GA, JA, ABA, and salicylic acid (SA) regulated the blooming process by involving the circadian clock. Cytokinin (CTK), ethylene (ETH), and SA were key regulators that affected flower senescence. Additionally, several floral integrators (HaLFY, HaSOC1-2, HaAP1, HaFULL, HaAGL24, HaFLC, etc.) were dominant contributors to the development of H. arborescens flowers. Overall, this research provides a comprehensive understanding of the dynamic mechanism underlying the entire process of current-year flower development, thereby offering valuable insights for further studies on the flower development of H. arborescens 'Annabelle'.

Genome-wide identification of the NAC gene family and its functional analysis in Liriodendron

As one of the largest plant specific transcription factor families, NAC family members play an important role in plant growth, development and stress resistance. To investigate the function of NAC transcription factors during abiotic stress, as well as during somatic embryogenesis, we identified and characterized the NAC gene family in Liriodendron chinense. We found that most LcNAC members contain more than three exons, with a relatively conserved gene and motif structure, especially at the N-terminus. Interspecies collinearity analysis revealed a closer relationship between the L. chinense NACs and the P. trichocarpa NACs. We analyzed the expression of LcNAC in different tissues and under three abiotic stresses. We found that 12 genes were highly expressed during the ES3 and ES4 stages of somatic embryos, suggesting that they are involved in the development of somatic embryos. 6 LcNAC genes are highly expressed in flower organs. The expression pattern analysis of LcNACs based on transcriptome data and RT-qPCR obtained from L. chinense leaves indicated differential expression responses to drought, cold, and heat stress. Genes in the NAM subfamily expressed differently during abiotic stress, and LcNAC6/18/41/65 might be the key genes in response to abiotic stress. LcNAC6/18/41/65 were cloned and transiently transformed into Liriodendron protoplasts, where LcNAC18/65 was localized in cytoplasm and nucleus, and LcNAC6/41 was localized only in nucleus. Overall, our findings suggest a role of the NAC gene family during environmental stresses in L. chinense. This research provides a basis for further study of NAC genes in Liriodendron chinense.

Characterization of the Molecular Events Underlying the Establishment of Axillary Meristem Region in Pepper

Plant architecture is a major motif of plant diversity, and shoot branching patterns primarily determine the aerial architecture of plants. In this study, we identified an inbred pepper line with fewer lateral branches, 20C1734, which was free of lateral branches at the middle and upper nodes of the main stem with smooth and flat leaf axils. Successive leaf axil sections confirmed that in normal pepper plants, for either node n, Pn (Primordium n) < 1 cm and Pn+1 < 1 cm were the critical periods between the identification of axillary meristems and the establishment of the region, whereas Pn+3 < 1 cm was fully developed and formed a completely new organ. In 20C1734, the normal axillary meristematic tissue region establishment and meristematic cell identity confirmation could not be performed on the axils without axillary buds. Comparative transcriptome analysis revealed that "auxin-activated signaling pathway", "response to auxin", "response to abscisic acid", "auxin biosynthetic process", and the biosynthesis of the terms/pathways, such as "secondary metabolites", were differentially enriched in different types of leaf axils at critical periods of axillary meristem development. The accuracy of RNA-seq was verified using RT-PCR for some genes in the pathway. Several differentially expressed genes (DEGs) related to endogenous phytohormones were targeted, including several genes of the PINs family. The endogenous hormone assay showed extremely high levels of IAA and ABA in leaf axils without axillary buds. ABA content in particular was unusually high. At the same time, there is no regular change in IAA level in this type of leaf axils (normal leaf axils will be accompanied by AM formation and IAA content will be low). Based on this, we speculated that the contents of endogenous hormones IAA and ABA in 20C1734 plant increased sharply, which led to the abnormal expression of genes in related pathways, which affected the formation of Ams in leaf axils in the middle and late vegetative growth period, and finally, nodes without axillary buds and side branches appeared.

Characterization of PsmiR319 during flower development in early- and late-flowering tree peonies cultivars

The flowering period is the most important ornamental trait of tree peony, while industrial development of tree peony has been limited by short flowering period. miR319 plays an important regulatory role in plant flowering. In the current study, the expression characteristics and evolution of PsmiR319 in tree peony flowering was explored using 'Feng Dan' and 'Lian He', which are early-flowering and late-flowering varieties of tree peony, respectively. The structure, evolution, and target(s) of PsmiR319 were analyzed by bioinformatics. Evolution analysis showed that pre-PsmiR319 was distributed in 41 plant species, among which the length of the precursor sequence exhibited marked differences (between 52 and 308 bp). Pre-PsmiR319 of tree peony was located close to the corresponding sequences of Linum usitatissimum and Picea abies in the phylogenetic tree, and in addition, could form a typical hairpin structure including a mature body with a length of 20 bp located on the 3p arm and part of the loop sequence. The mature sequence of miR319 was highly conserved among different species. Target genes of PsmiR319 include MYB-related transcription factor in tree peony. Expression of PsmiR319, assayed by qRT-PCR, differed between 'Feng Dan' and 'Lian He' during different flower development periods. PsmiR319 and its target gene showed a negative expression regulation relationship during the periods of CE (color exposure), BS (blooming stage), IF (initial flowering), and HO (half opening) in the early-flowering 'Feng Dan', and the same in FB (Full blooming) periods of late-flowering 'Lian He'. Findings from this study provide a reference for further investigation into the mechanism of miR319 in the development of different varieties of tree peony.

Genetic studies on continuous flowering in woody plant Osmanthus fragrans

Continuous flowering is a key horticultural trait in ornamental plants, whereas the specific molecular regulation mechanism remains largely unknown. In sweet osmanthus (Osmanthus fragrans Lour.), plants based on their flowering characteristics are divided into once-flowering (OF) habit and continuous flowering (CF) habit. Here, we first described the flowering phenology shifts of OF and CF habits in sweet osmanthus through paraffin section and microscope assay. Phenotypic characterization showed that CF plants had constant new shoot growth, floral transition, and blooming for 1 year, which might lead to a continuous flowering trait. We performed the transcriptome sequencing of OF and CF sweet osmanthus and analyzed the transcriptional activity of flowering-related genes. Among the genes, three floral integrators, OfFT, OfTFL1, and OfBFT, had a differential expression during the floral transition process in OF and CF habits. The expression patterns of the three genes in 1 year were revealed. The results suggested that their accumulations corresponded to the new shoots occurring and the floral transition process. Function studies suggested that OfFT acted as a flowering activator, whereas OfBFT was a flowering inhibitor. Yeast one-hybrid assay indicated that OfSPL8 was a common upstream transcription factor of OfFT and OfBFT, suggesting the vital role of OfSPL8 in continuous flowering regulation. These results provide a novel insight into the molecular mechanism of continuous flowering.

From Floral Induction to Blooming: The Molecular Mysteries of Flowering in Woody Plants

Flowering is a pivotal developmental process in response to the environment and determines the start of a new life cycle in plants. Woody plants usually possess a long juvenile nonflowering phase followed by an adult phase with repeated flowering cycles. The molecular mechanism underlying flowering regulation in woody plants is believed to be much more complex than that in annual herbs. In this review, we briefly describe the successive but distinct flowering processes in perennial trees, namely the vegetative phase change, the floral transition, floral organogenesis, and final blooming, and summarize in detail the most recent advances in understanding how woody plants regulate flowering through dynamic gene expression. Notably, the florigen gene FLOWERING LOCUS T(FT) and its antagonistic gene TERMINAL FLOWER 1 (TFL1) seem to play a central role in various flowering transition events. Flower development in different taxa requires interactions between floral homeotic genes together with AGL6 conferring floral organ identity. Finally, we illustrate the issues and corresponding measures of flowering regulation investigation. It is of great benefit to the future study of flowering in perennial trees.

Genome-Wide Identification and Expression Analysis of SnRK2 Gene Family in Dormant Vegetative Buds of Liriodendron chinense in Response to Abscisic Acid, Chilling, and Photoperiod

Protein kinases play an essential role in plants’ responses to environmental stress signals. SnRK2 (sucrose non-fermenting 1-related protein kinase 2) is a plant-specific protein kinase that plays a crucial role in abscisic acid and abiotic stress responses in some model plant species. In apple, corn, rice, pepper, grapevine, Arabidopsis thaliana, potato, and tomato, a genome-wide study of the SnRK2 protein family was performed earlier. The genome-wide comprehensive investigation was first revealed to categorize the SnRK2 genes in the Liriodendron chinense (L. chinense). The five SnRK2 genes found in the L. chinense genome were highlighted in this study. The structural gene variants, 3D structure, chromosomal distributions, motif analysis, phylogeny, subcellular localization, cis-regulatory elements, expression profiles in dormant buds, and photoperiod and chilling responses were all investigated in this research. The five SnRK2 genes from L. chinense were grouped into groups (I–IV) based on phylogeny analysis, with three being closely related to other species. Five hormones-, six stress-, two growths and biological process-, and two metabolic-related responsive elements were discovered by studying the cis-elements in the promoters. According to the expression analyses, all five genes were up- and down-regulated in response to abscisic acid (ABA), photoperiod, chilling, and chilling, as well as photoperiod treatments. Our findings gave insight into the SnRK2 family genes in L. chinense and opened up new study options.

Genome-wide identification of the Liriodendron chinense WRKY gene family and its diverse roles in response to multiple abiotic stress

BACKGROUND

Liriodendron chinense (Lchi) is a tree species within the Magnoliaceae family and is considered a basal angiosperm. The too low or high temperature or soil drought will restrict its growth as the adverse environmental conditions, thus improving L. chinense abiotic tolerance was the key issues to study. WRKYs are a major family of plant transcription factors known to often be involved in biotic and abiotic stress responses. So far, it is still largely unknown if and how the LchiWRKY gene family is tied to regulating L. chinense stress responses. Therefore, studying the involvement of the WRKY gene family in abiotic stress regulation in L. chinense could be very informative in showing how this tree deals with such stressful conditions.

RESULTS

In this research, we performed a genome-wide analysis of the Liriodendron chinense (Lchi) WRKY gene family, studying their classification relationships, gene structure, chromosomal locations, gene duplication, cis-element, and response to abiotic stress. The 44 members of the LchiWRKY gene family contain a significant amount of sequence diversity, with their lengths ranging from 525 bp to 40,981 bp. Using classification analysis, we divided the 44 LchiWRKY genes into three phylogenetic groups (I, II, II), with group II then being further divided into five subgroups (IIa, IIb, IIc, IId, IIe). Comparative phylogenetic analysis including the WRKY families from 17 plant species suggested that LchiWRKYs are closely related to the Magnolia Cinnamomum kanehirae WRKY family, and has fewer family members than higher plants. We found the LchiWRKYs to be evenly distributed across 15 chromosomes, with their duplication events suggesting that tandem duplication may have played a major role in LchiWRKY gene expansion model. A Ka/Ks analysis indicated that they mainly underwent purifying selection and distributed in the group IId. Motif analysis showed that LchiWRKYs contained 20 motifs, and different phylogenetic groups contained conserved motif. Gene ontology (GO) analysis showed that LchiWRKYs were mainly enriched in two categories, i.e., biological process and molecular function. Two group IIc members (LchiWRKY10 and LchiWRKY37) contain unique WRKY element sequence variants (WRKYGKK and WRKYGKS). Gene structure analysis showed that most LchiWRKYs possess 3 exons and two different types of introns: the R- and V-type which are both contained within the WRKY domain (WD). Additional promoter cis-element analysis indicated that 12 cis-elements that play different functions in environmental adaptability occur across all LchiWRKY groups. Heat, cold, and drought stress mainly induced the expression of group II and I LchiWRKYs, some of which had undergone gene duplication during evolution, and more than half of which had three exons. LchiWRKY33 mainly responded to cold stress and LchiWRKY25 mainly responded to heat stress, and LchiWRKY18 mainly responded to drought stress, which was almost 4-fold highly expressed, while 5 LchiWRKYs (LchiWRKY5, LchiWRKY23, LchiWRKY14, LchiWRKY27, and LchiWRKY36) responded equally three stresses with more than 6-fold expression. Subcellular localization analysis showed that all LchiWRKYs were localized in the nucleus, and subcellular localization experiments of LchiWRKY18 and 36 also showed that these two transcription factors were expressed in the nucleus.

CONCLUSIONS

This study shows that in Liriodendron chinense, several WRKY genes like LchiWRKY33, LchiWRKY25, and LchiWRKY18, respond to cold or heat or drought stress, suggesting that they may indeed play a role in regulating the tree's response to such conditions. This information will prove a pivotal role in directing further studies on the function of the LchiWRKY gene family in abiotic stress response and provides a theoretical basis for popularizing afforestation in different regions of China.

Insight Into the Multiple Branches Traits of a Mutant in Larix olgensis by Morphological, Cytological, and Transcriptional Analyses

Larix olgensis is a tall deciduous tree species that has many applications in the wood fiber industry. Bud mutations are somatic mutations in plants and are considered an ideal material to identify and describe the molecular mechanism of plant mutation. However, the molecular regulatory mechanisms of bud mutations in L. olgensis remain unknown. In this study, dwarfed (or stunted), short-leaved, and multi-branched mutants of L. olgensis were found and utilized to identify crucial genes and regulatory networks controlling the multiple branch structure of L. olgensis. The physiological data showed that the branch number, bud number, fresh and dry weight, tracheid length, tracheid length-width ratio, inner tracheid diameter, and epidermal cell area of mutant plants were higher than that of wild-type plants. Hormone concentration measurements found that auxin, gibberellin, and abscisic acid in the mutant leaves were higher than that in wild-type plants. Moreover, the transcriptome sequencing of all samples using the Illumina Hiseq sequencing platform. Transcriptome analysis identified, respectively, 632, 157, and 199 differentially expressed genes (DEGs) in buds, leaves, and stems between mutant plants and wild type. DEGs were found to be involved in cell division and differentiation, shoot apical meristem activity, plant hormone biosynthesis, and sugar metabolism. Furthermore, bZIP, WRKY, and AP2/ERF family transcription factors play a role in bud formation. This study provides new insights into the molecular mechanisms of L. olgensis bud and branch formation and establishes a fundamental understanding of the breeding of new varieties in L. olgensis.