An updated phylogenetic and biogeographic analysis based on genome skimming data reveals convergent evolution of shrubby habit in Clematis in the Pliocene and Pleistocene

Convergent evolution, often viewed as the inevitable outcome of natural selection, has received special attention since the time of Darwin. Clematis is well known for its climbing habit, but it has some shrubby species, known as sect. Fruticella s.l. The shrubby Clematis species are distributed in the dry habitats of Central Asia and adjacent areas showing possible convergent evolution. In this study, we assembled the complete plastome and nuclear ribosomal DNA (nrDNA) sequences of 56 Clematis species, representing most sections and covering most of the shrubby species, to reconstruct their evolutionary histories. Using both maximum likelihood and Bayesian methods, the plastome and nrDNA datasets generated similar, but not identical, phylogenetic relationships, which are better resolved than in previous studies. Then, molecular dating, historical range reconstruction, and character optimization analyses were conducted based on this updated phylogenetic framework. All the morphological characters widely used for taxonomy were shown to have evolved multiple times. Molecular dating inferred that Clematis diverged from its sister in the mid Miocene, and all six major clades of Clematis originated during the late Miocene, with a species radiation during the Pliocene to Pleistocene. The results clearly showed that the shrubby habit evolved independently in four lineages of Clematis in Asia. We also revealed that the shrubby lineages have emerged since the very beginning of Pliocene. Asian monsoon variation in the Pliocene and glacial period fluctuation in the Pleistocene may be the driving forces for the origin and diversification of the shrubby Clematis in Central Asia and adjacent dry areas.

The complete chloroplast genome of Artocarpus altilis (Moraceae) and phylogenetic relationships

The chloroplast (cp) is an essential organelle in higher plants. The genes of the plastome are well suited to infer phylogenetic relationships among taxa. In this study, we report the assembly of the cp genome of Artocarpus altilis and its phylogeny among species from Moraceae family. The cp genome of A. altilis was 160,822 base pair (bp) in length, comprising one large single-copy region of 88,692 bp, one small single-copy region of 19,290 bp, and a pair of inverted repeat regions (IRs) of 26,420 bp. A total of 113 different genes were predicted, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic analysis of 19 species belonging to the Moraceae family confirmed the phylogenetic proximity of the genus Artocarpus and Morus and the genetic similarity of A. camansi and A. altilis.

Discordant evolution of organellar genomes in peas (Pisum L.)

Plastids and mitochondria have their own small genomes, which do not undergo meiotic recombination and may have evolutionary fates different from each other and that of the nuclear genome. For the first time, we sequenced mitochondrial genomes of pea (Pisum L.) from 42 accessions mostly representing diverse wild germplasm from throughout the wild pea geographical range. Six structural types of the pea mitochondrial genome were revealed. From the same accessions, plastid genomes were sequenced. Phylogenetic trees based on the plastid and mitochondrial genomes were compared. The topologies of these trees were highly discordant, implying not less than six events of hybridisation between diverged wild peas in the past, with plastids and mitochondria differently inherited by the descendants. Such discordant inheritance of organelles could have been driven by plastid-nuclear incompatibility, which is known to be widespread in crosses involving wild peas and affects organellar inheritance. The topology of the phylogenetic tree based on nucleotide sequences of a nuclear gene, His5, encoding a histone H1 subtype, corresponded to the current taxonomy and resembled that based on the plastid genome. Wild peas (Pisum sativum subsp. elatius s.l.) inhabiting Southern Europe were shown to be of hybrid origin, resulting from crosses of peas related to those presently inhabiting the eastern Mediterranean in a broad sense. These results highlight the roles of hybridisation and cytonuclear conflict in shaping plant microevolution.

Complete plastid genome of Cumathamnion serrulatum (Ceramiales, Rhodophyta)

