The Complete Chloroplast Genome Sequence of the Medicinal Plant Salvia miltiorrhiza

Salvia miltiorrhiza is an important medicinal plant with great economic and medicinal value. The complete chloroplast (cp) genome sequence of Salvia miltiorrhiza, the first sequenced member of the Lamiaceae family, is reported here. The genome is 151,328 bp in length and exhibits a typical quadripartite structure of the large (LSC, 82,695 bp) and small (SSC, 17,555 bp) single-copy regions, separated by a pair of inverted repeats (IRs, 25,539 bp). It contains 114 unique genes, including 80 protein-coding genes, 30 tRNAs and four rRNAs. The genome structure, gene order, GC content and codon usage are similar to the typical angiosperm cp genomes. Four forward, three inverted and seven tandem repeats were detected in the Salvia miltiorrhiza cp genome. Simple sequence repeat (SSR) analysis among the 30 asterid cp genomes revealed that most SSRs are AT-rich, which contribute to the overall AT richness of these cp genomes. Additionally, fewer SSRs are distributed in the protein-coding sequences compared to the non-coding regions, indicating an uneven distribution of SSRs within the cp genomes. Entire cp genome comparison of Salvia miltiorrhiza and three other Lamiales cp genomes showed a high degree of sequence similarity and a relatively high divergence of intergenic spacers. Sequence divergence analysis discovered the ten most divergent and ten most conserved genes as well as their length variation, which will be helpful for phylogenetic studies in asterids. Our analysis also supports that both regional and functional constraints affect gene sequence evolution. Further, phylogenetic analysis demonstrated a sister relationship between Salvia miltiorrhiza and Sesamum indicum. The complete cp genome sequence of Salvia miltiorrhiza reported in this paper will facilitate population, phylogenetic and cp genetic engineering studies of this medicinal plant.     

Complete Chloroplast Genome Sequence of <Emphasis Type="Italic">Glycine max</Emphasis> and Comparative Analyses with other Legume Genomes

Lack of complete chloroplast genome sequences is still one of the major limitations to extending chloroplast genetic engineering technology to useful crops. Therefore, we sequenced the soybean chloroplast genome and compared it to the other completely sequenced legumes, Lotus and Medicago. The chloroplast genome of Glycine is 152,218 basepairs (bp) in length, including a pair of inverted repeats of 25,574 bp of identical sequence separated by a small single copy region of 17,895 bp and a large single copy region of 83,175 bp. The genome contains 111 unique genes, and 19 of these are duplicated in the inverted repeat (IR). Comparisons of Glycine, Lotus and Medicago confirm the organization of legume chloroplast genomes based on previous studies. Gene content of the three legumes is nearly identical. The rpl22 gene is missing from all three legumes, and Medicago is missing rps16 and one copy of the IR. Gene order in Glycine, Lotus, and Medicago differs from the usual gene order for angiosperm chloroplast genomes by the presence of a single, large inversion of 51 kilobases (kb). Detailed analyses of repeated sequences indicate that many of the Glycine repeats that are located in the intergenic spacer regions and introns occur in the same location in the other legumes and in Arabidopsis, suggesting that they may play some functional role. The presence of small repeats of psbA and rbcL in legumes that have lost one copy of the IR indicate that this loss has only occurred once during the evolutionary history of legumes.     

The complete chloroplast genome sequence and phylogenetic analysis of <Emphasis Type="Italic">Chuanminshen</Emphasis> (<Emphasis Type="Italic">Chuanminshenviolaceum</Emphasis> Sheh et Shan)

Chloroplast genome sequences are very useful for species identification and phylogenetics. Chuanminshen (Chuanminshen violaceum Sheh et Shan) is an important traditional Chinese medicinal plant, for which the phylogenetic position is still controversial. In this study, the complete chloroplast genome of Chuanminshen violaceum Sheh et Shan was determined. The total size of Chuanminshen chloroplast genome was 154,529 bp with 37.8% GC content. It has the typical quadripartite structure, a large single copy (17,800 bp) and a small single copy (84,171 bp) and a pair of inverted repeats (26,279 bp). The whole genome harbors 132 genes, which includes 85 protein coding genes, 37 tRNA genes, eight rRNA genes, and two pseudogenes. Thirty-nine SSR loci, 32 tandem repeats and 49 dispersed repeats were found. Phylogenetic analyses results with the help of MEGA showed a new insight for the Chuanminshen phylogenetic relationship with the reported chloroplast genomes in Apiales plants.     

