Gene duplication, transfer, and evolution in the chloroplast genome

Phylogeny, biogeography, and character evolution of Anaphalis (Gnaphalieae, Asteraceae)

The HAP clade, mainly including Helichrysum Mill, Anaphalis DC., and Pseudognaphalium Kirp., is a major component of tribe Gnaphalieae (Asteraceae). In this clade, Anaphalis represents the largest genus of Asian Gnaphalieae. The intergeneric relationships among Anaphalis and its related genera and the infrageneric taxonomy of this genus are complex and remain controversial. However, there are few studies that have focused on these issues. Herein, based on the current most comprehensive sampling of the HAP clade, especially Anaphalis, we conducted phylogenetic analyses using chloroplast (cp) genome and nuclear ribosomal DNA (nrDNA) to evaluate the relationships within HAP clade, test the monophyly of Anaphalis, and examine the infrageneric taxonomy of this genus. Meanwhile, the morphological characters were verified to determine the circumscription and infrageneric taxonomy system of Anaphalis. Additionally, the biogeographical history, diversification processes, and evolution of crucial morphological characters were estimated and inferred. Our phylogenetic analyses suggested that Anaphalis is polyphyletic because it nested with Helichrysum and Pseudognaphalium. Two and four main clades of Anaphalis were identified in cp genome and nrDNA trees, respectively. Compared with nrDNA trees, the cp genome trees were more effective for phylogenetic resolution. After comprehensively analyzing morphological and phylogenetic evidence, it was concluded that the achene surface ornamentation and leaf base showed less homoplasy and supported the two Anaphalis lineages that were inferred from cp genome. Our biogeographical analyses based on cp genome indicated that HAP clade underwent rapid diversification from late Miocene to Pliocene. The two Anaphalis lineages appeared to have originated in Africa, then spread to Western and Southern Asia, and subsequently moved into Southwestern China forming a diversity center. The dispersal patterns of the two Anaphalis lineages were different. One dispersed around the world, except in Africa and South America. The other one dispersed to Eastern and Southeastern Asia from the ancestral origin region.

Comparative analyses of complete chloroplast genomes reveal interspecific difference and intraspecific variation of Tripterygium genus

The genus Tripterygium was of great medicinal value and attracted much attention on the taxonomic study using morphological and molecular methods. In this study, we assembled 12 chloroplast genomes of Tripterygium to reveal interspecific difference and intraspecific variation. The sequence length (156,692-157,061 bp) and structure of Tripterygium were conserved. Comparative analyses presented abundant variable regions for further study. Meanwhile, we determined the ndhB gene under positive selection through adaptive evolution analysis. And the phylogenetic analyses based on 15 chloroplast genomes supported the monophyly of Tripterygium hypoglaucum and the potential sister relationship between Tripterygium wilfordii and Tripterygium regelii. Molecular dating analysis indicated that the divergence time within Tripterygium was approximately 5.99 Ma (95% HPD = 3.11-8.68 Ma). The results in our study provided new insights into the taxonomy, evolution process, and phylogenetic construction of Tripterygium using complete plastid genomes.

Chloroplast Genome Comparison and Phylogenetic Analysis of the Commercial Variety Actinidia chinensis 'Hongyang'

Actinidia chinensis 'Hongyang', also known as red yangtao (red heart kiwifruit), is a vine fruit tree native to China possessing significant nutritional and economic value. However, information on its genetic diversity and phylogeny is still very limited. The first chloroplast (cp) genome of A. chinensis 'Hongyang' cultivated in China was sequenced using de novo technology in this study. A. chinensis 'Hongyang' possesses a cp genome that spans 156,267 base pairs (bp), exhibiting an overall GC content of 37.20%. There were 132 genes that were annotated, with 85 of them being protein-coding genes, 39 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. A total of 49 microsatellite sequences (SSRs) were detected, mainly single nucleotide repeats, mostly consisting of A or T base repeats. Compared with 14 other species, the cp genomes of A. chinensis 'Hongyang' were biased towards the use of codons containing A/U, and the non-protein coding regions in the A. chinensis 'Hongyang' cpDNA showed greater variation than the coding regions. The nucleotide polymorphism analysis (Pi) yielded nine highly variable region hotspots, most in the large single copy (LSC) region. The cp genome boundary analysis revealed a conservative order of gene arrangement in the inverted repeats (IRs) region of the cp genomes of 15 Actinidia plants, with small expansions and contractions of the boundaries. Furthermore, phylogenetic tree indicated that A. chinensis 'Hongyang' was the closest relative to A. indochinensis. This research provides a useful basis for future genetic and evolutionary studies of A. chinensis 'Hongyang', and enriches the biological information of Actinidia species.

New Insights Into The Evolution of Chloroplast Genomes in Ochna Species (Ochnaceae, Malpighiales)

Ochnaceae DC. includes more than 600 species that exhibit potential values for environmental ecology, ornamental, pharmaceutical, and timber industries. Although studies on phylogeny and phytochemicals have been intensively conducted, chloroplast genome data of Ochnaceae species have not been fully explored. In this study, the next-generation sequencing method was used to sequence the chloroplast genomes of Ochna integerrima and Ochna serrulata which were 157 329 and 157 835 bp in length, respectively. These chloroplast genomes had a quadripartite structure and contained 78 protein-coding genes, 30 tRNAs, and 4 rRNAs. Comparative analysis revealed 8 hypervariable regions, including trnK_UUU-trnQ_UUG, rpoB-psbM, trnS_GGA-rps4, accD-psaI, rpl33-rps18, rpl14-rpl16, ndhF-trnL_UAG, and rps15-ycf1 among 6 Ochnaceae taxa. Additionally, there were shared and unique repeats among 6 examined chloroplast genomes. The notable changes were the loss of rpl32 in Ochna species and the deletion of rps16 exon 2 in O. integerrima compared to other taxa. This study is the first comprehensive comparative genomic analysis of complete chloroplast genomes of Ochna species and related taxa in Ochnaceae. Consequently, the current study provides initial results for further research on genomic evolution, population genetics, and developing molecular markers in Ochnaceae and related taxa.

The complete chloroplast genome sequences of six Hylotelephium species: Comparative genomic analysis and phylogenetic relationships