We report the complete plastid genome of Cumathamnion serrulatum, also known as Delesseria serrulata. The plastid genome was 174,192 bp in size. Annotation showed there were 193 protein coding genes, three ribosomal RNAs, and 29 transfer RNAs. One intron was found, and the GC content was 27.2%. The maximum likelihood tree with the concatenated 177 plastid coding genes showed a strong monophyletic relationship to Membranoptera spp. within the Ceramiales.

Elatostema qinzhouense (Urticaceae), a new species from limestone karst in Guangxi, China

Elatostema qinzhouense L.F. Fu, A.K. Monro & Y.G. Wei, a new species from Guangxi, China is described and illustrated. Morphologically, E. qinzhouense is most similar to E. hezhouense from which it differs by having smaller size of leaf laminae, fewer and smaller staminate peduncle bracts, longer pistillate peduncle bracts and a larger achene. This result is supported by the molecular evidence. The phylogenetic position of the new species within Elatostema is evaluated using three DNA regions, ITS, trnH-psbA and psbM-trnD, for 107 taxa of Elatostema s.l. (including E. qinzhouense). Bayesian inference (BI) and maximum likelihood (ML) analyses each recovered the same strongly supported tree topologies, indicating that E. qinzhouense is a member of the core Elatostema clade and sister to E. hezhouense. Along with the phylogenetic studies, plastid genome and ribosomal DNA (rDNA) sequences of the new species are assembled and annotated. The plastid genome is 150,398 bp in length and comprises two inverted repeats (IRs) of 24,688 bp separated by a large single-copy of 83,919 bp and a small single-copy of 17,103 bp. A total of 113 functional genes are recovered, comprising 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The rDNA is 5,804 bp in length and comprised the 18S ribosomal RNA partial sequence (1,809 bp), internal transcribed spacer 1 (213 bp), 5.8S ribosomal RNA (164 bp), internal transcribed spacer 2 (248 bp) and 26S ribosomal RNA partial sequence (3,370 bp). In addition, the chromosome number of E. qinzhouense is observed to be 2n = 26, suggesting that the species is diploid. Given a consistent relationship between ploidy level and reproductive system in Elatostema, the new species is also considered to be sexually reproducing. Our assessment of the extinction threat for E. qinzhouense is that it is Endangered (EN) according to the criteria of the International Union for Conservation of Nature.

Comparative analysis of complete plastid genomes from Lilium lankongense Franchet and its closely related species and screening of Lilium-specific primers

Lilium lankongense Franchet is a lily species found on the Qinghai-Tibet Plateau. It is pink with deep red spots, has a high ornamental value, and is used in hybrid breeding of horticultural lily varieties. We have insufficient knowledge of the genetic resources of L. lankongense and its phylogenetic relationships with related species. Recent molecular phylogenetic studies have shown a very close phylogenetic relationship between L. lankongense and the five species L. duchartrei, L. stewartianum, L. matangense, L. lophophorum, and L. nanum. However, molecular markers still lack sufficient signals for population-level research of the genus Lilium. We sequenced and compared the complete plastid sequences of L. lankongense and its five related species. The genomes ranged from 152,307 bp to 152,611 bp. There was a slight inconsistency detected in inverted repeat and single copy boundaries and there were 53 to 63 simple sequence repeats in the six species. Two of the 12 highly variable regions (trnC-petN and rpl32-trnL) were verified in 11 individuals and are promising for population-level studies. We used the complete sequence of 33 plastid genomes, the protein-coding region sequence, and the nuclear ITS sequence to reconstruct the phylogenetic tree of Lilium species. Our results showed that the plastid gene tree and nuclear gene tree were not completely congruent, which may be caused by hybridization, insufficient information contained in the nuclear ITS, or the small number of samples. The results of phylogenetic analysis based on plastid genomes indicated that the six Lilium species were closely related. Our study provides a preliminarily rebuilt backbone phylogeny that is significant for future molecular and morphological studies of Lilium.

Complete chloroplast genome sequence of Passiflora serrulata Jacq. (Passifloraceae)

