The complete chloroplast genome sequence of Morus mongolica and a comparative analysis within the Fabidae clade

Comparative and phylogenetic analyses of the chloroplast genome reveal the taxonomy of the Morus genus

Mulberry (genus Morus) is an economically important woody plant with an altered ploidy level. The variable number of Morus species recognized by different studies indicates that the genus is in need of revision. In this study, the chloroplast (CP) genomes of 123 Morus varieties were de novo assembled and systematically analyzed. The 123 varieties represented six Morus species, namely, Morus alba, Morus nigra, Morus notabilis, Morus rubra, Morus celtidifolia, and Morus serrata. The Morus CP genome was found to be 158,969~159,548 bp in size with 125 genes, including 81 protein coding, 36 tRNA, and 8 rRNA genes. The 87 out of 123 mulberry accessions were assigned to 14 diverse groups with identical CP genome, which indicated that they are maternally inherited and share 14 common ancestors. Then 50 diverse CP genomes occurred in 123 mulberry accessions for further study. The CP genomes of the Morus genus with a quadripartite structure have two inverted repeat (IR) regions (25,654~25,702 bp) dividing the circular genome into a large single-copy (LSC) region (87,873~88,243 bp) and small single-copy (SSC) region (19,740~19,994 bp). Analysis of the phylogenetic tree constructed using the complete CP genome sequences of Morus revealed a monophyletic genus and that M. alba consisted of two clades, M. alba var. alba and M. alba var. multicaulis. The Japanese cultivated germplasms were derived from M. alba var. multicaulis. We propose that the Morus genus be classified into six species, M. nigra, M. notabilis, M. serrata, M. celtidifolia, M. rubra, and M. alba with two subspecies, M. alba var. alba and M. alba var. multicaulis. Our findings provide a valuable resource for the classification, domestication, and breeding improvement of mulberry.

The complete chloroplast genome sequences of three Broussonetia species and comparative analysis within the Moraceae

Background

Species of Broussonetia (family Moraceae) are commonly used to make textiles and high-grade paper. The distribution of Broussonetia papyrifera L. is considered to be related to the spread and location of humans. The complete chloroplast (cp) genomes of B. papyrifera, Broussonetia kazinoki Sieb., and Broussonetia kaempferi Sieb. were analyzed to better understand the status and evolutionary biology of the genus Broussonetia.

Methods

The cp genomes were assembled and characterized using SOAPdenovo2 and DOGMA. Phylogenetic and molecular dating analysis were performed using the concatenated nucleotide sequences of 35 species in the Moraceae family and were based on 66 protein-coding genes (PCGs). An analysis of the sequence divergence (pi) of each PCG among the 35 cp genomes was conducted using DnaSP v6. Codon usage indices were calculated using the CodonW program.

Results

All three cp genomes had the typical land plant quadripartite structure, ranging in size from 160,239 bp to 160,841 bp. The ribosomal protein L22 gene (RPL22) was either incomplete or missing in all three Broussonetia species. Phylogenetic analysis revealed two clades. Clade 1 included Morus and Artocarpus, whereas clade 2 included the other seven genera. Malaisia scandens Lour. was clustered within the genus Broussonetia. The differentiation of Broussonetia was estimated to have taken place 26 million years ago. The PCGs' pi values ranged from 0.0005 to 0.0419, indicating small differences within the Moraceae family. The distribution of most of the genes in the effective number of codons plot (ENc-plot) fell on or near the trend line; the slopes of the trend line of neutrality plots were within the range of 0.0363-0.171. These results will facilitate the identification, taxonomy, and utilization of the Broussonetia species and further the evolutionary studies of the Moraceae family.

Complete chloroplast genome sequencing and comparative analysis reveals changes to the chloroplast genome after allopolyploidization in Cucumis

