Comparison of nuclear, plastid, and mitochondrial phylogenies and the origin of wild octoploid strawberry species

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Phylogenomic conflict analyses of the plastid and mitochondrial genomes via deep genome skimming highlight their independent evolutionary histories: A case study in the cinquefoil genus Potentilla sensu lato (Potentilleae, Rosaceae)

Phylogenomic conflicts are widespread among genomic data, with most previous studies primarily focusing on nuclear datasets instead of organellar genomes. In this study, we investigate phylogenetic conflict analyses within and between plastid and mitochondrial genomes using Potentilla as a case study. We generated three plastid datasets (coding, noncoding, and all-region) and one mitochondrial dataset (coding regions) to infer phylogenies based on concatenated and multispecies coalescent (MSC) methods. Conflict analyses were then performed using PhyParts and Quartet Sampling (QS). Both plastid and mitochondrial genomes divided the Potentilla into eight highly supported clades, two of which were newly identified in this study. While most organellar loci were uninformative for the majority of nodes (bootstrap value < 70%), PhyParts and QS detected conflicting signals within the two organellar genomes. Regression analyses revealed that conflict signals mainly occurred among shorter loci, whereas longer loci tended to be more concordant with the species tree. In addition, two significant disagreements between the two organellar genomes were detected, likely attributed to hybridization and/or incomplete lineage sorting. Our results demonstrate that mitochondrial genes can fully resolve the phylogenetic relationships among eight major clades of Potentilla and are not always linked with plastome in evolutionary history. Stochastic inferences appear to be the primary source of observed conflicts among the gene trees. We recommend that the loci with short sequence length or containing limited informative sites should be used cautiously in MSC analysis, and suggest the joint application of concatenated and MSC methods for phylogenetic inference using organellar genomes.

Phylogenomics and topological conflicts in the tribe Anthospermeae (Rubiaceae)

Genome skimming (shallow whole-genome sequencing) offers time- and cost-efficient production of large amounts of DNA data that can be used to address unsolved evolutionary questions. Here we address phylogenetic relationships and topological incongruence in the tribe Anthospermeae (Rubiaceae), using phylogenomic data from the mitochondrion, the nuclear ribosomal cistron, and the plastome. All three genomic compartments resolve relationships in the Anthospermeae; the tribe is monophyletic and consists of three major subclades. Carpacoce Sond. is sister to the remaining clade, which comprises an African subclade and a Pacific subclade. Most results, from all three genomic compartments, are statistically well supported; however, not fully consistent. Intergenomic topological incongruence is most notable in the Pacific subclade but present also in the African subclade. Hybridization and introgression followed by organelle capture may explain these conflicts but other processes, such as incomplete lineage sorting (ILS), can yield similar patterns and cannot be ruled out based on the results. Whereas the null hypothesis of congruence among all sequenced loci in the individual genomes could not be rejected for nuclear and mitochondrial data, it was rejected for plastid data. Phylogenetic analyses of three subsets of plastid loci identified using the hierarchical likelihood ratio test demonstrated statistically supported intragenomic topological incongruence. Given that plastid genes are thought to be fully linked, this result is surprising and may suggest modeling or sampling error. However, biological processes such as biparental inheritance and inter-plastome recombination have been reported and may be responsible for the observed intragenomic incongruence. Mitochondrial insertions into the plastome are rarely documented in angiosperms. Our results indicate that a mitochondrial insertion event in the plastid trnS GGA - rps4 IGS region occurred in the common ancestor of the Pacific clade of Anthospermeae. Exclusion/inclusion of this locus in phylogenetic analyses had a strong impact on topological results in the Pacific clade.

Target capture and genome skimming for plant diversity studies

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

"Omics" insights into plastid behavior toward improved carotenoid accumulation

Plastids are a group of diverse organelles with conserved carotenoids synthesizing and sequestering functions in plants. They optimize the carotenoid composition and content in response to developmental transitions and environmental stimuli. In this review, we describe the turbulence and reforming of transcripts, proteins, and metabolic pathways for carotenoid metabolism and storage in various plastid types upon organogenesis and external influences, which have been studied using approaches including genomics, transcriptomics, proteomics, and metabonomics. Meanwhile, the coordination of plastid signaling and carotenoid metabolism including the effects of disturbed carotenoid biosynthesis on plastid morphology and function are also discussed. The "omics" insight extends our understanding of the interaction between plastids and carotenoids and provides significant implications for designing strategies for carotenoid-biofortified crops.

Fragaria mitogenomes evolve rapidly in structure but slowly in sequence and incur frequent multinucleotide mutations mediated by microinversions

