CAM: an alignment-free method to recover phylogenies using codon aversion motifs

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.

Intraspecific and Intrageneric Genomic Variation across Three Sedum Species (Crassulaceae): A Plastomic Perspective

Sedum is the largest succulent genus in Crassulaceae. Because of predominant maternal inheritance, little recombination, and slow evolution, plastomes can serve as powerful super barcodes for inter- or intra-species phylogenetic analyses. While previous research has focused on plastomes between Sedum species, intra-species studies are scarce. Here, we sequenced plastomes from three Sedum species (Sedum alfredii, Sedum plumbizincicola, and Sedum japonicum) to understand their evolutionary relationships and plastome structural evolution. Our analyses revealed minimal size and GC content variation across species. However, gene distribution at IR boundaries, repeat structures, and codon usage patterns showed diversity at both inter-specific and intra-specific levels. Notably, an rps19 gene expansion and a bias toward A/T-ending codons were observed. Codon aversion motifs also varied, potentially serving as markers for future studies. Phylogenetic analyses confirmed the non-monophyly of Sedum and divided the Acre clade into two groups. Individuals from the same species clustered together, with strong support for the relationships between S. alfredii, S. tricarpum, and S. plumbizincicola. Additionally, S. japonicum clearly affiliates with the Acre clade. This study provides valuable insights into both intra-specific and intra-generic plastome variation in Sedum, as well as overall plastome evolution within the genus.

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.

Plastomes of Bletilla (Orchidaceae) and Phylogenetic Implications

The genus Bletilla is a small genus of only five species distributed across Asia, including B. chartacea, B. foliosa, B. formosana, B. ochracea and B. striata, which is of great medicinal importance. Furthermore, this genus is a member of the key tribe Arethuseae (Orchidaceae), harboring an extremely complicated taxonomic history. Recently, the monophyletic status of Bletilla has been challenged, and the phylogenetic relationships within this genus are still unclear. The plastome, which is rich in both sequence and structural variation, has emerged as a powerful tool for understanding plant evolution. Along with four new plastomes, this work is committed to exploring plastomic markers to elucidate the phylogeny of Bletilla. Our results reveal considerable plastomic differences between B. sinensis and the other three taxa in many aspects. Most importantly, the specific features of the IR junction patterns, novel pttRNA structures and codon aversion motifs can serve as useful molecular markers for Bletilla phylogeny. Moreover, based on maximum likelihood and Bayesian inference methods, our phylogenetic analyses based on two datasets of Arethuseae strongly imply that Bletilla is non-monophyletic. Accordingly, our findings from this study provide novel potential markers for species identification, and shed light on the evolution of Bletilla and Arethuseae.

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

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.

A comprehensive analysis of the phylogenetic signal in ramp sequences in 211 vertebrates

Ramp sequences increase translational speed and accuracy when rare, slowly-translated codons are found at the beginnings of genes. Here, the results of the first analysis of ramp sequences in a phylogenetic construct are presented. Ramp sequences were compared from 247 vertebrates (114 Mammalian and 133 non-mammalian), where the presence and absence of ramp sequences was analyzed as a binary character in a parsimony and maximum likelihood framework. Additionally, ramp sequences were mapped to the Open Tree of Life synthetic tree to determine the number of parallelisms and reversals that occurred, and those results were compared to random permutations. Parsimony and maximum likelihood analyses of the presence and absence of ramp sequences recovered phylogenies that are highly congruent with established phylogenies. Additionally, 81% of vertebrate mammalian ramps and 81.2% of other vertebrate ramps had less parallelisms and reversals than the mean from 1000 randomly permuted trees. A chi-square analysis of completely orthologous ramp sequences resulted in a p-value < 0.001 as compared to random chance. Ramp sequences recover comparable phylogenies as other phylogenomic methods. Although not all ramp sequences appear to have a phylogenetic signal, more ramp sequences track speciation than expected by random chance. Therefore, ramp sequences may be used in conjunction with other phylogenomic approaches if many orthologs are taken into account. However, phylogenomic methods utilizing few orthologs should be cautious in incorporating ramp sequences because individual ramp sequences may provide conflicting signals.

CUBAP: an interactive web portal for analyzing codon usage biases across populations

Synonymous codon usage significantly impacts translational and transcriptional efficiency, gene expression, the secondary structure of both mRNA and proteins, and has been implicated in various diseases. However, population-specific differences in codon usage biases remain largely unexplored. Here, we present a web server, https://cubap.byu.edu, to facilitate analyses of codon usage biases across populations (CUBAP). Using the 1000 Genomes Project, we calculated and visually depict population-specific differences in codon frequencies, codon aversion, identical codon pairing, co-tRNA codon pairing, ramp sequences, and nucleotide composition in 17,634 genes. We found that codon pairing significantly differs between populations in 35.8% of genes, allowing us to successfully predict the place of origin for African and East Asian individuals with 98.8% and 100% accuracy, respectively. We also used CUBAP to identify a significant bias toward decreased CTG pairing in the immunity related GTPase M (IRGM) gene in East Asian and African populations, which may contribute to the decreased association of rs10065172 with Crohn's disease in those populations. CUBAP facilitates in-depth gene-specific and codon-specific visualization that will aid in analyzing candidate genes identified in genome-wide association studies, identifying functional implications of synonymous variants, predicting population-specific impacts of synonymous variants and categorizing genetic biases unique to certain populations.

Codon Pairs are Phylogenetically Conserved: A comprehensive analysis of codon pairing conservation across the Tree of Life

Identical codon pairing and co-tRNA codon pairing increase translational efficiency within genes when two codons that encode the same amino acid are translated by the same tRNA before it diffuses from the ribosome. We examine the phylogenetic signal in both identical and co-tRNA codon pairing across 23 428 species using alignment-free and parsimony methods. We determined that conserved codon pairing typically has a smaller window size than the length of a ribosome, and codon pairing tracks phylogenies across various taxonomic groups. We report a comprehensive analysis of codon pairing, including the extent to which each codon pairs. Our parsimony method generally recovers phylogenies that are more congruent with the established phylogenies than our alignment-free method. However, four of the ten taxonomic groups did not have sufficient orthologous codon pairings and were therefore analyzed using only the alignment-free methods. Since the recovered phylogenies using only codon pairing largely match phylogenies from the Open Tree of Life and the NCBI taxonomy, and are comparable to trees recovered by other algorithms, we propose that codon pairing biases are phylogenetically conserved and should be considered in conjunction with other phylogenomic techniques.

Molecular epidemiology of carbapenem-resistance plasmids using publicly available sequences

Carbapenem-resistant bacteria have quickly become a worldwide concern in nosocomial infections. Of the seven known carbapenemases, four have been shown to be particularly problematic: KPC, NDM, IMP, and VIM. To date, many local and species- or carbapenemase-specific epidemiological studies have been performed, which often focus on the organism itself. This report attempts to perform an inclusive (encompass both species and carbapenemase) epidemiologic study using publicly available plasmid sequences from NCBI. In this report, the gene content of these various plasmids has been characterized, replicon types of the plasmids identified, and the global spread and species promiscuity of the plasmids analyzed. Additionally, support to several groups targeting plasmid maintenance and transfer mechanisms to slow the spread of resistance plasmids is given.