The plastid genome of twenty-two species from Ferula, Talassia, and Soranthus: comparative analysis, phylogenetic implications, and adaptive evolution

BACKGROUND

The Ferula genus encompasses 180-185 species and is one of the largest genera in Apiaceae, with many of Ferula species possessing important medical value. The previous studies provided more information for Ferula, but its infrageneric relationships are still confusing. In addition, its genetic basis of its adaptive evolution remains poorly understood. Plastid genomes with more variable sites have the potential to reconstruct robust phylogeny in plants and investigate the adaptive evolution of plants. Although chloroplast genomes have been reported within the Ferula genus, few studies have been conducted using chloroplast genomes, especially for endemic species in China.

RESULTS

Comprehensively comparative analyses of 22 newly sequenced and assembled plastomes indicated that these plastomes had highly conserved genome structure, gene number, codon usage, and repeats type and distribution, but varied in plastomes size, GC content, and the SC/IR boundaries. Thirteen mutation hotspot regions were detected and they would serve as the promising DNA barcodes candidates for species identification in Ferula and related genera. Phylogenomic analyses with high supports and resolutions showed that Talassia transiliensis and Soranthus meyeri were nested in the Ferula genus, and thus they should be transferred into the Ferula genus. Our phylogenies also indicated the monophyly of subgenera Sinoferula and subgenera Narthex in Ferula genus. Twelve genes with significant posterior probabilities for codon sites were identified in the positively selective analysis, and their function may relate to the photosystem II, ATP subunit, and NADH dehydrogenase. Most of them might play an important role to help Ferula species adapt to high-temperatures, strong-light, and drought habitats.

CONCLUSION

Plastome data is powerful and efficient to improve the support and resolution of the complicated Ferula phylogeny. Twelve genes with significant posterior probabilities for codon sites were helpful for Ferula to adapt to the harsh environment. Overall, our study supplies a new perspective for comprehending the phylogeny and evolution of Ferula.

Putting small and big pieces together: a genome assembly approach reveals the largest Lamiid plastome in a woody vine

The plastid genome of flowering plants generally shows conserved structural organization, gene arrangement, and gene content. While structural reorganizations are uncommon, examples have been documented in the literature during the past years. Here we assembled the entire plastome of Bignonia magnifica and compared its structure and gene content with nine other Lamiid plastomes. The plastome of B. magnifica is composed of 183,052 bp and follows the canonical quadripartite structure, synteny, and gene composition of other angiosperms. Exceptionally large inverted repeat (IR) regions are responsible for the uncommon length of the genome. At least four events of IR expansion were observed among the seven Bignoniaceae species compared, suggesting multiple expansions of the IRs over the SC regions in the family. A comparison with 6,231 other complete plastomes of flowering plants available on GenBank revealed that the plastome of B. magnifica is the longest Lamiid plastome described to date. The newly generated plastid genome was used as a source of selected genes. These genes were combined with orthologous regions sampled from other species of Bignoniaceae and all gene alignments concatenated to infer a phylogeny of the family. The tree recovered is consistent with known relationships within the Bignoniaceae.

Towards the plastome evolution and phylogeny of Cycas L. (Cycadaceae): molecular-morphology discordance and gene tree space analysis

BACKGROUND

Plastid genomes (plastomes) present great potential in resolving multiscale phylogenetic relationship but few studies have focused on the influence of genetic characteristics of plastid genes, such as genetic variation and phylogenetic discordance, in resolving the phylogeny within a lineage. Here we examine plastome characteristics of Cycas L., the most diverse genus among extant cycads, and investigate the deep phylogenetic relationships within Cycas by sampling 47 plastomes representing all major clades from six sections.

RESULTS

All Cycas plastomes shared consistent gene content and structure with only one gene loss detected in Philippine species C. wadei. Three novel plastome regions (psbA-matK, trnN-ndhF, chlL-trnN) were identified as containing the highest nucleotide variability. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection except ndhB. Phylogenomic analyses that alternatively included concatenated and coalescent methods, both identified four clades but with conflicting topologies at shallow nodes. Specifically, we found three species-rich Cycas sections, namely Stangerioides, Indosinenses and Cycas, were not or only weakly supported as monophyly based on plastomic phylogeny. Tree space analyses based on different tree-inference methods both revealed three gene clusters, of which the cluster with moderate genetic properties showed the best congruence with the favored phylogeny.