Allopolyploids undergo "genomic shock" leading to significant genetic and epigenetic modifications. Previous studies have mainly focused on nuclear changes, while little is known about the inheritance and changes of organelle genome in allopolyploidization. The synthetic allotetraploid Cucumis ×hytivus, which is generated via hybridization between C. hystrix and C. sativus, is a useful model system for studying cytonuclear variation. Here, we report the chloroplast genome of allotetraploid C. ×hytivus and its diploid parents via sequencing and comparative analysis. The size of the obtained chloroplast genomes ranged from 154 673 to 155 760 bp, while their gene contents, gene orders, and GC contents were similar to each other. Comparative genome analysis supports chloroplast maternal inheritance. However, we identified 51 indels and 292 SNP genetic variants in the chloroplast genome of the allopolyploid C. ×hytivus relative to its female parent C. hystrix. Nine intergenic regions with rich variation were identified through comparative analysis of the chloroplast genomes within the subgenus Cucumis. The phylogenetic network based on the chloroplast genome sequences clarified the evolution and taxonomic position of the synthetic allotetraploid C. ×hytivus. The results of this study provide us with an insight into the changes of organelle genome after allopolyploidization, and a new understanding of the cytonuclear evolution.

The Complete Plastid Genome of Artocarpus camansi: A High Degree of Conservation of the Plastome Structure in the Family Moraceae

Understanding the plastid genome is extremely important for the interpretation of the genetic mechanisms associated with essential physiological and metabolic functions, the identification of possible marker regions for phylogenetic or phylogeographic analyses, and the elucidation of the modes through which natural selection operates in different regions of this genome. In the present study, we assembled the plastid genome of Artocarpus camansi, compared its repetitive structures with Artocarpus heterophyllus, and searched for evidence of synteny within the family Moraceae. We also constructed a phylogeny based on 56 chloroplast genes to assess the relationships among three families of the order Rosales, that is, the Moraceae, Rhamnaceae, and Cannabaceae. The plastid genome of A. camansi has 160,096 bp, and presents the typical circular quadripartite structure of the Angiosperms, comprising a large single copy (LSC) of 88,745 bp and a small single copy (SSC) of 19,883 bp, separated by a pair of inverted repeat (IR) regions each with a length of 25,734 bp. The total GC content was 36.0%, which is very similar to Artocarpus heterophyllus (36.1%) and other moraceous species. A total of 23,068 codons and 80 SSRs were identified in the A. camansi plastid genome, with the majority of the SSRs being mononucleotide (70.0%). A total of 50 repeat structures were observed in the A. camansi plastid genome, in contrast with 61 repeats in A. heterophyllus. A purifying selection signal was found in 70 of the 79 protein-coding genes, indicating that they have all been highly conserved throughout the evolutionary history of the genus. The comparative analysis of the structural characteristics of the chloroplast among different moraceous species found a high degree of similarity in the sequences, which indicates a highly conserved evolutionary model in these plastid genomes. The phylogenetic analysis also recovered a high degree of similarity between the chloroplast genes of A. camansi and A. heterophyllus, and reconfirmed the hypothesis of the intense conservation of the plastome in the family Moraceae.

Molecular phylogeny of mulberries reconstructed from ITS and two cpDNA sequences

Background

Species in the genus Morus (Moraceae) are deciduous woody plants of great economic importance. The classification and phylogenetic relationships of Morus, especially the abundant mulberry resources in China, is still undetermined. Internal transcribed spacer (ITS) regions are among the most widely used molecular markers in phylogenetic analyses of angiosperms. However, according to the previous phylogenetic analyses of ITS sequences, most of the mulberry accessions collected in China were grouped into the largest clade lacking for phylogenetic resolution. Compared with functional ITS sequences, ITS pseudogenes show higher sequence diversity, so they can provide useful phylogenetic information.

Methods

We sequenced the ITS regions and the chloroplast DNA regions TrnL-TrnF and TrnT-TrnL from 33 mulberry accessions, and performed phylogenetic analyses to explore the evolution of mulberry.

Results

We found ITS pseudogenes in 11 mulberry accessions. In the phylogenetic tree constructed from ITS sequences, clade B was separated into short-type sequence clades (clades 1 and 2), and a long-type sequence clade (clade 3). Pseudogene sequences were separately clustered into two pseudogroups, designated as pseudogroup 1 and pseudogroup 2. The phylogenetic tree generated from cpDNA sequences also separated clade B into two clades.

Conclusions

Two species were separated in clade B. The existence of three connection patterns and incongruent distribution patterns between the phylogenetic trees generated from cpDNA and ITS sequences suggested that the ITS pseudogene sequences connect with genetic information from the female progenitor. Hybridization has played important roles in the evolution of mulberry, resulting in low resolution of the phylogenetic analysis based on ITS sequences. An evolutionary pattern illustrating the evolution history of mulberry is proposed. These findings have significance for the conservation of local mulberry resources. Polyploidy, hybridization, and concerted evolution have all played the roles in the evolution of ITS sequences in mulberry. This study will expand our understanding of mulberry evolution.

