Comparison of the chloroplast genomics of nine endangered Habenaria species and phylogenetic analysis.

BACKGROUND: Habenaria, a genus in the family Orchidaceae, are the nearly cosmopolitan orchids, and most species have significant medicinal and ornamental values. Despite the morphological and molecular data that have been studied in recent years, the phylogenetic relationship is still unclear. RESULTS: We sequenced, assembled, and annotated the chloroplast (cp) genomes of two species (Habenaria aitchisonii Rchb.f. and Habenaria tibetica Schltr.ex Limpricht) of Habenaria grown on the Qinghai-Tibetan Plateau (QTP), and compared them with seven previously published cp genomes which may aid in the genomic profiling of these species. The two genomes ranged from 155,259-155,269 bp in length and both included 132 genes, encoding 86 proteins, 38 tRNAs and 8 rRNAs. In the cp genomes, the tandem repeats (797), SSRs (2195) and diverse loci (3214) were identified. Comparative analyses of codon usage, amino frequency, microsatellite, oligo repeats and transition and transversion substitutions revealed similarities between the species. Moreover, we identified 16 highly polymorphic regions with a nucleotide diversity above 0.02, which may be suitable for robust authentic barcoding and inferring in the phylogeny of Habenaria species. Among the polymorphic regions, positive selection was significantly exerted on several genes, such as cemA, petA, and ycf1. This finding may suggest an important adaptation strategy for the two Habenaria species on the QTP. The phylogenetic relationship revealed that H. aitchisonii and H. tibetica were more closely related to each other than to the other species, and the other seven species were clustered in three groups. In addition, the estimated divergence time suggested that the two species separated from the others approximately 0.39 Mya in the Neogene period. Our findings also suggest that Habenaria can be divided into different sections. CONCLUSIONS: The results of this study enriched the genomics resources of Habenaria, and SSR marker may aid in the conservation management of two endangered species.

Chloroplast genome evolution in the Dracunculus clade (Aroideae, Araceae).

Chloroplast (cp) genomes are considered important for the study of lineage-specific molecular evolution, population genetics, and phylogenetics. Our aim here was to elucidate the molecular evolution in cp genomes of species in the Dracunculus clade (Aroideae, Araceae). We report de novo assembled cp genomes for eight species from eight genera and also retrieved cp genomes of four species from the National Center for Biotechnology Information (NCBI). The cp genomes varied in size from 162,424 bp to 176,835 bp. Large Single Copy (LSC) region ranged in size from 87,141 bp to 95,475 bp; Small Single Copy (SSC) from 14,338 bp to 23,981 bp; and Inverted Repeats (IRa and IRb) from 25,131 bp to 32,708 bp. The expansion in inverted repeats led to duplication of ycf1 genes in four species. The genera showed high similarity in gene content and yielded 113 unique genes (79 protein-coding, 4 rRNA, and 30 tRNA genes). Codon usage, amino acid frequency, RNA editing sites, microsatellites repeats, transition and transversion substitutions, and synonymous and non-synonymous substitutions were also similar across the clade. A previous study reported deletion of ycf1, accD, psbE, trnL-CAA, and trnG-GCC genes in four Amorphophallus species. Our study supports conservative structure of cp genomes in the Dracunculus clade including Amorphophallus species and does not support gene deletion mentioned above. We also report suitable polymorphic loci based on comparative analyses of Dracunculus clade species, which could be useful for phylogenetic inference. Overall, the current study broad our knowledge about the molecular evolution of chloroplast genome in aroids.

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

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

Characterization of Withania somnifera chloroplast genome and its comparison with other selected species of Solanaceae.

Withania somnifera (L) Dunal, a wonder herb of family Solanaceae, has multiple medicinal properties. Here, we reported the chloroplast genome sequence of Withania somnifera (154,386 bp) which comprises of a large single copy region (85,688 bp), and a small single copy region (18,464 bp), separated by a pair of large inverted repeats (25,117 bp). The chloroplast genome has 132 genes including 86 protein-coding, 37 tRNAs and 8 rRNAs. Comparison of chloroplast genomes of Withania somnifera with four other Solanaceae species revealed similarities in genomic features, including structure, nucleotide content, codon usage, RNA editing sites, simple sequence repeats (SSRs), oligonucleotide repeats, and tandem repeats. We identified 147 simple sequence repeats in protein-coding, and 229 in non-protein-coding regions. We observed numerous post-transcriptional substitutions of Serine to Leucine, specifically at the second nucleotide position of the codon. Maximum likelihood and maximum parsimony tree reconstructed displayed Withania somnifera a sister taxon of Physalis peruviana.

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

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

