Phylogenomics and plastomics offer new evolutionary perspectives on Kalanchoideae (Crassulaceae).

BACKGROUND AND AIMS: Kalanchoideae is one of three subfamilies within Crassulaceae and contains four genera. Despite previous efforts, the phylogeny of Kalanchoideae remains inadequately resolved with persistent issues including low support, unstructured topologies and polytomies. This study aimed to address two central objectives: (1) resolving the pending phylogenetic questions within Kalanchoideae by using organelle-scale 'barcodes' (plastomes) and nuclear data; and (2) investigating interspecific diversity patterns among Kalanchoideae plastomes. METHODS: To explore the plastome evolution in Kalanchoideae, we newly sequenced 38 plastomes representing all four constituent genera (Adromischus, Cotyledon, Kalanchoe and Tylecodon). We performed comparative analyses of plastomic features, including GC and gene contents, gene distributions at the IR (inverted repeat) boundaries, nucleotide divergence, plastomic tRNA (pttRNA) structures and codon aversions. Additionally, phylogenetic inferences were inferred using both the plastomic dataset (79 genes) and nuclear dataset (1054 genes). KEY RESULTS: Significant heterogeneities were observed in plastome lengths among Kalanchoideae, strongly correlated with LSC (large single copy) lengths. Informative diversities existed in the gene content at SSC/IRa (small single copy/inverted repeat a), with unique patterns individually identified in Adromischus leucophyllus and one major Kalanchoe clade. The ycf1 gene was assessed as a shared hypervariable region among all four genera, containing nine lineage-specific indels. Three pttRNAs exhibited unique structures specific to Kalanchoideae and the genera Adromischus and Kalanchoe. Moreover, 24 coding sequences revealed a total of 41 lineage-specific unused codons across all four constituent genera. The phyloplastomic inferences clearly depicted internal branching patterns in Kalanchoideae. Most notably, by both plastid- and nuclear-based phylogenies, our research offers the first evidence that Kalanchoe section Eukalanchoe is not monophyletic. CONCLUSIONS: This study conducted comprehensive analyses on 38 newly reported Kalanchoideae plastomes. Importantly, our results not only reconstructed well-resolved phylogenies within Kalanchoideae, but also identified highly informative unique markers at the subfamily, genus and species levels. These findings significantly enhance our understanding of the evolutionary history of Kalanchoideae.

Genome evolution of Buchnera aphidicola (Gammaproteobacteria): Insights into strand compositional asymmetry, codon usage bias, and phylogenetic implications.

Taxa of Buchnera aphidicola (hereafter "Buchnera") are mutualistic intracellular symbionts of aphids, known for their remarkable biological traits such as genome reduction, strand compositional asymmetry, and symbiont-host coevolution. With the growing availability of genomic data, we performed a comprehensive analysis of 103 genomes of Buchnera strains from 12 host subfamilies, focusing on the genomic characterizations, codon usage patterns, and phylogenetic implications. Our findings revealed consistent features among all genomes, including small genome sizes, low GC contents, and gene losses. We also identified strong strand compositional asymmetries in all strains at the genome level. Further investigation suggested that mutation pressure may have played a crucial role in shaping codon usage of Buchnera. Moreover, the genomic asymmetries were reflected in asymmetric codon usage preferences within chromosomal genes. Notably, the levels of these asymmetries were varied among strains and were significantly influenced by the degrees of genome shrinkages. Lastly, our phylogenetic analyses presented an alternative topology of Aphididae, based on the Buchnera symbionts, providing robust confirmation of the paraphylies of Eriosomatinae, and Macrosiphini. Our objectives are to further understand the strand compositional asymmetry and codon usage bias of Buchnera taxa, and provide new perspectives for phylogenetic studies of Aphididae.

Comparative genomics of the medicinal plants Lonicera macranthoides and L. japonica provides insight into genus genome evolution and hederagenin-based saponin biosynthesis.