To evaluate the phylogenetic relationships between Hylotelephium and Orostachys, and to provide important information for further studies, we analyzed the complete chloroplast genomes of six Hylotelephium species and compared the sequences to those of published chloroplast genomes of congeneric species and species of the closely related genus, Orostachys. The total chloroplast genome length of nineteen species, including the six Hylotelephium species analyzed in this study and the thirteen Hylotelephium and Orostachys species analyzed in previous studies, ranged from 150,369 bp (O. minuta) to 151,739 bp (H. spectabile). Their overall GC contents were almost identical (37.7-37.8%). The chloroplast genomes of the nineteen species contained 113 unique genes comprising 79 protein-coding genes (PCGs), 30 transfer RNA genes (tRNAs), and four ribosomal RNA genes (rRNAs). Among the annotated genes, fourteen genes contained one intron, and two genes contained two introns. The chloroplast genomes of the nineteen Hylotelephium and Orostachys species had identical structures. Additionally, the large single copy (LSC), inverted repeat (IR), and small single copy (SSC) junction regions were conserved in the Hylotelephium and Orostachys species. The nucleotide diversity between the Hylotelephium chloroplast genomes was extremely low in all regions, and only one region showed a high Pi value (>0.03). In all nineteen chloroplast genomes, six regions had a high Pi value (>0.03). The phylogenetic analysis showed that the genus delimitation could not be clearly observed even in this study because Hylotelephium formed a paraphyly with subsect. Orostachys of the genus Orostachys. Additionally, the data supported the taxonomic position of Sedum taqeutii, which was treated as a synonym for H. viridescens in previous studies, as an independent taxon.

The complete chloroplast genome of Diplodiscus trichospermus and phylogenetic position of Brownlowioideae within Malvaceae

BACKGROUND

Malvaceae is an economically important plant family of 4,225 species in nine subfamilies. Phylogenetic relationships among the nine subfamilies have always been controversial, especially for Brownlowioideae, whose phylogenetic position remains largely unknown due to the lack of samples in previous analysis datasets. To greatly clarify the phylogenetic relationship of Malvaceae, we newly sequenced and assembled the plastome of Diplodiscus trichospermus taxonomically located in Brownlowioideae, and downloaded the allied genomes from public database to build a dataset covering all subfamily members of Malvaceae.

RESULTS

The annotation results showed that the plastome of Diplodiscus trichospermus has a typical quadripartite structure, comprising 112 unique genes, namely 78 protein-coding genes, 30 tRNA genes and 4 rRNA genes. The total length was 158,570 bp with 37.2% GC content. Based on the maximum likelihood method and Bayesian inference, a robust phylogenetic backbone of Malvaceae was reconstructed. The topology showed that Malvaceae was divided distinctly into two major branches which were previously recognized as Byttneriina and Malvadendrina. In the Malvadendrina clade, Malvoideae and Bombacoideae formed, as always, a close sister clade named as Malvatheca. Subfamily Helicteroideae occupied the most basal position and was followed by Sterculioideae which was sister to the alliance of Malvatheca, Brownlowioideae, Dombeyoideae, and Tilioideae. Brownlowioideae together with the clade comprising Dombeyoideae and Tilioideae formed a sister clade to Malvatheca. In addition, one specific conservation SSR and three specific palindrome sequences were observed in Brownlowioideae.

CONCLUSIONS

In this study, the phylogenetic framework of subfamilies in Malvaceae has been resolved clearly based on plastomes, which may contribute to a better understanding of the classification and plastome evolution for Malvaceae.

LysR-type transcriptional regulators: state of the art

The LysR-type transcriptional regulators (LTTRs) are DNA-binding proteins present in bacteria, archaea, and in algae. Knowledge about their distribution, abundance, evolution, structural organization, transcriptional regulation, fundamental roles in free life, pathogenesis, and bacteria-plant interaction has been generated. This review focuses on these aspects and provides a current picture of LTTR biology.

Phylogeny of Leontopodium (Asteraceae) in China-with a reference to plastid genome and nuclear ribosomal DNA

The infrageneric taxonomy system, species delimitation, and interspecies systematic relationships of Leontopodium remain controversial and complex. However, only a few studies have focused on the molecular phylogeny of this genus. In this study, the characteristics of 43 chloroplast genomes of Leontopodium and its closely related genera were analyzed. Phylogenetic relationships were inferred based on chloroplast genomes and nuclear ribosomal DNA (nrDNA). Finally, together with the morphological characteristics, the relationships within Leontopodium were identified and discussed. The results showed that the chloroplast genomes of Filago, Gamochaeta, and Leontopodium were well-conserved in terms of gene number, gene order, and GC content. The most remarkable differences among the three genera were the length of the complete chloroplast genome, large single-copy region, small single-copy region, and inverted repeat region. In addition, the chloroplast genome structure of Leontopodium exhibited high consistency and was obviously different from that of Filago and Gamochaeta in some regions, such as matk, trnK (UUU)-rps16, petN-psbM, and trnE (UUC)-rpoB. All the phylogenetic trees indicated that Leontopodium was monophyletic. Except for the subgeneric level, our molecular phylogenetic results were inconsistent with the previous taxonomic system, which was based on morphological characteristics. Nevertheless, we found that the characteristics of the leaf base, stem types, and carpopodium base were phylogenetically correlated and may have potential value in the taxonomic study of Leontopodium. In the phylogenetic trees inferred using complete chloroplast genomes, the subgen. Leontopodium was divided into two clades (Clades 1 and 2), with most species in Clade 1 having herbaceous stems, amplexicaul, or sheathed leaves, and constricted carpopodium; most species in Clade 2 had woody stems, not amplexicaul and sheathed leaves, and not constricted carpopodium.

Phylogeny of Trigonotis in China-with a special reference to its nutlet morphology and plastid genome

The genus Trigonotis comprises nearly 60 species mainly distributed in East and Southeast Asia. China has the largest number of Trigonotis species in the world, with a total of 44 species, of which 38 are endemic. Nutlet morphology is useful for the taxonomic delimitation of Trigonotis. However, there are still controversial circumscriptions of nutlet shape in some species. In previous studies, interspecies phylogenetic relationships were inferred using few DNA markers and very few taxa, which possibly led to erroneous or incomplete conclusions. In this study, the nutlet morphology of 39 Trigonotis taxa and the characteristics of 34 complete chloroplast genomes (29 taxa) were investigated and analyzed. Then, the phylogenetic relationships were discussed within this genus based on complete chloroplast genomes. To the best of our knowledge, this study is the first comprehensive analysis of nutlet morphology and complete chloroplast genome of Trigonotis. Based on nutlet morphology, Trigonotis can be divided into two groups: Group 1, hemispherical or oblique tetrahedron with carpopodiums, and Group 2, inverted tetrahedron without carpopodiums. The chloroplast genome of Trigonotis exhibited a typical quadripartite structure, including 84-86 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes, with a total length of 147,247-148,986 bp. Genes in the junctions were well conserved in Trigonotis, similar to those in other Boraginaceae s.str. species. Furthermore, Trigonotis chloroplast genomes showed relatively high diversity, with more conserved genic regions than intergenic regions; in addition, we detected 14 hot spots (Pi > 0.005) in non-coding regions. Phylogenetic analyses based on chloroplast genome data identified highly resolved relationships between Trigonotis species. Specifically, Trigonotis was divided into two clades with strong support: one clade included species with hemispherical or oblique tetrahedron nutlets with carpopodiums and bracts, whereas the other clade included species with inverted tetrahedron nutlets without carpopodiums or bracts. Our results may inform future taxonomic, phylogenetic, and evolutionary studies on Boraginaceae.

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.

