Complete sequences of organelle genomes from the medicinal plant Rhazya stricta (Apocynaceae) and contrasting patterns of mitochondrial genome evolution across asterids

Mitochondrial genome of Hancornia speciosa gomes: intergenic regions containing retrotransposons and predicted genes

BACKGROUND

Plant mitochondrial genomes are characterized by high homologous recombination, extensive intergenic spacers, conservation in DNA sequences, and gene content. The Hancornia genus belongs to the Apocynaceae family, with H. speciosa Gomes being the sole species in the genus. It is an siganificant commercial fruit crop; however, only a number of studies have been conducted. In this study, we present the mitochondrial genome of H. speciosa and compare it with other mitochondrial genomes within the Apocynaceae family.

METHODS AND RESULTS

A total of 2.8 Gb of Illumina paired-end reads were used to obtain the mitogenome, resulting in 22 contigs that were merged using 6.1 Gb of Illumina mate-pair reads to obtain a circular chromosome. The mitochondrial genome of H. speciosa is circular, containing 63 predicted functional genes, spanning a length of 741,811 bp, with a CG content of 44%. Within the mitogenome, 50 chloroplast DNA sequences, equivalent to 1.72% of the genome, were detected. However, intergenic spaces accounted for 703,139 bp (94.79% of the genome), and 287 genes were predicted, totaling 173,721 bp.

CONCLUSION

This suggests the incorporation of nuclear DNA into the mitogenome of H. speciosa and self duplication. Comparative analysis among the mitogenomes in the Apocynaceae family revealed a diversity in the structure mediated by recombination, with similar gene content and large intergenic spaces.

Intrageneric structural variation in organelle genomes from the genus Dystaenia (Apiaceae): genome rearrangement and mitochondrion-to-plastid DNA transfer

Introduction

During plant evolution, intracellular DNA transfer (IDT) occurs not only from organelles to the nucleus but also between organelles. To further comprehend these events, both organelle genomes and transcriptomes are needed.

Methods

In this study, we constructed organelle genomes and transcriptomes for two Dystaenia species and described their dynamic IDTs between their nuclear and mitochondrial genomes, or plastid and mitochondrial genomes (plastome and mitogenome).

Results and Discussion

We identified the putative functional transfers of the mitochondrial genes 5' rpl2, rps10, rps14, rps19, and sdh3 to the nucleus in both Dystaenia species and detected two transcripts for the rpl2 and sdh3 genes. Additional transcriptomes from the Apicaceae species also provided evidence for the transfers and duplications of these mitochondrial genes, showing lineage-specific patterns. Intrageneric variations of the IDT were found between the Dystaenia organelle genomes. Recurrent plastid-to-mitochondrion DNA transfer events were only identified in the D. takeshimana mitogenome, and a pair of mitochondrial DNAs of plastid origin (MIPTs) may generate minor alternative isoforms. We only found a mitochondrion-to-plastid DNA transfer event in the D. ibukiensis plastome. This event may be linked to inverted repeat boundary shifts in its plastome. We inferred that the insertion region involved an MIPT that had already acquired a plastid sequence in its mitogenome via IDT. We propose that the MIPT acts as a homologous region pairing between the donor and recipient sequences. Our results provide insight into the evolution of organelle genomes across the family Apiaceae.

Comparative analysis of the complete mitogenome of Geoffroea decorticans: a native tree surviving in the Atacama Desert

Chañar (Geoffroea decorticans (Gill., ex Hook. & Arn.) Burkart) has been highly significant for indigenous people in the Atacama Desert for over 3,000 years. Through evolutionary processes, the G. decorticans mitogenome likely underwent changes facilitating its adaptation to the extreme conditions of the Atacama Desert. Here, we compare the mitochondrial genome of G. decorticans with those of other Papilionoideae family species. The complete mitogenome of G. decorticans was sequenced and assembled, making it the first in the genus Geoffroea. The mitogenome contained 383,963 base pairs, consisting of 33 protein coding genes, 21 transfer RNA genes, and 3 ribosomal RNA genes. The Chañar mitogenome is relatively compact, and has two intact genes (sdh4 and nad1) which were not observed in most other species. Additionally, Chañar possessed the highest amount of mitochondrial DNA of plastid origin among angiosperm species. The phylogenetic analysis of the mitogenomes of Chañar and 12 other taxa displayed a high level of consistency in taxonomic classification, when compared to those of the plastid genome. Atp8 was subjected to positive selection, while the ccmFc and rps1 were subjected to neutral selection. This study provides valuable information regarding its ability to survive the extreme environmental conditions of the Atacama Desert.

Rate accelerations in plastid and mitochondrial genomes of Cyperaceae occur in the same clades

Cyperaceae, the second largest family in the monocot order Poales, comprises >5500 species and includes the genus Eleocharis with ∼ 250 species. A previous study of complete plastomes of two Eleocharis species documented extensive structural heteroplasmy, gene order changes, high frequency of dispersed repeats along with gene losses and duplications. To better understand the phylogenetic distribution of gene and intron content as well as rates and patterns of sequence evolution within and between mitochondrial and plastid genomes of Eleocharis and Cyperaceae, an additional 29 Eleocharis organelle genomes were sequenced and analyzed. Eleocharis experienced extensive gene loss in both genomes while loss of introns was mitochondria-specific. Eleocharis has higher rates of synonymous (dS) and nonsynonymous (dN) substitutions in the plastid and mitochondrion than most sampled angiosperms, and the pattern was distinct from other eudicot lineages with accelerated rates. Several clades showed higher dS and dN in mitochondrial genes than in plastid genes. Furthermore, nucleotide substitution rates of mitochondrial genes were significantly accelerated on the branch leading to Cyperaceae compared to most angiosperms. Mitochondrial genes of Cyperaceae exhibited dramatic loss of RNA editing sites and a negative correlation between RNA editing and dS values was detected among angiosperms. Mutagenic retroprocessing and dysfunction of DNA replication, repair and recombination genes are the most likely cause of striking rate accelerations and loss of edit sites and introns in Eleocharis and Cyperaceae organelle genomes.

Apostasia Mitochondrial Genome Analysis and Monocot Mitochondria Phylogenomics

Apostasia shenzhenica belongs to the subfamily Apostasioideae and is a primitive group located at the base of the Orchidaceae phylogenetic tree. However, the A. shenzhenica mitochondrial genome (mitogenome) is still unexplored, and the phylogenetic relationships between monocots mitogenomes remain unexplored. In this study, we discussed the genetic diversity of A. shenzhenica and the phylogenetic relationships within its monocotyledon mitogenome. We sequenced and assembled the complete mitogenome of A. shenzhenica, resulting in a circular mitochondrial draft of 672,872 bp, with an average read coverage of 122× and a GC content of 44.4%. A. shenzhenica mitogenome contained 36 protein-coding genes, 16 tRNAs, two rRNAs, and two copies of nad4L. Repeat sequence analysis revealed a large number of medium and small repeats, accounting for 1.28% of the mitogenome sequence. Selection pressure analysis indicated high mitogenome conservation in related species. RNA editing identified 416 sites in the protein-coding region. Furthermore, we found 44 chloroplast genomic DNA fragments that were transferred from the chloroplast to the mitogenome of A. shenzhenica, with five plastid-derived genes remaining intact in the mitogenome. Finally, the phylogenetic analysis of the mitogenomes from A. shenzhenica and 28 other monocots showed that the evolution and classification of most monocots were well determined. These findings enrich the genetic resources of orchids and provide valuable information on the taxonomic classification and molecular evolution of monocots.

Comparative analysis of the organelle genomes of three Rhodiola species provide insights into their structural dynamics and sequence divergences

BACKGROUND

Plant organelle genomes are a valuable resource for evolutionary biology research, yet their genome architectures, evolutionary patterns and environmental adaptations are poorly understood in many lineages. Rhodiola species is a type of flora mainly distributed in highland habitats, with high medicinal value. Here, we assembled the organelle genomes of three Rhodiola species (R. wallichiana, R. crenulata and R. sacra) collected from the Qinghai-Tibet plateau (QTP), and compared their genome structure, gene content, structural rearrangements, sequence transfer and sequence evolution rates.

RESULTS

The results demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes in three Rhodiola species, with the former possessing more conserved genome structure but faster evolutionary rates of sequence, while the latter exhibiting structural diversity but slower rates of sequence evolution. Some lineage-specific features were observed in Rhodiola mitogenomes, including chromosome fission, gene loss and structural rearrangement. Repeat element analysis shows that the repeats occurring between the two chromosomes may mediate the formation of multichromosomal structure in the mitogenomes of Rhodiola, and this multichromosomal structure may have recently formed. The identification of homologous sequences between plastomes and mitogenomes reveals several unidirectional protein-coding gene transfer events from chloroplasts to mitochondria. Moreover, we found that their organelle genomes contained multiple fragments of nuclear transposable elements (TEs) and exhibited different preferences for TEs insertion type. Genome-wide scans of positive selection identified one gene matR from the mitogenome. Since the matR is crucial for plant growth and development, as well as for respiration and stress responses, our findings suggest that matR may participate in the adaptive response of Rhodiola species to environmental stress of QTP.

