Extensive chloroplast haplotype variation indicates Pleistocene hybridization and radiation of North American Arabis drummondii, A. x divaricarpa, and A. holboellii (Brassicaceae)

A Hyb-Seq phylogeny of Boechera and related genera using a combination of Angiosperms353 and Brassicaceae-specific bait sets

PREMISE

Although Boechera (Boechereae, Brassicaceae) has become a plant model system for both ecological genomics and evolutionary biology, all previous phylogenetic studies have had limited success in resolving species relationships within the genus. The recent effective application of sequence data from target enrichment approaches to resolve the evolutionary relationships of several other challenging plant groups prompted us to investigate their usefulness in Boechera and Boechereae.

METHODS

To resolve the phylogeny of Boechera and closely related genera, we utilized the Hybpiper pipeline to analyze two combined bait sets: Angiosperms353, with broad applicability across flowering plants; and a Brassicaceae-specific bait set designed for use in the mustard family. Relationships for 101 samples representing 81 currently recognized species were inferred from a total of 1114 low-copy nuclear genes using both supermatrix and species coalescence methods.

RESULTS

Our analyses resulted in a well-resolved and highly supported phylogeny of the tribe Boechereae. Boechereae is divided into two major clades, one comprising all western North American species of Boechera, the other encompassing the eight other genera of the tribe. Our understanding of relationships within Boechera is enhanced by the recognition of three core clades that are further subdivided into robust regional species complexes.

CONCLUSIONS

This study presents the first broadly sampled, well-resolved phylogeny for most known sexual diploid Boechera. This effort provides the foundation for a new phylogenetically informed taxonomy of Boechera that is crucial for its continued use as a model system.

Identifying Climatic Drivers of Hybridization with a New Ancestral Niche Reconstruction Method

Applications of molecular phylogenetic approaches have uncovered evidence of hybridization across numerous clades of life, yet the environmental factors responsible for driving opportunities for hybridization remain obscure. Verbal models implicating geographic range shifts that brought species together during the Pleistocene have often been invoked, but quantitative tests using paleoclimatic data are needed to validate these models. Here, we produce a phylogeny for Heuchereae, a clade of 15 genera and 83 species in Saxifragaceae, with complete sampling of recognized species, using 277 nuclear loci and nearly complete chloroplast genomes. We then employ an improved framework with a coalescent simulation approach to test and confirm previous hybridization hypotheses and identify one new intergeneric hybridization event. Focusing on the North American distribution of Heuchereae, we introduce and implement a newly developed approach to reconstruct potential past distributions for ancestral lineages across all species in the clade and across a paleoclimatic record extending from the late Pliocene. Time calibration based on both nuclear and chloroplast trees recovers a mid- to late-Pleistocene date for most inferred hybridization events, a timeframe concomitant with repeated geographic range restriction into overlapping refugia. Our results indicate an important role for past episodes of climate change, and the contrasting responses of species with differing ecological strategies, in generating novel patterns of range contact among plant communities and therefore new opportunities for hybridization. The new ancestral niche method flexibly models the shape of niche while incorporating diverse sources of uncertainty and will be an important addition to the current comparative methods toolkit. [Ancestral niche reconstruction; hybridization; paleoclimate; pleistocene.].

Moving beyond assumptions: Polyploidy and environmental effects explain a geographical parthenogenesis scenario in European plants

Geographical parthenogenesis (GP) describes the phenomenon that apomicts tend to have larger distribution areas and/or occur at higher altitudes or latitudes compared to sexual relatives. However, the complex effects of genome-wide heterozygosity, ploidy, reproduction mode (sexual vs. apomictic), and environment shaping GP of plants are still not well understood. We ascertained ploidy and reproduction mode by flow cytometry of 221 populations, and added genomic RADseq data (maximum 33,165 loci) of 80 taxa of the Ranunculus auricomus polyploid plant complex in temperate Europe. We observed 7% mainly diploid sexual, 28% facultative apomictic (mean sexuality 7.1%), and 65% obligate apomictic populations. Sexuals occupied a more southern, smaller distribution area, whereas apomicts expanded their range to higher latitudes. Within the complex, we detected three main genetic clusters and highly reticulate relationships. A genetically-informed path analysis using GLMMs revealed several significant relationships. Sexuality of populations (percent of sexual seeds) was higher in diploids compared to polyploids, associated with more petals, and similar between forests and open habitats. In contrast to other apomictic plant complexes, sexuality was mainly positively correlated to solar radiation and isothermality, which fits the southern distribution. We found up to three times higher heterozygosity in polyploids compared to diploids, and generally more heterozygous individuals in forests compared with open habitats. Interestingly, we revealed a previously unknown positive association between heterozygosity and temperature seasonality, suggesting a higher resistance of polyploids to more extreme climatic conditions. We provide empirical evidence for intrinsic and extrinsic factors shaping the GP pattern in a polyploid plant complex.

The Evolution of Sex is Tempered by Costly Hybridization in Boechera (Rock Cress)

Despite decades of research, the evolution of sex remains an enigma in evolutionary biology. Typically, research addresses the costs of sex and asexuality to characterize the circumstances favoring one reproductive mode. Surprisingly few studies address the influence of common traits that are, in many organisms, obligately correlated with asexuality, including hybridization and polyploidy. These characteristics have substantial impacts on traits under selection. In particular, the fitness consequences of hybridization (i.e., reduced fitness due to interspecific reproductive isolation) will influence the evolution of sex. This may comprise a cost of either sex or asexuality due to the link between hybridity and asexuality. We examined reproductive isolation in the formation of de novo hybrid lineages between 2 widespread species in the ecological model system Boechera. Seventeen percent of 664 crosses produced F1 fruits, and only 10% of these were viable, suggesting that postmating prezygotic and postzygotic barriers inhibit hybrid success in this system. The postmating prezygotic barrier was asymmetrical, with 110 of 115 total F1 fruits produced when Boechera stricta acted as maternal parent. This asymmetry was confirmed in wild-collected lineages, using a chloroplast phylogeny of wild-collected B. stricta, Boechera retrofracta, and hybrids. We next compared fitness of F2 hybrids and selfed parental B. stricta lines, finding that F2 fitness was reduced by substantial hybrid sterility. Multiple reproductively isolating barriers influence the formation and fitness of hybrid lineages in the wild, and the costs of hybridization likely have profound impacts on the evolution of sex in the natural environment.

Comparative population genetic analyses suggest hybrid origin of R hododendron pubicostatum, an endangered plant species with extremely small populations endemic to Yunnan, China

Gene flow between sympatric congeneric plants is thought to be very common and may pose serious threats to endangered species. In the present study, we evaluate the genetic diversity and divergence of three sympatric Rhododendron species in Jiaozi Mountain using newly developed microsatellites through the Illumina MiSeq sequencing approach. Genetic diversity of all three Rhododendron species studied was moderate in comparison to genetic parameters previously reported from species of this genus. Interestingly, genetic structure analysis of the three species identified a possible hybrid origin of the threatened Rh. pubicostatum. This sympatry should be considered a unimodal hybrid zone, since R h. pubicostatum is predominant here. Unimodal hybrid zones are uncommon in Rhododendron, despite the fact that hybridization frequently occurs in the genus. Issues pertaining to the conservation of R h. pubicostatum resulting from admixture of genetic material from its parental species are discussed.

The Quaternary evolutionary history of Bristol rock cress (Arabis scabra, Brassicaceae), a Mediterranean element with an outpost in the north-western Atlantic region

BACKGROUND AND AIMS

Bristol rock cress is among the few plant species in the British Isles considered to have a Mediterranean-montane element. Spatiotemporal patterns of colonization of the British Isles since the last interglacial and after the Last Glacial Maximum (LGM) from mainland Europe are underexplored and have not yet included such floristic elements. Here we shed light on the evolutionary history of a relic and outpost metapopulation of Bristol rock cress in the south-western UK.

METHODS

Amplified fragment length polymorphisms (AFLPs) were used to identify distinct gene pools. Plastome assembly and respective phylogenetic analysis revealed the temporal context. Herbarium material was largely used to exemplify the value of collections to obtain a representative sampling covering the entire distribution range.

KEY RESULTS