The complete chloroplast genome sequences of eight Orostachys species: Comparative analysis and assessment of phylogenetic relationships

We analyzed the complete chloroplast genomes of eight Orostachys species and compared the sequences to those of published chloroplast genomes of the congeneric and closely related genera, Meterostachys and Hylotelephium. The total chloroplast genome length of thirteen species, including the eight species analyzed in this study and the five species analyzed in previous studies, ranged from 149,860 (M. sikokianus) to 151,707 bp (H. verticillatum). The overall GC contents of the genomes were almost identical (37.6 to 37.8%). The thirteen chloroplast genomes each contained 113 unique genes comprising 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Among the annotated genes, sixteen genes contained one or two introns. Although the genome structures of all Orostachys and Hylotelephium species were identical, Meterostachys differed in structure due to a relatively large gene block (trnS-GCU-trnS-GGA) inversion. The nucleotide diversity among the subsect. Orostachys chloroplast genomes was extremely low in all regions, and among the subsect. Appendiculatae, genus Orostachys, and all thirteen chloroplast genomes showed high values of Pi (>0.03) in one, five, or three regions. The phylogenetic analysis showed that Orostachys formed polyphyly, and subsect. Orostachys and Appendiculatae were clustered with Hylotelephium and Meterostachys, respectively, supporting the conclusion that each subsection should be considered as an independent genus. Furthermore, the data supported the taxonomic position of O. margaritifolia and O. iwarenge f. magnus, which were treated as synonyms for O. iwarenge in a previous study, as independent taxa. Our results suggested that O. ramosa and O. japonica f. polycephala were individual variations of O. malacophylla and O. japonica, respectively. The exact taxonomic position of O. latielliptica and the phylogenetic relationship among the three species, O. chongsunensis, O. malacophylla and O. ramosa, should be a topic of future study.

Comparative Analysis of Bacillariophyceae Chloroplast Genomes Uncovers Extensive Genome Rearrangements Associated with Speciation

The Bacillariophyceae is a species-rich, ecologically significant class of Bacillariophyta. Despite their critical importance in marine ecosystems as primary producers and in the development of harmful algal blooms (HABs), taxonomic research on Bacillariophyceae species has been hindered because of their limited morphological features, plasticity of morphologies, and the low resolution of common molecular markers. Hence molecular markers with improved resolution are urgently needed. Organelle genomes, which can be constructed efficiently with the recent development of high throughput DNA sequencing technologies and the advancement of bioinformatics tools, have been proposed as super barcodes for their higher resolution for distinguishing different species and intra-species genomic variations. In this study, we tested the value of full-length chloroplast genomes (cpDNAs) as super barcodes for distinguishing diatom species, by constructing cpDNAs of 11 strains of the class Bacillariophyceae, including Nitzschia ovalis, Nitzschia traheaformis, Cylindrotheca spp., Psammodictyon constrictum, Bacillaria paxillifer, two strains of Haslea tsukamotoi, Haslea avium, Navicula arenaria, and Pleurosigma sp. Comparative analysis of cpDNAs revealed that cpDNAs were not only adequate for resolving different species, but also for enabling recognition of high levels of genome rearrangements between cpDNAs of different species, especially for species of the genera Nitzschia, Cylindrotheca, Navicula and Haslea. Additionally, comparative analysis suggested that the positioning of species in the genus Haslea should be transferred to the genus Navicula. Chloroplast genome-based evolutionary analysis suggested that the Bacillariophyceae species first appeared during the Cretaceous period and the diversity of species rose after the mass extinction about 65 Mya. This study highlighted the value of cpDNAs in research on the biodiversity and evolution of Bacillariophyceae species, and, with the construction of more cpDNAs representing additional genera, deeper insight into the biodiversity and evolutionary relationships of Bacillariophyceae species will be gained.

Characteristics and comparative analysis of Mesona chinensis Benth chloroplast genome reveals DNA barcode regions for species identification

Mesona chinensis Benth (MCB) is an important medicinal and edible plant in Southern China and Southeast Asian countries. Chloroplast (cp) genome is usually used for plant phylogeny, species identification, and chloroplast genetic engineering. To characterize the cp genome and determine the evolutionary position and perform the genetic diversity analysis of MCB, we sequence and characterize the MCB cp genome. The results show that the cp genome of MCB is a single circular molecule with a length of 152,635 bp. It is a typical quadripartite structure, comprising a large single-copy region (LSC, 83,514 bp) and a small single-copy region (SSC, 17,751 bp) separated by two inverted repeat regions (IRs, 51,370 bp). It encodes 129 unique genes, including 84 protein-coding genes (PCGs), 37 transfer RNAs (tRNAs), and 8 ribosomal RNAs (rRNAs). Altogether 127 simple sequence repeats (SSRs) are identified in the MCB cp genome with 86.61% of mononucleotide repeats. Phylogenetic analysis reveals that MCB is most closely related to Ocimum basilicum based on the whole cp genomes. Several highly divergent regions are found, such as trnH_psbA, rps16_trnQ, trnS_trnG, trnE_trnT, psaA_ycf3, rpl32_trnL, ccsA_ndhD, ndhG_ndhI, and rps15_ycf1, which can be proposed for use as DNA barcode regions. Genetic diversity analysis unveils a relatively narrow genetic basis of MCB germplasm resources. Therefore, the innovative breeding of MCB is very urgent and necessary in future research.

A Chromosome-Level Genome Assembly of the European Beech (Fagus sylvatica) Reveals Anomalies for Organelle DNA Integration, Repeat Content and Distribution of SNPs

The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.

The complete chloroplast genome sequence of Clerodendranthus spicatus, a medicinal plant for preventing and treating kidney diseases from Lamiaceae family

BACKGROUND

Clerodendranthus spicatus (Thunb.) C. Y. Wu ex H. W. Li is one of the most important medicines for the treatment of nephrology in the southeast regions of China. To understand the taxonomic classification of Clerodendranthus species and identify species discrimination markers, we sequenced and characterized its chloroplast genome in the current study.

METHODS AND RESULTS