Plant mitochondrial DNA has been described as evolving rapidly in structure but slowly in sequence. However, many of the noncoding portions of plant mitogenomes are not homologous among species, raising questions about the rate and spectrum of mutations in noncoding regions. Recent studies have suggested that the lack of homology in noncoding regions could be due to increased sequence divergence. We compared 30 kb of coding and 200 kb of noncoding DNA from 13 sequenced Fragaria mitogenomes, followed by analysis of the rate of sequence divergence, microinversion events and structural variations. Substitution rates in synonymous sites and nongenic sites are nearly identical, suggesting that the genome-wide point mutation rate is generally consistent. A surprisingly high number of large multinucleotide substitutions were detected in Fragaria mitogenomes, which may have resulted from microinversion events and could affect phylogenetic signal and local rate estimates. Fragaria mitogenomes preferentially accumulate deletions relative to insertions and substantial genomic arrangements, whereas mutation rates could positively associate with these sequence and structural changes among species. Together, these observations suggest that plant mitogenomes exhibit low point mutations genome-wide but exceptionally high structural variations, and our results favour a gain-and-loss model for the rapid loss of homology among plant mitogenomes.

ORTHOSKIM: in silico sequence capture from genomic and transcriptomic libraries for phylogenomic and barcoding applications

Low-coverage whole genome shotgun sequencing (or genome skimming) has emerged as a cost-effective method for acquiring genomic data in non-model organisms. This method provides sequence information on chloroplast genome (cpDNA), mitochondrial genome (mtDNA) and nuclear ribosomal regions (rDNA), which are over-represented within cells. However, numerous bioinformatic challenges remain to accurately and rapidly obtain such data in organisms with complex genomic structures and rearrangements, in particular for mtDNA in plants or for cpDNA in some plant families. Here we introduce the pipeline ORTHOSKIM, which performs in silico capture of targeted sequences from genomic and transcriptomic libraries without assembling whole organelle genomes. ORTHOSKIM proceeds in three steps: 1) global sequence assembly, 2) mapping against reference sequences, and 3) target sequence extraction; importantly it also includes a range of quality control tests. Different modes are implemented to capture both coding and non-coding regions of cpDNA, mtDNA and rDNA sequences, along with predefined nuclear sequences (e.g. ultra-conserved elements) or collections of single-copy ortholog genes. Moreover, aligned DNA matrices are produced for phylogenetic reconstructions, by performing multiple alignments of the captured sequences. While ORTHOSKIM is suitable for any eukaryote, a case study is presented here, using 114 genome-skimming libraries and 4 RNAseq libraries obtained for two plant families, Primulaceae and Ericaceae, the latter being a well-known problematic family for cpDNA assemblies. ORTHOSKIM recovered with high success rates cpDNA, mtDNA and rDNA sequences, well suited to accurately infer evolutionary relationships within these families. ORTHOSKIM is released under a GPL-3 license and is available at: https://github.com/cpouchon/ORTHOSKIM.

A mitochondrial genome phylogeny of voles and lemmings (Rodentia: Arvicolinae): Evolutionary and taxonomic implications

Arvicolinae is one of the most impressive placental radiations with over 150 extant and numerous extinct species that emerged since the Miocene in the Northern Hemisphere. The phylogeny of Arvicolinae has been studied intensively for several decades using morphological and genetic methods. Here, we sequenced 30 new mitochondrial genomes to better understand the evolutionary relationships among the major tribes and genera within the subfamily. The phylogenetic and molecular dating analyses based on 11,391 bp concatenated alignment of protein-coding mitochondrial genes confirmed the monophyly of the subfamily. While Bayesian analysis provided a high resolution across the entire tree, Maximum Likelihood tree reconstruction showed weak support for the ordering of divergence and interrelationships of tribal level taxa within the most ancient radiation. Both the interrelationships among tribes Lagurini, Ellobiusini and Arvicolini, comprising the largest radiation and the position of the genus Dinaromys within it also remained unresolved. For the first time complex relationships between genus level taxa within the species-rich tribe Arvicolini received full resolution. Particularly Lemmiscus was robustly placed as sister to the snow voles Chionomys in the tribe Arvicolini in contrast with a long-held belief of its affinity with Lagurini. Molecular dating of the origin of Arvicolinae and early divergences obtained from the mitogenome data were consistent with fossil records. The mtDNA estimates for putative ancestors of the most genera within Arvicolini appeared to be much older than it was previously proposed in paleontological studies.

Plastid phylogenomics improves resolution of phylogenetic relationship in the Cheirostylis and Goodyera clades of Goodyerinae (Orchidoideae, Orchidaceae)

