Thomas Marcussen

Ancient hybridizations among the ancestral genomes of bread wheat

The allohexaploid bread wheat genome consists of three closely related subgenomes (A, B, and D), but a clear understanding of their phylogenetic history has been lacking. We used genome assemblies of bread wheat and five diploid relatives to analyze genome-wide samples of gene trees, as well as to estimate evolutionary relatedness and divergence times. We show that the A and B genomes diverged from a common ancestor ~7 million years ago and that these genomes gave rise to the D genome through homoploid hybrid speciation 1 to 2 million years later. Our findings imply that the present-day bread wheat genome is a product of multiple rounds of hybrid speciation (homoploid and polyploid) and lay the foundation for a new framework for understanding the wheat genome as a multilevel phylogenetic mosaic.

Inferring Species Networks from Gene Trees in High-Polyploid North American and Hawaiian Violets (Viola, Violaceae)

Abstract The phylogenies of allopolyploids take the shape of networks and cannot be adequately represented as bifurcating trees. Especially for high polyploids (i.e., organisms with more than six sets of nuclear chromosomes), the signatures of gene homoeolog loss, deep coalescence, and polyploidy may become confounded, with the result that gene trees may be congruent with more than one species network. Herein, we obtained the most parsimonious species network by objective comparison of competing scenarios involving polyploidization and homoeolog loss in a high-polyploid lineage of violets (Viola, Violaceae) mostly or entirely restricted to North America, Central America, or Hawaii. We amplified homoeologs of the low-copy nuclear gene, glucose-6-phosphate isomerase (GPI), by single-molecule polymerase chain reaction (PCR) and the chloroplast trnL-F region by conventional PCR for 51 species and subspecies. Topological incongruence among GPI homoeolog subclades, owing to deep coalescence and two instances of putative loss (or lack of detection) of homoeologs, were reconciled by applying the maximum tree topology for each subclade. The most parsimonious species network and the fossil-based calibration of the homoeolog tree favored monophyly of the high polyploids, which has resulted from allodecaploidization 9–14 Ma, involving sympatric ancestors from the extant Viola sections Chamaemelanium (diploid), Plagiostigma (paleotetraploid), and Viola (paleotetraploid). Although two of the high-polyploid lineages (Boreali-Americanae, Pedatae) remained decaploid, recurrent polyploidization with tetraploids of section Plagiostigma within the last 5 Ma has resulted in two 14-ploid lineages (Mexicanae, Nosphinium) and one 18-ploid lineage (Langsdorffianae). This implies a more complex phylogenetic and biogeographic origin of the Hawaiian violets (Nosphinium) than that previously inferred from rDNA data and illustrates the necessity of considering polyploidy in phylogenetic and biogeographic reconstruction.

Allozymic variation and relationships within Viola subsection Viola (Violaceae)

Allozyme markers from ten European taxa ofViola subsectionViola suggest that this group is allotetraploid, based on x = 5. All taxa had distinct multilocus phenotypes exceptV. alba subspp.alba andscotophylla, which were identical and different from subsp.dehnhardtii. Variation was consistently higher in Mediterranean populations than in North European ones. Hybridisation seems extensive but putative F1 hybrids were distinctly less fertile than the parental species. Nevertheless, increased fertility in later-generation hybrids and shared band patterns among taxa indicate an important role of hybridisation and introgression in past and present evolution within the subsection. The octoploidV. ambigua shows affinity toV. hirta (tetraploid). The octoploidV. suavis probably originated fromV. pyrenaica and other unidentified tetraploids, and high variability suggests polytopy or even polyphyly. The stoloniferous condition (seriesFlagellatae) seems to be primitive in the subsection but the reduction of stolons (seriesEflagellatae) may have originated multiple times.

From Gene Trees to a Dated Allopolyploid Network: Insights from the Angiosperm Genus Viola (Violaceae)