Total genomic DNA were isolated from dried leaves of C. spicatus and sequenced using an Illumina sequencing platform. The data were assembled and annotated by the NOVOPlasty software and CpGAVAS2 web service. The complete chloroplast genome of C. spicatus was 152,155 bp, including a large single-copy region of 83,098 bp, a small single-copy region of 17,665 bp, and a pair of inverted repeat regions of 25,696 bp. The Isoleucine codons are the most abundant, accounting for 4.17% of all codons. The codons of AUG, UUA, and AGA demonstrated a high degree of usage bias. Twenty-eight simple sequence repeats, thirty-six tandem repeats, and forty interspersed repeats were identified. The distribution of the specific rps19, ycf1, rpl2, trnH, psbA genes were analyzed. Analysis of the genetic distance of the intergenic spacer regions shows that ndhG-ndhI, accD-psaI, rps15-ycf1, rpl20-clpP, ccsA-ndhD regions have high K2p values. Phylogenetic analysis showed that C. spicatu is closely related to two Lamiaceae species, Tectona grandis, and Glechoma longituba.

CONCLUSIONS

In this study, we sequenced and characterized the chloroplast genome of C. spicatus. Phylogenomic analysis has identified species closely related to C. spicatus, which represent potential candidates for the development of drugs improving renal functions.

The complete chloroplast genome sequence of Rubus hirsutus Thunb. and a comparative analysis within Rubus species

Rubus hirsutus is a type of tonifying kidney-essence herb that belongs to the Rosaceae family, and has been commonly used to treat multiple diseases, such as polyuria, impotence, and infertility. In this study, we determined the complete chloroplast sequence of R. hirsutus and conduced a comparative analysis within the genus Rubus. The assembled chloroplast (cp.) genome is 156,380 bp in length with a GC content of 37.0% and shares a conserved quadripartite structure within the other cp. genomes in this genus. A total of 132 unique genes were annotated in the cp. genome of R. hirsutus, which contained 87 protein-coding genes, 37 tRNAs, and eight rRNAs. Seventeen duplicated genes were identified in the inverted repeats region. Furthermore, 70 simple sequence repeats and 35 long repeats were detected in total in the R. hirsutus chloroplast genome. Eight mutational hotspots were identified in the cp. genome of this species with higher nucleotide variations in non-coding regions than those of coding regions. Furthermore, the gene order, codon usage, and repeat sequence distribution were highly consistent in Rubus according to the results of a comparative analysis. A phylogenetic analysis indicated that there was a sister relationship between R. hirsutus and R. chingii. Overall, the complete chloroplast genome of R. hirsutus and the comparative analysis will help to further the evolutionary study, conservation, phylogenetic reconstruction, and development of molecular barcodes for the genus Rubus.

Codon usage pattern and genetic diversity in chloroplast genomes of Panicum species

Exploring the molecular identities and the genetic diversity of a plant species is crucial in figuring out the evolutionary pressure of genes as well as in molecular breeding application. Nineteen chloroplast genomes of Panicum species in the National Center for Biotechnology Information database were downloaded and analyzed. The base composition, the effective number of codons, the relative synonymous codon usage, the codon bias index and the codon adaptation index of all genes in all chloroplast genomes, as well as the correlation coefficient among them, were calculated and discussed. The correspondence analysis and the clustering characteristics among nineteen genomes base on the relative synonymous codon usage values of nineteen chloroplast genomes were calculated and analyzed. In order to figuring out the evolutionary diversity of certain genes, the codon usage pattern of forty-one typical genes were separately counted and compared. Summations of their standard deviations were considered to evaluate their genetic diversities. The results of codon usage pattern showed that all genes were obvious AU-rich ones in chloroplast genomes of Panicum species, revealing that the natural selection was the main factor that influenced their evolutionary process. The correspondence and clustering analysis among nineteen chloroplast genomes showed that the overall evolutionary differences among them were not significant. However, the analysis on the genetic diversity of tyical genes showed that the degrees of diversity are different, and that the shorter sequences are more prone to instability. These findings would improve our understanding on the evolution of chloroplast genomes of Panicum species and be useful for further study on their evolutionary phenomenon.

The complete chloroplast genome sequence of Casearia glomerata Roxb. (Malpighiales; Salicaceae) from Fujian, China

Casearia glomerata Roxb. is classified in Salicaceae and has a high economic value. Herein, we report the complete chloroplast (cp) genome of C. glomerata using Illumina sequence data. The cp genome is 156,809 bp in length and contains a large single-copy (LSC) region of 84,888 bp and a small single-copy (SSC) region of 17,039 bp separated by two inverted repeat (IR) regions of 27,441 bp each. It contained a total of 123 genes, with an overall GC content of 36.81%. The phylogenetic analysis showed that C. glomerata is closely related to Casearia velutina. This study provides important sequence information for species identification and its phylogenetic relationship in the Salicaceae.

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.

Comparative analysis of chloroplast genomes indicated different origin for Indian tea (Camellia assamica cv TV1) as compared to Chinese tea

Based upon the morphological characteristics, tea is classified botanically into 2 main types i.e. Assam and China, which are morphologically very distinct. Further, they are so easily pollinated among themselves, that a third category, Cambod type is also described. Although the general consensus of origin of tea is India, Burma and China adjoining area, yet specific origin of China and Assam type tea are not yet clear. Thus, we made an attempt to understand the origin of Indian tea through the comparative analysis of different chloroplast (cp) genomes under the Camellia genus by performing evolutionary study and comparing simple sequence repeats (SSRs) and codon usage distribution patterns among them. The Cp genome based phylogenetic analysis indicated that Indian Tea, TV1 formed a different group from that of China tea, indicating that TV1 might have undergone different domestications and hence owe different origins. The simple sequence repeats (SSRs) analysis and codon usage distribution patterns also supported the clustering order in the cp genome based phylogenetic tree.

The complete plastome of Sorbaria kirilowii: genome structure, comparative analysis, and phylogenetic implications

Sorbaria kirilowii is a deciduous perennial admired for its showy white blossoms. Though of importance for horticultural purposes, the plastomic study concerning this species is still lacking. Here, the plastome of S. kirilowii was de novo assembled using the high-throughput sequencing data. The complete plastome assembly of S. kirilowii was 160,810 bp in length, with a GC content of 36.03%. It featured a typical quadripartite structure, containing a pair of inverted repeats (IRs; 26,338 bp) separated by a large single-copy (LSC; 88,762 bp) and a small single-copy (SSC, 19,372 bp). In total, 132 genes were annotated in the plastome, including 87 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. Furthermore, 63 SSRs, most of which were AT-rich, were identified in the cp genome of S. kirilowii. 71.7% of the cpSSRs were shown to be located in the intergenic regions. In addition, 49 repeats of varying sizes and types were also identified in the plastome. Through comparison, eight divergence hotspots were identified between the plastome of S. kirilowii and S. sorbifolia var. stellipila. These variable regions could potentially be developed into molecular markers for species delimitation or phylogenetics in future studies. We re-investigated the relationship among 17 Rosaceae species using the plastomic sequences, and S. kirilowii was shown to be a sister to S. sorbifolia var. stellipila. Overall, this study provides plastomic resources which could facilitate marker development and phylogenomics of Rosaceae.

