Complete mitochondrial genomes (mitogenomes) of plants are valuable resources for nucleocytoplasmic interactions, plant evolution, and plant cytoplasmic male sterile line breeding. However, the complete assembly of plant mitogenomes is challenging due to frequent recombination events and horizontal gene transfers. Previous studies have adopted Illumina, PacBio, and Nanopore sequencing data to assemble plant mitogenomes, but the poor assembly completeness, low sequencing accuracy, and high cost limit the sampling capacity. Here, we present an efficient assembly toolkit (PMAT) for de novo assembly of plant mitogenomes using low-coverage HiFi sequencing data. PMAT has been applied to the de novo assembly of 13 broadly representative plant mitogenomes, outperforming existing organelle genome assemblers in terms of assembly accuracy and completeness. By evaluating the assembly of plant mitogenomes from different sequencing data, it was confirmed that PMAT only requires 1× HiFi sequencing data to obtain a complete plant mitogenome. The source code for PMAT is available at https://github.com/bichangwei/PMAT. The developed PMAT toolkit will indeed accelerate the understanding of evolutionary variation and breeding application of plant mitogenomes.

MADS-box genes encode transcription factors that affect plant growth and development. Camellia chekiangoleosa is an oil tree species with ornamental value, but there have been few molecular biological studies on the developmental regulation of this species. To explore their possible role in C. chekiangoleosa and lay a foundation for subsequent research, 89 MADS-box genes were identified across the whole genome of C. chekiangoleosa for the first time. These genes were present on all the chromosomes and were found to have expanded by tandem duplication and fragment duplication. Based on the results of a phylogenetic analysis, the 89 MADS-box genes could be divided into either type I (38) or type II (51). Both the number and proportion of the type II genes were significantly greater than those of Camellia sinensis and Arabidopsis thaliana, indicating that C. chekiangoleosa type II genes experienced a higher duplication rate or a lower loss rate. The results of both a sequence alignment and a conserved motif analysis suggest that the type II genes are more conserved, meaning that they may have originated and differentiated earlier than the type I genes did. At the same time, the presence of extra-long amino acid sequences may be an important feature of C. chekiangoleosa. Gene structure analysis revealed the number of introns of MADS-box genes: twenty-one type I genes had no introns, and 13 type I genes contained only 1~2 introns. The type II genes have far more introns and longer introns than the type I genes do. Some MIKC^C genes have super large introns (≥15 kb), which are rare in other species. The super large introns of these MIKC^C genes may indicate richer gene expression. Moreover, the results of a qPCR expression analysis of the roots, flowers, leaves and seeds of C. chekiangoleosa showed that the MADS-box genes were expressed in all those tissues. Overall, compared with that of the type I genes, the expression of the type II genes was significantly higher. The CchMADS31 and CchMADS58 genes (type II) were highly expressed specifically in the flowers, which may in turn regulate the size of the flower meristem and petals. CchMADS55 was expressed specifically in the seeds, which might affect seed development. This study provides additional information for the functional characterization of the MADS-box gene family and lays an important foundation for in-depth study of related genes, such as those involved in the development of the reproductive organs of C. chekiangoleosa.

The gibberellic acid-stimulated Arabidopsis (GASA) gene family plays an important regulatory role in the growth and development of plants. In this study, we identified 19 GASA genes using bioinformatics-based methods in Populus trichocarpa, and these PtGASA genes could be divided into three categories based on their phylogenetic relationships. Based on an analysis of the structure and motifs of these genes, it was concluded that PtGASA class II members are more conserved than class I and class III members are, and the results of collinearity analysis showed that members of class II are collinearly related in poplar. Expression analysis of Populus trichocarpa roots, stems, and leaves showed that most of the PtGASA genes are expressed at higher levels in the stems or roots than in the leaves; a similar expression pattern was found in Vitis vinifera, indicating that the GASA-family members mainly play a role in the morphogenesis of poplar. Considering the phenomenon of gene amplification, we found that the higher the similarity of homologous genes was, the more similar the expression patterns. This study represents the first whole-genome identification and expression-profile analysis of the GASA-gene family in poplar, a model species, laying a foundation for functional studies of poplar GASA genes and serving as a reference for related research on other woody plant species.