This study was the first report for the complete chloroplast genome of Passiflora serrulata Jacq. (Passifloraceae). The cp genome was 149,683 bp in length contained two inverted repeats (IRs) of 25,470 bp, which were separated by large single-copy (LSC) and small single-copy (SSC) of 86,252 bp and 13,491 bp, respectively. A total of 110 functional genes were encoded, comprised 76 protein-coding genes, 30 tRNA genes, and four rRNA genes. The GC content was 37.0%. The maximum likelihood phylogenetic tree indicated that P. serrulata was recovered as the member of subg. Passiflora and most closely related to the clade formed by P. serratodigitata and P. ligularis.

Comparative plastid genomics of four Pilea (Urticaceae) species: insight into interspecific plastid genome diversity in Pilea

BACKGROUND

Pilea is a genus of perennial herbs from the family Urticaceae, and some species are used as courtyard ornamentals or for medicinal purposes. At present, there is no information about the plastid genome of Pilea, which limits our understanding of this genus. Here, we report 4 plastid genomes of Pilea taxa (Pilea mollis, Pilea glauca 'Greizy', Pilea peperomioides and Pilea serpyllacea 'Globosa') and performed comprehensive comparative analysis.

RESULTS

The four plastid genomes all have a typical quartile structure. The lengths of the plastid genomes ranged from 150,398 bp to 152,327 bp, and each genome contained 113 unique genes, including 79 protein-coding genes, 4 rRNA genes, and 30 tRNA genes. Comparative analysis showed a rather high level of sequence divergence in the four genomes. Moreover, eight hypervariable regions were identified (petN-psbM, psbZ-trnG-GCC, trnT-UGU-trnL-UAA, accD-psbI, ndhF-rpl32, rpl32-trnL-UAG, ndhA-intron and ycf1), which are proposed for use as DNA barcode regions. Phylogenetic relationships based on the plastid genomes of 23 species of 14 genera of Urticaceae resulted in the placement of Pilea in the middle and lower part of the phylogenetic tree, with 100% bootstrap support within Urticaceae.

CONCLUSION

Our results enrich the resources concerning plastid genomes. Comparative plastome analysis provides insight into the interspecific diversity of the plastid genome of Pilea. The identified hypervariable regions could be used for developing molecular markers applicable in various research areas.

Changes in plastid biogenesis leading to the formation of albino regenerants in barley microspore culture

BACKGROUND

Microspore embryogenesis is potentially the most effective method of obtaining doubled haploids (DH) which are utilized in breeding programs to accelerate production of new cultivars. However, the regeneration of albino plants significantly limits the exploitation of androgenesis for DH production in cereals. Despite many efforts, the precise mechanisms leading to development of albino regenerants have not yet been elucidated. The objective of this study was to reveal the genotype-dependent molecular differences in chloroplast differentiation that lead to the formation of green and albino regenerants in microspore culture of barley.

RESULTS

We performed a detailed analysis of plastid differentiation at successive stages of androgenesis in two barley cultivars, 'Jersey' and 'Mercada' that differed in their ability to produce green regenerants. We demonstrated the lack of transition from the NEP-dependent to PEP-dependent transcription in plastids of cv. 'Mercada' that produced mostly albino regenerants in microspore culture. The failed NEP-to-PEP transition was associated with the lack of activity of Sig2 gene encoding a sigma factor necessary for transcription of plastid rRNA genes. A very low level of 16S and 23S rRNA transcripts and impaired plastid translation machinery resulted in the inhibition of photomorphogenesis in regenerating embryos and albino regenerants. Furthermore, the plastids present in differentiating 'Mercada' embryos contained a low number of plastome copies whose replication was not always completed. Contrary to 'Mercada', cv. 'Jersey' that produced 90% green regenerants, showed the high activity of PEP polymerase, the highly increased expression of Sig2, plastid rRNAs and tRNAGlu, which indicated the NEP inhibition. The increased expression of GLKs genes encoding transcription factors required for induction of photomorphogenesis was also observed in 'Jersey' regenerants.