Identification and Characterization of γ-Ray-Induced Mutations in Rice Cytoplasmic Genomes by Whole-Genome Sequencing

Chloroplasts and mitochondria are semi-autonomous organelles and have their own genomes (cytoplasmic genomes). Physical radiations (e.g., γ-rays) have been widely used in artificial mutation induction for plant germplasm enhancement and for breeding new cultivars. However, little is known at the genomic level about which kind of cytoplasmic mutations and/or characteristics could be induced in plants. The present study aimed to investigate the type, number, and distribution of inheritable cytoplasmic mutations induced by γ-rays in rice (Oryza sativa L.). Six plants were selected from the 2nd generation (M2) populations after γ-ray (137Cs) irradiation of the rice cultivar Nipponbare, 2 each for the 3 irradiation doses (150, 250, and 350 Gy), and their genomes were sequenced on an Illumina platform. Together with the whole-genome sequencing data of 3 external Nipponbare control plants, single-base substitutions (SBSs) and insertions/deletions (InDels) in chloroplast (cp) and mitochondrial (mt) genomes were identified and analyzed in-depth using bioinformatic tools. The majority of SBSs and InDels identified were background mutations in the 6 M2 plants, and the number of induced mutations varied greatly among the plants. Most induced mutations were present in a heterogeneous state, reflecting the fact that multiple cp and mt copies existed in the progenitor cells. The induced mutations were distributed in different genomic regions in the 6 M2 plants, including exonic regions, but none of them was predicted to cause nonsynonymous mutations or frameshifts. Our study thus revealed, at the genomic level, characteristics of cytoplasmic mutations induced by γ-rays in rice.

The chloroplast genome of the lichen-symbiont microalga Trebouxia sp. Tr9 (Trebouxiophyceae, Chlorophyta) shows short inverted repeats with a single gene and loss of the rps4 gene, which is encoded by the nucleus

The Trebouxiophyceae is the class of Chlorophyta algae from which the highest number of chloroplast genome (cpDNA) sequences has been obtained. Several species in this class participate in symbioses with fungi to form lichens. However, no cpDNA has been obtained from any Trebouxia lichen-symbiont microalgae, which are present in approximately half of all lichens. Here, we report the sequence of the completely assembled cpDNA from Trebouxia sp. TR9 and a comparative study with other Trebouxio-phyceae. The organization of the chloroplast genome of Trebouxia sp. TR9 has certain features that are unusual in the Trebouxiophyceae and other green algae. The most remarkable characteristics are the presence of long intergenic spacers, a quadripartite structure with short inverted repeated sequences (IRs), and the loss of the rps4 gene. The presence of long intergenic spacers accounts for a larger cpDNA size in comparison to other closely related Trebouxiophyceae. The IRs, which were thought to be lost in the Trebouxiales, are distinct from most of cpDNAs since they lack the rRNA operon and uniquely includes the rbcL gene. The functional transfer of the rps4 gene to the nuclear genome has been confirmed by sequencing and examination of the gene architecture, which includes three spliceosomal introns as well as the verification of the presence of the corresponding transcript. This is the first documented transfer of the rps4 gene from the chloroplast to the nucleus among Viridiplantae. Additionally, a fairly well-resolved phylogenetic reconstruction, including Trebouxia sp. TR9 along with other Trebouxiophyceae, was obtained based on a set of conserved chloroplast genes.

Intergeneric Relationships within the Family Salicaceae s.l. based on Plastid Phylogenomics

Many Salicaceae s.l. plants are recognized for their important role in the production of products such as wood, oils, and medicines, and as a model organism in life studies. However, the difference in plastid sequence, phylogenetic relationships, and lineage diversification of the family Salicaceae s.l. remain poorly understood. In this study, we compare 24 species representing 18 genera of the family. Simple sequence repeats (SSRs) are considered effective molecular markers for plant species identification and population genetics. Among them, a total of 1798 SSRs were identified, among which mononucleotide repeat was the most common with 1455 accounts representing 80.92% of the total. Most of the SSRs are located in the non-coding region. We also identified five other types of repeats, including 1750 tandems, 434 forward, 407 palindromic, 86 reverse, and 30 complementary repeats. The species in Salicaceae s.l. have a conserved plastid genome. Each plastome presented a typical quadripartite structure and varied in size due to the expansion and contraction of the inverted repeat (IR) boundary, lacking major structural variations, but we identified six divergence hotspot regions. We obtained phylogenetic relationships of 18 genera in Salicaceae s.l. and the 24 species formed a highly supported lineage. Casearia was identified as the basal clade. The divergence time between Salicaceae s.l. and the outgroup was estimated as ~93 Mya; Salix, and Populus diverged around 34 Mya, consistent with the previously reported time. Our research will contribute to a better understanding of the phylogenetic relationships among the members of the Salicaceae s.l.

The complete chloroplast genome sequences of four Viola species (Violaceae) and comparative analyses with its congeneric species

We report the complete chloroplast genomes of four Viola species (V. mirabilis, V. phalacrocarpa, V. raddeana, and V. websteri) and the results of a comparative analysis between these species and the published plastid genome of the congeneric species V. seoulensis. The total genome length of the five Viola species, including the four species analyzed in this study and the species analyzed in the previous study, ranged from 156,507 (V. seoulensis) to 158,162 bp (V. mirabilis). The overall GC contents of the genomes were almost identical (36.2-36.3%). The five Viola plastomes each contained 111 unique genes comprising 77 protein-coding genes, 30 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes. Among the annotated genes, 16 contained one or two introns. Based on the results of a chloroplast genome structure comparison using MAUVE, all five Viola plastomes were almost identical. Additionally, the large single copy (LSC), inverted repeat (IR), and small single copy (SSC) junction regions were conserved among the Viola species. A total of 259 exon, intron, and intergenic spacer (IGS) fragments were compared to verify the divergence hotspot regions. The nucleotide diversity (Pi) values ranged from 0 to 0.7544. The IR region was relatively more conserved than the LSC and SSC regions. The Pi values in ten noncoding regions were relatively high (>0.03). Among these regions, all but rps19-trnH, petG-trnW, rpl16-rps3, and rpl2-rpl23 represent useful molecular markers for phylogenetic studies and will be helpful to resolve the phylogenetic relationships of Viola. The phylogenetic tree, which used 76 protein-coding genes from 21 Malpighiales species and one outgroup species (Averrhoa carambola), revealed that Malpighiales is divided into five clades at the family level: Erythroxylaceae, Chrysobalanaceae, Euphorbiaceae, Salicaceae, and Violaceae. Additionally, Violaceae was monophyletic, with a bootstrap value of 100% and was divided into two subclades.