Complete chloroplast genome of Myracrodruon urundeuva and its phylogenetics relationships in Anacardiaceae family

Continuous exploratory use of tree species is threatening the existence of several plants in South America. One of these threatened species is Myracroduron urundeuva, highly exploited due to the high quality and durability of its wood. The chloroplast (cp) has been used for several evolutionary studies as well traceability of timber origin, based on its gene sequences and simple sequence repeats (SSR) variability. Cp genome organization is usually consisting of a large single copy and a small single copy region separated by two inverted repeats regions. We sequenced the complete cp genome from M. urundeuva based on Illumina next-generation sequencing. Our results show that the cp genome is 159,883 bp in size. The 36 SSR identified ranging from mono- to hexanucleotides. Positive selection analysis revealed nine genes related to photosystem, protein synthesis, and DNA replication, and protease are under positive selection. Genome comparison a other Anacardiaceae chloroplast genomes showed great variability in the family. The phylogenetic analysis using complete chloroplast genome sequences of other Anacardiaceae family members showed a close relationship with two other economically important genera, Pistacia and Rhus. These results will help future investigations of timber monitoring and population and evolutionary studies. Supplementary InformationThe online version contains supplementary material available at 10.1007/s12298-021-00989-1.     

Chloroplast genome sequences of Artemisia maritima and Artemisia absinthium: Comparative analyses,mutational hotspots in genus Artemisia and phylogeny in family Asteraceae

Artemisia L. is a complex genus of medicinal importance. Publicly available chloroplast genomes of few Artemisia species are insufficient to resolve taxonomic discrepancies at species level. We report chloroplast genome sequences of two further Artemisia species: A. maritima (151,061 bp) and A. absinthium (151,193 bp). Both genomes possess typical quadripartite structure comprising of a large single copy, a small single copy and a pair of long inverted repeats. The two genomes exhibited high similarities in genome sizes, gene synteny, GC content, synonymous and non-synonymous substitutions, codon usage, amino acids frequencies, RNA editing sites, microsatellites, and oligonucleotide repeats. Transition to transversion ratio was <1. Maximum likelihood tree showed Artemisia a monophyletic genus, sister to genus Chrysanthemum. We also identified 20 highly polymorphic regions including rpoC2-rps2, trnR-UCU-trnG-UCC, rps18-rpl20, and trnL-UAG-rpl32 that could be used to develop authentic and cost-effective markers to resolve taxonomic discrepancies and infer phylogenetic relationships among Artemisia species.     

Chloroplast genome of Hibiscus rosa-sinensis (Malvaceae): Comparative analyses and identification of mutational hotspots

Previous studies to resolve phylogenetic and taxonomic discrepancies of Hibiscus remained inconclusive. Here, we report chloroplast genome sequence of Hibiscus rosa-sinensis. Hibiscus rosa-sinensis chloroplast genome was 160,951 bp, comprising of large single copy (89,509 bp) and small single copy (20,246 bp) regions, separated by IRa and IRb (25,598 bp each). The genome contained 130 genes including 85 protein-coding genes, 37 transfer RNAs and 8 ribosomal RNAs. Comparative analyses of chloroplast genomes revealed similar structure among 12 species within family Malvaceae. Evolutionary rates of 77 protein-coding genes showed 95% similarities. Analyses of codon usage, amino acid frequency, putative RNA editing sites, and repeats showed a great extent of similarities between Hibiscus rosa-sinensis and Hibiscus syriacus. We identified 30 mutational hotpots including psbZ-trnG, trnK-rps16, trnD-trnY, trnW-trnP, rpl33-rps18, petG-trnW, trnS-trnG, trnH-psbA, atpB-rbcL, and rpl32-trnL that might be used as polymorphic and robust markers to resolve phylogenetic discrepancies in genus Hibiscus.     

The complete chloroplast genome sequence of the endangered species Syringa pinnatifolia (Oleaceae)