The AFLPs recognized two distinct gene pools, with the Iberian Peninsula as the primary centre of genetic diversity and the origin of lineages expanding before and after the LGM towards mountain areas in France and Switzerland. No present-day lineages are older than 51 ky, which is in sharp contrast to the species stem group age of nearly 2 My, indicating severe extinction and bottlenecks throughout the Pleistocene. The British Isles were colonized after the LGM and feature high genetic diversity.

CONCLUSIONS

The short-lived perennial herb Arabis scabra, which is restricted to limestone, has expanded its distribution range after the LGM, following corridors within an open landscape, and may have reached the British Isles via the desiccated Celtic Sea at about 16 kya. This study may shed light on the origin of other rare and peculiar species co-occurring in limestone regions in the south-western British Isles.

Differential activity of F-box genes and E3 ligases distinguishes sexual versus apomictic germline specification in Boechera

Germline specification is the first step during sexual and apomictic plant reproduction, and takes place in the nucellus of the ovule, a specialized domain of the reproductive flower tissues. In each case, a sporophytic cell is determined to form the sexual megaspore mother cell (MMC) or an apomictic initial cell (AIC). These differ in their developmental fates: while the MMC undergoes meiosis, the AIC modifies or omits meiosis to form the female gametophyte. Despite great interest in these distinct developmental processes, little is known about their gene regulatory basis. To elucidate the gene regulatory networks underlying germline specification, we conducted tissue-specific transcriptional profiling using laser-assisted microdissection and RNA sequencing to compare the transcriptomes of nucellar tissues between different sexual and apomictic Boechera accessions representing four species and two ploidy levels. This allowed us to distinguish between expression differences caused by genetic background or reproductive mode. Statistical data analysis revealed 45 genes that were significantly differentially expressed, and which potentially play a role for determination of the reproductive mode. Based on annotations, these included F-box genes and E3 ligases that most likely relate to genes previously described as regulators important for germline development. Our findings provide novel insights into the transcriptional basis of sexual and apomictic reproduction.

The Boechera Genus as a Resource for Apomixis Research

The genera Boechera (A. Löve et D. Löve) and Arabidopsis, the latter containing the model plant Arabidopsis thaliana, belong to the same clade within the Brassicaceae family. Boechera is the only among the more than 370 genera in the Brassicaceae where apomixis is well documented. Apomixis refers to the asexual reproduction through seed, and a better understanding of the underlying mechanisms has great potential for applications in agriculture. The Boechera genus currently includes 110 species (of which 38 are reported to be triploid and thus apomictic), which are distributed mostly in the North America. The apomictic lineages of Boechera occur at both the diploid and triploid level and show signs of a hybridogenic origin, resulting in a modification of their chromosome structure, as reflected by alloploidy, aneuploidy, substitutions of homeologous chromosomes, and the presence of aberrant chromosomes. In this review, we discuss the advantages of the Boechera genus to study apomixis, consider its modes of reproduction as well as the inheritance and possible mechanisms controlling apomixis. We also consider population genetic aspects and a possible role of hybridization at the origin of apomixis in Boechera. The molecular tools available to study Boechera, such as transformation techniques, laser capture microdissection, analysis of transcriptomes etc. are also discussed. We survey available genome assemblies of Boechera spp. and point out the challenges to assemble the highly heterozygous genomes of apomictic species. Due to these challenges, we argue for the application of an alternative reference-free method for the comparative analysis of such genomes, provide an overview of genomic sequencing data in the genus Boechera suitable for such analysis, and provide examples of its application.

Calochortus gunnisonii furthers evidence for the complex genetic legacy of historical climate change in the southern Rocky Mountains

PREMISE OF THE STUDY

Climate cycles of the Quaternary have impacted plants at a global scale, leaving behind a complex genetic legacy. Species of the northern Rocky Mountains of North America were exposed to more uniform glacial patterns than the central and southern ranges, where synergistic relationships between temperature and precipitation caused differences in the timing and extent of glacier onset. We examined the genetic impacts of climate oscillations on Calochortus gunnisonii (Liliaceae) in the central and southern Rocky Mountains.

METHODS

Populations were sampled from disjunct mountain ranges across the basins of Wyoming and northern and central Colorado. Allelic data from nuclear microsatellites and plastid sequences (trnV-ndhC, petA-psbJ, and rpl16) were used to examine patterns of genetic structure between and among populations along the southern Rocky Mountain corridor.

KEY RESULTS

We infer considerable population structure concordant with mountain range of origin. Clustering analysis supports separate north and south genetic clusters on either side of major basins in Wyoming, suggesting that populations were maintained in two distinct refugia. Additionally, populations within the Sierra Madre Range of southern Wyoming show localized, divergent genetic signal indicative of a third potential glacial refugium. By contrast, recent genetic admixture is observed in the Laramie, Medicine Bow, and Front ranges, where population expansion from glacial refugia has likely occurred.

CONCLUSIONS

We conclude that during climate cycles of the Quaternary, C. gunnisonii experienced periods of population expansion and reduction, habitat fragmentation, isolation in three or more refugia, and admixture mirroring genetic impacts of other southern Rocky Mountains organisms.

Ecological differentiation facilitates fine-scale coexistence of sexual and asexual Boechera

PREMISE OF THE STUDY

Ecological differentiation (ED) between sexual and asexual organisms may permit the maintenance of reproductive polymorphism. Several studies of sexual/asexual ED in plants have shown that the geographic ranges of asexuals extend beyond those of sexuals, often in areas of higher latitude or elevation. But very little is known about ED at fine scales, wherein coexistence of sexuals and asexuals may be permitted by differential niche occupation.

METHODS

We used 149 populations of sexual and apomictic lineages in the genus Boechera (rock cress) collected across a portion of this mustard's vast range. We characterized reproductive mode, ploidy, and species identity or hybrid parentage of each individual, and then used a multipronged statistical approach to (1) identify ED between sexuals and asexuals; (2) investigate the impacts of two confounding factors, polyploidy and hybridization, on ED; and (3) determine the environmental variables underlying ED.

KEY RESULTS

We found that sexuals and asexuals are significantly ecologically differentiated across the landscape, despite fine-scale interdigitation of these two reproductive forms. Asexual reproduction was strongly associated with greater disturbance, reduced slope, and greater environmental variability. Although ploidy had little effect on the patterns observed, hybridization has a unique impact on the relationships between asexual reproduction and specific environmental variables.

CONCLUSIONS

Ecological differentiation along the axes of disturbance, slope, and climatic variability, as well as the effects of heterozygosity, may contribute to the maintenance of sexuality and asexuality across the landscape, ultimately impacting the establishment and spread of asexual lineages.

Establishing the cell biology of apomictic reproduction in diploid Boechera stricta (Brassicaceae)

Background and aims

In the Brassicaceae family, apomictic development is characteristic of the genus Boechera. Hybridization, polyploidy and environmental adaptation that arose during the evolution of Boechera may serve as (epi)genetic regulators of apomictic initiation in this genus. Here we focus on Boechera stricta, a predominantly diploid species that reproduces sexually. However, apomictic development in this species has been reported in several studies, indicating non-obligate sexuality.

Methods

A progressive investigation of flower development was conducted using three accessions to assess the reproductive system of B. stricta. We employed molecular and cyto-embryological identification using histochemistry, transmission electron microscopy and Nomarski and epifluorescence microscopy.

Key Results

Data from internal transcribed spacer (ITS) and chloroplast haplotype sequencing, in addition to microsatellite variation, confirmed the B. stricta genotype for all lines. Embryological data indicated irregularities in sexual reproduction manifested by heterochronic ovule development, longevity of meiocyte and dyad stages, diverse callose accumulation during meiocyte-to-gametophyte development, and the formation of triads and tetrads in several patterns. The arabinogalactan-related sugar epitope recognized by JIM13 immunolocalized to one or more megaspores. Furthermore, pollen sterility and a high frequency of seed abortion appeared to accompany reproduction of the accession ES512, along with the initiation of parthenogenesis. Data from flow cytometric screening revealed both sexual and apomictic seed formation.

Conclusion

These results imply that B. stricta is a species with an underlying ability to initiate apomixis, at least with respect to the lines examined here. The existence of apomixis in an otherwise diploid sexual B. stricta may provide the genomic building blocks for establishing highly penetrant apomictic diploids and hybrid relatives. Our findings demonstrate that apomixis per se is a variable trait upon which natural selection could act.