CONCLUSION

The study analyzed the organelle genomes of three Rhodiola species and demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes. Signals of positive selection were detected in the matR gene of Rhodiola mitogenomes, suggesting the potential role of this gene in Rhodiola adaptation to QTP. Together, the study is expected to enrich the genomic resources and provide valuable insights into the structural dynamics and sequence divergences of Rhodiola species.

Characterization of the complete chloroplast genome of Amsonia elliptica (Apocynaceae)

Amsonia elliptica (Apocynaceae), endangered species in Korea, is a perennial herb that is economically important as traditional medicine and used as ornamentals. Natural populations of this species are facing extinction due to small population size and isolated distribution. Here, we report the complete chloroplast (cp) genome of A. elliptica using Illumina HiSeq sequencing and its phylogenetic position in subfamily Rauvolfioideae based on 20 Apocynaceae cp genomes. The cp genome of A. elliptica was 154,242 bp in length with a pair of inverted repeats of 25,711 bp, separated by large single-copy and small single-copy regions of 85,382 bp and 17,438 bp, respectively. Our phylogenomic analyses revealed that A. elliptica was closely related to Rhazya stricta in Rauvolfioideae (Apocynaceae).

Complete mitochondrial genome assembly and comparison of Camellia sinensis var. Assamica cv. Duntsa

Camellia sinensis var. Assamica cv. Duntsa (C.duntsa), a valuable Theaceae from Hunan Province, has been looked at as a precious tea resource by local farmers because of its economic and ecological value. Genomics study on C.duntsa is essential for the domestication and enhancement of tea tree varieties. In the present study, we used a hybrid approach based on Illumina and PacBio data to sequence and assemble the mitochondrial genome of C.duntsa. The mitochondrial genome of C.duntsa was estimated to be 1,081,996 base pairs (bp) and eighty-one genes consisting of one pseudogene, three ribosomal RNA (rRNA) genes, thirty transfer RNA (tRNA) genes, and forty-seven protein-coding genes (PCGs). Tetramer repetitions made up 43.90% of simple sequence repeats (SSRs). The codon usage bias of the Theaceae mitochondrial gene atp9 was altered by mutation, but the codon usage of other genes was shaped by natural selection. Besides, there are eighteen gene-containing homologous regions between the chloroplast and mitochondrial genomes of C. duntsa.Some genomes including atp8, cox1, cox3, nad7, nad9, rpl16, rpl2, rps19, rps4, and sdh4 are absent in the mitochondrial genome of several Theaceae plant. However, C. duntsa maintains these genes integrity and functionality. Another gene, rps16, is either lacking from the mitochondrial genome of C. duntsa or is present as a pseudogene. C. duntsa and C. sinensis (OM809792) are very similar, as shown by a collinear match across four species of Theaceae; the most conservative genes are nad5, atp9, cox2, rps3, trnA-TGC, trnI-GAT, rrn18, trnV-GAC, and ccmFN. Similarly, the genome's phylogenetic trees revealed that C. duntsa was the sister species to C. sinensis. The results confirmed that the C. duntsa and C. sinensis (OM809792) mitochondrial genome underwent gene rearrangement.In general, our results shows that genomic information from organelles can help us understand plant phylogeny and can also be used to make molecular markers and study how genetic traits change over time. Our research will contribute to the population genetics and evolution of tea plant.

The complete mitochondrial genome of Camellia nitidissima (Theaceae)

The mitochondrial genome of Camellia nitidissima was sequenced by Illumina and Pacbio sequencing. The results of sequences showed that a total length was 949,915 bp, and the GC content was 45.7% in assembled mitochondrial genome of C. nitidissima. 71 unigenes had been found, including 36 coding proteins and 35 non-coding proteins. Subsequently, the phylogenetic tree was built on 24 plants with the maximum-likelihood method, which had high bootstrap value and fited to the angiosperm phylogeny group classification (APG IV). The study's findings unravel the taxonomic status of C. nitidissima and benefit the evolution study.

Comparative analysis of mitochondrial genomes of two alpine medicinal plants of Gentiana (Gentianaceae)

Gentiana crassicaulis and G. straminea are alpine plants of Gentiana with important medicinal value and complex genetic backgrounds. In this study, the mitochondrial genomes (mtDNAs) of these two species were sequenced. The mtDNAs of G. crassicaulis and G. straminea are 368,808 and 410,086 bp long, respectively, 52 and 49 unique genes are annotated in the two species, and the gene arrangement varies widely. Compared to G. crassicaulis, G. straminea loses three effective genes, namely atp6, trnG-GCC and trnV-GAC. As a pseudogene, the atp6 gene of G. straminea is incomplete, which is rare in higher plants. We detected 1696 and 1858 pairs of long repeats and 213 SSRs and 250 SSs in the mtDNAs of G. crassicaulis and G. straminea, respectively. There are 392 SNPs and 18 InDels between the two genomes, and syntenic sequence and structural variation analysis show low collinearity between the two genomes. Chloroplast DNA transferring to mtDNA is observed in both species, and 46,511 and 55,043 bp transferred segments containing three tRNA genes are identified, respectively. Comparative analysis of mtDNAs of G. crassicaulis, G. straminea and four species of Gentianales determined 18 core genes, and there is no specific gene in G. crassicaulis and G. straminea. The phylogenetic tree based on mtDNAs places Gentianaceae in a branch of Gentianales. This study is the first to analyze the mtDNAs of Gentianaceae, which could provide information for analysis of the structure of mtDNAs of higher plants and phylogenetic research of Gentianaceae and Gentianales.

Plastid phylogenomics and plastome evolution in the morning glory family (Convolvulaceae)

Convolvulaceae, the morning glories or bindweeds, is a large family containing species of economic value, including crops, traditional medicines, ornamentals, and vegetables. However, not only are the phylogenetic relationships within this group still debated at the intertribal and intergeneric levels, but also plastid genome (plastome) complexity within Convolvulaceae is not well surveyed. We gathered 78 plastomes representing 17 genera across nine of the 12 Convolvulaceae tribes. Our plastid phylogenomic trees confirm the monophyly of Convolvulaceae, place the genus Jacquemontia within the subfamily Dicranostyloideae, and suggest that the tribe Merremieae is paraphyletic. In contrast, positions of the two genera Cuscuta and Erycibe are uncertain as the bootstrap support of the branches leading to them is moderate to weak. We show that nucleotide substitution rates are extremely variable among Convolvulaceae taxa and likely responsible for the topological uncertainty. Numerous plastomic rearrangements are detected in Convolvulaceae, including inversions, duplications, contraction and expansion of inverted repeats (IRs), and losses of genes and introns. Moreover, integrated foreign DNA of mitochondrial origin was found in the Jacquemontia plastome, adding a rare example of gene transfer from mitochondria to plastids in angiosperms. In the IR of Dichondra, we discovered an extra copy of rpl16 containing a direct repeat of ca. 200 bp long. This repeat was experimentally demonstrated to trigger effective homologous recombination, resulting in the coexistence of intron-containing and -lacking rpl16 duplicates. Therefore, we propose a hypothetical model to interpret intron loss accompanied by invasion of direct repeats at appropriate positions. Our model complements the intron loss model driven by retroprocessing when genes have lost introns but contain abundant RNA editing sites adjacent to former splicing sites.

Characterization, comparative phylogenetic, and gene transfer analyses of organelle genomes of Rhododendron × pulchrum

Rhododendron × pulchrum, an important horticultural species, is widely distributed in Europe, Asia, and North America. To analyze the phylogenetic and organelle genome information of R. × pulchrum and its related species, the organelle genome of R. × pulchrum was sequenced and assembled. The complete mitochondrial genome showed lineage DNA molecules, which were 816,410 bp long and contained 64 genes, namely 24 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 37 protein-coding genes. The chloroplast genome of R. × pulchrum was reassembled and re-annotated; the results were different from those of previous studies. There were 42 and 46 simple sequence repeats (SSR) identified from the mitochondrial and chloroplast genomes of R. × pulchrum, respectively. Five genes (nad1, nad2, nad4, nad7, and rps3) were potentially useful molecular markers. The R. × pulchrum mitochondrial genome collinear alignment among five species of the Ericaceae showed that the mitochondrial genomes of these related species have a high degree of homology with R. × pulchrum in this gene region, and the most conservative genes were trnC-GCA, trnD-GUC, trnM-CAU, trnN-GUU, trnY-GUA, atp4, nad4, nad2, nad5, ccmC, and rrn26. The phylogenetic trees of mitochondrial genome showed that R. simsii was a sister to R. × pulchrum. The results verified that there was gene rearrangement between R. × pulchrum and R. simsii mitochondrial genomes. The codon usage bias of 10 Ericaceae mitochondrial genes and 7 Rhododendron chloroplast genes were influenced by mutation, while other genes codon usages had undergone selection. The study identified 13 homologous fragments containing gene sequences between the chloroplast and mitochondrial genomes of R. × pulchrum. Overall, our results illustrate the organelle genome information could explain the phylogenetics of plants and could be used to develop molecular markers and genetic evolution. Our study will facilitate the study of population genetics and evolution in Rhododendron and other genera in Ericaceae.

