The complete chloroplast and mitochondrial genome sequences of Boea hygrometrica: insights into the evolution of plant organellar genomes

The mitochondrial genome of Lavandula angustifolia Mill. (Lamiaceae) sheds light on its genome structure and gene transfer between organelles

BACKGROUND

Lavandula angustifolia holds importance as an aromatic plant with extensive applications spanning the fragrance, perfume, cosmetics, aromatherapy, and spa sectors. Beyond its aesthetic and sensory applications, this plant offers medicinal benefits as a natural herbal remedy and finds use in household cleaning products. While extensive genomic data, inclusive of plastid and nuclear genomes, are available for this species, researchers have yet to characterize its mitochondrial genome. This gap in knowledge hampers deeper understanding of the genome organization and its evolutionary significance.

RESULTS

Through the course of this study, we successfully assembled and annotated the mitochondrial genome of L. angustifolia, marking a first in this domain. This assembled genome encompasses 61 genes, which comprise 34 protein-coding genes, 24 transfer RNA genes, and three ribosomal RNA genes. We identified a chloroplast sequence insertion into the mitogenome, which spans a length of 10,645 bp, accounting for 2.94% of the mitogenome size. Within these inserted sequences, there are seven intact tRNA genes (trnH-GUG, trnW-CCA, trnD-GUC, trnS-GGA, trnN-GUU, trnT-GGU, trnP-UGG) and four complete protein-coding genes (psbA, rps15, petL, petG) of chloroplast derivation. Additional discoveries include 88 microsatellites, 15 tandem repeats, 74 palindromic repeats, and 87 forward long repeats. An RNA editing analysis highlighted an elevated count of editing sites in the cytochrome c oxidase genes, notably ccmB with 34 editing sites, ccmFN with 32, and ccmC with 29. All protein-coding genes showed evidence of cytidine-to-uracil conversion. A phylogenetic analysis, utilizing common protein-coding genes from 23 Lamiales species, yielded a tree with consistent topology, supported by high confidence values.

CONCLUSIONS

Analysis of the current mitogenome resource revealed its typical circular genome structure. Notably, sequences originally from the chloroplast genome were found within the mitogenome, pointing to the occurrence of horizontal gene transfer between organelles. This assembled mitogenome stands as a valuable resource for subsequent studies on mitogenome structures, their evolution, and molecular biology.

Assembly and comparative analysis of the first complete mitochondrial genome of Salix psammophila, a good windbreak and sand fixation shrub

Salix psammophila, commonly known as the sandlive willow, is a vital shrub species within the Salicaceae family, particularly significant for its ecological role in regions susceptible to desertification and sandy soils. In this study, we assembled the complete S. psammophila mitochondrial genome using Pacbio HiFi third-generation sequencing data. The genome was found to be a typical single circular structure, with a total length of 715,555 bp and a GC content of 44.89%. We annotated 33 unique protein-coding genes (PCGs), which included 24 core mitochondrial genes and 9 variable genes, as well as 18 tRNA genes (5 of which were multicopy genes) and 3 rRNA genes. Comparative analysis of the PCGs from the mitochondrial genomes of S. psammophila, Populus deltoides, Populus simonii, Salix wilsonii, and Salix suchowensis revealed that these genes are relatively conserved within the Salicaceae family, with variability primarily occurring in the ribosomal protein genes. The absence of the rps14, which encodes a ribosomal protein, may have played a role in the evolution of stress tolerance in Salicaceae plants. Additionally, we identified 232 SSRs, 19 tandem repeat sequences, and 236 dispersed repeat sequences in the S. psammophila mitochondrial genome, with palindromic and forward repeats being the most abundant. The longest palindromic repeat measured 260 bp, while the longest forward repeat was 86,068 bp. Furthermore, 324 potential RNA editing sites were discovered, all involving C-to-U edits, with the nad4 having the highest number of edits. These findings provide valuable insights into the phylogenetic and genetic research of Salicaceae plants.

Characterization of the complete chloroplast genome of Rhodiola sachalinensis and comparative analysis with its congeneric plants

Rhodiola, belonging to the Crassulaceae family, is a perennial herbaceous plant genus. There are about 90 Rhodiola species worldwide, some of which have been reported to have medicinal properties. Rhodiola sachalinensis is a perennial medicinal herb within this genus and, in the present study, its chloroplast genome was sequenced, assembled, annotated and compared with 24 other Rhodiola species. The results obtained show that the chloroplast genome of R. sachalinensis is 151 595 bp long and has a CG content of 37.7%. The inverted repeats (IR) region of the Rhodiola chloroplast genome is the most conserved region, with the main differences being observed in the ycf1 and ndhF genes at the IRb-small single copy boundary, and rps19 and trnH genes at the IRa-large single copy boundary. Phylogenetic analysis showed that Rhodiola species form two major clades, and species with recorded medicinal properties, clustered together in one branch except for R. dumulosa. Within the genus, R. sachalinensis is most closely related to Rhodiola rosea, although comparative analyses showed that only R. sachalinensis and Rhodiola subopposita contained the psbZ gene, which encodes a highly conserved protein subunit of the Photosystem II core complex. Overall, the present study contributes to the understanding of the chloroplast genome of Rhodiola species, and provides a theoretical basis for the study of their genetic diversity and possible use as medicinal plants.

Comparative Analysis and Phylogeny of the Complete Chloroplast Genomes of Nine Cynanchum (Apocynaceae) Species

Cynanchum belongs to the Apocynaceae family and is a morphologically diverse genus that includes around 200 shrub or perennial herb species. Despite the utilization of CPGs, few molecular phylogenetic studies have endeavored to elucidate infrafamilial relationships within Cynanchum through extensive taxon sampling. In this research, we constructed a phylogeny and estimated divergence time based on the chloroplast genomes (CPGs) of nine Cynanchum species. We sequenced and annotated nine chloroplast (CP) genomes in this study. The comparative analysis of these genomes from these Cynanchum species revealed a typical quadripartite structure, with a total sequence length ranging from 158,283 to 161,241 base pairs (bp). The CP genome (CPG) was highly conserved and moderately differentiated. Through annotation, we identified a total of 129-132 genes. Analysis of the boundaries of inverted repeat (IR) regions showed consistent positioning: the rps19 gene was located in the IRb region, varying from 46 to 50 bp. IRb/SSC junctions were located between the trnN and ndhF genes. We did not detect major expansions or contractions in the IR region or rearrangements or insertions in the CPGs of the nine Cynanchum species. The results of SSR analysis revealed a variation in the number of SSRs, ranging from 112 to 150. In five types of SSRs, the largest number was mononucleotide repeats, and the smallest number was hexanucleotide repeats. The number of long repeats in the cp genomes of nine Cynanchum species was from 35 to 80. In nine species of Cynanchum, the GC3s values ranged from 26.80% to 27.00%, indicating a strong bias towards A/U-ending codons. Comparative analyses revealed four hotspot regions in the CPG, ndhA-ndhH, trnI-GAU-rrn16, psbI-trnS-GCU, and rps7-ndhB, which could potentially serve as molecular markers. In addition, phylogenetic tree construction based on the CPG indicated that the nine Cynanchum species formed a monophyletic group. Molecular dating suggested that Cynanchum diverged from its sister genus approximately 18.87 million years ago (Mya) and species diversification within the Cynanchum species primarily occurred during the recent Miocene epoch. The divergence time estimation presented in this study will facilitate future research on Cynanchum, aid in species differentiation, and facilitate diverse investigations into this economically and ecologically important genus.

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

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

Twelve newly assembled jasmine chloroplast genomes: unveiling genomic diversity, phylogenetic relationships and evolutionary patterns among Oleaceae and Jasminum species

BACKGROUND

Jasmine (Jasminum), renowned for its ornamental value and captivating fragrance, has given rise to numerous species and accessions. However, limited knowledge exists regarding the evolutionary relationships among various Jasminum species.

RESULTS

In the present study, we sequenced seven distinct Jasminum species, resulting in the assembly of twelve high-quality complete chloroplast (cp) genomes. Our findings revealed that the size of the 12 cp genomes ranged from 159 to 165 kb and encoded 134-135 genes, including 86-88 protein-coding genes, 38-40 tRNA genes, and 8 rRNA genes. J. nudiflorum exhibited a larger genome size compared to other species, mainly attributed to the elevated number of forward repeats (FRs). Despite the typically conservative nature of chloroplasts, variations in the presence or absence of accD have been observed within J. sambac. The calculation of nucleotide diversity (Pi) values for 19 cp genomes indicated that potential mutation hotspots were more likely to be located in LSC regions than in other regions, particularly in genes ycf2, rbcL, atpE, ndhK, and ndhC (Pi > 0.2). Ka/Ks values revealed strong selection pressure on the genes rps2, atpA, rpoA, rpoC1, and rpl33 when comparing J. sambac with the three most closely related species (J. auriculatum, J. multiflorum, and J. dichotomum). Additionally, SNP identification, along with the results of Structure, PCA, and phylogenetic tree analyses, divided the Jasminum cp genomes into six groups. Notably, J. polyanthum showed gene flow signals from both the G5 group (J. nudiflorum) and the G3 group (J. tortuosum and J. fluminense). Phylogenetic tree analysis reflected that most species from the same genus clustered together with robust support in Oleaceae, strongly supporting the monophyletic nature of cp genomes within the genus Jasminum.

