Complete genome sequence of mulberry dwarf phytoplasma strain MDAK-01: comparative genomics and evolutionary analyses
- 1Shaanxi Key Laboratory of Sericulture, Ankang University, Ankang, 725000 China
- 2School of Modern Agriculture & Biotechnology, Ankang University, Ankang, 725000 China
- 3The Sericultural and Silk Research Institute, College of Animal Science and Technology, North West Agriculture and Forestry University, Yangling, 712100 China
Abstract
Background: Phytoplasmas are phloem-limited plant pathogens that cause substantial economic losses in global agriculture. Complete genome sequences are essential for clarifying their taxonomic status, evolutionary relationships, and pathogenic mechanisms. In this study, we generated a high-quality, complete circular genome of Mulberry dwarf phytoplasma (MDP) strain MDAK-01—the first formally described species in the genus Candidatus Phytoplasma—using a hybrid sequencing strategy combining Illumina and Oxford Nanopore Technology. Furthermore, we performed a comprehensive comparative analysis between the MDAK-01 genome and 36 publicly available phytoplasma genomes.
Results: The MDAK-01 genome is 782,702 bp in length, with 755 coding sequences (CDSs), 2 rRNA operons and 32 tRNA genes. Functional annotation based on homology searches identified 753 CDSs. In addition, 10 types of Potential Mobile Units with a total of 61 copies were detected, and 12 secreted proteins were predicted. Phylogenetic analyses based on 151 single-copy homologous proteins and 16 S rRNA gene sequences from 37 phytoplasma strains bothassigned MDAK-01 to the 16SrI group. Average amino acid identity values within 16Sr groups were consistently higher than those between groups. However, inter-group ANI between 16SrI and 16SrXII (97.65%) surpassed the intra-group ANI of the 16SrXII group (96.77%). For the 16SrI group, Relative Synonymous Codon Usage analysis revealed a preference for AT-ending codons. Further analyses of Parity Rule 2-bias, GC12-GC3, and Effective Number of Codons-GC3s indicated that codon usage patterns were mainly shaped by natural selection. All 37 phytoplasma strains possessed the dnaA gene (a potential replication origin), located near the trough of typical or atypical V-shaped cumulative GC-skew plot.
Conclusions: This study presents the complete genome of the MDP strain MDAK-01, analyzes the codon usage bias of the 16SrI group, clarifies the phylogenetic relationships and replication characteristics of 37 phytoplasma strains, and provides valuable genomic resources for the taxonomic revision and evolutionary research of phytoplasmas.
Supplementary Information: The online version contains supplementary material available at https://doi.org/10.1186/s12866-026-05512-8.
Related
- Complete genome sequence of Comamonas testosteroni TA441 and comparative analysis of steroid degradation gene clusters across Proteobacteria and Actinomycetota
- Comparative evolutionary pharmacogenomics of human prostaglandin-endoperoxide synthase paralogs identifies population-structured coding variation and protein-contextual candidates in N-terminal leader-sequence and catalytic-channel regions
- The complete genome of ‘Candidatus Phytoplasma phoenicium’ highlights carboxylic acids as primary carbon and energy substrates in phytoplasmas
- Evolutionary Dynamics of the Complete Chemosensory Repertoire in Kissing Bugs of the Genus Rhodnius : Divergent Odorant Receptors Contrast with Conserved Gene Families
- The complete genome sequence of a poplar, Populus davidiana.
BibTeX
@article{kong2026complete,
title = {{Complete genome sequence of mulberry dwarf phytoplasma strain MDAK-01: comparative genomics and evolutionary analyses}},
author = {Weiqing Kong and Kerong Lv and Fang Yang and Yunwu Peng and Wei Jiang and Jiajun Liang and Jinhong Yang and Lijun Bao},
journal = {BMC Microbiology},
year = {2026},
month = aug,
doi = {10.1186/s12866-026-05512-8},
url = {https://doi.org/10.1186/s12866-026-05512-8},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - JOUR
TI - Complete genome sequence of mulberry dwarf phytoplasma strain MDAK-01: comparative genomics and evolutionary analyses
AU - Weiqing Kong
AU - Kerong Lv
AU - Fang Yang
AU - Yunwu Peng
AU - Wei Jiang
AU - Jiajun Liang
AU - Jinhong Yang
AU - Lijun Bao
PY - 2026
DA - 2026/08/08
JO - BMC Microbiology
DO - 10.1186/s12866-026-05512-8
UR - https://doi.org/10.1186/s12866-026-05512-8
AB - Background: Phytoplasmas are phloem-limited plant pathogens that cause substantial economic losses in global agriculture. Complete genome sequences are essential for clarifying their taxonomic status, evolutionary relationships, and pathogenic mechanisms. In this study, we generated a high-quality, complete circular genome of Mulberry dwarf phytoplasma (MDP) strain MDAK-01—the first formally described species in the genus Candidatus Phytoplasma—using a hybrid sequencing strategy combining Illumina and Oxford Nanopore Technology. Furthermore, we performed a comprehensive comparative analysis between the MDAK-01 genome and 36 publicly available phytoplasma genomes. Results: The MDAK-01 genome is 782,702 bp in length, with 755 coding sequences (CDSs), 2 rRNA operons and 32 tRNA genes. Functional annotation based on homology searches identified 753 CDSs. In addition, 10 types of Potential Mobile Units with a total of 61 copies were detected, and 12 secreted proteins were predicted. Phylogenetic analyses based on 151 single-copy homologous proteins and 16 S rRNA gene sequences from 37 phytoplasma strains bothassigned MDAK-01 to the 16SrI group. Average amino acid identity values within 16Sr groups were consistently higher than those between groups. However, inter-group ANI between 16SrI and 16SrXII (97.65%) surpassed the intra-group ANI of the 16SrXII group (96.77%). For the 16SrI group, Relative Synonymous Codon Usage analysis revealed a preference for AT-ending codons. Further analyses of Parity Rule 2-bias, GC12-GC3, and Effective Number of Codons-GC3s indicated that codon usage patterns were mainly shaped by natural selection. All 37 phytoplasma strains possessed the dnaA gene (a potential replication origin), located near the trough of typical or atypical V-shaped cumulative GC-skew plot. Conclusions: This study presents the complete genome of the MDP strain MDAK-01, analyzes the codon usage bias of the 16SrI group, clarifies the phylogenetic relationships and replication characteristics of 37 phytoplasma strains, and provides valuable genomic resources for the taxonomic revision and evolutionary research of phytoplasmas. Supplementary Information: The online version contains supplementary material available at https://doi.org/10.1186/s12866-026-05512-8.
SN - 1471-2180
ER -