Complete Chloroplast Genomes of Anthurium huixtlense and Pothos scandens (Pothoideae, Araceae): Unique Inverted Repeat Expansion and Contraction Affect Rate of Evolution.

The subfamily Pothoideae belongs to the ecologically important plant family Araceae. Here, we report the chloroplast genomes of two species of the subfamily Pothoideae: Anthurium huixtlense (size: 163,116 bp) and Pothos scandens (size: 164,719 bp). The chloroplast genome of P. scandens showed unique contraction and expansion of inverted repeats (IRs), thereby increasing the size of the large single-copy region (LSC: 102,956 bp) and decreasing the size of the small single-copy region (SSC: 6779 bp). This led to duplication of many single-copy genes due to transfer to IR regions from the small single-copy (SSC) region, whereas some duplicate genes became single copy due to transfer to large single-copy regions. The rate of evolution of protein-coding genes was affected by the contraction and expansion of IRs; we found higher mutation rates for genes that exist in single-copy regions as compared to those in IRs. We found a 2.3-fold increase of oligonucleotide repeats in P. scandens when compared with A. huixtlense, whereas amino acid frequency and codon usage revealed similarities. The ratio of transition to transversion mutations was 2.26 in P. scandens and 2.12 in A. huixtlense. Transversion mutations mostly translated in non-synonymous substitutions. The phylogenetic inference of the limited species showed the monophyly of the Araceae subfamilies. Our study provides insight into the molecular evolution of chloroplast genomes in the subfamily Pothoideae and family Araceae.

Comparative Plastomics of Plantains (Plantago, Plantaginaceae) as a Tool for the Development of Species-Specific DNA Barcodes.

Plantago (plantains, Plantaginaceae) is a cosmopolitan genus including over 250 species used as functional foods, forage, and traditional medicine. Among them, Plantago lanceolata is commonly used as an ingredient of herbal products, but the close similarity to other Plantago species can cause misidentifications with potentially serious consequences for product safety/quality. To test the possibility of developing species-specific barcoding markers, we de novo assembled plastome sequences of individuals of Plantago argentea, Plantago atrata, P. lanceolata, and Plantago maritima. These genomes were characterized in comparison with both previously sequenced conspecific accessions and other publicly available plastomes, thus providing an assessment of both intraspecific and interspecific genetic variation in Plantago plastomes. Additionally, molecular evolutionary analyses indicated that eleven protein-coding genes involved in different plastid functions in Plantago plastomes underwent positive selection, suggesting they might have contributed to enhancing species' adaptation during the evolutionary history of Plantago. While the most variable mutational hotspots in Plantago plastomes were not suitable for the development of species-specific molecular markers, species-specific polymorphisms could discriminate P. lanceolata from its closest relatives. Taken together, these results highlight the potential of plastome sequencing for the development of molecular markers to improve the identification of species with relevance in herbal products.

Characterization of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence.

Anisodus tanguticus Maxim. (Solanaceae), a traditional endangered Tibetan herb, is endemic to the Qinghai-Tibet Plateau. Here, we report the de novo assembled chloroplast (cp) genome sequences of A. tanguticus (155,765 bp). The cp contains a pair of inverted repeated (IRa and IRb) regions of 25,881 bp that are separated by a large single copy (LSC) region (86,516 bp) and a small single copy SSC (17,487 bp) region. A total of 132 functional genes were annotated in the cp genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Moreover, 199 simple sequence repeats (SSR) and 65 repeat structures were detected. Comparative plastome analyses revealed a conserved gene order and high similarity of protein-coding sequences. The A. tanguticus cp genome exhibits contraction and expansion, which differs from Przewalskia tangutica and other related Solanaceae species. We identified 30 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which may be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Anisodus and family Solanaceae. Analysis of the Ka/Ks ratios of the Hyoscyameae tribe revealed significant positive selection exerted on the cemA, rpoC2, and clpP genes, which suggests that protein metabolism may be an important strategy for A. tanguticus and other species in Hyoscyameae in adapting to the adverse environment on the Qinghai-Tibetan Plateau. Phylogenetic analysis revealed that A. tanguticus clustered closer with Hyoscyamus niger than P. tangutica. Our results provide reliable genetic information for future exploration of the taxonomy and phylogenetic evolution of the Hyoscyameae tribe and related species.

Comparative plastome analysis of Blumea, with implications for genome evolution and phylogeny of Asteroideae.