CONCLUSIONS

Proplastids present in microspore-derived embryos of albino-producing genotypes did not pass the early checkpoints of their development that are required for induction of further light-dependent differentiation of chloroplasts. The failed activation of plastid-encoded RNA polymerase during differentiation of embryos was associated with the genotype-dependent inability to regenerate green plants in barley microspore culture. The better understanding of molecular mechanisms underlying formation of albino regenerants may be helpful in overcoming the problem of albinism in cereal androgenesis.

DNA barcoding and molecular phylogeny of Dumasia (Fabaceae: Phaseoleae) reveals a cryptic lineage

Dumasia taxonomy and classification have long been problematic. Species within this genus have few morphological differences and plants without flowers or fruits are difficult to accurately identify. In this study, we evaluated the ability of six DNA barcoding sequences, one nuclear (ITS) and five chloroplast regions (trnH-psbA, matK, rbcL, trnL-trnF, psbB-psbF), to efficiently identify Dumasia species. Most single markers or their combinations identify obvious barcoding gaps between intraspecific and interspecific genetic variation. Most combined analyses including ITS showed good species resolution and identification efficiency. We therefore suggest that ITS alone or a combination of ITS with any cpDNA marker are most suitable for DNA barcoding of Dumasia. The phylogenetic analyses clearly indicated that Dumasia yunnanensis is not monophyletic and is separated as two independent branches, which may result from cryptic differentiation. Our results demonstrate that molecular data can deepen the comprehension of taxonomy of Dumasia and provide an efficient approach for identification of the species.

The complete plastid genome of Ficus erecta (Moraceae)

Ficus erecta is a wild relative of F. carica, which is an important economic plant. Here, we determined the complete plastid genome of F. erecta using the Illumina paired reads to provide genomic feature resources. The whole plastid genome of F. erecta is 160,603 bp in length, containing two inverted repeats (IRs) of 25,899 bp separated by a large single-copy (LSC 88,640 bp) region and a small single-copy (SSC 20,165 bp) region. The complete plastome sequence of Ficus erecta will provide a useful resource for the evolutionary biology study as well as for the phylogenetic studies.

The complete chloroplast genome of Akebia trifoliata subsp. australis (Lardizabalaceae), a medicinal plant in China

Akebia trifoliata subsp. australis (Diels) T.Shimizu is a medicinal plant in China. Here, the complete chloroplast (cp) genome sequence of A. trifoliata subsp. australis was assembled and characterized as a resource for future genetic studies. The whole cp genome was 157,952 bp in length, containing of a large single-copy (LSC) region of 86,596 bp, a small single-copy (SSC) region of 19,060 bp, and two inverted repeat (IR) regions of 26,148 bp. The new sequence possessed total 113 unique genes, including 79 protein-coding genes, 30 tRNA genes and 4 rRNA genes. The nucleotide composition was asymmetric (30.3% A, 19.7% C, 19.0% G and 31.0% T) with an overall GC content of 38.7%. The maximum likelihood phylogenetic analysis based on 10 cp genomes indicated that A. trifoliata subsp. australis was closely related to Akebia trifoliata subsp. trifoliata. However, Akebia quinata was closely related to Stauntonia obovatifoliola.

Phylogeny of Euclidieae (Brassicaceae) based on plastome and nuclear ribosomal DNA data

Euclidieae, a morphologically diverse tribe in the family Brassicaceae (Cruciferae), consists of 29 genera and more than 150 species distributed mainly in Asia. Prior phylogenetic analyses on Euclidieae are inadequate. In this study, sequence data from the plastid genome and nuclear ribosomal DNA of 72 species in 27 genera of Euclidieae were used to infer the inter- and intra-generic relationships within. The well-resolved and strongly supported plastome phylogenies revealed that Euclidieae could be divided into five clades. Both Cymatocarpus and Neotorularia are polyphyletic in nuclear and plastome phylogenies. Besides, the conflicts of systematic positions of three species of Braya and two species of Solms-laubachia s.l. indicated that hybridization and or introgression might have happened during the evolutionary history of the tribe. Results from divergence-time analyses suggested an early Miocene origin of Euclidieae, and it probably originated from the Central Asia, Pamir Plateau and West Himalaya. In addition, multiple ndh genes loss and pseudogenization were detected in eight species based on comparative genomic study.