CONCLUSION

Overall, this study provides comprehensive insights into the genomic composition, variation, and phylogenetic relationships among various Jasminum species. These findings enhance our understanding of the genetic diversity and evolutionary history of Jasminum.

Comparative Analysis of Codon Usage Patterns in Nuclear and Chloroplast Genome of Dalbergia (Fabaceae)

The Dalbergia plants are widely distributed across more than 130 tropical and subtropical countries and have significant economic and medicinal value. Codon usage bias (CUB) is a critical feature for studying gene function and evolution, which can provide a better understanding of biological gene regulation. In this study, we comprehensively analyzed the CUB patterns of the nuclear genome, chloroplast genome, and gene expression, as well as systematic evolution of Dalbergia species. Our results showed that the synonymous and optimal codons in the coding regions of both nuclear and chloroplast genome of Dalbergia preferred ending with A/U at the third codon base. Natural selection was the primary factor affecting the CUB features. Furthermore, in highly expressed genes of Dalbergia odorifera, we found that genes with stronger CUB exhibited higher expression levels, and these highly expressed genes tended to favor the use of G/C-ending codons. In addition, the branching patterns of the protein-coding sequences and the chloroplast genome sequences were very similar in the systematic tree, and different with the cluster from the CUB of the chloroplast genome. This study highlights the CUB patterns and features of Dalbergia species in different genomes, explores the correlation between CUB preferences and gene expression, and further investigates the systematic evolution of Dalbergia, providing new insights into codon biology and the evolution of Dalbergia plants.

Decoding the complete chloroplast genome of Cissus quadrangularis: insights into molecular structure, comparative genome analysis and mining of mutational hotspot regions

Cissus quadrangularis L., a member of the Vitaceae family, is an important medicinal plant with widespread application in Indian traditional medicines. C. quadrangularis L. whole chloroplast genome of 160,404 bp was assembled using a genome skimming approach from the whole genome library. The assembled chloroplast genome contained a large single-copy region (88,987 bp), a small single-copy region (18,621 bp), and pairs of inverted repeat regions (26,398 bp). It also comprised 133 genes, including 37 tRNAs, eight rRNAs, and 88 protein-coding genes. Aside from that, we annotated three genes atpH, petB, and psbL, as well as one duplicated copy of the ycf1 gene in C. quadrangularis L. that had previously been missing from the annotation of compared Cissus chloroplast genomes. Five divergent hotspot regions such as petA_psbJ (0.1237), rps16_trnQ-UUG (0.0913), psbC_trnS-UGA (0.0847), rps15_ycf1 (0.0788), and rps2_rpoC2 (0.0788) were identified in the investigation that could aid in future species discrimination. Surprisingly, we found the overlapping genes ycf1 and ndhF on the IRb/SSC junction, rarely seen in angiosperms. The results of the phylogenetic study showed that the genomes of the Cissus species under study formed a single distinct clade. The detailed annotations given in this study could be useful in the future for genome annotations of Cissus species. The current findings of the study have the potential to serve as a useful resource for future research in the field of population genetics and the evolutionary relationships in the Cissus genus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01312-w.

Comparative genomics analysis of endangered wild Egyptian Moringa peregrina (Forssk.) Fiori plastome, with implications for the evolution of Brassicales order

Moringa is a mono-genus belonging to the Moringaceae family, which includes 13 species. Among them, Moringa peregrina is plant species native to the Arabian Peninsula, Southern Sinai in Egypt, and the Horn of Africa, and comprehensive studies on its nutritional, industrial, and medicinal values have been performed. Herein, we sequenced and analyzed the initial complete chloroplast genome of Moringa peregrina. Concurrently, we analyzed the new chloroplast genome along with 25 chloroplast genomes related to species representing eight families in the Brassicales order. The results indicate that the plastome sequence of M. peregrina consists of 131 genes, with an average GC content of 39.23%. There is a disparity in the IR regions of the 26 species ranging from 25,804 to 31,477 bp. Plastome structural variations generated 20 hotspot regions that could be considered prospective DNA barcode locations in the Brassicales order. Tandem repeats and SSR structures are reported as significant evidence of structural variations among the 26 tested specimens. Furthermore, selective pressure analysis was performed to estimate the substitution rate within the Moringaceae family, which revealing that the ndhA and accD genes are under positive selective pressure. The phylogenetic analysis of the Brassicales order produced an accurate monophyletic annotation cluster of the Moringaceae and Capparaceae species, offering unambiguous identification without overlapping groups between M. oleifera and M. peregrina, which are genetically strongly associated. Divergence time estimation suggests that the two Moringa species recently diversified, 0.467 Ma. Our findings highlight the first complete plastome of the Egyptian wild-type of M. peregrina, which can be used for determining plastome phylogenetic relationships and systematic evolution history within studies on the Moringaceae family.

Chloroplast genomic comparison provides insights into the evolution of seagrasses

BACKGROUND

Seagrasses are a polyphyletic group of monocotyledonous angiosperms that have evolved to live entirely submerged in marine waters. Thus, these species are ideal for studying plant adaptation to marine environments. Herein, we sequenced the chloroplast (cp) genomes of two seagrass species (Zostera muelleri and Halophila ovalis) and performed a comparative analysis of them with 10 previously published seagrasses, resulting in various novel findings.

RESULTS

The cp genomes of the seagrasses ranged in size from 143,877 bp (Zostera marina) to 178,261 bp (Thalassia hemprichii), and also varied in size among different families in the following order: Hydrocharitaceae > Cymodoceaceae > Ruppiaceae > Zosteraceae. The length differences between families were mainly related to the expansion and contraction of the IR region. In addition, we screened out 2,751 simple sequence repeats and 1,757 long repeat sequence types in the cp genome sequences of the 12 seagrass species, ultimately finding seven hot spots in coding regions. Interestingly, we found nine genes with positive selection sites, including two ATP subunit genes (atpA and atpF), three ribosome subunit genes (rps4, rps7, and rpl20), one photosystem subunit gene (psbH), and the ycf2, accD, and rbcL genes. These gene regions may have played critical roles in the adaptation of seagrasses to diverse environments. In addition, phylogenetic analysis strongly supported the division of the 12 seagrass species into four previously recognized major clades. Finally, the divergence time of the seagrasses inferred from the cp genome sequences was generally consistent with previous studies.

CONCLUSIONS

In this study, we compared chloroplast genomes from 12 seagrass species, covering the main phylogenetic clades. Our findings will provide valuable genetic data for research into the taxonomy, phylogeny, and species evolution of seagrasses.

Complete chloroplast genomes and comparative analysis of Ligustrum species

In this study, we assembled and annotated the chloroplast (cp) genomes of four Ligustrum species, L. sinense, L. obtusifolium, L. vicaryi, and L. ovalifolium 'Aureum'. Including six other published Ligustrum species, we compared various characteristics such as gene structure, sequence alignment, codon preference, and nucleic acid diversity, and performed positive-selection genes screening and phylogenetic analysis. The results showed that the cp genome of Ligustrum was 162,185-166,800 bp in length, with a circular tetrad structure, including a large single-copy region (86,885-90,106 bp), a small single-copy region (11,446-11,499 bp), and a pair of IRa and IRb sequences with the same coding but in opposite directions (31,608-32,624 bp). This structure is similar to the cp genomes of most angiosperms. We found 132-137 genes in the cp genome of Ligustrum, including 89-90 protein-coding genes, 35-39 tRNAs, and 8 rRNAs. The GC content was 37.93-38.06% and varied among regions, with the IR region having the highest content. The single-nucleotide (A/T)n was dominant in simple-sequence repeats of the Ligustrum cp genome, with an obvious A/T preference. Six hotspot regions were identified from multiple sequence alignment of Ligustrum; the ycf1 gene region and the clpP1 exon region can be used as potential DNA barcodes for the identification and phylogeny of the genus Ligustrum. Branch-site model and Bayes empirical Bayes (BEB) analysis showed that four protein-coding genes (accD, clpP, ycf1, and ycf2) were positively selected, and BEB analysis showed that accD and rpl20 had positively selected sites. A phylogenetic tree of Oleaceae species was constructed based on the whole cp genomes, and the results were consistent with the traditional taxonomic results. The phylogenetic results showed that genus Ligustrum is most closely related to genus Syringa. Our study provides important genetic information to support further investigations of the phylogenetic development and adaptive evolution of Ligustrum species.

Analysis of the chloroplast genome and phylogenetic evolution of three species of Syringa

BACKGROUND

By the time our study was completed, the chloroplast genomes of Syringa oblata, S. pubescents subsp. Microphylla, and S. reticulate subsp. Amurensis had not been sequenced, and their genetic background was not clear.

THE RESEARCH CONTENT

In this study, the chloroplast genomes of Syringa oblata, S. pubescents subsp. Microphylla, S. reticulate subsp. Amurensis, and five other species of Syringa were sequenced for a comparative genomics analysis, inverted repeat (IR) boundary analysis, collinearity analysis, codon preference analysis and a nucleotide variability analysis. Differences in the complete chloroplast genomes of 30 species of Oleaceae were compared with that of S. oblata as the reference species, and Ginkgo biloba was used as the out group to construct the phylogenetic tree.

RESULTS