Complete chloroplast genome sequence and phylogenetic analysis of Spathiphyllum 'Parrish'

Spathiphyllum is a very important tropical plant used as a small, potted, ornamental plant in South China, with an annual output value of hundreds of millions of yuan. In this study, we sequenced and analyzed the complete nucleotide sequence of the Spathiphyllum 'Parrish' chloroplast genome. The whole chloroplast genome is 168,493 bp in length, and includes a pair of inverted repeat (IR) regions (IRa and IRb, each 31,600 bp), separated by a small single-copy (SSC, 15,799 bp) region and a large single-copy (LSC, 89,494 bp) region. Our annotation revealed that the S. 'Parrish' chloroplast genome contained 132 genes, including 87 protein coding genes, 37 transfer RNA genes, and 8 ribosomal RNA genes. In the repeat structure analysis, we detected 281 simple sequence repeats (SSRs) which included mononucleotides (223), dinucleotides (28), trinucleotides (12), tetranucleotides (11), pentanucleotides (6), and hexanucleotides (1), in the S. 'Parrish' chloroplast genome. In addition, we identified 50 long repeats, comprising 18 forward repeats, 13 reverse repeats, 17 palindromic repeats, and 2 complementary repeats. Single nucleotide polymorphism (SNP) and insertion/deletion (indel) analyses of the chloroplast genome of the S. 'Parrish' relative S. cannifolium revealed 962 SNPs in S. 'Parrish'. There were 158 indels (90 insertions and 68 deletions) in the S. 'Parrish' chloroplast genome relative to the S. cannifolium chloroplast genome. Phylogenetic analysis of five species found S. 'Parrish' to be more closely related to S. kochii than to S. cannifolium. This study identified the characteristics of the S. 'Parrish' chloroplast genome, which will facilitate species identification and phylogenetic analysis within the genus Spathiphyllum.

Determination of the evolutionary pressure on Camellia oleifera on Hainan Island using the complete chloroplast genome sequence

Camellia oleifera is one of the four largest woody edible oil plants in the world with high ecological and medicinal values. Due to frequent interspecific hybridization, it was difficult to study its genetics and evolutionary history. This study used C. oleifera that was collected on Hainan Island to conduct our research. The unique island environment makes the quality of tea oil higher than that of other species grown in the mainland. Moreover, a long-term geographic isolation might affect gene structure. In order to better understand the molecular biology of this species, protect excellent germplasm resources, and promote the population genetics and phylogenetic studies of Camellia plants, high-throughput sequencing technology was used to obtain the chloroplast genome sequence of Hainan C. oleifera. The results showed that the whole chloroplast genome of C. oleifera in Hainan was 156,995 bp in length, with a typical quadripartite structure of a large single copy (LSC) region of 86,648 bp, a small single copy (SSC) region of 18,297 bp, and a pair of inverted repeats (IRs) of 26,025 bp. The whole genome encoded a total of 141 genes (115 different genes), including 88 protein-coding genes, 45 tRNA genes, and eight rRNA genes. Among these genes, nine genes contained one intron, two genes contained two introns, and four overlapping genes were also detected. The total GC content of Hainan C. oleifera's chloroplast genome was 37.29%. The chloroplast genome structure characteristics of Hainan C. oleifera were compared with mainland C. oleifera and those of the other eight closely related Theaceae species; it was found that the contractions and expansions of the IR/LSC and IR/SSC regions affected the length of chloroplast genome. The chloroplast genome sequences of these Theaceae species were highly similar. A comparative analysis indicated that the Theaceae species were conserved in structure and evolution. A total of 51 simple sequence repeat (SSR) loci were detected in the chloroplast genome of Hainan C. oleifera, and all Camellia plants did not have pentanucleotide repeats, which could be used as a good marker in phylogenetic studies. We also detected seven long repeats, the base composition of all repeats was biased toward A/T, which was consistent with the codon bias. It was found that Hainan C. oleifera had a similar evolutionary relationship with C. crapnelliana, through the use of codons and phylogenetic analysis. This study can provide an effective genomic resource for the evolutionary history of Theaceae family.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Aster tataricus