Admixture, evolution, and variation in reproductive isolation in the Boechera puberula clade

BACKGROUND

Hybridization is very common in plants, and the incorporation of new alleles into existing lineages (i.e. admixture) can blur species boundaries. However, admixture also has the potential to increase standing genetic variation. With new sequencing methods, we can now study admixture and reproductive isolation at a much finer scale than in the past. The genus Boechera is an extraordinary example of admixture, with over 400 hybrid derivates of varying ploidy levels. Yet, few studies have assessed admixture in this genus on a genomic scale.

RESULTS

In this study, we used Genotyping-by-Sequencing (GBS) to clarify the evolution of the Boechera puberula clade, whose six members are scattered across the western United States. We further assessed patterns of admixture and reproductive isolation within the group, including two additional species (B. stricta and B. retrofracta) that are widespread across North America. Based on 14,815 common genetic variants, we found evidence for some cases of hybridization. We find evidence of both recent and more ancient admixture, and that levels of admixture vary across species.

CONCLUSIONS

We present evidence for a monophyletic origin of the B. puberula group, and a split of B. puberula into two subspecies. Further, when inferring reproductive isolation on the basis of presence and absence of admixture, we found that the accumulation of reproductive isolation between species does not seem to occur linearly with time since divergence in this system. We discuss our results in the context of sexuality and asexuality in Boechera.

Asymmetric reproductive interference: The consequences of cross-pollination on reproductive success in sexual-apomictic populations of Potentilla puberula (Rosaceae)

Apomixis evolves from a sexual background and usually is linked to polyploidization. Pseudogamous gametophytic apomicts, which require a fertilization to initiate seed development, of various ploidy levels frequently co‐occur with their lower‐ploid sexual ancestors, but the stability of such mixed populations is affected by reproductive interferences mediated by cross‐pollination. Thereby, reproductive success of crosses depends on the difference in ploidy levels of mating partners, that is, on tolerance of deviation from the balanced ratio of maternal versus paternal genomes. Quality of pollen can further affect reproductive success in intercytotype pollinations. Cross‐fertilization, however, can be avoided by selfing which may be induced upon pollination with mixtures of self‐ and cross‐pollen (i.e., mentor effects). We tested for reproductive compatibility of naturally co‐occurring tetraploid sexuals and penta‐ to octoploid apomicts in the rosaceous species Potentilla puberula by means of controlled crosses. We estimated the role of selfing as a crossing barrier and effects of self‐ and cross‐pollen quality as well as maternal: paternal genomic ratios in the endosperm on reproductive success. Cross‐fertilization of sexuals by apomicts was not blocked by selfing, and seed set was reduced in hetero‐ compared to homoploid crosses. Thereby, seed set was negatively related to deviations from balanced parental genomic ratios in the endosperm. In contrast, seed set in the apomictic cytotypes was not reduced in hetero‐ compared to homoploid crosses. Thus, apomictic cytotypes either avoided intercytotype cross‐fertilization through selfing, tolerated intercytotype cross‐fertilizations without negative effects on reproductive success, or even benefitted from higher pollen quality in intercytotype pollinations. Our experiment provides evidence for asymmetric reproductive interference, in favor of the apomicts, with significantly reduced seed set of sexuals in cytologically mixed populations, whereas seed set in apomicts was not affected. Incompleteness of crossing barriers further indicated at least partial losses of a parental genomic endosperm balance requirement.

Pathways to polyploidy: indications of a female triploid bridge in the alpine species Ranunculus kuepferi (Ranunculaceae)

Polyploidy is one of the most important evolutionary processes in plants. In natural populations, polyploids usually emerge from unreduced gametes which either fuse with reduced ones, resulting in triploid offspring (triploid bridge), or with other unreduced gametes, resulting in tetraploid embryos. The frequencies of these two pathways, and male versus female gamete contributions, however, are largely unexplored. Ranunculus kuepferi occurs with diploid, triploid and autotetraploid cytotypes in the Alps, whereby diploids are mostly sexual, while tetraploids are facultative apomicts. To test for the occurrence of polyploidization events by triploid bridge, we investigated 551 plants of natural populations via flow cytometric seed screening. We assessed ploidy shifts in the embryo to reconstruct female versus male gamete contributions to polyploid embryo and/or endosperm formation. Seed formation via unreduced egg cells (BIII hybrids) occurred in all three cytotypes, while only in one case both gametes were unreduced. Polyploids further formed seeds with reduced, unfertilized egg cells (polyhaploids and aneuploids). Pollen was highly variable in diameter, but only pollen >27 μm was viable, whereby diploids produced higher proportions of well-developed pollen. Pollen size was not informative for the formation of unreduced pollen. These results suggest that a female triploid bridge via unreduced egg cells is the major pathway toward polyploidization in R. kuepferi, maybe as a consequence of constraints of endosperm development. Triploids resulting from unreduced male gametes were not observed, which explains the lack of obligate sexual tetraploid individuals and populations. Unreduced egg cell formation in diploids represents the first step toward apomixis.

Divergence of annual and perennial species in the Brassicaceae and the contribution of cis-acting variation at FLC orthologues

Variation in life history contributes to reproductive success in different environments. Divergence of annual and perennial angiosperm species is an extreme example that has occurred frequently. Perennials survive for several years and restrict the duration of reproduction by cycling between vegetative growth and flowering, whereas annuals live for 1 year and flower once. We used the tribe Arabideae (Brassicaceae) to study the divergence of seasonal flowering behaviour among annual and perennial species. In perennial Brassicaceae, orthologues of FLOWERING LOCUS C (FLC), a floral inhibitor in Arabidopsis thaliana, are repressed by winter cold and reactivated in spring conferring seasonal flowering patterns, whereas in annuals, they are stably repressed by cold. We isolated FLC orthologues from three annual and two perennial Arabis species and found that the duplicated structure of the A. alpina locus is not required for perenniality. The expression patterns of the genes differed between annuals and perennials, as observed among Arabidopsis species, suggesting a broad relevance of these patterns within the Brassicaceae. Also analysis of plants derived from an interspecies cross of A. alpina and annual A. montbretiana demonstrated that cis-regulatory changes in FLC orthologues contribute to their different transcriptional patterns. Sequence comparisons of FLC orthologues from annuals and perennials in the tribes Arabideae and Camelineae identified two regulatory regions in the first intron whose sequence variation correlates with divergence of the annual and perennial expression patterns. Thus, we propose that related cis-acting changes in FLC orthologues occur independently in different tribes of the Brassicaceae during life history evolution.

Pathways to polyploidy: indications of a female triploid bridge in the alpine species Ranunculus kuepferi (Ranunculaceae)

Polyploidy is one of the most important evolutionary processes in plants. In natural populations, polyploids usually emerge from unreduced gametes which either fuse with reduced ones, resulting in triploid offspring (triploid bridge), or with other unreduced gametes, resulting in tetraploid embryos. The frequencies of these two pathways, and male versus female gamete contributions, however, are largely unexplored. Ranunculus kuepferi occurs with diploid, triploid and autotetraploid cytotypes in the Alps, whereby diploids are mostly sexual, while tetraploids are facultative apomicts. To test for the occurrence of polyploidization events by triploid bridge, we investigated 551 plants of natural populations via flow cytometric seed screening. We assessed ploidy shifts in the embryo to reconstruct female versus male gamete contributions to polyploid embryo and/or endosperm formation. Seed formation via unreduced egg cells (B_III hybrids) occurred in all three cytotypes, while only in one case both gametes were unreduced. Polyploids further formed seeds with reduced, unfertilized egg cells (polyhaploids and aneuploids). Pollen was highly variable in diameter, but only pollen >27 μm was viable, whereby diploids produced higher proportions of well-developed pollen. Pollen size was not informative for the formation of unreduced pollen. These results suggest that a female triploid bridge via unreduced egg cells is the major pathway toward polyploidization in R. kuepferi, maybe as a consequence of constraints of endosperm development. Triploids resulting from unreduced male gametes were not observed, which explains the lack of obligate sexual tetraploid individuals and populations. Unreduced egg cell formation in diploids represents the first step toward apomixis.