Syringa pinnatifolia is an endangered endemic species in China with important ornamental and medicinal value, and it needs urgent protection. Here, we report the complete chloroplast (cp) genome structure of S. pinnatifolia and its evolution is inferred through comparative studies with related species. The S. pinnatifolia cp genome was 155 326 bp and contained a large single copy region (LSC) of 86 167 bp and a small single copy region (SSC) of 17 775 bp, as well as a pair of inverted repeat regions (IRs) of 25 692 bp. A total of 113 unique genes were annotated, including 79 protein‐coding genes, 30 tRNA genes and four rRNA genes. The GC content of the S. pinnatifolia cp genome was 37.9%, and the corresponding values in the LSC, SSC and IR regions were 36.0, 32.1, 43.2% respectively. Repetitive sequences analysis revealed that the S. pinnatifolia cp genome contained 38 repeats. Microsatellite marker detection analysis identified 253 simple sequence repeats (SSRs), which provides opportunities for future studies of the population genetics and phylogenetic relationships of Syringa. Phylogenetic analysis of 29 selected cp genomes revealed that S. pinnatifolia is closely related to Syringa vulgaris and all 27 Lamiales species formed a clade separate from the two outgroup species. This newly characterized S. pinnatifolia chloroplast genome will provide a useful genomic resource of phylogenetic inference and the development of more genetic markers for species discrimination and population studies in the genus Syringa.     

The First Complete Chloroplast Genome Sequences in Actinidiaceae: Genome Structure and Comparative Analysis

Actinidia chinensis is an important economic plant belonging to the basal lineage of the asterids. Availability of a complete Actinidia chloroplast genome sequence is crucial to understanding phylogenetic relationships among major lineages of angiosperms and facilitates kiwifruit genetic improvement. We report here the complete nucleotide sequences of the chloroplast genomes for Actinidia chinensis and A. chinensis var deliciosa obtained through de novo assembly of Illumina paired-end reads produced by total DNA sequencing. The total genome size ranges from 155,446 to 157,557 bp, with an inverted repeat (IR) of 24,013 to 24,391 bp, a large single copy region (LSC) of 87,984 to 88,337 bp and a small single copy region (SSC) of 20,332 to 20,336 bp. The genome encodes 113 different genes, including 79 unique protein-coding genes, 30 tRNA genes and 4 ribosomal RNA genes, with 16 duplicated in the inverted repeats, and a tRNA gene (trnfM-CAU) duplicated once in the LSC region. Comparisons of IR boundaries among four asterid species showed that IR/LSC borders were extended into the 5’ portion of the psbA gene and IR contraction occurred in Actinidia. The clap gene has been lost from the chloroplast genome in Actinidia, and may have been transferred to the nucleus during chloroplast evolution. Twenty-seven polymorphic simple sequence repeat (SSR) loci were identified in the Actinidia chloroplast genome. Maximum parsimony analyses of a 72-gene, 16 taxa angiosperm dataset strongly support the placement of Actinidiaceae in Ericales within the basal asterids.     

Complete chloroplast genome sequences from Korean ginseng (Panax schinseng Nees) and comparative analysis of sequence evolution among 17 vascular plants.

The nucleotide sequence of Korean ginseng (Panax schinseng Nees) chloroplast genome has been completed (AY582139). The circular double-stranded DNA, which consists of 156,318 bp, contains a pair of inverted repeat regions (IRa and IRb) with 26,071 bp each, which are separated by small and large single copy regions of 86,106 bp and 18,070 bp, respectively. The inverted repeat region is further extended into a large single copy region which includes the 5' parts of the rpsl9 gene. Four short inversions associated with short palindromic sequences that form stem-loop structures were also observed in the chloroplast genome of P. schinseng compared to that of Nicotiana tabacum. The genome content and the relative positions of 114 genes (75 peptide-encoding genes, 30 tRNA genes, 4 rRNA genes, and 5 conserved open reading frames [ycfs]), however, are identical with the chloroplast DNA of N. tabacum. Sixteen genes contain one intron while two genes have two introns. Of these introns, only one (trnL-UAA) belongs to the self-splicing group I; all remaining introns have the characteristics of six domains belonging to group II. Eighteen simple sequence repeats have been identified from the chloroplast genome of Korean ginseng. Several of these SSR loci show infra-specific variations. A detailed comparison of 17 known completed chloroplast genomes from the vascular plants allowed the identification of evolutionary modes of coding segments and intron sequences, as well as the evaluation of the phylogenetic utilities of chloroplast genes. Furthermore, through the detailed comparisons of several chloroplast genomes, evolutionary hotspots predominated by the inversion end points, indel mutation events, and high frequencies of base substitutions were identified. Large-sized indels were often associated with direct repeats at the end of the sequences facilitating intra-molecular recombination.     