Comparison of eight complete plastid genomes from three moss families Amblystegiaceae, Calliergonaceae and Pylaisiaceae

We sequenced and assembled eight complete plastid genomes from three closely related pleurocarpous moss families: Amblystegium serpens, Campyliadelphus stellatus, Cratoneuron filicinum, Drepanocladus aduncus, and Leptodictyum humile (Amblystegiaceae), Calliergon sarmentosum and Warnstorfia exannulata (Calliergonaceae), and Calliergonella cuspidata (Pylaisiaceae). The newly generated plastid genomes range from 124,256 to 124,819 bp, with two inverted repeat regions (9,624–9,696 bp) separated by a large single-copy region (86,422–86,924 bp) and a small single-copy region (18,430–18,514 bp). All these plastid genomes encode 116 unique genes including 82 protein-coding genes, 30 tRNA genes and four rRNAgenes. The overall GC content is between 28.6%–29.3%. Phylogenetic analysis showed that all Amblystegiaceae species Amblystegium serpens, Campyliadelphus stellatus, Cratoneuron filicinum, Drepanocladus aduncus, Leptodictyum humile, and Sanionia uncinata clustered in one clade, which is sister to the Pylaisiaceae species Calliergonella cuspidata. The two Calliergonaceae species Calliergon sarmentosum and Warnstorfia exannulata form a clade and is sister to Amblystegiaceae and Pylaisiaceae.

Degradation of key photosynthetic genes in the critically endangered semi-aquatic flowering plant Saniculiphyllum guangxiense (Saxifragaceae)

BACKGROUND

Plastid gene loss and pseudogenization has been widely documented in parasitic and mycoheterotrophic plants, which have relaxed selective constraints on photosynthetic function. More enigmatic are sporadic reports of pseudogenization and loss of important photosynthesis genes in lineages thought to be fully photosynthetic. Here we report the complete plastid genome of Saniculiphyllum guangxiense, a critically endangered and phylogenetically isolated plant lineage, along with genomic evidence of reduced chloroplast function. We also report 22 additional plastid genomes representing the diversity of its containing clade Saxifragales, characterizing gene content and placing variation in a broader phylogenetic context.

RESULTS

We find that the plastid genome of Saniculiphyllum has experienced pseudogenization of five genes of the ndh complex (ndhA, ndhB, ndhD, ndhF, and ndhK), previously reported in flowering plants with an aquatic habit, as well as the surprising pseudogenization of two genes more central to photosynthesis (ccsA and cemA), contrasting with strong phylogenetic conservatism of plastid gene content in all other sampled Saxifragales. These genes participate in photooxidative protection, cytochrome synthesis, and carbon uptake. Nuclear paralogs exist for all seven plastid pseudogenes, yet these are also unlikely to be functional.

CONCLUSIONS

Saniculiphyllum appears to represent the greatest degree of plastid gene loss observed to date in any fully photosynthetic lineage, perhaps related to its extreme habitat specialization, yet plastid genome length, structure, and substitution rate are within the variation previously reported for photosynthetic plants. These results highlight the increasingly appreciated dynamism of plastid genomes, otherwise highly conserved across a billion years of green plant evolution, in plants with highly specialized life history traits.

Plastome of the mycoheterotrophic eudicot Exacum paucisquama (Gentianaceae) exhibits extensive gene loss and a highly expanded inverted repeat region