Correlations of polyploidy and apomixis with elevation and associated environmental gradients in an alpine plant

Apomictic plants expand their geographical distributions more to higher elevations compared to their sexual progenitors. It was so far unclear whether this tendency is related to mode of reproduction itself or represents a side effect of polyploidy. Apomixis is advantageous for range expansions as no mating partners and pollinators are needed (Baker's rule). Polyploidy is thought to infer fitness advantages and a higher vigour that would enable plants to adjust better to more extreme climatic conditions. However, little is known about actual performance of plants at higher elevations. We analyzed 81 populations of Ranunculus kuepferi from the whole distribution area in the European Alps to quantify apomictic versus sexual seed formation via flow cytometric seed screening. Seed set and vegetative growth were measured as fitness parameters. All parameters were correlated to geographical distribution, elevation, temperature and precipitation. Flow cytometric seed screening revealed predominantly obligate sexuality (88.9 %) and facultative apomixis in diploid populations, while tetraploid populations are predominantly facultative (65.4 %) to obligate apomictic. Apomictic seed formation correlated significantly to higher elevations, which explains also the observed niche shift to lower temperatures. However, within the tetraploid range, there is no apparent correlation of degree of facultative apomixis to geographical distance. Apomixis appeared in diploids three times independently in separated, otherwise sexual populations in the southwestern refugial areas of the Alps. Diploid apomixis was not successful in range expansions, and obligate sexual polyploids were not observed. Polyploidy may relate to cold tolerance as an adaptation to conditions at high elevations, where diploid sexuals have no fitness advantage. Instead, facultative apomixis may have aided colonization of higher elevations and range expansions in the Alps without mate and pollinator limitation, but did not necessarily involve long-distance dispersal. A direct influence of low temperatures on unreduced gamete formation cannot be ruled out.

Molecular phylogeny, divergence time estimates and historical biogeography within one of the world's largest monocot genera

A primary aim of historical biogeography is to identify the causal factors or processes that have shaped the composition and distribution of biotas over time. Another is to infer the evolution of geographic ranges of species and clades in a phylogenetic context. To this end, historical biogeography addresses important questions such as: Where were ancestors distributed? Where did lineages originate? Which processes cause geographic ranges to evolve through time? Allium subgenus Anguinum comprises approximately twelve taxa with a disjunct distribution in the high mountains from south-western Europe to eastern Asia and in northeastern North America. Although both the systematic position and the geographical limits of Anguinum have been identified, to date no molecular systematic study has been performed utilizing a comprehensive sampling of these species. With an emphasis on the Anguinum eastern Asian geographical group, the goals of the present study were: (i) to infer species-level phylogenetic relationships within Anguinum, (ii) to assess molecular divergence and estimated the times of the major splits in Anguinum and (iii) to trace the biogeographic history of the subgenus. Four DNA sequences (ITS, matK, trnH-psbA, rps16) were used to reconstruct the phylogeny of Allium subgen. Anguinum RbcL sequences were used to estimate divergences time for Allium, and sequences of ITS were used to estimate the divergence times for Anguinum and its main lineages and to provide implications for the evolutionary history of the subgenus. Phylogenetic analyses for all Allium corroborate that Anguinum is monophyletic and indicate that Anguinum is composed of two sister groups: one with a Eurasian-American distribution, and the other restricted to eastern Asia. In the eastern Asian geographical group, incongruence between gene trees and morphology-based taxonomies was recovered as was incongruence between data from plastid and nuclear sequences. This incongruence is likely due to the combined effects of a recent radiation, incomplete lineage sorting, and hybridization/introgression. Divergence time estimates suggest that the crown group of Anguinum originated during the late Miocene (ca. 7.16 Mya) and then diverged and dispersed. Biogeographic analyses using statistical dispersal-vicariance analysis (S-DIVA) and a likelihood method support an eastern Asia origin of Anguinum It is inferred that in the late Pliocene/Early Pleistocene, with cooling climates and the uplift of the Himalayas and Hengduan Mountains, the ancestor of the eastern Asian alliance clade underwent a very recent radiation.

Phylogeography and modes of reproduction in diploid and tetraploid halophytes of Limonium species (Plumbaginaceae): evidence for a pattern of geographical parthenogenesis

BACKGROUND AND AIMS

The genus Limonium (Plumbaginaceae) has long been recognized to have sexual and apomictic (asexual seed formation) modes of reproduction. This study aimed to elucidate phylogeographical patterns and modes of reproduction in diploid and tetraploid Limonium species, namely three putative sexual diploid species with morphological affinities (L. nydeggeri, L. ovalifolium, L. lanceolatum) and three related, probably apomict tetraploid species (L. binervosum, L. dodartii, L. multiflorum).

METHODS

cpDNA diversity and differentiation between natural populations of the species were investigated using two chloroplast sequence regions (trnL intron and trnL-trnF intergenic spacer). Floral heteromorphies, ovule cytoembryological analyses and pollination and crossing tests were performed in representative species of each ploidy group, namely diploid L. ovalifolium and tetraploid L. multiflorum, using plants from greenhouse collections.

KEY RESULTS AND CONCLUSIONS

Genetic analyses showed that diploid species have a higher haplotype diversity and a higher number of unique (endemic) haplotypes than tetraploid species. Network analysis revealed correlations between cpDNA haplotype distribution and ploidy groups, species groups and geographical origin, and haplotype sharing within and among species with distinct ploidy levels. Reproductive biology analyses showed that diploid L. ovalifolium mainly forms meiotically reduced tetrasporic embryo sacs of Gagea ova, Adoxa and Drusa types. Limonium multiflorum, however, has only unreduced, diplosporic (apomictic) embryo sacs of Rudbeckia type, and autonomous apomictic development seems to occur. Taken together, the findings provide evidence of a pattern of 'geographical parthenogenesis' in which quaternary climatic oscillations appear to be involved in the geographical patterns of coastal diploid and tetraploid Limonium species.

Beringian sub-refugia revealed in blackfish (Dallia): implications for understanding the effects of Pleistocene glaciations on Beringian taxa and other Arctic aquatic fauna

Background

Pleistocene climatic instability had profound and diverse effects on the distribution and abundance of Arctic organisms revealed by variation in phylogeographic patterns documented in extant Arctic populations. To better understand the effects of geography and paleoclimate on Beringian freshwater fishes, we examined genetic variability in the genus Dallia (blackfish: Esociformes: Esocidae). The genus Dallia groups between one and three nominal species of small, cold- and hypoxia-tolerant freshwater fishes restricted entirely in distribution to Beringia from the Yukon River basin near Fairbanks, Alaska westward including the Kuskokwim River basin and low-lying areas of Western Alaska to the Amguema River on the north side of the Chukotka Peninsula and Mechigmen Bay on the south side of the Chukotka Peninsula. The genus has a non-continuous distribution divided by the Bering Strait and the Brooks Range. We examined the distribution of genetic variation across this range to determine the number and location of potential sub-refugia within the greater Beringian refugium as well as the roles of the Bering land bridge, Brooks Range, and large rivers within Beringia in shaping the current distribution of populations of Dallia. Our analyses were based on DNA sequence data from two nuclear gene introns (S7 and RAG1) and two mitochondrial genome fragments from nineteen sampling locations. These data were examined under genetic clustering and coalescent frameworks to identify sub-refugia within the greater Beringia refugium and to infer the demographic history of different populations of Dallia.

Results

We identified up to five distinct genetic clusters of Dallia. Four of these genetic clusters are present in Alaska: (1) Arctic Coastal Plain genetic cluster found north of the Brooks Range, (2) interior Alaska genetic cluster placed in upstream locations in the Kuskokwim and Yukon river basins, (3) a genetic cluster found on the Seward Peninsula, and (4) a coastal Alaska genetic cluster encompassing downstream Kuskokwim River and Yukon River basin sample locations and samples from Southwest Alaska not in either of these drainages. The Chukotka samples are assigned to their own genetic cluster (5) similar to the coastal Alaska genetic cluster. The clustering and ordination analyses implemented in Discriminant Analysis of Principal Components (DAPC) and STRUCTURE showed mostly concordant groupings and a high degree of differentiation among groups. The groups of sampling locations identified as genetic clusters correspond to geographic areas divided by likely biogeographic barriers including the Brooks Range and the Bering Strait. Estimates of sequence diversity (θ) are highest in the Yukon River and Kuskokwim River drainages near the Bering Sea. We also infer asymmetric migration rates between genetic clusters. The isolation of Dallia on the Arctic Coastal Plain of Alaska is associated with very low estimated migration rates between the coastal Alaska genetic cluster and the Arctic Coastal Plain genetic cluster.

Conclusions

Our results support a scenario with multiple aquatic sub-refugia in Beringia during the Pleistocene and the preservation of that structure in extant populations of Dallia. An inferred historical presence of Dallia across the Bering land bridge explains the similarities in the genetic composition of Dallia in West Beringia and western coastal Alaska. In contrast, historic and contemporary isolation across the Brooks Range shaped the distinctiveness of present day Arctic Coastal Plain Dallia. Overall this study uncovered a high degree of genetic structuring among populations of Dallia supporting the idea of multiple Beringian sub-refugia during the Pleistocene and which appears to be maintained to the present due to the strictly freshwater nature and low dispersal ability of this genus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0413-2) contains supplementary material, which is available to authorized users.