The results showed that the chloroplast genomes of S. oblata, S. pubescents subsp. Microphylla, and S. reticulate subsp. Amurensis were similar to those of other angiosperms and showed a typical four-segment structure, with full lengths of 155,569, 160,491, 155,419, and protein codes of 88, 95, and 87, respectively. Because the IR boundary of S. pubescents subsp. Microphylla was significantly expanded to the large single copy (LSC) region, resulting in complete replication of some genes in the IR region, the LSC region of S. pubescents subsp. Microphylla was significantly shorter than those of S. oblate and S. reticulate subsp. Amurensis. Similar to most higher plants, these three species have a preference for their codons ending with A/T.

CONCLUSIONS

We consider the genus Syringa to be a synphyletic group. The nucleotide variability and phylogenetic analyses showed that Syringa differentiated before Ligustrum and Ligustrum developed from Syringa. We propose removing the existing section division and directly dividing Syringa into five series.

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.

Plastome phylogenomics provide new perspective into the phylogeny and evolution of Betulaceae (Fagales)

BACKGROUND

Betulaceae is a relatively small but morphologically diverse family, with many species having important economic and ecological values. Although plastome structure of Betulaceae has been reported sporadically, a comprehensive exploration for plastome evolution is still lacking. Besides, previous phylogenies had been constructed based on limited gene fragments, generating unrobust phylogenetic framework and hindering further studies on divergence ages, biogeography and character evolution. Here, 109 plastomes (sixteen newly assembled and 93 previously published) were subject to comparative genomic and phylogenomic analyses to reconstruct a robust phylogeny and trace the diversification history of Betulaceae.

RESULTS

All Betulaceae plastomes were highly conserved in genome size, gene order, and structure, although specific variations such as gene loss and IR boundary shifts were revealed. Ten divergent hotspots, including five coding regions (P_i > 0.02) and five noncoding regions (P_i > 0.035), were identified as candidate DNA barcodes for phylogenetic analysis and species delimitation. Phylogenomic analyses yielded high-resolution topology that supported reciprocal monophyly between Betula and Alnus within Betuloideae, and successive divergence of Corylus, Ostryopsis, and Carpinus-Ostrya within Coryloideae. Incomplete lineage sorting and hybridization may be responsible for the mutual paraphyly between Ostrya and Carpinus. Betulaceae ancestors originated from East Asia during the upper Cretaceous; dispersals and subsequent vicariance accompanied by historical environment changes contributed to its diversification and intercontinental disjunction. Ancestral state reconstruction indicated the acquisition of many taxonomic characters was actually the results of parallel or reversal evolution.

CONCLUSIONS

Our research represents the most comprehensive taxon-sampled and plastome-level phylogenetic inference for Betulaceae to date. The results clearly document global patterns of plastome structural evolution, and established a well-supported phylogeny of Betulaceae. The robust phylogenetic framework not only provides new insights into the intergeneric relationships, but also contributes to a perspective on the diversification history and evolution of the family.

The Complete Chloroplast Genome Sequence of Laportea bulbifera (Sieb. et Zucc.) Wedd. and Comparative Analysis with Its Congeneric Species

Laportea bulbifera (L. bulbifera) is an important medicinal plant of Chinese ethnic minorities, with high economic and medicinal value. However, the medicinal materials of the genus Laportea are prone to be misidentified due to the similar morphological characteristics of the original plants. Thus, it is crucial to discover their molecular marker points and to precisely identify these species for their exploitation and conservation. Here, this study reports detailed information on the complete chloroplast (cp) of L. bulbifera. The result indicates that the cp genome of L. bulbifera of 150,005 bp contains 126 genes, among them, 37 tRNA genes and 81 protein-coding genes. The analysis of repetition demonstrated that palindromic repeats are more frequent. In the meantime, 39 SSRs were also identified, the majority of which were mononucleotides Adenine-Thymine (A-T). Furthermore, we compared L. bulbifera with eight published Laportea plastomes, to explore highly polymorphic molecular markers. The analysis identified four hypervariable regions, including rps16, ycf1, trnC-GCA and trnG-GCC. According to the phylogenetic analysis, L. bulbifera was most closely related to Laportea canadensis (L. canadensis), and the molecular clock analysis speculated that the species originated from 1.8216 Mya. Overall, this study provides a more comprehensive analysis of the evolution of L. bulbifera from the perspective of phylogenetic and intrageneric molecular variation in the genus Laportea, which is useful for providing a scientific basis for further identification, taxonomic, and evolutionary studies of the genus.

Plastome structure, phylogenomic analyses and molecular dating of Arecaceae

Arecaceae is a species-rich clade of Arecales, while also being regarded as a morphologically diverse angiosperm family with numerous species having significant economic, medicinal, and ornamental value. Although in-depth studies focused on the chloroplast structure of Arecaceae, as well as inferring phylogenetic relationships using gene fragments, have been reported in recent years, a comprehensive analysis of the chloroplast structure of Arecaceae is still needed. Here we perform a comprehensive analysis of the structural features of the chloroplast genome of Arecaceae, compare the variability of gene sequences, infer phylogenetic relationships, estimate species divergence times, and reconstruct ancestral morphological traits. In this study, 74 chloroplast genomes of Arecaceae were obtained, covering five subfamilies. The results show that all chloroplast genomes possess a typical tetrad structure ranging in size between 153,806-160,122 bp, with a total of 130-137 genes, including 76-82 protein-coding genes, 29-32 tRNA genes, and 4 rRNA genes. Additionally, the total GC content was between 36.9-37.7%. Analysis of the SC/IR boundary indicated that the IR region underwent expansion or contraction. Phylogenetic relationships indicate that all five subfamilies in Arecaceae are monophyletic and that Ceroxyloideae and Arecoideae are sister groups (BS/PP = 100/1). The results of molecular dating indicate that the age of the crown group of Arecaceae is likely to be 96.60 [84.90-107.60] Ma, while the age of the stem group is 102.40 [93.44-111.17] Ma. Reconstruction of ancestral traits indicate that the ancestral characteristics of the family include monoecious plants, one seed, six stamens, and a smooth pericarp.

Comparative chloroplast genome and phylogenetic analyses of Chinese Polyspora

Polyspora Sweet (Theaceae) are winter ornamental landscape plants native to southern and southeastern Asia, some of which have medicinal value. The chloroplast (cp) genome data of Polyspora are scarce, and the gene evolution and interspecific relationship are still unclear. In this study, we sequenced and annotated Polyspora chrysandra cp genome and combined it with previously published genomes for other Chinese Polyspora species. The results showed that cp genomes of six Chinese Polyspora varied in length between 156,452 bp (P. chrysandra) and 157,066 bp (P. speciosa), but all contained 132 genes, with GC content of 37.3%, and highly similar genes distribution and codon usage. A total of eleven intergenic spacer regions were found having the highest levels of divergence, and eight divergence hotspots were identified as molecular markers for Phylogeography and genetic diversity studies in Polyspora. Gene selection pressure suggested that five genes were subjected to positive selection. Phylogenetic relationships among Polyspora species based on the complete cp genomes were supported strongly, indicating that the cp genomes have the potential to be used as super barcodes for further analysis of the phylogeny of the entire genus. The cp genomes of Chinese Polyspora species will provide valuable information for species identification, molecular breeding and evolutionary analysis of genus Polyspora.

Complete chloroplast genome of Lilium ledebourii (Baker) Boiss and its comparative analysis: lights into selective pressure and adaptive evolution

Lilium ledebourii (Baker) Boiss is a rare species, which exhibits valuable traits. However, before its genetic diversity and evolutionary were uncovered, its wild resources were jeopardized. Moreover, some ambiguities in phylogenetic relationships of this genus remain unresolved. Therefore, obtaining the whole chloroplast sequences of L. ledebourii and its comparative analysis along with other Lilium species is crucial and pivotal to understanding the evolution of this genus as well as the genetic populations. A multi-scale genome-level analysis, especially selection pressure, was conducted. Detailed third‑generation sequencing and analysis revealed a whole chloroplast genome of 151,884 bp, with an ordinary quadripartite and protected structure comprising 37.0% GC. Overall, 113 different genes were recognized in the chloroplast genome, consisting of 30 distinct tRNA genes, four distinct ribosomal RNAs genes, and 79 unique protein-encoding genes. Here, 3234 SSRs and 2053 complex repeats were identified, and a comprehensive analysis was performed for IR expansion and contraction, and codon usage bias. Moreover, genome-wide sliding window analysis revealed the variability of rpl32-trnL-ccsA, petD-rpoA, ycf1, psbI-trnS-trnG, rps15-ycf1, trnR, trnT-trnL, and trnP-psaJ-rpl33 were higher among the 48 Lilium cp genomes, displaying higher variability of nucleotide in SC regions. Following 1128 pairwise comparisons, ndhB, psbJ, psbZ, and ycf2 exhibit zero synonymous substitution, revealing divergence or genetic restriction. Furthermore, out of 78 protein-coding genes, we found that accD and rpl36 under positive selection: however, at the entire-chloroplast protein scale, the Lilium species have gone through a purifying selection. Also, a new phylogenetic tree for Lilium was rebuilt, and we believe that the Lilium classification is clearer than before. The genetic resources provided here will aid future studies in species identification, population genetics, and Lilium conservation.

