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.

Highly Diverse Shrub Willows (Salix L.) Share Highly Similar Plastomes

Plastome phylogenomics is used in a broad range of studies where single markers do not bear enough information. Phylogenetic reconstruction in the genus Salix is difficult due to the lack of informative characters and reticulate evolution. Here, we use a genome skimming approach to reconstruct 41 complete plastomes of 32 Eurasian and North American Salix species representing different lineages, different ploidy levels, and separate geographic regions. We combined our plastomes with published data from Genbank to build a comprehensive phylogeny of 61 samples (50 species) using RAxML (Randomized Axelerated Maximum Likelihood). Additionally, haplotype networks for two observed subclades were calculated, and 72 genes were tested to be under selection. The results revealed a highly conserved structure of the observed plastomes. Within the genus, we observed a variation of 1.68%, most of which separated subg. Salix from the subgeneric Chamaetia/Vetrix clade. Our data generally confirm previous plastid phylogenies, however, within Chamaetia/Vetrix phylogenetic results represented neither taxonomical classifications nor geographical regions. Non-coding DNA regions were responsible for most of the observed variation within subclades and 5.6% of the analyzed genes showed signals of diversifying selection. A comparison of nuclear restriction site associated DNA (RAD) sequencing and plastome data on a subset of 10 species showed discrepancies in topology and resolution. We assume that a combination of (i) a very low mutation rate due to efficient mechanisms preventing mutagenesis, (ii) reticulate evolution, including ancient and ongoing hybridization, and (iii) homoplasy has shaped plastome evolution in willows.

Plastome and Nuclear Phylogenies of Dwarf Mistletoes (Arceuthobium: Viscaceae)

Abstract— Dwarf mistletoes are a lineage of morphologically-reduced stem parasites inhabiting Pinaceae and Cupressaceae throughout the northern hemisphere and equatorial East Africa. Though diagnosable by a suite of morphological traits, phylogenetic knowledge of species relationships has been limited to studies employing either comprehensive taxonomic sampling of one or two genes, or more sequence data from a limited number of individuals. We used data from genome skimming to assemble 3kb of the nuclear ribosomal cistron and up to 45kb of the plastome to clarify the phylogenetic root of the genus, monophyly of species, and relationships among infraspecific taxa. Genomic differentiation among terminal taxa was variable; however, we found strong support for reciprocally monophyletic New World and Old World lineages, congruent nrDNA and plastome topologies at the species level and below, and monophyly of most taxonomic sections and species. Plastome gene content was stable across the genus with minimal pseudogenization or loss, as in other hemiparasites, with the notable exception of cemA. These findings form the basis of our re-evaluation of historical biogeographical hypotheses, species- and subspecies-level taxonomy, and plastome evolution in Arceuthobium. More broadly, this work provides a foundation for future clade-focused comparative and biosystematics studies of Arceuthobium.

Plastome Evolution in the Hyperdiverse Genus Euphorbia (Euphorbiaceae) Using Phylogenomic and Comparative Analyses: Large-Scale Expansion and Contraction of the Inverted Repeat Region

With c. 2,000 species, Euphorbia is one of the largest angiosperm genera, yet a lack of chloroplast genome (plastome) resources impedes a better understanding of its evolution. In this study, we assembled and annotated 28 plastomes from Euphorbiaceae, of which 15 were newly sequenced. Phylogenomic and comparative analyses of 22 plastome sequences from all four recognized subgenera within Euphorbia revealed that plastome length in Euphorbia is labile, presenting a range of variation c. 42 kb. Large-scale expansions of the inverted repeat (IR) region were identified, and at the extreme opposite, the near-complete loss of the IR region (with only 355 bp left) was detected for the first time in Euphorbiaceae. Other structural variations, including gene inversion and duplication, and gene loss/pseudogenization, were also observed. We screened the most promising molecular markers from both intergenic and coding regions for phylogeny-based utilities, and estimated maximum likelihood and Bayesian phylogenies from four datasets including whole plastome sequences. The monophyly of Euphorbia is supported, and its four subgenera are recovered in a successive sister relationship. Our study constitutes the first comprehensive investigation on the plastome structural variation in Euphorbia and it provides resources for phylogenetic research in the genus, facilitating further studies on its taxonomy, evolution, and conservation.

Plastome evolution and organisation in the Hoya group (Apocynaceae)

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

The chloroplast genome evolution of Venus slipper (Paphiopedilum): IR expansion, SSC contraction, and highly rearranged SSC regions