Lonicera macranthoides (LM) and L. japonica (LJ) are medicinal plants widely used in treating viral diseases, such as COVID-19. Although the two species are morphologically similar, their secondary metabolite profiles are significantly different. Here, metabolomics analysis showed that LM contained ~86.01 mg/g hederagenin-based saponins, 2000-fold higher than LJ. To gain molecular insights into its secondary metabolite production, a chromosome-level genome of LM was constructed, comprising 9 pseudo-chromosomes with 40 097 protein-encoding genes. Genome evolution analysis showed that LM and LJ were diverged 1.30-2.27 million years ago (MYA). The two plant species experienced a common whole-genome duplication event that occurred ∼53.9-55.2 MYA before speciation. Genes involved in hederagenin-based saponin biosynthesis were arranged in clusters on the chromosomes of LM and they were more highly expressed in LM than in LJ. Among them, oleanolic acid synthase (OAS) and UDP-glycosyltransferase 73 (UGT73) families were much more highly expressed in LM than in LJ. Specifically, LmOAS1 was identified to effectively catalyse the C-28 oxidation of ß-Amyrin to form oleanolic acid, the precursor of hederagenin-based saponin. LmUGT73P1 was identified to catalyse cauloside A to produce α-hederin. We further identified the key amino acid residues of LmOAS1 and LmUGT73P1 for their enzymatic activities. Additionally, comparing with collinear genes in LJ, LmOAS1 and LmUGT73P1 had an interesting phenomenon of 'neighbourhood replication' in LM genome. Collectively, the genomic resource and candidate genes reported here set the foundation to fully reveal the genome evolution of the Lonicera genus and hederagenin-based saponin biosynthetic pathway.

Structural Diversities and Phylogenetic Signals in Plastomes of the Early-Divergent Angiosperms: A Case Study in Saxifragales.

As representative of the early-divergent groups of angiosperms, Saxifragales is extremely divergent in morphology, comprising 15 families. Within this order, our previous case studies observed significant structural diversities among the plastomes of several lineages, suggesting a possible role in elucidating their deep phylogenetic relationships. Here, we collected 208 available plastomes from 11 constituent families to explore the evolutionary patterns among Saxifragales. With thorough comparisons, the losses of two genes and three introns were found in several groups. Notably, 432 indel events have been observed from the introns of all 17 plastomic intron-containing genes, which could well play an important role in family barcoding. Moreover, numerous heterogeneities and strong intrafamilial phylogenetic implications were revealed in pttRNA (plastomic tRNA) structures, and the unique structural patterns were also determined for five families. Most importantly, based on the well-supported phylogenetic trees, evident phylogenetic signals were detected in combinations with the identified pttRNAs features and intron indels, demonstrating abundant lineage-specific characteristics for Saxifragales. Collectively, the results reported here could not only provide a deeper understanding into the evolutionary patterns of Saxifragales, but also provide a case study for exploring the plastome evolution at a high taxonomic level of angiosperms.

Ten Plastomes of Crassula (Crassulaceae) and Phylogenetic Implications.

The genus Crassula is the second-largest genus in the family Crassulaceae, with about 200 species. As an acknowledged super-barcode, plastomes have been extensively utilized for plant evolutionary studies. Here, we first report 10 new plastomes of Crassula. We further focused on the structural characterizations, codon usage, aversion patterns, and evolutionary rates of plastomes. The IR junction patterns-IRb had 110 bp expansion to rps19-were conservative among Crassula species. Interestingly, we found the codon usage patterns of matK gene in Crassula species are unique among Crassulaceae species with elevated ENC values. Furthermore, subgenus Crassula species have specific GC-biases in the matK gene. In addition, the codon aversion motifs from matK, pafI, and rpl22 contained phylogenetic implications within Crassula. The evolutionary rates analyses indicated all plastid genes of Crassulaceae were under the purifying selection. Among plastid genes, ycf1 and ycf2 were the most rapidly evolving genes, whereas psaC was the most conserved gene. Additionally, our phylogenetic analyses strongly supported that Crassula is sister to all other Crassulaceae species. Our findings will be useful for further evolutionary studies within the Crassula and Crassulaceae.