Evaluation of cytotoxicity and antiviral activity of Rhazya stricta Decne leaves extract against influenza A/PR/8/34 (H1N1)

Influenza viruses have developed resistance to the current classes of drugs, which means they could eventually become more virulent and cause more mortality and hospitalization. Our study aims to investigate the antiviral activity of Rhazya stricta Decne leaves extract in vitro and search for new promising drugs from R. stricta identified compounds in silico. The study was performed in vitro by utilizing Madin-Darby Canine Kidney cell line (MDCK) as a substrate for the influenza virus and estimating the inhibition performance of the plant leaves extract. Additionally, in silico screening was conducted to explore the antiviral activity of R. stricta phytochemicals. We investigated the cytotoxicity of R. stricta leaves extract and its antiviral activity against influenza virus (A/Puerto Rico/8/34 (H1N1)) using the MTT assay. The mode of action of the plant leaves extract during the viral life cycle was tested using time-of-addition assay. In silico analyses were performed, including molecular docking, drug-likeness analysis, and toxicity risk assessment, to state the leading compounds to be developed into an anti-influenza virus drug. The 50% cytotoxicity concentration of the leaves extract was CC50: 184.6 µg/mL, and the 50% inhibition concentration was CI50: 19.71 µg\mL. The time of addition assay revealed that R. stricta leaves extract exerted its activity in the late step of the influenza virus replication cycle. In comparison to Oseltamivir, the leading compounds showed better binding affinity and can be developed into oral drugs with low toxicity risk. Isolation and purification of the leading compounds and testing their antiviral activity in vitro and in vivo are required.

Graph-based models of the Oenothera mitochondrial genome capture the enormous complexity of higher plant mitochondrial DNA organization

Plant mitochondrial genomes display an enormous structural complexity, as recombining repeat-pairs lead to the generation of various sub-genomic molecules, rendering these genomes extremely challenging to assemble. We present a novel bioinformatic data-processing pipeline called SAGBAC (Semi-Automated Graph-Based Assembly Curator) that identifies recombinogenic repeat-pairs and reconstructs plant mitochondrial genomes. SAGBAC processes assembly outputs and applies our novel ISEIS (Iterative Sequence Ends Identity Search) algorithm to obtain a graph-based visualization. We applied this approach to three mitochondrial genomes of evening primrose (Oenothera), a plant genus used for cytoplasmic genetics studies. All identified repeat pairs were found to be flanked by two alternative and unique sequence-contigs defining so-called 'double forks', resulting in four possible contig-repeat-contig combinations for each repeat pair. Based on the inferred structural models, the stoichiometry of the different contig-repeat-contig combinations was analyzed using Illumina mate-pair and PacBio RSII data. This uncovered a remarkable structural diversity of the three closely related mitochondrial genomes, as well as substantial phylogenetic variation of the underlying repeats. Our model allows predicting all recombination events and, thus, all possible sub-genomes. In future work, the proposed methodology may prove useful for the investigation of the sub-genome organization and dynamics in different tissues and at various developmental stages.

The chloroplast genome sequence and phylogenetic analysis of Apocynum venetum L

Apocynum venetum L. (Apocynaceae) is valuable for its medicinal compounds and fiber content. Native A. venetum populations are threatened and require protection. Wild A. venetum resources are limited relative to market demand and a poor understanding of the composition of A. venetum at the molecular level. The chloroplast genome contains genetic markers for phylogenetic analysis, genetic diversity evaluation, and molecular identification. In this study, the entire genome of the A. venetum chloroplast was sequenced and analyzed. The A. venetum cp genome is 150,878 bp, with a pair of inverted repeat regions (IRA and IRB). Each inverted repeat region is 25,810 bp, which consist of large (LSC, 81,951 bp) and small (SSC, 17,307 bp) single copy areas. The genome-wide GC content was 38.35%, LSC made up 36.49%, SSC made up 32.41%, and IR made up 43.3%. The A. venetum chloroplast genome encodes 131 genes, including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. This study identified the unique characteristics of the A. venetum chloroplast genome, which will help formulate effective conservation and management strategies as well as molecular identification approaches for this important medicinal plant.

A Review of Rhazya stricta Decne Phytochemistry, Bioactivities, Pharmacological Activities, Toxicity, and Folkloric Medicinal Uses

The local medicinal plant Rhazya stricta Decne is reviewed for its folkloric medicinal, phytochemical, pharmacological, biological, and toxicological features. R. stricta has been used widely in different cultures for various medical disorders. The phytochemical studies performed on the R. stricta extract revealed many alkaloidal and fatty acid compounds. Moreover, several flavonoid and terpenoid compounds were also detected. Pharmacological activates of R. stricta extracts are approved to possess antimicrobial, antioxidant, anticancer, antidiabetic, and antihypertensive activities. Additionally, R. stricta extract was found to hold biological activates such as larvicidal and phytoremediation activates R. stricta extract was found to be toxic, genotoxic, and mutagenic. R. stricta contains novel phytochemical compounds that have not been investigated pharmacologically. Further research is needed through in vitro and in vivo experiments to pave the road for these compounds for medical, veterinary, and ecological uses.

Targeted designing functional markers revealed the role of retrotransposon derived miRNAs as mobile epigenetic regulators in adaptation responses of pistachio

We developed novel miRNA-based markers based on salt responsive miRNA sequences to detect polymorphisms in miRNA sequences and locations. The validation of 76 combined miRNA + miRNA and miRNA + ISSR markers in the three extreme pistachio populations led to the identification of three selected markers that could link salt tolerance phenotype to genotype and divided pistachio genotypes and Pistacia species into three clusters. This novel functional marker system, in addition to more efficient performance, has higher polymorphisms than previous miRNA-based marker systems. The functional importance of the target gene of five miRNAs in the structure of the three selected markers in regulation of different genes such as ECA2, ALA10, PFK, PHT1;4, PTR3, KUP2, GRAS, TCP, bHLH, PHD finger, PLATZ and genes involved in developmental, signaling and biosynthetic processes shows that the polymorphism associated with these selected miRNAs can make a significant phenotypic difference between salt sensitive and tolerant pistachio genotypes. The sequencing results of selected bands showed the presence of conserved miRNAs in the structure of the mitochondrial genome. Further notable findings of this study are that the sequences of PCR products of two selected markers were annotated as Gypsy and Copia retrotransposable elements. The transposition of retrotransposons with related miRNAs by increasing the number of miRNA copies and changing their location between nuclear and organellar genomes can affect the regulatory activity of these molecules. These findings show the crucial role of retrotransposon-derived miRNAs as mobile epigenetic regulators between intracellular genomes in regulating salt stress responses as well as creating new and tolerant phenotypes for adaptation to environmental conditions.

The chicken or the egg? Plastome evolution and an independent loss of the inverted repeat in papilionoid legumes

The plastid genome (plastome), while surprisingly constant in gene order and content across most photosynthetic angiosperms, exhibits variability in several unrelated lineages. During the diversification history of the legume family Fabaceae, plastomes have undergone many rearrangements, including inversions, expansion, contraction and loss of the typical inverted repeat (IR), gene loss and repeat accumulation in both shared and independent events. While legume plastomes have been the subject of study for some time, most work has focused on agricultural species in the IR-lacking clade (IRLC) and the plant model Medicago truncatula. The subfamily Papilionoideae, which contains virtually all of the agricultural legume species, also comprises most of the plastome variation detected thus far in the family. In this study three non-papilioniods were included among 34 newly sequenced legume plastomes, along with 33 publicly available sequences, to assess plastome structural evolution in the subfamily. In an effort to examine plastome variation across the subfamily, approximately 20% of the sampling represents the IRLC with the remainder selected to represent the early-branching papilionoid clades. A number of IR-related and repeat-mediated changes were identified and examined in a phylogenetic context. Recombination between direct repeats associated with ycf2 resulted in intraindividual plastome heteroplasmy. Although loss of the IR has not been reported in legumes outside of the IRLC, one genistoid taxon was found to completely lack the typical plastome IR. The role of the IR and non-IR repeats in the progression of plastome change is discussed.

Both Conifer II and Gnetales are characterized by a high frequency of ancient mitochondrial gene transfer to the nuclear genome

BACKGROUND

Mitochondrial gene transfer/loss is common in land plants, and therefore the fate of missing mitochondrial genes has attracted more and more attention. The gene content of gymnosperm mitochondria varies greatly, supplying a system for studying the evolutionary fate of missing mitochondrial genes.

RESULTS