CONCLUSIONS

Our exploration in plastomic data for Cycas supports the idea that plastid protein-coding genes may exhibit discordance in phylogenetic signals. The incongruence between molecular phylogeny and morphological classification reported here may largely be attributed to the uniparental attribute of plastid, which cannot offer sufficient information to resolve the phylogeny. Contrasting to a previous consensus that genes with longer sequences and a higher proportion of variances are superior for phylogeny reconstruction, our result implies that the most effective phylogenetic signals could come from loci that own moderate variation, GC content, sequence length, and underwent modest selection.

The complete plastomes of seven Peucedanum plants: comparative and phylogenetic analyses for the Peucedanum genus

BACKGROUND

The Peucedanum genus is the backbone member of Apiaceae, with many economically and medically important plants. Although the previous studies on Peucedanum provide us with a good research basis, there are still unclear phylogenetic relationships and many taxonomic problems in Peucedanum, and a robust phylogenetic framework of this genus still has not been obtained, which severely hampers the improvement and revision of taxonomic system for this genus. The plastid genomes possessing more variable characters have potential for reconstructing a robust phylogeny in plants.

RESULTS

In the current study, we newly sequenced and assembled seven Peucedanum plastid genomes. Together with five previously published plastid genomes of Peucedanum, we performed a comprehensively comparative analyses for this genus. Twelve Peucedanum plastomes were similar in terms of genome structure, codon bias, RNA editing sites, and SSRs, but varied in genome size, gene content and arrangement, and border of SC/IR. Fifteen mutation hotspot regions were identified among plastid genomes that can serve as candidate DNA barcodes for species identification in Peucedanum. Our phylogenetic analyses based on plastid genomes generated a phylogeny with high supports and resolutions for Peucedanum that robustly supported the non-monophyly of genus Peucedanum.

CONCLUSION

The plastid genomes of Peucedanum showed both conservation and diversity. The plastid genome data were efficient and powerful for improving the supports and resolutions of phylogeny for the complex Peucedanum genus. In summary, our study provides new sights into the plastid genome evolution, taxonomy, and phylogeny for Peucedanum species.

Extensive genomic rearrangements mediated by repetitive sequences in plastomes of Medicago and its relatives

BACKGROUND

Although plastomes are highly conserved with respect to gene content and order in most photosynthetic angiosperms, extensive genomic rearrangements have been reported in Fabaceae, particularly within the inverted repeat lacking clade (IRLC) of Papilionoideae. Two hypotheses, i.e., the absence of the IR and the increased repeat content, have been proposed to affect the stability of plastomes. However, this is still unclear for the IRLC species. Here, we aimed to investigate the relationships between repeat content and the degree of genomic rearrangements in plastomes of Medicago and its relatives Trigonella and Melilotus, which are nested firmly within the IRLC.

RESULTS

We detected abundant repetitive elements and extensive genomic rearrangements in the 75 newly assembled plastomes of 20 species, including gene loss, intron loss and gain, pseudogenization, tRNA duplication, inversion, and a second independent IR gain (IR ~ 15 kb in Melilotus dentata) in addition to the previous first reported cases in Medicago minima. We also conducted comparative genomic analysis to evaluate plastome evolution. Our results indicated that the overall repeat content is positively correlated with the degree of genomic rearrangements. Some of the genomic rearrangements were found to be directly linked with repetitive sequences. Tandem repeated sequences have been detected in the three genes with accelerated substitution rates (i.e., accD, clpP, and ycf1) and their length variation could be explained by the insertions of tandem repeats. The repeat contents of the three localized hypermutation regions around these three genes with accelerated substitution rates are also significantly higher than that of the remaining plastome sequences.

CONCLUSIONS