The Mitogenome of Sedum plumbizincicola (Crassulaceae): Insights into RNA Editing, Lateral Gene Transfer, and Phylogenetic Implications.

As the largest family within the order Saxifragales, Crassulaceae contains about 34 genera with 1400 species. Mitochondria play a critical role in cellular energy production. Since the first land plant mitogenome was reported in Arabidopsis, more than 400 mitogenomic sequences have been deposited in a public database. However, no entire mitogenome data have been available for species of Crassulaceae to date. To better understand the evolutionary history of the organelles of Crassulaceae, we sequenced and performed comprehensive analyses on the mitogenome of Sedum plumbizincicola. The master mitogenomic circle is 212,159 bp in length, including 31 protein-coding genes (PCGs), 14 tRNA genes, and 3 rRNA genes. We further identified totally 508 RNA editing sites in PCGs, and demonstrated that the second codon positions of mitochondrial genes are most prone to RNA editing events. Notably, by neutrality plot analyses, we observed that the mitochondrial RNA editing events have large effects on the driving forces of plant evolution. Additionally, 4 MTPTs and 686 NUMTs were detected in the mitochondrial and nuclear genomes of S. plumbizincicola, respectively. Additionally, we conducted further analyses on gene transfer, secondary structures of mitochondrial RNAs, and phylogenetic implications. Therefore, the findings presented here will be helpful for future investigations on plant mitogenomes.

Plastome evolution of Aeonium and Monanthes (Crassulaceae): insights into the variation of plastomic tRNAs, and the patterns of codon usage and aversion.

MAIN CONCLUSION: This study reported 13 new plastomes from Aeonium and Monanthes, and observed new markers for phylogeny and DNA barcoding, such as novel tRNA structures and codon usage bias and aversion. The Macaronesian clade of Crassulaceae consists of three genera: Aichryson, with about 15 species; Monanthes, with about 10 species; Aeonium, with about 40 species. Within this clade, Aeonium, known as "the botanical equivalent of Darwin's finches", is regarded as an excellent model plant for researching adaptive evolution. Differing from the well-resolved relationships among three genera of the Macaronesian clade, the internal branching patterns within the genus Aeonium are largely unclear. In this study, we first reported 13 new plastomes from genus Aeonium and the closely related genus Monanthes. We further performed comprehensive analyses of the plastomes, with focuses on the secondary structures of pttRNAs and the patterns of codon usage and aversion. With a typical circular and quadripartite structure, the 13 plastomes ranged from 149,900 to 151,030 bp in size, and the unique pattern in IR junctions might become a family-specific marker for Crassulaceae species. Surprisingly, the π values of plastomes from Monanthes were almost twice those from Aeonium. Most importantly, we strongly recommend that highly polymorphic regions, novel putative pttRNA structures, patterns of codon usage bias and aversion derived from plastomes might have phylogenetic implications, and could act as new markers for DNA barcoding of plants. The results of phylogenetic analyses strongly supported a clear internal branching pattern in Macaronesian clade (represented by Aeonium and Monanthes), with higher nodal support values. The findings reported here will provide new insights into the variation of pttRNAs, and the patterns of codon usage and aversion of the family Crassulaceae.

Mitogenomic Codon Usage Patterns of Superfamily Certhioidea (Aves, Passeriformes): Insights into Asymmetrical Bias and Phylogenetic Implications.