Allopolyploidization accounts for a significant fraction of speciation events in many eukaryotic lineages. However, existing phylogenetic and dating methods require tree-like topologies and are unable to handle the network-like phylogenetic relationships of lineages containing allopolyploids. No explicit framework has so far been established for evaluating competing network topologies, and few attempts have been made to date phylogenetic networks. We used a four-step approach to generate a dated polyploid species network for the cosmopolitan angiosperm genus Viola L. (Violaceae Batch.). The genus contains ca 600 species and both recent (neo-) and more ancient (meso-) polyploid lineages distributed over 16 sections. First, we obtained DNA sequences of three low-copy nuclear genes and one chloroplast region, from 42 species representing all 16 sections. Second, we obtained fossil-calibrated chronograms for each nuclear gene marker. Third, we determined the most parsimonious multilabeled genome tree and its corresponding network, resolved at the section (not the species) level. Reconstructing the “correct” network for a set of polyploids depends on recovering all homoeologs, i.e., all subgenomes, in these polyploids. Assuming the presence of Viola subgenome lineages that were not detected by the nuclear gene phylogenies (“ghost subgenome lineages”) significantly reduced the number of inferred polyploidization events. We identified the most parsimonious network topology from a set of five competing scenarios differing in the interpretation of homoeolog extinctions and lineage sorting, based on (i) fewest possible ghost subgenome lineages, (ii) fewest possible polyploidization events, and (iii) least possible deviation from expected ploidy as inferred from available chromosome counts of the involved polyploid taxa. Finally, we estimated the homoploid and polyploid speciation times of the most parsimonious network. Homoploid speciation times were estimated by coalescent analysis of gene tree node ages. Polyploid speciation times were estimated by comparing branch lengths and speciation rates of lineages with and without ploidy shifts. Our analyses recognize Viola as an old genus (crown age 31 Ma) whose evolutionary history has been profoundly affected by allopolyploidy. Between 16 and 21 allopolyploidizations are necessary to explain the diversification of the 16 major lineages (sections) of Viola, suggesting that allopolyploidy has accounted for a high percentage—between 67% and 88%—of the speciation events at this level. The theoretical and methodological approaches presented here for (i) constructing networks and (ii) dating speciation events within a network, have general applicability for phylogenetic studies of groups where allopolyploidization has occurred. They make explicit use of a hitherto underexplored source of ploidy information from chromosome counts to help resolve phylogenetic cases where incomplete sequence data hampers network inference. Importantly, the coalescent-based method used herein circumvents the assumption of tree-like evolution required by most techniques for dating speciation events.

Evolution, phylogeography, and taxonomy within the Viola alba complex (Violaceae)

Abstract. Taxa of the Viola alba complex were investigated using allozymes and morphometry. A taxonomic revision is presented. A wide delimitation of V. alba with only three morphological and geographical subspecies is suggested: (1) ssp. dehnhardtii distributed in the Mediterranean eastwards to Turkey; (2) ssp. alba flanking ssp. dehnhardtii in the north and east; and (3) ssp. cretica endemic to Crete. Ssp. cretica, up to now treated as a separate species, is particularly close to ssp. dehnhardtii. Viola cadevallii (NW Mediterranean) is included in the synonymy of ssp. dehnhardtii. Ssp. scotophylla (S Europe), ssp. thessala (Balkan), V. armena (Turkey), and V. besseri (Caucasus) are reduced to synonyms of V. alba ssp. alba. Viola pentelica (Greece) might represent transitional forms between ssp. alba and ssp. dehnhardtii. Glacial refugia for ssp. alba are suggested from the eastern Mediterranean via Turkey to the Caucasus, for ssp. dehnhardtii in the Mediterranean area in general, and for ssp. cretica in Crete. A key to the subspecies is provided. Taxonomic recombination: Viola alba Bess. ssp. cretica (Boiss. & Heldr.) Marcussen, comb. nov.

Marginal Likelihood Estimate Comparisons to Obtain Optimal Species Delimitations in Silene sect. Cryptoneurae (Caryophyllaceae)