Mycoheterotrophic plants are highly specialized species able to acquire organic carbon from symbiotic fungi, with relaxed dependence on photosynthesis for carbon fixation. The relaxation of the functional constraint of photosynthesis and thereby the relaxed selective pressure on functional photosynthetic genes usually lead to substantial gene loss and a highly degraded plastid genome in heterotrophs. In this study, we sequenced and analyzed the plastome of the eudicot Exacum paucisquama, providing the first plastid genome of a mycoheterotroph in the family Gentianaceae to date. The E. paucisquama plastome was 44,028 bp in length, which is much smaller than the plastomes of autotrophic eudicots. Although the E. paucisquama plastome had a quadripartite structure, a distinct boundary shift was observed in comparison with the plastomes of other eudicots. We detected extensive gene loss and only 21 putative functional genes (15 protein-coding genes, four rRNA genes and two tRNA genes). Our results provide valuable information for comparative evolutionary analyses of plastomes of heterotrophic species belonging to different phylogenetic groups.

Complete plastome sequencing resolves taxonomic relationships among species of Calligonum L. (Polygonaceae) in China

BACKGROUND

Calligonum (Polygonaceae) is distributed from southern Europe through northern Africa to central Asia, and is typically found in arid, desert regions. Previous studies have revealed that standard DNA barcodes fail to discriminate Calligonum species. In this study, the complete plastid genomes (plastome) for 32 accessions of 21 Calligonum species is sequenced to not only generate the first complete plastome sequence for the genus Calligonum but to also 1) Assess the ability of the complete plastome sequence to discern species within the group, and 2) screen the plastome sequence for a cost-effective DNA barcode that can be used in future studies to resolve taxonomic relationships within the group.

RESULTS

The whole plastomes of Calligonum species possess a typical quadripartite structure. The size of the Calligonum plastome is approximately 161 kilobase pairs (kbp), and encodes 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Based on ML phylogenetic tree analyses, the complete plastome has higher species identification (78%) than combinations of standard DNA barcodes (rbcL + matK + nrITS, 56%). Five newly screened gene regions (ndhF, trnS-G, trnC-petN, ndhF-rpl32, rpl32-trnL) had high species resolution, where ndhF and trnS-G were able to distinguish the highest proportion of Calligonum species (56%).

CONCLUSIONS

The entire plastid genome was the most effective barcode for the genus Calligonum, although other gene regions showed great potential as taxon-specific barcodes for species identification in Calligonum.

The loss of photosynthesis pathway and genomic locations of the lost plastid genes in a holoparasitic plant Aeginetia indica

BACKGROUND

With three origins of holoparasitism, Orobanchaceae provides an ideal system to study the evolution of holoparasitic lifestyle in plants. The evolution of holoparasitism can be revealed by plastid genome degradation and coordinated changes in the nuclear genome, since holoparasitic plants lost the capability of photosynthesis. Among the three clades with holoparasitic plants in Orobanchaceae, only Clade VI has no available plastid genome sequences for holoparasitic plants. In this study, we sequenced the plastome and transcriptome of Aeginetia indica, a holoparasitic plant in Clade VI of Orobanchaceae, to study its plastome evolution and the corresponding changes in the nuclear genome as a response of the loss of photosynthetic function.

RESULTS

The plastome of A. indica is reduced to 86,212 bp in size, and almost all photosynthesis-related genes were lost. Massive fragments of the lost plastid genes were transferred into the mitochondrial and/or nuclear genomes. These fragments could not be detected in its transcriptomes, suggesting that they were non-functional. Most protein coding genes in the plastome showed the signal of relaxation of purifying selection. Plastome and transcriptome analyses indicated that the photosynthesis pathway is completely lost, and that the porphyrin and chlorophyll metabolism pathway is partially retained, although chlorophyll synthesis is not possible.

CONCLUSIONS

Our study suggests the loss of photosynthesis-related functions in A. indica in both the nuclear and plastid genomes. The lost plastid genes are transferred into its nuclear and/or mitochondrial genomes, and exist in very small fragments with no expression and are thus non-functional. The Aeginetia indica plastome also provides a resource for comparative studies on the repeated evolution of holoparasitism in Orobanchaceae.