The superfamily Certhioidea currently comprises five families. Due to the rapid diversification, the phylogeny of Certhioidea is still controversial. The advent of next generation sequencing provides a unique opportunity for a mitogenome-wide study. Here, we first provided six new complete mitogenomes of Certhioidea (Certhia americana, C. familiaris, Salpornis spilonota, Cantorchilus leucotis, Pheugopedius coraya, and Pheugopedius genibarbis). We further paid attention to the genomic characteristics, codon usages, evolutionary rates, and phylogeny of the Certhioidea mitogenomes. All mitogenomes we analyzed displayed typical ancestral avian gene order with 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, and one control region (CR). Our study indicated the strand-biased compositional asymmetry might shape codon usage preferences in mitochondrial genes. In addition, natural selection might be the main factor in shaping the codon usages of genes. Additionally, evolutionary rate analyses indicated all mitochondrial genes were under purifying selection. Moreover, MT-ATP8 and MT-CO1 were the most rapidly evolving gene and conserved genes, respectively. According to our mitophylogenetic analyses, the monophylies of Troglodytidae and Sittidae were strongly supported. Importantly, we suggest that Salpornis should be separated from Certhiidae and put into Salpornithidae to maintain the monophyly of Certhiidae. Our findings are useful for further evolutionary studies within Certhioidea.

Codon usage patterns and evolution of HSP60 in birds.

Heat shock protein 60 (HSP60) is highly conserved from prokaryotic to eukaryotic organisms, acting as molecular chaperone and other vital biological functions. In spite of increasing knowledge of HSP60, its evolutionary mechanism on functional adaption is still far from completely understood. Moreover, analysis of codon usage bias (CUB) is a powerful tool to understand evolutionary association studies. However, so far, as we know, no scientific work on CUB of HSP60 in birds has been reported. In this study, we provide a comprehensive analysis on the codon usage and molecular evolution of HSP60 across birds. The results indicated that HSP60 had a weak codon usage bias with high ENC values (range from 52.66 to 61), low RSCU, and A/T-ending codons were mostly preferred. Meanwhile, it was considered that mutation, natural selection, and genetic drift combined to shape codon usage patterns with different strength proportions among various birds for HSP60. Then, the LRT tests suggested that different lineages of birds might be under similar selective pressures. Besides, the two positive selection sites (151 and 131) were detected and might undergo radical changes. These findings would contribute to understand function diversity and molecular evolution of HSP60 in birds.

Quantification and evolution of mitochondrial genome rearrangement in Amphibians.

BACKGROUND: Rearrangement is an important topic in the research of amphibian mitochondrial genomes ("mitogenomes" hereafter), whose causes and mechanisms remain enigmatic. Globally examining mitogenome rearrangements and uncovering their characteristics can contribute to a better understanding of mitogenome evolution. RESULTS: Here we systematically investigated mitogenome arrangements of 232 amphibians including four newly sequenced Dicroglossidae mitogenomes. The results showed that our new sequenced mitogenomes all possessed a trnM tandem duplication, which was not exclusive to Dicroglossidae. By merging the same arrangements, the mitogenomes of ~ 80% species belonged to the four major patterns, the major two of which were typical vertebrate arrangement and typical neobatrachian arrangement. Using qMGR for calculating rearrangement frequency (RF) (%), we found that the control region (CR) (RF = 45.04) and trnL2 (RF = 38.79) were the two most frequently rearranged components. Forty-seven point eight percentage of amphibians possessed rearranged mitogenomes including all neobatrachians and their distribution was significantly clustered in the phylogenetic trees (p < 0.001). In addition, we argued that the typical neobatrachian arrangement may have appeared in the Late Jurassic according to possible occurrence time estimation. CONCLUSION: It was the first global census of amphibian mitogenome arrangements from the perspective of quantity statistics, which helped us to systematically understand the type, distribution, frequency and phylogenetic characteristics of these rearrangements.

Characterization of deoxyribonucleic methylation and transcript abundance of sex-related genes during tempera ture-dependent sex determination in Mauremys reevesii†.