Here, we studied the tempo and pattern of mitochondrial gene transfer/loss in gymnosperms represented by all 13 families, using high-throughput sequencing of both DNA and cDNA. All 41 mitochondrial protein-coding genes were found in cycads, Ginkgo and Pinaceae, whereas multiple mitochondrial genes were absent in Conifer II and Gnetales. In Conifer II, gene transfer from mitochondria to the nucleus followed by loss of the mitochondrial copy was common, but complete loss of a gene in both mitochondrial and nuclear genomes was rare. In contrast, both gene transfer and loss were commonly found in Gnetales. Notably, in Conifer II and Gnetales, the same five mitochondrial genes were transferred to the nuclear genome, and these gene transfer events occurred, respectively, in ancestors of the two lineages. A two-step transfer mechanism (retroprocessing and subsequent DNA-mediated gene transfer) may be responsible for mitochondrial gene transfer in Conifer II and Gnetales. Moreover, the mitochondrial gene content variation is correlated with gene length, GC content, hydrophobicity, and nucleotide substitution rates in land plants.

CONCLUSIONS

This study reveals a complete evolutionary scenario for variations of mitochondrial gene transferring in gymnosperms, and the factors responsible for mitochondrial gene content variation in land plants.

Plastome evolution and organisation in the Hoya group (Apocynaceae)

The genus Hoya is highly diverse and many of its species are popular ornamental plants. However, the relationships between Hoya and related genera (the Hoya group) are not fully resolved. In this study, we report 20 newly sequenced plastomes of species in the Hoya group. The complete plastomes vary in length from 175,405 to 178,525 bp while the LSCs vary from 90,248 to 92,364 bp and the complete SSCs vary from 2,285 to 2,304 bp, making the SSC in the Hoya group one of the shortest known in the angiosperms. The plastome structure in the Hoya group is characterised by a massive increase in the size of the inverted repeats as compared to the outgroups. In all ingroup species, the IR/SSC boundary moved from ycf1 to ndhF while this was not observed in outgroup taxa, making it a synapomorphy for the Hoya group. We have also assembled the mitogenome of Hoya lithophytica, which, at 718,734 bp, is the longest reported in the family. The phylogenetic analysis using exons from 42 taxa in the Hoya group and three outgoups confirms that the earliest divergent genus in the Hoya group is Papuahoya, followed by Dischidia. The relationship between Dischidia and the clade which includes all Hoya and Oreosparte taxa, is not fully supported. Oreosparte is nested in Hoya making it paraphyletic unless Clemensiella is recognised as a separate genus.

Comparative Mitogenomic Analysis Reveals Gene and Intron Dynamics in Rubiaceae and Intra-Specific Diversification in Damnacanthus indicus

The dynamic evolution of mitochondrial gene and intron content has been reported across the angiosperms. However, a reference mitochondrial genome (mitogenome) is not available in Rubiaceae. The phylogenetic utility of mitogenome data at a species level is rarely assessed. Here, we assembled mitogenomes of six Damnacanthus indicus (Rubiaceae, Rubioideae) representing two varieties (var. indicus and var. microphyllus). The gene and intron content of D. indicus was compared with mitogenomes from representative angiosperm species and mitochondrial contigs from the other Rubiaceae species. Mitogenome structural rearrangement and sequence divergence in D. indicus were analyzed in six individuals. The size of the mitogenome in D. indicus varied from 417,661 to 419,435 bp. Comparing the number of intact mitochondrial protein-coding genes in other Gentianales taxa (38), D. indicus included 32 genes representing several losses. The intron analysis revealed a shift from cis to trans splicing of a nad1 intron (nad1i728) in D. indicus and it is a shared character with the other four Rubioideae taxa. Two distinct mitogenome structures (type A and B) were identified. Two-step direct repeat-mediated recombination was proposed to explain structural changes between type A and B mitogenomes. The five individuals from two varieties in D. indicus diverged well in the whole mitogenome-level comparison with one exception. Collectively, our study elucidated the mitogenome evolution in Rubiaceae along with D. indicus and showed the reliable phylogenetic utility of the whole mitogenome data at a species-level evolution.

Comparative analysis of the complete plastid genomes of Mangifera species and gene transfer between plastid and mitochondrial genomes

Mango is an important commercial fruit crop belonging to the genus Mangifera. In this study, we reported and compared four newly sequenced plastid genomes of the genus Mangifera, which showed high similarities in overall size (157,780–157,853 bp), genome structure, gene order, and gene content. Three mutation hotspots (trnG-psbZ, psbD-trnT, and ycf4-cemA) were identified as candidate DNA barcodes for Mangifera. These three DNA barcode candidate sequences have high species identification ability. We also identified 12 large fragments that were transferred from the plastid genome to the mitochondrial genome, and found that the similarity was more than 99%. The total size of the transferred fragment was 35,652 bp, accounting for 22.6% of the plastid genome. Fifteen intact chloroplast genes, four tRNAs and numerous partial genes and intergenic spacer regions were identified. There are many of these genes transferred from mitochondria to the chloroplast in other species genomes. Phylogenetic analysis based on whole plastid genome data provided a high support value, and the interspecies relationships within Mangifera were resolved well.

Characterization of Chloroplast Genomes From Two Salvia Medicinal Plants and Gene Transfer Among Their Mitochondrial and Chloroplast Genomes

Salvia species have been widely used as medicinal plants and have played an important role in the treatment and recovery of individuals with COVID-19. In this study, we reported two newly identified whole chloroplast genome sequences of Salvia medicinal plants (Salvia yangii and Salvia miltiorrhiza f. alba) and compared them with those of seven other reported Salvia chloroplast genomes. These were proven to be highly similar in terms of overall size, genome structure, gene content, and gene order. We identified 10 mutation hot spots (trnK-rps16, atpH-atpI, psaA-ycf3, ndhC-trnV, ndhF, rpl32-trnL, ndhG-ndhI, rps15-ycf1, ycf1a, and ycf1b) as candidate DNA barcodes for Salvia. Additionally, we observed the transfer of nine large-sized chloroplast genome fragments, with a total size of 49,895 bp (accounting for 32.97% of the chloroplast genome), into the mitochondrial genome as they shared >97% sequence similarity. Phylogenetic analyses of the whole chloroplast genome provided a high resolution of Salvia. This study will pave the way for the identification and breeding of Salvia medicinal plants and further phylogenetic evolutionary research on them as well.

Deciphering the organelle genomes and transcriptomes of a common ornamental plant Ligustrum quihoui reveals multiple fragments of transposable elements in the mitogenome

Ligustrum quihoui (L. quihoui) is an important hedge material for landscaping and also possesses medicinal value. To generate genomic resources for better understanding the evolutionary history of this important plant, the organelle genomes of L. quihoui are de novo assembled and functionally annotated. Compared with other Oleaceae species, the 163,069 bp chloroplast genome of L. quihoui exhibits a typical quadripartite structure with highly conserved gene content and gene order, while the 848,451 bp mitochondrial genome of L. quihoui exhibits highly divergent genome size and gene content. Codon usage analyses show that genes related with photosynthesis and mitochondrial respiratory chain show inconsistent codon biases. A total of 48,760 bp transposable elements (TEs) fragments and 41,887 bp chloroplast-like sequences are found in the L. quihoui mitochondrial genome. A striking discrepancy of RNA editing between the two organelle genomes is found in L. quihoui, in which 146 mitochondrial editing sites coexist with only 43 such sites in chloroplast. Based on DNA and RNA-Seq data, we propose that GTG may act as the start codon of mitochondrial rpl16 in Oleaceae species. Phylogenetic analysis based on chloroplast genome shows that L. quihoui and L. japonicum form a sister clade within the genus Ligustrum.

Extensive variation in nucleotide substitution rate and gene/intron loss in mitochondrial genomes of Pelargonium

Geraniaceae organelle genomes have been shown to exhibit several highly unusual features compared to most other photosynthetic angiosperms. This includes massively rearranged plastomes with considerable size variation, extensive gene and intron loss, accelerated rates of nucleotide substitutions in both mitogenomes and plastomes, and biparental inheritance and cytonuclear incompatibility of the plastome. Most previous studies have focused on plastome evolution with mitogenome comparisons limited to only a few taxa or genes. In this study, mitogenomes and transcriptomes were examined for 27 species of Geraniales, including 13 species of Pelargonium. Extensive gene and intron losses were detected across the Geraniales with Pelargonium representing the most gene depauperate lineage in the family. Plotting these events on the Geraniaceae phylogenetic tree showed that gene losses occurred multiple times, whereas intron losses more closely reflected the relationships among taxa. In addition, P. australe acquired an intron by horizontal transfer. Comparisons of nucleotide substitution rates in Pelargonium showed that synonymous changes in nuclear genes were much lower than in mitochondrial genes. This is in contrast to the previously published studies that indicated that nuclear genes have 16 fold higher rates than mitochondrial genes across angiosperms. Elevated synonymous substitutions occurred for each mitochondrial gene in Pelargonium with the highest values 783 and 324 times higher than outgroups and other Geraniaceae, respectively. Pelargonium is one of four unrelated genera of angiosperms (Ajuga, Plantago and Silene) that have experienced highly accelerated nucleotide substitutions in mitogenomes. It is distinct from most angiosperms in also having elevated substitution rates in plastid genes but the cause of rate accelerations in Pelargonium plastomes and mitogenomes may be different.

