http://orcid.org/0000-0002-7016-7485
Figures
Abstract
In Rubiaceae phylogenetics, the number of markers often proved a limitation with authors failing to provide well-supported trees at tribal and generic levels. A robust phylogeny is a prerequisite to study the evolutionary patterns of traits at different taxonomic levels. Advances in next-generation sequencing technologies have revolutionized biology by providing, at reduced cost, huge amounts of data for an increased number of species. Due to their highly conserved structure, generally recombination-free, and mostly uniparental inheritance, chloroplast DNA sequences have long been used as choice markers for plant phylogeny reconstruction. The main objectives of this study are: 1) to gain insight in chloroplast genome evolution in the Rubiaceae (Ixoroideae) through efficient methodology for de novo assembly of plastid genomes; and, 2) to test the efficiency of mining SNPs in the nuclear genome of Ixoroideae based on the use of a coffee reference genome to produce well-supported nuclear trees. We assembled whole chloroplast genome sequences for 27 species of the Rubiaceae subfamily Ixoroideae using next-generation sequences. Analysis of the plastid genome structure reveals a relatively good conservation of gene content and order. Generally, low variation was observed between taxa in the boundary regions with the exception of the inverted repeat at both the large and short single copy junctions for some taxa. An average of 79% of the SNP determined in the Coffea genus are transferable to Ixoroideae, with variation ranging from 35% to 96%. In general, the plastid and the nuclear genome phylogenies are congruent with each other. They are well-resolved with well-supported branches. Generally, the tribes form well-identified clades but the tribe Sherbournieae is shown to be polyphyletic. The results are discussed relative to the methodology used and the chloroplast genome features in Rubiaceae and compared to previous Rubiaceae phylogenies.
Citation: Ly SN, Garavito A, De Block P, Asselman P, Guyeux C, Charr J-C, et al. (2020) Chloroplast genomes of Rubiaceae: Comparative genomics and molecular phylogeny in subfamily Ixoroideae. PLoS ONE 15(4): e0232295. https://doi.org/10.1371/journal.pone.0232295
Editor: Shilin Chen, Chinese Academy of Medical Sciences and Peking Union Medical College, CHINA
Received: December 20, 2019; Accepted: April 11, 2020; Published: April 30, 2020
Copyright: © 2020 Ly et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All the files are available from the Genbank database (accession numbers: MN851267-MN851274 and MK577905-MK577918).
Funding: The author(s) received no specific funding for this work.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Rubiaceae (coffee family) belongs to Gentianales in the eudicots. It is the fourth most species-rich and diverse family in the flowering plants [1, 2, https://stateoftheworldsplants.org/2017/], comprising ca. 13,600 species grouped in ca. 620 genera and ca. 60 tribes [2, 3]. Rubiaceae are mainly tropical trees and shrubs, and less often annual or perennial herbs [4]. They occupy a large range of ecological niches from desert to evergreen humid forests and from sea level to high altitudes (above 4,000 m [5]). While some herbaceous species reached the temperate regions, Rubiaceae are especially abundant (species diversity and biomass) in lowland humid tropical forest, where they often are the most species-abundant of the woody plant families [2]. The Rubiaceae are divided into two subfamilies, Rubioideae and Cinchonoideae by [1], whereas Bremer and Eriksson [6] recognized three subfamilies, splitting the Cinchonoideae into Ixoroideae and Cinchonoideae.
The pantropical Ixoroideae subfamily comprises ca. 4,000 species [7], distributed into 27 tribes [7, 8, 9], and several well-known genera, i.e. the economically important Coffea and the horticulturally important Gardenia and Ixora [10], besides other less economically important genera such as Vangueria, Alibertia and Duroia L.f.
Molecular phylogenetic analyses of Rubiaceae have been carried out using either nuclear sequences (ETS, ITS, 5S-NTS, pep-C large, pep-V small, PI, Tpi), plastid DNA sequences (accD-psa1, atpB-rbcL, ndhF, matK, petD, rbcL, rpl16, rpl32-trnL, rps16, trnG, trnH-psbA, trnL-F, trnT-L, trnS-G) or a combination of both [7, 11, 12, 13, 14]. Altogether, more than twenty markers (fourteen from cpDNA and seven nuclear) have been used for Rubiaceae phylogeny reconstruction, the most popular being ITS and rbcL. However, the actual number of amplicons used in individual studies is much lower, e.g. for the dating of the family, subfamily and tribes based on fossils, only five plastid sequences were used [6]. The number of markers used often proved a limitation at tribal and generic levels as authors failed to provide well-supported trees [15, 16]. For instance, Maurin and coworkers [15], using four plastid regions and the internal transcribed spacer (ITS) region of nuclear rDNA (ITS 1/5.8S/ITS 2), failed to get a robust molecular tree for Coffea. Similarly, using one plastid and one nuclear marker, Khan and coworkers [17] re-circumscribed the Sabiceeae tribe but were unable to perform a proper biogeographic analysis, given that the molecular tree was largely unresolved. The availability of a robust phylogeny is a prerequisite to accurately study trait evolution at different taxonomic levels (family, subfamily, tribe or genus). This is the case, for instance, when mapping morphological and functional traits in Gardenieae [9], when investigating the evolution of sexual systems and growth habit in Mussaenda [18], or when studying the evolution of caffeine content in Coffea [19]. Since two decades, advances in next-generation sequencing (NGS) technologies have revolutionized the field of biology by providing, at reduced cost, huge amounts of data for an increased number of plant species. Among them, short-read sequencing technologies occupy an important place as they need relatively small amounts of DNA (from 600 ng to 1 μg), which allows the use of limited quantities of initial material, such as from herbarium samples [20]. Sequencing on total DNA permits to reconstruct whole chloroplast (cp) genome sequences of around 150–170 kb [21, 22, 23, 24], which can be used to construct robust phylogenies.
Chloroplasts are derived from endosymbiosis between independent living cyanobacteria and a non-photosynthetic eukaryotic host [25, 26]. Most flowering plants, including Coffea species [24] and Emmenopterys henryi [23, 27, 28] have a quadripartite circular chloroplast structure with two copies of Inverted Repeat (IR) regions (further called IRA and IRB) separating two regions of unique DNA sequence named large single copy (LSC) and small single copy (SSC) according to their length and gene composition [21]. The comparison of the structure and gene composition in cp genomes in broad sets of organisms permits to better understand their origin and function [29]. Due to their highly conserved structure, generally recombination-free, and mostly uniparental inheritance, cp DNA sequences have long been used as choice markers for plant phylogeny [30, 31, 32]. However, the low degree of polymorphism among the regular DNA markers used for Rubiaceae phylogenetics often does not resolve relationships at genus level in case of recent speciation [10]. In such conditions, the use of whole cp genome sequences could be a good alternative. In Rubiaceae, complete cp genomes are available for three species of two tribes in subfamily Cinchonoideae (tribe Naucleeae: Mitragyna speciosa Korth. [33], Neolamarckia cadamba (Roxb.) Bosser [34]; tribe Guettardeae: Antirhea chinensis (Champ. ex Benth.) Benth. & Hook.f. ex F.B.Forbes & Hemsl. [35]), five species belonging to at least three tribes in subfamily Rubioideae (tribe Spermacoceae: Hedyotis ovata Thunb. ex Maxim. [36]; insertae sedis: Paralasianthus hainanensis (Merr.) H.Zhu (as Saprosma merrillii H.S.Lo; [37]; tribe Rubieae: Galium mollugo L. (NC_028009); tribe Morindeae: Gynochthodes officinalis (F.C.How) Razafim. & B.Bremer (as Morinda officinalis F.C.How; NC_028009), Gynochthodes nanlingensis (Y.Z.Ruan) Razafim. & B.Bremer (NC_028614)], and two species belonging to subfamily Ixoroideae (tribe Condamineae: Emmenopterys henryi [23] and tribe Scyphiphoreae: Scyphiphora hydrophyllaceae [38]). However, large projects aiming to develop a library of plastid genomes (including Rubiaceae) are ongoing [33 (GenomeTrakrCP project), 39]).
Nuclear genomic raw data can be assembled into short contigs and used to mine Single Nucleotide Polymorphisms (SNPs) to study the genetic diversity within and between populations and species [40], the evolution of traits of interest [19] or the dynamics of transposable elements [41]. Methodologies based upon short read sequencing such as Genotyping-By-Sequencing (GBS) using a reference genome, permit to define sets of nuclear SNPs for high numbers of genotypes (convenient for multiples of 96 well-plates) as was done for Coffea species [19]. This is possible even with low nuclear genome coverage sequencing (about 10 x coverage). The combination of independent whole genome short read sequencing and bioinformatics tools permit to search these SNPs in different sets of species.
The main objectives of this study were i)- to develop efficient methodology to obtain complete de novo assembled cp genomes permitting comparative genomics and a robust molecular phylogenetic tree, ii)- to test the efficiency of mining SNPs in the nuclear genome of non-coffee Rubiaceae based on the use of a coffee reference genome in order to produce a well-supported nuclear tree, and, iii)- to gain insight in chloroplast genome evolution in the Rubiaceae.
Material and methods
Material
For this study, we have limited the sampling to subfamily Ixoroideae, to which also Coffea belongs. We included 27 taxa representing 10 tribes (Coffeeae, Condamineae, Cordiereae, Gardenieae, Ixoreae, Mussaendeae, Octotropideae, Pavetteae, Sherbournieae and Vanguerieae) plus Emmenopterys henryi [23], the complete cp sequence of which was retrieved from NCBI. Detailed information on sampling is given in Table 1.
Our analyses resulted in a single sample, Sherbournia, with a phylogenetic position different from what was expected. In order to verify the identity of this sample, TrnL-F and rps16 sequences were blasted in GenBank. Blasting was then repeated with sequences from other Sherbournia samples obtained with Sanger sequencing. Samples used were S. bignoniiflora (Welw.) Hua [Boyekoli Ebale 283 (BR)], S. buccularia [Lachenaud et al. 730 (BR)] and S. zenkeri Hua [Dessein et al. 1428 (BR)].
Authors of genus and species names of the studied taxa are given in Table 1; for other taxa they are given in the text upon first use.
Genome sequencing
Whole genomic DNA was isolated from silica or living plant material following a modified cetyltrimethylammonium bromide (CTAB) method [42]. A total of 25 mg (silica dried) or 100 mg (fresh) leaf material was ground into a fine powder. To eliminate secondary metabolites, two consecutive chloroform cleaning steps were carried out. DNA was lysed either in elution buffer (10 mM Tris-HCl, pH 8.0–8.5) or sterile PCR-grade water. The use of EDTA in elution buffer should be avoided to circumvent possible enzymatic inhibition in downstream applications which may lead to lower library quality. The short-read sequencing was done using the BGI-seq 500 platform, 2x100 bp paired-end. The quality of reads was verified using the Java software FASTQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). The raw data were checked in order to detect potential contamination using Kraken and Krona tools [43]. Raw reads were cleaned when necessary using Trimmomatic [44].
Chloroplast genome assembly and annotation
Cp genomes were de novo assembled using NOVOplasty software [45] from the sorted cp raw sequences obtained. Good quality raw reads were split into two sets corresponding to forward (F) and reverse (R) reads. With the aim to sort only cp data, the two sets of data (F and R) were mapped against the Coffea arabica cp reference genome [46] using Bowtie2 [47]. The choice of C. arabica (the genetically best-known species among Rubiaceae) at this step is justified since one objective of this study is to test the transferability of tools and methodology developed for Coffea to other members of Rubiaceae. Then, in silico filtered reads were considered for the cp genome de novo assembly using NOVOplasty. The recalcitrant cp genomes were re-assembled using Abyss [48]. At the end of the assembly process, the cp genomes were compared with the reference C. arabica genome using Gepard [49] in order to check for incongruences of the assembly. To confirm the overall structure, the pair-end illumina reads are mapped back to the assembly using Bowtie2 [47] and BAM files are used to display the read coverage using Artemis (https://www.sanger.ac.uk/science/tools/artemis-comparison-tool-act). The cp genomes were then annotated using Geseq [50] as recommended by Guyuex et al., [24]. The circular visualization of each cp genome was obtained using the OrganellarGenomeDRAW tool (OGDRAW) [51]. The linear gene order comparison was obtained using ACT (Artemis Comparison Tool [52]).
Sequence divergence and junction sequences divergence
Given that Coffea arabica could not be used as outgroup in the phylogeny, we decided to use the annotated genomes of Antirhea chinensis [35], Mitragyna speciosa [33] and Neolamarckia cadamba [34] belonging to the Cinchonoideae subfamily as reference genomes and outgroup taxa for the cp genome analyses. The alignments of the complete chloroplast genome sequences of the 28 studied Rubiaceae were visualized using mVISTA [53] in order to show global interspecific variation and variation within the tribes.
Taking into account data obtained from for each taxon (length of regions LSC, SSC, IR and gene annotation), we calculated the distance between boundaries and the nearest gene to visualize junction sequence divergence between species and within tribes.
Phylogenetic relationships
The plastid phylogeny was produced using a total of 28 cp sequences and one of three outgroups belonging to the Cinchonoideae subfamily (Antirhea chinensis, Mitragyna speciosa and Neolamarckia cadamba, retrieved from GenBank). The sequences were first aligned using MAFFT 7.305 with the following parameters BLOSUM62 and 200PAM/ k = 2 [54].
The nuclear tree was produced using a total of 27 taxa (no data available for Emmenopterys henryi). No non-Ixoroideae data were available, so the trees were rooted midpoint. The 28,800 SNPs used for Coffea [19] were mined according to the methodology developed by these authors. In a second step, in order to reduce the amount of missing data, rare or too common sites were removed using Tassel ver. 5.0 [55] with the following parameters: minimum count = 18, minimum frequency = 0.2, maximum frequency = 0.8. A total of 1,726 sites (SNPs) were kept.
MAFFT alignment of cp sequences and nuclear SNP concatenation were used to infer the phylogenetic trees. The phylogenetic reconstructions were done using RAxML version 8 [56] under the General Time Reversible nucleotide substitution model with gamma distributed rate variation among sites (GTR+G), ML estimate of alpha-parameter, BFGS method to optimize GTR rate parameters and Felsenstein’s bootstraps option autoMRE as recommended by the author). The trees were then edited with FigTree ver. 1.3.1 [57] and Inkscape (https://inkscape.org/fr/release/0.91/).
Results and discussion
Chloroplast genome features in Ixoroideae
Among the 28 studied samples 25 exhibited the classical quadripartite structure but three had an apparent tripartite structure with only one IR (Mussaenda pubescens, Feretia aeruginescens and Pavetta schumanniana). These latter belong to three different tribes (Mussaendeae, Octotropideae and Pavetteae, respectively) but the tripartite structure is not present in all representatives of these tribes.
Regarding the quadripartite genomes, total length ranges from 153,056 bp for Bertiera breviflora to 155,328 bp for Sherbournia buccularia. Similar length differences are observed among the tripartite genomes (from 127,396 bp for Pavetta schumanniana to 129,508 bp for Mussaenda pubescens). Length variations were also noted for the different regions: from 83,406 bp (Vangueria infausta) to 85,461 bp (Belonophora coffeoides) for LSC, from 17,915 bp (Mitriostigma axillare) to 18,245 bp (Emmenopterys henryi) for SSC and from 24,855 bp (Bertiera laxa) to 25,978 bp (Mussaenda pubescens) for IR. In all species, GC content was similar for the complete cp genome (ca. 37%) as well as for each of the cp subregions (LSC: ca. 35%; SSC: ca. 31%; IR: ca. 43%). Individual information is summarized in Table 2.
We annotated a total of 118 different genes belonging to 14 functional categories and present in all the genomes with the exception of Tarenna grevei which has lost trnH-GUG. One hundred genes were present as single copy, 17 were duplicated and one (rps12) was triplicated in IR. The LSC, SSC and IR regions contained 87, 13 and 18 genes respectively. Among the 80 protein-coding genes identified, only nine include introns. Seven of these contain one intron (atpF, ndhA & B, rpoC1, rps12, rps16, rpl2) whereas clpP and ycf3 have two introns. Complete infA and pbf1 genes were present in all studied species. For the three taxa showing only one IR, the corresponding genes were present in only one copy.
The annotated cp genome sequences permitted to compare the length of the junctions of the main regions LSC, IR and SSC among the studied Rubiaceae (Table 3). Generally, low variation was observed between the taxa in the boundary regions. However, while the distance for LSC/IRB junctions was 88 bp for most species, variation was noted in Emmenopterys henryi (Condamineae) with 30 bp, in Psilanthus ebracteolatus (Coffeeae) with 358 bp, in Coffea sessiliflora (Coffeeae) with 157 bp and in Sherbournia buccularia (Gardenieae) with 148 bp. A similar variation was obtained for the IRB/SSC junction. (S1 Table)
Estimations were done from data extracted from [33, 34, 35, 59] or calculated from data extracted from GenBank for Morinda officinalis (NC_028009) and Gynochtodes nanlingensis (NC_028614).
Data obtained for other Rubiaceae such as Hedyotis ovata [36] and Paralasianthus hainanensis (as Saprosma merrillii; [37]) from the Rubioideae subfamily and Antirhea chinensis [35] from the Cinchonoideae subfamily indicate a quadripartite structure and a total of 114 genes (eight duplicated genes counted once and eight genes missing in Rubiaceae, see below) of which 80 are unique protein-coding genes. For Neolamarckia cadamba (Cinchonoideae), [34] revealed a total of 130 genes, 79 of which are protein-coding. Data obtained from GenBank for three Rubioideae (accession numbers NC_036970 for Galium mollugo, NC_028009 for Morinda officinalis and NC_028614 for Gynochtodes nanlenginsis) and from literature for two Cinchonoideae [34, 35] permitted us to determine their plastid features (Table 3). All have the classical quadripartite structure. In the Rubioideae, the total cp length varies from 153,398 to 154,086 bp; LSC from 84,011 to 84,471 bp; SSC from 17,056 to 18,115 bp and IR from 25,766 to 26,075 bp. In the Cinchonoideae, total cp length varies from 154,999 to 155,616 bp; LSC from 85,880 to 86,252 bp; SSC from 17,851 to 17,984 bp and IR from 25,634 to 25,690 bp. The GC content for the total sequence and for the different regions are similar to our results in Ixoroideae. The chloroplast genome features of the Ixoroideae, and of the Rubiaceae as a whole, are in the ranges reported for most flowering plants [29, 30, 58].
The tripartite genome structure was not yet reported in the Rubiaceae but was recorded for Fabaceae [60], Geraniaceae [61], Pinaceae [62], Cactaceae [63], Arecaceae [64] and Passifloraceae [65]. Within the Ixoroideae, the chloroplast assembly of three species showed a tripartite genome structure. Besides frequent inversions, duplications, or losses of fragments [65, 66], IR expansion/contraction and even IR absence contributed to substantial variation in cp genome length [67]. However, the robustness of our assemblies of species showing only one IR were tested. All reads were mapped on the assembly and the read coverage was displayed. The read coverage showed an increase at the IR region of the assembly suggesting that two IR regions may be present but assembled into only one IR (S1 Fig). So, the assembly of these chloroplast sequences (Mussaenda pubescens, Feretia aeruginescens and Pavetta schumanniana) and the IR absence should be considered with caution since this event is very rare in most other plant families [68] and since we cannot exclude that the assembly process collapsed the IR regions into only one. Increasing the number of taxa investigated and using long read sequencing techniques, such as PacBio and Oxford Nanopore may demonstrate it to be a not so rare event in the Rubiaceae, which would lead to questions on the role of two or one IR in the evolution of land plants.
Circular visualization of the Bertiera breviflora cp genome is given in Fig 1 as an example. Gene order and orientation from pairwise comparisons were generally well-conserved although some gene orientations were different (S2 Fig).
Sequence divergence was visualized using mVISTA with Coffea arabica as the reference annotated genome. The choice of Coffea arabica instead of Antirhea chinensis (outgroup used for the plastid phylogeny) was justified by the level of divergence between Cinchonideae and Ixoroideae. Globally, sequence divergence among all taxa was relatively high and mainly concentrated in conserved non-coding sequences and in Untranslated Transcribed Regions (UTR). However, variation among species seemed to be negligible for UTRs located in the IR region (rpl2, ndhB and rps12 genes). Substitutions were more frequent but indels were observed as well, even in the ycf2 exon. Four (the conserved non-coding regions between matK and atpA, rpoB and psbD, rps4 and ndhJ and ndhC and atpE), and two (ndhF—ccsA and ycf1) hypervariable regions were identified in the LSC and SSC regions respectively. A representation of sequence divergence is given for a selected set of taxa (Fig 2). In total, 31.5% sites of the complete alignment included indels.
The location and orientation of the genes are indicated on the top. Exons and UTRs are in purple and turquoise respectively. Conserved non-coding regions are in orange. The y-axis ranges from 100% (top) to 50% identity between each sequence and the reference. The order of the taxa used from top to bottom is: Bertiera iturensis (Bertiereae), Psilanthus ebracteolatus (Coffeeae), Empogona congesta (Coffeeae), Euclinia longiflora (Gardenieae), Pavetta schumanniana (Pavetteae), Pseudomussaenda stenocarpa (Mussaendeae), Feretia aeruginescens (Octotropideae), Mitriostigma axillare (Sherbournieae) and Alibertia edulis (Cordiereae). Feretia aeruginescens and Pavetta schumanniana showing only one IR in the current assembly.
In Ixoroideae, plotting the length variation of the different regions relative to the smallest cp genome (here Vangueria infausta and considering only the quadripartite cp genomes), showed the pattern of variation given in Fig 3. The length of the different regions did not increase simultaneously to the total cp length except for the smallest four cp genomes. The increase in size seems to be mainly due to increase in length of LSC and possibly to gene and/or intron length increases. For Bertiera breviflora and Bertiera laxa, the increase in cp size is mainly due to the increase in length of LSC associated with a decrease in length of IR. Variation in the length of SSC has only limited impact on cp size variation. Regarding the four Cinchonoideae and the five Rubioideae species, similar patterns of variation are observed with the exception of Galium mollugo for which a decrease in length of SSC is notable. Therefore, with the exception of a few species, it seems that length variation in LSC is the main contributor to cp size variation. Among eudicots, the progressive expansion of the IR has been documented in Pelargonium L'Hér. ex Aiton [69] and Passiflora L. [65], and a similar molecular mechanism driving the IR evolution in these two unrelated lineages could be a possibility.
The taxa are ordered in increasing total cp genome size in each subfamily (24 Ixoroideae, four Cinchonideae in blue box marked with one asterisk and five Rubioideae in red box marked with two asterisks). Data for taxa indicated with asterisk(s) was retrieved from literature [23, 33, 34, 35, 36, 37] or calculated from data extracted from GenBank for Morinda officinalis (NC_028009), Gallium mollugo (NC_036970), Gynochtodes nanlingensis (NC_028614).
Angiosperm cp genomes exhibit a remarkably conserved gene content and order as observed for instance within Fagaceae [21, 70, 71, 72] and more specifically for Quercus L. [22]. Likewise, gene content and order were nearly identical in the Ixoroideae representatives studied as well as in representatives of the two other Rubiaceae subfamilies.
Recorded in tobacco and in most others members of Solanaceae as a pseudogene [46], infA was intact in all Ixoroideae species of this study and in the Rubioideae and Cinchonoideae species for which whole cp genomes are available. Similarly, putatively involved in photosystem I and II biogenesis, pbf1 (psbN in Coffea arabica, [46]) was present in all Ixoroideae as well as in the Rubioideae and Cinchonoideae. In all Ixoroideae of this study and in the species of the other subfamilies, a fragment of rps19 appeared duplicated at the IR/LSC boundaries as reported in Solanaceae with the exception of tobacco [73]. Eight genes (CHLB, CHLL, CHLN, CYSA, CYST, MBPX, PSAM, and RPL21) were absent in the study of 16 wild coffee trees [24]. This study showed their absence in all Ixoroideae and the other Rubioideae and Cinchonoideae tested. Finally, despite minor changes in gene content, orientation and order, Ixoroideae plastid genomes are well conserved within and between tribes. This was also the case in the available Cinchonoideae and Rubioideae species and, therefore, could be true for the whole family. However, sequence divergence within and between tribes was observed and at much higher level (Fig 2) than reported in Quercus [22].
Plastid molecular phylogeny and comparison to previous Rubiaceae phylogenies
The complete cp genome-based phylogeny included 28 Ixoroideae taxa and Antirhea chinensis (Cinchonoideae subfamily) as outgroup. Maximum Likelihood analyses resulted in a generally well-resolved topology with highly supported branches, except for four lineages: the branch between Empogona and the Belonophora/Tricalysia clade, the branch between Leptactina and Pavetta/Tarenna within the Pavetteae tribe, the branch towards the Cordiereae/Octotropideae clade and the branch towards the Mussaendeae/Condamineae clade (BS < 80%, Fig 4). The ingroup has three main clades: Mussaendeae (Fig 4; in green), Condamineae (Fig 4; in red) and a large clade comprising all other taxa. The Mussaendeae and Condamineae are well-supported as distinct monophyletic lineages (BS = 100) but their mutual relationship and their relationship with the rest of the ingroup remain unclear. The rest of the ingroup forms a well-supported clade (BS = 100) and comprises two well-supported subclades (BS = 100) that correspond to the Vanguerieae alliance (Ixora and Vangueria) and the Coffeeae alliance. Within the Coffeeae alliance, the tribe Pavetteae (Fig 4; in pink), the Coffeeae/Bertiereae lineage with Bertiera sister to the Coffeeae (Fig 4; in blue) and the clade comprising Schumanniophyton, Gardenia, Sherbournia, Euclina and Atractocarpus (Gardenieae, Fig 4; in brown) are supported as monophyletic groups (BS = 100). All tribes represented by at least two representatives are retrieved as monophyletic with the exception of the Sherbournieae. Sherbournia does not form a clade with Mitriostigma but is firmly embedded in Gardenieae. The same phylogeny was obtained with Mitragyna speciosa (Cinchonoideae) as outgroup. However, when using Neolamarckia cadamba as outgroup, a slightly different phylogenetic tree was obtained (data not shown).
Four well-supported clades are marked in green for Mussaendeae, red for Condamineae, pink for Pavetteae, brown for Gardenieae and blue for Coffeeae/Bertiereae.
This chloroplast phylogeny concurs well with previously published phylogenetic trees based on Sanger sequencing of several markers [1, 6, 7]. Within the Ixoroideae (Ixoridinae sensu Robbrecht and Manen), Robbrecht and Manen [1] recognized a basal clade (basal Ixoridinae; not represented in our analysis) and two main evolutionary lineages Ixoridinae I and Ixoridinae II. Ixoridineae I is essentially neotropical and represented here by the tribe Condamineae. Ixoridineae II is mainly paleotropical and includes all other members of the ingroup. Unlike the super-tree of Robbrecht and Manen [1], our plastid phylogeny is not resolved at the base and does not clearly separate Ixoridineae I and II, since Mussaendeae is considered part of Ixoridineae II by [1]. Bremer and Eriksson [6] did not distinguish lineages within the subfamily Ixoroideae, probably because the base of their phylogenetic tree is unresolved. Kainulainen et al., [7] recognized within the Ixoroideae a basal grade (here represented by Condamineae and Mussaendeae) and a clade of core Ixoroideae. Our analysis confirms the subfamilial classification of Kainulainen et al., [7] rather than that of Robbrecht and Manen [1]. Within the core Ixoroideae, Kainulainen et al., [7] differentiated between the Vanguerieae alliance and the Coffeeae alliance. These two clades are also retrieved in our analysis. Our analysis confirms the sister relationships between Ixoreae and Vanguerieae [1, 6, 7] and between Coffeeae and Bertiereae [1, 6, 7, 9].
Gardenieae is retrieved as monophyletic only with the inclusion of Sherbournia, which has been considered part of the tribe Sherbournieae [9]. The Sherbournieae were recently instated [9] to include the former Gardenieae genera Sherbournia, Mitriostigma, Atractogyne and Oxyanthus, the last two of which are not included in our analysis. With the exception of Sherbournia, these genera are characterized by pollen grains in tetrads [73]. Persson and more recently Bremer and Eriksson [74, 6] also retrieved this group of three genera with pollen in tetrads as monophyletic. However, the inclusion of Sherbournia makes the tribe morphologically heterogeneous as regards to pollen characters (pollen in monads). In order to check the identity of our Sherbournia sample we separated TrnL-F and rps16 sequences from the whole genomes sequence and blasted them in GenBank, where they showed more similarity with Rothmannia Thunb. than with the Sherbournia sequences present there. This was repeated with sequences from other Sherbournia species obtained with Sanger sequencing with the same results. We are therefore confident that Sherbournia does not form part of the tribe Sherbournieae as delimited by Darwin [9], but belongs to the Gardenieae. It should be noted that also in the phylogeny of Persson [74], Sherbournia groups with Rothmannia. The tribe Gardenieae has been demonstrated in several studies to be polyphyletic [1, 7, 9]. The fact that it is not so in our analysis is the result of the small number of representatives included, notably five genera out of over fifty [9]. The five genera making up the Gardenieae clade in our analysis are not generally considered closely related. Atractocarpus is part of Gardenieae IV in [1] and of the Porterandia group in [9], Gardenia is part of Gardenieae II [1] and the Gardenia group [9], Euclinia belongs to Gardenieae III [1] and the Randia group [9], Schumanniophyton belongs to Gardenieae I [1] and remains unplaced in [9] and Sherbournia is unplaced in [1] and belongs to the Sherbournieae tribe in [9].
Nuclear SNP mining and Efficiency of transferability of methods from Coffea to Ixoroideae
The genome of Coffea canephora Pierre ex A.Froehner was used as reference genome to mine SNPs as described in Hamon et al. (2017). This methodology was efficient despite unequal results between taxa. No outgroup from another subfamily of the Rubiaceae was available so the tree was rooted midpoint. In this analysis Coffea canephora was added but Emmenopterys henryi could not be included since no nuclear genome data was available.
An average of 22,906 SNPs was sorted with the extremes ranging from 10,335 in Pentagonia tinajita to 27,642 in Tricalysia lasiodelphys. Among the 806,400 individual data expected (28 x 28,800), excluding all Coffea species, the average percentage of missing data was 31% ranging from 10% in Tricalysia hensii to 77% in Ixora chinensis and Pentagonia tinajita. The percentage of heterozygotes was 0.6% on average but varied from 0.25% in Atractocarpus fitzalanii to 3.2% in Psilanthus ebracteolatus. The nucleotide percentage was 29.1% for A, 29.3% for T, 20.5% for G and 20.2% for C (S1 Appendix). So, the SNP transferability from Coffea [19] to non-coffee Rubiaceae belonging to ten tribes of subfamily Ixoroideae can be considered as successful. Interestingly, the phylogenetically most distant species are those with the fewest orthologous sequences.
The species relationships obtained with the complete dataset (28,800 sites) are shown in Fig 5. The Maximum Likelihood tree shows a majority of well-supported branches (BS of 86–100%). The ingroup shows two main, well-supported clades, the first comprising the Coffeeae alliance and the second comprising the Vanguerieae alliance, the Condamineae and the Mussaendeae. The Mussaendeae (Fig 5; in green), Pavetteae (Fig 5; in pink), Coffeeae/Bertiereae (Fig 5; in blue) and Gardenieae (Fig 5; in brown) are supported as monophyletic with high branch support values. Bertiera is sister to the Coffeeae. The following results are in contrast to the results of the plastid phylogeny: the Gardenieae clade (Fig 5; in brown) does not include Schumanniophyton; Ixora and Vangueria (Vanguerieae alliance) do not form a monophyletic group. The monophyly of the Condamineae cannot be evaluated because only a single representative is present in this analysis (no data for Emmenopterys henryi). The tribe Sherbournieae (Sherbournia and Mitriostigma) is not retrieved as monophyletic and it is embedded in the Gardenieae clade.
Colored clades indicated well-defined tribes. The tree is rooted midpoint as no outgroup is available. Green for Mussaendeae; pink for Pavetteae; brown for Gardenieae; blue for Coffeeae and turquoise for Coffeeae/Bertiereae.
The phylogenetic tree resulting from the SNP mining of the nuclear genome is similar to the chloroplast based phylogenetic tree with the same clades (Mussaendeae, Pavetteae, Bertiereae, Gardenieae, Coffeeae and Coffeeae/Bertiereae) being retrieved and highly supported even though the position of individual taxa within the clades may be different. Other relationships, such as the sister relationship between Ixoreae and Vanguerieae and between the Vanguerieae alliance and the Coffeeae alliance, are not retrieved in the nuclear phylogeny.
With the aim to use a dataset with less missing data, SNPs were filtered leading to a total of 1,726 sites (SNPs) retained for further analysis. The resulting tree (S3 Fig) shows a long branch for the Psilanthus-Coffea clade that may indicate a highly divergent evolution between these species and the rest of the ingroup. The tree further differs from the one based on 28.800 SNPs in that Bertiera is not sister to the Coffeeae but to the ingroup clade consisting of all species except for Coffeeae. Similarly to the 28,800 SNPs-based tree, branch support values are generally high. The clades Pavetteae, Mussaendeae, Gardenieae (excluding Schumanniophyton), Bertiereae and Coffeeae are supported as monophyletic. However, while the main clades are similar, their relative position is not the same in the two analyses. Sherbournieae are not retrieved as monophyletic, indicating that the reduction of the number of SNPs should be done with care due to possible bias in markers genomic distribution.
Conclusions
In this study we reported and analyzed the chloroplast genome sequences for 27 species of the Rubiaceae subfamily Ixoroideae using next-generation sequences (NGS. Plastid and nuclear genome phylogenies are well congruent with each other with an overall well-supported branch. Generally, the tribes form well-identified clades but the tribe Sherbournieae is shown to be polyphyletic. With continuously dropping prices and an increasing output and efficiency of bioinformatic tools, NGS appears to be now the best choice to study difficult or neglected plant families, tribes or genera. Our methodology used here combined plastid genome reconstruction and SNP mining of the nuclear genome and was successful for Ixoroideae. The same methodology should be extended to the two other Rubiaceae subfamilies (Cinchonoideae and Rubioideae). This would permit to clarify the relationships between Rubiaceae taxa and to better understand genome evolution in the family in relation to adaptive traits. The increased availability of more reference genomes other than Coffea genomes will facilitate and speed up this process.
Supporting information
S1 Table. Junction sequence divergence among 28 Rubiaceae.
Taxa are ordered alphabetically within tribes. The genes considered at the border between the main regions LSC, IR and SSC are those identified in this study. The distances between genes and junctions are given in bp. IRA is not confirmed the assembly for three species (Mussaenda pubescens, Feretia aeruginescens and Pavetta schumanniana).
https://doi.org/10.1371/journal.pone.0232295.s001
(XLSX)
S1 Fig. Illumina read coverage in the IR region of Mussaenda pubescens.
https://doi.org/10.1371/journal.pone.0232295.s002
(TIFF)
S2 Fig. Gene order and orientation visualized in some pairwise comparisons using Artemis Comparison Tool.
https://doi.org/10.1371/journal.pone.0232295.s003
(TIFF)
S3 Fig. Maximum likelihood nuclear tree of 28 Ixoroideae based on 1,726 nuclear SNPs (RaxML with GTR model of substitution) and bootstrap values to estimate branch supports.
Colored clades indicate well-defined tribes. The tree is rooted midpoint since no outgroup outside Ixoroideae is available.
https://doi.org/10.1371/journal.pone.0232295.s004
(TIFF)
S1 Appendix. Fasta sequences of assembled cp genomes.
https://doi.org/10.1371/journal.pone.0232295.s005
(TXT)
Acknowledgments
We thank the gardeners, Mr. H. Leqeux, Mr. J. Van Eeckhoudt, and the scientific curators, Ms. V. leyman and Dr. M. Reynders, for their care for the living Rubiaceae collection at Meise Botanic Garden. We thank Mr. F. Van Caekenberghe for help with the collection of plant material. We are grateful to the manager of the molecular laboratory at Meise Botanic Garden, Mr. W. Baert, for his goodwill towards this project. We thank Meise Botanic Garden for funding sequencing work, travel and publication fees. We thank the Institut de Recherche pour le Développement (IRD) for funding sequencing work and for providing a fellowship for S.N. Ly. The Laboratory Chrono-Environment, Université de Bourgogne Franche-Comté, also funded part of the sequencing.
References
- 1. Robbrecht E, Manen J-F. The major evolutionary lineages of the coffee family (Rubiaceae, angiosperms). Combined analysis (nDNA and cpDNA) to infer the position of Coptosapelta and Luculia, and supertree construction based on rbcL, rps16, trnL-trnF and atpB-rbcL data. A new classification in two subfamilies, Cinchonoideae and Rubioideae. Syst Geogr Plants. 2006; 76: 85–145.
- 2. Davis AP, Govaerts R, Bridson DM, Rushsam M, Moat J, Brummitt NA. A global assessment of distribution, diversity, endemism, and taxonomic effort in the Rubiaceae. Ann Mo Bot Gard. 2009; 96: 68–78.
- 3. Govaerts R, Ruhsam M, Andersson L, Robbrecht E, Bridson DM, Davis AP, et al. World checklist of Rubiaceae. 2016. Available from http://apps.kew.org/wcsp/
- 4. Robbrecht E. Tropical woody Rubiaceae. Characteristic features and progressions. Contributions to a new subfamilial classification. Opera Bot Belg. 1988; 1: 1–271.
- 5. Wikström N, Kainulainen K, Razafimandimbison SG, Smedmark JE, Bremer B. A revised time tree of the asterids: Establishing a temporal framework for evolutionary studies of the coffee family (Rubiaceae). PLOS ONE. 2015; 10, e0126690. pmid:25996595
- 6. Bremer B, Eriksson T. Time tree of Rubiaceae: phylogeny and dating the family, subfamilies, and tribes. Int J Plant Sci. 2009; 170: 766–793.
- 7. Kainulainen K, Razafimandimbison SG, Bremer B. Phylogenetic relationships and new tribal delimitations in subfamily Ixoroideae (Rubiaceae). Bot J Linn Soc. 2013; 173: 387–406.
- 8. Darwin SP. The subfamilial, tribal and subtribal nomenclature of the Rubiaceae. Taxon. 1976; 25: 595–610.
- 9. Mouly A, Kainulainen K, Persson C, Davis AP, Wong KM, Razafimandimbison SG, et al. Phylogenetic structure and clade circumscriptions in the Gardenieae complex (Rubiaceae). Taxon. 2014; 63: 801–818.
- 10. Andreasen K, Bremer B. Combined phylogenetic analysis in the Rubiaceae-Ixoroideae: Morphology, nuclear and chloroplast DNA data. Am J Bot. 2000; 87: 1731–1748. pmid:11080124
- 11. Bremer B. A review of molecular phylogenetic studies of Rubiaceae. Ann Mo Bot Gard. 2009; 96: 4–26.
- 12. Tosh J, Davis AP, Dessein S, De Block P, Chester M, Maurin O, et al. Phylogeny of Tricalysia A.Rich. (Rubiaceae) and its relationships with allied genera based on cp DNA data. Ann Mo Bot Gard. 2009; 96: 194–213.
- 13. Cristians S, Bye R, Nieto-Sotelo J. Molecular markers associated with chemical analysis: A powerful tool for quality control assessment of copalchi medicinal plant complex. Front Pharmacol. 2018; 9: 666. pmid:29988415
- 14. De Block P, Rakotonasolo F, Ntore S, Razafimandimbison SG, Janssens S. Four new endemic genera of Rubiaceae (Pavetteae) from Madagascar represent multiple radiations into drylands. Phytokeys. 2018; 99: 1–66. pmid:29861651
- 15. Maurin O, Davis AP, Chester M, Mvungi EF, Jaufeerally-Fakim Y, Fay MF. Towards a phylogeny for Coffea (Rubiaceae): Identifying well-supported lineages based on nuclear and plastid DNA sequences. Ann Bot. 2007; 100: 1565–1583. pmid:17956855
- 16. Mouly A, Razafimandimbison SG, Florence J, Jérémie J, Bremer B. Paraphyly of Ixora and new tribal delimitation of Ixoreae (Rubiaceae): Inference from combined chloroplast (rps16, rbcl, and trnt-f) sequence data. Ann Mo Bot Gard. 2009; 96: 146–160.
- 17. Khan S, Razafimandimbison SG, Bremer B, Liede-Schumann S. Sabiceeae and Virectarieae (Rubiaceae, Ixoroideae): One or two tribes?–New tribal and generic circumscriptions of Sabiceeae and biogeography of Sabicea s.l. Taxon. 2008; 57: 7–23.
- 18. Duan T, Deng X, Chen S, Luo Z, Zhao Z, Tu T, et al. Evolution of sexual systems and growth habit in Mussaenda (Rubiaceae): Insights into the evolutionary pathways of dioecy. Mol Phylogenet Evol. 2018; 123: 113–122. pmid:29454889
- 19. Hamon P, Grover CE, Davis AP, Rakotomalala J-J, Raharimalala NE, Albert VA, et al. Genotyping-by-sequencing provides the first well-resolved phylogeny for coffee (Coffea) and insights into the evolution of caffeine content in its species: GBS coffee phylogeny and the evolution of caffeine content. Mol Phylogenet Evol. 2017; 109: 351–361. pmid:28212875
- 20. Staats M, Erkens RHJ, van de Vossenberg B, Wieringa JJ, Kraaijeveld K, et al. Genomic treasure troves: Complete genome sequencing of herbarium and insect museum specimens. PLOS ONE. 2013; 8: e69189. pmid:23922691
- 21. Jansen RK, Saski C, Lee S-B, Hansen AK, Daniell H. Complete plastid genome sequences of three rosids (Castanea, Prunus, Theobroma): Evidence for at least two independent transfers of rpl22 to the nucleus. Mol Biol Evol. 2011; 28: 835–847. pmid:20935065
- 22. Yang Y, Zhou T, Duan D, Yang J, Feng L, Zhao G. Comparative analysis of the complete chloroplast genomes of five Quercus species. Front Plant Sci. 2016; 7: 959. pmid:27446185
- 23. Duan R, Huang M, Yang L, Liu Z. Characterization of the complete chloroplast genome of Emmenopterys henryi (Gentianales: Rubiaceae), an endangered relict tree species endemic to China. Conserv Genet Resour. 2017; 9: 459–461.
- 24. Guyeux J, Charr JC, Hue TMT, Furtado A, Henry RB, Crouzillat D, et al. Evaluation of chloroplast genome annotation tools and application to analysis of the evolution of coffee species. PLOS ONE. 2019; 14: e0216347. pmid:31188829
- 25. McFadden GI. Chloroplast origin and integration. Plant Physiol. 2001; 125: 503. pmid:11154294
- 26. Keeling PJ. Diversity and evolutionary history of plastids and their hosts. Am J Bot. 2004; 91: 1481–1493. pmid:21652304
- 27. Jansen RK, Raubeson LA, Boore JL, dePamphilis CW, Chumley TW, et al. Methods for obtaining and analyzing whole chloroplast genome sequences. Methods Enzymol. 2005; 395: 348–384. pmid:15865976
- 28.
Jansen RK, Ruhlman TA. Plastid genomes of seed plants. In: Bock R., Knoop V, editors. Genomics of chloroplasts and mitochondria. Advances in photosynthesis and respiration (including bioenergy and related processes). Dordrecht: Springer; 2012; 103–126.
- 29. Wicke S, Schneeweiss GM, dePamphilis CW, Müller KF, Quandt D. The evolution of the plastid chromosome in land plants: gene content, gene order, gene function. Plant Mol Biol. 2011; 76: 273–297. pmid:21424877
- 30. Jansen RK, Cai Z, Raubeson LA, Daniell H, dePamphillis CW, Leebens-Mack J, et al. Analysis of 81 genes from 64 plastid genomes resolves relationships in angiosperms and identifies genome scale evolutionary patterns. Proc Natl Acad Sci USA. 2007; 104: 19369–19374. pmid:18048330
- 31. Parks M, Cronn R, Liston A. Increasing phylogenetic resolution at low taxonomic levels using massively parallel sequencing of chloroplast genomes. BMC Biol. 2009; 7: 84. pmid:19954512
- 32. Moore MJ, Soltis PS, Bell CD, Burleigh JG, Soltis DE. Phylogenetic analysis of 83 plastid genes further resolves the early diversification of eudicots. Proc Natl Acad Sci USA. 2010; 107: 4623–4628. pmid:20176954
- 33. Zhang N, Ramachandran P, Wen J, Duke JA, Metzman H, McLaughlin W, et al. Development of a reference standard library of chloroplast genome sequences, GenomeTrakrCP. Planta med 2017; 83: 1256–1263.
- 34. Li J, Zhang D, Ouyang K, Chen X. The complete chloroplast genome of the miracle tree Neolamarckia cadamba and its comparison in Rubiaceae family. Biotechnol. Biotechnol. Equip. 2018; 32: 1087–1097.
- 35. Fan W-W, Wang J-H, Zhao K-K, Wang H-X, Zhu Z-X, Wang H-F. Complete plastome sequence of Antirhea chinensis (Champ. ex Benth.) Forbes et Hemst: An endemic species in South China. Mitochondrial DNA B. 2019; 4: 538–539.
- 36. Zhang X-F, Wang J-H, Zhao K-K, Fan W-W, Wang H-X, Zhu Z-X, et al. Complete plastome sequence of Hedyotis ovata Thunb. ex Maxim (Rubiaceae): an endemic shrub in Hainan, China. Mitochondrial DNA B. 2019; 4: 675–676.
- 37. Zhu Z-X, Wang J-H, Zhao K-K, Fan W-W, Wang H-X, Wang H-F. Complete chloroplast genome of Saprosma merrillii Lo (Rubiaceae): A Near Threaten (NT) shrub species endemic to Hainan province, China. Mitochondrial DNA B. 2019; 4: 742–743.
- 38. Zhang Y, Zhang J-W, Yang Y, Li X-N. Structural and Comparative Analysis of the Complete Chloroplast Genome of a Mangrove Plant: Scyphiphora hydrophyllacea Gaertn. f. and Related Rubiaceae Species. Forests 2019, 10(11), 1000.
- 39. Rydin C, Wikström N, Bremer B. Conflicting results from mitochondrial genomic data challenge current views of Rubiaceae phylogeny. Am J Bot. 2017; 104: 1522–1532. pmid:29885222
- 40. Curk F, Ancillo G, Perrier X, Jacquemoud-Collet J-P, Garcia-Lor A, Navarro L, et al. Nuclear species-diagnostic SNP markers mined from 454 amplicon sequencing reveal admixture genomic structure of modern Citrus varieties. PLOS ONE. 2015; 10, e0125628. pmid:25973611
- 41. Guyot R, Darré T, Dupeyron M, de Kochko A, Hamon S, Couturon E, et al. Partial sequencing reveals the transposable element composition of Coffea genomes and provides evidence for distinct evolutionary stories. Mol Genet Genomics. 2016; 291: 1979–1990. pmid:27469896
- 42. Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987; 19: 11–15.
- 43. Ondov BD, Bergman NH, Phillippy AM. Interactive metagenomic visualization in a Web browser. BMC Bioinformatics. 2011; 12: 385. pmid:21961884
- 44. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114–2120. pmid:24695404
- 45. Dierckxsens N, Mardulyn P, Smits G. NOVOPlasty: de novo assembly of organelle genomes from whole genome data. Nucleic Acids Res. 2017; 45, e18. pmid:28204566
- 46. Samson N, Bausher MG, Lee S-B, Jansen RK, Daniell H. The complete nucleotide sequence of the coffee (Coffea arabica L.) chloroplast genome: organization and implications for biotechnology and phylogenetic relationships amongst angiosperms. Plant Biotechnol J. 2007; 5: 339–353. pmid:17309688
- 47. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat. Methods. 2012; 9: 357–359. pmid:22388286
- 48. Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJM, Birol I. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009; 19: 1117–1123. pmid:19251739
- 49. Krumsiek J, Arnold R, Rattei T. Gepard: a rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics. 2007; 23: 1026–1028. pmid:17309896
- 50. Tillich M, Lehwark P, Pellizzer T, Ulbricht-Jones ES, Fischer A, Bock R, et al. GeSeq - versatile and accurate annotation of organelle genomes. Nucleic Acids Res. 2017; 45, W6–W11. pmid:28486635
- 51. Greiner S, Lehwark P, Bock R. OrganellarGenomeDRAW (OGDRAW) version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes. Nucleic Acids Res. 2019. 47: W59–W64. pmid:30949694
- 52. Carver TJ, Rutherford KM, Berriman M, Rajandream M-A, Barrell BG, Parkhill J. ACT: the Artemis Comparison Tool. Bioinformatics. 2005; 21: 3422–3423. pmid:15976072
- 53. Frazer KA, Pachter L, Poliakov A, Rubin EM, Dubchak I. VISTA: computational tools for comparative genomics. Nucleic acids res. 2004; 32: W273–W279. pmid:15215394
- 54. Kuraku S, Zmasek CM, Nishimura O, Katoh K. aLeaves facilitates on-demand exploration of metazoan gene family trees on MAFFT sequence alignment server with enhanced interactivity. Nucleic Acids Res. 2013; 41: W22–W28. pmid:23677614
- 55. Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007; 23: 2633–2635. pmid:17586829
- 56. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014; 30: 1312–1313. pmid:24451623
- 57.
Rambaut A. FigTree. 2007. Available from http://tree.bio.ed.ac.uk/software/figtree/
- 58. Kim J-H, Lee S-I, Kim B-R, Choi I-Y, Ryser P, Kim N-S. Chloroplast genomes of Lilium lancifolium, L. amabile, L. callosum, and L. philadelphicum: Molecular characterization and their use in phylogenetic analysis in the genus Lilium and other allied genera in the order Liliales. PLOS ONE. 2017; 12: e0186788. pmid:29065181
- 59.
Dann M, Bellot S, Schepella S, Schaefer H, Tellier A. National Center for Biotechnology Information, NIH, Bethesda, MD 20894, USA. Galium mollugo NC_036970.1. 2018
- 60. Cai Z, Guisinger M, Kim H-G, Ruck E, Blazier JC, McMurtry V, et al. Extensive reorganization of the plastid genome of Trifolium subterraneum (Fabaceae) is associated with numerous repeated sequences and novel dna insertions. J Mol Evol. 2008; 67: 696–704. pmid:19018585
- 61. Guisinger MM, Kuehl JV, Boore JL, Jansen RK. Extreme reconfiguration of plastid genomes in the angiosperm family Geraniaceae: rearrangements, repeats, and codon usage. Mol Biol Evol. 2010; 28: 583–600. pmid:20805190
- 62. Lin CP, Huang JP, Wu CS, Hsu CY, Chaw SM. Comparative chloroplast genomics reveals the evolution of Pinaceae genera and subfamilies. Genome Biol Evol. 2010; 2: 504–517. pmid:20651328
- 63. Sanderson MJ, Copetti D, Búrquez A, Bustamante E, Charboneau JL, Eguiarte LE, et al. Exceptional reduction of the plastid genome of saguaro cactus (Carnegiea gigantea): loss of the ndh gene suite and inverted repeat. Am J Bot. 2015; 102: 1115–1127. pmid:26199368
- 64. Barrett CF, Baker WJ, Comer JR, Conran JG, Lahmeyer SC, Leebens-Mack , et al. Plastid genomes reveal support for deep phylogenetic relationships and extensive rate variation among palms and other commelinid monocots. New Phytol. 2016; 209: 855–870. pmid:26350789
- 65. Shrestha B, Weng M-L, Theriot EC, Gilbert LE, Ruhlman TA, Krosnick SE, et al. Highly accelerated rates of genomic rearrangements and nucleotide substitutions in plastid genomes of Passiflora subgenus Decaloba. Mol Phylogenet Evol. 2019; 138: 53–64. pmid:31129347
- 66. Zhai W, Duanb X, Zhang R, Guo C, Li L, Xu G, et al. Chloroplast genomic data provide new and robust insights into the phylogeny and evolution of the Ranunculaceae. Mol Phylogenet Evol. 2019; 135: 12–21. pmid:30826488
- 67.
Ruhlman TA, Jansen RK. The plastid genomes of flowering plants. In: Maliga P, editor. Chloroplast biotechnology: methods and protocols. Totowa, NJ: Humana Press; 2014. pp. 3–38.
- 68. Zheng X-M, Junrui W, Li F, Sha L, Hongbo P, Lan Q, et al. Inferring the evolutionary mechanism of the chloroplast genome size by comparing whole chloroplast genome sequences in seed plants. Sci Rep. 2017; 7: 1555. pmid:28484234
- 69. Weng M-L, Ruhlman TA, Jansen RK. Expansion of inverted repeat does not decrease substitution rates in Pelargonium plastid genomes. New Phytol. 2017; 214: 842–851. pmid:27991660
- 70. Alexander LW, Woeste KE. Pyrosequencing of the northern red oak (Quercus rubra L.) chloroplast genome reveals high quality polymorphisms for population management. Tree Genet. Genomes. 2014; 10: 803–812.
- 71. Dane F, Wang Z, Goertzen L. Analysis of the complete chloroplast genome of Castanea pumila var. pumila, the Allegheny chinkapin. Tree Genet. Genomes. 2015; 11: 1–6.
- 72. Lu S, Hou M, Du FK, Li J, Yin K. Complete chloroplast genome of the Oriental whiteoak: Quercus aliena Blume. Mitochondrial DNA. 2015; 27: 2802–2804. pmid:26114324
- 73. Chung H-J, Jung JD, Park H-W, Kim J-H, Cha HW, Min SR, et al. The complete chloroplast genome sequence of Solanum tuberosum and comparative analysis with Solanaceae species identified the presence of a 241-bp deletion in cultivated potato chloroplast DNA sequence. Plant Cell Rep. 2006; 25: 1369–1379. pmid:16835751
- 74. Persson C. Phylogeny of Gardenieae (Rubiaceae) based on chloroplast DNA sequences from the rps16 intron and trnL(UAA)‐F(GAA) intergenic spacer. Nordic J Bot. 2000; 20: 257–269.