Hybrid apomicts trapped in the ecological niches of their sexual ancestors

Asexual reproduction is expected to reduce the adaptive potential to novel or changing environmental conditions, restricting or altering the ecological niche of asexual lineages. Asexual lineages of plants and animals are typically polyploid, an attribute that may influence their genetic variation, plasticity, adaptive potential, and niche breadth. The genus Boechera (Brassicaceae) represents an ideal model to test the relative ecological and biogeographic impacts of reproductive mode and ploidy because it is composed of diploid sexual and both diploid and polyploid asexual (i.e., apomictic) lineages. Here, we demonstrate a strong association between a transcriptionally conserved allele and apomictic seed formation. We then use this allele as a proxy apomixis marker in 1,649 accessions to demonstrate that apomixis is likely to be a common feature across the Boechera phylogeny. Phylogeographic analyses of these data demonstrate (i) species-specific niche differentiation in sexuals, (ii) extensive niche conservation between differing reproductive modes of the same species, (iii) ploidy-specific niche differentiation within and among species, and (iv) occasional niche drift between apomicts and their sexual ancestors. We conclude that ploidy is a substantially stronger and more common driver of niche divergence within and across Boechera species although variation in both traits may not necessarily lead to niche evolution on the species scale.

Taming the wild: resolving the gene pools of non-model Arabidopsis lineages

BACKGROUND

Wild relatives in the genus Arabidopsis are recognized as useful model systems to study traits and evolutionary processes in outcrossing species, which are often difficult or even impossible to investigate in the selfing and annual Arabidopsis thaliana. However, Arabidopsis as a genus is littered with sub-species and ecotypes which make realizing the potential of these non-model Arabidopsis lineages problematic. There are relatively few evolutionary studies which comprehensively characterize the gene pools across all of the Arabidopsis supra-groups and hypothesized evolutionary lineages and none include sampling at a world-wide scale. Here we explore the gene pools of these various taxa using various molecular markers and cytological analyses.

RESULTS

Based on ITS, microsatellite, chloroplast and nuclear DNA content data we demonstrate the presence of three major evolutionary groups broadly characterized as A. lyrata group, A. halleri group and A. arenosa group. All are composed of further species and sub-species forming larger aggregates. Depending on the resolution of the marker, a few closely related taxa such as A. pedemontana, A. cebennensis and A. croatica are also clearly distinct evolutionary lineages. ITS sequences and a population-based screen based on microsatellites were highly concordant. The major gene pools identified by ITS sequences were also significantly differentiated by their homoploid nuclear DNA content estimated by flow cytometry. The chloroplast genome provided less resolution than the nuclear data, and it remains unclear whether the extensive haplotype sharing apparent between taxa results from gene flow or incomplete lineage sorting in this relatively young group of species with Pleistocene origins.

CONCLUSIONS

Our study provides a comprehensive overview of the genetic variation within and among the various taxa of the genus Arabidopsis. The resolved gene pools and evolutionary lineages will set the framework for future comparative studies on genetic diversity. Extensive population-based phylogeographic studies will also be required, however, in particular for A. arenosa and their affiliated taxa and cytotypes.

Prospects for discriminating Zingiberaceae species in India using DNA barcodes

We evaluated nine plastid (matK, rbcL, rpoC1, rpoB, rpl36-rps8, ndhJ, trnL-F, trnH-psbA, accD) and two nuclear (ITS and ITS2) barcode loci in family Zingiberaceae by analyzing 60 accessions of 20 species belonging to seven genera from India. Bidirectional sequences were recovered for every plastid locus by direct sequencing of polymerase chain reaction (PCR) amplicons in all the accessions tested. However, only 35 (58%) and 40 accessions (66%) yielded ITS and ITS2 sequences, respectively, by direct sequencing. In different bioinformatics analyses, matK and rbcL consistently resolved 15 species (75%) into monophyletic groups and five species into two paraphyletic groups. The 173 ITS sequences, including 138 cloned sequences from 23 accessions, discriminated only 12 species (60%), and the remaining species were entered into three paraphyletic groups. Phylogenetic and genealogic analyses of plastid and ITS sequences imply the possible occurrence of natural hybridizations in the evolutionary past in giving rise to species paraphyly and intragenomic ITS heterogeneity in the species tested. The results support using matK and rbcL loci for barcoding Zingiberaceae members and highlight the poor utility of ITS and the highly regarded ITS2 in barcoding this family, and also caution against proposing ITS loci for barcoding taxa based on limited sampling.

ITS polymorphisms shed light on hybrid evolution in apomictic plants: a case study on the Ranunculus auricomus complex

The reconstruction of reticulate evolutionary histories in plants is still a major methodological challenge. Sequences of the ITS nrDNA are a popular marker to analyze hybrid relationships, but variation of this multicopy spacer region is affected by concerted evolution, high intraindividual polymorphism, and shifts in mode of reproduction. The relevance of changes in secondary structure is still under dispute. We aim to shed light on the extent of polymorphism within and between sexual species and their putative natural as well as synthetic hybrid derivatives in the Ranunculus auricomus complex to test morphology-based hypotheses of hybrid origin and parentage of taxa. We employed direct sequencing of ITS nrDNA from 68 individuals representing three sexuals, their synthetic hybrids and one sympatric natural apomict, as well as cloning of ITS copies in four representative individuals, RNA secondary structure analysis, and landmark geometric morphometric analysis on leaves. Phylogenetic network analyses indicate additivity of parental ITS variants in both synthetic and natural hybrids. The triploid synthetic hybrids are genetically much closer to their maternal progenitors, probably due to ploidy dosage effects, although exhibiting a paternal-like leaf morphology. The natural hybrids are genetically and morphologically closer to the putative paternal progenitor species. Secondary structures of ITS1-5.8S-ITS2 were rather conserved in all taxa. The observed similarities in ITS polymorphisms suggest that the natural apomict R. variabilis is an ancient hybrid of the diploid sexual species R. notabilis and the sexual species R. cassubicifolius. The additivity pattern shared by R. variabilis and the synthetic hybrids supports an evolutionary and biogeographical scenario that R. variabilis originated from ancient hybridization. Concerted evolution of ITS copies in R. variabilis is incomplete, probably due to a shift to asexual reproduction. Under the condition of comprehensive inter- and intraspecific sampling, ITS polymorphisms are powerful for elucidating reticulate evolutionary histories.

On the origin and evolution of apomixis in Boechera

The genetic mechanisms causing seed development by gametophytic apomixis in plants are predominantly unknown. As apomixis is consistently associated with hybridity and polyploidy, these confounding factors may either (a) be the underlying mechanism for the expression of apomixis, or (b) obscure the genetic factors which cause apomixis. To distinguish between these hypotheses, we analyzed the population genetic patterns of diploid and triploid apomictic lineages and their sexual progenitors in the genus Boechera (Brassicaceae). We find that while triploid apomixis is associated with hybridization, the majority of diploid apomictic lineages are likely the product of intra-specific crosses. We then show that these diploid apomicts are more likely to sire triploid apomictic lineages than conspecific sexuals. Combined with flow cytometric seed screen phenotyping for male and female components of apomixis, our analyses demonstrate that hybridization is an indirect correlate of apomixis in Boechera.