Our results suggest that IR reemergence in the IRLC species does not ensure their plastome stability. Instead, repeat-mediated illegitimate recombination is the major mechanism leading to genome instability, a pattern in agreement with recent findings in other angiosperm lineages. The plastome data generated herein provide valuable genomic resources for further investigating the plastome evolution in legumes.

Phylogenomic and evolutionary dynamics of inverted repeats across Angelica plastomes

BACKGROUND

Angelica L. (family Apiaceae) is an economically important genus comprising ca. One hundred ten species. Angelica species are found on all continents of the Northern Hemisphere, and East Asia hosts the highest number of species. Morphological characters such as fruit anatomy, leaf morphology and subterranean structures of Angelica species show extreme diversity. Consequently, the taxonomic classification of Angelica species is complex and remains controversial, as the classifications proposed by previous studies based on morphological data and molecular data are highly discordant. In addition, the phylogenetic relationships of major clades in the Angelica group, particularly in the Angelica s. s. clade, remain unclear. Chloroplast (cp) genome sequences have been widely used in phylogenetic studies and for evaluating genetic diversity.

RESULTS

In this study, we sequenced and assembled 28 complete cp genomes from 22 species, two varieties and two cultivars of Angelica. Combined with 36 available cp genomes in GenBank from representative clades of the subfamily Apioideae, the characteristics and evolutionary patterns of Angelica cp genomes were studied, and the phylogenetic relationships of Angelica species were resolved. The Angelica cp genomes had the typical quadripartite structure including a pair of inverted repeats (IRs: 5836-34,706 bp) separated by a large single-copy region (LSC: 76,657-103,161 bp) and a small single-copy region (SSC: 17,433-21,794 bp). Extensive expansion and contraction of the IR region were observed among cp genomes of Angelica species, and the pattern of the diversification of cp genomes showed high consistency with the phylogenetic placement of Angelica species. Species of Angelica were grouped into two major clades, with most species grouped in the Angelica group and A. omeiensis and A. sinensis grouped in the Sinodielsia with Ligusticum tenuissimum.

CONCLUSIONS

Our results further demonstrate the power of plastid phylogenomics in enhancing the phylogenetic reconstructions of complex genera and provide new insights into plastome evolution across Angelica L.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae

Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

Plastome comparison and evolution within the tribes of Plantaginaceae: Insights from an Asian gypsyweed

In spite of availability of several plastomes representing different tribes of Plantaginaceae, sparse attempts have been made to understand the plastome structure, evolution, and phylogenomics. In the present study, we have made an effort to understand the gene content and plastome evolution in the family Plantaginaceae using the newly generated plastome sequence of Veronica ovata subsp. kiusiana, a taxon native to SE Asia. In the first-ever attempt, plastomes of seven out of 10 tribes of Plantaginaceae have been compared to understand the evolution across the tribes of Plantaginaceae. The size of the plastome of V. ovata subsp. kiusiana is 152,249 bp, showing a typical quadripartite structure containing LSC, SSC, and two IRs with the sizes of 83,187, 17,704, and 25,679 respectively. The plastome comparison revealed the unique deletions in ycf2 and ndhF genes of members of different tribes, and also revealed high nucleotide variable hotspots. The study also revealed six highly variable genes and intergenic spacer viz. rps16, rps15-ycf1, ccsA-ndhD, ndhC-trnV, petN-psbM, and ycf1-trnN as potential DNA barcodes for the genus Veronica. The phylogenomic study revealed the sister relationship between V. ovata subsp. kiusiana and V. persica and also suggested the tentative placement of seven tribes in the family Plantaginaceae.

The plastomes of Astrocaryum aculeatum G. Mey. and A. murumuru Mart. show a flip-flop recombination between two short inverted repeats

