Transcriptomic and microstructural analyses in Liriodendron tulipifera Linn. reveal candidate genes involved in nectary development and nectar secretion

Predicting the impact of climate change on the distribution of two relict Liriodendron species by coupling the MaxEnt model and actual physiological indicators in relation to stress tolerance

Climate change has a crucial impact on the distributions of plants, especially relict species. Hence, predicting the potential impact of climate change on the distributions of relict plants is critical for their future conservation. Liriodendron plants are relict trees, and only two natural species have survived: L. chinense and L. tulipifera. However, the extent of the impact of future climate change on the distributions of these two Liriodendron species remains unclear. Therefore, we predicted the suitable habitat distributions of two Liriodendron species under present and future climate scenarios using MaxEnt modeling. The results showed that the area of suitable habitats for two Liriodendron species would significantly decrease. However, the two relict species presented different habitat shift patterns, with a local contraction of suitable habitat for L. chinense and a northward shift in suitable habitat for L. tulipifera, indicating that changes in environmental factors will affect the distributions of these species. Among the environmental factors assessed, May precipitation induced the largest impact on the L. chinense distribution, while L. tulipifera was significantly affected by precipitation in the driest quarter. Furthermore, to explore the relationship between habitat suitability and Liriodendron stress tolerance, we analyzed six physiological indicators of stress tolerance by sampling twelve provenances of L. chinense and five provenances of L. tulipifera. The composite index of six physiological indicators was significantly negatively correlated with the habitat suitability of the species. The stress tolerance of Liriodendron plants in highly suitable areas was lower than that in areas with moderate or low suitability. Overall, these findings improve our understanding of the ecological impacts of climate change, informing future conservation efforts for Liriodendron species.

Identification, Phylogenetic and Expression Analyses of the AAAP Gene Family in Liriodendron chinense Reveal Their Putative Functions in Response to Organ and Multiple Abiotic Stresses

In this study, 52 AAAP genes were identified in the L. chinense genome and divided into eight subgroups based on phylogenetic relationships, gene structure, and conserved motif. A total of 48 LcAAAP genes were located on the 14 chromosomes, and the remaining four genes were mapped in the contigs. Multispecies phylogenetic tree and codon usage bias analysis show that the LcAAAP gene family is closer to the AAAP of Amborella trichopoda, indicating that the LcAAAP gene family is relatively primitive in angiosperms. Gene duplication events revealed six pairs of segmental duplications and one pair of tandem duplications, in which many paralogous genes diverged in function before monocotyledonous and dicotyledonous plants differentiation and were strongly purification selected. Gene expression pattern analysis showed that the LcAAAP gene plays a certain role in the development of Liriodendron nectary and somatic embryogenesis. Low temperature, drought, and heat stresses may activate some WRKY/MYB transcription factors to positively regulate the expression of LcAAAP genes to achieve long-distance transport of amino acids in plants to resist the unfavorable external environment. In addition, the GAT and PorT subgroups could involve gamma-aminobutyric acid (GABA) transport under aluminum poisoning. These findings could lay a solid foundation for further study of the biological role of LcAAAP and improvement of the stress resistance of Liriodendron.

Identification and Expression Analysis of R2R3-MYB Family Genes Associated with Salt Tolerance in Cyclocarya paliurus

R2R3-MYB transcription factors are most abundant in the MYB superfamily, while the R2R3-MYB genes play an important role in plant growth and development, especially in response to environmental stress. Cyclocarya paliurus is a multifunction tree species, and the existing resources cannot meet the requirement for its leaf production and medical use. Therefore, lands with some environmental stresses would be potential sites for developing C. paliurus plantations. However, the function of R2R3-MYB genes in C.paliurus in response to environmental stress remains unknown. In this study, to identify the roles of R2R3-MYB genes associated with salt stress response, 153 CpaMYB genes and their corresponding protein sequences were identified from the full-length transcriptome. Based on the comparison with MYB protein sequences of Arabidopsis thaliana, 69 R2R3-MYB proteins in C. paliurus were extracted for further screening combined with conserved functional domains. Furthermore, the MYB family members were analyzed from the aspects of protein sequences alignment, evolution, motif prediction, promoter cis-acting element analysis, and gene differential expression under different salt treatments using both a pot experiment and hydroponic experiment. The results showed that the R2R3-MYB genes of C.paliurus conserved functional domains, whereas four R2R3-MYB genes that might respond to salt stress via regulating plant hormone signals were identified in this study. This work provides a basis for further functional characterization of R2R3-MYB TFs in C. paliurus.

Transcriptomic profiling of the floral fragrance biosynthesis pathway of Liriodendron and functional characterization of the LtuDXR gene