Enlarged and highly repetitive plastome of Lagarostrobos and plastid phylogenomics of Podocarpaceae

Podocarpaceae is the largest family in cupressophytes (conifers II), but its plastid genomes (plastomes) are poorly studied, with plastome data currently existing for only four of the 19 Podocarpaceous genera. In this study, we sequenced and assembled the complete plastomes from representatives of eight additional genera, including Afrocarpus, Dacrydium, Lagarostrobos, Lepidothamnus, Pherosphaera, Phyllocladus, Prumnopitys, and Saxegothaea. We found that Lagarostrobos, a monotypic genus native to Tasmania, has the largest plastome (151,496 bp) among any cupressophytes studied to date. Plastome enlargement in Lagarostrobos coincides with increased intergenic spacers, repeats, and duplicated genes. Among the Podocarpaceae, Lagarostrobos has the most rearranged plastome, but its substitution rates are modest. Plastid phylogenomic analyses based on 81 plastid genes clarify the positions of previously conflicting Podocarpaceous genera. Tree topologies firmly support the division of Podocarpaceae into two sister clades: (1) the Prumnopityoid clade and (2) the clade containing Podocarpoid, Dacrydioid, Pherosphaera, and Saxegothaea. The Phyllocladus is nested within the Podocarpaceae, thus familial status of the monotypic Phyllocladaceae is not supported.

Identification of plant species using variable length chloroplast DNA sequences

The correct identification of species in the highly divergent group of plants is crucial for several forensic investigations. Previous works had difficulties in the establishment of a rapid and robust method for the identification of plants. For instance, DNA barcoding requires the analysis of two or three different genomic regions to attain reasonable levels of discrimination. Therefore, new methods for the molecular identification of plants are clearly needed. Here we tested the utility of variable-length sequences in the chloroplast DNA (cpDNA) as a way to identify plant species. The SPInDel (Species Identification by Insertions/Deletions) approach targets hypervariable genomic regions that contain multiple insertions/deletions (indels) and length variability, which are found interspersed with highly conserved regions. The combination of fragment lengths defines a unique numeric profile for each species, allowing its identification. We analysed more than 44,000 sequences retrieved from public databases belonging to 206 different plant families. Four target regions were identified as suitable for the SPInDel concept: atpF-atpH, psbA-trnH, trnL CD and trnL GH. When considered alone, the discrimination power of each region was low, varying from 5.18% (trnL GH) to 42.54% (trnL CD). However, the discrimination power reached more than 90% when the length of some of these regions is combined. We also observed low diversity in intraspecific data sets for all target regions, suggesting they can be used for identification purposes. Our results demonstrate the utility of the SPInDel concept for the identification of plants.

Chloroplast genomic resources for phylogeny and DNA barcoding: a case study on Fritillaria

The genus Fritillaria comprises approximately 130 perennial herbaceous species. In the Pharmacopoeia of the People’s Republic of China, the bulbs of 11 Fritillaria species are used in Chinese herbal medicines. However, the traditional methods of morphological classification cannot accurately identify closely related species of Fritillaria. Previous studies have attempted to identify these species with universal molecular markers, but insufficient phylogenetic signal was available. In this study, the complete chloroplast genomes of eight Fritillaria species were compared. The length of the eight Fritillaria chloroplast genomes ranges from 151,009 bp to 152,224 bp. A total of 136 SSR loci were identified, including 124 polymorphic SSR loci. For large repeat sequences, 108 repeat loci and four types of repeats were observed. Ten highly variable regions were identified as potential molecular markers. These SSRs, large repeat sequences and highly variable regions provide important information for the development of genetic markers and DNA fingerprints. Phylogenetic analyses showed that the topological structures of all data sets (except the IR regions) were in complete agreement and well resolved. Overall, this study provides comprehensive chloroplast genomic resources, which will be valuable for future studies of evolution and species identification in Fritillaria.

A rare case of plastid protein-coding gene duplication in the chloroplast genome of Euglena archaeoplastidiata (Euglenophyta)

Gene duplication is an important evolutionary process that allows duplicate functions to diverge, or, in some cases, allows for new functional gains. However, in contrast to the nuclear genome, gene duplications within the chloroplast are extremely rare. Here, we present the chloroplast genome of the photosynthetic protist Euglena archaeoplastidiata. Upon annotation, it was found that the chloroplast genome contained a novel tandem direct duplication that encoded a portion of RuBisCO large subunit (rbcL) followed by a complete copy of ribosomal protein L32 (rpl32), as well as the associated intergenic sequences. Analyses of the duplicated rpl32 were inconclusive regarding selective pressures, although it was found that substitutions in the duplicated region, all non-synonymous, likely had a neutral functional effect. The duplicated region did not exhibit patterns consistent with previously described mechanisms for tandem direct duplications, and demonstrated an unknown mechanism of duplication. In addition, a comparison of this chloroplast genome to other previously characterized chloroplast genomes from the same family revealed characteristics that indicated E. archaeoplastidiata was probably more closely related to taxa in the genera Monomorphina, Cryptoglena, and Euglenaria than it was to other Euglena taxa. Taken together, the chloroplast genome of E. archaeoplastidiata demonstrated multiple characteristics unique to the euglenoid world, and has justified the longstanding curiosity regarding this enigmatic taxon.

Molecular polymorphism in Pistacia vera L. using non-coding regions of chloroplast DNA

The present study describes plastid DNA polymorphism and reports a comparative analysis of two non-coding cpDNA regions (trnC-trnD and atpB-rbcL) in pistachio. Seventeen different genotypes of domestic and wild pistachio from Iran, Syria, Turkey and America were sampled. Total genomic DNA was extracted and amplified with trnC-trnD and atpB-rbcL specific primers and then were sequenced. Phylogenetic relationships and depiction of phylogenetic trees were conducted. Cultivated genotypes of Pistacia vera were classified in a group regardless of their geographic location. P. vera was isolated from Sarakhs but they placed in the two close groups. Among cultivated genotypes, Jalab was separated from other cultivated genotypes. Pistacia Khinjuk was classified with Pistacia atlantica subsp. mutica. The findings confirm the common splitting hypothesis for commercial pistachio genotypes of the P. vera wild-type and also indicated the direct impact of Iranian genotypes in the evolutionary process of cultivated pistachios in other parts of the world. In conclusion it can be inferred that cultivated varieties of pistachio and P. vera var. sarakhs have the same origin, moreover genomic chloroplast could appropriately identify the interspecies relationships of pistachios.

