Full-length transcriptome characterization and comparative analysis of Gleditsia sinensis

Analysis of differential mRNA and miRNA expression induced by heterogeneous grafting in Gleditsia sinensis

Gleditsia sinensis Lam. is a multifaceted plant with medicinal, edible, chemical, timber, and ornamental applications. However, the effect of rootstocks on scions after grafting is still unclear. This study examined the mRNA and miRNA transcriptome among homografts, heterografts, and seedlings. GO enrichment analysis between seedlings and homograft/heterograft combinations revealed that biosynthesis, degradation, and transport were enriched. The KEGG enrichment results showed that plant hormone signal transduction and the plant MAPK signaling pathway were enriched in both seedlings and heterograft combinations. Through weighted correlation network analysis (WGCNA), the hub genes related to the content of plant hormones were obtained. Taking G. sinensis as the scion, there were 4594, 2887, 3429, and 5959 mRNAs that were specifically expressed in the grafted plants of G. sinensis/G. fera, G. sinensis/G. delavayi, G. sinensis/G. microphylla, and G. sinensis/G. japonica, respectively. The specifically expressed mRNA genes may participate in such processes and pathways as the rhythmic process, circadian rhythm, gibberellic-acid-mediated signaling pathway, and peptide-based amino acid modification. Additionally, 3, 16, 2, and 15 specifically expressed miRNAs were identified. This study examines the impact of grafting on gene expression in Gleditsia plants and establishes a foundation for the development of new resources and rootstock breeding.

Characterization of the chloroplast genome of Gleditsia species and comparative analysis

The genus Gleditsia has significant medicinal and economic value, but information about the chloroplast genomic characteristics of Gleditsia species has been limited. Using the Illumina sequencing, we assembled and annotated the whole chloroplast genomes of seven Gleditsia species (Gleditsia sinensis, Gleditsia japonica var. delavayi (G. delavayi), G. fera, G. japonica, G. microphylla, Fructus Gleditsiae Abnormalis (Zhū Yá Zào), G. microphylla mutant). The assembled genomes revealed that Gleditsia species have a typical circular tetrad structure, with genome sizes ranging from 162,746 to 170,907 bp. Comparative genomic analysis showed that most (65.8-75.8%) of the abundant simple sequence repeats in Gleditsia and Gymnocladus species were located in the large single copy region. The Gleditsia chloroplast genome prefer T/A-ending codons and avoid C/G-ending codons, positive selection was acting on the rpoA, rpl20, atpB, ndhA and ycf4 genes, most of the chloroplast genes of Gleditsia species underwent purifying selection. Expansion and contraction of the inverted repeat (IR)/single copy (SC) region showed similar patterns within the Gleditsia genus. Polymorphism analysis revealed that coding regions were more conserved than non-coding regions, and the IR region was more conserved than the SC region. Mutational hotspots were mostly found in intergenic regions such as "rps16-trnQ", "trnT-trnL", "ndhG-ndhI", and "rpl32-trnL" in Gleditsia. Phylogenetic analysis showed that G. fera is most closely related to G. sinensis,G. japonica and G. delavayi are relatively closely related. Zhū Yá Zào can be considered a bud mutation of the G. sinensis. The albino phenotype of G. microphylla mutant is not caused by variations in the chloroplast genome, and that the occurrence of the albino phenotype may be due to mutations in chloroplast-related genes involved in splicing or localization functions. This study will help us enhance our exploration of the genetic evolution and geographical origins of the Gleditsia genus.