Research Article |
|
Corresponding author: Lu Tan ( tanlu19910222@163.com ) Academic editor: Brecht Verstraete
© 2025 Lu Tan, Meng Hu, Dan-Dan Wu, Yi-Ran Cheng, Li-Na Sha, Xing Fan, Hou-Yang Kang, Yi Wang, Ana Valdés-Florido, Hai-Qin Zhang, Yong-Hong Zhou.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Tan L, Hu M, Wu D-D, Cheng Y-R, Sha L-N, Fan X, Kang H-Y, Wang Y, Valdés-Florido A, Zhang H-Q, Zhou Y-H (2025) Genome origin and phylogenetic relationships of Campeiostachys (Triticeae: Poaceae) based on nuclear and chloroplast DNA regions. Plant Ecology and Evolution 158(3): 337-349. https://doi.org/10.5091/plecevo.153974
|
Background and aims – Campeiostachys is an allohexaploid perennial genus of the Triticeae tribe (Poaceae). The allopolyploids of Triticeae are produced by interspecific hybridization of different genera. In this study, we investigate the genome origin of Campeiostachys and the relationships of some species based on phylogenetic analyses.
Material and methods – Two nuclear (Acc1 and DMC1) and two chloroplast (matK and rps16) DNA regions of the species of Campeiostachys and its related genera were used for phylogenetic analyses.
Key results – The Acc1 and DMC1 sequences revealed that the genome composition of all Campeiostachys species in our study is StYH, suggesting that Campeiostachys may have originated by the natural hybridization between species with StY and H genomes, as no species with Y or HY genomes have been found in the wild. The results from the chloroplast regions indicated that the maternal donor of the Campeiostachys species contains the St subgenome. In addition, phylogenetic analysis of the nuclear sequences showed that C. purpuraristata always groups with the species of the C. dahurica complex in the St, Y, or H clade, distinct from other species in the genus. Also, C. calcicola, C. kamoji, and C. tsukushiensis var. transiens are distinct yet closely related species.
Conclusion – Campeiostachys species originated from the natural hybridization of the tetraploid species of Roegneria (StY) with the diploid species of Hordeum (H), with Roegneria (StY) acting as the maternal donor. Campeiostachys purpuraristata should be classified into the C. dahurica complex and treated as C. dahurica var. purpuraristata.
allohexaploid, Campeiostachys, genome constitution, maternal donor, natural hybridization, phylogeny
Hybridization and polyploidization play a key role in plant evolution and speciation (
The genus Campeiostachys Drobow was established by
The relationship between some species within Campeiostachys remains unclear. For example, the taxa included in the Campeiostachys dahurica (Turcz. ex Griseb.) B.R.Baum, J.L.Yang & C.Yen (= Elymus dahuricus Turcz. ex Griseb.) complex have been debated. Because the morphological difference is small,
Phylogenetic analyses have been proven as a fast and effective way for identifying genome composition, species relationships, and progenitor species of allopolyploid taxa, revealing the origin and evolutionary history of polyploid plants (
In this context, the objectives of this study on 15 Campeiostachys polyploids and the diploid and polyploids of related genera are: (1) to elucidate the genome origin of Campeiostachys; (2) to investigate the maternal donor of Campeiostachys species; (3) to explore the phylogenetic relationships among Campeiostachys species.
Most of the material was collected by the authors’ research team, except for material of Campeiostachys drobovii (Nevski) B.R.Baum, J.L.Yang, & C.Yen (PI 314203), C. tsukushiensis var. transiens (PI 276396), and Bromus inermis Leyss. (PI 618974), which was kindly provided by the USDA National Plant Germplasm System (https://www.ars-grin.gov). Specimens of the material listed in Table
| Species | Genome | Accession No. | Origin |
| Campeiostachys aristiglumis | StYH | Y 0614 | Xinjiang |
| Campeiostachys calcicola | StYH | ZY 1005 | Sichuan |
| Campeiostachys drobovii | StYH | PI 314203 | Russian |
| Campeiostachys dahurica | StYH | ZY 11033 | Inner Mongolia |
| Campeiostachys dahurica var. cylindrica | StYH | Y 0750 | Xinjiang |
| Campeiostachys dahurica var. excelsis | StYH | ZY 11034 | Inner Mongolia |
| Campeiostachys dahurica var. tangutorum | StYH | Y 2092 | Sichuan |
| Campeiostachys kamoji | StYH | ZY 1007 | Sichuan |
| Campeiostachys nutans | StYH | Y 2235 | - |
| Campeiostachys purpuraristata | StYH | ZY 11075 | Inner Mongolia |
| Campeiostachys schrenkiana | StYH | Y 2426 | - |
| Campeiostachys tsukushiensis var. transiens | StYH | PI 276396 | Sweden |
| Elymus atratus | StYH | ZY 15005 | Sichuan |
| Elymus breviaristatus | StYH | ZY 17008 | Sichuan |
| Elymus sinosubmuticus | StYH | ZY 17004 | Sichuan |
| Hordeum bogdanii | H | ZY 11066 | Inner Mongolia |
| Bromus inermis subsp. inermis | - | PI 618974 | Xinjiang |
In addition to the material mentioned in Table
The total genomic DNA was extracted from fresh leaves using the CTAB method (
| Gene | Primer | Sequence of primer (5’-3’) | PCR profiles |
| Acc1 | F | CCCAATATTTATCATGAGACTTGCA | 1 cycle: 5 min 95°C; 35 cycles: 30 s 95°C, 30 s 56°C, 2 min 30 s 68°C; 1 cycle: 10 min 68°C |
| R | CAACATTTGAATGAAThCTCCACG | ||
| DMC1 | F | TGCCAATTGCTGAGAGATTTG | 1 cycle: 4 min 95°C; 35 cycles: 1 min 95°C, 1 min 52°C, 1 min 72°C; 1 cycle: 10 min 72°C |
| R | AGCCACCTGTTGTAATCTGG | ||
| matK | F | CGATCTATTCATTCAATATTTC | 1 cycle: 4 min 95°C; 35 cycles: 1 min 95°C, 1 min 50°C, 1 min 30 s 72°C; 1 cycle: 10 min 72°C |
| R | TCTAGCACACGAAAGTCGAAGT | ||
| rps16 | F | AAACGATGTGGTAGAAAGCAAC | 1 cycle: 4 min 95°C; 35 cycles: 1 min 95°C, 1 min 53°C, 1 min 72°C; 1 cycle: 10 min 72°C |
| R | AAACGATGTGGTAGAAAGCAAC |
DNA sequences were confirmed through BLAST (
The length of the Acc1 sequences of the Campeiostachys species ranges from 1423 to 1448 bp. The data matrix contains 1827 characters, of which 288 are parsimony uninformative and 134 are parsimony informative. The Acc1 data matrix of 88 sequences was analysed with ML using the TIM1+I+G model (-Ln likelihood = 8078.0741). The assumed nucleotide frequencies were A = 0.2546, C = 0.1827, G = 0.2161, T = 0.3467. The tree topology generated by the BI analysis is similar to that inferred by the ML analysis. The ML tree with bootstrap support values (BS, above the branches) and Bayesian posterior probability (PP, below the branches) is displayed in Fig.
Maximum likelihood tree derived from Acc1 sequences of Campeiostachys and related species. The capital letters in brackets after the species name indicate the genome composition of the species. The numbers above and below the branches indicate bootstrap values > 50% and Bayesian posterior probability values > 0.90, respectively.
All Campeiostachys species have three copies of the Acc1 sequence, which are grouped in the St, Y, and H clades (Fig.
The length of the DMC1 sequences of the Campeiostachys species ranges from 1013 to 1087 bp. The data matrix contains 1266 characters, of which 249 are parsimony uninformative and 103 are parsimony informative. The DMC1 data matrix was analysed with ML using the TIM3+G model (-Ln likelihood = 5162.2108). The assumed nucleotide frequencies were A = 0.3219, C = 0.2140, G = 0.2088, T = 0.2553. The phylogenetic analysis of 103 DMC1 sequences was performed using Bromus inermis as the outgroup (Fig.
Maximum likelihood tree derived from DMC1 sequences of Campeiostachys and related species. The capital letters in brackets after the species name indicate the genome composition of the species. The numbers above and below the branches indicate bootstrap values > 50% and Bayesian posterior probability values > 0.90, respectively.
Three DMC1 sequence copies of the Campeiostachys species are divided into three well-supported clades, which are named the St, Y, and H clades (Fig.
The phylogenetic tree constructed by combining Acc1 and DMC1 sequences is consistent with the one constructed by the single regions. All Campeiostachys species are divided into three clades (Fig.
Maximum likelihood tree derived from Acc1+DMC1 sequences of Campeiostachys and related species. The capital letters in brackets after the species name indicate the genome composition of the species. The numbers above and below the branches indicate bootstrap values > 50% and Bayesian posterior probability values > 0.90, respectively.
In the Y clade, the Campeiostachys species are clustered together (BS = 93%, PP = 1.00). Of which, C. dahurica, C. dahurica var. cylindrica, C. dahurica var. excelsis, C. dahurica var. tangutorum, C. purpuraristata cluster together (BS = 80%, PP = 0.99). Besides, C. kamoji, C. calcicola, and C. tsukushiensis var. transiens cluster into one group (BS = 99%, PP = 1.00). The H clade (BS = 100%, PP = 1.00) includes the species of Hordeum (H), Elymus (StH), and Campeiostachys (StYH). Among them, C. kamoji, C. calcicola, and C. tsukushiensis var. transiens are grouped together (BS = 100%, PP = 1.00). Besides, C. dahurica, C. dahurica var. cylindrica, C. dahurica var. excelsis, and C. purpuraristata cluster together (BS = 94%, PP = 1.00).
The matK matrix contains 60 taxa and 844 characters, including 99 variable information loci and 38 parsimony informative loci. The phylogenetic analysis was based on maximum likelihood (ML) using GTR+I+G as the best-fit model (-Ln likelihood = 2142.3028). The assumed nucleotide frequencies were A = 0.3094, C = 0.1808, G = 0.1523, T = 0.3575. Both ML and BI trees show the matK sequences of Campeiostachys species divided into the St+V+E clade (BS = 51%) (Fig.
Maximum likelihood tree derived from chloroplast regions of Campeiostachys and related species. A. matK. B. rps16. C. matK+rps16. The capital letters in brackets after the species name indicate the genome composition of the species. The numbers above and below the branches indicate bootstrap values > 50% and Bayesian posterior probability values > 0.90, respectively.
A total of 53 rps16 sequences were used for ML analysis. The rps16 sequences matrix contains 706 characters, of which 47 informative loci and 22 parsimony informative loci. The phylogenetic analysis based on the rps16 sequences was conducted using TIM1+G, which was identified as the best-fit model (-Ln likelihood = 1688.6674). The assumed nucleotide frequencies were A = 0.2991, C = 0.1925, G = 0.1478, T = 0.3606. In addition to Campeiostachys species, the St+V+E clade (BS = 50%) also included the diploid species of Pseudoroegneria (St), Lophopyrum (Ee), Thinopyrum (Eb), and Dasypyrum (V) (Fig.
The BI tree and ML tree based on concatenated gene sequences exhibit highly similar topologies. All Campeiostachys species cluster in the same clade, named the St+V+E clade (BS = 58%, PP = 0.91) (Fig.
Traditionally, the species in Triticeae with the same genome or genome combination have been classified into the same genus (
Allopolyploids in Triticeae arise by interspecific hybridization involving different genera with different genome compositions (
The cpDNA is maternally inherited in grasses (
Based on the analysis of the single-copy nuclear gene and cpDNA sequences (Figs
Elymus purpuraristatus (= Campeiostachys purpuraristata) is a perennial grass, mainly distributed in Inner Mongolia (
Elymus purpuraristatus
C.P.Wang & H.L.Yang, Bulletin of Botanical Research Harbin 4(4): 83. 1984. (
Campeiostachys purpuraristata
(C.P.Wang & X.L.Yang) Y.H.Zhou, H.Q.Zhang & Wei Huan Chen (
CHINA • Inner Mongolia, Daqing Mountains; 6 Aug. 1965; C.P. Wang 278; holotype: NMAC!.
Roegneria calcicola Keng (= Campeiostachys calcicola) and Roegneria kamoji (Ohwi) Ohwi ex Keng (= C. kamoji) are perennial herbs within Triticeae (
Campeiostachys species originated from natural hybridization between the tetraploid species of Roegneria (StY) and the diploid species of Hordeum (H), with Roegneria (StY) acting as the maternal donor. Campeiostachys purpuraristata should be classified into the Campeiostachys dahurica complex and treated as Campeiostachys dahurica var. purpuraristata (C.P.Wang & H.L.Yang) Y.H.Zhou, H.Q.Zhang, W.H.Chen & L.Tan.
We would like to express our appreciation to Xiaoxia Zhu, Yang Song, and Qingxiang Huang for their management of the experimental plants. This project was supported by the National Natural Science Foundation of China (No. 32200180), the Science and Technology Bureau of Sichuan Province (2023NSFSC1995).
The Acc1 sequences used in the phylogenetic analyses.
The DMC1 sequences used in the phylogenetic analyses.
The matK and rps16 sequences used in the phylogenetic analyses.