A number of genes relevant for sex determination have been found in species with temperature-dependent sex determination. Epigenetics play a key role in sex determination, but characterization of deoxyribonucleic acid methylation of sex-related genes on temperature-dependent sex determination remains unclear. Mauremys reevesii is a typical species with temperature-dependent sex determination. In this study, we analyzed the Cytosine Guanine (CpG) methylation status of the proximal promoters, the messenger ribonucleic acid expression patterns and the correlation between methylation and expression levels of Aromatase, Forkhead box protein L2, Doublesex and mab3-related transcription factor 1, sex-determining region on Y chromosome-box 9, and anti-Müllerian hormone, which are key genes in sex determination in other species. We also analyzed the expression level of genes that encode enzymes involved in methylation and demethylation. The expression levels of Aromatase and Forkhead box protein L2 at the female producing temperature were higher than those at the male producing temperature; the expression levels of Doublesex and mab3-related transcription factor 1, sex-determining region on Y chromosome-box 9, and anti-Müllerian hormone were higher at MPT. The expression of some genes involved in methylation and demethylation is significantly different between male producing temperature and female producing temperature. The expression of messenger ribonucleic acid of genes involved in deoxyribonucleic acid methylation and demethylation affected by temperature, together with other factors, may change the methylation level of the regulatory regions of sex-related genes, which may further lead to temperature-specific expression of sex-related genes, and eventually affect the differentiation of the gonads.

Next-Generation Genome Sequencing of Sedum plumbizincicola Sheds Light on the Structural Evolution of Plastid rRNA Operon and Phylogenetic Implications within Saxifragales.

The genus Sedum, with about 470 recognized species, is classified in the family Crassulaceae of the order Saxifragales. Phylogenetic relationships within the Saxifragales are still unresolved and controversial. In this study, the plastome of S. plumbizincicola was firstly presented, with a focus on the structural analysis of rrn operon and phylogenetic implications within the order Saxifragaceae. The assembled complete plastome of S. plumbizincicola is 149,397 bp in size, with a typical circular, double-stranded, and quadripartite structure of angiosperms. It contains 133 genes, including 85 protein-coding genes (PCGs), 36 tRNA genes, 8 rRNA genes, and four pseudogenes (one ycf1, one rps19, and two ycf15). The predicted secondary structure of S. plumbizincicola 16S rRNA includes three main domains organized in 74 helices. Further, our results confirm that 4.5S rRNA of higher plants is associated with fragmentation of 23S rRNA progenitor. Notably, we also found the sequence of putative rrn5 promoter has some evolutionary implications within the order Saxifragales. Moreover, our phylogenetic analyses suggested that S. plumbizincicola had a closer relationship with S. sarmentosum than S. oryzifolium, and supported the taxonomic revision of Phedimus. Our findings of the present study will be useful for further investigation of the evolution of plastid rRNA operon and phylogenetic relationships within Saxifragales.

Two new mitogenomes of Picidae (Aves, Piciformes): Sequence, structure and phylogenetic analyses.

The family Picidae is classified in the order Piciformes with 236 recognized species in 35 genera. Phylogenetic relationships within the Piciformes are still unresolved by major cause of plumage convergence. This study presented two new mitogenomes of Picidae, Darjeeling woodpecker (Dendrocopos darjellensis) and Gray-capped woodpecker (Yungipicus canicapillus), with the size of 16,727 and 16,847 bp, respectively. They have the identical gene order: 37 genes, including 13 PCGs, two rRNA genes, 22 tRNA genes, and two control regions. Furthermore, an extra cytosine insertion at position 174 of MT-ND3 was also found in both species. The MT-RNR1 of Y. canicapillus shared 89.0% nucleotide sequence identity and a similar secondary structure with that of D. darjellensis, consisting of 46 helices and three main domains. There is more divergence in the MT-RNR2 secondary structures than those of MT-RNR1 between the two species. We found the highest dN/dS for the MT-ATP8 (0.08326) among Piciformes, while the lowest for the MT-CO1 (0.00802). Our phylogenetic analyses suggested that the subfamily Picumninae is paraphyletic, and D. darjellensis had a closer relationship with D. major than D. leucotos. The findings of the present study will be useful for further investigating the evolutionary relationships within Piciformes.

Systematic Identification and Evolution Analysis of Sox Genes in Coturnix japonica Based on Comparative Genomics.