Floral fragrance, which has the function of attracting pollinators, is a class of volatile secondary metabolites mainly released by the secretory tissue of petals. Terpenoids are key components of floral volatile substances. Previous studies have shown that there are significant differences in the concentration and composition of volatile floral fragrances, especially terpenoids, between Liriodendron chinense and L. tulipifera. At present, the mechanism by which the synthesis of floral fragrance is regulated in Liriodendron remains unexplored. In this study, we analyzed the differentially expressed genes (DEGs) of L. chinense and L. tulipifera, and identified 130 DEGs related to terpenoid synthesis. A KEGG enrichment analysis of DEGs related to terpenoid biosynthesis revealed that the monoterpenoid biosynthesis pathway was the most significant. We cloned the LtuDXR gene from L. tulipifera using RACE technology. RT-qPCR results showed that the expression of the LtuDXR gene was the highest in the early florescence petals, indicating that the LtuDXR gene may play a role in the synthesis of volatile terpenoids. Subcellular localization showed that the LtuDXR protein is mainly localized in the chloroplast. Overexpression of LtuDXR in Arabidopsis thaliana significantly increased the plant height, DXR enzyme activity, and carotenoid content. In this study, we identified and functionally characterized LtuDXR, which is involved in terpenoid synthesis in Liriodendron. Our work lays the foundation for further exploration of the molecular mechanism by which terpenoid biosynthesis is regulated in Liriodendron.

Genome-Wide Identification and Expression Analysis of R2R3-MYB Family Genes Associated with Petal Pigment Synthesis in Liriodendron

The MYB transcription factor family is one of the largest families in plants, and its members have various biological functions. R2R3-MYB genes are involved in the synthesis of pigments that yield petal colors. Liriodendron plants are widely cultivated as ornamental trees owing to their peculiar leaves, tulip-like flowers, and colorful petals. However, the mechanism underlying petal coloring in this species is unknown, and minimal information about MYB genes in Liriodendron is available. Herein, this study aimed to discern gene(s) involved in petal coloration in Liriodendron via genome-wide identification, HPLC, and RT-qPCR assays. In total, 204 LcMYB superfamily genes were identified in the Liriodendron chinense genome, and 85 R2R3-MYB genes were mapped onto 19 chromosomes. Chromosome 4 contained the most (10) R2R3-MYB genes, and chromosomes 14 and 16 contained the fewest (only one). MEME analysis showed that R2R3-MYB proteins in L. chinense were highly conserved and that their exon-intron structures varied. The HPLC results showed that three major carotenoids were uniformly distributed in the petals of L. chinense, while lycopene and β-carotene were concentrated in the orange band region in the petals of Liriodendron tulipifera. Furthermore, the expression profiles via RT-qPCR assays revealed that four R2R3-MYB genes were expressed at the highest levels at the S3P/S4P stage in L. tulipifera. This result combined with the HPLC results showed that these four R2R3-MYB genes might participate in carotenoid synthesis in the petals of L. tulipifera. This work laid a cornerstone for further functional characterization of R2R3-MYB genes in Liriodendron plants.

Isolation, expression, and functional analysis of the geranylgeranyl pyrophosphate synthase (GGPPS) gene from Liriodendron tulipifera

Terpenoids are important secondary metabolites in plants and are involved in stress responses and pollinator attraction. Geranylgeranyl pyrophosphate synthase (GGPPS) is a key synthase in the 2C-methyl-D-erythritol-4-phosphate (MEP) pathway of terpenoid synthesis, catalyzing the synthesis of diterpenoids. Liriodendron tulipifera is a nectar plant in North America. Little is known about the key genes involved in the biosynthetic pathways of terpenoids, the precursors of most compounds related to nectar, fragrance and coloring in flowers in L. tulipifera. In this study, the LtuGGPPS2 gene and its promoter (LtuGGPPS2-pro) were cloned from L. tulipifera. The results of sequence alignment showed that the LtuGGPPS2 gene is highly homologous to GGPPS genes of other plants. Subcellular localization analysis showed that the LtuGGPPS2 protein localizes to chloroplasts, suggesting that the LtuGGPPS2 gene is probably related to carotenoid and chlorophyll synthesis. Based on tissue expression profiles revealed by RT-qPCR, the expression level of the LtuGGPPS2 gene was highest in petals. These results were consistent with the changes in volatile and nonvolatile terpenoids in the flowers of L. tulipifera. GUS staining to examine the LtuGGPPS2 promoter indicated that it is responsive to hormones. Overexpression of the LtuGGPPS2 gene increased the carotenoid content and GGPPS enzyme activity in Arabidopsis thaliana, indicating that LtuGGPPS2 is the key terpenoid synthase in the flowers of L. tulipifera. Our findings lay a foundation for further functional analysis of the LtuGGPPS2 gene and deeper investigation of the terpenoid biosynthetic pathway in L. tulipifera.

Integrative analysis of transcriptome and proteome revealed nectary and nectar traits in the plant-pollinator interaction of Nitraria tangutorum Bobrov

BACKGROUND

Nitraria tangutorum is an important desert shrub that shows resistance to drought, salt and wind erosion stresses. It is a central ecological species in its area. Here, we have studied how N. tangutorum has adapted to achieve a successful reproduction strategy.

RESULTS

