Complete Chloroplast Genome Sequences and Comparative Analysis of Chenopodium quinoa and C. album

An updated phylogeny and adaptive evolution within Amaranthaceae s.l. inferred from multiple phylogenomic datasets

Amaranthaceae s.l. is a widely distributed family consisting of over 170 genera and 2000 species. Previous molecular phylogenetic studies have shown that Amaranthaceae s.s. and traditional Chenopodiaceae form a monophyletic group (Amaranthaceae s.l.), however, the relationships within this evolutionary branch have yet to be fully resolved. In this study, we assembled the complete plastomes and full-length ITS of 21 Amaranthaceae s.l. individuals and compared them with 38 species of Amaranthaceae s.l. Through plastome structure and sequence alignment analysis, we identified a reverse complementary region approximately 5200 bp long in the genera Atriplex and Chenopodium. Adaptive evolution analysis revealed significant positive selection in eight genes, which likely played a driving role in the evolution of Amaranthaceae s.l., as demonstrated by partitioned evolutionary analysis. Furthermore, we found that about two-thirds of the examined species lack the ycf15 gene, potentially associated with natural selection pressures from their adapted habitats. The phylogenetic tree indicated that some genera (Chenopodium, Halogeton, and Subtr. Salsolinae) are paraphyletic lineages. Our results strongly support the clustering of Amaranthaceae s.l. with monophyletic traditional Chenopodiaceae (Clades I and II) and Amaranthaceae s.s. After a comprehensive analysis, we determined that cytonuclear conflict, gene selection by adapted habitats, and incomplete lineage sorting (ILS) events were the primary reasons for the inconsistent phylogeny of Amaranthaceae s.l. During the last glacial period, certain species within Amaranthaceae s.l. underwent adaptations to different environments and began to differentiate rapidly. Since then, these species may have experienced morphological and genetic changes distinct from those of other genera due to intense selection pressure.

Phylogenetic placement of the monotypic Baolia (Amaranthaceae s.l.) based on morphological and molecular evidence

BACKGROUND

Baolia H.W.Kung & G.L.Chu is a monotypic genus only known in Diebu County, Gansu Province, China. Its systematic position is contradictory, and its morphoanatomical characters deviate from all other Chenopodiaceae. Recent study has regarded Baolia as a sister group to Corispermoideae. We therefore sequenced and compared the chloroplast genomes of this species, and resolved its phylogenetic position based on both chloroplast genomes and marker sequences.

RESULTS

We sequenced 18 chloroplast genomes of 16 samples from two populations of Baolia bracteata and two Corispermum species. These genomes of Baolia ranged in size from 152,499 to 152,508 bp. Simple sequence repeats (SSRs) were primarily located in the LSC region of Baolia chloroplast genomes, and most of them consisted of single nucleotide A/T repeat sequences. Notably, there were differences in the types and numbers of SSRs between the two populations of B. bracteata. Our phylogenetic analysis, based on both complete chloroplast genomes from 33 species and a combination of three markers (ITS, rbcL, and matK) from 91 species, revealed that Baolia and Corispermoideae (Agriophyllum, Anthochlamys, and Corispermum) form a well-supported clade and sister to Acroglochin. According to our molecular dating results, a major divergence event between Acroglochin, Baolia, and Corispermeae occurred during the Middle Eocene, approximately 44.49 mya. Ancestral state reconstruction analysis showed that Baolia exhibited symplesiomorphies with those found in core Corispermoideae characteristics including pericarp and seed coat.

CONCLUSIONS

Comparing the chloroplast genomes of B. bracteata with those of eleven typical Chenopodioideae and Corispermoideae species, we observed a high overall similarity and a one notable noteworthy case of inversion of approximately 3,100 bp. of DNA segments only in two Atriplex and four Chenopodium species. We suggest that Corispermoideae should be considered in a broader sense, it includes Corispermeae (core Corispermoideae: Agriophyllum, Anthochlamys, and Corispermum), as well as two new monotypic tribes, Acroglochineae (Acroglochin) and Baolieae (Baolia).

Complete plastid genome of Iris orchioides and comparative analysis with 19 Iris plastomes

Iris is a cosmopolitan genus comprising approximately 280 species distributed throughout the Northern Hemisphere. Although Iris is the most diverse group in the Iridaceae, the number of taxa is debatable owing to various taxonomic issues. Plastid genomes have been widely used for phylogenetic research in plants; however, only limited number of plastid DNA markers are available for phylogenetic study of the Iris. To understand the genomic features of plastids within the genus, including its structural and genetic variation, we newly sequenced and analyzed the complete plastid genome of I. orchioides and compared it with those of 19 other Iris taxa. Potential plastid markers for phylogenetic research were identified by computing the sequence divergence and phylogenetic informativeness. We then tested the utility of the markers with the phylogenies inferred from the markers and whole-plastome data. The average size of the plastid genome was 152,926 bp, and the overall genomic content and organization were nearly identical among the 20 Iris taxa, except for minor variations in the inverted repeats. We identified 10 highly informative regions (matK, ndhF, rpoC2, ycf1, ycf2, rps15-ycf, rpoB-trnC, petA-psbJ, ndhG-ndhI and psbK-trnQ) and inferred a phylogeny from each region individually, as well as from their concatenated data. Remarkably, the phylogeny reconstructed from the concatenated data comprising three selected regions (rpoC2, ycf1 and ycf2) exhibited the highest congruence with the phylogeny derived from the entire plastome dataset. The result suggests that this subset of data could serve as a viable alternative to the complete plastome data, especially for molecular diagnoses among closely related Iris taxa, and at a lower cost.

Phylogeny and diversification of genus Sanicula L. (Apiaceae): novel insights from plastid phylogenomic analyses

BACKGROUND

The genus Sanicula L. is a unique perennial herb that holds important medicinal values. Although the previous studies on Sanicula provided us with a good research basis, its taxonomic system and interspecific relationships have not been satisfactorily resolved, especially for those endemic to China. Moreover, the evolutionary history of this genus also remains inadequately understood. The plastid genomes possessing highly conserved structure and limited evolutionary rate have proved to be an effective tool for studying plant phylogeny and evolution.

RESULTS

In the current study, we newly sequenced and assembled fifteen Sanicula complete plastomes. Combined with two previously reported plastomes, we performed comprehensively plastid phylogenomics analyses to gain novel insights into the evolutionary history of this genus. The comparative results indicated that the seventeen plastomes exhibited a high degree of conservation and similarity in terms of their structure, size, GC content, gene order, IR borders, codon bias patterns and SSRs profiles. Such as all of them displayed a typical quadripartite structure, including a large single copy region (LSC: 85,074-86,197 bp), a small single copy region (SSC: 17,047-17,132 bp) separated by a pair of inverted repeat regions (IRs: 26,176-26,334 bp). And the seventeen plastomes had similar IR boundaries and the adjacent genes were identical. The rps19 gene was located at the junction of the LSC/IRa, the IRa/SSC junction region was located between the trnN gene and ndhF gene, the ycf1 gene appeared in the SSC/IRb junction and the IRb/LSC boundary was located between rpl12 gene and trnH gene. Twelve specific mutation hotspots (atpF, cemA, accD, rpl22, rbcL, matK, ycf1, trnH-psbA, ycf4-cemA, rbcL-accD, trnE-trnT and trnG-trnR) were identified that can serve as potential DNA barcodes for species identification within the genus Sanicula. Furthermore, the plastomes data and Internal Transcribed Spacer (ITS) sequences were performed to reconstruct the phylogeny of Sanicula. Although the tree topologies of them were incongruent, both provided strong evidence supporting the monophyly of Saniculoideae and Apioideae. In addition, the sister groups between Saniculoideae and Apioideae were strongly suggested. The Sanicula species involved in this study were clustered into a clade, and the Eryngium species were also clustered together. However, it was clearly observed that the sections of Sanicula involved in the current study were not respectively recovered as monophyletic group. Molecular dating analysis explored that the origin of this genus was occurred during the late Eocene period, approximately 37.84 Ma (95% HPD: 20.33-52.21 Ma) years ago and the diversification of the genus was occurred in early Miocene 18.38 Ma (95% HPD: 10.68-25.28 Ma).

CONCLUSION

The plastome-based tree and ITS-based tree generated incongruences, which may be attributed to the event of hybridization/introgression, incomplete lineage sorting (ILS) and chloroplast capture. Our study highlighted the power of plastome data to significantly improve the phylogenetic supports and resolutions, and to efficiently explore the evolutionary history of this genus. Molecular dating analysis explored that the diversification of the genus occurred in the early Miocene, which was largely influenced by the prevalence of the East Asian monsoon and the uplift of the Hengduan Mountains (HDM). In summary, our study provides novel insights into the plastome evolution, phylogenetic relationships, taxonomic framework and evolution of genus Sanicula.