The first genome for the Cape Primrose Streptocarpus rexii (Gesneriaceae), a model plant for studying meristem-driven shoot diversity

Cape Primroses (Streptocarpus, Gesneriaceae) are an ideal study system for investigating the genetics underlying species diversity in angiosperms. Streptocarpus rexii has served as a model species for plant developmental research for over five decades due to its unusual extended meristem activity present in the leaves. In this study, we sequenced and assembled the complete nuclear, chloroplast, and mitochondrial genomes of S. rexii using Oxford Nanopore Technologies long read sequencing. Two flow cells of PromethION sequencing resulted in 32 billion reads and were sufficient to generate a draft assembly including the chloroplast, mitochondrial and nuclear genomes, spanning 776 Mbp. The final nuclear genome assembly contained 5,855 contigs, spanning 766 Mbp of the 929-Mbp haploid genome with an N50 of 3.7 Mbp and an L50 of 57 contigs. Over 70% of the draft genome was identified as repeats. A genome repeat library of Gesneriaceae was generated and used for genome annotation, with a total of 45,045 genes annotated in the S. rexii genome. Ks plots of the paranomes suggested a recent whole genome duplication event, shared between S. rexii and Primulina huaijiensis. A new chloroplast and mitochondrial genome assembly method, based on contig coverage and identification, was developed, and successfully used to assemble both organellar genomes of S. rexii. This method was developed into a pipeline and proved widely applicable. The nuclear genome of S. rexii and other datasets generated and reported here will be invaluable resources for further research to aid in the identification of genes involved in morphological variation underpinning plant diversification.

Complete Chloroplast Genome Sequence of Fortunella venosa (Champ. ex Benth.) C.C.Huang (Rutaceae): Comparative Analysis, Phylogenetic Relationships, and Robust Support for Its Status as an Independent Species

Fortunella venosa (Rutaceae) is an endangered species endemic to China and its taxonomic status has been controversial. The genus Fortunella contains a variety of important economic plants with high value in food, medicine, and ornamental. However, the placement of Genus Fortunella into Genus Citrus has led to controversy on its taxonomy and Systematics. In this present research, the Chloroplast genome of F. venosa was sequenced using the second-generation sequencing, and its structure and phylogenetic relationship analyzed. The results showed that the Chloroplast genome size of F. venosa was 160,265 bp, with a typical angiosperm four-part ring structure containing a large single copy region (LSC) (87,597 bp), a small single copy region (SSC) (18,732 bp), and a pair of inverted repeat regions (IRa\IRb) (26,968 bp each). There are 134 predicted genes in Chloroplast genome, including 89 protein-coding genes, 8 rRNAs, and 37 tRNAs. The GC-content of the whole Chloroplast genome was 43%, with the IR regions having a higher GC content than the LSC and the SSC regions. There were no rearrangements present in the Chloroplast genome; however, the IR regions showed obvious contraction and expansion. A total of 108 simple sequence repeats (SSRs) were present in the entire chloroplast genome and the nucleotide polymorphism was high in LSC and SSC. In addition, there is a preference for codon usage with the non-coding regions being more conserved than the coding regions. Phylogenetic analysis showed that species of Fortunella are nested in the genus of Citrus and the independent species status of F. venosa is supported robustly, which is significantly different from F. japonica. These findings will help in the development of DNA barcodes that can be useful in the study of the systematics and evolution of the genus Fortunella and the family Rutaceae.

Characterization of the complete chloroplast genome of Sanicula rubriflora F. Schmidt ex Maxim

Sanicula rubriflora refers to a kind of edible wild herbs, which possesses reputed medicinal value. Chloroplasts (cp) is essential organelle in plant cells and has a genome that can be inherited autonomously. The complete chloroplast genome of S. rubriflora was assemble and annotated in the present study. It was a circular molecular genome with a size of 155,700 bp in length, which has a typical quadripartite structure. For instance, a large single-copy (LSC) of 85,979 bp and a pair of inverted repeats (IRs) of 26,333 bp were disconnected by a small single-copy (SSC) of 17,053bp. A total of 130 genes have been annotated, including 86 protein-coding genes (PCGs), 36 transfer RNA genes (tRNAs) and 8 ribosomal RNA genes (rRNAs). The total GC content of the complete chloroplast genome reached 37.9%. A maximum likelihood phylogenetic analysis with the reported chloroplast genomes revealed that S. rubriflora is most closely related to Saniculachinensis in the phylogenetic relationships.

Structural variation and evolution of chloroplast tRNAs in green algae

As one of the important groups of the core Chlorophyta (Green algae), Chlorophyceae plays an important role in the evolution of plants. As a carrier of amino acids, tRNA plays an indispensable role in life activities. However, the structural variation of chloroplast tRNA and its evolutionary characteristics in Chlorophyta species have not been well studied. In this study, we analyzed the chloroplast genome tRNAs of 14 species in five categories in the green algae. We found that the number of chloroplasts tRNAs of Chlorophyceae is maintained between 28-32, and the length of the gene sequence ranges from 71 nt to 91 nt. There are 23-27 anticodon types of tRNAs, and some tRNAs have missing anticodons that are compensated for by other types of anticodons of that tRNA. In addition, three tRNAs were found to contain introns in the anti-codon loop of the tRNA, but the analysis scored poorly and it is presumed that these introns are not functional. After multiple sequence alignment, the Ψ-loop is the most conserved structural unit in the tRNA secondary structure, containing mostly U-U-C-x-A-x-U conserved sequences. The number of transitions in tRNA is higher than the number of transversions. In the replication loss analysis, it was found that green algal chloroplast tRNAs may have undergone substantial gene loss during the course of evolution. Based on the constructed phylogenetic tree, mutations were found to accompany the evolution of the Green algae chloroplast tRNA. Moreover, chloroplast tRNAs of Chlorophyceae are consistent with those of monocotyledons and gymnosperms in terms of evolutionary patterns, sharing a common multi-phylogenetic pattern and rooted in a rich common ancestor. Sequence alignment and systematic analysis of tRNA in chloroplast genome of Chlorophyceae, clarified the characteristics and rules of tRNA changes, which will promote the evolutionary relationship of tRNA and the origin and evolution of chloroplast.

Comparative analysis of chloroplast genome structure and molecular dating in Myrtales

BACKGROUND

Myrtales is a species rich branch of Rosidae, with many species having important economic, medicinal, and ornamental value. At present, although there are reports on the chloroplast structure of Myrtales, a comprehensive analysis of the chloroplast structure of Myrtales is lacking. Phylogenetic and divergence time estimates of Myrtales are mostly constructed by using chloroplast gene fragments, and the support for relationships is low. A more reliable method to reconstruct the species divergence time and phylogenetic relationships is by using whole chloroplast genomes. In this study, we comprehensively analyzed the structural characteristics of Myrtales chloroplasts, compared variation hotspots, and reconstructed the species differentiation time of Myrtales with four fossils and one secondary calibration point.

RESULTS

A total of 92 chloroplast sequences of Myrtales, representing six families, 16 subfamilies and 78 genera, were obtained including nine newly sequenced chloroplasts by whole genome sequencing. Structural analyses showed that the chloroplasts range in size between 152,214-171,315 bp and exhibit a typical four part structure. The IR region is between 23,901-36,747 bp, with the large single copy region spanning 83,691-91,249 bp and the small single copy region spanning 11,150-19,703 bp. In total, 123-133 genes are present in the chloroplasts including 77-81 protein coding genes, four rRNA genes and 30-31 tRNA genes. The GC content was 36.9-38.9%, with the average GC content being 37%. The GC content in the LSC, SSC and IR regions was 34.7-37.3%, 30.6-36.8% and 39.7-43.5%, respectively. By analyzing nucleotide polymorphism of the chloroplast, we propose 21 hypervariable regions as potential DNA barcode regions for Myrtales. Phylogenetic analyses showed that Myrtales and its corresponding families are monophyletic, with Combretaceae and the clade of Onagraceae + Lythraceae (BS = 100%, PP = 1) being sister groups. The results of molecular dating showed that the crown of Myrtales was most likely to be 104.90 Ma (95% HPD = 87.88-114.18 Ma), and differentiated from the Geraniales around 111.59 Ma (95% HPD = 95.50-118.62 Ma).

CONCLUSIONS

The chloroplast genome structure of Myrtales is similar to other angiosperms and has a typical four part structure. Due to the expansion and contraction of the IR region, the chloroplast genome sizes in this group are slightly different. The variation of noncoding regions of the chloroplast genome is larger than those of coding regions. Phylogenetic analysis showed that Combretaceae and Onagraceae + Lythraceae were well supported as sister groups. Molecular dating indicates that the Myrtales crown most likely originated during the Albian age of the Lower Cretaceous. These chloroplast genomes contribute to the study of genetic diversity and species evolution of Myrtales, while providing useful information for taxonomic and phylogenetic studies of Myrtales.

Decoding first complete chloroplast genome of toothbrush tree (Salvadora persica L.): insight into genome evolution, sequence divergence and phylogenetic relationship within Brassicales

BACKGROUND

Salvadora persica L. (Toothbrush tree - Miswak; family-Salvadoraceae) grows in the arid-land ecosystem and possesses economic and medicinal importance. The species, genus and the family have no genomic datasets available specifically on chloroplast (cp) genomics and taxonomic evolution. Herein, we have sequenced the complete chloroplast genome of S. persica for the first time and compared it with 11 related specie's cp genomes from the order Brassicales.