Insight into infrageneric circumscription through complete chloroplast genome sequences of two Trillium species

Genomic events including gene loss, duplication, pseudogenization and rearrangement in plant genomes are valuable sources for exploring and understanding the process of evolution in angiosperms. The family Melanthiaceae is distributed in temperate regions of the Northern Hemisphere and divided into five tribes (Heloniadeae, Chionographideae, Xerophylleae, Melanthieae and Parideae) based on the molecular phylogenetic analyses. At present, complete chloroplast genomes of the Melanthiaceae have been reported from three species. In the previous genomic study of Liliales, atrnI-CAU gene duplication event was reported fromParis verticillata, a member of Parideae. To clarify the significant genomic events of the tribe Parideae, we analysed the complete chloroplast genome sequences of twoTrilliumspecies representing two subgenera:TrilliumandPhyllantherum InTrillium tschonoskii(subgenusTrillium), the circular double-stranded cpDNA sequence of 156 852 bp consists of two inverted repeat (IR) regions of 26 501 bp each, a large single-copy (LSC) region of 83 981 bp and a small single-copy (SSC) region of 19 869 bp. The chloroplast genome sequence ofT. maculatum(subgenusPhyllantherum) is 157 359 bp in length, consisting of two IRs (25 535 bp), one SSC (19 949 bp) and one LSC (86 340 bp), and is longer than that ofT. tschonoskii The results showed that the cpDNAs of Parideae are highly conserved across genome structure, gene order and contents. However, the chloroplast genome ofT. maculatumcontained a 3.4-kb inverted sequence betweenndhCandrbcLin the LSC region, and it was a unique feature for subgeneraPhyllantherum In addition, we found three different types of gene duplication in the intergenic spacer betweenrpl23andycf2containingtrnI-CAU, which were in agreement with the circumscription of subgenera and sections in Parideae excludingT. govanianum These genomic features provide informative molecular markers for identifying the infrageneric taxa ofTrilliumand improve our understanding of the evolution patterns of Parideae in Melanthiaceae.

A trnI_CAU triplication event in the complete chloroplast genome of Paris verticillata M.Bieb. (Melanthiaceae, Liliales)

The chloroplast is an essential plant organelle responsible for photosynthesis. Gene duplication, relocation, and loss in the chloroplast genome (cpDNA) are useful for exploring the evolution and phylogeny of plant species. In this study, the complete chloroplast genome of Paris verticillata was sequenced using the 454 sequencing system and Sanger sequencing method to trace the evolutionary pattern in the tribe Parideae of the family Melanthiaceae (Liliales). The circular double-stranded cpDNA of P. verticillata (157,379 bp) consists of two inverted repeat regions each of 28,373 bp, a large single copy of 82,726 bp, and a small single copy of 17,907 bp. Gene content and order are generally similar to the previously reported cpDNA sequences within the order Liliales. However, we found that trnI_CAU was triplicated in P. verticillata. In addition, cemA is suspected to be a pseudogene due to the presence of internal stop codons created by poly(A) insertion and single small CA repeats. Such changes were not found in previously examined cpDNAs of the Melanthiaceae or other families of the Liliales, suggesting that such features are unique to the tribe Parideae of Melanthiaceae. The characteristics of P. verticillata cpDNA will provide useful information for uncovering the evolution within Paris and for further research of plastid genome evolution and phylogenetic studies in Liliales.

Establishing oleaginous microalgae research models for consolidated bioprocessing of solar energy

Algal feedstock is the foundation of the emerging algal biofuel industry. However, few algae found in nature have demonstrated the combination of high biomass accumulation rate, robust oil yield and tolerance to environmental stresses, all complex traits that a large-scale, economically competitive production scheme demands. Therefore, untangling the intricate sub-cellular networks underlying these complex traits, in one or a series of carefully selected algal research models, has become an urgent research mission, which can take advantage of the emerging model oleaginous microalgae that have already demonstrated small, simple and tackleable genomes and the potential for large-scale open-pond cultivation. The revolutions in whole-genome-based technologies, coupled with systems biology, metabolic engineering and synthetic biology approaches, would enable the rational design and engineering of algal feedstock and help to fill the gaps between the technical and economical reality and the enormous potential of algal biofuels.

Insights from the complete chloroplast genome into the evolution of Sesamum indicum L

Sesame (Sesamum indicum L.) is one of the oldest oilseed crops. In order to investigate the evolutionary characters according to the Sesame Genome Project, apart from sequencing its nuclear genome, we sequenced the complete chloroplast genome of S. indicum cv. Yuzhi 11 (white seeded) using Illumina and 454 sequencing. Comparisons of chloroplast genomes between S. indicum and the 18 other higher plants were then analyzed. The chloroplast genome of cv. Yuzhi 11 contains 153,338 bp and a total of 114 unique genes (KC569603). The number of chloroplast genes in sesame is the same as that in Nicotiana tabacum, Vitis vinifera and Platanus occidentalis. The variation in the length of the large single-copy (LSC) regions and inverted repeats (IR) in sesame compared to 18 other higher plant species was the main contributor to size variation in the cp genome in these species. The 77 functional chloroplast genes, except for ycf1 and ycf2, were highly conserved. The deletion of the cp ycf1 gene sequence in cp genomes may be due either to its transfer to the nuclear genome, as has occurred in sesame, or direct deletion, as has occurred in Panax ginseng and Cucumis sativus. The sesame ycf2 gene is only 5,721 bp in length and has lost about 1,179 bp. Nucleotides 1–585 of ycf2 when queried in BLAST had hits in the sesame draft genome. Five repeats (R10, R12, R13, R14 and R17) were unique to the sesame chloroplast genome. We also found that IR contraction/expansion in the cp genome alters its rate of evolution. Chloroplast genes and repeats display the signature of convergent evolution in sesame and other species. These findings provide a foundation for further investigation of cp genome evolution in Sesamum and other higher plants.

Rates of molecular evolution and diversification in plants: chloroplast substitution rates correlate with species-richness in the Proteaceae

BACKGROUND