Comparative Analysis of Six Chloroplast Genomes in Chenopodium and Its Related Genera (Amaranthaceae): New Insights into Phylogenetic Relationships and the Development of Species-Specific Molecular Markers

Species within the genus Chenopodium hold significant research interest due to their nutritional richness and salt tolerance. However, the morphological similarities among closely related species and a dearth of genomic resources have impeded their comprehensive study and utilization. In the present research, we conduct the sequencing and assembly of chloroplast (cp) genomes from six Chenopodium and related species, five of which were sequenced for the first time. These genomes ranged in length from 151,850 to 152,215 base pairs, showcased typical quadripartite structures, and encoded 85 protein-coding genes (PCGs), 1 pseudogene, 37 tRNA genes, and 8 rRNA genes. Compared with the previously published sequences of related species, these cp genomes are relatively conservative, but there are also some interspecific differences, such as inversion and IR region contraction. We discerned 929 simple sequence repeats (SSRs) and a series of highly variable regions across 16 related species, predominantly situated in the intergenic spacer (IGS) region and introns. The phylogenetic evaluations revealed that Chenopodium is more closely related to genera such as Atriplex, Beta, Dysphania, and Oxybase than to other members of the Amaranthaceae family. These lineages shared a common ancestor approximately 60.80 million years ago, after which they diverged into distinct genera. Based on InDels and SNPs between species, we designed 12 pairs of primers for species identification, and experiments confirmed that they could completely distinguish 10 related species.

Loss of plastid ndh genes in an autotrophic desert plant

Plant plastid genomes are highly conserved with most flowering plants having the same complement of essential plastid genes. Here, we report the loss of five of the eleven NADH dehydrogenase subunit genes (ndh) in the plastid of a desert plant jojoba (Simmondsia chinensis). The plastid genome of jojoba was 156,496 bp with one large single copy region (LSC), a very small single copy region (SSC) and two expanded inverted repeats (IRA + IRB). The NADH dehydrogenase (NDH) complex is comprised of several protein subunits, encoded by the ndh genes of the plastome and the nucleus. The ndh genes are critical to the proper functioning of the photosynthetic electron transport chain and protection of plants from oxidative stress. Most plants are known to contain all eleven ndh genes. Plants with missing or defective ndh genes are often heterotrophs either due to their complete or holo- or myco- parasitic nature. Plants with a defective NDH complex, caused by the deletion/pseudogenisation of some or all the ndh genes, survive in milder climates suggesting the likely extinction of plant lineages lacking these genes under harsh climates. Interestingly, some autotrophic plants do exist without ndh gene/s and can cope with high or low light. This implies that these plants are protected from oxidative stress by mechanisms excluding ndh genes. Jojoba has evolved mechanisms to cope with a non-functioning NDH complex and survives in extreme desert conditions with abundant sunlight and limited water.

Genetic resources and breeding approaches for improvement of amaranth (Amaranthus spp.) and quinoa (Chenopodium quinoa)

Nowadays, the human population is more concerned about their diet and very specific in choosing their food sources to ensure a healthy lifestyle and avoid diseases. So people are shifting to more smart nutritious food choices other than regular cereals and staple foods they have been eating for a long time. Pseudocereals, especially, amaranth and quinoa, are important alternatives to traditional cereals due to comparatively higher nutrition, essential minerals, amino acids, and zero gluten. Both Amaranchaceae crops are low-input demanding and hardy plants tolerant to stress, drought, and salinity conditions. Thus, these crops may benefit developing countries that follow subsistence agriculture and have limited farming resources. However, these are underutilized orphan crops, and the efforts to improve them by reducing their saponin content remain ignored for a long time. Furthermore, these crops have very rich variability, but the progress of their genetic gain for getting high-yielding genotypes is slow. Realizing problems in traditional cereals and opting for crop diversification to tackle climate change, research should be focused on the genetic improvement for low saponin, nutritionally rich, tolerant to biotic and abiotic stresses, location-specific photoperiod, and high yielding varietal development of amaranth and quinoa to expand their commercial cultivation. The latest technologies that can accelerate the breeding to improve yield and quality in these crops are much behind and slower than the already established major crops of the world. We could learn from past mistakes and utilize the latest trends such as CRISPR/Cas, TILLING, and RNA interference (RNAi) technology to improve these pseudocereals genetically. Hence, the study reviewed important nutrition quality traits, morphological descriptors, their breeding behavior, available genetic resources, and breeding approaches for these crops to shed light on future breeding strategies to develop superior genotypes.

Development of InDels markers for the identification of cytoplasmic male sterility in Sorghum by complete chloroplast genome sequences analysis

Cytoplasmic male sterility (CMS) is predominantly used for F1 hybrid breeding and seed production in Sorghum. DNA markers to distinguish between normal fertile (CMS-N) and sterile (CMS-S) male cytoplasm can facilitate F1 hybrid cultivar development in Sorghum breeding programs. In this study, the complete chloroplast (cp) genome sequences of CMS-S and Korean Sorghum cultivars were obtained using next-generation sequencing. The de novo assembled genome size of ATx623, the CMS-S line of the chloroplast, was 140,644bp. When compared to the CMS-S and CMS-N cp genomes, 19 single nucleotide polymorphisms (SNPs) and 142 insertions and deletions (InDels) were identified, which can be used for marker development for breeding, population genetics, and evolution studies. Two InDel markers with sizes greater than 20 bp were developed to distinguish cytotypes based on the copy number variation of lengths as 28 and 22 bp tandem repeats, respectively. Using the newly developed InDel markers with five pairs of CMS-S and their near isogenic maintainer line, we were able to easily identify their respective cytotypes. The InDel markers were further examined and applied to 1,104 plants from six Korean Sorghum cultivars to identify variant cytotypes. Additionally, the phylogenetic analysis of seven Sorghum species with complete cp genome sequences, including wild species, indicated that CMS-S and CMS-N contained Milo and Kafir cytotypes that might be hybridized from S. propinquum and S. sudanese, respectively. This study can facilitate F1 hybrid cultivar development by providing breeders with reliable tools for marker-assisted selection to breed desirable Sorghum varieties.

Genome-wide identification and characterization of glutathione S-transferase gene family in quinoa (Chenopodium quinoa Willd.)

The present investigation was envisaged for large scale in-silico genome wide identification and characterization of glutathione S-transferases (GSTs) in Chenopodium quinoa. In this study, a total of 120 GST genes (CqGSTs) were identified and divided into 11 classes of which tau and phi were highest in numbers. The average protein length of protein was found to be 279.06 with their corresponding average molecular weight of 31,819.4 kDa. The subcellular localization analysis results showed that proteins were centrally localized in the cytoplasm followed by chloroplast, mitochondria and plastids. Structural analysis revealed the presence of 2 -14 exons in CqGST genes. Most of the proteins possessed two exon one intron organization. MEME analysis identified 15 significantly conserved motifs with a width of 6-50 amino acids. Motifs 1, 3, 2, 5, 6, 8, 9 and 13 were found specifically in tau class family; motifs 3, 4, 5, 6, 7 and 9 were found in phi class gene family, while motifs 3, 4, 13 and 14 were found in metaxin class. Multiple sequence alignment revealed highly conserved N-terminus with active site serine (Ser; S) or cysteine (Cys; C) residue for the activation of GSH binding and GST catalytic activity. The gene loci were found to be unevenly distributed across 18 different chromosomes with a maximum of 17 genes located on chromosome number 7. Dominance of alpha helix was followed by coil, extended strand and beta turns. Gene duplication analysis revealed that segmental duplication and purifying type selection were highest in number and found to be main source of expansion of GST gene family. Cis acting regulatory elements analysis showed the presence of 21 different elements involved in stress, hormone and light response and cellular development. The evolutionary relationship of CqGST proteins carried out using maximum likelihood method revealed that all the tau and phi class GSTs were closely associated with those of G. max, O. sativa and A. thaliana. Molecular docking of GST molecules with the fungicide metalaxyl showed that the CqGSTF1 had the lowest binding energy. The comprehensive study of CqGST gene family in quinoa provides groundwork for further functional analysis of CqGST genes in the species at molecular level and has potential applications in plant breeding.

Molecular Characterization of Dehydrin in Azraq Saltbush among Related Atriplex Species