RESULTS

The S. persica cp genome was 153,379 bp in length containing a sizeable single-copy region (LSC) of 83,818 bp which separated from the small single-copy region (SSC) of 17,683 bp by two inverted repeats (IRs) each 25,939 bp. Among these genomes, the largest cp genome size (160,600 bp) was found in M. oleifera, while in S. persica it was the smallest (153,379 bp). The cp genome of S. persica encoded 131 genes, including 37 tRNA genes, eight rRNA genes and 86 protein-coding genes. Besides, S. persica contains 27 forward, 36 tandem and 19 palindromic repeats. The S. persica cp genome had 154 SSRs with the highest number in the LSC region. Complete cp genome comparisons showed an overall high degree of sequence resemblance between S. persica and related cp genomes. Some divergence was observed in the intergenic spaces of other species. Phylogenomic analyses of 60 shared genes indicated that S. persica formed a single clade with A. tetracantha with high bootstrap values. The family Salvadoraceae is closely related to Capparaceae and Petadiplandraceae rather than to Bataceae and Koberliniacaea.

CONCLUSION

The current genomic datasets provide pivotal genetic resources to determine the phylogenetic relationships, genome evolution and future genetic diversity-related studies of S. persica in complex angiosperm families.

The complete mitogenome sequence of the coral lily (Lilium pumilum) and the Lanzhou lily (Lilium davidii) in China

Background

The mitogenomes of higher plants are conserved. This study was performed to complete the mitogenome of two China Lilium species (Lilium pumilum Redouté and Lilium davidii var. unicolor (Hoog) cotton).

Methods

Genomic DNA was separately extracted from the leaves of L. pumilum and L. davidii in triplicate and used for sequencing. The mitogenome of Allium cepa was used as a reference. Genome assembly, annotation and phylogenetic tree were analyzed.

Results

The mitogenome of L. pumilum and L. davidii was 988,986 bp and 924,401 bp in length, respectively. There were 22 core protein-coding genes (including atp1, atp4, atp6, atp9, ccmB, ccmC, ccmFc, ccmFN1, ccmFN2, cob, cox3, matR, mttB, nad1, nad2, nad3, nad4, nad4L, nad5, nad6, nad7 and nad9), one open reading frame and one ribosomal protein-coding gene (rps12) in the mitogenomes. Compared with the A. cepa mitogenome, the coding sequence of the 24 genes and intergenic spacers in L. pumilum and L. davidii mitogenome contained 1,621 and 1,617 variable sites, respectively. In the phylogenetic tree, L. pumilum and L. davidii were distinct from A. cepa (NC_030100).

Conclusions

L. pumilum and L. davidii mitogenomes have far distances from other plants. This study provided additional information on the species resources of China Lilium.

Generation of Chloroplast Molecular Markers to Differentiate Sophora toromiro and Its Hybrids as a First Approach to Its Reintroduction in Rapa Nui (Easter Island)

Sophora toromiro is an endemic tree of Rapa Nui with religious and cultural relevance that despite being extinct in the wild, still persists in botanical gardens and private collections around the world. The authenticity of some toromiro trees has been questioned because the similarities among hybrid lines leads to misclassification of the species. The conservation program of toromiro has the objective of its reinsertion into Rapa Nui, but it requires the exact genotyping and certification of the selected plants in order to efficiently reintroduce the species. In this study, we present for the first time the complete chloroplast genome of S. toromiro and four other Sophora specimens, which were sequenced de-novo and assembled after mapping the raw reads to a chloroplast database. The length of the chloroplast genomes ranges from 154,239 to 154,473 bp. A total of 130-143 simple sequence repeats (SSR) loci and 577 single nucleotide polymorphisms (SNPs) were identified.

Comparative plastomic analysis and insights into the phylogeny of Salvia (Lamiaceae)

Salvia is the largest genus of Lamiaceae, with almost 1000 species, and has been divided into 11 subgenera. Salvia subg. Glutinaria, native to East Asia, is particularly important because of its potential medicinal value. However, the interspecific relationships of this subgenus have not been resolved and the plastomes of Salvia have rarely been studied. In the current study, we compared plastid genome structure and organization of 19 species of Salvia (14 newly sequenced and 5 previously published). Our comparative analysis showed that all Salvia plastomes examined have a quadripartite structure typical of most angiosperms and contain an identical set of 114 unique genes (80 protein-coding genes, 4 rRNA genes, and 30 tRNA genes). The plastome structure of all Salvia species is highly conserved like other Lamiaceae plastomes. Gene content, gene order, and GC content were highly similar in these plastomes. The inverted repeats/single copy region (IR/SC) boundaries of Salvia are highly conserved, and IR contraction only occurred in two species (Salvia mekongensis and S. rosmarinus). In Salvia, sequence divergence was higher in non-coding regions than in coding regions. We found that using large single copy (LSC) and small single copy regions (SSC) with exclusion of the rapidly evolving sites produced the highest resolution in phylogenetic analysis of Salvia, suggesting that using suitable informative sites to build trees is more conducive in phylogenetic research. This study assembled a powerful matrix data set for studying the phylogeny of Salvia, resolving the interspecific relationship of Salvia subg. Glutinaria. The newly sequenced plastid genomes will also enrich the plastome database of Salvia, providing the scientific basis for the development and utilization of germplasm resources of this large and important genus.

Characterization and Analysis of the Mitochondrial Genome of Common Bean (Phaseolus vulgaris) by Comparative Genomic Approaches

The common bean (Phaseolus vulgaris) is a major source of protein and essential nutrients for humans. To explore the genetic diversity and phylogenetic relationships of P. vulgaris, its complete mitochondrial genome (mitogenome) was sequenced and assembled. The mitogenome is 395,516 bp in length, including 31 unique protein-coding genes (PCGs), 15 transfer RNA (tRNA) genes, and 3 ribosomal RNA (rRNA) genes. Among the 31 PCGs, four genes (mttB, nad1, nad4L, and rps10) use ACG as initiation codons, which are altered to standard initiation codons by RNA editing. In addition, the termination codon CGA in the ccmFC gene is converted to UGA. Selective pressure analysis indicates that the ccmB, ccmFC, rps1, rps10, and rps14 genes were under evolutionary positive selection. The proportions of five amino acids (Phe, Leu, Pro, Arg, and Ser) in the whole amino acid profile of the proteins in each mitogenome can be used to distinguish angiosperms from gymnosperms. Phylogenetic analyses show that P. vulgaris is evolutionarily closer to the Glycininae than other leguminous plants. The results of the present study not only provide an important opportunity to conduct further genomic breeding studies in the common bean, they also provide valuable information for future evolutionary and molecular studies of leguminous plants.

The First Glimpse of Streptocarpus ionanthus (Gesneriaceae) Phylogenomics: Analysis of Five Subspecies' Chloroplast Genomes

Streptocarpus ionanthus (Gesneriaceae) comprise nine herbaceous subspecies, endemic to Kenya and Tanzania. The evolution of Str. ionanthus is perceived as complex due to morphological heterogeneity and unresolved phylogenetic relationships. Our study seeks to understand the molecular variation within Str. ionanthus using a phylogenomic approach. We sequence the chloroplast genomes of five subspecies of Str. ionanthus, compare their structural features and identify divergent regions. The five genomes are identical, with a conserved structure, a narrow size range (170 base pairs (bp)) and 115 unique genes (80 protein-coding, 31 tRNAs and 4 rRNAs). Genome alignment exhibits high synteny while the number of Simple Sequence Repeats (SSRs) are observed to be low (varying from 37 to 41), indicating high similarity. We identify ten divergent regions, including five variable regions (psbM, rps3, atpF-atpH, psbC-psbZ and psaA-ycf3) and five genes with a high number of polymorphic sites (rps16, rpoC2, rpoB, ycf1 and ndhA) which could be investigated further for phylogenetic utility in Str. ionanthus. Phylogenomic analyses here exhibit low polymorphism within Str. ionanthus and poor phylogenetic separation, which might be attributed to recent divergence. The complete chloroplast genome sequence data concerning the five subspecies provides genomic resources which can be expanded for future elucidation of Str. ionanthus phylogenetic relationships.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Quercus bawanglingensis Huang, Li et Xing, a Vulnerable Oak Tree in China