The Mitogenome of Norway Spruce and a Reappraisal of Mitochondrial Recombination in Plants

Plant mitogenomes can be difficult to assemble because they are structurally dynamic and prone to intergenomic DNA transfers, leading to the unusual situation where an organelle genome is far outnumbered by its nuclear counterparts. As a result, comparative mitogenome studies are in their infancy and some key aspects of genome evolution are still known mainly from pregenomic, qualitative methods. To help address these limitations, we combined machine learning and in silico enrichment of mitochondrial-like long reads to assemble the bacterial-sized mitogenome of Norway spruce (Pinaceae: Picea abies). We conducted comparative analyses of repeat abundance, intergenomic transfers, substitution and rearrangement rates, and estimated repeat-by-repeat homologous recombination rates. Prompted by our discovery of highly recombinogenic small repeats in P. abies, we assessed the genomic support for the prevailing hypothesis that intramolecular recombination is predominantly driven by repeat length, with larger repeats facilitating DNA exchange more readily. Overall, we found mixed support for this view: Recombination dynamics were heterogeneous across vascular plants and highly active small repeats (ca. 200 bp) were present in about one-third of studied mitogenomes. As in previous studies, we did not observe any robust relationships among commonly studied genome attributes, but we identify variation in recombination rates as a underinvestigated source of plant mitogenome diversity.

The Mitochondrial Genome Assembly of Fennel (Foeniculum vulgare) Reveals Two Different atp6 Gene Sequences in Cytoplasmic Male Sterile Accessions

Cytoplasmic male sterility (CMS) has always aroused interest among researchers and breeders, being a valuable resource widely exploited not only to breed F1 hybrid varieties but also to investigate genes that control stamen and pollen development. With the aim of identifying candidate genes for CMS in fennel, we adopted an effective strategy relying on the comparison between mitochondrial genomes (mtDNA) of both fertile and sterile genotypes. mtDNA raw reads derived from a CMS genotype were assembled in a single molecule (296,483 bp), while a draft mtDNA assembly (166,124 nucleotides, 94 contigs) was performed using male fertile sample (MF) sequences. From their annotation and alignment, two atp6-like sequences were identified. atp6-, the putative mutant copy with a 300 bp truncation at the 5'-end, was found only in the mtDNA of CMS samples, while the wild type copy (atp6+) was detected only in the MF mtDNA. Further analyses (i.e., reads mapping and Sanger sequencing), revealed an atp6+ copy also in CMS samples, probably in the nuclear DNA. However, qPCRs performed on different tissues proved that, despite its availability, atp6+ is expressed only in MF samples, while apt6- mRNA was always detected in CMS individuals. In the light of these findings, the energy deficiency model could explain the pollen deficiency observed in male sterile flower. atp6- could represent a gene whose mRNA is translated into a not-fully functional protein leading to suboptimal ATP production that guarantees essential cellular processes but not a high energy demand process such as pollen development. Our study provides novel insights into the fennel mtDNA genome and its atp6 genes, and paves the way for further studies aimed at understanding their functional roles in the determination of male sterility.

Mitochondrial plastid DNA can cause DNA barcoding paradox in plants

The transfer of ancestral plastid genomes into mitochondrial genomes to generate mitochondrial plastid DNA (MTPT) is known to occur in plants, but its impacts on mitochondrial genome complexity and the potential for causing a false-positive DNA barcoding paradox have been underestimated. Here, we assembled the organelle genomes of Cynanchum wilfordii and C. auriculatum, which are indigenous medicinal herbs in Korea and China, respectively. In both species, it is estimated that 35% of the ancestral plastid genomes were transferred to mitochondrial genomes over the past 10 million years and remain conserved in these genomes. Some plastid barcoding markers co-amplified the conserved MTPTs and caused a barcoding paradox, resulting in mis-authentication of botanical ingredients and/or taxonomic mis-positioning. We identified dynamic and lineage-specific MTPTs that have contributed to mitochondrial genome complexity and might cause a putative barcoding paradox across 81 plant species. We suggest that a DNA barcoding guidelines should be developed involving the use of multiple markers to help regulate economically motivated adulteration.

Comprehensive comparative analysis of chloroplast genomes from seven Panax species and development of an authentication system based on species-unique single nucleotide polymorphism markers

BACKGROUND

Panax species are important herbal medicinal plants in the Araliaceae family. Recently, we reported the complete chloroplast genomes and 45S nuclear ribosomal DNA sequences from seven Panax species, two (P . quinqu e folius and P . trifolius) from North America and five (P . ginseng, P . notoginseng, P . japonicus, P . vietnamensis, and P . stipuleanatus) from Asia.

METHODS

We conducted phylogenetic analysis of these chloroplast sequences with 12 other Araliaceae species and comprehensive comparative analysis among the seven Panax whole chloroplast genomes.

RESULTS

We identified 1,128 single nucleotide polymorphisms (SNP) in coding gene sequences, distributed among 72 of the 79 protein-coding genes in the chloroplast genomes of the seven Panax species. The other seven genes (including psaJ, psbN, rpl23, psbF, psbL, rps18, and rps7) were identical among the Panax species. We also discovered that 12 large chloroplast genome fragments were transferred into the mitochondrial genome based on sharing of more than 90% sequence similarity. The total size of transferred fragments was 60,331 bp, corresponding to approximately 38.6% of chloroplast genome. We developed 18 SNP markers from the chloroplast genic coding sequence regions that were not similar to regions in the mitochondrial genome. These markers included two or three species-specific markers for each species and can be used to authenticate all the seven Panax species from the others.

CONCLUSION

The comparative analysis of chloroplast genomes from seven Panax species elucidated their genetic diversity and evolutionary relationships, and 18 species-specific markers were able to discriminate among these species, thereby furthering efforts to protect the ginseng industry from economically motivated adulteration.

Fluctuations in Fabaceae mitochondrial genome size and content are both ancient and recent

BACKGROUND

Organelle genome studies of Fabaceae, an economically and ecologically important plant family, have been biased towards the plastid genome (plastome). Thus far, less than 15 mitochondrial genome (mitogenome) sequences of Fabaceae have been published, all but four of which belong to the subfamily Papilionoideae, limiting the understanding of size variation and content across the family. To address this, four mitogenomes were sequenced and assembled from three different subfamilies (Cercidoideae, Detarioideae and Caesalpinioideae).

RESULTS

Phylogenetic analysis based on shared mitochondrial protein coding regions produced a fully resolved and well-supported phylogeny that was completely congruent with the plastome tree. Comparative analyses suggest that two kinds of mitogenome expansions have occurred in Fabaceae. Size expansion of four genera (Tamarindus, Libidibia, Haematoxylum, and Leucaena) in two subfamilies (Detarioideae and Caesalpinioideae) occurred in relatively deep nodes, and was mainly caused by intercellular gene transfer and/or interspecific horizontal gene transfer (HGT). The second, more recent expansion occurred in the Papilionoideae as a result of duplication of native mitochondrial sequences. Family-wide gene content analysis revealed 11 gene losses, four (rps2, 7, 11 and 13) of which occurred in the ancestor of Fabaceae. Losses of the remaining seven genes (cox2, rpl2, rpl10, rps1, rps19, sdh3, sdh4) were restricted to specific lineages or occurred independently in different clades. Introns of three genes (cox2, ccmFc and rps10) showed extensive lineage-specific length variation due to large sequence insertions and deletions. Shared DNA analysis among Fabaceae mitogenomes demonstrated a substantial decay of intergenic spacers and provided further insight into HGT between the mimosoid clade of Caesalpinioideae and the holoparasitic Lophophytum (Balanophoraceae).

CONCLUSION

This study represents the most exhaustive analysis of Fabaceae mitogenomes so far, and extends the understanding the dynamic variation in size and gene/intron content. The four newly sequenced mitogenomes reported here expands the phylogenetic coverage to four subfamilies. The family has experienced multiple mitogenome size fluctuations in both ancient and recent times. The causes of these size variations are distinct in different lineages. Fabaceae mitogenomes experienced extensive size fluctuation by recruitment of exogenous DNA and duplication of native mitochondrial DNA.

A draft genome and transcriptome of common milkweed (Asclepias syriaca) as resources for evolutionary, ecological, and molecular studies in milkweeds and Apocynaceae