Background Paphiopedilum is the largest genus of slipper orchids. Previous studies showed that the phylogenetic relationships of this genus are not well resolved, and sparse taxon sampling documented inverted repeat (IR) expansion and small single copy (SSC) contraction of the chloroplast genomes of Paphiopedilum. Results Here, we sequenced, assembled, and annotated 77 plastomes of Paphiopedilum species (size range of 152,130 – 164,092 bp). The phylogeny based on the plastome resolved the relationships of the genus except for the phylogenetic position of two unstable species. We used phylogenetic and comparative genomic approaches to elucidate the plastome evolution of Paphiopedilum. The plastomes of Paphiopedilum have a conserved genome structure and gene content except in the SSC region. The large single copy/inverted repeat (LSC/IR) boundaries are relatively stable, while the boundaries of the inverted repeat and small single copy region (IR/SSC) varied among species. Corresponding to the IR/SSC boundary shifts, the chloroplast genomes of the genus experienced IR expansion and SSC contraction. The IR region incorporated one to six genes of the SSC region. Unexpectedly, great variation in the size, gene order, and gene content of the SSC regions was found, especially in the subg. Parvisepalum. Furthermore, Paphiopedilum provides evidence for the ongoing degradation of the ndh genes in the photoautotrophic plants. The estimated substitution rates of the protein coding genes show accelerated rates of evolution in clpP, psbH, and psbZ. Genes transferred to the IR region due to the boundary shift also have higher substitution rates. Conclusions We found IR expansion and SSC contraction in the chloroplast genomes of Paphiopedilum with dense sampling, and the genus shows variation in the size, gene order, and gene content of the SSC region. This genus provides an ideal system to investigate the dynamics of plastome evolution.

Comparison and Phylogenetic Analyses of Nine Complete Chloroplast Genomes of Zingibereae

Zingibereae is a large tribe in the family Zingiberaceae, which contains plants with important medicinal, edible, and ornamental values. Although tribes of Zingiberaceae are well circumscribed, the circumscription of many genera within Zingibereae and the relationships among them remain elusive, especially for the genera of Boesenbergia, Curcuma, Kaempferia and Pyrgophyllum. In this study, we investigated the plastome variation in nine species representing five genera of Zingibereae. All plastomes showed a typical quadripartite structure with lengths ranging from 162,042 bp to 163,539 bp and contained 132–134 genes, consisting of 86–88 coding genes, 38 transfer RNA genes and eight ribosomal RNA genes. Moreover, the characteristics of the long repeats sequences and simple sequence repeats (SSRs) were detected. In addition, we conducted phylogenomic analyses of the Zingibereae and related taxa with plastomes data from additional 32 species from Genbank. Our results confirmed that Stahlianthus is closely related to Curcuma, supporting the idea of merging it into Curcuma. Kaempferia, Boesenbergia and Zingiber were confirmed as close relatives and grouped together as the Kaempferia group. Pyrgophyllum is not allied with the Curcuma clade but instead is embedded within the Hedychium clade. Our results demonstrate the power of plastid phylogenomics in improving the phylogenetic relationships within Zingibereae and provide a new insight into plastome evolution in Zingibereceae.

Plastome structure and phylogenetic relationships of Styracaceae (Ericales)

Background The Styracaceae are a woody, dicotyledonous family containing 12 genera and an estimated 160 species. Recent studies have shown that Styrax and Sinojackia are monophyletic, Alniphyllum and Bruinsmia cluster into a clade with an approximately 20-kb inversion in the Large Single-Copy (LSC) region. Halesia and Pterostyrax are not supported as monophyletic, while Melliodendron and Changiostyrax always form sister clades. Perkinsiodendron and Changiostyrax are newly established genera of Styracaceae. However, the phylogenetic relationship of Styracaceae at the generic level needs further research. Results We collected 28 complete plastomes of Styracaceae, including 12 sequences newly reported here and 16 publicly available sequences, comprising 11 of the 12 genera of Styracaceae. All species possessed the typical quadripartite structure of angiosperm plastomes, with sequence differences being minor, except for a large 20-kb (14 genes) inversion found in Alniphyllum and Bruinsmia. Seven coding sequences (rps4, rpl23, accD, rpoC1, psaA, rpoA and ndhH) were identified to possess positively selected sites. Phylogenetic reconstructions based on seven data sets (i.e., LSC, SSC, IR, Coding, Non-coding, combination of LSC + SSC and concatenation of LSC + SSC + one IR) produced similar topologies. In our analyses, all genera were strongly supported as monophyletic. Styrax was sister to the remaining genera. Alniphyllum and Bruinsmia form a clade. Halesia diptera does not cluster with Perkinsiodendron, while Perkinsiodendron and Rehderodendron form a clade. Changiostyrax is sister to a clade of Pterostyrax and Sinojackia. Conclusion Overall, our results demonstrate the power of plastid phylogenomics in improving estimates of phylogenetic relationships among genera. This study also provides insight into plastome evolution across Styracaceae.