We found that N. tangutorum is mainly pollinated by insects of the Hymenoptera, Diptera and Coleoptera orders. Nitraria tangutorum has very small flowers, with the nectary composed of secretive epidermal cells from which nectar is secreted, located within the inner petals. In addition, analyzing the transcriptome of four successive flower developmental stages revealed that mainly differentially expressed genes associated with flower and nectary development, nectar biosynthesis and secretion, flavonoid biosynthesis, plant hormone signal transduction and plant-pathogen interaction show dynamic expression. From the nectar, we could identify seven important proteins, of which the L-ascorbate oxidase protein was first found in plant nectar. Based on the physiological functions of these proteins, we predict that floral nectar proteins of N. tangutorum play an important role in defending against microbial infestation and scavenging active oxygen.

CONCLUSIONS

This study revealed that N. tangutorum is an insect-pollinated plant and its nectary is composed of secretive epidermal cells that specialized into secretive trichomes. We identified a large number of differentially expressed genes controlling flower and nectary development, nectar biosynthesis and secretion, flavonoid biosynthesis, plant hormone signal transduction and plant-pathogen interaction. We suggest that proteins present in N. tangutorum nectar may have both an antibacterial and oxygen scavenging effect. These results provide a scientific basis for exploring how the reproductive system of N. tangutorum and other arid-desert plants functions.

The PIN gene family in relic plant L. chinense: Genome-wide identification and gene expression profiling in different organizations and abiotic stress responses

The auxin efflux carrier PIN-FORMED (PIN) proteins are required for the polar transport of auxin between cells through their asymmetric distribution on the plasma membrane, thus mediating the differential distribution of auxin in plants, finally, affecting plant growth and developmental processes. In this study, 11 LcPIN genes were identified. The structural characteristics and evolutionary status of LcPIN genes were thoroughly investigated and interpreted combining physicochemical property analysis, evolutionary analysis, gene structure analysis, chromosomal localization, etc. Multi-species protein sequence analysis showed that angiosperm PIN genes have strong purification options and some functional sites were predicted about PIN protein polarity, trafficking and activity in L. chinense. Further qRT-PCR and transcriptome data analysis indicated that the long LcPINs have highly expressed from globular embryo to plantlet, and the LcPIN6a started upregulated in cotyledon embryo. The LcPIN3 and LcPIN6a are both highly expressed during the development of stamens and petals and the expression of LcPIN2 is related to root elongation, suggesting that they may play an important role in these processes. Experiment data indicates that LcPIN5 and LcPIN8 might play a key role in auxin transport in Liriodendron stems and leaves under abiotic stress. Analyzed the response of LcPIN genes to abiotic stress and as a basis for uncovering the biological role of LcPIN genes in development and adaption to adverse environments. This study provides a foundation for further genetic and functional analyses.

Genome-wide identification of MIKC-type genes related to stamen and gynoecium development in Liriodendron

The organogenesis and development of reproductive organs, i.e., stamen and gynoecium, are important floral characteristics that are closely related to pollinators and reproductive fitness. As a genus from Magnoliaceae, Liriodendron has only two relict species: L. chinense and L. tulipifera. Despite the similar flower shapes of these species, their natural seed-setting rates differ significantly, implying interspecies difference in floral organogenesis and development. MADS-box genes, which participate in floral organogenesis and development, remain unexplored in Liriodendron. Here, to explore the interspecies difference in floral organogenesis and development and identify MADS-box genes in Liriodendron, we examined the stamen and gynoecium primordia of the two Liriodendron species by scanning electron microscopy combined with paraffin sectioning, and then collected two types of primordia for RNA-seq. A total of 12 libraries were constructed and 42,268 genes were identified, including 35,269 reference genes and 6,999 new genes. Monoterpenoid biosynthesis was enriched in L. tulipifera. Genome-wide analysis of 32 MADS-box genes was conducted, including phylogenetic trees, exon/intron structures, and conserved motif distributions. Twenty-six genes were anchored on 17 scaffolds, and six new genes had no location information. The expression profiles of MIKC-type genes via RT-qPCR acrossing six stamen and gynoecium developmental stages indicates that the PI-like, AG/STK-like, SEP-like, and SVP-like genes may contribute to the species-specific differentiation of the organogenesis and development of reproductive organs in Liriodendron. Our findings laid the groundwork for the future exploration of the mechanism underlying on the interspecific differences in reproductive organ development and fitness in Liriodendron.

The Genetic Control of Nectary Development

Nectar is the most important reward offered by flowering plants to pollinators for pollination services. Since pollinator decline has emerged as a major threat for agriculture, and the food demand is growing globally, studying the nectar gland is of utmost importance. Although the genetic mechanisms that control the development of angiosperm flowers have been quite well understood for many years, the development and maturation of the nectar gland and the secretion of nectar in synchrony with the maturation of the sexual organs appears to be one of the flower's best-kept secrets. Here we review key findings controlling these processes. We also raise key questions that need to be addressed to develop crop ecological functions that take into consideration pollinators' needs.