Alternative Splicing Enhances the Transcriptome Complexity of Liriodendron chinense

Alternative splicing (AS) plays pivotal roles in regulating plant growth and development, flowering, biological rhythms, signal transduction, and stress responses. However, no studies on AS have been performed in Liriodendron chinense, a deciduous tree species that has high economic and ecological value. In this study, we used multiple tools and algorithms to analyze transcriptome data derived from seven tissues via hybrid sequencing. Although only 17.56% (8,503/48,408) of genes in L. chinense were alternatively spliced, these AS genes occurred in 37,844 AS events. Among these events, intron retention was the most frequent AS event, producing 1,656 PTC-containing and 3,310 non-PTC-containing transcripts. Moreover, 183 long noncoding RNAs (lncRNAs) also underwent AS events. Furthermore, weighted gene coexpression network analysis (WGCNA) revealed that there were great differences in the activities of transcription and post-transcriptional regulation between pistils and leaves, and AS had an impact on many physiological and biochemical processes in L. chinense, such as photosynthesis, sphingolipid metabolism, fatty acid biosynthesis and metabolism. Moreover, our analysis showed that the features of genes may affect AS, as AS genes and non-AS genes had differences in the exon/intron length, transcript length, and number of exons/introns. In addition, the structure of AS genes may impact the frequencies and types of AS because AS genes with more exons or introns tended to exhibit more AS events, and shorter introns tended to be retained, whereas shorter exons tended to be skipped. Furthermore, eight AS genes were verified, and the results were consistent with our analysis. Overall, this study reveals that AS and gene interaction are mutual-on one hand, AS can affect gene expression and translation, while on the other hand, the structural characteristics of the gene can also affect AS. This work is the first to comprehensively report on AS in L. chinense, and it can provide a reference for further research on AS in L. chinense.

Positive feedback regulation of a Lycium chinense-derived VDE gene by drought-induced endogenous ABA, and over-expression of this VDE gene improve drought-induced photo-damage in Arabidopsis

Violaxanthin de-epoxidase (VDE) plays an important role in protecting the photosynthetic apparatus from photo-damage by dissipating excessively absorbed light energy as heat, via the conversion of violaxanthin (V) to intermediate product antheraxanthin (A) and final product zeaxanthin (Z) under light stress. We have cloned a VDE gene (LcVDE) from Lycium chinense, a deciduous woody perennial halophyte, which can grow in a large variety of soil types. The amino acid sequence of LcVDE has high homology with VDEs in other plants. Under drought stress, relative expression of LcVDE and the de-epoxidation ratio (Z+0.5A)/(V+A+Z) increased rapidly, and non-photochemical quenching (NPQ) also rose. Interestingly, these elevations induced by drought stress were reduced by the topical administration of abamine SG, a potent ABA inhibitor via inhibition of NCED in the ABA synthesis pathway. Until now, little has been done to explore the relationship between endogenous ABA and the expression of VDE genes. Since V serves as a common precursor for ABA, these data support the possible involvement of endogenous ABA in the positive feedback regulation of LcVDE gene expression in L. chinense under drought stress. Moreover, the LcVDE may be involved in modulating the level of photosynthesis damage caused by drought stress. Furthermore, the ratio of (Z+0.5A)/(V+A+Z) and NPQ increased more in transgenic Arabidopsis over-expressing LcVDE gene than the wild types under drought stress. The maximum quantum yield of primary photochemistry of PSII (Fv/Fm) in transgenic Arabidopsis decreased more slowly during the stressed period than that in wild types under the same conditions. Furthermore, transgenic Arabidopsis over-expressing LcVDE showed increased tolerance to drought stress.

A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense

Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.

Gibberellin Oxidase Gene Family in L. chinense: Genome-Wide Identification and Gene Expression Analysis