The complete chloroplast genome sequences of Solanum tuberosum and comparative analysis with Solanaceae species identified the presence of a 241-bp deletion in cultivated potato chloroplast DNA sequence

The complete nucleotide sequence of the chloroplast genome of potato Solanum tuberosum L. cv. Desiree was determined. The circular double-stranded DNA, which consists of 155,312 bp, contains a pair of inverted repeat regions (IRa, IRb) of 25,595 bp each. The inverted repeat regions are separated by small and large single copy regions of 18,373 and 85,749 bp, respectively. The genome contains 79 proteins, 30 tRNAs, 4 rRNAs, and unidentified genes. A comparison of chloroplast genomes of seven Solanaceae species revealed that the gene content and their relative positions of S. tuberosum are similar to the other six Solanaceae species. However, undefined open reading frames (ORFs) in LSC region were highly diverged in Solanaceae species except N. sylvestris. Detailed comparison was identified by numerous indels in the intergenic regions that were mostly located in the LSC region. Among them, a single large 241-bp deletion, was not associated with direct repeats and found in only S. tuberosum, clearly discriminates a cultivated potato from wild potato species Solanum bulbocastanum. The extent of sequence divergence may provide the basis for evaluating genetic diversity within the Solanaceae species, and will be useful to examine the evolutionary processes in potato landraces.     

Evolutionary dynamics of chloroplast genomes in subfamily Aroideae (Araceae)

Aroideae is the largest and most diverse subfamily of the plant family Araceae. Despite its agricultural and horticultural importance, the genomic resources are sparse for this subfamily. Here, we report de novo assembled and fully annotated chloroplast genomes of 13 Aroideae species. The quadripartite chloroplast genomes (size range of 158,177–170,037 bp) are comprised of a large single copy (LSC; 75,594–94,702 bp), a small single copy (SSC; 12,903–23,981 bp) and a pair of inverted repeats (IRs; 25,266–34,840 bp). Notable gene rearrangements and IRs contraction / expansions were found for Anchomanes hookeri and Zantedeschia aethiopica. Codon usage, amino acid frequencies, oligonucleotide repeats, GC contents, and gene features revealed similarities among the 13 species. The number of oligonucleotide repeats was uncorrelated with genome size or phylogenetic position of the species. Phylogenetic analyses corroborated the monophyly of Aroideae but were unable to resolve the positions of Calla and Schismatoglottis.     

The chloroplast genome sequence of Syzygium cumini (L.) and its relationship with other angiosperms

Only a few studies to date have conducted comparative genomics in the Myrtaceae family. Here, we report the complete sequence and bioinformatics analysis of the chloroplast genome of Syzygium cumini (L.), one of the family members. The size of S. cumini cp genome was within the range of reported angiosperm chloroplast genomes. Comparison of S. cumini cpDNA sequence with previously reported partial sequences of S. cumini revealed several SNPs that resulted in non-synonymous mutations in maturase K and NADH-plastoquinone oxidoreductase subunit-5. These polymorphic characters might serve as intra-specific markers to address whether lineage sorting from polymorphic ancestry has occurred. Comparison of the S. cumini chloroplast genome with related dicots revealed an expansion in the intergenic spacer located between IRA/large single copy (LSC) border and the first gene of LSC region, driven by sequence of 54 bp. This type of variation in the intergenic regions can be utilized in the development of species-specific vectors for chloroplast genetic engineering. Several of the longer (30–40 bp) repeats were found to be conserved in other dicot species, suggesting that they might be widespread in angiosperm chloroplast genomes.     

藏波罗花叶绿体基因组分析