The genus Blumea (Asteroideae, Asteraceae) comprises about 100 species, including herbs, shrubs, and small trees. Previous studies have been unable to resolve taxonomic issues and the phylogeny of the genus Blumea due to the low polymorphism of molecular markers. Therefore, suitable polymorphic regions need to be identified. Here, we de novo assembled plastomes of the three Blumea species B. oxyodonta, B. tenella, and B. balsamifera and compared them with 26 other species of Asteroideae after correction of annotations. These species have quadripartite plastomes with similar gene content, genome organization, and inverted repeat contraction and expansion comprising 113 genes, including 80 protein-coding, 29 transfer RNA, and 4 ribosomal RNA genes. The comparative analysis of codon usage, amino acid frequency, microsatellite repeats, oligonucleotide repeats, and transition and transversion substitutions has revealed high resemblance among the newly assembled species of Blumea. We identified 10 highly polymorphic regions with nucleotide diversity above 0.02, including rps16-trnQ, ycf1, ndhF-rpl32, petN-psbM, and rpl32-trnL, and they may be suitable for the development of robust, authentic, and cost-effective markers for barcoding and inference of the phylogeny of the genus Blumea. Among these highly polymorphic regions, five regions also co-occurred with oligonucleotide repeats and support use of repeats as a proxy for the identification of polymorphic loci. The phylogenetic analysis revealed a close relationship between Blumea and Pluchea within the tribe Inuleae. At tribe level, our phylogeny supports a sister relationship between Astereae and Anthemideae rooted as Gnaphalieae, Calenduleae, and Senecioneae. These results are contradictory to recent studies which reported a sister relationship between "Senecioneae and Anthemideae" and "Astereae and Gnaphalieae" or a sister relationship between Astereae and Gnaphalieae rooted as Calenduleae, Anthemideae, and then Senecioneae using nuclear genome sequences. The conflicting phylogenetic signals observed at the tribal level between plastidt and nuclear genome data require further investigation.

Pseudogenization of the chloroplast threonine (trnT-GGU) gene in the sunflower family (Asteraceae).

The chloroplast genome evolves through the course of evolution. Various types of mutational events are found within the chloroplast genome, including insertions-deletions (InDels), substitutions, inversions, gene rearrangement, and pseudogenization of genes. The pseudogenization of the chloroplast threonine (trnT-GGU) gene was previously reported in Cryptomeria japonica (Cupressaceae), Pelargonium × hortorum (Geraniaceae), and Anaphalis sinica and Leontopodium leiolepis of the tribe Gnaphalieae (Asteroideae, Asteraceae). Here, we performed a broad analysis of the trnT-GGU gene among the species of 13 subfamilies of Asteraceae and found this gene as a pseudogene in core Asteraceae (Gymnarrhenoideae, Cichorioideae, Corymbioideae, and Asteroideae), which was linked to an insertion event within the 5' acceptor stem and is not associated with ecological factors such as habit, habitat, and geographical distribution of the species. The pseudogenization of trnT-GGU was not predicted in codon usage, indicating that the superwobbling phenomenon occurs in core Asteraceae in which a single transfer RNA (trnT-UGU) decodes all four codons of threonine. To the best of our knowledge, this is the first evidence of a complete clade of a plant species using the superwobbling phenomenon for translation.

Comparative Plastomics of Ashwagandha (Withania, Solanaceae) and Identification of Mutational Hotspots for Barcoding Medicinal Plants.

Within the family Solanaceae, Withania is a small genus belonging to the Solanoideae subfamily. Here, we report the de novo assembled chloroplast genome sequences of W. coagulans, W. adpressa, and W. riebeckii. The length of these genomes ranged from 154,162 to 154,364 base pairs (bp). These genomes contained a pair of inverted repeats (IRa and IRb) ranging from 25,029 to 25,071 bp that were separated by a large single-copy (LSC) region of 85,635-85,765 bp and a small single-copy (SSC) region of 18,457-18,469 bp. We analyzed the structural organization, gene content and order, guanine-cytosine content, codon usage, RNA-editing sites, microsatellites, oligonucleotide and tandem repeats, and substitutions of Withania plastomes, which revealed high similarities among the species. Comparative analysis among the Withania species also highlighted 10 divergent hotspots that could potentially be used for molecular marker development, phylogenetic analysis, and species identification. Furthermore, our analyses showed that even three mutational hotspots (rps4-trnT, trnM-atpE, and rps15) were sufficient to discriminate the Withania species included in current study.

Plastid genomics of Nicotiana (Solanaceae): insights into molecular evolution, positive selection and the origin of the maternal genome of Aztec tobacco (Nicotiana rustica).