GAox is a key enzyme for the transformation of gibberellins, and belongs to the 2-ketoglutarate dependent dioxygenase gene family (2ODD). However, a systematic analysis of GAox in the angiosperm L. chinense has not yet been reported. Here, we identified all LcGAox gene family members in L. chinense, which were classified into the three subgroups of GA20ox, C19GA2ox, and C20GA2ox. Comparison of the gene structure, conserve motifs, phylogenetic relationships, and syntenic relationships of gibberellin oxidase gene families in different species indicated that the gene functional differences may be due to the partial deletion of their domains during evolution. Furthermore, evidence for purifying selection was detected between orthologous GAox genes in rice, grape, Arabidopsis, and L. chinense. Analysis of the codon usage patterns showed that mutation pressure and natural selection might have induced codon usage bias in angiosperms; however, the LcGAox genes in mosses, lycophytes, and ambarella plants exhibited no obvious codon usage preference. These results suggested that the gibberellin oxidase genes were more primitive. The gene expression pattern was analyzed in different organs subjected to multiple abiotic stresses, including GA, abscisic acid (ABA), and chlormequat (CCC) treatment, by RNA-seq and qRT-PCR, and the stress- and phytohormone-responsive cis-elements were counted. The results showed that the synthesis and decomposition of GA were regulated by different LcGAox genes in the vegetative and reproductive organs of L. chinense, and only LcGA2ox1,4, and 7 responded to the NaCl, polyethylene glycol, 4 °C, GA, ABA, and CCC treatment in the roots, stems, and leaves of seedlings at different time periods, revealing the potential role of LcGAox in stress resistance.

Genome-Wide Identification and Expression Analysis of CCO Gene Family in Liriodendron chinense

Carotenoid cleavage oxygenase (CCO) is an enzyme that can catalyze carotenoids to volatile aromatic substances and participate in the biosynthesis of two important phytohormones, i.e., abscisic acid (ABA) and strigolactone (SL). However, the genome-wide identification and analysis of the CCO gene family in the rare and endangered woody plant Liriodendron chinense has not been reported. Here, we performed a genome-wide analysis of the CCO gene family in the L. chinense genome and examined its expression pattern during different developmental processes and in response to various abiotic stresses. A total of 10 LcCCO genes were identified and divided into 6 subfamilies according to the phylogenetic analysis. Subcellular localization prediction showed that most of the LcCCO proteins were located in the cytoplasm. Gene replication analysis showed that segmental and tandem duplication contributed to the expansion of this gene family in the L. chinense genome. Cis-element prediction showed that cis-elements related to plant hormones, stress and light response were widely distributed in the promoter regions of LcCCO genes. Gene expression profile analysis showed that LcNCED3b was extensively involved in somatic embryogenesis, especially the somatic embryo maturation, as well as in response to heat and cold stress in leaves. Furthermore, qRT-PCR analysis showed that LcNCED3b obviously responded to drought stress in roots and leaves. This study provides a comprehensive overview of the LcCCO gene family and a potential gene target for the optimization of the somatic embryogenesis system and resistance breeding in the valuable forest tree L. chinense.

The Identification and Expression Analysis of the Liriodendron chinense F-Box Gene Family

The F-box gene family is one of the largest gene families in plants, and it plays a crucial role in regulating plant development, reproduction, cellular protein degradation, and response to biotic and abiotic stresses. Despite their significance, a comprehensive analysis of the F-box gene family in Liriodendron chinense and other magnoliaceae species has not been reported. In this study, we report for the first time the identification of 144 full-length F-box genes in L. chinense. Based on specific domains and phylogenetic analyses, these genes were divided into 10 distinct subfamilies. We further analyzed their gene structure, conserved domain and chromosome distribution, genome-wide replication events, and collinearity. Additionally, based on GO analysis, we found that F-box genes exhibit functional specificity, with a significant proportion of them being involved in protein binding (GO:0005515), suggesting that F-box genes may play an important role in gene regulation in L. chinense. Transcriptome data and q-PCR results also showed that F-box genes are involved in the development of multiple tissues in L. chinense, regulate the somatic embryogenesis of Liriodendron hybrids, and play a pivotal role in abiotic stress. Altogether, these findings provide a foundation for understanding the biological function of F-box genes in L. chinense and other plant species.