Coalescent-based inference of phylogenetic relationships among species takes into account gene tree incongruence due to incomplete lineage sorting, but for such methods to make sense species have to be correctly delimited. Because alternative assignments of individuals to species result in different parametric models, model selection methods can be applied to optimise model of species classification. In a Bayesian framework, Bayes factors (BF), based on marginal likelihood estimates, can be used to test a range of possible classifications for the group under study. Here, we explore BF and the Akaike Information Criterion (AIC) to discriminate between different species classifications in the flowering plant lineage Silene sect. Cryptoneurae (Caryophyllaceae). We estimated marginal likelihoods for different species classification models via the Path Sampling (PS), Stepping Stone sampling (SS), and Harmonic Mean Estimator (HME) methods implemented in BEAST. To select among alternative species classification models a posterior simulation-based analog of the AIC through Markov chain Monte Carlo analysis (AICM) was also performed. The results are compared to outcomes from the software BP&P. Our results agree with another recent study that marginal likelihood estimates from PS and SS methods are useful for comparing different species classifications, and strongly support the recognition of the newly described species S. ertekinii.

Evolution of plant RNA polymerase IV/V genes: evidence of subneofunctionalization of duplicated NRPD2/NRPE2-like paralogs in Viola (Violaceae)

BackgroundDNA-dependent RNA polymerase IV and V (Pol IV and V) are multi-subunit enzymes occurring in plants. The origin of Pol V, specific to angiosperms, from Pol IV, which is present in all land plants, is linked to the duplication of the gene encoding the largest subunit and the subsequent subneofunctionalization of the two paralogs (NRPD1 and NRPE1). Additional duplication of the second-largest subunit, NRPD2/NRPE2, has happened independently in at least some eudicot lineages, but its paralogs are often subject to concerted evolution and gene death and little is known about their evolution nor their affinity with Pol IV and Pol V.ResultsWe sequenced a ~1500 bp NRPD2/E2-like fragment from 18 Viola species, mostly paleopolyploids, and 6 non-Viola Violaceae species. Incongruence between the NRPD2/E2-like gene phylogeny and species phylogeny indicates a first duplication of NRPD2 relatively basally in Violaceae, with subsequent sorting of paralogs in the descendants, followed by a second duplication in the common ancestor of Viola and Allexis. In Viola, the mutation pattern suggested (sub-) neofunctionalization of the two NRPD2/E2-like paralogs, NRPD2/E2-a and NRPD2/E2-b. The dN/dSratios indicated that a 54 bp region exerted strong positive selection for both paralogs immediately following duplication. This 54 bp region encodes a domain that is involved in the binding of the Nrpd2 subunit with other Pol IV/V subunits, and may be important for correct recognition of subunits specific to Pol IV and Pol V. Across all Viola taxa 73 NRPD2/E2-like sequences were obtained, of which 23 (32%) were putative pseudogenes - all occurring in polyploids. The NRPD2 duplication was conserved in all lineages except the diploid MELVIO clade, in which NRPD2/E2-b was lost, and its allopolyploid derivates from hybridization with the CHAM clade, section Viola and section Melanium, in which NRPD2/E2-a occurred in multiple copies while NRPD2/E2-b paralogs were either absent or pseudogenized.ConclusionsFollowing the relatively recent split of Pol IV and Pol V, our data indicate that these two multi-subunit enzymes are still in the process of specialization and each acquiring fully subfunctionalized copies of their subunit genes. Even after specialization, the NRPD2/E2-like paralogs are prone to pseudogenization and gene conversion and NRPD2 and NRPE2 copy number is a highly dynamic process modulated by allopolyploidy and gene death.

A Phylogeny of the Violaceae (Malpighiales) Inferred from Plastid DNA Sequences: Implications for Generic Diversity and Intrafamilial Classification

Abstract The Violaceae consist of 1,000–1,100 species of herbs, shrubs, lianas, and trees that are placed in 22 recognized genera. In this study we tested the monophyly of genera with a particular focus on the morphologically heterogeneous Rinorea and Hybanthus, the second and third most species-rich genera in the family, respectively. We also investigated intrafamilial relationships in the Violaceae with taxon sampling which included all described genera and several unnamed generic segregates. Phylogenetic inference was based on maximum parsimony, maximum likelihood, and Bayesian analyses of DNA sequences from the trnL/trnL-F and rbcL plastid regions for 102 ingroup accessions. Results from phylogenetic analyses showed Rinorea and Hybanthus to be polyphyletic, with each genus represented by three and nine clades, respectively. Results also showed that most intrafamilial taxa from previous classifications of the Violaceae were not supported. The phylogenetic inferences presented in this study illustrate the need to describe new generic segregates and to reinstate other genera, as well as to revise the traditionally accepted intrafamilial classification, which is artificial and principally based on the continuous and homoplasious character state of floral symmetry.