Species of the genus Nicotiana (Solanaceae), commonly referred to as tobacco plants, are often cultivated as non-food crops and garden ornamentals. In addition to the worldwide production of tobacco leaves, they are also used as evolutionary model systems due to their complex development history tangled by polyploidy and hybridization. Here, we assembled the plastid genomes of five tobacco species: N. knightiana, N. rustica, N. paniculata, N. obtusifolia and N. glauca. De novo assembled tobacco plastid genomes had the typical quadripartite structure, consisting of a pair of inverted repeat (IR) regions (25,323-25,369 bp each) separated by a large single-copy (LSC) region (86,510-86,716 bp) and a small single-copy (SSC) region (18,441-18,555 bp). Comparative analyses of Nicotiana plastid genomes with currently available Solanaceae genome sequences showed similar GC and gene content, codon usage, simple sequence and oligonucleotide repeats, RNA editing sites, and substitutions. We identified 20 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which could be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Nicotiana and family Solanaceae. Our comparative plastid genome analysis revealed that the maternal parent of the tetraploid N. rustica was the common ancestor of N. paniculata and N. knightiana, and the later species is more closely related to N. rustica. Relaxed molecular clock analyses estimated the speciation event between N. rustica and N. knightiana appeared 0.56 Ma (HPD 0.65-0.46). Biogeographical analysis supported a south-to-north range expansion and diversification for N. rustica and related species, where N. undulata and N. paniculata evolved in North/Central Peru, while N. rustica developed in Southern Peru and separated from N. knightiana, which adapted to the Southern coastal climatic regimes. We further inspected selective pressure on protein-coding genes among tobacco species to determine if this adaptation process affected the evolution of plastid genes. These analyses indicate that four genes involved in different plastid functions, including DNA replication (rpoA) and photosynthesis (atpB, ndhD and ndhF), came under positive selective pressure as a result of specific environmental conditions. Genetic mutations in these genes might have contributed to better survival and superior adaptations during the evolutionary history of tobacco species.

Comparison of Chloroplast Genomes among Species of Unisexual and Bisexual Clades of the Monocot Family Araceae.

The chloroplast genome provides insight into the evolution of plant species. We de novo assembled and annotated chloroplast genomes of four genera representing three subfamilies of Araceae: Lasia spinosa (Lasioideae), Stylochaeton bogneri, Zamioculcas zamiifolia (Zamioculcadoideae), and Orontium aquaticum (Orontioideae), and performed comparative genomics using these chloroplast genomes. The sizes of the chloroplast genomes ranged from 163,770 bp to 169,982 bp. These genomes comprise 113 unique genes, including 79 protein-coding, 4 rRNA, and 30 tRNA genes. Among these genes, 17-18 genes are duplicated in the inverted repeat (IR) regions, comprising 6-7 protein-coding (including trans-splicing gene rps12), 4 rRNA, and 7 tRNA genes. The total number of genes ranged between 130 and 131. The infA gene was found to be a pseudogene in all four genomes reported here. These genomes exhibited high similarities in codon usage, amino acid frequency, RNA editing sites, and microsatellites. The oligonucleotide repeats and junctions JSB (IRb/SSC) and JSA (SSC/IRa) were highly variable among the genomes. The patterns of IR contraction and expansion were shown to be homoplasious, and therefore unsuitable for phylogenetic analyses. Signatures of positive selection were seen in three genes in S. bogneri, including ycf2, clpP, and rpl36. This study is a valuable addition to the evolutionary history of chloroplast genome structure in Araceae.

Neuroprotective, antidiabetic and antioxidant effect of Hedera nepalensis and lupeol against STZ + AlCl3 induced rats model.

PURPOSE: This study was aimed to evaluate the effect of Hedera nepalensis crude extract (HNC) and its isolated compound lupeol on antioxidant defence system, biochemical parameters and behavioural indices of Alzheimer disease generated in diabetic rats. METHODS: To evaluate the effect of the plant extract and lupeol, symptoms of Alzheimer and diabetes were induced in rats by STZ + AlCl3 treatment. Glucose level was measured with glucometer followed by antioxidant and biochemical assessment of the treated and untreated animals. Behavioural response of the rats was determined by Elevated Plus Maze (EPM) test and Morris Water Maze (MWM) test followed by determination of brain neurotransmitters by HPLC. RESULTS: HNC significantly reduced blood glucose level in a time dependent manner and elevated liver function markers were significantly (P < 0.05) reinstated to normal levels. HNC showed increase in level of catalase (CAT), superoxide dismutase (SOD) and reduced glutathione (GSH). HPLC quantification revealed that HNC treatment led to significant (p < 0.001) elevation in the level of neurotransmitters (dopamine and serotonin) in the midbrain region as compared to Alzheimer control (AC) group. EPM and MWM test showed decrease in cognitive and memory impairment in a rat group treated with HNC as compared to AC group. CONCLUSION: Overall, results showed that H. nepalensis has therapeutic potential for the treatment of diseases like Alzheimer and diabetes. Graphical abstract Therapeutic effect of Hedera nepalensis K. Koch and lupeol against STZ + AICI3 induced diabetic rats model.