Whole plastome sequences from five ginger species facilitate marker development and define limits to barcode methodology

Highly divergent regions in the complete plastome sequences of Aquilaria are suitable for DNA barcoding applications including identifying species origin of agarwood products

Members of Aquilaria Lam. (Thymelaeaceae) are evergreen trees that are widely distributed in the Indomalesia region. Aquilaria is highly prized for its unique scented resin, agarwood, which is often the subject of unlawful trade activities. Survival of the tree is heavily threatened by destructive harvesting and agarwood poaching, leading to its protection under the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES). Unfortunately, an efficient species identification method, which is crucial to aid in the conservation efforts of Aquilaria is lacking. Here, we described our search for a suitable specific DNA barcode for Aquilaria species using eight complete plastome sequences. We identified five highly variable regions (HVR) (matK-rps16, ndhF-rpl32, psbJ-petA, trnD, and trnT-trnL) in the plastomes. These regions were further analyzed using the neighbor-joining (NJ) method to assess their ability at discriminating the eight species. Coupled with in silico primer design, two potential barcoding regions, psbJ-petA and trnT-trnL, were identified. Their strengths in species delimitation were evaluated individually and in combination, via DNA barcoding analysis. Our findings showed that the combined dataset, psbJ-petA + trnT-trnL, effectively resolved members of the genus Aquilaria by clustering all species into their respective clades. In addition, we demonstrated that the newly proposed DNA barcode was capable at identifying the species of origin of six commercial agarwood samples that were included as unknown samples. Such achievement offers a new technical advancement, useful in the combat against illicit agarwood trades and in assisting the conservation of these valuable species in natural populations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03479-1.

Deoxyribonucleic Acid Barcoding for the Identification of Botanicals

The Natural Herbal Products industry uses botanicals or herbs as raw materials for production of herbal products or dietary supplements. Recently, the demand for natural herbal products has increased tremendously and this has led to adulteration and to counterfeit herbal products. The present chapter deals with currently used molecular methods from "simple" single genomic regions to high-throughput whole genome or transcriptome sequencing methods used in the identification of botanicals.

The whole chloroplast genome of Mesua ferrea: Insight into the dynamic pattern of evolution and its comparison with species from recently diverged families

Mesua ferrea is an important source of timber, oil and herbal medicines. In the present investigation, we assembled the whole chloroplast genome of M. ferrea of size 161.4 kb. The genome contained 86 protein-coding genes, 38 tRNAs, 8 rRNA genes and exhibited a characteristic quadripartite structural orientation, with two inverted repeats (27,614 bp) separated by an LSC (88,746 bp) region and an SSC (27,614 bp) (17,470 bp). Interestingly, no gene loss was identified in the M. ferrea genome, contrary to what has been observed in other Clusioid species. We compared the chloroplast genome of M. ferrea with the chloroplast genome of Bonnetia and Garcinia belonging to Bonnetiaceae and Clusiaceae families. Overall, the compared genomes possess a similar synteny of gene order except for a small inversion in Garcinia species. M. ferrea has the largest chloroplast genome size in Clusioid clade owing to the lengthening of the LSC, IR, and non-coding regions. Substantial differences were observed in population of simple sequence repeats (SSRs) and RNA editing sites among the studied genomes. A comparative assessment of chloroplast genomes revealed five highly divergence regions: rpl32, trnS-GCU_trnG-UCC, petN-psbM, psbZ_trnG-GCC and ccsA_ndhD among the analyzed sequences. Phylogenetic analyses and sequence homology search indicate that M. ferrea is closely related to the Garcinia species.

Phylogenomics, plastome structure and species identification in Mahonia (Berberidaceae)

BACKGROUND

Elucidating the phylogenetic relationships within species-rich genera is essential but challenging, especially when lineages are assumed to have been going through radiation events. Mahonia Nutt. (Berberidaceae) is a genus with cosmopolitan distribution, comprising approximately 100 species, two of which are known as Caulis Mahoniae (M. bealei and M. fortunei) with crucial pharmacological significance in Chinese herbal medicine. Mahonia is a taxonomically challenging genus, and intrageneric phylogenetic relationships still need to be explored using genome data. Universal DNA barcodes and floral morphological attributes have limited discriminatory power in Mahonia.

RESULTS