Lateral organ boundaries domain (LBD) proteins are plant-specific transcription factors involved in various transcriptional regulation processes. We identified a total of 37 GbLBD genes in ginkgo, and based on gene structure and phylogenetic analysis, the GbLBD gene family was classified into class I (33, with the largest number of Id genes (16)) and class II (4). The ginkgo LBD gene was also analyzed regarding its chromosomal distributions, gene duplications, promoters, and introns/exons. In addition, gene expression profiling and real-time quantitative PCR analysis showed that the expression of 14 GbLBD genes differed in six different tissues and three developmental stages. The GbLBD gene of class II were highly expressed relative to the class I gene in all tissues and developmental stages, while class Id gene were generally at low levels or were not expressed, especially in seed developmental stages. The expression pattern analysis of cold/drought treatment and IAA/ABA hormone treatment showed that abiotic stress treatment could significantly induce the expression of GbLBD gene, of which class II genes played a key role in stress treatment. Our study provides a solid foundation for further evolutionary and functional analysis of the ginkgo LBD gene family.

Color variation in plant leaves has a significant impact on their photosynthesis and plant growth. Camellia chekiangoleosa yellow-leaf mutants are ideal materials for studying the mechanisms of pigment synthesis and photosynthesis, but their mechanism of leaf variation is not clear. We systematically elucidated the intrinsic causes of leaf yellowing in the new Camellia chekiangoleosa variety ‘Diecui Liuji’ in terms of changes in its cell structure, pigment content, and transcript levels. This study indicates that the incomplete structure of chloroplast-like vesicles, the decrease in blue-green chlorophyll a, and the increase in yellow-green chlorophyll b in yellowing leaves are the direct causes of yellowing-leaf formation. The high expression of genes that catalyze the degradation of chlorophyll a (PAO and RCCR) and its conversion to chlorophyll b (CAO) in yellowing leaves leads to a decrease in the chlorophyll a content, while the low expression of CLH genes is the main reason for the increase in the chlorophyll b content. We also found transcription factors such as ERF, E2F, WRKY, MYB, TPC, TGA, and NFYC may regulate their expression. RT-qPCR assays of 12 DEGs confirm the RNA-seq results. This study will provide a foundation for investigating the transcriptional and regulatory mechanisms of leaf color changes.

The WRKY gene family is essential for controlling a variety of plant physiological functions, yet the involvement of specific WRKY members in pigment biosynthesis and accumulation in Camellia chekiangoleosa remains unexplored, particularly in anthocyanins and carotenoids, which play crucial roles in the pigmentation of C. chekiangoleosa. This study systematically identified 87 WRKY genes across 15 chromosomes in C. chekiangoleosa through bioinformatic approaches. Further structural and phylogenetic analyses of these TFs enabled their classification into six different subgroups. WRKY family expansion was shown to be mostly driven by tandem duplication. W-box elements, which can be binding sites for WRKY transcription factors, were present in a number of biosynthetic genes in the pigment production pathway. Yeast one-hybrid assay confirmed that five WRKY transcription factors (CchWRKY15/24/33/47/76) directly bind to the promoter regions of two key biosynthetic genes, CchPSY1 and Cch4CL1. Intriguingly, among the five WRKYs tested, the expression levels of CchWRKY15, CchWRKY33, and CchWRKY47 showed the strongest positive associations with flavonoid accumulation (p < 0.05, Pearson correlation analysis).These findings provide novel insights into the evolutionary patterns, transcriptional regulation, and functional characteristics of CchWRKYs, while elucidating their possible regulatory mechanisms in the fruit coloration of C. chekiangoleosa.