Chloroplast genome sequences have been used to understand evolutionary events and to infer efficiently phylogenetic relationships. Callitropsis funebris (Cupressaceae) is an endemic species in China. Its phylogenetic position is controversial due to morphological characters similar to those of Cupressus, Callitropsis, and Chamaecyparis. This study used next‐generation sequencing technology to sequence the complete chloroplast genome of Ca. funebris and then constructed the phylogenetic relationship between Ca. funebris and its related species based on a variety of data sets and methods. Simple sequence repeats (SSRs) and adaptive evolution analysis were also conducted. Our results showed that the monophyletic branch consisting of Ca. funebris and Cupressus tonkinensis is a sister to Cupressus, while Callitropsis is not monophyletic; Ca. nootkatensis and Ca. vietnamensis are nested in turn at the base of the monophyletic group Hesperocyparis. The statistical results of SSRs supported the closest relationship between Ca. funebris and Cupressus. By performing adaptive evolution analysis under the phylogenetic background of Cupressales, the Branch model detected three genes and the Site model detected 10 genes under positive selection; and the Branch‐Site model uncovered that rpoA has experienced positive selection in the Ca. funebries branch. Molecular analysis from the chloroplast genome highly supported that Ca. funebris is at the base of Cupressus. Of note, SSR features were found to be able to shed some light on phylogenetic relationships. In short, this chloroplast genomic study has provided new insights into the phylogeny of Ca. funebris and revealed multiple chloroplast genes possibly undergoing adaptive evolution.

Simple sequence repeats (SSRs), also known as microsatellites, are short tandem repeat sequences with a motif length of 1-6 bp characterized by high variability and codominant inheritance and have been widely used in species identification, genetic diversity studies, and phylogenetic relationship determination (Chmielewski et al., 2015;Dashnow et al., 2015). SSRs are caused by slipped strand mispairing and subsequent errors during DNA replication, repair, and recombination (Levinson & Gutman, 1987). SSRs are mainly found in intergenic and noncoding regions, with a few present in introns (Li et al., 2004;Liu et al., 2021;Su et al., 2018). Previous studies have shown that the characteristics of genomic SSRs in different taxa (such as their distribution patterns) reflect their phylogenetic relationships (Manee et al., 2020;Srivastava et al., 2019).The distribution of SSRs in some chloroplast (cp) genomes is nonrandom and dominated by mononucleotides, where A/T bases account for the majority (Ellegren, 2004;George et al., 2015;

Cloning the interferon regulatory factor 1 gene in lungfish (Protopterus annectens) and its molecular evolution among sarcopterygians Chinese Science Bulletin 56, 1782 (2011); Molecular evolution of stress-response gene Leptin in high-altitude Chinese snub-nosed monkeys ( Rhinopithecus genus) Chinese Science Bulletin 55, 4132 (2010); Structural characterization and expression pattern analysis of the rice PLT gene family Acta Biochimica et Biophysica Sinica 43, 688 (2011); Characterization and expression pattern of pouII1,a novel class Ⅱ POU gene in zebrafish

No abstract

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein‐coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non‐conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non‐conifers. It is more obvious that in non‐conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6‐ to 3.1‐fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2‐ to 3.6‐fold of non‐conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein‐coding genes due to gene function, plant habitat, or newly acquired IRs.