Quercus bawanglingensis Huang, Li et Xing, an endemic evergreen oak of the genus Quercus (Fagaceae) in China, is currently listed in the Red List of Chinese Plants as a vulnerable (VU) plant. No chloroplast (cp) genome information is currently available for Q. bawanglingensis, which would be essential for the establishment of guidelines for its conservation and breeding. In the present study, the cp genome of Q. bawanglingensis was sequenced and assembled into double-stranded circular DNA with a length of 161,394 bp. Two inverted repeats (IRs) with a total of 51,730 bp were identified, and the rest of the sequence was separated into two single-copy regions, namely, a large single-copy (LSC) region (90,628 bp) and a small single-copy (SSC) region (19,036 bp). The genome of Q. bawanglingensis contains 134 genes (86 protein-coding genes, 40 tRNAs and eight rRNAs). More forward (29) than inverted long repeats (21) are distributed in the cp genome. A simple sequence repeat (SSR) analysis showed that the genome contains 82 SSR loci, involving 84.15% A/T mononucleotides. Sequence comparisons among the nine complete cp genomes, including the genomes of Q. bawanglingensis, Q. tarokoensis Hayata (NC036370), Q. aliena var. acutiserrata Maxim. ex Wenz. (KU240009), Q. baronii Skan (KT963087), Q. aquifolioides Rehd. et Wils. (KX911971), Q. variabilis Bl. (KU240009), Fagus engleriana Seem. (KX852398), Lithocarpus balansae (Drake) A. Camus (KP299291) and Castanea mollissima Bl. (HQ336406), demonstrated that the diversity of SC regions was higher than that of IR regions, which might facilitate identification of the relationships within this extremely complex family. A phylogenetic analysis showed that Fagus engleriana and Trigonobalanus doichangensis form the basis of the produced evolutionary tree. Q. bawanglingensis and Q. tarokoensis, which belong to the group Ilex, share the closest relationship. The analysis of the cp genome of Q. bawanglingensis provides crucial genetic information for further studies of this vulnerable species and the taxonomy, phylogenetics and evolution of Quercus.

Comparative RNA editing profile of mitochondrial transcripts in cytoplasmic male sterile and fertile pigeonpea reveal significant changes at the protein level

RNA editing is a process which leads to post-transcriptional alteration of the nucleotide sequence of the corresponding mRNA molecule which may or may not lead to changes at the protein level. Apart from its role in providing variability at the transcript and protein levels, sometimes, such changes may lead to abnormal expression of the mitochondrial gene leading to a cytoplasmic male sterile phenotype. Here we report the editing status of 20 major mitochondrial transcripts in both male sterile (AKCMS11) and male fertile (AKPR303) pigeonpea genotypes. The validation of the predicted editing sites was done by mapping RNA-seq reads onto the amplified mitochondrial genes, and 165 and 159 editing sites were observed in bud tissues of the male sterile and fertile plant respectively. Among the resulting amino acid alterations, the most frequent one was the conversion of hydrophilic amino acids to hydrophobic. The alterations thus detected in our study indicates differential editing, but no major change in terms of the abnormal protein structure was detected. However, the above investigation provides an insight into the behaviour of pigeonpea mitochondrial genome in native and alloplasmic state and could hold clues in identification of editing factors and their role in adaptive evolution in pigeonpea.

The comparison of four mitochondrial genomes reveals cytoplasmic male sterility candidate genes in cotton

BACKGROUND

The mitochondrial genomes of higher plants vary remarkably in size, structure and sequence content, as demonstrated by the accumulation and activity of repetitive DNA sequences. Incompatibility between mitochondrial genome and nuclear genome leads to non-functional male reproductive organs and results in cytoplasmic male sterility (CMS). CMS has been used to produce F₁ hybrid seeds in a variety of plant species.

RESULTS

Here we compared the mitochondrial genomes (mitogenomes) of Gossypium hirsutum sterile male lines CMS-2074A and CMS-2074S, as well as their restorer and maintainer lines. First, we noticed the mitogenome organization and sequences were conserved in these lines. Second, we discovered the mitogenomes of 2074A and 2074S underwent large-scale substitutions and rearrangements. Actually, there were five and six unique chimeric open reading frames (ORFs) in 2074A and 2074S, respectively, which were derived from the recombination between unique repetitive sequences and nearby functional genes. Third, we found out four chimeric ORFs that were differentially transcribed in sterile line (2074A) and fertile-restored line.

CONCLUSIONS

These four novel and recombinant ORFs are potential candidates that confer CMS character in 2074A. In addition, our observations suggest that CMS in cotton is associated with the accelerated rates of rearrangement, and that novel expression products are derived from recombinant ORFs.

Characterization and Comparative Analysis of the Complete Chloroplast Genome of the Critically Endangered Species Streptocarpus teitensis (Gesneriaceae)

Streptocarpus teitensis (Gesneriaceae) is an endemic species listed as critically endangered in the International Union for Conservation of Nature (IUCN) red list of threatened species. However, the sequence and genome information of this species remains to be limited. In this article, we present the complete chloroplast genome structure of Streptocarpus teitensis and its evolution inferred through comparative studies with other related species. S. teitensis displayed a chloroplast genome size of 153,207 bp, sheltering a pair of inverted repeats (IR) of 25,402 bp each split by small and large single-copy (SSC and LSC) regions of 18,300 and 84,103 bp, respectively. The chloroplast genome was observed to contain 116 unique genes, of which 80 are protein-coding, 32 are transfer RNAs, and four are ribosomal RNAs. In addition, a total of 196 SSR markers were detected in the chloroplast genome of Streptocarpus teitensis with mononucleotides (57.1%) being the majority, followed by trinucleotides (33.2%) and dinucleotides and tetranucleotides (both 4.1%), and pentanucleotides being the least (1.5%). Genome alignment indicated that this genome was comparable to other sequenced members of order Lamiales. The phylogenetic analysis suggested that Streptocarpus teitensis is closely related to Lysionotus pauciflorus and Dorcoceras hygrometricum.

Application of Chloroplast Phylogenomics to Resolve Species Relationships Within the Plant Genus Amaranthus

Amaranthus species are an emerging and promising nutritious traditional vegetable food source. Morphological plasticity and poorly resolved dendrograms have led to the need for well resolved species phylogenies. We hypothesized that whole chloroplast phylogenomics would result in more reliable differentiation between closely related amaranth species. The aims of the study were therefore: to construct a fully assembled, annotated chloroplast genome sequence of Amaranthus tricolor; to characterize Amaranthus accessions phylogenetically by comparing barcoding genes (matK, rbcL, ITS) with whole chloroplast sequencing; and to use whole chloroplast phylogenomics to resolve deeper phylogenetic relationships. We generated a complete A. tricolor chloroplast sequence of 150,027 bp. The three barcoding genes revealed poor inter- and intra-species resolution with low bootstrap support. Whole chloroplast phylogenomics of 59 Amaranthus accessions increased the number of parsimoniously informative sites from 92 to 481 compared to the barcoding genes, allowing improved separation of amaranth species. Our results support previous findings that two geographically independent domestication events of Amaranthus hybridus likely gave rise to several species within the Hybridus complex, namely Amaranthus dubius, Amaranthus quitensis, Amaranthus caudatus, Amaranthus cruentus and Amaranthus hypochondriacus. Poor resolution of species within the Hybridus complex supports the recent and ongoing domestication within the complex, and highlights the limitation of chloroplast data for resolving recent evolution. The weedy Amaranthus retroflexus and Amaranthus powellii was found to share a common ancestor with the Hybridus complex. Leafy amaranth, Amaranthus tricolor, Amaranthus blitum, Amaranthus viridis and Amaranthus graecizans formed a stable sister lineage to the aforementioned species across the phylogenetic trees. This study demonstrates the power of next-generation sequencing data and reference-based assemblies to resolve phylogenies, and also facilitated the identification of unknown Amaranthus accessions from a local genebank. The informative phylogeny of the Amaranthus genus will aid in selecting accessions for breeding advanced genotypes to satisfy global food demand.

Phylogenic study of Lemnoideae (duckweeds) through complete chloroplast genomes for eight accessions

Background

Phylogenetic relationship within different genera of Lemnoideae, a kind of small aquatic monocotyledonous plants, was not well resolved, using either morphological characters or traditional markers. Given that rich genetic information in chloroplast genome makes them particularly useful for phylogenetic studies, we used chloroplast genomes to clarify the phylogeny within Lemnoideae.

Methods

DNAs were sequenced with next-generation sequencing. The duckweeds chloroplast genomes were indirectly filtered from the total DNA data, or directly obtained from chloroplast DNA data. To test the reliability of assembling the chloroplast genome based on the filtration of the total DNA, two methods were used to assemble the chloroplast genome of Landoltia punctata strain ZH0202. A phylogenetic tree was built on the basis of the whole chloroplast genome sequences using MrBayes v.3.2.6 and PhyML 3.0.

Results

Eight complete duckweeds chloroplast genomes were assembled, with lengths ranging from 165,775 bp to 171,152 bp, and each contains 80 protein-coding sequences, four rRNAs, 30 tRNAs and two pseudogenes. The identity of L. punctata strain ZH0202 chloroplast genomes assembled through two methods was 100%, and their sequences and lengths were completely identical. The chloroplast genome comparison demonstrated that the differences in chloroplast genome sizes among the Lemnoideae primarily resulted from variation in non-coding regions, especially from repeat sequence variation. The phylogenetic analysis demonstrated that the different genera of Lemnoideae are derived from each other in the following order: Spirodela, Landoltia, Lemna, Wolffiella, and Wolffia.

Discussion

This study demonstrates potential of whole chloroplast genome DNA as an effective option for phylogenetic studies of Lemnoideae. It also showed the possibility of using chloroplast DNA data to elucidate those phylogenies which were not yet solved well by traditional methods even in plants other than duckweeds.

The complete chloroplast genome of Primulina and two novel strategies for development of high polymorphic loci for population genetic and phylogenetic studies

Background

Primulina Hance is an emerging model for studying evolutionary divergence, adaptation and speciation of the karst flora. However, phylogenetic relationships within the genus have not been resolved due to low variation detected in the cpDNA regions. Chloroplast genomes can provide important information for phylogenetic and population genetic studies. Recent advances in next-generation sequencing (NGS) techniques greatly facilitate sequencing whole chloroplast genomes for multiple individuals. Consequently, novel strategies for development of highly polymorphic loci for population genetic and phylogenetic studies based on NGS data are needed.