A world-wide perspective on crucifer speciation and evolution: phylogenetics, biogeography and trait evolution in tribe Arabideae

BACKGROUND AND AIMS

Tribe Arabideae are the most species-rich monophyletic lineage in Brassicaceae. More than 500 species are distributed in the majority of mountain and alpine regions worldwide. This study provides the first comprehensive phylogenetic analysis for the species assemblage and tests for association of trait and characters, providing the first explanations for the enormous species radiation since the mid Miocene.

METHODS

Phylogenetic analyses of DNA sequence variation of nuclear encoded loci and plastid DNA are used to unravel a reliable phylogenetic tree. Trait and ancestral area reconstructions were performed and lineage-specific diversification rates were calculated to explain various radiations in the last 15 Myr in space and time.

KEY RESULTS

A well-resolved phylogenetic tree demonstrates the paraphyly of the genus Arabis and a new systematic concept is established. Initially, multiple radiations involved a split between lowland annuals and mountain/alpine perennial sister species. Subsequently, increased speciation rates occur in the perennial lineages. The centre of origin of tribe Arabideae is most likely the Irano-Turanian region from which the various clades colonized the temperate mountain and alpine regions of the world.

CONCLUSIONS

Mid Miocene early diversification started with increased speciation rates due to the emergence of various annual lineages. Subsequent radiations were mostly driven by diversification within perennial species during the Pliocene, but increased speciation rates also occurred during that epoch. Taxonomic concepts in Arabis are still in need of a major taxonomic revision to define monophyletic groups.

A continental-wide perspective: the genepool of nuclear encoded ribosomal DNA and single-copy gene sequences in North American Boechera (Brassicaceae)

74 of the currently accepted 111 taxa of the North American genus Boechera (Brassicaceae) were subject to pyhlogenetic reconstruction and network analysis. The dataset comprised 911 accessions for which ITS sequences were analyzed. Phylogenetic analyses yielded largely unresolved trees. Together with the network analysis confirming this result this can be interpreted as an indication for multiple, independent, and rapid diversification events. Network analyses were superimposed with datasets describing i) geographical distribution, ii) taxonomy, iii) reproductive mode, and iv) distribution history based on phylogeographic evidence. Our results provide first direct evidence for enormous reticulate evolution in the entire genus and give further insights into the evolutionary history of this complex genus on a continental scale. In addition two novel single-copy gene markers, orthologues of the Arabidopsis thaliana genes At2g25920 and At3g18900, were analyzed for subsets of taxa and confirmed the findings obtained through the ITS data.

Systematics and evolution of Arctic-Alpine Arabis alpina (Brassicaceae) and its closest relatives in the eastern Mediterranean

PREMISE OF THE STUDY

The high mountains in southern Anatolia and the eastern Mediterranean are assumed to play a major role as a primary center of genetic diversity and species richness in Eurasia. We tested this hypothesis by focusing on the widespread perennial arctic-alpine Arabis alpina and its sympatrically distributed closest relatives in the eastern Mediterranean.

METHODS

Plastid (trnL intron, trnL-F intergenic spacer) and nuclear (ITS) DNA sequence analysis was used for phylogenetic reconstruction. Broad-scale plastid haplotype analyses were conducted to infer ancestral biogeographic patterns.

KEY RESULTS

Five Arabis species, identified from the eastern Mediterranean (Turkey mainland and Cyprus), evolved directly and independently from A. alpina, leaving Arabis alpina as a paraphyletic taxon. These species are not affected by hybridization or introgression, and species divergence took place at the diploid level during the Pleistocene.

CONCLUSIONS

Pleistocene climate fluctuations produced local altitudinal range-shifts among mountain glacial survival areas, resulting not only in the accumulation of intraspecific genotype diversity but also in the formation of five local species. We also show that the closest sister group of Arabis alpina consists exclusively of annuals/winter annuals and diverged prior to Pleistocene climatic fluctuations during the colonization of the lowland Mediterranean landscape. These findings highlight that Anatolia is not only a center of species richness but also a center for life-history diversification.

Boechera, a model system for ecological genomics

The selection and development of a study system for evolutionary and ecological functional genomics (EEFG) depend on a variety of factors. Here, we present the genus Boechera as an exemplary system with which to address ecological and evolutionary questions. Our focus on Boechera is based on several characteristics as follows: (i) native populations in undisturbed habitats where current environments reflect historical conditions over several thousand years; (ii) functional genomics benefitting from its close relationship to Arabidopsis thaliana; (iii) inbreeding tolerance enabling development of recombinant inbred lines, near-isogenic lines and positional cloning; (iv) interspecific crosses permitting mapping for genetic analysis of speciation; (v) apomixis (asexual reproduction by seeds) in a genetically tractable diploid; and (vi) broad geographic distribution in North America, permitting ecological genetics for a large research community. These characteristics, along with the current sequencing of three Boechera species by the Joint Genome Institute, position Boechera as a rapidly advancing system for EEFG studies.

Testing species delimitation in sympatric species complexes: the case of an African tropical tree, Carapa spp. (Meliaceae)

Plant species delimitation within tropical ecosystems is often difficult because of the lack of diagnostic morphological characters that are clearly visible. The development of an integrated approach, which utilizes several different types of markers (both morphological and molecular), would be extremely useful in this context. Here we have addressed species delimitation of sympatric tropical tree species that belong to Carapa spp. (Meliaceae) in Central Africa. We adopted a population genetics approach, sampling numerous individuals from three locations where sympatric Carapa species are known to exist. Comparisons between morphological markers (the presence or absence of characters, leaf-shape traits) and molecular markers (chloroplast sequences, ribosomal internal transcribed spacer region (ITS) sequences, and nuclear microsatellites) demonstrated the following: (i) a strong correlation between morphological and nuclear markers; (ii) despite substantial polymorphism, the inability of chloroplast DNA to discriminate between species, suggesting that cytoplasmic markers represent ineffective DNA barcodes; (iii) lineage sorting effects when using ITS sequences; and (iv) a complex evolutionary history within the genus Carapa, which includes frequent inter-specific gene flow. Our results support the use of a population genetics approach, based on ultra-polymorphic markers, to address species delimitation within complex taxonomic groups.

CpDNA-based species identification and phylogeography: application to African tropical tree species

Despite the importance of the African tropical rainforests as a hotspot of biodiversity, their history and the processes that have structured their biodiversity are understood poorly. With respect to past demographic processes, new insights can be gained through characterizing the distribution of genetic diversity. However, few studies of this type have been conducted in Central Africa, where the identification of species in the field can be difficult. We examine here the distribution of chloroplast DNA (cpDNA) diversity in Lower Guinea in two tree species that are difficult to distinguish, Erythrophleum ivorense and Erythrophleum suaveolens (Fabaceae). By using a blind-sampling approach and comparing molecular and morphological markers, we first identified retrospectively all sampled individuals and determined the limits of the distribution of each species. We then performed a phylogeographic study using the same genetic data set. The two species displayed essentially parapatric distributions that were correlated well with the rainfall gradient, which indicated different ecological requirements. In addition, a phylogeographic structure was found for E. suaveolens and, for both species, substantially higher levels of diversity and allelic endemism were observed in the south (Gabon) than in the north (Cameroon) of the Lower Guinea region. This finding indicated different histories of population demographics for the two species, which might reflect different responses to Quaternary climate changes. We suggest that a recent period of forest perturbation, which might have been caused by humans, favoured the spread of these two species and that their poor recruitment at present results from natural succession in their forest formations.

Quantitative variation for apomictic reproduction in the genus Boechera (Brassicaceae)

PREMISE OF THE STUDY

The evolution of asexual seed production (apomixis) from sexual relatives is a great enigma of plant biology. The genus Boechera is ideal for studying apomixis because of its close relation to Arabidopsis and the occurrence of sexual and apomictic species at low ploidy levels (diploid and triploid). Apomixis is characterized by three components: unreduced embryo-sac formation (apomeiosis), fertilization-independent embryogenesis (parthenogenesis), and functional endosperm formation (pseudogamy or autonomous endosperm formation). Understanding the variation in these traits within and between species has been hindered by the laborious histological analyses required to analyze large numbers of samples. •