We sequenced and analyzed the complete chloroplast genome of Aster tataricus (family Asteraceae), a Chinese herb used medicinally to relieve coughs and reduce sputum. The A. tataricus chloroplast genome was 152,992 bp in size, and harbored a pair of inverted repeat regions (IRa and IRb, each 24,850 bp) divided into a large single-copy (LSC, 84,698 bp) and a small single-copy (SSC, 18,250 bp) region. Our annotation revealed that the A. tataricus chloroplast genome contained 115 genes, including 81 protein-coding genes, 4 ribosomal RNA genes, and 30 transfer RNA genes. In addition, 70 simple sequence repeats (SSRs) were detected in the A. tataricus chloroplast genome, including mononucleotides (36), dinucleotides (1), trinucleotides (23), tetranucleotides (1), pentanucleotides (8), and hexanucleotides (1). Comparative chloroplast genome analysis of three Aster species indicated that a higher similarity was preserved in the IR regions than in the LSC and SSC regions, and that the differences in the degree of preservation were slighter between A. tataricus and A. altaicus than between A. tataricus and A. spathulifolius. Phylogenetic analysis revealed that A. tataricus was more closely related to A. altaicus than to A. spathulifolius. Our findings offer valuable information for future research on Aster species identification and selective breeding.

Mulberry (Morus spp.) has the features to treat as a potential perennial model system

Mulberry (Morus spp.), a commercially exploited tree species as the host of monophagous pest silk worm (Bombyx mori), belongs to the family Moraceae. The domesticated tree has diverse beneficial characters such as traits associated with rapid growth and biomass production, plant insect/microbe interaction, abiotic stress tolerance and the traits associated with nutritional and medicinal values; some of which have been exploited. Draft genome of Morus notabilis has been sequenced and a large volume of transcriptome and genomic resources have been generated. In this review an attempt has been made to examine the options for considering mulberry as another tree model system to study unique traits associated with perennial systems. The diverse traits and features in mulberry suggest that the system can be a "comprehensive trait integrated tree system" quite different from other model tree systems.

The Complete Chloroplast Genome Sequences of Fritillaria ussuriensis Maxim. and Fritillaria cirrhosa D. Don, and Comparative Analysis with Other Fritillaria Species

The genus Fritillaria belongs to the widely distributed Liliaceae. The bulbs of Fritillaria, F. ussuriensis and F. cirrhosa are valuable herbaceous medicinal ingredients. However, they are still used indiscriminately in herbal medicine. Identification and molecular phylogenic analysis of Fritillaria species are therefore required. Here, we report the complete chloroplast (CP) genome sequences of F. ussuriensis and F. cirrhosa. The two Fritillaria CP genomes were 151,524 and 151,083 bp in length, respectively, and each included a pair of inverted repeated regions (52,678 and 52,156 bp) that was separated by a large single copy region (81,732 and 81,390 bp), and a small single copy region (17,114 and 17,537 bp). A total of 111 genes in F. ussuriensis and 112 in F. cirrhosa comprised 77 protein-coding regions in F. ussuriensis and 78 in F. cirrhosa, 30 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes. The gene order, content, and orientation of the two Fritillaria CP genomes exhibited the general structure of flowering plants, and were similar to those of other Fritillaria species. Comparison of the six Fritillaria species’ CP genomes indicated seven highly divergent regions in intergenic spacers and in the matK, rpoC1, rpoC2, ycf1, ycf2, ndhD, and ndhF coding regions. We established the position of the six species through phylogenic analysis. The complete chloroplast genome sequences of the two Fritillaria species and a comparison study are useful genomic information for identifying and for studying the phylogenetic relationship among Fritillaria species within the Liliaceae.

The complete chloroplast genome sequence of Morus cathayana and Morus multicaulis, and comparative analysis within genus Morus L

Trees in the Morus genera belong to the Moraceae family. To better understand the species status of genus Morus and to provide information for studies on evolutionary biology within the genus, the complete chloroplast (cp) genomes of M. cathayana and M. multicaulis were sequenced. The plastomes of the two species are 159,265 bp and 159,103 bp, respectively, with corresponding 83 and 82 simple sequence repeats (SSRs). Similar to the SSRs of M. mongolica and M. indica cp genomes, more than 70% are mononucleotides, ten are in coding regions, and one exhibits nucleotide content polymorphism. Results for codon usage and relative synonymous codon usage show a strong bias towards NNA and NNT codons in the two cp genomes. Analysis of a plot of the effective number of codons (ENc) for five Morus spp. cp genomes showed that most genes follow the standard curve, but several genes have ENc values below the expected curve. The results indicate that both natural selection and mutational bias have contributed to the codon bias. Ten highly variable regions were identified among the five Morus spp. cp genomes, and 154 single-nucleotide polymorphism mutation events were accurately located in the gene coding region.