Atriplex spp. (saltbush) is known to survive extremely harsh environmental stresses such as salinity and drought. It mitigates such conditions based on specialized physiological and biochemical characteristics. Dehydrin genes (DHNs) are considered major players in this adaptation. In this study, a novel DHN gene from Azrak (Jordan) saltbush was characterized along with other Atriplex species from diverse habitats. Intronless DHN-expressed sequence tags (495-761 bp) were successfully cloned and sequenced. Saltbush dehydrins contain one S-segment followed by three K-segments: an arrangement called SK3-type. Two substantial insertions were detected including three copies of the K2-segemnet in A. canescens. New motif variants other than the six-serine standard were evident in the S-segment. AhaDHN1 (A. halimus) has a cysteine residue (SSCSSS), while AgaDHN1 (A. gardneri var. utahensis) has an isoleucine residue (SISSSS). In contrast to the conserved K1-segment, both the K2- and K3-segment showed several substitutions, particularly in AnuDHN1 (A. nummularia). In addition, a parsimony phylogenetic tree based on homologs from related genera was constructed. The phylogenetic tree resolved DHNs for all of the investigated Atriplex species in a superclade with an 85% bootstrap value. Nonetheless, the DHN isolated from Azraq saltbush was uniquely subclustred with a related genera Halimione portulacoides. The characterized DHNs revealed tremendous diversification among the Atriplex species, which opens a new venue for their functional analysis.

The complete chloroplast genome of Pseudostellaria davidii (franch.) Pax, 1934

Pseudostellaria davidii (Franch.) Pax belongs to subseries distancs of Pseudostellaria (Caryophyllaceae), and is mainly distributed in north-eastern Asia. The complete chloroplast (cp) genome of P. davidii was assembled and annotated for the first time in this study. The cp genome of P. davidii is 149,732 bp in length with the GC content of 36.57%, and it consists of four subregions: a large single-copy (LSC) region of 81,156 bp, a small single-copy (SSC) region of 16,894 bp and two inverted repeats (IR) regions of 25,841 bp each. The cp genome of P. davidii encodes a total of 111 unique genes, which are 77 protein-coding genes, four rRNA genes, and 30 tRNA genes. The results of phylogenetic analysis strongly suggested that Pseudostellaria was a monophyletic group and P. davidii forms an independent sister clade to other species of Pseudostellaria.

Plastome variations reveal the distinct evolutionary scenarios of plastomes in the subfamily Cereoideae (Cactaceae)

BACKGROUND

The cactus family (Cactaceae) has been reported to have evolved a minimal photosynthetic plastome size, with the loss of inverted-repeat (IR) regions and NDH gene suites. However, there are very limited genomic data on the family, especially Cereoideae, the largest subfamily of cacti.

RESULTS

In the present study, we assembled and annotated 35 plastomes, 33 of which were representatives of Cereoideae, alongside 2 previously published plastomes. We analyzed the organelle genomes of 35 genera in the subfamily. These plastomes have variations rarely observed in those of other angiosperms, including size differences (with ~ 30 kb between the shortest and longest), dramatic dynamic changes in IR boundaries, frequent plastome inversions, and rearrangements. These results suggested that cacti have the most complex plastome evolution among angiosperms.

CONCLUSION

These results provide unique insight into the dynamic evolutionary history of Cereoideae plastomes and refine current knowledge of the relationships within the subfamily.

CgMYB1, an R2R3-MYB transcription factor, can alleviate abiotic stress in an annual halophyte Chenopodium glaucum

MYB transcription factors (TFs) are important regulators of the stress response in plants. In the present study, we characterized the CgMYB1 gene in Chenopodium glaucum, a member of the R2R3-MYB TF family. CgMYB1 was located in the nucleus with an activating domain at the C terminus. The CgMYB1 gene could be induced by salt and cold stress in C. glaucum. Overexpressing CgMYB1 in Arabidopsis significantly enhanced salt and cold tolerance, probably by improving physiological performance and stress-related gene expression. Further analysis suggests that the positive response of CgMYB1 to abiotic stress may partially be attributed to the interaction between CgMYB1 and the CgbHLH001 promoter followed by activation of downstream stress-responsive genes, which mediates stress tolerance. Our findings should contribute to further understanding of the function of R2R3 MYB TF in response to abiotic stress.

Promoter activity and transcriptome analyses decipher functions of CgbHLH001 gene (Chenopodium glaucum L.) in response to abiotic stress

BACKGROUND

Our previous studies revealed that CgbHLH001 transcription factor (TF) played an important role in abiotic stress tolerance, suggesting that its promoter was a potential target in response to stress signals. In addition, the regulatory mechanism of CgbHLH001 TF is still limited.

RESULTS

In the present study, a 1512 bp of 5'-flanking sequence of CgbHLH001 gene was identified, and the sequence carried quite a few of cis-acting elements. The gene promoter displayed strong activity and was induced by multiple abiotic stress. A series of 5'-deletions of the promoter sequence resulted in a gradual decrease in its activity, especially, the 5' untranslated region (UTR) was necessary to drive promoter activity. Further, CgbHLH001 promoter drove its own gene overexpression ectopically at the transcriptional and translational levels, which in turn conferred the stress tolerance to transgenic Arabidopsis. Transcriptome analysis showed that salt stress induced a large number of genes involved in multiple biological regulatory processes. Differentially expressed genes (DEGs) that mediate phytohormone signal transduction and mitogen-activated protein kinase (MAPK) signaling pathway were widely induced and mostly upregulated under salt stress, and the transcription levels in PbHLH::bHLH-overexpressing transgenic lines were higher than that of 35S::bHLH overexpression.

CONCLUSIONS

The CgbHLH001 promoter exhibited a positive response to abiotic stress and its 5' UTR sequence enhanced the regulation of gene expression to stress. A few important pathways and putative key genes involved in salt tolerance were identified, which can be used to elucidate the mechanism of salt tolerance and decipher the regulatory mechanism of promoters to develop an adaptation strategy for desert halophytes.

Genomic Sequence of Canadian Chenopodium berlandieri: A North American Wild Relative of Quinoa