Establishing the phylogenetic origin, history, and age of the narrow endemic Viola guadalupensis (Violaceae)

PREMISE OF THE STUDY Climate change and shifts in land use are two major threats to biodiversity and are likely to disproportionately impact narrow endemics. Understanding their origins and the extent of their genetic diversity will enable land managers to better conserve these unique, highly localized gene pools. Viola guadalupensis is a narrow endemic of the Guadalupe Mountains (west Texas, USA). Its affinities within Viola section Chamaemelanium have been the subject of some debate. Furthermore, the polyploid and presumably reticulate relationships within this section remain largely unknown. METHODS We counted chromosomes for V. guadalupensis. Phylogenies for the chloroplast trnL-F region and the low-copy nuclear gene GPI for 24 Viola taxa were generated and used to produce a polyploid phylogenetic network. Divergence dates were obtained by fossil calibration. KEY RESULTS Meiotic chromosome counts revealed that V. guadalupensis is tetraploid (n = 12), and the presence of two GPI homoeologs further suggested allotetraploidy. Phylogenetic reconstructions showed that it originated through hybridization between unidentified members of subsection Canadenses (paternal parent) and subsection Nuttallianae (maternal parent). A fossil-calibrated relaxed clock dating analysis of GPI estimated the maximum age of V. guadalupensis to be 8.6 (5.7-11.6) Myr, suggesting the species evolved after the Guadalupe Mountains formed 12-13 Ma. CONCLUSIONS Viola guadalupensis originated by intersubsectional hybridization followed by polyploidization. Within section Chamaemelanium, this phenomenon has occurred repeatedly in the last 9 Myr (at least for V. bakeri, V. douglasii, V. glabella, and V. sempervirens). Consequences for the systematics of the section are discussed.

Assignment of Homoeologs to Parental Genomes in Allopolyploids for Species Tree Inference, with an Example from Fumaria (Papaveraceae)

There is a rising awareness that species trees are best inferred from multiple loci while taking into account processes affecting individual gene trees, such as substitution model error (failure of the model to account for the complexity of the data) and coalescent stochasticity (presence of incomplete lineage sorting [ILS]). Although most studies have been carried out in the context of dichotomous species trees, these processes operate also in more complex evolutionary histories involving multiple hybridizations and polyploidy. Recently, methods have been developed that accurately handle ILS in allopolyploids, but they are thus far restricted to networks of diploids and tetraploids. We propose a procedure that improves on this limitation by designing a workflow that assigns homoeologs to hypothetical diploid ancestral genomes prior to genome tree construction. Conflicting assignment hypotheses are evaluated against substitution model error and coalescent stochasticity. Incongruence that cannot be explained by stochastic mechanisms needs to be explained by other processes (e.g., homoploid hybridization or paralogy). The data can then be filtered to build multilabeled genome phylogenies using inference methods that can recover species trees, either in the face of substitution model error and coalescent stochasticity alone, or while simultaneously accounting for hybridization. Methods are already available for folding the resulting multilabeled genome phylogeny into a network. We apply the workflow to the reconstruction of the reticulate phylogeny of the plant genus Fumaria (Papaveraceae) with ploidal levels ranging from 2[Formula: see text] to 14[Formula: see text]. We describe the challenges in recovering nuclear NRPB2 homoeologs in high ploidy species while combining in vivo cloning and direct sequencing techniques. Using parametric bootstrapping simulations we assign nuclear homoeologs and chloroplast sequences (four concatenated loci) to their common hypothetical diploid ancestral genomes. As these assignments hinge on effective population size assumptions, we investigate how varying these assumptions impacts the recovered multilabeled genome phylogeny.