Alfalfa is often called the queen of forage crops, but one of its wilder relatives may hold the keys to its climate-proof future. Medicago sativa subsp. falcata, the yellow-flowered alfalfa native to the steppes of Inner Mongolia and Central Asia, tolerates cold and drought far better than the cultivated alfalfa grown on millions of hectares worldwide. A new study published in BMC Plant Biology has now sequenced and compared the complete chloroplast genomes of diploid and tetraploid forms of this hardy subspecies, placing them within a genus-wide framework of 30 Medicago species. The results offer a detailed portrait of how the plastid genome behaves when a plant’s nuclear genome doubles, and they deliver a trove of molecular markers that breeders can use to build tougher forage cultivars.
The research team, led by Ying Xue and Yunpeng Gai of Beijing Forestry University together with colleagues including Guy Smagghe of the Vrije Universiteit Brussel, focused on a natural experiment that has fascinated evolutionary biologists for decades. Within M. sativa subsp. falcata, both diploid individuals with 16 chromosomes and tetraploid individuals with 32 chromosomes occur in nature. Whole-genome duplication, or polyploidization, is one of the most dramatic events that can happen to a plant lineage, reshaping nuclear gene dosage, silencing duplicate genes, and often triggering sweeping changes in cell size and physiology. Whether and how these nuclear upheavals ripple outward into the chloroplast genome, however, has remained poorly resolved, particularly within a single subspecies where cytotypes can be directly compared.
Chloroplasts carry their own small circular genomes, typically around 120,000 to 160,000 base pairs in flowering plants, descended from free-living cyanobacteria that were engulfed by an ancestral eukaryotic cell more than a billion years ago. Because plastid genomes are maternally inherited in most legumes, relatively compact, and present in thousands of copies per cell, they are workhorses of plant evolutionary biology, used to trace lineages, identify species, and reconstruct geographic histories. The team assembled the complete chloroplast genome of the diploid falcata cytotype at 125,934 base pairs and the tetraploid cytotype at 125,770 base pairs, a difference of just 164 base pairs. Both genomes showed nearly identical guanine-cytosine content of approximately 33.80 percent, and gene counts differed only marginally between the two ploidy levels.
That near-identity is itself a finding. If polyploidization imposed strong selective pressure on the plastid compartment, one might expect shifts in genome size, gene content, or structural organization between cytotypes. Instead, the comparison revealed that the chloroplast genomes of diploid and tetraploid falcata are essentially interchangeable in size, composition, and architecture. The authors conclude that polyploidization exerts limited influence on chloroplast genome size, gene content, and structural organization, and that strong functional constraints on the plastid genome are preserved across ploidy levels. In other words, whatever turbulence whole-genome duplication stirs up in the nucleus, the chloroplast keeps its house in order.
Beyond the ploidy comparison, the study dug into the fine-grained mechanics of chloroplast protein evolution through codon usage bias, the phenomenon by which genes preferentially use certain synonymous codons over others. The genetic code is degenerate, meaning most amino acids are specified by multiple three-letter codons, and organisms are far from random in choosing among them. In the Medicago chloroplast genomes, the analysis revealed a dominant preference for codons ending in adenine or uracil, the A/U-ending codons characteristic of AT-rich plastid genomes. Crucially, the statistical signatures indicated that this bias is shaped primarily by natural selection rather than by mutation pressure or random drift alone, suggesting that translational efficiency and accuracy continue to exert measurable selective forces on chloroplast coding sequences.
The team also examined the evolutionary forces acting on individual protein-coding genes by calculating the ratio of nonsynonymous substitution rates to synonymous substitution rates, abbreviated Ka/Ks. A ratio below one indicates purifying selection, the process by which natural selection removes harmful mutations and preserves protein function; a ratio near or above one would suggest neutral evolution or positive selection. Across the Medicago chloroplast gene set, 77.1 percent of genes showed Ka/Ks values below 0.5, a hallmark of strong purifying selection. The constraint was especially pronounced in genes involved in photosynthesis and energy metabolism, the core functions that the chloroplast exists to perform. This deep conservation underscores how little tolerance the photosynthetic machinery has for change, even across a genus spanning dozens of species and millions of years of divergence.
Not every region of the chloroplast genome is locked down, however, and that variability is precisely what makes the plastome useful for breeding and identification. The researchers mapped simple sequence repeats, short tandemly repeated DNA motifs that mutate rapidly and are abundant in the falcata chloroplast genomes. Most of these SSR loci were mononucleotide repeats, runs of a single base such as poly-A or poly-T tracts, and the vast majority sat in non-coding regions rather than within genes. Because these hypervariable stretches evolve quickly while flanking sequences remain conserved, they are ideal raw material for developing chloroplast-derived molecular markers that can distinguish accessions, assess genetic diversity in germplasm collections, and fingerprint cultivars without the complications of nuclear heterozygosity.
To place falcata in its evolutionary context, the team reconstructed phylogenetic relationships across 30 Medicago species and performed collinearity analysis, which compares the order and orientation of genes along the chloroplast genome. Both approaches demonstrated high structural conservation among closely related species, with gene order rearrangements rare and localized. In the resulting trees, the diploid and tetraploid falcata chloroplast genomes clustered together with other M. sativa accessions, confirming the subspecies-level relationships inferred from morphology and nuclear data and reinforcing the view of falcata as the wild gene pool most closely allied to cultivated alfalfa. The tight clustering also validates chloroplast genomes as reliable barcodes for untangling relationships within a genus famous for taxonomic complexity, driven in part by autopolyploidy and hybridization.
The practical implications extend well beyond evolutionary theory. Yellow-flowered alfalfa is regarded as an essential genetic resource for forage improvement, biofuel production, and sustainable agriculture, and its superior cold tolerance and drought resistance make it suitable for cultivation in marginal and climate-stressed environments where common alfalfa fails. As global agriculture confronts hotter summers, colder snaps, and increasingly unreliable rainfall, the wild cytotypes of falcata represent a reservoir of adaptive alleles that breeders are eager to tap. The chloroplast genomic resources and candidate SSR loci delivered by this study can accelerate germplasm evaluation, phylogenetic inference, and marker-assisted breeding of Medicago forage cultivars, allowing breeders to track valuable cytoplasms through crossing programs and to certify the maternal lineage of elite hybrids.
The study also contributes to a broader scientific conversation about how organellar genomes respond to nuclear genome doubling. Because chloroplasts and mitochondria encode proteins that must interface with thousands of nuclear-encoded partners, polyploidization could in principle disrupt these co-adapted systems and force compensatory evolution in the organelle. The falcata comparison suggests that, at least for the chloroplast, such disruption is minimal: the plastome’s size, GC content, gene inventory, and codon preferences ride out whole-genome duplication essentially unchanged, buffered by intense purifying selection. For a crop genus that feeds livestock across the globe and anchors grassland agriculture from China to North America, that stability is good news. It means the chloroplast markers developed here should remain robust across ploidy levels, and it means the wild, cold-hardy cytoplasms of falcata can be introduced into breeding lines without unpredictable plastid side effects. As sequencing costs continue to fall, comparative plastomics of this kind is poised to become a routine tool in the effort to future-proof the world’s forage crops, one chloroplast genome at a time.
Subject of Research: Comparative chloroplast genome evolution and polyploidy in Medicago sativa subsp. falcata
Article Title: Comparative chloroplast genomics of Medicago: codon usage bias, purifying selection, and phylogenomics
Article References: Xue, Y., Liang, X., Zhao, W., Smagghe, G., Shen, L., Li, S., Wang, Z., Gao, F., Li, D., Zhang, T., & Gai, Y. (2026). Comparative chloroplast genomics of Medicago: codon usage bias, purifying selection, and phylogenomics. BMC Plant Biology. https://doi.org/10.1186/s12870-026-10073-z
Image Credits: AI Generated
DOI: 10.1186/s12870-026-10073-z
Keywords: Medicago, chloroplast genome, polyploidy, codon usage bias, purifying selection, simple sequence repeats, phylogenomics, alfalfa, forage breeding, plant molecular biology, genomics, Comparative