Chenopodium berlandieri (pitseed goosefoot) is a widespread native North American plant, which was cultivated and consumed by indigenous peoples prior to the arrival of European colonists. Chenopodium berlandieri is closely related to, and freely hybridizes with the domesticated South American food crop C. quinoa. As such it is a potential source of wild germplasm for breeding with C. quinoa, for improved quinoa production in North America. The C. berlandieri genome sequence could also be a useful source of information for improving quinoa adaptation. To this end, we first optimized barcode markers in two chloroplast genes, rbcL and matK. Together these markers can distinguish C. berlandieri from the morphologically similar Eurasian invasive C. album (lamb's quarters). Second, we performed whole genome sequencing and preliminary assembly of a C. berlandieri accession collected in Manitoba, Canada. Our assembly, while fragmented, is consistent with the expected allotetraploid structure containing diploid Chenopodium sub-genomes A and B. The genome of our accession is highly homozygous, with only one variant site per 3-4000 bases in non-repetitive sequences. This is consistent with predominant self-fertilization. As previously reported for the genome of a partly domesticated Mexican accession of C. berlandieri, our genome assembly is similar to that of C. quinoa. Somewhat unexpectedly, the genome of our accession had almost as many variant sites when compared to the Mexican C. berlandieri, as compared to C. quinoa. Despite the overall similarity of our genome sequence to that of C. quinoa, there are differences in genes known to be involved in the domestication or genetics of other food crops. In one example, our genome assembly appears to lack one functional copy of the SOS1 (salt overly sensitive 1) gene. SOS1 is involved in soil salinity tolerance, and by extension may be relevant to the adaptation of C. berlandieri to the wet climate of the Canadian region where it was collected. Our genome assembly will be a useful tool for the improved cultivation of quinoa in North America.

Genome-wide identification and expression pattern analysis of quinoa BBX family

BBX is a transcription factor encoding zinc finger protein that plays a key role in plant growth and development as well as in responding to abiotic stresses. However, in quinoa, which is known as a "super grain" and has extremely high nutritional value, this gene family has not yet been thoroughly studied. In this study, in order to fully understand the family function of the BBX in quinoa, a total of 31 BBX members were identified by bioinformatics methods. These BBX members were mainly acidic proteins, and most of their secondary structures were random coil s, 31 CqBBX members were unevenly distributed on 17 chromosomes, and the analysis of replication events found that quinoa BBX genes produced a total of 14 pairs of gene replication. The BBX genes were divided into five subfamilies according to phylogenetics, and its gene structure and conserved motif were basically consistent with the classification of its phylogenetic tree. In addition, a total of 43 light response elements, hormone response elements, tissue-specific expression response elements, and abiotic stress response elements were found in the promoter region, involving stress elements such as drought and low temperature. Finally, the expression patterns of CqBBX genes in different tissues and abiotic stresses were studied by combining transcriptome data and qRT-PCR , and all 13 genes responded to drought, salt, and low-temperature stress to varying degrees. This study is the first comprehensive study of the BBX family of quinoa, and its results provide important clues for further analysis of the function of the abiotic stress response.

The complete chloroplast genome sequences of three lilies: genome structure, comparative genomic and phylogenetic analyses

We sequenced and analyzed the complete chloroplast genomes of Lilium amoenum, Lilium souliei, and Nomocharis forrestii in detail, including the first sequence and structural comparison of Nomocharis forrestii. We found that the lengths and nucleotide composition of the three chloroplast genes showed little variation. The chloroplast genomes of the three Lilium species contain 87 protein coding genes (PCGs), 38 tRNAs, and 8 rRNA genes. The only difference is that Nomocharis forrestii had an additional infA pseudogene. In the sequence analysis of the Lilium chloroplast genomes, 216 SSRs, 143 pairs of long repeats, 571 SNPs, and 202 indels were detected. In addition, we identified seven hypervariable regions that can be used as potential molecular markers and DNA barcodes of Lilium through complete sequence alignment. The phylogenetic tree was constructed from the three chloroplast genome sequences of Lilium obtained here and 40 chloroplast genome sequences from the NCBI database (including 35 Lilium species, 4 Fritillaria species, and one species of Smilax). The analysis showed that the species clustering of the genus Lilium essentially conformed to the classical morphological classification system of Comber, but differences in the classification of individual species remained. In our report, we support the reclassification of Lilium henryi and Lilium rosthorniiy in the genus Lilium. In general, this study not only provides genome data for three Lilium species, but also provides a comparative analysis of the Lilium chloroplast genomes. These advances will help to identify Lilium species, clarify the phylogenetic analysis of the Lilium genus, and help to solve and improve the disputes and deficiencies in the traditional morphological classification.

Complete chloroplast genome sequences of Phlomis fruticosa and Phlomoides strigosa and comparative analysis of the genus Phlomis sensu lato (Lamiaceae)

The taxonomic terms "Phlomis" and "Phlomoides" had been used to describe two sections within the genus Phlomis belonging to the Lamiaceae family. Recently, phylogenetic analyses using molecular markers showed that Phlomis and Phlomoides formed two monophyletic clades, and thus they are generally accepted as separate genera. In this study, we assembled the complete chloroplast genome of Phlomis fruticosa, which is the first reported chloroplast genome belonging to Phlomis genus, as well as the complete chloroplast genome of Phlomoides strigosa belonging to Phlomoides genus. The results showed that the length of chloroplast genome was 151,639 bp (Phlomis fruticosa) and 152,432 bp (Phlomoides strigosa), with conserved large single copy regions, small single copy regions, and inverted repeat regions. 121 genes in Phlomis fruticosa and 120 genes in Phlomoides strigosa were annotated. The chloroplast genomes of Phlomis fruticosa, Phlomoides strigosa, and three reported Phlomoides species, as well as those of 51 species from the Lamiaceae family, which covered 12 subfamilies, were subjected to phylogenetic analyses. The Phlomis and Phlomoides species were split into two groups, which were well supported by both maximum likelihood and Bayesian inference tree analyses. Our study provided further evidence to recognize Phlomis and Phlomoides as independent genera.

Complete chloroplast genome sequencing and comparative analysis of threatened dragon trees Dracaena serrulata and Dracaena cinnabari

Dracaena (Asparagaceae family) tree is famous for producing "dragon blood"—a bioactive red-colored resin. Despite its long history of use in traditional medicine, little knowledge exists on the genomic architecture, phylogenetic position, or evolution. Hence, in this study, we sequenced the whole chloroplast (cp) genomes of D. serrulata and D. cinnabari and performed comparative genomics of nine genomes of the genus Dracaena. The results showed that the genome sizes range from 155,055 (D. elliptica) to 155,449 (D. cochinchinensis). The cp genomes of D. serrulata and D. cinnabari encode 131 genes, each including 85 and 84 protein-coding genes, respectively. However, the D. hokouensis had the highest number of genes (133), with 85 protein coding genes. Similarly, about 80 and 82 repeats were identified in the cp genomes of D. serrulata and D. cinnabari, respectively, while the highest repeats (103) were detected in the cp genome of D. terniflora. The number of simple sequence repeats (SSRs) was 176 and 159 in D. serrulata and D. cinnabari cp genomes, respectively. Furthermore, the comparative analysis of complete cp genomes revealed high sequence similarity. However, some sequence divergences were observed in accD, matK, rpl16, rpoC2, and ycf1 genes and some intergenic spacers. The phylogenomic analysis revealed that D. serrulata and D. cinnabari form a monophyletic clade, sister to the remaining Dracaena species sampled in this study, with high bootstrap values. In conclusion, this study provides valuable genetic information for studying the evolutionary relationships and population genetics of Dracaena, which is threatened in its conservation status.

Trends and Limits for Quinoa Production and Promotion in Pakistan

Quinoa is known as a super food due to its extraordinary nutritional qualities and has the potential to ensure future global food and nutritional security. As a model plant with halophytic behavior, quinoa has potential to meet the challenges of climate change and salinization due to its capabilities for survival in harsh climatic conditions. The quinoa crop has received worldwide attention due to its adoption and production expanded in countries out of the native Andean region. Quinoa was introduced to Pakistan in 2009 and it is still a new crop in Pakistan. The first quinoa variety was registered in 2019, then afterward, its cultivation started on a larger scale. Weed pressure, terminal heat stress, stem lodging, bold grain size, and an unstructured market are the major challenges in the production and promotion of the crop. The potential of superior features of quinoa has not been fully explored and utilized. Hence, there is a need to acquire more diverse quinoa germplasm and to establish a strong breeding program to develop new lines with higher productivity and improved crop features for the Pakistan market. Mechanized production, processing practices, and a structured market are needed for further scaling of quinoa production in Pakistan. To achieve these objectives, there is a dire need to create an enabling environment for quinoa production and promotion through the involvement of policymakers, research institutions, farmers associations, and the private sector.

A dPCR Method for Quantitative Authentication of Wild Lingonberry (Vaccinium vitis-idaea) versus Cultivated American Cranberry (V. macrocarpon)

Berries of the genus Vaccinium are highly valued health-beneficial superfoods, which are commonly subjected to adulteration and mixed with each other, or with other common berry species. A quantitative DNA-based method utilizing a chip-based digital polymerase chain reaction (dPCR) technique was developed for identifying and quantifying wild lingonberry (V. vitis-idaea) and cultivated American cranberry (V. macrocarpon). The dPCR method with species-specific primers for mini-barcoding was designed based on the indel regions found in the trnI-CAU-trnL-CAA locus in the chloroplast genome. The designed primers were able to amplify only target species, enabling to distinguish the two closely related species with good sensitivity. Our results illustrated the ability of the method to identify lingonberry and American cranberry DNA using PCR without the need for probes or further sequencing. The dPCR method could also quantify the DNA copy number in mixed samples. Based on this study, the method provides a basis for a simple, fast, and sensitive quantitative authentication analysis of lingonberry and American cranberry by dPCR. Moreover, it can also provide a platform for authentication analyses of other plant species as well by utilizing the indel regions of chloroplast genomes.

Characterisation of the Complete Chloroplast Genomes of Seven Hyacinthus orientalis L. Cultivars: Insights into Cultivar Phylogeny

To improve agricultural performance and obtain potential economic benefits, an understanding of phylogenetic relationships of Hyacinthus cultivars is needed. This study aims to revisit the phylogenetic relationships of Hyacinthus cultivars using complete chloroplast genomes. Nine chloroplast genomes were de novo sequenced, assembled and annotated from seven cultivars of Hyacinthus orientalis and two Scilloideae species including Bellevalia paradoxa and Scilla siberica. The chloroplast genomes of Hyacinthus cultivars ranged from 154,458 bp to 154,641 bp, while those of Bellevalia paradoxa and Scilla siberica were 154,020 bp and 154,943 bp, respectively. Each chloroplast genome was annotated with 133 genes, including 87 protein-coding genes, 38 transfer RNA genes and 8 ribosomal RNA genes. Simple sequence repeats AAGC/CTTG and ACTAT/AGTAT were identified only in ‘Eros’, while AAATC/ATTTG were identified in all cultivars except ‘Eros’. Five haplotypes were identified based on 460 variable sites. Combined with six other previously published chloroplast genomes of Scilloideae, a sliding window analysis and a phylogenetic analysis were performed. Divergence hotspots ndhA and trnG-UGC were identified with a nucleotide diversity threshold at 0.04. The phylogenetic positions of Hyacinthus cultivars were different from the previous study using ISSR. Complete chloroplast genomes serve as new evidence in Hyacinthus cultivar phylogeny, contributing to cultivar identification, preservation and breeding.

The complete chloroplast genome of Atriplex gmelinii C. A. Mey. ex Bong. (Amaranthaceae)

Atriplex gmelinii C. A. Mey. Ex Bong._1838 is an annual halophytic herb found in East Asia and North America. The chloroplast genome of A. gmelinii was successfully sequenced. The assembled genome (151,852 bp; GC ratio, 37.3%) is composed of four subregions, with the large single copy (LSC; 83,872 bp; 35.4%) and small single copy (SSC; 17,812 bp; 30.9%) regions separated by two regions of inverted repeat regions (25,084 bp; 42.8%). A total of 130 genes were predicted with 85 protein-coding genes, 8 rRNAs, and 37 tRNAs. The phylogenetic analyses inferred from whole chloroplast genomes of 35 species, including 34 species in Amaranthaceae and one outgroup species, suggest a close relationship between A. gmelinii and A. centralasiatica.

Characterization of the complete chloroplast genome of Carya laciniosa (F.Michx.) G.Don (Juglandaceae)

Carya laciniosa (F.Michx.) G.Don is a slow-growing valuable tree in the family Juglandaceae. In this research, we first reported the chloroplast genome sequence of C. laciniosa. The chloroplast genome was a circular form with 160,832 bp in size, comprising four subregions: a large single-copy (LSC) region of 90,065 bp, a small single-copy (SSC) region of 18,791 bp, and two copies of inverted repeats (IRs; IRa and IRb) of 25,988 bp each. A total of 132 genes were annotated in the chloroplast genome of C. laciniosa, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The overall GC content of the whole genome was 36.14%. Phylogenetic analysis suggested that C. laciniosa was closely related to C. ovata.

Comparative chloroplast genome analyses of cultivated spinach and two wild progenitors shed light on the phylogenetic relationships and variation

Spinacia is a genus of important leafy vegetable crops worldwide and includes cultivated Spinacia oleracea and two wild progenitors, Spinacia turkestanica and Spinacia tetrandra. However, the chloroplast genomes of the two wild progenitors remain unpublished, limiting our knowledge of chloroplast genome evolution among these three Spinacia species. Here, we reported the complete chloroplast genomes of S. oleracea, S. turkestanica, and S. tetrandra obtained via Illumina sequencing. The three chloroplast genomes exhibited a typical quadripartite structure and were 150,739, 150,747, and 150,680 bp in size, respectively. Only three variants were identified between S. oleracea and S. turkestanica, whereas 690 variants were obtained between S. oleracea and S. tetrandra, strongly demonstrating the close relationship between S. turkestanica and S. oleracea. This was further supported by phylogenetic analysis. We reported a comprehensive variant dataset including 503 SNPs and 83 Indels using 85 Spinacia accessions containing 61 S. oleracea, 16 S. turkestanica, and eight S. tetrandra accessions. Thirteen S. oleracea accessions were derived through introgression from S. turkestanica that acts as the maternal parent. Together, these results provide a valuable resource for spinach breeding programs and improve our understanding of the phylogenetic relationships within Amaranthaceae.

New Insight into the Phylogeny and Taxonomy of Cultivated and Related Species of Crataegus in China, Based on Complete Chloroplast Genome Sequencing

Hawthorns (Crataegus L.) are one of the most important processing and table fruits in China, due to their medicinal properties and health benefits. However, the interspecific relationships and evolution history of cultivated Crataegus in China remain unclear. Our previously published data showed C. bretschneideri may be derived from the hybridization of C. pinnatifida with C. maximowiczii, and that introgression occurs between C. hupehensis, C. pinnatifida, and C. pinnatifida var. major. In the present study, chloroplast sequences were used to further elucidate the phylogenetic relationships of cultivated Crataegus native to China. The chloroplast genomes of three cultivated species and one related species of Crataegus were sequenced for comparative and phylogenetic analyses. The four chloroplast genomes of Crataegus exhibited typical quadripartite structures and ranged from 159,607 bp (C. bretschneideri) to 159,875 bp (C. maximowiczii) in length. The plastomes of the four species contained 113 genes consisting of 79 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Six hypervariable regions (ndhC-trnV(UAC)-trnM(CAU), ndhA, atpH-atpI, ndhF, trnR(UCU)-atpA, and ndhF-rpl32), 196 repeats, and a total of 386 simple sequence repeats were detected as potential variability makers for species identification and population genetic studies. In the phylogenomic analyses, we also compared the entire chloroplast genomes of three published Crataegus species: C. hupehensis (MW201730.1), C. pinnatifida (MN102356.1), and C. marshallii (MK920293.1). Our phylogenetic analyses grouped the seven Crataegus taxa into two main clusters. One cluster included C. bretschneideri, C. maximowiczii, and C. marshallii, whereas the other included C. hupehensis, C. pinnatifida, and C. pinnatifida var. major. Taken together, our findings indicate that C. maximowiczii is the maternal origin of C. bretschneideri. This work provides further evidence of introgression between C. hupehensis, C. pinnatifida, and C. pinnatifida var. major, and suggests that C. pinnatifida var. major might have been artificially selected and domesticated from hybrid populations, rather than evolved from C. pinnatifida.

Characterization of the complete chloroplast genome of an annual herb, Chenopodium album (Amaranthaceae)

Chenopodium album is an annual herb from Amaranthaceae with worldwide distribution. It is a leafy vegetable as well as an important subsidiary grain crop with high nutritional value and medicinal value. In this study, we reported the complete chloroplast genome of C. album. The total chloroplast genome was 152,167 bp in length, containing a large single-copy region (LSC, 83,676 bp), a small single-copy region (SSC, 18,105 bp), and a pair of inverted repeat regions (IRs, 25,193 bp). The complete chloroplast genome contains 110 genes, including 78 protein-coding genes, 28 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes with an overall GC content of 37.3%. Phylogenetic analysis showed that C. album was sister to C. acuminatum within Chenopodioideae. The complete chloroplast genome of C. album will provide useful resources for the development and utilization of this species and the phylogenetic study of Amaranthaceae.

Plastid Phylogenomic Data Offers Novel Insights Into the Taxonomic Status of the Trichosanthes kirilowii Complex (Cucurbitaceae) in South Korea

Trichosanthes is a genus in Cucurbitaceae comprising 90-100 species. Trichosanthes species are valuable as herbaceous medicinal ingredients. The fruits, seeds, and roots of species such as T. kirilowii and T. rosthornii are used in Korean traditional herbal medicines. T. rosthornii is only found in China, whereas in South Korea two varieties, T. kirilowii var. kirilowii and T. kirilowii var. japonica, are distributed. T. kirilowii var. kirilowii and T. kirilowii var. japonica have different fruit and leaf shapes but are recognized as belonging to the same species. Furthermore, although its members have herbal medicine applications, genomic information of the genus is still limited. The broad goals of this study were (i) to evaluate the taxonomy of Trichosanthes using plastid phylogenomic data and (ii) provide molecular markers specific for T. kirilowii var. kirilowii and T. kirilowii var. japonica, as these have differences in their pharmacological effectiveness and thus should not be confused and adulterated. Comparison of five Trichosanthes plastid genomes revealed locally divergent regions, mainly within intergenic spacer regions (trnT-UGU-trnL-UAA: marker name Tri, rrn4.5-rrn5: TRr, trnE-UUC-trnT-GGU: TRtt). Using these three markers as DNA-barcodes for important herbal medicine species in Trichosanthes, the identity of Trichosanthes material in commercial medicinal products in South Korea could be successfully determined. Phylogenetic analysis of the five Trichosanthes species revealed that the species are clustered within tribe Sicyoeae. T. kirilowii var. kirilowii and T. rosthornii formed a clade with T. kirilowii var. japonica as their sister group. As T. kirilowii in its current circumscription is paraphyletic and as the two varieties can be readily distinguished morphologically (e.g., in leaf shape), T. kirilowii var. japonica should be treated (again) as an independent species, T. japonica.

Genome-Wide Identification and Expression Analysis of the β-Amylase Gene Family in Chenopodium quinoa

β-Amylase (BAM) is an important starch hydrolase, playing a role in a variety of plant growth and development processes. In this study, 22 BAM gene family members (GFMs) were identified in quinoa (Chenopodium quinoa), an ancient crop gaining modern consumer acceptance because of its nutritional qualities. The genetic structure, phylogenetic and evolutionary relationships, and expression patterns of CqBAM GFMs in different tissues, were analyzed. Phylogenetic analyses assigned the CqBAMs, AtBAMs, and OsBAMs into four clades. The CqBAM gene family had expanded due to segmental duplication. RNA-seq analysis revealed expression of the duplicated pairs to be similar, with the expression of CqBAM GFM pairs showing a degree of tissue specificity that was confirmed by reverse transcription quantitative PCR (RT-qPCR). Several CqBAM GFMs were also responsive to abiotic stresses in shoots and/or roots. In conclusion, the BAM gene family in quinoa was identified and systematically analyzed using bioinformatics and experimental methods. These results will help to elucidate the evolutionary relationship and biological functions of the BAM gene family in quinoa.

Characterization of the Complete Chloroplast Genome of Buddleja Lindleyana

BACKGROUND

Buddleja lindleyana Fort., which belongs to the Loganiaceae with a distribution throughout the tropics, is widely used as an ornamental plant in China. Buddleja contains several morphologically similar species, which need to be identified by molecular identification. But there is little molecular research on the genus Buddleja.

OBJECTIVE

Using molecular biology techniques to sequence and analyze the complete chloroplast (cp) genome of B. lindleyana.

METHODS

According to next-generation sequencing to sequence the genome data, a series of bioinformatics software were used to assembly and analysis the molecular structure of cp genome of B. lindleyana.

RESULTS

The complete cp genome of B. lindleyana is a circular 154,487-bp-long molecule with a GC content of 38.1%. It has a familiar quadripartite structure, including a large single-copy region (LSC; 85,489 bp), a small single-copy region (SSC; 17,898bp) and a pair of inverted repeats (IRs; 25,550 bp). A total of 133 genes were identified in the genome, including 86 protein-coding genes, 37 tRNA genes, 8 rRNA genes and 2 pseudogenes.

CONCLUSIONS

These results suggested that B. lindelyana cp genome could be used as a potential genomic resource to resolve the phylogenetic positions and relationships of Loganiaceae, and will offer valuable information for future research in the identification of Buddleja species and will conduce to genomic investigations of these species.

The Comparative Analyses of Six Complete Chloroplast Genomes of Morphologically Diverse Chenopodium album L. (Amaranthaceae) Collected in Korea

Chenopodium album sensu stricto belonging to C. album aggregate is an annual cosmopolitan weed displaying the diversity of morphologies. We completed the six chloroplast genomes of C. album s. str. collected in Korea to understand the relationship between the diversity of chloroplast genomes and their morphological variations. All six C. album chloroplast genomes have a typical quadripartite structure with length ranging from 151,906 bp to 152,199 bp, similar to the previously sequenced C. album chloroplast genome (NC_034950). In total, 56 single nucleotide polymorphisms (SNPs) and 26 insertion and deletion (INDEL) regions (308 bp in total) were identified from the six chloroplast genomes, presenting a low level of intraspecific variations in comparison to the other angiosperm species. 376 normal simple sequence repeats were identified in all seven C. album chloroplast genomes. The phylogenetic analysis based on all available complete Amaranthaceae chloroplast genomes presents phylogenetic positions of six C. album samples as well as correlation with one of C. album morphological features. Our results provide the way to investigate intraspecific features of C. album chloroplast genomes and also the insights of understanding various intraspecific characteristics including morphological features.

Chloroplast genome sequence of Chongming lima bean (Phaseolus lunatus L.) and comparative analyses with other legume chloroplast genomes

Background

Lima bean (Phaseolus lunatus L.) is a member of subfamily Phaseolinae belonging to the family Leguminosae and an important source of plant proteins for the human diet. As we all know, lima beans have important economic value and great diversity. However, our knowledge of the chloroplast genome level of lima beans is limited.

Results

The chloroplast genome of lima bean was obtained by Illumina sequencing technology for the first time. The Cp genome with a length of 150,902 bp, including a pair of inverted repeats (IRA and IRB 26543 bp each), a large single-copy (LSC 80218 bp) and a small single-copy region (SSC 17598 bp). In total, 124 unique genes including 82 protein-coding genes, 34 tRNA genes, and 8 rRNA genes were identified in the P. lunatus Cp genome. A total of 61 long repeats and 290 SSRs were detected in the lima bean Cp genome. It has a typical 50 kb inversion of the Leguminosae family and an 70 kb inversion to subtribe Phaseolinae. rpl16, accD, petB, rsp16, clpP, ndhA, ndhF and ycf1 genes in coding regions was found significant variation, the intergenic regions of trnk-rbcL, rbcL-atpB, ndhJ-rps4, psbD-rpoB, atpI-atpA, atpA-accD, accD-psbJ, psbE-psbB, rsp11-rsp19, ndhF-ccsA was found in a high degree of divergence. A phylogenetic analysis showed that P. lunatus appears to be more closely related to P. vulgaris, V.unguiculata and V. radiata.

Conclusions

The characteristics of the lima bean Cp genome was identified for the first time, these results will provide useful insights for species identification, evolutionary studies and molecular biology research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07467-8.

The complete chloroplast genome of Chenopodium acuminatum Willd. (Amaranthaceae)

The complete chloroplast genome (plastome) of the Chenopodium acuminatum was assembled and annotated in this study. The complete plastome was composed of circular DNA molecules with a total length of 152,200 bp, comprising a large single-copy region (83,683 bp), a small single-copy region (18,131 bp), and two inverted repeat regions (25,193 bp). GC content of this complete plastome was 37.2%. In total, 113 unique genes were annotated, including 79 protein-coding genes (PCGs), 30 transfer RNAs, and 4 ribosomal RNAs. Phylogenomic analysis showed that C. acuminatum was closely related to C. album.

Characterization of the complete plastome of Cyperus rotundus L. (Cyperaceae)

Cyperus rotundus L. (C. rotundus) is a sedge belonging to the family Cyperaceae and is widely distributed in tropical and warmer temperate regions worldwide. It is one of the oldest traditional medicinal herbs in China, India, Japan, and Korea. In this study, we sequenced the complete chloroplast genome of C. rotundus on the Illumina HiSeq Platform. The chloroplast genome is 182,986 bp in length, with a typical quadripartite structure and consisting of a pair of inverted repeat (IR) regions (35,969 bp) separated by a large single-copy (LSC) region (100,733 bp) and a small single-copy (SSC) region (10,315 bp). It was predicted to contain a total of 133 genes, with an overall GC content of 33.26%. Phylogenetic analysis suggested C. rotundus is sister to Eleocharis celluosa and Eleocharis dulcis.

Complete chloroplast genome of the grain Chenopodium quinoa Willd., an important economical and dietary plant

The grain Chenopodium quinoa Willd. is the main traditional food of Inca aboriginal, which was a native grain in South American Andes Mountains, the edible and cultivation history of which has been more than five thousand years. In this study, we sequenced the complete chloroplast genome of C. quinoa on the Illumina platform by shotgun genome skimming method. The complete chloroplast genome of C. quinoa was 152,087 bp in length with the GC content 37.2%, which was comprised of a large single copy (LSC) region of 83,570 bp, a small single copy (SSC) region of 18,107 bp, and a pair of inverted repeats (IRA/IRB) of 25,205 bp. The chloroplast genome encoded 133 genes including 88 protein-coding genes, 37 tRNA genes and eight rRNA genes. Phylogenetic analysis constructed using the maximum likelihood (ML) method strongly supported the monophyly of each genus in the family Chenopodiaceae, and the genus Chenopodium is sister to Spinacia as a cluster, which closely grouped to Dysphania.

Dissection for Floral Micromorphology and Plastid Genome of Valuable Medicinal Borages Arnebia and Lithospermum (Boraginaceae)

The genera Arnebia and Lithospermum (Lithospermeae-Boraginaceae) comprise 25-30 and 50-60 species, respectively. Some of them are economically valuable, as their roots frequently contain a purple-red dye used in the cosmetic industry. Furthermore, dried roots of Arnebia euchroma, A. guttata, and Lithospermum erythrorhizon, which have been designated Lithospermi Radix, are used as traditional Korean herbal medicine. This study is the first report on the floral micromorphology and complete chloroplast (cp) genome sequences of A. guttata (including A. tibetana), A. euchroma, and L. erythrorhizon. We reveal great diversity in floral epidermal cell patterns, gynoecium, and structure of trichomes. The cp genomes were 149,361-150,465 bp in length, with conserved quadripartite structures. In total, 112 genes were identified, including 78 protein-coding regions, 30 tRNA genes, and four rRNA genes. Gene order, content, and orientation were highly conserved and were consistent with the general structure of angiosperm cp genomes. Comparison of the four cp genomes revealed locally divergent regions, mainly within intergenic spacer regions (atpH-atpI, petN-psbM, rbcL-psaI, ycf4-cemA, ndhF-rpl32, and ndhC-trnV-UAC). To facilitate species identification, we developed molecular markers psaA- ycf3 (PSY), trnI-CAU- ycf2 (TCY), and ndhC-trnV-UAC (NCTV) based on divergence hotspots. High-resolution phylogenetic analysis revealed clear clustering and a close relationship of Arnebia to its Lithospermum sister group, which was supported by strong bootstrap values and posterior probabilities. Overall, gynoecium characteristics and genetic distance of cp genomes suggest that A. tibetana, might be recognized as an independent species rather than a synonym of A. guttata. The present morphological and cp genomic results provide useful information for future studies, such as taxonomic, phylogenetic, and evolutionary analysis of Boraginaceae.

The complete plastome sequences of five Aponogeton species (Aponogetonaceae): Insights into the structural organization and mutational hotspots

Members of the aquatic plant genus Aponogeton are widely used commercially in aquariums because of their variable leaf shape and unique inflorescences. However, due to extensive similarity between species in this genus, morphological characters are generally inadequate for taxonomic classification. Currently, molecular makers available for taxonomic and phylogenetic studies of Aponogeton are limited. One approach to clarifying relationships between species in these complex groups is to use divergence hotspot regions within the genome. Here, we sequenced and analyzed the plastomes of five Aponogeton species collected from China, Zambia, and Kenya, and subsequently screened these plastomes for divergent DNA hotspots. The five plastomes are circular structures with sizes ranging from 154,167 bp to 154,860 bp. The Large and the Small Single Copies are separated by two Inverted Repeats. One hundred and thirteen unique genes were identified including 79 protein-coding, 30 tRNA, and four rRNA genes. We found that the most abundant repeats in all but one species were mononucleotide repeats (A/T) and that there were 23 potential RNA ending sites. Interestingly, a ~3 kb inversion, which includes the accD gene, was detected within the Asian species of Aponogeton. The inversion may be related to more frequent exchanges between this region and the nuclear genome. Furthermore, we detected mutational hotspot sites among the five Aponogeton species. Three of these hotspots are intergenic spacer regions (accD-psaI, rbcL-accD and trnH-GUG-psbA) that might be suitable for use as barcodes to resolve intra-generic relationships. We also identified four highly variable protein-coding genes (ccsA, rpl22, rps16 and ycf1) may be used as barcodes to resolve the higher-level phylogenies. Our study will provide valuable molecular resources for the taxonomic and phylogenomic study of the complex genus Aponogeton.

Adaptation Evolution and Phylogenetic Analyses of Species in Chinese Allium Section Pallasia and Related Species Based on Complete Chloroplast Genome Sequences

The section Pallasia is one of the components of the genus Allium subgenus Allium (Amaryllidaceae), and species relationship in this section is still not resolved very well, which hinders further evolutionary and adaptive studies. Here, the complete chloroplast genomes of five sect. Pallasia species were reported, and a comparative analysis was performed with other three related Allium species. The genome size of the eight species ranged from 151,672 bp to 153,339 bp in length, GC content changed from 36.7% to 36.8%, and 130 genes (except Allium pallasii), 37 tRNA, and 8 rRNA were identified in each genome. By analyzing the IR/LSC and IR/SSC boundary, A. pallasii exhibited differences compared with other seven species. Phylogenetic analysis achieved high supports in each branch, seven of the eight Allium species cluster into a group, and A. pallasii exhibit a close relationship with A. obliquum. Higher pairwise Ka/Ks ratios were found in A. schoenoprasoides compared to A. caeruleum and A. macrostemon while a lower value of Ka/Ks ratios was detected between A. caeruleum and A. macrostemon. This study will be a great contribution to the future phylogenetic and adaptive research in Allium.

Chloroplast genomic comparison of two sister species Allium macranthum and A. fasciculatum provides valuable insights into adaptive evolution

BACKGROUND

Allium macranthum and Allium fasciculatum are two sister species and their natural populations are separated by high mountains and deep valleys with exact opposite habitat. The chloroplast genome in angiosperms has showed useful for investigating plant evolution and systematic studies.

OBJECTIVE

Comparative analysis of these genomes revealed potential markers and phylogenetic analysis, and discuss the influence of positive selected sites on adaptive evolution.

METHODS

Here, we sequenced the complete chloroplast genomes of these two species and analyzed the repeat sequences components, nucleotide diversity, selection pressure and the phylogeny relationships with related species.

RESULTS

A typical quadripartite structure was detected with a genome size changed from 152,148 to 152,931 bp. We identified 67 and 79 simple sequence repeats in A. macranthum and A. fasciculatum, in which the mono-nucleotide repeats A/T possess the highest percentage. Three mutational hotspots (rpl32, rps16 and matK) at the SSC and LSC regions were observed, which showed remarkably higher Pi value (> 0.03). Additionally, eight genes (rpoA, atpF, cemA, rps4, ccsA, rpoC2, rpl14 and clpP) exhibited elevated pairwise Ka/Ks ratios in alpine species. Phylogenetic analyses based on the CDS sequences and the whole complete genomes showed same topologies with high support, and A. macranthum was closely clustered with A. fasciculatum within the fourteen Amaryllidaceae species.

CONCLUSION

Their coding proteins of these genes often functioned in chloroplast protein synthesis, gene transcription, energy transformation and regulation and photosynthesis. These results provide valuable insights into the alpine species adaptation and evolution.

Comparative Analysis of Actaea Chloroplast Genomes and Molecular Marker Development for the Identification of Authentic Cimicifugae Rhizoma

Actaea (Ranunculaceae; syn. Cimicifuga) is a controversial and complex genus. Dried rhizomes of Actaea species are used as Korean traditional herbal medicine. Although Actaea species are valuable, given their taxonomic classification and medicinal properties, sequence information of Actaea species is limited. In this study, we determined the complete chloroplast (cp) genome sequences of three Actaea species, including A. simplex, A. dahurica, and A. biternata. The cp genomes of these species varied in length from 159,523 to 159,789 bp and contained 112 unique functional genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. Gene order, orientation, and content were well conserved in the three cp genomes. Comparative sequence analysis revealed the presence of hotspots, including ndhC-trnV-UAC, in Actaea cp genomes. High-resolution phylogenetic relationships were established among Actaea species based on cp genome sequences. Actaea species were clustered into each Actaea section, consistent with the Angiosperm Phylogeny Group (APG) IV system of classification. We also developed a novel indel marker, based on copy number variation of tandem repeats, to facilitate the authentication of the herbal medicine Cimicifugae Rhizoma. The availability Actaea cp genomes will provide abundant information for the taxonomic and phylogenetic analyses of Actaea species, and the Actaea (ACT) indel marker will be useful for the authentication of the herbal medicine.

Comparative Analysis of the Complete Chloroplast Genomes in Allium Subgenus Cyathophora (Amaryllidaceae): Phylogenetic Relationship and Adaptive Evolution

Recent advances in molecular phylogenetics provide us with information of Allium L. taxonomy and evolution, such as the subgenus Cyathophora, which is monophyletic and contains five species. However, previous studies detected distinct incongruence between the nrDNA and cpDNA phylogenies, and the interspecies relationships of this subgenus need to be furtherly resolved. In our study, we newly assembled the whole chloroplast genome of four species in subgenus Cyathophora and two allied Allium species. The complete cp genomes were found to possess a quadripartite structure, and the genome size ranged from 152,913 to 154,174 bp. Among these cp genomes, there were subtle differences in the gene order, gene content, and GC content. Seven hotspot regions (infA, rps16, rps15, ndhF, trnG-UCC, trnC-GCA, and trnK-UUU) with nucleotide diversity greater than 0.02 were discovered. The selection analysis showed that some genes have elevated Ka/Ks ratios. Phylogenetic analysis depended on the complete chloroplast genome (CCG), and the intergenic spacer regions (IGS) and coding DNA sequences (CDS) showed same topologies with high support, which revealed that subgenus Cyathophora was a monophyletic group, containing four species, and A. cyathophorum var. farreri was sister to A. spicatum with 100% bootstrap value. Our study revealed selective pressure may exert effect on several genes of the six Allium species, which may be useful for them to adapt to their specific living environment. We have well resolved the phylogenetic relationship of species in the subgenus Cyathophora, which will contribute to future evolutionary studies or phylogeographic analysis of Allium.

The complete chloroplast genome of Microcycas calocoma (Miq.) A. DC. (Zamiaceae, Cycadales) and evolution in Cycadales

Cycadales is an extant group of seed plants occurring in subtropical and tropical regions comprising putatively three families and 10 genera. At least one complete plastid genome sequence has been reported for all of the 10 genera except Microcycas, making it an ideal plant group to conduct comprehensive plastome comparisons at the genus level. This article reports for the first time the plastid genome of Microcycas calocoma. The plastid genome has a length of 165,688 bp with 134 annotated genes including 86 protein-coding genes, 47 non-coding RNA genes (39 tRNA and eight rRNA) and one pseudogene. Using global sequence variation analysis, the results showed that all cycad genomes share highly similar genomic profiles indicating significant slow evolution and little variation. However, identity matrices coinciding with the inverted repeat regions showed fewer similarities indicating that higher polymorphic events occur at those sites. Conserved non-coding regions also appear to be more divergent whereas variations in the exons were less discernible indicating that the latter comprises more conserved sequences. Phylogenetic analysis using 81 concatenated protein-coding genes of chloroplast (cp) genomes, obtained using maximum likelihood and Bayesian inference with high support values (>70% ML and = 1.0 BPP), confirms that Microcycas is closest to Zamia and forms a monophyletic clade with Ceratozamia and Stangeria. While Stangeria joined the Neotropical cycads Ceratozamia, Zamia and Microcyas, Bowenia grouped with the Southern Hemisphere cycads Encephalartos, Lepidozamia and Macrozamia. All Cycas species formed a distinct clade separated from the other genera. Dioon, on the other hand, was outlying from the rest of Zamiaceae encompassing two major clades—the Southern Hemisphere cycads and the Neotropical cycads. Analysis of the whole cp genomes in phylogeny also supports that the previously recognized family—Stangeriaceae—which contained Bowenia and Stangeria, is not monophyletic. Thus, the cp genome topology obtained in our study is congruent with other molecular phylogenies recognizing only a two-family classification (Cycadaceae and Zamiaceae) within extant Cycadales.

Analysis of complete chloroplast genomes of Curcuma and the contribution to phylogeny and adaptive evolution

Curcuma is an important member of Zingiberaceae. Many species of this genus are widely used in traditional medicine and have important cultural value in East Asia. Among them, C. longa is considered to be the main source of curcumin and has a very wide range of uses. The rapid development of molecular phylogeny has deepened our understanding of taxonomy and evolution of Curcuma. However, little is known about the chloroplast genome phylogeny and the genetic bases of adaptative evolution. In this work, we sequenced the complete chloroplast genome of 4 Curcuma species. Curcuma chloroplast genomes showed highly conserved structures and the length ranged from 159,423 bp to 152,723 bp. A total of 133 genes were observed. Multiple repeats and simple sequence repeats (SSRs) were detected. By comparing with related species, 7 highly variable regions were identified as potential specific DNA barcodes for species identification. Phylogenetic analysis of complete plastome sequences and specific data sets revealed discordance with expected genus boundary. Chloroplast phylogenetic relationships were better predicted by geography than by morphological and nuclear DNA, indicating a substantial existence of introgression. 9 genes were proved to have high posteriori probability in positive selection analysis, and 4 of them (psbA, psbD, PetA and rbcL) closely related to photosynthesis, implying that chloroplast genes may had undergone positive selection pressure in evolution. These results are of great significance for us to understand the genetic basis, phylogeny and adaptive evolution of Curcuma chloroplast.

Comparative analysis of the complete chloroplast genome sequences of six species of Pulsatilla Miller, Ranunculaceae

BACKGROUND

Baitouweng is a traditional Chinese medicine with a long history of different applications. Although referred to as a single medicine, Baitouweng is actually comprised of many closely related species. It is therefore critically important to identify the different species that are utilized in these medicinal applications. Knowledge about their phylogenetic relationships can be derived from their chloroplast genomes and may provide additional insights into development of molecular markers.

METHODS

Genomic DNA was extracted from six species of Pulsatilla and then sequenced on an Illumina HiSeq 4000. Sequences were assembled into contigs by SOAPdenovo 2.04, aligned to the reference genome using BLAST, and then manually corrected. Genome annotation was performed by the online DOGMA tool. General characteristics of the cp genomes of the six species were analyzed and compared with closely related species. Additionally, phylogenetic trees were constructed, based on single nucleotide polymorphisms (SNPs) and 51 shared protein-coding gene sequences in the cp genome among all 31 species via maximum likelihood.

RESULTS

The size of cp genomes of P. chinensis (Bge.) Regel, P. chinensis (Bge.) Regel var. kissii (Mandl) S. H. Li et Y. H. Huang, P. cernua (Thunb.) Bercht. et Opiz f. plumbea J. X. Ji et Y. T. zhao, P. dahurica (Fisch.) Spreng, P. turczaninovii Kryl. et Serg, and P. cernua (Thunb.) Bercht. et Opiz. were 163,851 bp, 163,756 bp, 162,481 bp, 162,450 bp, 162,795 bp, and 162,924 bp, respectively. Each species included two inverted repeat regions, a small single-copy region, and a large single-copy region. A total of 134 genes were annotated, including 90 protein-coding genes, 36 tRNAs, and eight rRNAs across all species. In simple sequence repeat analysis, only P. dahurica was found to contain hexanucleotide repeats. A total of 26, 39, 32, 37, 32 and 43 large repeat sequences were identified in the genic regions of the six Pulsatilla species. Nucleotide diversity analysis revealed that the rpl36 gene and ccsA-ndhD region have the highest Pi value. In addition, two phylogenetic trees of the cp genomes were constructed, which laced all Pulsatilla species into one branch within Ranunculaceae.

CONCLUSIONS

We identified and analyzed the cp genome features of six species of P. Miller, with implications for species identification and phylogenetic analysis.

Characterization of the Complete Chloroplast Genome of Acer truncatum Bunge (Sapindales: Aceraceae): A New Woody Oil Tree Species Producing Nervonic Acid

Acer truncatum, which is a new woody oil tree species, is an important ornamental and medicinal plant in China. To assess the genetic diversity and relationships of A. truncatum, we analyzed its complete chloroplast (cp) genome sequence. The A. truncatum cp genome comprises 156,492 bp, with the large single-copy, small single-copy, and inverted repeat (IR) regions consisting of 86,010, 18,050, and 26,216 bp, respectively. The A. truncatum cp genome contains 112 unique functional genes (i.e., 4 rRNA, 30 tRNA, and 78 protein-coding genes) as well as 78 simple sequence repeats, 9 forward repeats, 1 reverse repeat, 5 palindromic repeats, and 7 tandem repeats. We analyzed the expansion/contraction of the IR regions in the cp genomes of six Acer species. A comparison of these cp genomes indicated the noncoding regions were more diverse than the coding regions. A phylogenetic analysis revealed that A. truncatum is closely related to A. miaotaiense. Moreover, a novel ycf4-cemA indel marker was developed for distinguishing several Acer species (i.e., A. buergerianum, A. truncatum, A. henryi, A. negundo, A. ginnala, and A. tonkinense). The results of the current study provide valuable information for future evolutionary studies and the molecular barcoding of Acer species.

Comparison of the Complete Eragrostis pilosa Chloroplast Genome with Its Relatives in Eragrostideae (Chloridoideae; Poaceae)

Eragrostis of the tribe Eragrostideae is a taxonomically complex genus, because of its polyploid nature and the presence of similar morphological characters among its species. However, the relationship between these morphologically indistinguishable species at the genomic level has not yet been investigated. Here, we report the complete chloroplast genome of E. pilosa and compare its genome structures, gene contents, simple sequence repeats (SSRs), sequence divergence, codon usage bias, and Kimura 2-parameter (K2P) interspecific genetic distances with those of other Eragrostideae species. The E. pilosa chloroplast genome was 134,815 bp in length and contained 132 genes and four regions, including a large single-copy region (80,100 bp), a small single-copy region (12,661 bp), and a pair of inverted repeats (21,027 bp). The average nucleotide diversity between E. pilosa and E. tef was estimated to be 0.011, and 0.01689 among all species. The minimum and maximum K2P interspecific genetic distance values were identified in psaA (0.007) and matK (0.029), respectively. Of 45 SSRs, eight were shared with E. tef, all of which were in the LSC region. Phylogenetic analysis resolved the monophyly of the sampled Eragrostis species and confirmed the close relationship between E. pilosa and E. tef. This study provides useful chlorophyll genomic information for further species identification and phylogenetic reconstruction of Eragrostis species.

Complete chloroplast genomes of all six Hosta species occurring in Korea: molecular structures, comparative, and phylogenetic analyses

BACKGROUND

The genus Hosta is a group of economically appreciated perennial herbs consisting of approximately 25 species that is endemic to eastern Asia. Due to considerable morphological variability, the genus has been well recognized as a group with taxonomic problems. Chloroplast is a cytoplasmic organelle with its own genome, which is the most commonly used for phylogenetic and genetic diversity analyses for land plants. To understand the genomic architecture of Hosta chloroplasts and examine the level of nucleotide and size variation, we newly sequenced four (H. clausa, H. jonesii, H. minor, and H. venusta) and analyzed six Hosta species (including the four, H. capitata and H. yingeri) distributed throughout South Korea.

RESULTS

The average size of complete chloroplast genomes for the Hosta taxa was 156,642 bp with a maximum size difference of ~ 300 bp. The overall gene content and organization across the six Hosta were nearly identical with a few exceptions. There was a single tRNA gene deletion in H. jonesii and four genes were pseudogenized in three taxa (H. capitata, H. minor, and H. jonesii). We did not find major structural variation, but there were a minor expansion and contractions in IR region for three species (H. capitata, H. minor, and H. venusta). Sequence variations were higher in non-coding regions than in coding regions. Four genic and intergenic regions including two coding genes (psbA and ndhD) exhibited the largest sequence divergence showing potential as phylogenetic markers. We found compositional codon usage bias toward A/T at the third position. The Hosta plastomes had a comparable number of dispersed and tandem repeats (simple sequence repeats) to the ones identified in other angiosperm taxa. The phylogeny of 20 Agavoideae (Asparagaceae) taxa including the six Hosta species inferred from complete plastome data showed well resolved monophyletic clades for closely related taxa with high node supports.

CONCLUSIONS

Our study provides detailed information on the chloroplast genome of the Hosta taxa. We identified nucleotide diversity hotspots and characterized types of repeats, which can be used for developing molecular markers applicable in various research area.