Methods

For development of high polymorphic loci for population genetic and phylogenetic studies, two novel strategies are proposed here. The first protocol develops lineage-specific highly variable markers from the true high variation regions (Con_Seas) across whole cp genomes, instead of traditional noncoding regions. The pipeline has been integrated into a single perl script, and named "Con_Sea_Identification_and_PIC_Calculation". The second method assembles chloroplast fragments (poTs) and sub-super-marker (CpContigs) through our "SACRing" pipeline. This approach can fundamentally alter the strategies used in phylogenetic and population genetic studies based on cp markers, facilitating a transition from traditional Sanger sequencing to RAD-Seq. Both of these scripts are available at https://github.com/scbgfengchao/.

Results

Three complete Primulina chloroplast genomes were assembled from genome survey data, and then two novel strategies were developed to yield highly polymorphic markers. For experimental evaluation of the first protocol, a set of Primulina species were used for PCR amplification. The results showed that these newly developed markers are more variable than traditional ones, and seem to be a better choice for phylogenetic and population studies in Primulina. The second method was also successfully applied in population genetic studies of 21 individuals from three natural populations of Primulina.

Conclusions

These two novel strategies may provide a pathway for similar research in other non-model species. The newly developed high polymorphic loci in this study will promote further the phylogenetic and population genetic studies in Primulina and other genera of the family Gesneriaceae.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1067-z) contains supplementary material, which is available to authorized users.

Efficient modulation of photosynthetic apparatus confers desiccation tolerance in the resurrection plant Boea hygrometrica

Boea hygrometrica (B. hygrometrica) can tolerate severe desiccation and resume photosynthetic activity rapidly upon water availability. However, little is known about the mechanisms by which B. hygrometrica adapts to dehydration and resumes competence upon rehydration. Here we determine how B. hygrometrica deals with oxidative stress, excessive excitation/electron pressures as well as photosynthetic apparatus modulation during dehydration/rehydration. By measuring ROS generation and scavenging efficiency, we found that B. hygrometrica possesses efficient strategies to maintain cellular redox homeostasis. Transmission electron microscopy (TEM) analysis revealed a remarkable alteration of chloroplast architecture and plastoglobules (PGs) accumulation during dehydration/rehydration. Pulse-amplitude modulated (PAM) chlorophyll fluorescence measurements, P700 redox assay as well as chlorophyll fluorescence emission spectra analysis on leaves of B. hygrometrica during dehydration/rehydration were also performed. Results showed that the photochemical activity of PSII as well as photoprotective energy dissipation in PSII undergo gradual inactivation/activation during dehydration/rehydration in B. hygrometrica; PSI activity is relatively induced upon water deficit, and dehydration leads to physical interaction between PSI and LHCII. Furthermore, blue-native polyacrylamide gel electrophoresis (BN-PAGE) and immunoblot analysis revealed that the protein abundance of light harvesting complexes decrease markedly along with internal water deficit to restrict light absorption and attenuate electron transfer, resulting in limited light excitation and repressed photosynthesis. In contrast, many thylakoid proteins remain at a basal level even after full dehydration. Taken together, our study demonstrated that efficient modulation of cellular redox homeostasis and photosynthetic activity confers desiccation tolerance in B. hygrometrica.

Complete chloroplast genome sequences contribute to plant species delimitation: A case study of the Anemopaegma species complex

PREMISE OF THE STUDY

Bignoniaceae is an important component of neotropical forests and a model for evolutionary and biogeographical studies. A previous combination of molecular markers and morphological traits improved the phylogeny of the group. Here we demonstrate the value of next-generation sequencing (NGS) to assemble the chloroplast genome of eight Anemopaegma species and solve taxonomic problems.

METHODS

Three NGS platforms were used to sequence total DNA of Anemopaegma species. After genome assembly and annotation, we compared chloroplast genomes within Anemopaegma, with other Lamiales species, and the evolutionary rates of protein-coding genes using Tanaecium tetragonolobum as the outgroup. Phylogenetic analyses of Anemopaegma with different data sets were performed.

KEY RESULTS

Chloroplast genomes of Anemopaegma species ranged from 167,413 bp in A. foetidum to 168,987 bp in A. acutifolium ("typical" form). They exhibited a characteristic quadripartite structure with a large single-copy region (75,070-75,761 bp), a small single-copy region (12,766-12,817 bp) and a pair of inverted repeat regions (IRs) (39,480-40,481) encoding an identical set of 112 genes. An inversion of a fragment with ca. 8 kb, located in the IRs and containing the genes trnI-AAU, ycf2, and trnL-CAA, was observed in these chloroplast genomes when compared with those of other Lamiales.

CONCLUSIONS

Anemopaegma species have the largest genomes within the Lamiales possibly due to the large amount of repetitive sequences and IR expansion. Variation was higher in coding regions than in noncoding regions, and some genes were identified as markers for differentiation between species. The use of the entire chloroplast genome gave better phylogenetic resolution of the taxonomic groups. We found that two forms of A. acutifolium result from different maternal lineages.

Diversity of mitochondrial plastid DNAs (MTPTs) in seed plants

Mitochondrial plastid DNAs (MTPTs) refer to plastid-derived DNA fragments in mitochondrial genomes. While the MTPTs have been described for numerous species, its overall patterns have not been examined in details. Here, we carried out a systematic analysis of MTPTs among 73 plant species, including 28 algae, 1 liverwort, 2 moss, 1 lycophyte, 1 gymnosperm, 1 magnoliid, 12 monocots, 26 eudicots and 1 relic angiosperm Amborella trichopoda. A total of 300 MTPT gene clusters were found in 39 seed plants, which represented 144 MTPT gene cluster types. The detected MTPT gene clusters were evaluated in seven aspects, and they were found to be enriched particularly in monocots and asterids of eudicots. Some MTPT gene clusters were found to be shared by closely related species. All chloroplast genes were found in MTPTs, suggesting that there is no functional relevancy for genes that were transferred. However, after calculation of the frequency of the 115 chloroplast genes, five hot spots and three cold spots were discovered in chloroplast genome. In summary, this study demonstrated the high degree of diversity in MTPTs. The discovered MTPTs would facilitate the accurate assembly of chloroplast and mitochondrial genomes as well as the understanding of organelle genome evolution.

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.

Chloroplast Genome Analysis of Resurrection Tertiary Relict Haberlea rhodopensis Highlights Genes Important for Desiccation Stress Response

Haberlea rhodopensis is a paleolithic tertiary relict species, best known as a resurrection plant with remarkable tolerance to desiccation. When exposed to severe drought stress, H. rhodopensis shows an ability to maintain the structural integrity of its photosynthetic apparatus, which re-activates easily upon rehydration. We present here the results from the assembly and annotation of the chloroplast (cp) genome of H. rhodopensis, which was further subjected to comparative analysis with the cp genomes of closely related species. H. rhodopensis showed a cp genome size of 153,099 bp, harboring a pair of inverted repeats (IR) of 25,415 bp separated by small and large copy regions (SSC and LSC) of 17,826 and 84,443 bp. The genome structure, gene order, GC content and codon usage are similar to those of the typical angiosperm cp genomes. The genome hosts 137 genes representing 70.66% of the plastome, which includes 86 protein-coding genes, 36 tRNAs, and 4 rRNAs. A comparative plastome analysis with other closely related Lamiales members revealed conserved gene order in the IR and LSC/SSC regions. A phylogenetic analysis based on protein-coding genes from 33 species defines this species as belonging to the Gesneriaceae family. From an evolutionary point of view, a site-specific selection analysis detected positively selected sites in 17 genes, most of which are involved in photosynthesis (e.g., rbcL, ndhF, accD, atpE, etc.). The observed codon substitutions may be interpreted as being a consequence of molecular adaptation to drought stress, which ensures an evolutionary advantage to H. rhodopensis.

Comparative analysis of single nucleotide polymorphisms in the nuclear, chloroplast, and mitochondrial genomes in identification of phylogenetic association among seven melon (Cucumis melo L.) cultivars

A variety of melons are cultivated worldwide, and their specific biological properties make them an attractive model for molecular studies. This study aimed to investigate the single nucleotide polymorphisms (SNPs) from the mitochondrial, chloroplast, and nuclear genomes of seven melon accessions (Cucumis melo L.) to identify the phylogenetic relationships among melon cultivars with the Illumina HiSeq 2000 platform and bioinformatical analyses. The data showed that there were a total of 658 mitochondrial SNPs (207-295 in each), while there were 0-60 chloroplast SNPs among these seven melon cultivars, compared to the reference genome. Bioinformatical analysis showed that the mitochondrial tree topology was unable to separate the melon features, whereas the maximum parsimony/neighbor joining (MP/NJ) tree of the chloroplast SNPs could define melon features such as seed length, width, thickness, 100-seed weight, and type. SNPs of the nuclear genome were better than the mitochondrial and chloroplast SNPs in the identification of melon features. The data demonstrated the usefulness of mitochondrial, chloroplast, and nuclear SNPs in identification of phylogenetic associations among these seven melon cultivars.