The plastomes of Astrocaryum murumuru and A. aculeatum revealed a lineage-specific structural feature originated by flip-flop recombination, non-synonymous substitutions in conserved genes and several molecular markers. Astrocaryum murumuru Mart. and A. aculeatum G.Mey. are two palm species of Amazon forest that are economically important as source of food, oil and raw material for several applications. Genetic studies aiming to establish strategies for conservation and domestication of both species are still in the beginning given that the exploitation is mostly by extractive activity. The identification and characterization of molecular markers are essential to assess the genetic diversity of natural populations of both species. Therefore, we sequenced and characterized in detail the plastome of both species. We compared both species and identified 32 polymorphic SSR loci, 150 SNPs, 46 indels and eight hotspots of nucleotide diversity. Additionally, we reported a specific RNA editing site found in the ccsA gene, which is exclusive to A. murumuru. Moreover, the structural analysis in the plastomes of both species revealed a 4.6-kb inversion encompassing a set of genes involved in chlororespiration and plastid translation. This 4.6-kb inversion is a lineage-specific structural feature of the genus Astrocaryum originated by flip-flop recombination between two short inverted repeats. Furthermore, our phylogenetic analysis using whole plastomes of 39 Arecaceae species placed the Astrocaryum species sister to Acrocomia within the tribe Cocoseae. Finally, our data indicated substantial changes in the plastome structure and sequence of both species of the genus Astrocaryum, bringing new molecular markers, several structural and evolving features, which can be applied in several areas such as genetic, evolution, breeding, phylogeny and conservation strategies for both species.

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.

Plastomes of the green algae Hydrodictyon reticulatum and Pediastrum duplex (Sphaeropleales, Chlorophyceae)

Background

Comparative studies of chloroplast genomes (plastomes) across the Chlorophyceae are revealing dynamic patterns of size variation, gene content, and genome rearrangements. Phylogenomic analyses are improving resolution of relationships, and uncovering novel lineages as new plastomes continue to be characterized. To gain further insight into the evolution of the chlorophyte plastome and increase the number of representative plastomes for the Sphaeropleales, this study presents two fully sequenced plastomes from the green algal family Hydrodictyaceae (Sphaeropleales, Chlorophyceae), one from Hydrodictyon reticulatum and the other from Pediastrum duplex.

Methods

Genomic DNA from Hydrodictyon reticulatum and Pediastrum duplex was subjected to Illumina paired-end sequencing and the complete plastomes were assembled for each. Plastome size and gene content were characterized and compared with other plastomes from the Sphaeropleales. Homology searches using BLASTX were used to characterize introns and open reading frames (orfs) ≥ 300 bp. A phylogenetic analysis of gene order across the Sphaeropleales was performed.

Results

The plastome of Hydrodictyon reticulatum is 225,641 bp and Pediastrum duplex is 232,554 bp. The plastome structure and gene order of H. reticulatum and P. duplex are more similar to each other than to other members of the Sphaeropleales. Numerous unique open reading frames are found in both plastomes and the plastome of P. duplex contains putative viral protein genes, not found in other Sphaeropleales plastomes. Gene order analyses support the monophyly of the Hydrodictyaceae and their sister relationship to the Neochloridaceae.

Discussion

The complete plastomes of Hydrodictyon reticulatum and Pediastrum duplex, representing the largest of the Sphaeropleales sequenced thus far, once again highlight the variability in size, architecture, gene order and content across the Chlorophyceae. Novel intron insertion sites and unique orfs indicate recent, independent invasions into each plastome, a hypothesis testable with an expanded plastome investigation within the Hydrodictyaceae.

Discovery of novel plastid phenylalanine (trnF) pseudogenes defines a distinctive clade in Solanaceae

BACKGROUND

The plastome of embryophytes is known for its high degree of conservation in size, structure, gene content and linear order of genes. The duplication of entire tRNA genes or their arrangement in a tandem array composed by multiple pseudogene copies is extremely rare in the plastome. Pseudogene repeats of the trnF gene have rarely been described from the chloroplast genome of angiosperms.

FINDINGS

We report the discovery of duplicated copies of the original phenylalanine (trnFGAA) gene in Solanaceae that are specific to a larger clade within the Solanoideae subfamily. The pseudogene copies are composed of several highly structured motifs that are partial residues or entire parts of the anticodon, T- and D-domains of the original trnF gene.

CONCLUSIONS

The Pseudosolanoid clade consists of 29 genera and includes many economically important plants such as potato, tomato, eggplant and pepper.