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.

Cinnamomum camphora has been cultivated as an economically important tree for its medicinal and aromatic properties. Selective breeding has produced Cinnamomum plants for special uses, including spice strains with characteristic flavors and aromas and high-potency medicinal cultivars. The molecular biology underlying terpenoid biosynthesis is still unexplored. Gas chromatography-mass spectrometry was used to analyze the differences in contents and compositions of essential oil terpenoids in linalool- and borneol-type chemotypes of C. camphora. The data revealed that the essential oils consist primarily of monoterpenes with only very minor quantities of sesquiterpenes and diterpenes and that the essential oil differs in different chemotypes of C. camphora, with higher yields of (−)-borneol from the borneol-type than from the linalool-type. To study the terpenoid biosynthesis of signature compounds of the major monoterpenes, we performed RNA sequencing to profile the leaf transcriptomes of the two chemotypes of C. camphora. A total of 23.76 Gb clean data was generated from two chemotypes and assembled into 156,184 unigenes. The total length, average length, N50 and GC content of unigenes were 155,645,929 bp, 997 bp, 1430 bp, and 46.5%, respectively. Among them, 76,421 unigenes were annotated by publicly available databases, of which 67 candidate unigenes were identified to be involved in terpenoid biosynthesis in C. camphora. A total of 2863 unigenes were identified to be differentially expression between borneol-type and linalool-type, including 1714 up-regulated and 1149 down-regulated unigenes. Most genes encoding proteins involved in terpenoid precursor MVA and MEP pathways were expressed in similar levels in both chemotypes of C. camphora. In addition, 10 and 17 DEGs were significantly enriched in the terpene synthase activity and oxidoreductase activity terms of their directed acyclic graphs (DAG), respectively. Three monoterpene synthase genes, TPS14-like1, TPS14-like2 and TPS14-like3 were up-regulated in the borneol-type compared to the linalool-type, and their expression levels were further verified using quantitative real-time PCR. This study provides a global overview of gene expression patterns related to terpenoid biosynthesis in C. camphora, and could contribute to a better understanding of the differential accumulation of terpenoids in different C. camphora chemotypes.

As one of the most ancient tree species, the codon usage pattern analysis of Ginkgo biloba is a useful way to understand its evolutionary and genetic mechanisms. Several studies have been conducted on angiosperms, but seldom on gymnosperms. Based on RNA-Seq data of the G. biloba transcriptome, amount to 17,579 unigenes longer than 300 bp were selected and analyzed from 68,547 candidates. The codon usage pattern tended towards more frequently use of A/U-ending codons, which showed an obvious gradient progressing from gymnosperms to dicots to monocots. Meanwhile, analysis of high/low-expression unigenes revealed that high-expression unigenes tended to use G/C-ending codons together with more codon usage bias. Variation of unigenes with different functions suggested that unigenes involving in environment adaptation use G/C-ending codons more frequently with more usage bias, and these results were consistent with the conclusion that the formation of G. biloba codon usage bias was dominated by natural selection.

Cinnamomum camphora, a member of theLauraceaefamily, is a valuable aromatic and timber tree that is indigenous to the south of China and Japan. All parts ofCinnamomum camphorahave secretory cells containing different volatile chemical compounds that are utilized as herbal medicines and essential oils. Here, we reported the complete sequencing of the chloroplast genome ofCinnamomum camphorausing illumina technology. The chloroplast genome ofCinnamomum camphorais 152,570 bp in length and characterized by a relatively conserved quadripartite structure containing a large single copy region of 93,705 bp, a small single copy region of 19,093 bp and two inverted repeat (IR) regions of 19,886 bp. Overall, the genome contained 123 coding regions, of which 15 were repeated in the IR regions. An analysis of chloroplast sequence divergence revealed that the small single copy region was highly variable among the different genera in theLauraceaefamily. A total of 40 repeat structures and 83 simple sequence repeats were detected in both the coding and non-coding regions. A phylogenetic analysis indicated thatCalycanthusis most closely related toLauraceae, both being members ofLaurales, which forms a sister group toMagnoliids. The complete sequence of the chloroplast ofCinnamomum camphorawill aid in in-depth taxonomical studies of theLauraceaefamily in the future. The genetic sequence information will also have valuable applications for chloroplast genetic engineering.

Camellia oil extracted from Camellia seeds is rich in unsaturated fatty acids (UFAs) and secondary metabolites beneficial to human health. However, no oil-tea tree genome has yet been published, which is a major obstacle to investigating the heredity improvement of oil-tea trees. Here, using both Illumina and PicBio sequencing technologies, we present the first chromosome-level genome sequence of the oil-tea tree species Camellia chekiangoleosa Hu. (CCH). The assembled genome consists of 15 pseudochromosomes with a genome size of 2.73 Gb and a scaffold N50 of 185.30 Mb. At least 2.16 Gb of the genome assembly consists of repetitive sequences, and the rest involves a high-confidence set of 64 608 protein-coding gene models. Comparative genomic analysis revealed that the CCH genome underwent a whole-genome duplication (WGD) event shared across the Camellia genus at ~57.48 MYA and a γ-WGT event shared across all core eudicot plants at ~120 MYA. Gene family clustering revealed that the genes involved in terpenoid biosynthesis have undergone rapid expansion. Furthermore, we determined the expression patterns of oleic acid accumulation- and terpenoid biosynthesis-associated genes in six tissues. We found that these genes tend to be highly expressed in leaves, pericarp tissues, roots, and seeds. The first chromosome-level genome of oil-tea trees will provide valuable resources for determining Camellia evolution and utilizing the germplasm of this taxon.