Many factors have been identified as correlates of the rate of molecular evolution, such as body size and generation length. Analysis of many molecular phylogenies has also revealed correlations between substitution rates and clade size, suggesting a link between rates of molecular evolution and the process of diversification. However, it is not known whether this relationship applies to all lineages and all sequences. Here, in order to investigate how widespread this phenomenon is, we investigate patterns of substitution in chloroplast genomes of the diverse angiosperm family Proteaceae. We used DNA sequences from six chloroplast genes (6278bp alignment with 62 taxa) to test for a correlation between diversification and the rate of substitutions.

RESULTS

Using phylogenetically-independent sister pairs, we show that species-rich lineages of Proteaceae tend to have significantly higher chloroplast substitution rates, for both synonymous and non-synonymous substitutions.

CONCLUSIONS

We show that the rate of molecular evolution in chloroplast genomes is correlated with net diversification rates in this large plant family. We discuss the possible causes of this relationship, including molecular evolution driving diversification, speciation increasing the rate of substitutions, or a third factor causing an indirect link between molecular and diversification rates. The link between the synonymous substitution rate and clade size is consistent with a role for the mutation rate of chloroplasts driving the speed of reproductive isolation. We find no significant differences in the ratio of non-synonymous to synonymous substitutions between lineages differing in net diversification rate, therefore we detect no signal of population size changes or alteration in selection pressures that might be causing this relationship.

Plant organelle proteomics: collaborating for optimal cell function

Organelle proteomics describes the study of proteins present in organelle at a particular instance during the whole period of their life cycle in a cell. Organelles are specialized membrane bound structures within a cell that function by interacting with cytosolic and luminal soluble proteins making the protein composition of each organelle dynamic. Depending on organism, the total number of organelles within a cell varies, indicating their evolution with respect to protein number and function. For example, one of the striking differences between plant and animal cells is the plastids in plants. Organelles have their own proteins, and few organelles like mitochondria and chloroplast have their own genome to synthesize proteins for specific function and also require nuclear-encoded proteins. Enormous work has been performed on animal organelle proteomics. However, plant organelle proteomics has seen limited work mainly due to: (i) inter-plant and inter-tissue complexity, (ii) difficulties in isolation of subcellular compartments, and (iii) their enrichment and purity. Despite these concerns, the field of organelle proteomics is growing in plants, such as Arabidopsis, rice and maize. The available data are beginning to help better understand organelles and their distinct and/or overlapping functions in different plant tissues, organs or cell types, and more importantly, how protein components of organelles behave during development and with surrounding environments. Studies on organelles have provided a few good reviews, but none of them are comprehensive. Here, we present a comprehensive review on plant organelle proteomics starting from the significance of organelle in cells, to organelle isolation, to protein identification and to biology and beyond. To put together such a systematic, in-depth review and to translate acquired knowledge in a proper and adequate form, we join minds to provide discussion and viewpoints on the collaborative nature of organelles in cell, their proper function and evolution.

Molecular evolution of rbcL in three gymnosperm families: identifying adaptive and coevolutionary patterns

Background

The chloroplast-localized ribulose-1, 5-biphosphate carboxylase/oxygenase (Rubisco), the primary enzyme responsible for autotrophy, is instrumental in the continual adaptation of plants to variations in the concentrations of CO₂. The large subunit (LSU) of Rubisco is encoded by the chloroplast rbcL gene. Although adaptive processes have been previously identified at this gene, characterizing the relationships between the mutational dynamics at the protein level may yield clues on the biological meaning of such adaptive processes. The role of such coevolutionary dynamics in the continual fine-tuning of RbcL remains obscure.

Results

We used the timescale and phylogenetic analyses to investigate and search for processes of adaptive evolution in rbcL gene in three gymnosperm families, namely Podocarpaceae, Taxaceae and Cephalotaxaceae. To understand the relationships between regions identified as having evolved under adaptive evolution, we performed coevolutionary analyses using the software CAPS. Importantly, adaptive processes were identified at amino acid sites located on the contact regions among the Rubisco subunits and on the interface between Rubisco and its activase. Adaptive amino acid replacements at these regions may have optimized the holoenzyme activity. This hypothesis was pinpointed by evidence originated from our analysis of coevolution that supported the correlated evolution between Rubisco and its activase. Interestingly, the correlated adaptive processes between both these proteins have paralleled the geological variation history of the concentration of atmospheric CO₂.

Conclusions

The gene rbcL has experienced bursts of adaptations in response to the changing concentration of CO₂ in the atmosphere. These adaptations have emerged as a result of a continuous dynamic of mutations, many of which may have involved innovation of functional Rubisco features. Analysis of the protein structure and the functional implications of such mutations put forward the conclusion that this evolutionary scenario has been possible through a complex interplay between adaptive mutations, often structurally destabilizing, and compensatory mutations. Our results unearth patterns of evolution that have likely optimized the Rubisco activity and uncover mutational dynamics useful in the molecular engineering of enzymatic activities.

Reviewers

This article was reviewed by Prof. Christian Blouin (nominated by Dr W Ford Doolittle), Dr Endre Barta (nominated by Dr Sandor Pongor), and Dr Nicolas Galtier.

Identification of new polymorphic regions and differentiation of cultivated olives (Olea europaea L.) through plastome sequence comparison

BACKGROUND

The cultivated olive (Olea europaea L.) is the most agriculturally important species of the Oleaceae family. Although many studies have been performed on plastid polymorphisms to evaluate taxonomy, phylogeny and phylogeography of Olea subspecies, only few polymorphic regions discriminating among the agronomically and economically important olive cultivars have been identified. The objective of this study was to sequence the entire plastome of olive and analyze many potential polymorphic regions to develop new inter-cultivar genetic markers.

RESULTS

The complete plastid genome of the olive cultivar Frantoio was determined by direct sequence analysis using universal and novel PCR primers designed to amplify all overlapping regions. The chloroplast genome of the olive has an organisation and gene order that is conserved among numerous Angiosperm species and do not contain any of the inversions, gene duplications, insertions, inverted repeat expansions and gene/intron losses that have been found in the chloroplast genomes of the genera Jasminum and Menodora, from the same family as Olea.The annotated sequence was used to evaluate the content of coding genes, the extent, and distribution of repeated and long dispersed sequences and the nucleotide composition pattern. These analyses provided essential information for structural, functional and comparative genomic studies in olive plastids. Furthermore, the alignment of the olive plastome sequence to those of other varieties and species identified 30 new organellar polymorphisms within the cultivated olive.

CONCLUSIONS

In addition to identifying mutations that may play a functional role in modifying the metabolism and adaptation of olive cultivars, the new chloroplast markers represent a valuable tool to assess the level of olive intercultivar plastome variation for use in population genetic analysis, phylogenesis, cultivar characterisation and DNA food tracking.