METHODS

To quantify variability for the different components of apomictic seed development, we developed a high-throughput flow cytometric seed screen technique to measure embryo:endosperm ploidy in over 22000 single seeds derived from 71 accessions of diploid and triploid Boechera. •

KEY RESULTS

Three interrelated features were identified within and among Boechera species: (1) variation for most traits associated with apomictic seed formation, (2) three levels of apomeiosis expression (low, high, obligate), and (3) correlations between apomeiosis and parthenogenesis/pseudogamy. •

CONCLUSIONS

The data presented here provide a framework for choosing specific genotypes for correlations with large "omics" data sets being collected for Boechera to study population structure, gene flow, and evolution of specific traits. We hypothesize that low levels of apomeiosis represent an ancestral condition of Boechera, whereas high apomeiosis levels may have been induced by global gene regulatory changes associated with hybridization.

Geographic and phylogenetic patterns in Silene section Melandrium (Caryophyllaceae) as inferred from chloroplast and nuclear DNA sequences

The phylogenetic relationships between the five dioecious species in Silene section Melandrium (Caryophyllaceae) and their putative hermaphrodite relatives are investigated based on an extensive geographic and taxonomic sample, using DNA sequence data from the chloroplast genome and the nuclear ribosomal ITS region. The hermaphrodite S. noctiflora (the type species of section Elisanthe) is distantly related to the dioecious species. With the exception of chloroplast sequences in one S. latifolia population from Turkey, the dioecious taxa form a strongly supported monophyletic group (Silene section Melandrium). The phylogenetic structure within section Melandrium differs between chloroplast and nuclear sequences. While there is extensive sharing of chloroplast haplotypes among all the dioecious species (the observed patterns reflect geographic structure), the nuclear ITS phylogeny shows a higher degree of taxonomic structure. Chloroplast-sharing by the section Melandrium species is most plausibly explained by a history of hybridization and extensive backcrossing.

Review of the Application of Modern Cytogenetic Methods (FISH/GISH) to the Study of Reticulation (Polyploidy/Hybridisation)

The convergence of distinct lineages upon interspecific hybridisation, including when accompanied by increases in ploidy (allopolyploidy), is a driving force in the origin of many plant species. In plant breeding too, both interspecific hybridisation and allopolyploidy are important because they facilitate introgression of alien DNA into breeding lines enabling the introduction of novel characters. Here we review how fluorescence in situ hybridisation (FISH) and genomic in situ hybridisation (GISH) have been applied to: 1) studies of interspecific hybridisation and polyploidy in nature, 2) analyses of phylogenetic relationships between species, 3) genetic mapping and 4) analysis of plant breeding materials. We also review how FISH is poised to take advantage of nextgeneration sequencing (NGS) technologies, helping the rapid characterisation of the repetitive fractions of a genome in natural populations and agricultural plants.

Review of the Application of Modern Cytogenetic Methods (FISH/GISH) to the Study of Reticulation (Polyploidy/Hybridisation)

The convergence of distinct lineages upon interspecific hybridisation, including when accompanied by increases in ploidy (allopolyploidy), is a driving force in the origin of many plant species. In plant breeding too, both interspecific hybridisation and allopolyploidy are important because they facilitate introgression of alien DNA into breeding lines enabling the introduction of novel characters. Here we review how fluorescence in situ hybridisation (FISH) and genomic in situ hybridisation (GISH) have been applied to: 1) studies of interspecific hybridisation and polyploidy in nature, 2) analyses of phylogenetic relationships between species, 3) genetic mapping and 4) analysis of plant breeding materials. We also review how FISH is poised to take advantage of nextgeneration sequencing (NGS) technologies, helping the rapid characterisation of the repetitive fractions of a genome in natural populations and agricultural plants.

Colonizing the American continent: Systematics of the genus Arabis in North America (Brassicaceae)

The circumscription of the genus Arabis underwent many and drastic changes within the past. Using DNA sequence information from the nuclear ribosomal RNA and parts of the plastid genome (trnL-trnLF), as well as a critical evaluation of herbarium material from East Asia and North America, we circumscribe the various Arabis taxa of North America. The American and East Asian Arabis species are closely related and, contrary to what was previously believed, they are not closely related to the Eurasian A. hirsuta. Using cpDNA, we found five North American lineages of Arabis with distinct distribution patterns, of which only the purple/red-flowered lineage consists of proven diploids that evolved directly from East Asian progenitors. All other four lineages evolved via ancient hybridization either on the Asian continent prior to migration to North America or showed significant evidence for hybridization and reticulation while diversifying on the American continent. We also provide the first evidence for the systematic circumscription of East Asian Arabis taxa, which together with the North American taxa, form one clade distantly related to European A. ciliata and Eurasian A. hirsuta. The findings also represent the first record of A. pycnocarpa for the floras of China, Japan, and Russian Far East.

Genetic structure and phylogeography of a temperate-boreal herb, Cardamine scutata (Brassicaceae), in northeast Asia inferred from AFLPs and cpDNA haplotypes

PREMISE OF THE STUDY

Studies on genetic structure of plant populations help us understand the history of local flora and vegetation. In this study, we focus on the temperate-boreal herb Cardamine scutata from northeast Asia, an area with scarce phylogeographic studies. We explore patterns of genetic variation within this species, with an aim to infer its (post-) glacial history with reference to colonization routes and migrations via land bridges. •

METHODS

We analyzed 46 populations sampled in Japan, Kamchatka, and Korea using AFLP and cpDNA sequence data. •

KEY RESULTS

Two intraspecific genetic groups were resolved, distributed in the northeastern and southwestern part of the study area, most likely reflecting lineages isolated from each other during (at least) the last glaciation. A zone of secondary contacts was found in central/northern Honshu, and a few cases of long-distance dispersal were observed. We detected efficient gene flow across the marine straits, supporting the role of land bridges created by sea level decline during the last glacial period. The cpDNA data indicated extensive recent expansion and diversification within both lineages. We inferred recent colonization of Kamchatka from Hokkaido, associated with genetic impoverishment. •

CONCLUSIONS

The pattern of north-south genetic differentiation found in C. scutata is rather common among several other plant species studied in Japan, despite their distinct biological features. We assume that different processes and factors may have brought about this similarity. Overall, this study contributes to better understanding of the biogeography of northeast Asia.

The evolutionary history of the Arabidopsis lyrata complex: a hybrid in the amphi-Beringian area closes a large distribution gap and builds up a genetic barrier

BACKGROUND

The genomes of higher plants are, on the majority, polyploid, and hybridisation is more frequent in plants than in animals. Both polyploidisation and hybridisation contribute to increased variability within species, and may transfer adaptations between species in a changing environment. Studying these aspects of evolution within a diversified species complex could help to clarify overall spatial and temporal patterns of plant speciation. The Arabidopsis lyrata complex, which is closely related to the model plant Arabidopsis thaliana, is a perennial, outcrossing, herbaceous species complex with a circumpolar distribution in the Northern Hemisphere as well as a disjunct Central European distribution in relictual habitats. This species complex comprises three species and four subspecies, mainly diploids but also several tetraploids, including one natural hybrid. The complex is ecologically, but not fully geographically, separated from members of the closely related species complex of Arabidopsis halleri, and the evolutionary histories of both species compexes have largely been influenced by Pleistocene climate oscillations.

RESULTS

Using DNA sequence data from the nuclear encoded cytosolic phosphoglucoisomerase and Internal Transcribed Spacers 1 and 2 of the ribosomal DNA, as well as the trnL/F region from the chloroplast genome, we unravelled the phylogeography of the various taxonomic units of the A. lyrata complex. We demonstrate the existence of two major gene pools in Central Europe and Northern America. These two major gene pools are constructed from different taxonomic units. We also confirmed that A. kamchatica is the allotetraploid hybrid between A. lyrata and A. halleri, occupying the amphi-Beringian area in Eastern Asia and Northern America. This species closes the large distribution gap of the various other A. lyrata segregates. Furthermore, we revealed a threefold independent allopolyploid origin of this hybrid species in Japan, China, and Kamchatka.

CONCLUSIONS