Implications of plastome evolution in the true lilies (monocot order Liliales)

The families of the monocot order Liliales exhibit highly contrasting characteristic of photosynthetic and mycoheterotrophic life histories. Although previous phylogenetic and morphological studies of Liliales have been conducted, they have not examined molecular evolution associated with this contrasting phenomenon. Here, we conduct the first comparative plastome study of all ten families of Liliales using 29 newly sequenced plastid genomes analyzed together with previously published data. We also present a phylogenetic analysis for Liliales of 78 plastid genes combined with 22 genes from all three genomes (nuclear 18S rDNA and phyC; 17 plastid genes; and mitochondrial matR, atpA, and cob). Within the newly generated phylogenetic tree of Liliales, we evaluate the ancestral state changes of selected morphological traits in the order. There are no significant differences in plastid genome features among species that show divergent characteristics correlated with family circumscriptions. However, the results clearly differentiate between photosynthetic and mycoheterotrophic taxa of Liliales in terms of genome structure, and gene content and order. The newly sequenced plastid genomes and combined three-genome data revealed Smilacaceae as sister to Liliaceae instead of Philesiaceae and Ripogonaceae. Additionally, we propose a revised familial classification system of Liliales that consists of nine families, considering Ripogonaceae a synonym of Philesiaceae. The ancestral state reconstruction indicated synapomorphies for each family of Liliales, except Liliaceae, Melanthiaceae and Colchicaceae. A taxonomic key for all nine families of Liliales is also provided.

-Complete plastid genome sequences of two species of the Neotropical genus Brunellia (Brunelliaceae)

Here we present the first two complete plastid genomes for Brunelliaceae, a Neotropical family with a single genus, Brunellia. We surveyed the entire plastid genome in order to find variable cpDNA regions for further phylogenetic analyses across the family. We sampled morphologically different species, B. antioquensis and B. trianae, and found that the plastid genomes are 157,685 and 157,775 bp in length and display the typical quadripartite structure found in angiosperms. Despite the clear morphological distinction between both species, the molecular data show a very low level of divergence. The amount of nucleotide substitutions per site is one of the lowest reported to date among published congeneric studies (π = 0.00025). The plastid genomes have gene order and content coincident with other COM (Celastrales, Oxalidales, Malpighiales) relatives. Phylogenetic analyses of selected superrosid representatives show high bootstrap support for the ((C,M)O) topology. The N-fixing clade appears as the sister group of the COM clade and Zygophyllales as the sister to the rest of the fabids group.

Complete plastid genome of the Chinese medicinal herb Paeonia obovata subsp. Willmottiae (Paeoniaceae): characterization and phylogeny

The plastid genome (plastome) of the endemic Chinese medicinal herb Paeonia obovata subsp. Willmottiae (Paeoniaceae) was sequenced and investigated in this study. The complete plastome is 152,713 bp in length with the typical quadripartite structure, which consists of a large single-copy region (LSC, 84,419 bp), a small single-copy region (SSC, 16,982 bp), and a pair of inverted repeat regions (IRs, 25,656 bp). The overall GC content is 33.2%, and the IR regions are more GC rich (43.2%) than the LSC (36.7%) and SSC (32.8%) regions. A total of 114 unique genes, including 79 protein-coding genes, 31 tRNAs, and four rRNAs were identified. Phylogenetic reconstruction based on complete plastome sequences demonstrated that P. obovata subsp. Willmottiae is phylogenetically closest to P. obovata.

Large scale genome skimming from herbarium material for accurate plant identification and phylogenomics

BACKGROUND