The Chloroplast Genome of Passiflora edulis (Passifloraceae) Assembled from Long Sequence Reads: Structural Organization and Phylogenomic Studies in Malpighiales

The family Passifloraceae consists of some 700 species classified in around 16 genera. Almost all its members belong to the genus Passiflora. In Brazil, the yellow passion fruit (Passiflora edulis) is of considerable economic importance, both for juice production and consumption as fresh fruit. The availability of chloroplast genomes (cp genomes) and their sequence comparisons has led to a better understanding of the evolutionary relationships within plant taxa. In this study, we obtained the complete nucleotide sequence of the P. edulis chloroplast genome, the first entirely sequenced in the Passifloraceae family. We determined its structure and organization, and also performed phylogenomic studies on the order Malpighiales and the Fabids clade. The P. edulis chloroplast genome is characterized by the presence of two copies of an inverted repeat sequence (IRA and IRB) of 26,154 bp, each separating a small single copy region of 13,378 bp and a large single copy (LSC) region of 85,720 bp. The annotation resulted in the identification of 105 unique genes, including 30 tRNAs, 4 rRNAs, and 71 protein coding genes. Also, 36 repetitive elements and 85 SSRs (microsatellites) were identified. The structure of the complete cp genome of P. edulis differs from that of other species because of rearrangement events detected by means of a comparison based on 22 members of the Malpighiales. The rearrangements were three inversions of 46,151, 3,765 and 1,631 bp, located in the LSC region. Phylogenomic analysis resulted in strongly supported trees, but this could also be a consequence of the limited taxonomic sampling used. Our results have provided a better understanding of the evolutionary relationships in the Malpighiales and the Fabids, confirming the potential of complete chloroplast genome sequences in inferring evolutionary relationships and the utility of long sequence reads for generating very accurate biological information.

The Complete Chloroplast Genome of Ye-Xing-Ba (Scrophularia dentata; Scrophulariaceae), an Alpine Tibetan Herb

Scrophularia dentata is an important Tibetan medicinal plant and traditionally used for the treatment of exanthema and fever in Traditional Tibetan Medicine (TTM). However, there is little sequence and genomic information available for S. dentata. In this paper, we report the complete chloroplast genome sequence of S. dentata and it is the first sequenced member of the Sect. Tomiophyllum within Scrophularia (Scrophulariaceae). The gene order and organization of the chloroplast genome of S. dentata are similar to other Lamiales chloroplast genomes. The plastome is 152,553 bp in length and includes a pair of inverted repeats (IRs) of 25,523 bp that separate a large single copy (LSC) region of 84,058 bp and a small single copy (SSC) region of 17,449 bp. It has 38.0% GC content and includes 114 unique genes, of which 80 are protein-coding, 30 are transfer RNA, and 4 are ribosomal RNA. Also, it contains 21 forward repeats, 19 palindrome repeats and 41 simple sequence repeats (SSRs). The repeats and SSRs within S. dentata were compared with those of S. takesimensis and present certain discrepancies. The chloroplast genome of S. dentata was compared with other five publicly available Lamiales genomes from different families. All the coding regions and non-coding regions (introns and intergenic spacers) within the six chloroplast genomes have been extracted and analysed. Furthermore, the genome divergent hotspot regions were identified. Our studies could provide basic data for the alpine medicinal species conservation and molecular phylogenetic researches of Scrophulariaceae and Lamiales.

Comparative Analysis of the Complete Chloroplast Genomes of Five Quercus Species