藏波罗花（Incarvillea younghusbandii Sprague）是一种传统的补益类中药。其根作草药使用,用于滋补强壮,治产后少乳、久病虚弱、头晕、贫血等症。但目前关于藏波罗花分子遗传信息的研究很少。本研究基于高通量测序技术对藏波罗花叶绿体基因组进行测序、组装和注释,并对其序列特征、密码子偏好性、重复序列、系统发育和分化时间进行分析。结果表明,藏波罗花叶绿体基因组全长为159 323 bp,包含1个大单拷贝区（80 197 bp）、1个小单拷贝区（9 030 bp）和2个反向重复区（35 048 bp）;共注释出120个基因,包括77个蛋白编码基因、8个rRNA基因和35个tRNA基因;密码子偏好性分析显示,AAA是藏波罗花叶绿体基因组中使用最频繁的密码子;从藏波罗花叶绿体基因组中共检测到42个简单重复序列（simple sequence repeats,SSR）;系统发育分析表明,藏波罗花与密生波罗花（Incarvillea compacta）的亲缘关系最近,且在大概466万年前产生分化。本研究对藏波罗花相关资源的科学保护和开发具有重要的现实意义,也可以为后续角蒿属（Incarvillea）的物种鉴定、紫葳科（Bignoniaceae）的种群遗传多样性研究提供基本的遗传资源。     

极小种群植物川柿叶绿体基因组特征与系统发育分析（附录）

川柿(Diospyros sutchuensis)为极小种群和国家重点保护野生植物,分布范围狭窄,种群数量极少。目前,川柿基因组信息缺乏,在柿属(Diospyros)中的系统亲缘关系不明确。该研究通过Illumina平台对川柿叶绿体基因组进行测序,应用Getorganellev1.7.3.4和PGA软件对基因组进行组装和注释,使用DnaSP6.12.03软件进行多序列对比分析,并使用REPuter、Tandem Reapeats Finder和MISA软件进行重复序列分析,使用CodonW1.4和EasyCodemL软件分别进行密码子偏好性和选择压力分析。同时,基于4个不同的叶绿体基因组序列数据集,使用IQtree软件分析川柿与11个柿属物种的系统发育关系。结果表明:(1)川柿叶绿体基因组全长157 917 bp,包含1对26 111 bp的反向重复区、大单拷贝区(87 303 bp)和小单拷贝区(18 392 bp),GC碱基含量为37.4%。(2)川柿叶绿体基因组共注释到113个基因,包括79个蛋白编码基因、30个tRNA基因和4个rRNA基因; 共检测到49个长重复序列、27个串联重复序列和34个简单重复序列; 蛋白编码基因中高频密码子31个,多数密码子末位碱基为A或U,编码亮氨酸的密码子使用最多; 基因组编码区比非编码区更为保守,10个高变热点区域可作为潜在的分子标记; 蛋白编码基因中有8个基因(ndhB、ndhG、ndhI、rbcL、rpoB、petB、petD、rps12)受到正选择压力。(3)系统发育分析显示,川柿与老鸦柿(D. rhombifolia)和乌柿(D. cathayensis)亲缘关系最为密切,它们与海南柿(D. hainanensis)共同形成一个单系分支。该研究结果既为川柿及柿属种质资源鉴定、遗传多样性保护以及种群恢复等提供了叶绿体基因组资源,也为阐明川柿的系统进化提供了重要的分子信息。     

Complete Chloroplast Genome Sequences of Mongolia Medicine Artemisia frigida and Phylogenetic Relationships with Other Plants

Background

Artemisia frigida Willd. is an important Mongolian traditional medicinal plant with pharmacological functions of stanch and detumescence. However, there is little sequence and genomic information available for Artemisia frigida, which makes phylogenetic identification, evolutionary studies, and genetic improvement of its value very difficult. We report the complete chloroplast genome sequence of Artemisia frigida based on 454 pyrosequencing.

Methodology/Principal Findings

The complete chloroplast genome of Artemisia frigida is 151,076 bp including a large single copy (LSC) region of 82,740 bp, a small single copy (SSC) region of 18,394 bp and a pair of inverted repeats (IRs) of 24,971 bp. The genome contains 114 unique genes and 18 duplicated genes. The chloroplast genome of Artemisia frigida contains a small 3.4 kb inversion within a large 23 kb inversion in the LSC region, a unique feature in Asteraceae. The gene order in the SSC region of Artemisia frigida is inverted compared with the other 6 Asteraceae species with the chloroplast genomes sequenced. This inversion is likely caused by an intramolecular recombination event only occurred in Artemisia frigida. The existence of rich SSR loci in the Artemisia frigida chloroplast genome provides a rare opportunity to study population genetics of this Mongolian medicinal plant. Phylogenetic analysis demonstrates a sister relationship between Artemisia frigida and four other species in Asteraceae, including Ageratina adenophora, Helianthus annuus, Guizotia abyssinica and Lactuca sativa, based on 61 protein-coding sequences. Furthermore, Artemisia frigida was placed in the tribe Anthemideae in the subfamily Asteroideae (Asteraceae) based on ndhF and trnL-F sequence comparisons.