Coturnix japonica (Japanese quail) has been extensively used as a model animal for biological studies. The Sox gene family, which was systematically characterized by a high-mobility group (HMG-box) in many animal species, encodes transcription factors that play central roles during multiple developmental processes. However, genome-wide investigations on the Sox gene family in birds are scarce. In the current study, we first performed a genome-wide study to explore the Sox gene family in galliform birds. Based on available genomic sequences retrieved from the NCBI database, we focused on the global identification of the Sox gene family in C. japonica and other species in Galliformes, and the evolutionary relationships of Sox genes. In our result, a total of 35 Sox genes in seven groups were identified in the C. japonica genome. Our results also revealed that dispersed gene duplications contributed the most to the expansion of the Sox gene family in Galliform birds. Evolutionary analyses indicated that Sox genes are an ancient gene family, and strong purifying selections played key roles in the evolution of CjSox genes of C. japonica. More interestingly, we observed that most Sox genes exhibited highly embryo-specific expression in both gonads. Our findings provided new insights into the molecular function and phylogeny of Sox gene family in birds.

Rapid and recent diversification patterns in Anseriformes birds: Inferred from molecular phylogeny and diversification analyses.

The Anseriformes is a well-known and widely distributed bird order, with more than 150 species in the world. This paper aims to revise the classification, determine the phylogenetic relationships and diversification patterns in Anseriformes by exploring the Cyt b, ND2, COI genes and the complete mitochondrial genomes (mito-genomes). Molecular phylogeny and genetic distance analyses suggest that the Dendrocygna species should be considered as an independent family, Dendrocygnidae, rather than a member of Anatidae. Molecular timescale analyses suggests that the ancestral diversification occurred during the Early Eocene Climatic Optimum (58 ~ 50 Ma). Furthermore, diversification analyses showed that, after a long period of constant diversification, the median initial speciation rate was accelerated three times, and finally increased to approximately 0.3 sp/My. In the present study, both molecular phylogeny and diversification analyses results support that Anseriformes birds underwent rapid and recent diversification in their evolutionary history, especially in modern ducks, which show extreme diversification during the Plio-Pleistocene (~ 5.3 Ma). Therefore, our study support that the Plio-Pleistocene climate fluctuations are likely to have played a significant role in promoting the recent diversification for Anseriformes.

The Mitochondrial Genomes of Aquila fasciata and Buteo lagopus (Aves, Accipitriformes): Sequence, Structure and Phylogenetic Analyses.

The family Accipitridae is one of the largest groups of non-passerine birds, including 68 genera and 243 species globally distributed. In the present study, we determined the complete mitochondrial sequences of two species of accipitrid, namely Aquila fasciata and Buteo lagopus, and conducted a comparative mitogenome analysis across the family. The mitogenome length of A. fasciata and B. lagopus are 18,513 and 18,559 bp with an A + T content of 54.2% and 55.0%, respectively. For both the two accipitrid birds mtDNAs, obvious positive AT-skew and negative GC-skew biases were detected for all 12 PCGs encoded by the H strand, whereas the reverse was found in MT-ND6 encoded by the L strand. One extra nucleotide'C'is present at the position 174 of MT-ND3 gene of A. fasciata, which is not observed at that of B. lagopus. Six conserved sequence boxes in the Domain II, named boxes F, E, D, C, CSBa, and CSBb, respectively, were recognized in the CRs of A. fasciata and B. lagopus. Rates and patterns of mitochondrial gene evolution within Accipitridae were also estimated. The highest dN/dS was detected for the MT-ATP8 gene (0.32493) among Accipitridae, while the lowest for the MT-CO1 gene (0.01415). Mitophylogenetic analysis supported the robust monophyly of Accipitriformes, and Cathartidae was basal to the balance of the order. Moreover, we performed phylogenetic analyses using two other data sets (two mitochondrial loci, and combined nuclear and mitochondrial loci). Our results indicate that the subfamily Aquilinae and all currently polytypic genera of this subfamily are monophyletic. These two novel mtDNA data will be useful in refining the phylogenetic relationships and evolutionary processes of Accipitriformes.