Milkweeds (Asclepias) are used in wide-ranging studies including floral development, pollination biology, plant-insect interactions and co-evolution, secondary metabolite chemistry, and rapid diversification. We present a transcriptome and draft nuclear genome assembly of the common milkweed, Asclepias syriaca. This reconstruction of the nuclear genome is augmented by linkage group information, adding to existing chloroplast and mitochondrial genomic resources for this member of the Apocynaceae subfamily Asclepiadoideae. The genome was sequenced to 80.4× depth and the draft assembly contains 54,266 scaffolds ≥1 kbp, with N50 = 3,415 bp, representing 37% (156.6 Mbp) of the estimated 420 Mbp genome. A total of 14,474 protein-coding genes were identified based on transcript evidence, closely related proteins, and ab initio models, and 95% of genes were annotated. A large proportion of gene space is represented in the assembly, with 96.7% of Asclepias transcripts, 88.4% of transcripts from the related genus Calotropis, and 90.6% of proteins from Coffea mapping to the assembly. Scaffolds covering 75 Mbp of the Asclepias assembly formed 11 linkage groups. Comparisons of these groups with pseudochromosomes in Coffea found that six chromosomes show consistent stability in gene content, while one may have a long history of fragmentation and rearrangement. The progesterone 5β-reductase gene family, a key component of cardenolide production, is likely reduced in Asclepias relative to other Apocynaceae. The genome and transcriptome of common milkweed provide a rich resource for future studies of the ecology and evolution of a charismatic plant family.

Prospects of genomic resources available at the global databases for the flora of United Arab Emirates

This article emphasizes available genomic resources at the global databases National Center for Biotechnology Information (NCBI) GenBank, Gramene and Phytozome for the selected 378 plant taxa of the United Arab Emirates (UAE). Germplasm of these species was collected and banked at the Sharjah Seed Bank and Herbarium (SSBH) along with their related information on habit, habitat and occurrence. The occurrence statistics exhibits almost 19.84% species under rare-to-very rare category, the GenBank search statistics for this category indicates 17.72% species as studied and 2.11% as not studied. Overall, from the global search statistics for 378 plant species, it seems that about 40 (10.58%) species remained unstudied. Most of the unstudied species were herbaceous plants belonging to the mountainous habitat. Moreover, full genomes were recorded for 7 species at NCBI GenBank, 2 species at Phytozome and 1 species at Gramene database. The local search statistics (for UAE) exhibits about 10.58% of the flora that still remained unstudied and only 11 (2.90%) of the recorded species were having genomic information at NCBI GenBank. It is necessary to prioritize studies on such species that could provide valuable insight on their genetic composition in order to understand their adaptation to the natural environment. At present, the SSBH is cataloguing UAE's flora using core barcode and assisted markers that could provide a robust DNA barcode library for native plants of UAE. Our study appeals researchers to recognize and prioritize the species that need attention to enrich their genomic resources at the global databases by supporting nucleotide libraries with their conspecifics. At present, genomic resources for UAE plant taxa are limited, but with the advent of low-cost sequencing technologies these resources would flourish in the near future. Nevertheless, the information generated through genomic studies could be utilized for conservation and management of threatened and endangered plant species, Crop Wild Relatives and medicinal plants. We hope this article will promote interest in conducting additional studies in genomics of desert plants by encouraging researchers to participate in this emerging field.

Evolution and Diversification of Kiwifruit Mitogenomes through Extensive Whole-Genome Rearrangement and Mosaic Loss of Intergenic Sequences in a Highly Variable Region

Angiosperm mitochondrial genomes (mitogenomes) are notable for their extreme diversity in both size and structure. However, our current understanding of this diversity is limited, and the underlying mechanism contributing to this diversity remains unclear. Here, we completely assembled and compared the mitogenomes of three kiwifruit (Actinidia) species, which represent an early divergent lineage in asterids. We found conserved gene content and fewer genomic repeats, particularly large repeats (>1 kb), in the three mitogenomes. However, sequence transfers such as intracellular events are variable and dynamic, in which both ancestral shared and recently species-specific events as well as complicated transfers of two plastid-derived sequences into the nucleus through the mitogenomic bridge were detected. We identified extensive whole-genome rearrangements among kiwifruit mitogenomes and found a highly variable V region in which fragmentation and frequent mosaic loss of intergenic sequences occurred, resulting in greatly interspecific variations. One example is the fragmentation of the V region into two regions, V1 and V2, giving rise to the two mitochondrial chromosomes of Actinidia chinensis. Finally, we compared the kiwifruit mitogenomes with those of other asterids to characterize their overall mitogenomic diversity, which identified frequent gain/loss of genes/introns across lineages. In addition to repeat-mediated recombination and import-driven hypothesis of genome size expansion reported in previous studies, our results highlight a pattern of dynamic structural variation in plant mitogenomes through global genomic rearrangements and species-specific fragmentation and mosaic loss of intergenic sequences in highly variable regions on the basis of a relatively large ancestral mitogenome.

The plastid and mitochondrial genomes of Eucalyptus grandis

BACKGROUND

Land plant organellar genomes have significant impact on metabolism and adaptation, and as such, accurate assembly and annotation of plant organellar genomes is an important tool in understanding the evolutionary history and interactions between these genomes. Intracellular DNA transfer is ongoing between the nuclear and organellar genomes, and can lead to significant genomic variation between, and within, species that impacts downstream analysis of genomes and transcriptomes.

RESULTS

In order to facilitate further studies of cytonuclear interactions in Eucalyptus, we report an updated annotation of the E. grandis plastid genome, and the second sequenced and annotated mitochondrial genome of the Myrtales, that of E. grandis. The 478,813 bp mitochondrial genome shows the conserved protein coding regions and gene order rearrangements typical of land plants. There have been widespread insertions of organellar DNA into the E. grandis nuclear genome, which span 141 annotated nuclear genes. Further, we identify predicted editing sites to allow for the discrimination of RNA-sequencing reads between nuclear and organellar gene copies, finding that nuclear copies of organellar genes are not expressed in E. grandis.

CONCLUSIONS

The implications of organellar DNA transfer to the nucleus are often ignored, despite the insight they can give into the ongoing evolution of plant genomes, and the problems they can cause in many applications of genomics. Future comparisons of the transcription and regulation of organellar genes between Eucalyptus genotypes may provide insight to the cytonuclear interactions that impact economically important traits in this widely grown lignocellulosic crop species.

The Promise of Molecular and Genomic Techniques for Biodiversity Research and DNA Barcoding of the Arabian Peninsula Flora

The Arabian Peninsula is known to have a comprehensive and rich endowment of unique and genetically diverse plant genetic resources. Analysis and conservation of biological diversity is a crucial issue to the whole Arabian Peninsula. The rapid and accurate delimitation and identification of a species is crucial to genetic diversity analysis and the first critical step in the assessment of distribution, population abundance and threats related to a particular target species. During the last two decades, classical strategies of evaluating genetic variability, such as morphology and physiology, have been greatly complemented by phylogenetic, taxonomic, genetic diversity and breeding research molecular studies. At present, initiatives are taking place around the world to generate DNA barcode libraries for vascular plant flora and to make these data available in order to better understand, conserve and utilize biodiversity. The number of herbarium collection-based plant evolutionary genetics and genomics studies being conducted has been increasing worldwide. The herbaria provide a rich resource of already preserved and identified material, and these as well as freshly collected samples from the wild can be used for creating a reference DNA barcode library for the vascular plant flora of a region. This review discusses the main molecular and genomic techniques used in plant identification and biodiversity analysis. Hence, we highlight studies emphasizing various molecular techniques undertaken during the last 10 years to study the plant biodiversity of the Arabian Peninsula. Special emphasis on the role of DNA barcoding as a powerful tool for plant biodiversity analysis is provided, along with the crucial role of herbaria in creating a DNA barcode library.

Trends in herbgenomics

From Shen Nong's Herbal Classic (Shennong Bencao Jing) to the Compendium of Materia Medica (Bencao Gangmu) and the first scientific Nobel Prize for the mainland of China, each milestone in the historical process of the development of traditional Chinese medicine (TCM) involves screening, testing and integrating. After thousands of years of inheritance and development, herbgenomics (bencaogenomics) has bridged the gap between TCM and international advanced omics studies, promoting the application of frontier technologies in TCM. It is a discipline that uncovers the genetic information and regulatory networks of herbs to clarify their molecular mechanism in the prevention and treatment of human diseases. The main theoretical system includes genomics, functional genomics, proteomics, transcriptomics, metabolomics, epigenomics, metagenomics, synthetic biology, pharmacogenomics of TCM, and bioinformatics, among other fields. Herbgenomics is mainly applicable to the study of medicinal model plants, genomic-assisted breeding, herbal synthetic biology, protection and utilization of gene resources, TCM quality evaluation and control, and TCM drug development. Such studies will accelerate the application of cutting-edge technologies, revitalize herbal research, and strongly promote the development and modernization of TCM.

Towards a comprehensive picture of C-to-U RNA editing sites in angiosperm mitochondria