Conclusion

The chloroplast genome sequence of Artemisia frigida was assembled and analyzed in this study, representing the first plastid genome sequenced in the Anthemideae tribe. This complete chloroplast genome sequence will be useful for molecular ecology and molecular phylogeny studies within Artemisia species and also within the Asteraceae family.     

The Complete Sequence of the Acacia ligulata Chloroplast Genome Reveals a Highly Divergent clpP1 Gene

Legumes are a highly diverse angiosperm family that include many agriculturally important species. To date, 21 complete chloroplast genomes have been sequenced from legume crops confined to the Papilionoideae subfamily. Here we report the first chloroplast genome from the Mimosoideae, Acacia ligulata, and compare it to the previously sequenced legume genomes. The A. ligulata chloroplast genome is 158,724 bp in size, comprising inverted repeats of 25,925 bp and single-copy regions of 88,576 bp and 18,298 bp. Acacia ligulata lacks the inversion present in many of the Papilionoideae, but is not otherwise significantly different in terms of gene and repeat content. The key feature is its highly divergent clpP1 gene, normally considered essential in chloroplast genomes. In A. ligulata, although transcribed and spliced, it probably encodes a catalytically inactive protein. This study provides a significant resource for further genetic research into Acacia and the Mimosoideae. The divergent clpP1 gene suggests that Acacia will provide an interesting source of information on the evolution and functional diversity of the chloroplast Clp protease complex.     

A comparative analysis of the Lactuca and Helianthus (Asteraceae) plastid genomes: identification of divergent regions and categorization of shared repeats

We have sequenced two complete chloroplast genomes in the Asteraceae, Helianthus annuus (sunflower), and Lactuca sativa (lettuce), which belong to the distantly related subfamilies, Asteroideae and Cichorioideae, respectively. The Helianthus chloroplast genome is 151?104 bp and the Lactuca genome is 152?772 bp long, which is within the usual size range for chloroplast genomes in flowering plants. When compared to tobacco, both genomes have two inversions: a large 22.8-kb inversion and a smaller 3.3-kb inversion nested within it. Pairwise sequence divergence across all genes, introns, and spacers in Helianthus and Lactuca has resulted in the discovery of new, fast-evolving DNA sequences for use in species-level phylogenetics, such as the trnY-rpoB, trnL-rpl32, and ndhC-trnV spacers. Analysis and categorization of shared repeats resulted in seven classes useful for future repeat studies: double tandem repeats, three or more tandem repeats, direct repeats dispersed in the genome, repeats found in reverse complement orientation, hairpin loops, runs of A's or T's in excess of 12 bp, and gene or tRNA similarity. Results from BLAST searches of our genomic sequence against expressed sequence tag (EST) databases for both genomes produced eight likely RNA edited sites (C → U changes). These detailed analyses in Asteraceae contribute to a broader understanding of plastid evolution across flowering plants.     

Chloroplast genome evolution in the genus Avena

The genus Avena contains five different chloroplast genomes, I-V. A physical map of chloroplast (ct) DNA of Avena sativa (type I chloroplast genome) was constructed using three restriction endonucleases, PstI, SalI and SmaI. This genome is ca. 135.5 kbp in size, and contains two inverted repeats of ca. 22.5 kbp each, separated by a large (ca. 79.0 kbp) and small (ca. 12.5 kbp) single copy region. The rbcL gene which codes for the large subunit of ribulose 1,5-bisphosphate carboxylase, was located in the map. Restriction fragment patterns of all five chloroplast genomes were compared, and among them five fragment size and five restriction site mutations were disclosed. Four site mutations were found in two or more chloroplast genomes, the other site and five fragment size mutations were specific to one or another of the chloroplast genomes. A dendrogram showing phylogenetic relationships among the five chloroplast genomes, based on the distribution of the common and specific mutations among them, indicates that chloroplast genome divergence characterized by three restriction site mutations occurred first between two diploid groups, each carrying A and C genome (nuclear), respectively, followed by further speciation in each group.