Herbaria are valuable sources of extensive curated plant material that are now accessible to genetic studies because of advances in high-throughput, next-generation sequencing methods. As an applied assessment of large-scale recovery of plastid and ribosomal genome sequences from herbarium material for plant identification and phylogenomics, we sequenced 672 samples covering 21 families, 142 genera and 530 named and proposed named species. We explored the impact of parameters such as sample age, DNA concentration and quality, read depth and fragment length on plastid assembly error. We also tested the efficacy of DNA sequence information for identifying plant samples using 45 specimens recently collected in the Pilbara.

RESULTS

Genome skimming was effective at producing genomic information at large scale. Substantial sequence information on the chloroplast genome was obtained from 96.1% of samples, and complete or near-complete sequences of the nuclear ribosomal RNA gene repeat were obtained from 93.3% of samples. We were able to extract sequences for the core DNA barcode regions rbcL and matK from 96 to 93.3% of samples, respectively. Read quality and DNA fragment length had significant effects on sequencing outcomes and error correction of reads proved essential. Assembly problems were specific to certain taxa with low GC and high repeat content (Goodenia, Scaevola, Cyperus, Bulbostylis, Fimbristylis) suggesting biological rather than technical explanations. The structure of related genomes was needed to guide the assembly of repeats that exceeded the read length. DNA-based matching proved highly effective and showed that the efficacy for species identification declined in the order cpDNA >> rDNA > matK >> rbcL.

CONCLUSIONS

We showed that a large-scale approach to genome sequencing using herbarium specimens produces high-quality complete cpDNA and rDNA sequences as a source of data for DNA barcoding and phylogenomics.

Intra-individual heteroplasmy in the Gentiana tongolensis plastid genome (Gentianaceae)

Chloroplasts are typically inherited from the female parent and are haploid in most angiosperms, but rare intra-individual heteroplasmy in plastid genomes has been reported in plants. Here, we report an example of plastome heteroplasmy and its characteristics in Gentiana tongolensis (Gentianaceae). The plastid genome of G. tongolensis is 145,757 bp in size and is missing parts of petD gene when compared with other Gentiana species. A total of 112 single nucleotide polymorphisms (SNPs) and 31 indels with frequencies of more than 2% were detected in the plastid genome, and most were located in protein coding regions. Most sites with SNP frequencies of more than 10% were located in six genes in the LSC region. After verification via cloning and Sanger sequencing at three loci, heteroplasmy was identified in different individuals. The cause of heteroplasmy at the nucleotide level in plastome of G. tongolensis is unclear from the present data, although biparental plastid inheritance and transfer of plastid DNA seem to be most likely. This study implies that botanists should reconsider the heredity and evolution of chloroplasts and be cautious with using chloroplasts as genetic markers, especially in Gentiana.

The complete plastid genome of Iris domestica: a traditional Chinese medicine

Iris domestica has been used as Chinese traditional medicine to treat inflammation and throat disorders for many centuries. In this study, the complete plastid genome of I. domestica was first reported and characterized. The complete plastid genome is a typical quadripartite circular molecule of 153,729 bp in length, including a large single copy (LSC) region of 83,136 bp and a small single copy (SSC) region of 18,165 bp separated by two inverted repeat (IR) regions of 26,214 bp. A total of 132 genes including 86 protein-coding genes, 38 tRNA genes, and eight rRNA genes were identified. The phylogenetic analysis revealed that I. domestica was closer to I. gatesii.

Complete chloroplast genome of the tiny marine diatom Nanofrustulum shiloi (Bacillariophyta) from the Adriatic Sea

We report the chloroplast genome sequence of Nanofrustulum shiloi, a tiny araphid pennate diatom collected from the Adriatic Sea. The 160,994-bp N. shiloi genome displays a quadripartite structure and its gene repertoire resembles those of other diatom chloroplast genomes. Besides the genes located in the inverted repeat, psbY is duplicated. A gene-poor region in the large single-copy region contains multiple ORFs sharing sequence similarities with plasmids and chloroplast ORFs found in other diatom species. The genome features a single intron, a group II intron in petB. Phylogenomic analysis identified N. shiloi at a basal position within the araphid 2 clade.