Goodyerinae are one of phylogenetically unresolved groups of Orchidaceae. The lack of resolution achieved through the analyses of previous molecular sequences from one or a few markers has long confounded phylogenetic estimation and generic delimitation. Here, we present large-scale phylogenomic data to compare the plastome structure of the two main clades (Goodyera and Cheirostylis) in this subtribe and further adopt two strategies, combining plastid coding sequences and the whole plastome, to investigate phylogenetic relationships. A total of 46 species in 16 genera were sampled, including 39 species in 15 genera sequenced in this study. The plastomes of heterotrophic species are not drastically reduced in overall size, but display a pattern congruent with a loss of photosynthetic function. The plastomes of autotrophic species ranged from 147 to 165 kb and encoded from 132 to 137 genes. Three unusual structural features were detected: a 1.0-kb inversion in the large single-copy region of Goodyera schlechtendaliana; the loss and/or pseudogenization of ndh genes only in two species, Cheirostylis chinensis and C. montana; and the expansion of inverted repeat regions and contraction of small single-copy region in Hetaeria oblongifolia. Phylogenomic analyses provided improved resolution for phylogenetic relationships. All genera were recovered as monophyletic, except for Goodyera and Hetaeria, which were each recovered as non-monophyletic. Nomenclatural changes are needed until the broader sampling and biparental inherited markers. This study provides a phylogenetic framework of Goodyerinae and insight into plastome evolution of Orchidaceae.

Complete chloroplast genome sequencing of ten wild Fragaria species in China provides evidence for phylogenetic evolution of Fragaria