We sequenced 17 representative plastomes and integrated three published plastome data together to conduct comparative and phylogenetic analyses. We found that Mahonia and Berberis share a large IR expansion (~ 12 kb), which is recognized as a typical character of Berberideae. Repeated sequences are revealed in the species of Mahonia, which are valuable for further population genetic studies. Using a comparative plastome analysis, we determined eight hypervariable regions whose discriminative power is comparable to that of the whole plastid genomes. The incongruence of the ITS and the plastome tree topologies may be ascribed to ancestral hybridization events and/or to incomplete lineage sorting. In addition, we suggest that leaf epidermal characters could help to distinguish closely related species in Mahonia.

CONCLUSIONS

We propose an integrative approach combining special barcodes and micromorphological traits to circumscribe Mahonia species. The results cast a new light on the development of an integrative method for accurate species circumscription and provide abundant genetic resources for further research on Mahonia.

"Omics" insights into plastid behavior toward improved carotenoid accumulation

Plastids are a group of diverse organelles with conserved carotenoids synthesizing and sequestering functions in plants. They optimize the carotenoid composition and content in response to developmental transitions and environmental stimuli. In this review, we describe the turbulence and reforming of transcripts, proteins, and metabolic pathways for carotenoid metabolism and storage in various plastid types upon organogenesis and external influences, which have been studied using approaches including genomics, transcriptomics, proteomics, and metabonomics. Meanwhile, the coordination of plastid signaling and carotenoid metabolism including the effects of disturbed carotenoid biosynthesis on plastid morphology and function are also discussed. The "omics" insight extends our understanding of the interaction between plastids and carotenoids and provides significant implications for designing strategies for carotenoid-biofortified crops.

Comparative phylogenetic analysis of complete plastid genomes of Renanthera (Orchidaceae)

Owing to its attractive flower shape and color, Renanthera (Orchidaceae), comprising about 19 species, has significant ornamental value as a houseplant, in floral design and in landscape gardens. Two species of Renanthera are categorized as endangered and critically endangered in China’s Red List and international trade in these orchids is currently strictly monitored by the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES). This paper reports on the de novo assembled and annotated plastome of four species of Renanthera; R. citrina, R. coccinea, R. imschootiana, and R. philippinensis. The length of the plastome sequences ranged from 144,673 bp (R. imschootiana) to 149,007 bp (R. coccinea) with GC content of 36.6–36.7%. The plastomes showed a typical quadripartite structure, including a large single-copy (84,241–86,404 bp), a small single-copy (11,468–12,167 bp), and a pair of inverted repeats (24,482–25,715 bp) regions. Of the 120 genes detected, 74 were protein coding, 38 were tRNA, and eight were rRNA genes. The plastome of Renanthera is rather conserved, but nucleotide variations that could distinguish them apart are noticeable—the total number of tandem repeats ranged from 62 (in R. imschootiana) to 74 (in R. citrina); while the number of long repeats ranged from 21 (in R. imschootiana and R. philippinensis) to 43 (in R. citrina). Three hypervariable regions (psbI-trnS-GCU, trnG-GCC, rpl32) were identified. Phylogenetic analyses based on the CDS using maximum likelihood (ML) and Bayesian inference (BI) revealed that Renanthera is closely related to Holcoglossum, Neofinetia, Pendulorchis, and Vanda. The relationship between the four species of Renanthera was fully resolved; a monophyletic clade was formed and R. coccinea was recorded as the first to diverge from the rest. The genetic data obtained from this study could serve as a useful resource for species identification in Renanthera as well as contribute to future research on the phylogenomics of Orchidaceae.

Mitochondrial Genome of Fagus sylvatica L. as a Source for Taxonomic Marker Development in the Fagales

European beech, Fagus sylvatica L., is one of the most important and widespread deciduous tree species in Central Europe and is widely managed for its hard wood. The complete DNA sequence of the mitochondrial genome of Fagus sylvatica L. was assembled and annotated based on Illumina MiSeq reads and validated using long reads from nanopore MinION sequencing. The genome assembled into a single DNA sequence of 504,715 bp in length containing 58 genes with predicted function, including 35 protein-coding, 20 tRNA and three rRNA genes. Additionally, 23 putative protein-coding genes were predicted supported by RNA-Seq data. Aiming at the development of taxon-specific mitochondrial genetic markers, the tool SNPtax was developed and applied to select genic SNPs potentially specific for different taxa within the Fagales. Further validation of a small SNP set resulted in the development of four CAPS markers specific for Fagus, Fagaceae, or Fagales, respectively, when considering over 100 individuals from a total of 69 species of deciduous trees and conifers from up to 15 families included in the marker validation. The CAPS marker set is suitable to identify the genus Fagus in DNA samples from tree tissues or wood products, including wood composite products.

Applied Barcoding: The Practicalities of DNA Testing for Herbals