Assembly and analysis of the complete Salix purpurea L. (Salicaceae) mitochondrial genome sequence

Plant mitochondrial (mt) genomes possess several complex features, including a variable size, a dynamic genome structure, and complicated patterns of gene loss and gain throughout evolutionary history. Studies of plant mt genomes can, therefore, provide unique insights into organelle evolution. We assembled the complete Salix purpurea L. mt genome by screening genomic sequence reads generated by a Roche-454 pyrosequencing platform. The pseudo-molecule obtained has a typical circular structure 598,970 bp long, with an overall GC content of 55.06%. The S. purpurea mt genome contains 52 genes: 31 protein-coding, 18 tRNAs, and three rRNAs. Eighteen tandem repeats and 404 microsatellites are distributed unevenly throughout the S. purpurea mt genome. A phylogenetic tree of 23 representative terrestrial plants strongly supports S. purpurea inclusion in the Malpighiales clade. Our analysis contributes toward understanding the organization and evolution of organelle genomes in Salicaceae species.

Analysis of the Complete Mitochondrial Genome Sequence of the Diploid Cotton Gossypium raimondii by Comparative Genomics Approaches

Cotton is one of the most important economic crops and the primary source of natural fiber and is an important protein source for animal feed. The complete nuclear and chloroplast (cp) genome sequences of G. raimondii are already available but not mitochondria. Here, we assembled the complete mitochondrial (mt) DNA sequence of G. raimondii into a circular genome of length of 676,078 bp and performed comparative analyses with other higher plants. The genome contains 39 protein-coding genes, 6 rRNA genes, and 25 tRNA genes. We also identified four larger repeats (63.9 kb, 10.6 kb, 9.1 kb, and 2.5 kb) in this mt genome, which may be active in intramolecular recombination in the evolution of cotton. Strikingly, nearly all of the G. raimondii mt genome has been transferred to nucleus on Chr1, and the transfer event must be very recent. Phylogenetic analysis reveals that G. raimondii, as a member of Malvaceae, is much closer to another cotton (G. barbadense) than other rosids, and the clade formed by two Gossypium species is sister to Brassicales. The G. raimondii mt genome may provide a crucial foundation for evolutionary analysis, molecular biology, and cytoplasmic male sterility in cotton and other higher plants.

The Chloroplast Genome of Utricularia reniformis Sheds Light on the Evolution of the ndh Gene Complex of Terrestrial Carnivorous Plants from the Lentibulariaceae Family

Lentibulariaceae is the richest family of carnivorous plants spanning three genera including Pinguicula, Genlisea, and Utricularia. Utricularia is globally distributed, and, unlike Pinguicula and Genlisea, has both aquatic and terrestrial forms. In this study we present the analysis of the chloroplast (cp) genome of the terrestrial Utricularia reniformis. U. reniformis has a standard cp genome of 139,725bp, encoding a gene repertoire similar to essentially all photosynthetic organisms. However, an exclusive combination of losses and pseudogenization of the plastid NAD(P)H-dehydrogenase (ndh) gene complex were observed. Comparisons among aquatic and terrestrial forms of Pinguicula, Genlisea, and Utricularia indicate that, whereas the aquatic forms retained functional copies of the eleven ndh genes, these have been lost or truncated in terrestrial forms, suggesting that the ndh function may be dispensable in terrestrial Lentibulariaceae. Phylogenetic scenarios of the ndh gene loss and recovery among Pinguicula, Genlisea, and Utricularia to the ancestral Lentibulariaceae cladeare proposed. Interestingly, RNAseq analysis evidenced that U. reniformis cp genes are transcribed, including the truncated ndh genes, suggesting that these are not completely inactivated. In addition, potential novel RNA-editing sites were identified in at least six U. reniformis cp genes, while none were identified in the truncated ndh genes. Moreover, phylogenomic analyses support that Lentibulariaceae is monophyletic, belonging to the higher core Lamiales clade, corroborating the hypothesis that the first Utricularia lineage emerged in terrestrial habitats and then evolved to epiphytic and aquatic forms. Furthermore, several truncated cp genes were found interspersed with U. reniformis mitochondrial and nuclear genome scaffolds, indicating that as observed in other smaller plant genomes, such as Arabidopsis thaliana, and the related and carnivorous Genlisea nigrocaulis and G. hispidula, the endosymbiotic gene transfer may also shape the U. reniformis genome in a similar fashion. Overall the comparative analysis of the U. reniformis cp genome provides new insight into the ndh genes and cp genome evolution of carnivorous plants from Lentibulariaceae family.

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

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

The Complete Chloroplast Genome Sequences of the Medicinal Plant Pogostemon cablin

Pogostemon cablin, the natural source of patchouli alcohol, is an important herb in the Lamiaceae family. Here, we present the entire chloroplast genome of P. cablin. This genome, with 38.24% GC content, is 152,460 bp in length. The genome presents a typical quadripartite structure with two inverted repeats (each 25,417 bp in length), separated by one small and one large single-copy region (17,652 and 83,974 bp in length, respectively). The chloroplast genome encodes 127 genes, of which 107 genes are single-copy, including 79 protein-coding genes, four rRNA genes, and 24 tRNA genes. The genome structure, GC content, and codon usage of this chloroplast genome are similar to those of other species in the family, except that it encodes less protein-coding genes and tRNA genes. Phylogenetic analysis reveals that P. cablin diverged from the Scutellarioideae clade about 29.45 million years ago (Mya). Furthermore, most of the simple sequence repeats (SSRs) are short polyadenine or polythymine repeats that contribute to high AT content in the chloroplast genome. Complete sequences and annotation of P. cablin chloroplast genome will facilitate phylogenic, population and genetic engineering research investigations involving this particular species.

Mitochondrial Genome Analysis of Wild Rice (Oryza minuta) and Its Comparison with Other Related Species

Oryza minuta (Poaceae family) is a tetraploid wild relative of cultivated rice with a BBCC genome. O. minuta has the potential to resist against various pathogenic diseases such as bacterial blight (BB), white backed planthopper (WBPH) and brown plant hopper (BPH). Here, we sequenced and annotated the complete mitochondrial genome of O. minuta. The mtDNA genome is 515,022 bp, containing 60 protein coding genes, 31 tRNA genes and two rRNA genes. The mitochondrial genome organization and the gene content at the nucleotide level are highly similar (89%) to that of O. rufipogon. Comparison with other related species revealed that most of the genes with known function are conserved among the Poaceae members. Similarly, O. minuta mt genome shared 24 protein-coding genes, 15 tRNA genes and 1 ribosomal RNA gene with other rice species (indica and japonica). The evolutionary relationship and phylogenetic analysis revealed that O. minuta is more closely related to O. rufipogon than to any other related species. Such studies are essential to understand the evolutionary divergence among species and analyze common gene pools to combat risks in the current scenario of a changing environment.

Complete Plastid Genome of the Recent Holoparasite Lathraea squamaria Reveals Earliest Stages of Plastome Reduction in Orobanchaceae

Plants from the family Orobanchaceae are widely used as a model to study different aspects of parasitic lifestyle including host–parasite interactions and physiological and genomic adaptations. Among the latter, the most prominent are those that occurred due to the loss of photosynthesis; they include the reduction of the photosynthesis-related gene set in both nuclear and plastid genomes. In Orobanchaceae, the transition to non-photosynthetic lifestyle occurred several times independently, but only one lineage has been in the focus of evolutionary studies. These studies included analysis of plastid genomes and transcriptomes and allowed the inference of patterns and mechanisms of genome reduction that are thought to be general for parasitic plants. Here we report the plastid genome of Lathraea squamaria, a holoparasitic plant from Orobanchaceae, clade Rhinantheae. We found that in this plant the degree of plastome reduction is the least among non-photosynthetic plants. Like other parasites, Lathraea possess a plastome with elevated absolute rate of nucleotide substitution. The only gene lost is petL, all other genes typical for the plastid genome are present, but some of them–those encoding photosystem components (22 genes), cytochrome b₆/f complex proteins (4 genes), plastid-encoded RNA polymerase subunits (2 genes), ribosomal proteins (2 genes), ccsA and cemA–are pseudogenized. Genes for cytochrome b₆/f complex and photosystems I and II that do not carry nonsense or frameshift mutations have an increased ratio of non-synonymous to synonymous substitution rates, indicating the relaxation of purifying selection. Our divergence time estimates showed that transition to holoparasitism in Lathraea lineage occurred relatively recently, whereas the holoparasitic lineage Orobancheae is about two times older.

Infraspecific variation within and across complete organellar genomes and nuclear ribosomal repeats in a moss

Bryophytes (mosses, liverworts, and hornworts) are diverse and ecologically and evolutionarily significant yet genome scale data sets and analyses remain extremely sparse relative to other groups of plants, and are completely lacking at the infraspecific level. By sequencing the complete organellar genomes and nuclear ribosomal repeat from seven patches of a South American sub-Antarctic neo-endemic non-model moss, we present the first characterization of infraspecific polymorphism within and across the three genomic compartments for a bryophyte. Diversity within patches is accounted for by both intraindividual and interindividual variation for the nuclear ribosomal repeat and plastid genome, respectively. This represents the most extensive infraspecific genomic dataset generated for an early land plant lineage thus far and provides insight into relative rates of substitution between organellar genomes, including high rates of nonsynonymous to synonymous substitutions.