Our understanding of the dynamic and evolution of RNA editing in angiosperms is in part limited by the few editing sites identified to date. This study identified 10,217 editing sites from 17 diverse angiosperms. Our analyses confirmed the universality of certain features of RNA editing, and offer new evidence behind the loss of editing sites in angiosperms. RNA editing is a post-transcriptional process that substitutes cytidines (C) for uridines (U) in organellar transcripts of angiosperms. These substitutions mostly take place in mitochondrial messenger RNAs at specific positions called editing sites. By means of publicly available RNA-seq data, this study identified 10,217 editing sites in mitochondrial protein-coding genes of 17 diverse angiosperms. Even though other types of mismatches were also identified, we did not find evidence of non-canonical editing processes. The results showed an uneven distribution of editing sites among species, genes, and codon positions. The analyses revealed that editing sites were conserved across angiosperms but there were some species-specific sites. Non-synonymous editing sites were particularly highly conserved (~ 80%) across the plant species and were efficiently edited (80% editing extent). In contrast, editing sites at third codon positions were poorly conserved (~ 30%) and only partially edited (~ 40% editing extent). We found that the loss of editing sites along angiosperm evolution is mainly occurring by replacing editing sites with thymidines, instead of a degradation of the editing recognition motif around editing sites. Consecutive and highly conserved editing sites had been replaced by thymidines as result of retroprocessing, by which edited transcripts are reverse transcribed to cDNA and then integrated into the genome by homologous recombination. This phenomenon was more pronounced in eudicots, and in the gene cox1. These results suggest that retroprocessing is a widespread driving force underlying the loss of editing sites in angiosperm mitochondria.

Evolution on the backbone: Apocynaceae phylogenomics and new perspectives on growth forms, flowers, and fruits

PREMISE OF THE STUDY

We provide the largest phylogenetic analyses to date of Apocynaceae in terms of taxa and molecular data as a framework for analyzing the evolution of vegetative and reproductive traits.

METHODS

We produced maximum-likelihood phylogenies of Apocynaceae using 21 plastid loci sampled from 1045 species (nearly 25% of the family) and complete plastomes from 73 species. We reconstructed ancestral states and used model comparisons in a likelihood framework to analyze character evolution across Apocynaceae.

KEY RESULTS

We obtained a well-supported phylogeny of Apocynaceae, resolving poorly understood tribal and subtribal relationships (e.g., among Amsonieae and Hunterieae, within Asclepiadeae), rejecting monophyly of Melodineae and Odontadenieae, and placing previously unsampled and enigmatic taxa (e.g., Pycnobotrya). We provide new insights into the evolution of Apocynaceae, including frequent shifts between herbaceousness and woodiness, reversibility of twining, integrated evolution of the corolla and gynostegium, and ancestral baccate fruits.

CONCLUSIONS

Increased sampling and selection of best-fitting models of evolution provide more resolved and robust estimates of phylogeny and character evolution than obtained in previous studies. Evolutionary inferences are sensitive to choice of phylogenetic frameworks and models.

Authentication Markers for Five Major Panax Species Developed via Comparative Analysis of Complete Chloroplast Genome Sequences

Ginseng represents a set of high-value medicinal plants of different species: Panax ginseng (Asian ginseng), Panax quinquefolius (American ginseng), Panax notoginseng (Chinese ginseng), Panax japonicus (Bamboo ginseng), and Panax vietnamensis (Vietnamese ginseng). Each species is pharmacologically and economically important, with differences in efficacy and price. Accordingly, an authentication system is needed to combat economically motivated adulteration of Panax products. We conducted comparative analysis of the chloroplast genome sequences of these five species, identifying 34-124 InDels and 141-560 SNPs. Fourteen InDel markers were developed to authenticate the Panax species. Among these, eight were species-unique markers that successfully differentiated one species from the others. We generated at least one species-unique marker for each of the five species, and any of the species can be authenticated by selection among these markers. The markers are reliable, easily detectable, and valuable for applications in the ginseng industry as well as in related research.

Optimization and Comparative Analysis of Plant Organellar DNA Enrichment Methods Suitable for Next-generation Sequencing

Plant organellar genomes contain large, repetitive elements that may undergo pairing or recombination to form complex structures and/or sub-genomic fragments. Organellar genomes also exist in admixtures within a given cell or tissue type (heteroplasmy), and an abundance of subtypes may change throughout development or when under stress (sub-stoichiometric shifting). Next-generation sequencing (NGS) technologies are required to obtain deeper understanding of organellar genome structure and function. Traditional sequencing studies use several methods to obtain organellar DNA: (1) If a large amount of starting tissue is used, it is homogenized and subjected to differential centrifugation and/or gradient purification. (2) If a smaller amount of tissue is used (i.e., if seeds, material, or space is limited), the same process is performed as in (1), followed by whole-genome amplification to obtain sufficient DNA. (3) Bioinformatics analysis can be used to sequence the total genomic DNA and to parse out organellar reads. All these methods have inherent challenges and tradeoffs. In (1), it may be difficult to obtain such a large amount of starting tissue; in (2), whole-genome amplification could introduce a sequencing bias; and in (3), homology between nuclear and organellar genomes could interfere with assembly and analysis. In plants with large nuclear genomes, it is advantageous to enrich for organellar DNA to reduce sequencing costs and sequence complexity for bioinformatics analyses. Here, we compare a traditional differential centrifugation method with a fourth method, an adapted CpG-methyl pulldown approach, to separate the total genomic DNA into nuclear and organellar fractions. Both methods yield sufficient DNA for NGS, DNA that is highly enriched for organellar sequences, albeit at different ratios in mitochondria and chloroplasts. We present the optimization of these methods for wheat leaf tissue and discuss major advantages and disadvantages of each approach in the context of sample input, protocol ease, and downstream application.

The mitochondrial genome of the terrestrial carnivorous plant Utricularia reniformis (Lentibulariaceae): Structure, comparative analysis and evolutionary landmarks

The carnivorous plants of the family Lentibulariaceae have attained recent attention not only because of their interesting lifestyle, but also because of their dynamic nuclear genome size. Lentibulariaceae genomes span an order of magnitude and include species with the smallest genomes in angiosperms, making them a powerful system to study the mechanisms of genome expansion and contraction. However, little is known about mitochondrial DNA (mtDNA) sequences of this family, and the evolutionary forces that shape this organellar genome. Here we report the sequencing and assembly of the complete mtDNA from the endemic terrestrial Brazilian species Utricularia reniformis. The 857,234bp master circle mitochondrial genome encodes 70 transcriptionaly active genes (42 protein-coding, 25 tRNAs and 3 rRNAs), covering up to 7% of the mtDNA. A ltrA-like protein related to splicing and mobility and a LAGLIDADG homing endonuclease have been identified in intronic regions, suggesting particular mechanisms of genome maintenance. RNA-seq analysis identified properties with putative diverse and important roles in genome regulation and evolution: 1) 672kbp (78%) of the mtDNA is covered by full-length reads; 2) most of the 243kbp intergenic regions exhibit transcripts; and 3) at least 69 novel RNA editing sites in the protein-coding genes. Additional genomic features are hypothetical ORFs (48%), chloroplast insertions, including truncated plastid genes that have been lost from the chloroplast DNA (5%), repeats (5%), relics of transposable elements mostly related to LTR retrotransposons (5%), and truncated mitovirus sequences (0.4%). Phylogenetic analysis based on 32 different Lamiales mitochondrial genomes corroborate that Lentibulariaceae is a monophyletic group. In summary, the U. reniformis mtDNA represents the eighth largest plant mtDNA described to date, shedding light on the genomic trends and evolutionary characteristics and phylogenetic history of the family Lentibulariaceae.

Contrasting Patterns of Nucleotide Substitution Rates Provide Insight into Dynamic Evolution of Plastid and Mitochondrial Genomes of Geranium

Geraniaceae have emerged as a model system for investigating the causes and consequences of variation in plastid and mitochondrial genomes. Incredible structural variation in plastid genomes (plastomes) and highly accelerated evolutionary rates have been reported in selected lineages and functional groups of genes in both plastomes and mitochondrial genomes (mitogenomes), and these phenomena have been implicated in cytonuclear incompatibility. Previous organelle genome studies have included limited sampling of Geranium, the largest genus in the family with over 400 species. This study reports on rates and patterns of nucleotide substitutions in plastomes and mitogenomes of 17 species of Geranium and representatives of other Geraniaceae. As detected across other angiosperms, substitution rates in the plastome are 3.5 times higher than the mitogenome in most Geranium. However, in the branch leading to Geranium brycei/Geranium incanum mitochondrial genes experienced significantly higher dN and dS than plastid genes, a pattern that has only been detected in one other angiosperm. Furthermore, rate accelerations differ in the two organelle genomes with plastomes having increased dN and mitogenomes with increased dS. In the Geranium phaeum/Geranium reflexum clade, duplicate copies of clpP and rpoA genes that experienced asymmetric rate divergence were detected in the single copy region of the plastome. In the case of rpoA, the branch leading to G. phaeum/G. reflexum experienced positive selection or relaxation of purifying selection. Finally, the evolution of acetyl-CoA carboxylase is unusual in Geraniaceae because it is only the second angiosperm family where both prokaryotic and eukaryotic ACCases functionally coexist in the plastid.