Unglaciated parts of the Eastern Austrian Alps and arctic Eurasia, including Beringia, served as major glacial refugia of the Eurasian A. lyrata lineage, whereas A. halleri and its various subspecies probably survived in refuges in Central Europe and Eastern Asia with a large distribution gap in between. The North American A. lyrata lineage probably survived the glaciation in the southeast of North America. The dramatic climatic changes during glaciation and deglaciation cycles promoted not only secondary contact and formation of the allopolyploid hybrid A. kamchatica, but also provided the environment that allowed this species to fill a large geographic gap separating the two genetically different A. lyrata lineages from Eurasia and North America. With our example focusing on the evolutionary history of the A. lyrata species complex, we add substantial information to a broad evolutionary framework for future investigations within this emerging model system in molecular and evolutionary biology.

Genome wide analyses reveal little evidence for adaptive evolution in many plant species

The relative contribution of advantageous and neutral mutations to the evolutionary process is a central problem in evolutionary biology. Current estimates suggest that whereas Drosophila, mice, and bacteria have undergone extensive adaptive evolution, hominids show little or no evidence of adaptive evolution in protein-coding sequences. This may be a consequence of differences in effective population size. To study the matter further, we have investigated whether plants show evidence of adaptive evolution using an extension of the McDonald-Kreitman test that explicitly models slightly deleterious mutations by estimating the distribution of fitness effects of new mutations. We apply this method to data from nine pairs of species. Altogether more than 2,400 loci with an average length of approximately 280 nucleotides were analyzed. We observe very similar results in all species; we find little evidence of adaptive amino acid substitution in any comparison except sunflowers. This may be because many plant species have modest effective population sizes.

Molecular phylogeny and systematics of the genus Draba (Brassicaceae) and identification of its most closely related genera

Phylogenetic analyses of sequence data from the internal transcribed spacer regions of the nuclear ribosomal DNA (ITS) and the plastid trnL-F region were conducted to reconstruct the phylogenetic history of Draba and related genera. Out of the approximately 370 Draba species, 169 geographically and morphologically representative species are sampled here, including such "controversial" segregates as Abdra, Arabis, Athysanus, Drabopsis, Erophila, Graellsia, Heterodraba, Schivereckia and Tomostima. Phylogenetic analyses of the combined ITS and trnL-F markers indicate that Draba represents a monophyletic evolutionary lineage within the tribe Arabideae, but previously proposed infrageneric classification systems are mostly artificial (i.e. do not reflect true evolutionary history). Draba muralis formed an independent genus, possibly between Draba and Arabis, whereas D. hystrix fell outside Draba and was closely associated with Arabis. The New World annual, lowland Draba, D. platycarpa, D. reptans, D. cuneifolia, D. australis, D. arabidoides, (section Abdra) and D. brachycarpa and D. aspera (section Tomostima), appear to be independent genera that fall outside Draba and are monophyletic with the endemic North American Heterodraba unilateralis and Athysanus pusillus. Graellsia hederifolia and Erophila verna appear to be earlier diverging Draba species, with weak evidence of ancient hybridization in G. hederifolia. Core Draba species were organized into three major groups that encompass the segregate genera Drabopsis and Schivereckia. The three groups have geographic significance: Group I-Europe to Iran; Group II-North and South American Cordillera; Group III-Asian, Arctic, and Beringian. These three groups also have significant petal color and base chromosome level similarities. We also found that Arabis rimarum is a synonym for Draba aucheri, and is well within Group I of Core Draba. Therefore, we propose that many of the lowland and coastal temperate species often associated with Draba should be excluded from that genus. The compilation of these data indicate that true Draba are highly migrating arctic and higher alpine species that are most often perennial with white petals and n=8 as plesiomorphic characters.

Molecular signatures of apomictic and sexual ovules in the Boechera holboellii complex

Apomixis, a natural form of asexual seed production in plants, has evolved independently in various taxa, and represents an important potential technology for agriculture. The switch to apomixis is based on de-regulation of developmental pathways originally leading to sexual seed formation. Hybridization and polyploidy, both typical characteristics of asexual plants and animals, are mechanisms that could trigger de-regulation. Here we show that up-regulation of alleles in apomeiotic ovules is mediated by genomic duplication, heterochrony and the residual effects of ancient hybridization in diploid apomicts of the Boechera holboellii complex. Using SuperSAGE, we have identified over 4000 differentially expressed mRNA tags between micro-dissected ovules from two diploid sexual (Boechera stricta and B. holboellii) and two diploid apomictic (Boechera divaricarpa) accessions. Pairwise sequence comparisons between tags enabled identification of allelic variants of the same loci. Up-regulated candidate apomeiosis alleles consistently have more than three related alleles, thus demonstrating transcription from duplicated loci. A further 543 alleles were heterochronically expressed between sexual and apomeiotic ovules at developmental stages 2-II to 2-IV. Intriguingly, 69 B. holboellii specific alleles were preferentially up-regulated in apomeiotic ovules, thus showing a remnant'parent of origin' effect stemming from the Pleistocene origin of the hybrid B. divaricarpa from taxa related to B. holboellii and B. stricta. These data implicate polyploid gene dosage in the expression of asexual seed formation, and support hypotheses of de-regulation of the sexual pathway. The observed 'parent of origin' effect suggests that the genomic memory of hybridization has somehow been maintained after hundreds, if not thousands, of asexual generations.

Ecological genomics of Boechera stricta: identification of a QTL controlling the allocation of methionine- vs branched-chain amino acid-derived glucosinolates and levels of insect herbivory

In the Brassicaceae glucosinolates influence feeding, reproduction and development of many insect herbivores. Glucosinolate production and effects on herbivore feeding have been extensively studied in the model species Arabidopsis thaliana and Brassica crops, both of which constitutively produce leaf glucosinolates mostly derived from the amino acid methionine. Much less is known about the regulation or role in defense of glucosinolates derived from other aliphatic amino acids such as the branched-chain amino acids (BCAA) valine and isoleucine. We have identified a glucosinolate polymorphism in Boechera stricta controlling the allocation to BCAA- vs. methionine-derived glucosinolates in both leaves and seeds. Boechera stricta is a perennial species that grows in mostly undisturbed habitats of western North America. We have measured glucosinolate profiles and concentrations in 192 F₂ lines that have previously been used for genetic map construction. We also performed herbivory assays on six F₃ replicates per F₂ line using the generalist lepidopteran Trichoplusia ni. Quantitative Trait Locus (QTL) analysis identified a single locus controlling both glucosinolate profile and levels of herbivory, the Branched Chain-Methionine Allocation or BCMA QTL. We have delimited this QTL to a small genomic region with a 1.0 LOD confidence interval just 1.9 cM wide, which in A. thaliana contains ∼100 genes. We also found that methionine-derived glucosinolates provided significantly greater defense than the BCAA-derived glucosinolates against feeding by this generalist insect herbivore. The future positional cloning of this locus will allow for testing various adaptive explanations.

Multilocus patterns of nucleotide diversity, population structure and linkage disequilibrium in Boechera stricta, a wild relative of Arabidopsis

Information about polymorphism, population structure, and linkage disequilibrium (LD) is crucial for association studies of complex trait variation. However, most genomewide studies have focused on model systems, with very few analyses of undisturbed natural populations. Here, we sequenced 86 mapped nuclear loci for a sample of 46 genotypes of Boechera stricta and two individuals of B. holboellii, both wild relatives of Arabidopsis. Isolation by distance was significant across the species range of B. stricta, and three geographic groups were identified by structure analysis, principal coordinates analysis, and distance-based phylogeny analyses. The allele frequency spectrum indicated a genomewide deviation from an equilibrium neutral model, with silent nucleotide diversity averaging 0.004. LD decayed rapidly, declining to background levels in approximately 10 kb or less. For tightly linked SNPs separated by <1 kb, LD was dependent on the reference population. LD was lower in the specieswide sample than within populations, suggesting that low levels of LD found in inbreeding species such as B. stricta, Arabidopsis thaliana, and barley may result from broad geographic sampling that spans heterogeneous genetic groups. Finally, analyses also showed that inbreeding B. stricta and A. thaliana have approximately 45% higher recombination per kilobase than outcrossing A. lyrata.