Complete chloroplast genomes of ten wild Fragaria species native to China were sequenced. Phylogenetic analysis clustered Fragaria species into two clades: The south clade (F. iinumae, F. chinensis, F. pentaphylla, F. nilgerrensis, F. daltoniana, F. corymbosa, F. moupinensis, F. tibetica, F. nipponica, F. gracilis, and F. nubicola and north clade (F. viridis, F. orientalis, F. moschata, F. mandshurica, F. vesca, F. chiloensis, F. virginiana, and F. × ananassa), while F. iinumae is the oldest extant species. Molecular clock analysis suggested present Fragaria species share a common ancestor 3.57 million years ago (Ma), F. moschata and octoploid species evolve 0.89 and 0.97 Ma, respectively, but F. moschata be not directly involved in current octoploid species formation. Drastic global temperature change since the Palaeocene-Eocene, approx. 55 Ma, especially during uplifting of the Qinghai-Tibet plateau and quaternary glaciation may have driven the formation of Fragaria, separation of two groups and polyploidization.

Utility of the Mitochondrial Genome in Plant Taxonomic Studies

Size, structure, and sequence content lability of plant mitochondrial genome (mtDNA) across species has sharply limited its use in taxonomic studies. Historically, mtDNA variation has been first investigated with RFLPs, while the development of universal primers then allowed studying sequence polymorphisms within short genomic regions (<3 kb). The recent advent of NGS technologies now offers new opportunities by greatly facilitating the assembly of longer mtDNA regions, and even full mitogenomes. Phylogenetic works aiming at comparing signals from different genomic compartments (i.e., nucleus, chloroplast, and mitochondria) have been developed on a few plant lineages, and have been shown especially relevant in groups with contrasted inheritance of organelle genomes. This chapter first reviews the main characteristics of mtDNA and the application offered in taxonomic studies. It then presents tips for best sequencing protocol based on NGS data to be routinely used in mtDNA-based phylogenetic studies.

Resolving the Phylogeny of the Olive Family (Oleaceae): Confronting Information from Organellar and Nuclear Genomes

The olive family, Oleaceae, is a group of woody plants comprising 28 genera and ca. 700 species, distributed on all continents (except Antarctica) in both temperate and tropical environments. It includes several genera of major economic and ecological importance such as olives, ash trees, jasmines, forsythias, osmanthuses, privets and lilacs. The natural history of the group is not completely understood yet, but its diversification seems to be associated with polyploidisation events and the evolution of various reproductive and dispersal strategies. In addition, some taxonomical issues still need to be resolved, particularly in the paleopolyploid tribe Oleeae. Reconstructing a robust phylogenetic hypothesis is thus an important step toward a better comprehension of Oleaceae's diversity. Here, we reconstructed phylogenies of the olive family using 80 plastid coding sequences, 37 mitochondrial genes, the complete nuclear ribosomal cluster and a small multigene family encoding phytochromes (phyB and phyE) of 61 representative species. Tribes and subtribes were strongly supported by all phylogenetic reconstructions, while a few Oleeae genera are still polyphyletic (Chionanthus, Olea, Osmanthus, Nestegis) or paraphyletic (Schrebera, Syringa). Some phylogenetic relationships among tribes remain poorly resolved with conflicts between topologies reconstructed from different genomic regions. The use of nuclear data remains an important challenge especially in a group with ploidy changes (both paleo- and neo-polyploids). This work provides new genomic datasets that will assist the study of the biogeography and taxonomy of the whole Oleaceae.

The complete chloroplast genome sequence of the Dioscorea esculenta (Lour.) Burkill (Dioscoreaceae)

Dioscorea esculenta (Lour.) Burkill is an essential tuber crop with pharmacological effects in the family Dioscoreaceae. The complete chloroplast genome of D. esculenta was determined in this study. The total genome size is 153, 437 bp in length and demonstrates a typical quadripartite structure containing a large single copy (LSC, 83,628 bp) and a small single copy (SSC, 18,893 bp), separated by a pair of inverted repeats (IRa, IRb) of 25,458 bp. The GC content of the complete chloroplast genome sequence is 37.07%. A total of 131 genes were predicted including 86 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic tree analysis of 25 species belonging to the genus Dioscorea indicated that D. esculenta and D. sansibarensis were clustered into one branch.

Genome Synteny Has Been Conserved Among the Octoploid Progenitors of Cultivated Strawberry Over Millions of Years of Evolution

Allo-octoploid cultivated strawberry (Fragaria × ananassa) originated through a combination of polyploid and homoploid hybridization, domestication of an interspecific hybrid lineage, and continued admixture of wild species over the last 300 years. While genes appear to flow freely between the octoploid progenitors, the genome structures and diversity of the octoploid species remain poorly understood. The complexity and absence of an octoploid genome frustrated early efforts to study chromosome evolution, resolve subgenomic structure, and develop a single coherent linkage group nomenclature. Here, we show that octoploid Fragaria species harbor millions of subgenome-specific DNA variants. Their diversity was sufficient to distinguish duplicated (homoeologous and paralogous) DNA sequences and develop 50K and 850K SNP genotyping arrays populated with co-dominant, disomic SNP markers distributed throughout the octoploid genome. Whole-genome shotgun genotyping of an interspecific segregating population yielded 1.9M genetically mapped subgenome variants in 5,521 haploblocks spanning 3,394 cM in F. chiloensis subsp. lucida, and 1.6M genetically mapped subgenome variants in 3,179 haploblocks spanning 2,017 cM in F. × ananassa. These studies provide a dense genomic framework of subgenome-specific DNA markers for seamlessly cross-referencing genetic and physical mapping information and unifying existing chromosome nomenclatures. Using comparative genomics, we show that geographically diverse wild octoploids are effectively diploidized, nearly completely collinear, and retain strong macro-synteny with diploid progenitor species. The preservation of genome structure among allo-octoploid taxa is a critical factor in the unique history of garden strawberry, where unimpeded gene flow supported its origin and domestication through repeated cycles of interspecific hybridization.

Beyond balancing selection: frequent mitochondrial recombination contributes to high-female frequencies in gynodioecious Lobelia siphilitica (Campanulaceae)

Gynodioecy is a sexual system in which females and hermaphrodites co-occur. In most gynodioecious angiosperms, sex is determined by an interaction between mitochondrial male-sterility genes (CMS) that arise via recombination and nuclear restorer alleles that evolve to suppress them. In theory, gynodioecy occurs when multiple CMS types are maintained at equilibrium frequencies by balancing selection. However, some gynodioecious populations contain very high frequencies of females. High female frequencies are not expected under balancing selection, but could be explained by the repeated introduction of novel CMS types. To test for balancing selection and/or the repeated introduction of novel CMS, we characterised cytoplasmic haplotypes from 61 populations of Lobelia siphilitica that vary widely in female frequency. We confirmed that mitotype diversity and female frequency were positively correlated across populations, consistent with balancing selection. However, while low-female populations hosted mostly common mitotypes, high-female populations and female plants hosted mostly rare, recombinant mitotypes likely to carry novel CMS types. Our results suggest that balancing selection maintains established CMS types across this species, but extreme female frequencies result from frequent invasion by novel CMS types. We conclude that balancing selection alone cannot account for extreme population sex-ratio variation within a gynodioecious species.

Multi-gene incongruence consistent with hybridisation in Cladocopium (Symbiodiniaceae), an ecologically important genus of coral reef symbionts

Coral reefs rely on their intracellular dinoflagellate symbionts (family Symbiodiniaceae) for nutritional provision in nutrient-poor waters, yet this association is threatened by thermally stressful conditions. Despite this, the evolutionary potential of these symbionts remains poorly characterised. In this study, we tested the potential for divergent Symbiodiniaceae types to sexually reproduce (i.e. hybridise) within Cladocopium, the most ecologically prevalent genus in this family. With sequence data from three organelles (cob gene, mitochondrion; psbA^ncr region, chloroplast; and ITS2 region, nucleus), we utilised the Incongruence Length Difference test, Approximately Unbiased test, tree hybridisation analyses and visual inspection of raw data in stepwise fashion to highlight incongruences between organelles, and thus provide evidence of reticulate evolution. Using this approach, we identified three putative hybrid Cladocopium samples among the 158 analysed, at two of the seven sites sampled. These samples were identified as the common Cladocopium types C40 or C1 with respect to the mitochondria and chloroplasts, but the rarer types C3z, C3u and C1# with respect to their nuclear identity. These five Cladocopium types have previously been confirmed as evolutionarily distinct and were also recovered in non-incongruent samples multiple times, which is strongly suggestive that they sexually reproduced to produce the incongruent samples. A concomitant inspection of next generation sequencing data for these samples suggests that other plausible explanations, such as incomplete lineage sorting or the presence of co-dominance, are much less likely. The approach taken in this study allows incongruences between gene regions to be identified with confidence, and brings new light to the evolutionary potential within Symbiodiniaceae.

Employing hypothesis testing and data from multiple genomic compartments to resolve recalcitrant backbone nodes in Goodenia s.l. (Goodeniaceae)

Goodeniaceae is a primarily Australian flowering plant family with a complex taxonomy and evolutionary history. Previous phylogenetic analyses have successfully resolved the backbone topology of the largest clade in the family, Goodenia s.l., but have failed to clarify relationships within the species-rich and enigmatic Goodenia clade C, a prerequisite for taxonomic revision of the group. We used genome skimming to retrieve sequences for chloroplast, mitochondrial, and nuclear markers for 24 taxa representing Goodenia s.l., with a particular focus on Goodenia clade C. We performed extensive hypothesis tests to explore incongruence in clade C and evaluate statistical support for clades within this group, using datasets from all three genomic compartments. The mitochondrial dataset is comparable to the chloroplast dataset in providing resolution within Goodenia clade C, though backbone support values within this clade remain low. The hypothesis tests provided an additional, complementary means of evaluating support for clades. We propose that the major subclades of Goodenia clade C (C1-C3 + Verreauxia) are the result of a rapid radiation, and each represents a distinct lineage.

Genome skimming provides new insight into the relationships in Ludwigia section Macrocarpon, a polyploid complex

PREMISE OF THE STUDY

Interpreting relationships within groups containing polyploids, which are frequent in angiosperms, can be greatly assisted by genomic techniques. In this study, we used a genome-skimming approach to investigate the evolutionary relationships and origins of polyploids in the monophyletic group, Ludwigia section Macrocarpon (Onagraceae), which includes diploid, tetraploid, and hexaploid taxa.

METHODS

We sampled all known taxa and ploidy levels in the section and conducted shotgun sequencing. We assembled plastomes, mitochondrial sequences, and completed nuclear ribosomal regions, reconstructed phylogenies, and conducted comparative genomic analyses for plastomes to gain insights into the relationships among studied taxa.

KEY RESULTS

Within the section, results showed that the South American diploid taxa L. bonariensis and L. lagunae were closely related. We reported the first chromosome count (2n = 4× = 32) for L. neograndiflora, which is closely related to the two South American diploid taxa, although its exact origin remains unclear. The samples of the widespread, polyploid taxon L. octovalvis do not form a monophyletic group. Both tetraploid and hexaploid L. octovalvis lineages have originated more than once. At least one tetraploid in the L. octovalvis lineage may have been involved in the origins of hexaploids. One or more extinct/unsampled intermediate tetraploids in the L. octovalvis lineages had also likely been involved in the origins of hexaploids.

CONCLUSIONS

Genome skimming provided important insights into the complex evolutionary relationships within sect. Macrocarpon, but additional sampling and data from single-copy nuclear regions are necessary to further elucidate the origins of the polyploids in this section.

Plastid genomes reveal recurrent formation of allopolyploid Fragaria

PREMISE OF THE STUDY

Recurrent formation of polyploid taxa is a common observation in many plant groups. Haploid, cytoplasmic genomes like the plastid genome can be used to overcome the problem of homeologous genes and recombination in polyploid taxa. Fragaria (Rosaceae) contains several octo- and decaploid species. We use plastome sequences to infer the plastid ancestry of these taxa with special focus on the decaploid Fragaria cascadensis.

METHODS

We used genome skimming of 96 polyploid Fragaria samples on a single Illumina HiSeq 3000 lane to obtain whole plastome sequences. These sequences were used for phylogenetic reconstructions and dating analyses. Ploidy of all samples was inferred with flow cytometry, and plastid inheritance was examined in a controlled cross of F. cascadensis.

KEY RESULTS

The plastid genome phylogeny shows that only the octoploid F. chiloensis is monophyletic, all other polyploid taxa were supported to be para- or polyphyletic. The decaploid Fragaria cascadensis has biparental plastid inheritance and four different plastid donors. Diversification of the F. cascadensis clades occurred in the last 230,000 years. The southern part of its distribution range harbors considerably higher genetic diversity, suggestive of a potential refugium.

CONCLUSIONS

Fragaria cascadensis had at least four independent origins from parents with different plastomes. In contrast, para- and polyphyletic taxa of the octoploid Fragaria species are best explained by incomplete lineage sorting and/or hybridization. Biogeographic patterns in F. cascadensis are probably a result of range shift during the last glacial maximum.

New prospects in the detection and comparative analysis of hybridization in the tree of life

Assessing the relative importance of the various pathways to diversification is a central goal of biodiversity researchers. For plant biologists, and increasingly across the spectrum of biological sciences, among these pathways of interest is hybridization. New methodological developments are moving the field away from questions of whether natural hybridization occurs or hybrids can persist and toward more direct assessments of the long-term impact of hybridization on diversification and genome organization. Advances in theory and new data, especially phylogenomic data, have changed the face of this field, revealing extensive occurrences of hybridization at both shallow and deep levels, but lacking is a synthesis of these advancements. Here we provide an overview of methods that have been proposed for detecting hybridization with molecular data and advocate a time-extended, comparative view of reticulate evolution. In particular, we pose three overarching questions, newly placed within reach, that are critical for advancing our understanding of hybridization pattern and process: (1) How often is introgression biased toward certain genomes and loci, and is this bias selectively neutral? (2) What are the relative rates of formation of hybrid species and introgressants, and how does this compare to their subsequent fates? (3) Has the frequency of hybridization increased under historical periods of greater dynamism in climate and geographic range, such as the Pleistocene?

Prospects on the evolutionary mitogenomics of plants: A case study on the olive family (Oleaceae)

The mitogenome is rarely used to reconstruct the evolutionary history of plants, contrary to nuclear and plastid markers. Here, we evaluate the usefulness of mitochondrial DNA for molecular evolutionary studies in Oleaceae, in which cases of cytoplasmic male sterility (CMS) and of potentially contrasted organelle inheritance are known. We compare the diversity and the evolution of mitochondrial and chloroplast genomes by focusing on the olive complex and related genera. Using high-throughput techniques, we reconstructed complete mitogenomes (ca. 0.7 Mb) and plastomes (ca. 156 kb) for six olive accessions and one Chionanthus. A highly variable organization of mitogenomes was observed at the species level. In olive, two specific chimeric genes were identified in the mitogenome of lineage E3 and may be involved in CMS. Plastid-derived regions (mtpt) were observed in all reconstructed mitogenomes. Through phylogenetic reconstruction, we demonstrate that multiple integrations of mtpt regions have occurred in Oleaceae, but mtpt regions shared by all members of the olive complex derive from a common ancestor. We then assembled 52 conserved mitochondrial gene regions and complete plastomes of ten additional accessions belonging to tribes Oleeae, Fontanesieae and Forsythieae. Phylogenetic congruence between topologies based on mitochondrial regions and plastomes suggests a strong disequilibrium linkage between both organellar genomes. Finally, while phylogenetic reconstruction based on plastomes fails to resolve the evolutionary history of maternal olive lineages in the Mediterranean area, their phylogenetic relationships were successfully resolved with complete mitogenomes. Overall, our study demonstrates the great potential of using mitochondrial DNA in plant phylogeographic and metagenomic studies.

The complete chloroplast genome sequence of strawberry (Fragaria × ananassa Duch.) and comparison with related species of Rosaceae

Compared with other members of the family Rosaceae, the chloroplast genomes of Fragaria species exhibit low variation, and this situation has limited phylogenetic analyses; thus, complete chloroplast genome sequencing of Fragaria species is needed. In this study, we sequenced the complete chloroplast genome of F. × ananassa 'Benihoppe' using the Illumina HiSeq 2500-PE150 platform and then performed a combination of de novo assembly and reference-guided mapping of contigs to generate complete chloroplast genome sequences. The chloroplast genome exhibits a typical quadripartite structure with a pair of inverted repeats (IRs, 25,936 bp) separated by large (LSC, 85,531 bp) and small (SSC, 18,146 bp) single-copy (SC) regions. The length of the F. × ananassa 'Benihoppe' chloroplast genome is 155,549 bp, representing the smallest Fragaria chloroplast genome observed to date. The genome encodes 112 unique genes, comprising 78 protein-coding genes, 30 tRNA genes and four rRNA genes. Comparative analysis of the overall nucleotide sequence identity among ten complete chloroplast genomes confirmed that for both coding and non-coding regions in Rosaceae, SC regions exhibit higher sequence variation than IRs. The Ka/Ks ratio of most genes was less than 1, suggesting that most genes are under purifying selection. Moreover, the mVISTA results also showed a high degree of conservation in genome structure, gene order and gene content in Fragaria, particularly among three octoploid strawberries which were F. × ananassa 'Benihoppe', F. chiloensis (GP33) and F. virginiana (O477). However, when the sequences of the coding and non-coding regions of F. × ananassa 'Benihoppe' were compared in detail with those of F. chiloensis (GP33) and F. virginiana (O477), a number of SNPs and InDels were revealed by MEGA 7. Six non-coding regions (trnK-matK, trnS-trnG, atpF-atpH, trnC-petN, trnT-psbD and trnP-psaJ) with a percentage of variable sites greater than 1% and no less than five parsimony-informative sites were identified and may be useful for phylogenetic analysis of the genus Fragaria.

The complete chloroplast genome sequence of Aconitum coreanum and Aconitum carmichaelii and comparative analysis with other Aconitum species

Aconitum species (belonging to the Ranunculaceae) are well known herbaceous medicinal ingredients and have great economic value in Asian countries. However, there are still limited genomic resources available for Aconitum species. In this study, we sequenced the chloroplast (cp) genomes of two Aconitum species, A. coreanum and A. carmichaelii, using the MiSeq platform. The two Aconitum chloroplast genomes were 155,880 and 157,040 bp in length, respectively, and exhibited LSC and SSC regions separated by a pair of inverted repeat regions. Both cp genomes had 38% GC content and contained 131 unique functional genes including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. The gene order, content, and orientation of the two Aconitum cp genomes exhibited the general structure of angiosperms, and were similar to those of other Aconitum species. Comparison of the cp genome structure and gene order with that of other Aconitum species revealed general contraction and expansion of the inverted repeat regions and single copy boundary regions. Divergent regions were also identified. In phylogenetic analysis, Aconitum species positon among the Ranunculaceae was determined with other family cp genomes in the Ranunculales. We obtained a barcoding target sequence in a divergent region, ndhC–trnV, and successfully developed a SCAR (sequence characterized amplified region) marker for discrimination of A. coreanum. Our results provide useful genetic information and a specific barcode for discrimination of Aconitum species.

Genotyping-by-sequencing enables linkage mapping in three octoploid cultivated strawberry families

Genotyping-by-sequencing (GBS) was used to survey genome-wide single-nucleotide polymorphisms (SNPs) in three biparental strawberry (Fragaria × ananassa) populations with the goal of evaluating this technique in a species with a complex octoploid genome. GBS sequence data were aligned to the F. vesca ‘Fvb’ reference genome in order to call SNPs. Numbers of polymorphic SNPs per population ranged from 1,163 to 3,190. Linkage maps consisting of 30–65 linkage groups were produced from the SNP sets derived from each parent. The linkage groups covered 99% of the Fvb reference genome, with three to seven linkage groups from a given parent aligned to any particular chromosome. A phylogenetic analysis performed using the POLiMAPS pipeline revealed linkage groups that were most similar to ancestral species F. vesca for each chromosome. Linkage groups that were most similar to a second ancestral species, F. iinumae, were only resolved for Fvb 4. The quantity of missing data and heterogeneity in genome coverage inherent in GBS complicated the analysis, but POLiMAPS resolved F. × ananassa chromosomal regions derived from diploid ancestor F. vesca.

Comparative analysis of single nucleotide polymorphisms in the nuclear, chloroplast, and mitochondrial genomes in identification of phylogenetic association among seven melon (Cucumis melo L.) cultivars

A variety of melons are cultivated worldwide, and their specific biological properties make them an attractive model for molecular studies. This study aimed to investigate the single nucleotide polymorphisms (SNPs) from the mitochondrial, chloroplast, and nuclear genomes of seven melon accessions (Cucumis melo L.) to identify the phylogenetic relationships among melon cultivars with the Illumina HiSeq 2000 platform and bioinformatical analyses. The data showed that there were a total of 658 mitochondrial SNPs (207-295 in each), while there were 0-60 chloroplast SNPs among these seven melon cultivars, compared to the reference genome. Bioinformatical analysis showed that the mitochondrial tree topology was unable to separate the melon features, whereas the maximum parsimony/neighbor joining (MP/NJ) tree of the chloroplast SNPs could define melon features such as seed length, width, thickness, 100-seed weight, and type. SNPs of the nuclear genome were better than the mitochondrial and chloroplast SNPs in the identification of melon features. The data demonstrated the usefulness of mitochondrial, chloroplast, and nuclear SNPs in identification of phylogenetic associations among these seven melon cultivars.

Homomorphic ZW chromosomes in a wild strawberry show distinctive recombination heterogeneity but a small sex-determining region

Recombination in ancient, heteromorphic sex chromosomes is typically suppressed at the sex-determining region (SDR) and proportionally elevated in the pseudoautosomal region (PAR). However, little is known about recombination dynamics of young, homomorphic plant sex chromosomes. We examine male and female function in crosses and unrelated samples of the dioecious octoploid strawberry Fragaria chiloensis in order to map the small and recently evolved SDR controlling both traits and to examine recombination patterns on the incipient ZW chromosome. The SDR of this ZW system is located within a 280 kb window, in which the maternal recombination rate is lower than the paternal one. In contrast to the SDR, the maternal PAR recombination rate is much higher than the rates of the paternal PAR or autosomes, culminating in an elevated chromosome-wide rate. W-specific divergence is elevated within the SDR and a single polymorphism is observed in high species-wide linkage disequilibrium with sex. Selection for recombination suppression within the small SDR may be weak, but fluctuating sex ratios could favor elevated recombination in the PAR to remove deleterious mutations on the W. The recombination dynamics of this nascent sex chromosome with a modestly diverged SDR may be typical of other dioecious plants.

Congruent Deep Relationships in the Grape Family (Vitaceae) Based on Sequences of Chloroplast Genomes and Mitochondrial Genes via Genome Skimming

Vitaceae is well-known for having one of the most economically important fruits, i.e., the grape (Vitis vinifera). The deep phylogeny of the grape family was not resolved until a recent phylogenomic analysis of 417 nuclear genes from transcriptome data. However, it has been reported extensively that topologies based on nuclear and organellar genes may be incongruent due to differences in their evolutionary histories. Therefore, it is important to reconstruct a backbone phylogeny of the grape family using plastomes and mitochondrial genes. In this study,next-generation sequencing data sets of 27 species were obtained using genome skimming with total DNAs from silica-gel preserved tissue samples on an Illumina NextSeq 500 instrument [corrected]. Plastomes were assembled using the combination of de novo and reference genome (of V. vinifera) methods. Sixteen mitochondrial genes were also obtained via genome skimming using the reference genome of V. vinifera. Extensive phylogenetic analyses were performed using maximum likelihood and Bayesian methods. The topology based on either plastome data or mitochondrial genes is congruent with the one using hundreds of nuclear genes, indicating that the grape family did not exhibit significant reticulation at the deep level. The results showcase the power of genome skimming in capturing extensive phylogenetic data: especially from chloroplast and mitochondrial DNAs.

Multilocus Sex Determination Revealed in Two Populations of Gynodioecious Wild Strawberry, Fragaria vesca subsp. bracteata

Gynodioecy, the coexistence of females and hermaphrodites, occurs in 20% of angiosperm families and often enables transitions between hermaphroditism and dioecy. Clarifying mechanisms of sex determination in gynodioecious species can thus illuminate sexual system evolution. Genetic determination of gynodioecy, however, can be complex and is not fully characterized in any wild species. We used targeted sequence capture to genetically map a novel nuclear contributor to male sterility in a self-pollinated hermaphrodite of Fragaria vesca subsp. bracteata from the southern portion of its range. To understand its interaction with another identified locus and possibly additional loci, we performed crosses within and between two populations separated by 2000 km, phenotyped the progeny and sequenced candidate markers at both sex-determining loci. The newly mapped locus contains a high density of pentatricopeptide repeat genes, a class commonly involved in restoration of fertility caused by cytoplasmic male sterility. Examination of all crosses revealed three unlinked epistatically interacting loci that determine sexual phenotype and vary in frequency between populations. Fragaria vesca subsp. bracteata represents the first wild gynodioecious species with genomic evidence of both cytoplasmic and nuclear genes in sex determination. We propose a model for the interactions between these loci and new hypotheses for the evolution of sex determining chromosomes in the subdioecious and dioecious Fragaria.

A protocol for targeted enrichment of intron-containing sequence markers for recent radiations: A phylogenomic example from Heuchera (Saxifragaceae)

PREMISE OF THE STUDY

Phylogenetic inference is moving to large multilocus data sets, yet there remains uncertainty in the choice of marker and sequencing method at low taxonomic levels. To address this gap, we present a method for enriching long loci spanning intron-exon boundaries in the genus Heuchera.

METHODS

Two hundred seventy-eight loci were designed using a splice-site prediction method combining transcriptomic and genomic data. Biotinylated probes were designed for enrichment of these loci. Reference-based assembly was performed using genomic references; additionally, chloroplast and mitochondrial genomes were used as references for off-target reads. The data were aligned and subjected to coalescent and concatenated phylogenetic analyses to demonstrate support for major relationships.

RESULTS

Complete or nearly complete (>99%) sequences were assembled from essentially all loci from all taxa. Aligned introns showed a fourfold increase in divergence as opposed to exons. Concatenated analysis gave decisive support to all nodes, and support was also high and relationships mostly similar in the coalescent analysis. Organellar phylogenies were also well-supported and conflicted with the nuclear signal.

DISCUSSION

Our approach shows promise for resolving a recent radiation. Enrichment for introns is highly successful with little or no sequencing dropout at low taxonomic levels despite higher substitution and indel frequencies, and should be exploited in studies of species complexes.