Mitochondrial genome evolution in Alismatales: Size reduction and extensive loss of ribosomal protein genes

The order Alismatales is a hotspot for evolution of plant mitochondrial genomes characterized by remarkable differences in genome size, substitution rates, RNA editing, retrotranscription, gene loss and intron loss. Here we have sequenced the complete mitogenomes of Zostera marina and Stratiotes aloides, which together with previously sequenced mitogenomes from Butomus and Spirodela, provide new evolutionary evidence of genome size reduction, gene loss and transfer to the nucleus. The Zostera mitogenome includes a large portion of DNA transferred from the plastome, yet it is the smallest known mitogenome from a non-parasitic plant. Using a broad sample of the Alismatales, the evolutionary history of ribosomal protein gene loss is analyzed. In Zostera almost all ribosomal protein genes are lost from the mitogenome, but only some can be found in the nucleus.

The complete plastome sequence of Carissa macrocarpa (Eckl.) A. DC. (Apocynaceae)

In this study, we determined the complete plastome sequence of Carissa macrocarpa (Eckl.) A. DC. (Apocynaceae) (NCBI acc. no. KX364402). The gene order and structure of the C. macrocarpa plastome are similar to those of a typical angiosperm. The complete plastome is 155,297 bp in length, and consists of a large single-copy region of 85,586 bp and a small single-copy region of 18,131 bp, which are separated by two inverted repeats of 25,792 bp. The plastome contains 113 genes, of which 79 are protein-coding genes, 30 are tRNA genes and 4 are rRNA genes. Sixteen genes contained one intron and two genes have two introns. The average A–T content of the plastome is 62.0%. A total of 31 simple sequence repeat loci were identified within the genome. Phylogenetic analysis revealed that C. macrocarpa is a member of the paraphyletic subfamily Rauvolfioideae of Apocynaceae. The sister group relationship of C. macrocarpa to the Apocynoideae–Asclepiadoideae clade is supported by 100% bootstrap values.

The nuclear genome of Rhazya stricta and the evolution of alkaloid diversity in a medically relevant clade of Apocynaceae

Alkaloid accumulation in plants is activated in response to stress, is limited in distribution and specific alkaloid repertoires are variable across taxa. Rauvolfioideae (Apocynaceae, Gentianales) represents a major center of structural expansion in the monoterpenoid indole alkaloids (MIAs) yielding thousands of unique molecules including highly valuable chemotherapeutics. The paucity of genome-level data for Apocynaceae precludes a deeper understanding of MIA pathway evolution hindering the elucidation of remaining pathway enzymes and the improvement of MIA availability in planta or in vitro. We sequenced the nuclear genome of Rhazya stricta (Apocynaceae, Rauvolfioideae) and present this high quality assembly in comparison with that of coffee (Rubiaceae, Coffea canephora, Gentianales) and others to investigate the evolution of genome-scale features. The annotated Rhazya genome was used to develop the community resource, RhaCyc, a metabolic pathway database. Gene family trees were constructed to identify homologs of MIA pathway genes and to examine their evolutionary history. We found that, unlike Coffea, the Rhazya lineage has experienced many structural rearrangements. Gene tree analyses suggest recent, lineage-specific expansion and diversification among homologs encoding MIA pathway genes in Gentianales and provide candidate sequences with the potential to close gaps in characterized pathways and support prospecting for new MIA production avenues.

Renaissance in phytomedicines: promising implications of NGS technologies

Medicinal plant research is growing significantly in faith to discover new and more biologically compatible phytomedicines. Deposition of huge genome/trancriptome sequence data assisted by NGS technologies has revealed the new possibilities for producing upgraded bioactive molecules in medicinal plants. Growing interest of investors and consumers in the herbal drugs raises the need for extensive research to open the facts and details of every inch of life canvas of medicinal plants to produce improved quality of phytomedicines. As in agriculture crops, knowledge emergence from medicinal plant's genome/transcriptome, can be used to assure their amended quality and these improved varieties are then transported to the fields for cultivation. Genome studies generate huge sequence data which can be exploited further for obtaining information regarding genes/gene clusters involved in biosynthesis as well as regulation. This can be achieved rapidly at a very large scale with NGS platforms. Identification of new RNA molecules has become possible, which can lead to the discovery of novel compounds. Sequence information can be combined with advanced phytochemical and bioinformatics tools to discover functional herbal drugs. Qualitative and quantitative analysis of small RNA species put a light on the regulatory aspect of biosynthetic pathways for phytomedicines. Inter or intra genomic as well as transcriptomic interactive processes for biosynthetic pathways can be elucidated in depth. Quality management of herbal material will also become rapid and high throughput. Enrichment of sequence information will be used to engineer the plants to get more efficient phytopharmaceuticals. The present review comprises of role of NGS technologies to boost genomic studies of pharmaceutically important plants and further, applications of sequence information aiming to produce enriched phytomedicines. Emerging knowledge from the medicinal plants genome/transcriptome can give birth to deep understanding of the processes responsible for biosynthesis of medicinally important compounds.

Detecting and Characterizing the Highly Divergent Plastid Genome of the Nonphotosynthetic Parasitic Plant Hydnora visseri (Hydnoraceae)

Plastid genomes of photosynthetic flowering plants are usually highly conserved in both structure and gene content. However, the plastomes of parasitic and mycoheterotrophic plants may be released from selective constraint due to the reduction or loss of photosynthetic ability. Here we present the greatly reduced and highly divergent, yet functional, plastome of the nonphotosynthetic holoparasite Hydnora visseri (Hydnoraceae, Piperales). The plastome is 27 kb in length, with 24 genes encoding ribosomal proteins, ribosomal RNAs, tRNAs, and a few nonbioenergetic genes, but no genes related to photosynthesis. The inverted repeat and the small single copy region are only approximately 1.5 kb, and intergenic regions have been drastically reduced. Despite extreme reduction, gene order and orientation are highly similar to the plastome of Piper cenocladum, a related photosynthetic plant in Piperales. Gene sequences in Hydnora are highly divergent and several complementary approaches using the highest possible sensitivity were required for identification and annotation of this plastome. Active transcription is detected for all of the protein-coding genes in the plastid genome, and one of two introns is appropriately spliced out of rps12 transcripts. The whole-genome shotgun read depth is 1,400× coverage for the plastome, whereas the mitochondrial genome is covered at 40× and the nuclear genome at 2×. Despite the extreme reduction of the genome and high sequence divergence, the presence of syntenic, long transcriptionally active open-reading frames with distant similarity to other plastid genomes and a high plastome stoichiometry relative to the mitochondrial and nuclear genomes suggests that the plastome remains functional in H. visseri. A four-stage model of gene reduction, including the potential for complete plastome loss, is proposed to account for the range of plastid genomes in nonphotosynthetic plants.

Organellar Genomes of White Spruce (Picea glauca): Assembly and Annotation

The genome sequences of the plastid and mitochondrion of white spruce (Picea glauca) were assembled from whole-genome shotgun sequencing data using ABySS. The sequencing data contained reads from both the nuclear and organellar genomes, and reads of the organellar genomes were abundant in the data as each cell harbors hundreds of mitochondria and plastids. Hence, assembly of the 123-kb plastid and 5.9-Mb mitochondrial genomes were accomplished by analyzing data sets primarily representing low coverage of the nuclear genome. The assembled organellar genomes were annotated for their coding genes, ribosomal RNA, and transfer RNA. Transcript abundances of the mitochondrial genes were quantified in three developmental tissues and five mature tissues using data from RNA-seq experiments. C-to-U RNA editing was observed in the majority of mitochondrial genes, and in four genes, editing events were noted to modify ACG codons to create cryptic AUG start codons. The informatics methodology presented in this study should prove useful to assemble organellar genomes of other plant species using whole-genome shotgun sequencing data.

Genomics and Evolution in Traditional Medicinal Plants: Road to a Healthier Life

Medicinal plants have long been utilized in traditional medicine and ethnomedicine worldwide. This review presents a glimpse of the current status of and future trends in medicinal plant genomics, evolution, and phylogeny. These dynamic fields are at the intersection of phytochemistry and plant biology and are concerned with the evolution mechanisms and systematics of medicinal plant genomes, origin and evolution of the plant genotype and metabolic phenotype, interaction between medicinal plant genomes and their environment, the correlation between genomic diversity and metabolite diversity, and so on. Use of the emerging high-end genomic technologies can be expanded from crop plants to traditional medicinal plants, in order to expedite medicinal plant breeding and transform them into living factories of medicinal compounds. The utility of molecular phylogeny and phylogenomics in predicting chemodiversity and bioprospecting is also highlighted within the context of natural-product-based drug discovery and development. Representative case studies of medicinal plant genome, phylogeny, and evolution are summarized to exemplify the expansion of knowledge pedigree and the paradigm shift to the omics-based approaches, which update our awareness about plant genome evolution and enable the molecular breeding of medicinal plants and the sustainable utilization of plant pharmaceutical resources.