DNA barcoding is a widely accepted technique for the identification of plant materials, and its application to the authentication of commercial medicinal plants has attracted significant attention. The incorporation of DNA-based technologies into the quality testing protocols of international pharmacopoeias represents a step-change in status, requiring the establishment of standardized, reliable and reproducible methods. The process by which this can be achieved for any herbal medicine is described, using Hypericum perforatum L. (St John's Wort) and potential adulterant Hypericum species as a case study. A range of practical issues are considered including quality control of DNA sequences from public repositories and the construction of individual curated databases, choice of DNA barcode region(s) and the identification of informative polymorphic nucleotide sequences. A decision tree informs the structure of the manuscript and provides a template to guide the development of future DNA barcode tests for herbals.

Phoenix phylogeny, and analysis of genetic variation in a diverse collection of date palm (Phoenix dactylifera) and related species

Date palm (Phoenix dactylifera), one of the most ancient crops, is grown commercially in >30 countries. Using whole plastome assemblies, phylogenetic analyses revealed that cultivated date palm accessions share the same clade with P hoenix sylvestris, P hoenix pusilla and P hoenix acaulis, which are native to the Indian subcontinent, and Phoenix caespitosa that is native to the Arabian Peninsula and the deserts of Somalia. Analysis of genetic diversity and genetic relationships among date palm accessions from 13 producing countries involved 195 date palm accessions that were genotyped at 19 microsatellite loci. Extensive genetic diversity was observed, with many accessions heterozygous for most markers in this clonally propagated crop. The average number of alleles per locus (42.1), expected heterozygosity (0.8), observed heterozygosity (0.47) and fixation indices (F_ST = 0.42) demonstrated substantial genetic diversity and population structure. Iraqi accessions were found to have the richest allelic diversity, and the most private alleles. The model-based Bayesian method indicated that these accessions could be broadly divided into two structure groups, one group with predominantly African accessions and another predominantly Asian. Some germplasm, especially from Tunisia and Iraq, deviated from this generalization. Many accessions in the STRUCTURE-derived groups were found to be genetic admixtures, with gene flow between Asian and African groups. Indian and Pakistani date palms were found to be most closely related to North African germplasm.

Pan-plastome approach empowers the assessment of genetic variation in cultivated Capsicum species

Pepper species (Capsicum spp.) are widely used as food, spice, decoration, and medicine. Despite the recent old-world culinary impact, more than 50 commercially recognized pod types have been recorded worldwide from three taxonomic complexes (A, B, and P). The current study aimed to apply a pan-plastome approach to resolve the plastomic boundaries among those complexes and identify effective loci for the taxonomical resolution and molecular identification of the studied species/varieties. High-resolution pan-plastomes of five species and two varieties were assembled and compared from 321 accessions. Phyloplastomic and network analyses clarified the taxonomic position of the studied species/varieties and revealed a pronounced number of accessions to be the rare and endemic species, C. galapagoense, that were mistakenly labeled as C. annuum var. glabriusculum among others. Similarly, some NCBI-deposited plastomes were clustered differently from their labels. The rpl23-trnI intergenic spacer contained a 44 bp tandem repeat that, in addition to other InDels, was capable of discriminating the investigated Capsicum species/varieties. The rps16-trnQ/rbcL-accD/ycf3-trnS gene set was determined to be sufficiently polymorphic to retrieve the complete phyloplastomic signal among the studied Capsicum spp. The pan-plastome approach was shown to be useful in resolving the taxonomical complexes, settling the incomplete lineage sorting conflict and developing a molecular marker set for Capsicum spp. identification.

Phylogenomic Approaches to DNA Barcoding of Herbal Medicines: Developing Clade-Specific Diagnostic Characters for Berberis

DNA barcoding of herbal medicines has been mainly concerned with authentication of products in trade and has raised awareness of species substitution and adulteration. More recently DNA barcodes have been included in pharmacopoeias, providing tools for regulatory purposes. The commonly used DNA barcoding regions in plants often fail to resolve identification to species level. This can be especially challenging in evolutionarily complex groups where incipient or reticulate speciation is ongoing. In this study, we take a phylogenomic approach, analyzing whole plastid sequences from the evolutionarily complex genus Berberis in order to develop DNA barcodes for the medicinally important species Berberis aristata. The phylogeny reconstructed from an alignment of ∼160 kbp of chloroplast DNA for 57 species reveals that the pharmacopoeial species in question is polyphyletic, complicating development of a species-specific DNA barcode. Instead we propose a DNA barcode that is clade specific, using our phylogeny to define Operational Phylogenetic Units (OPUs). The plastid alignment is then reduced to small, informative DNA regions including nucleotides diagnostic for these OPUs. These DNA barcodes were tested on commercial samples, and shown to discriminate plants in trade and therefore to meet the requirement of a pharmacopoeial standard. The proposed method provides an innovative approach for inferring DNA barcodes for evolutionarily complex groups for regulatory purposes and quality control.

Green Biotechnology: A Brief Update on Plastid Genome Engineering

Plant genetic engineering has become an inevitable tool in the molecular breeding of crops. Significant progress has been made in the generation of novel plastid transformation vectors and optimized transformation protocols. There are several advantages of plastid genome engineering over conventional nuclear transformation. Some of the advantages include multigene engineering by expression of biosynthetic pathway genes as operons, extremely high-level expression of protein accumulation, lack of transgene silencing, etc. Transgene containment owing to maternal inheritance is another important advantage of plastid genome engineering. Chloroplast genome modification usually results in alteration of several thousand plastid genome copies in a cell. Several therapeutic proteins, edible vaccines, antimicrobial peptides, and industrially important enzymes have been successfully expressed in chloroplasts so far. Here, we critically recapitulate the latest developments in plastid genome engineering. Latest advancements in plastid genome sequencing are briefed. In addition, advancement of extending the toolbox for plastid engineering for selected applications in the area of molecular farming and production of industrially important enzyme is briefed.

Discovery of Lineage-Specific Genome Change in Rice Through Analysis of Resequencing Data

New mutations are rare, which makes their discovery laborious and time-consuming. Arthur and Bennetzen describe an approach for enriching recent mutations that relies only on a reference genome sequence and resequencing data for other...

Genome comparisons provide information on the nature of genetic change, but such comparisons are challenged to differentiate the importance of the actual sequence change processes relative to the role of selection. This problem can be overcome by identifying changes that have not yet had the time to undergo millions of years of natural selection. We describe a strategy to discover accession-specific changes in the rice genome using an abundant resource routinely provided for many genome analyses, resequencing data. The sequence of the fully sequenced rice genome from variety Nipponbare was compared to the pooled (∼114×) resequencing data from 126 japonica rice accessions to discover “Nipponbare-specific” sequences. Analyzing nonrepetitive sequences, 8504 “candidate” Nipponbare-specific changes were detected, of which around two-thirds are true novel sequence changes and the rest are predicted genome sequencing errors. Base substitutions outnumbered indels in this data set by > 28:1, with ∼8:5 bias toward transversions over transitions, and no transposable element insertions or excisions were observed. These results indicate that the strategy employed is effective for finding recent sequence changes, sequencing errors, and rare alleles in any organism that has both a reference genome sequence and a wealth of resequencing data.

Gene loss and genome rearrangement in the plastids of five Hemiparasites in the family Orobanchaceae

BACKGROUND

The chloroplast genomes (plastome) of most plants are highly conserved in structure, gene content, and gene order. Parasitic plants, including those that are fully photosynthetic, often contain plastome rearrangements. These most notably include gene deletions that result in a smaller plastome size. The nature of gene loss and genome structural rearrangement has been investigated in several parasitic plants, but their timing and contributions to the adaptation of these parasites requires further investigation, especially among the under-studied hemi-parasites.

RESULTS

De novo sequencing, assembly and annotation of the chloroplast genomes of five photosynthetic parasites from the family Orobanchaceae were employed to investigate plastome dynamics. Four had major structural rearrangements, including gene duplications and gene losses, that differentiated the taxa. The facultative parasite Aureolaria virginica had the most similar genome content to its close non-parasitic relative, Lindenbergia philippensis, with similar genome size and organization, and no differences in gene content. In contrast, the facultative parasite Buchnera americana and three obligate parasites in the genus Striga all had enlargements of their plastomes, primarily caused by expansion within the large inverted repeats (IRs) that are a standard plastome feature. Some of these IR increases were shared by multiple investigated species, but others were unique to particular lineages. Gene deletions and pseudogenization were also both shared and lineage-specific, with particularly frequent and independent loss of the ndh genes involved in electron recycling.

CONCLUSIONS

Five new plastid genomes were fully assembled and compared. The results indicate that plastome instability is common in parasitic plants, even those that retain the need to perform essential plastid functions like photosynthesis. Gene losses were slow and not identical across taxa, suggesting that different lineages had different uses or needs for some of their plastome gene content, including genes involved in some aspects of photosynthesis. Recent repeat region extensions, some unique to terminal species branches, were observed after the divergence of the Buchnera/Striga clade, suggesting that this otherwise rare event has some special value in this lineage.