Quercus is considered economically and ecologically one of the most important genera in the Northern Hemisphere. Oaks are taxonomically perplexing because of shared interspecific morphological traits and intraspecific morphological variation, which are mainly attributed to hybridization. Universal plastid markers cannot provide a sufficient number of variable sites to explore the phylogeny of this genus, and chloroplast genome-scale data have proven to be useful in resolving intractable phylogenetic relationships. In this study, the complete chloroplast genomes of four Quercus species were sequenced, and one published chloroplast genome of Quercus baronii was retrieved for comparative analyses. The five chloroplast genomes ranged from 161,072 bp (Q. baronii) to 161,237 bp (Q. dolicholepis) in length, and their gene organization and order, and GC content, were similar to those of other Fagaceae species. We analyzed nucleotide substitutions, indels, and repeats in the chloroplast genomes, and found 19 relatively highly variable regions that will potentially provide plastid markers for further taxonomic and phylogenetic studies within Quercus. We observed that four genes (ndhA, ndhK, petA, and ycf1) were subject to positive selection. The phylogenetic relationships of the Quercus species inferred from the chloroplast genomes obtained moderate-to-high support, indicating that chloroplast genome data may be useful in resolving relationships in this genus.

Complete chloroplast genome sequence of a major economic species, Ziziphus jujuba (Rhamnaceae)

Ziziphus jujuba is an important woody plant with high economic and medicinal value. Here, we analyzed and characterized the complete chloroplast (cp) genome of Z. jujuba, the first member of the Rhamnaceae family for which the chloroplast genome sequence has been reported. We also built a web browser for navigating the cp genome of Z. jujuba ( http://bio.njfu.edu.cn/gb2/gbrowse/Ziziphus_jujuba_cp/ ). Sequence analysis showed that this cp genome is 161,466 bp long and has a typical quadripartite structure of large (LSC, 89,120 bp) and small (SSC, 19,348 bp) single-copy regions separated by a pair of inverted repeats (IRs, 26,499 bp). The sequence contained 112 unique genes, including 78 protein-coding genes, 30 transfer RNAs, and four ribosomal RNAs. The genome structure, gene order, GC content, and codon usage are similar to other typical angiosperm cp genomes. A total of 38 tandem repeats, two forward repeats, and three palindromic repeats were detected in the Z. jujuba cp genome. Simple sequence repeat (SSR) analysis revealed that most SSRs were AT-rich. The homopolymer regions in the cp genome of Z. jujuba were verified and manually corrected by Sanger sequencing. One-third of mononucleotide repeats were found to be erroneously sequenced by the 454 pyrosequencing, which resulted in sequences of 1-4 bases shorter than that by the Sanger sequencing. Analyzing the cp genome of Z. jujuba revealed that the IR contraction and expansion events resulted in ycf1 and rps19 pseudogenes. A phylogenetic analysis based on 64 protein-coding genes showed that Z. jujuba was closely related to members of the Elaeagnaceae family, which will be helpful for phylogenetic studies of other Rosales species. The complete cp genome sequence of Z. jujuba will facilitate population, phylogenetic, and cp genetic engineering studies of this economic plant.

Complete chloroplast genomes of Aegilops tauschii Coss. and Ae. cylindrica Host sheds light on plasmon D evolution

Hexaploid wheat (Triticum aestivum L., genomes AABBDD) originated in South Caucasus by allopolyploidization of the cultivated Emmer wheat T. dicoccum (genomes AABB) with the Caucasian Ae. tauschii ssp strangulata (genomes DD). Genetic variation of Ae. tauschii is an important natural resource, that is why it is of particular importance to investigate how this variation was formed during Ae. tauschii evolutionary history and how it is presented through the species area. The D genome is also found in tetraploid Ae. cylindrica Host (2n = 28, CCDD). The plasmon diversity that exists in Triticum and Aegilops species is of great significance for understanding the evolution of these genera. In the present investigation the complete nucleotide sequence of plasmon D (chloroplast DNA) of nine accessions of Ae. tauschii and two accessions of Ae. cylindrica are presented. Twenty-eight SNPs are characteristic for both TauL1 and TauL2 accessions of Ae. tauschii using TauL3 as a reference. Four SNPs are additionally observed for TauL2 lineage. The longest (27 bp) indel is located in the intergenic spacer Rps15-ndhF of SSC. This indel can be used for simple determination of TauL3 lineage among Ae. tauschii accessions. In the case of Ae. cylindrica additionally 7 SNPs were observed. The phylogeny tree shows that chloroplast DNA of TauL1 and TauL2 diverged from the TauL3 lineage. TauL1 lineage is relatively older then TauL2. The position of Ae. cylindrica accessions on Ae. tauschii phylogeny tree constructed on chloroplast DNA variation data is intermediate between TauL1 and TauL2. The complete nucleotide sequence of chloroplast DNA of Ae. tauschii and Ae. cylindrica allows to refine the origin and evolution of D plasmon of genus Aegilops.