Spatial patterns and climatic drivers of phylogenetic structure of regional liverwort assemblages in China.

BACKGROUND AND AIMS: Latitudinal diversity gradients have been intimately linked to the tropical niche conservatism hypothesis, which posits that there has been a strong filter due to the challenges faced by ancestral tropical lineages to adapt to low temperatures and colonize extra-tropical regions. In liverworts, species richness is higher towards the tropics, but the centers of diversity of the basal lineages are distributed across extra-tropical regions, pointing to the colonization of tropical regions by phylogenetically clustered assemblages of species of temperate origin. Here, we test this hypothesis through analyses on the relationship between macroclimatic variation and phylogenetic diversity in Chinese liverworts. METHODS: Phylogenetic diversity metrics and their standardized effect sizes for liverworts in each of the 28 regional floras at the province level in China were related to latitude and six climate variables using regression analysis. We conducted variation partitioning analyses to determine the relative importance of each group of climatic variables. KEY RESULTS: We find that the number of species decreases with latitude, whereas phylogenetic diversity shows the reverse pattern, and that phylogenetic diversity is more strongly correlated with temperature-related variables, compared with precipitation-related variables. CONCLUSIONS: We interpret the opposite patterns observed in phylogenetic diversity and species richness in terms of a more recent origin of tropical diversity coupled with higher extinctions in temperate regions.

Disentangling climate change from air pollution effects on epiphytic bryophytes.

At the interface between atmosphere and vegetation, epiphytic floras have been largely used as indicators of air quality. The recovery of epiphytes from high levels of SO2 pollution has resulted in major range changes, whose interpretation has, however, been challenged by concomitant variation in other pollutants as well as climate change. Here, we combine historical and contemporary information on epiphytic bryophyte species distributions, climatic conditions, and pollution loads since the 1980s in southern Belgium to disentangle the relative impact of climate change and air pollution on temporal shifts in species composition. The relationship between the temporal variation of species composition, climatic conditions, SO2 , NO2 , O3 , and fine particle concentrations, was analyzed by variation partitioning. The temporal shift in species composition was such, that it was, on average, more than twice larger than the change in species composition observed today among communities scattered across the study area. The main driver, contributing to 38% of this temporal shift in species composition, was the variation of air quality. Climate change alone did not contribute to the substantial compositional shifts in epiphytic bryophyte communities in the course of the last 40 years. As a consequence of the substantial drop of N and S loads over the last decades, present-day variations of epiphytic floras were, however, better explained by the spatial variation of climatic conditions than by extant pollution loads. The lack of any signature of recolonization delays of formerly polluted areas in the composition of modern floras suggests that epiphytic bryophytes efficiently disperse at the landscape scale. We suggest that a monitoring of epiphyte communities at 10-year intervals would be desirable to assess the impact of raising pollution sources, and especially pesticides, whose impact on bryophytes remains poorly documented.

What do tropical cryptogams reveal? Strong genetic structure in Amazonian bryophytes.

Lowland tropical bryophytes have been perceived as excellent dispersers. In such groups, the inverse isolation hypothesis proposes that spatial genetic structure is erased beyond the limits of short-distance dispersal. Here, we determine the influence of environmental variation and geographic barriers on the spatial genetic structure of a widely dispersed and phylogenetically independent sample of Amazonian bryophytes. Single nucleotide polymorphism data were produced from a restriction site-associated DNA sequencing protocol for 10 species and analyzed through F-statistics and Mantel tests. Neither isolation-by-environment nor the impact of geographic barriers were recovered from the analyses. However, significant isolation-by-distance patterns were observed for 8 out of the 10 investigated species beyond the scale of short-distance dispersal (> 1 km), offering evidence contrary to the inverse isolation hypothesis. Despite a cadre of life-history traits and distributional patterns suggesting that tropical bryophytes are highly vagile, our analyses reveal spatial genetic structures comparable to those documented for angiosperms, whose diaspores are orders of magnitude larger. Dispersal limitation for tropical bryophytes flies in the face of traditional assumptions regarding their dispersal potential, and suggests that the plight of this component of cryptic biodiversity is more dire than previously considered in light of accelerated forest fragmentation in the Amazon.

Evidence for Gut-Associated Serratia symbiotica in Wild Aphids and Ants Provides New Perspectives on the Evolution of Bacterial Mutualism in Insects.

Many insects engage in symbiotic associations with diverse assemblages of bacterial symbionts that can deeply impact on their ecology and evolution. The intraspecific variation of symbionts remains poorly assessed while phenotypic effects and transmission behaviors, which are key processes for the persistence and evolution of symbioses, may differ widely depending on the symbiont strains. Serratia symbiotica is one of the most frequent symbiont species in aphids and a valuable model to assess this intraspecific variation since it includes both facultative and obligate symbiotic strains. Despite evidence that some facultative S. symbiotica strains exhibit a free-living capacity, the presence of these strains in wild aphid populations, as well as in insects with which they maintain regular contact, has never been demonstrated. Here, we examined the prevalence, diversity, and tissue tropism of S. symbiotica in wild aphids and associated ants. We found a high occurrence of S. symbiotica infection in ant populations, especially when having tended infected aphid colonies. We also found that the S. symbiotica diversity includes strains found located within the gut of aphids and ants. In the latter, this tissue tropism was found restricted to the proventriculus. Altogether, these findings highlight the extraordinary diversity and versatility of an insect symbiont and suggest the existence of novel routes for symbiont acquisition in insects.

Quantification of complex modular architecture in plants.

Morphometrics, the assignment of quantities to biological shapes, is a powerful tool to address taxonomic, evolutionary, functional and developmental questions. We propose a novel method for shape quantification of complex modular architecture in thalloid plants, whose extremely reduced morphologies, combined with the lack of a formal framework for thallus description, have long rendered taxonomic and evolutionary studies extremely challenging. Using graph theory, thalli are described as hierarchical series of nodes and edges, allowing for accurate, homologous and repeatable measurements of widths, lengths and angles. The computer program MorphoSnake was developed to extract the skeleton and contours of a thallus and automatically acquire, at each level of organization, width, length, angle and sinuosity measurements. Through the quantification of leaf architecture in Hymenophyllum ferns (Polypodiopsida) and a fully worked example of integrative taxonomy in the taxonomically challenging thalloid liverwort genus Riccardia, we show that MorphoSnake is applicable to all ramified plants. This new possibility of acquiring large numbers of quantitative traits in plants with complex modular architectures opens new perspectives of applications, from the development of rapid species identification tools to evolutionary analyses of adaptive plasticity.

Evolutionary origin of the latitudinal diversity gradient in liverworts.

A latitudinal diversity gradient towards the tropics appears as one most recurrent patterns in ecology, but the mechanisms underlying this pattern remain an area of controversy. In angiosperms, the tropical conservatism hypothesis proposes that most groups originated in the tropics and are adapted to a tropical climatic regime, and that relatively few species have evolved physiological adaptations to cold, dry or unpredictable climates. This mechanism is, however, unlikely to apply across land plants, and in particular, to liverworts, a group of about 7500 species, whose ability to withstand cold much better than their tracheophyte counterparts is at odds with the tropical conservatism hypothesis. Molecular dating, diversification rate analyses and ancestral area reconstructions were employed to explore the evolutionary mechanisms that account for the latitudinal diversity gradient in liverworts. As opposed to angiosperms, tropical liverwort genera are not older than their extra-tropical counterparts (median stem age of tropical and extra-tropical liverwort genera of 24.35 ±â€¯39.65 Ma and 39.57 ±â€¯49.07 Ma, respectively), weakening the 'time for speciation hypothesis'. Models of ancestral area reconstructions with equal migration rates between tropical and extra-tropical regions outperformed models with asymmetrical migration rates in either direction. The symmetry and intensity of migrations between tropical and extra-tropical regions suggested by the lack of resolution in ancestral area reconstructions towards the deepest nodes are at odds with the tropical niche conservatism hypothesis. In turn, tropical genera exhibited significantly higher net diversification rates than extra-tropical ones, suggesting that the observed latitudinal diversity gradient results from either higher extinction rates in extra-tropical lineages or higher speciation rates in the tropics. We discuss a series of experiments to help deciphering the underlying evolutionary mechanisms.

Range size heritability and diversification patterns in the liverwort genus Radula.

Why some species exhibit larger geographical ranges than others, and to what extent does variation in range size affect diversification rates, remains a fundamental, but largely unanswered question in ecology and evolution. Here, we implement phylogenetic comparative analyses and ancestral area estimations in Radula, a liverwort genus of Cretaceous origin, to investigate the mechanisms that explain differences in geographical range size and diversification rates among lineages. Range size was phylogenetically constrained in the two sub-genera characterized by their almost complete Australasian and Neotropical endemicity, respectively. The congruence between the divergence time of these lineages and continental split suggests that plate tectonics could have played a major role in their present distribution, suggesting that a strong imprint of vicariance can still be found in extant distribution patterns in these highly mobile organisms. Amentuloradula, Volutoradula and Metaradula species did not appear to exhibit losses of dispersal capacities in terms of dispersal life-history traits, but evidence for significant phylogenetic signal in macroecological niche traits suggests that niche conservatism accounts for their restricted geographic ranges. Despite their greatly restricted distribution to Australasia and Neotropics respectively, Amentuloradula and Volutoradula did not exhibit significantly lower diversification rates than more widespread lineages, in contrast with the hypothesis that the probability of speciation increases with range size by promoting geographic isolation and increasing the rate at which novel habitats are encountered. We suggest that stochastic long-distance dispersal events may balance allele frequencies across large spatial scales, leading to low genetic structure among geographically distant areas or even continents, ultimately decreasing the diversification rates in highly mobile, widespread lineages.

Increased diversification rates follow shifts to bisexuality in liverworts.

Shifts in sexual systems are one of the key drivers of species diversification. In contrast to angiosperms, unisexuality prevails in bryophytes. Here, we test the hypotheses that bisexuality evolved from an ancestral unisexual condition and is a key innovation in liverworts. We investigate whether shifts in sexual systems influence diversification using hidden state speciation and extinction analysis (HiSSE). This new method compares the effects of the variable of interest to the best-fitting latent variable, yielding robust and conservative tests. We find that the transitions in sexual systems are significantly biased toward unisexuality, even though bisexuality is coupled with increased diversification. Sexual systems are strongly conserved deep within the liverwort tree but become much more labile toward the present. Bisexuality appears to be a key innovation in liverworts. Its effects on diversification are presumably mediated by the interplay of high fertilization rates, massive spore production and long-distance dispersal, which may separately or together have facilitated liverwort speciation, suppressed their extinction, or both. Importantly, shifts in liverwort sexual systems have the opposite effect when compared to angiosperms, leading to contrasting diversification patterns between the two groups. The high prevalence of unisexuality among liverworts suggests, however, a strong selection for sexual dimorphism.

High migration rates shape the postglacial history of amphi-Atlantic bryophytes.

Paleontological evidence and current patterns of angiosperm species richness suggest that European biota experienced more severe bottlenecks than North American ones during the last glacial maximum. How well this pattern fits other plant species is less clear. Bryophytes offer a unique opportunity to contrast the impact of the last glacial maximum in North America and Europe because about 60% of the European bryoflora is shared with North America. Here, we use population genetic analyses based on approximate Bayesian computation on eight amphi-Atlantic species to test the hypothesis that North American populations were less impacted by the last glacial maximum, exhibiting higher levels of genetic diversity than European ones and ultimately serving as a refugium for the postglacial recolonization of Europe. In contrast with this hypothesis, the best-fit demographic model involved similar patterns of population size contractions, comparable levels of genetic diversity and balanced migration rates between European and North American populations. Our results thus suggest that bryophytes have experienced comparable demographic glacial histories on both sides of the Atlantic. Although a weak, but significant genetic structure was systematically recovered between European and North American populations, evidence for migration from and towards both continents suggests that amphi-Atlantic bryophyte population may function as a metapopulation network. Reconstructing the biogeographic history of either North American or European bryophyte populations therefore requires a large, trans-Atlantic geographic framework.

Is the sword moss (Bryoxiphium) a preglacial Tertiary relict?

The disjunction of floras between East Asia, Southeast North America, West North America, and Southwest Eurasia has been interpreted in terms of the fragmentation of a once continuous mixed mesophytic forest that occurred throughout the Northern Hemisphere due to the climatic and geological changes during the late Tertiary. The sword moss, Bryoxiphium, exhibits a distribution that strikingly resembles that of the mesophytic forest elements such as Liriodendron and is considered as the only living member of an early Tertiary flora in Iceland. These hypotheses are tested here using molecular dating analyses and ancestral area estimations. The results suggest that the extant range of Bryoxiphium results from the fragmentation of a formerly wider range encompassing North America and Southeast Asia about 10 million years ago. The split of continental ancestral populations is too recent to match with a continental drift scenario but is spatially and temporally remarkably congruent with that observed in Tertiary angiosperm relict species. The timing of the colonization of Iceland from Macaronesian ancestors, about two million years ago, is, however, incompatible with the hypothesis that Bryoxiphium is the only living member of an early Tertiary flora of the island. Alaska was recurrently colonized from East Asia. The ability of Bryoxiphium to overcome large oceanic barriers is further evidenced by its occurrence on remote oceanic archipelagos. In particular, Madeira was colonized twice independently from American and East Asian ancestors, respectively. The striking range disjunction of Bryoxiphium is interpreted in terms of its mating system, as the taxon exhibits a very singular pattern of spatial segregation of the sexes.

What are the evolutionary mechanisms explaining the similar species richness patterns in tropical mosses? Insights from the phylogeny of the pantropical genus Pelekium.

As opposed to angiosperms, moss species richness is similar among tropical regions of the world, in line with the hypothesis that tropical bryophytes are extremely good dispersers. Here, we reconstructed the phylogeny of the pantropical moss genus Pelekium to test the hypothesis that high migration rates erase any difference in species richness among tropical regions. In contrast with this hypothesis, several species considered to have a pantropical range were resolved as a complex of species with a strong geographic structure. Consequently, a significant phylogeographical signal was found in the data, evidencing that cladogenetic diversification within regions takes place at a faster rate than intercontinental migration. The shape of the Pelekium phylogeny, along with the selection of a constant-rate model of diversification among species in the genus, suggests, however, that the cladogenetic speciation patterns observed in Pelekium are not comparable to some of the spectacular examples of tropical radiations reported in angiosperms. Rather, the results presented here point to the constant accumulation of diversity through time in Pelekium. This, combined with evidence for long-distance dispersal limitations in the genus, suggests that the similar patterns of species richness among tropical areas are better explained in terms of comparable rates of diversification across tropical regions than by the homogenization of species richness by recurrent migrations.

Approximate Bayesian Computation Reveals the Crucial Role of Oceanic Islands for the Assembly of Continental Biodiversity.

The perceived low levels of genetic diversity, poor interspecific competitive and defensive ability, and loss of dispersal capacities of insular lineages have driven the view that oceanic islands are evolutionary dead ends. Focusing on the Atlantic bryophyte flora distributed across the archipelagos of the Azores, Madeira, the Canary Islands, Western Europe, and northwestern Africa, we used an integrative approach with species distribution modeling and population genetic analyses based on approximate Bayesian computation to determine whether this view applies to organisms with inherent high dispersal capacities. Genetic diversity was found to be higher in island than in continental populations, contributing to mounting evidence that, contrary to theoretical expectations, island populations are not necessarily genetically depauperate. Patterns of genetic variation among island and continental populations consistently fitted those simulated under a scenario of de novo foundation of continental populations from insular ancestors better than those expected if islands would represent a sink or a refugium of continental biodiversity. We, suggest that the northeastern Atlantic archipelagos have played a key role as a stepping stone for transoceanic migrants. Our results challenge the traditional notion that oceanic islands are the end of the colonization road and illustrate the significant role of oceanic islands as reservoirs of novel biodiversity for the assembly of continental floras.

Measuring spore settling velocity for an improved assessment of dispersal rates in mosses.

BACKGROUND AND AIMS: The settling velocity of diaspores is a key parameter for the measurement of dispersal ability in wind-dispersed plants and one of the most relevant parameters in explicit dispersal models, but remains largely undocumented in bryophytes. The settling velocities of moss spores were measured and it was determined whether settling velocities can be derived from spore diameter using Stokes' Law or if specific traits of spore ornamentation cause departures from theoretical expectations. METHODS: A fall tower design combined with a high-speed camera was used to document spore settling velocities in nine moss species selected to cover the range of spore diameters within the group. Linear mixed effect models were employed to determine whether settling velocity can be predicted from spore diameter, taking specific variation in shape and surface roughness into account. KEY RESULTS: Average settling velocity of moss spores ranged from 0·49 to 8·52 cm s(-1) There was a significant positive relationship between spore settling velocity and size, but the inclusion of variables of shape and texture of spores in the best-fit models provides evidence for their role in shaping spore settling velocities. CONCLUSIONS: Settling velocities in mosses can significantly depart from expectations derived from Stokes' Law. We suggest that variation in spore shape and ornamentation affects the balance between density and drag, and results in different dispersal capacities, which may be correlated with different life-history traits or ecological requirements. Further studies on spore ultrastructure would be necessary to determine the role of complex spore ornamentation patterns in the drag-to-mass ratio and ultimately identify what is the still poorly understood function of the striking and highly variable ornamentation patterns of the perine layer on moss spores.

The anagenetic world of spore-producing land plants.

A fundamental challenge to our understanding of biodiversity is to explain why some groups of species diversify, whereas others do not. On islands, the gradual evolution of a new species from a founder event has been called 'anagenetic speciation'. This process does not lead to rapid and extensive speciation within lineages and has received little attention. Based on a survey of the endemic bryophyte, pteridophyte and spermatophyte floras of nine oceanic archipelagos, we show that anagenesis, as measured by the proportion of genera with single endemic species within a genus, is much higher in bryophytes (73%) and pteridophytes (65%) than in spermatophytes (55%). Anagenesis contributed 49% of bryophyte and 40% of endemic pteridophyte species, but only 17% of spermatophytes. The vast majority of endemic bryophytes and pteridophytes are restricted to subtropical evergreen laurel forests and failed to diversify in more open environments, in contrast with the pattern exhibited by spermatophytes. We propose that the dominance of anagenesis in island bryophytes and pteridophytes is a result of a mixture of intrinsic factors, notably their strong preference for (sub)tropical forest environments, and extrinsic factors, including the long-term macro-ecological stability of these habitats and the associated strong phylogenetic niche conservatism of their floras.

Oceanic islands are not sinks of biodiversity in spore-producing plants.

Islands have traditionally been considered as migratory and evolutionary dead ends for two main reasons: island colonizers are typically assumed to lose their dispersal power, and continental back colonization has been regarded as unlikely because of niche preemption. The hypothesis that islands might actually represent dynamic refugia and migratory stepping stones for species that are effective dispersers, and in particular, for spore-producing plants, is formally tested here, using the archipelagos of the Azores, Canary Islands, and Madeira, as a model. Population genetic analyses based on nuclear microsatellite variation indicate that dispersal ability of the moss Platyhypnidium riparioides does not decrease in the island setting. The analyses further show that, unlike island populations, mainland (southwestern Europe and North Africa) populations underwent a severe bottleneck during the last glacial maximum (LGM). Our results thus refute the traditional view of islands as the end of the colonization road and point to a different perception of North Atlantic archipelagos as major sources of biodiversity for the postglacial recolonization of Europe by spore-producing plants.

To what extent can we predict variation of bryophyte and tracheophyte community composition at fine spatial scale along an elevation gradient?

Mounting evidence points to the need for high-resolution climatic data in biodiversity analyses under global change. As we move to finer resolution, other factors than climate, including other abiotic variables and biotic interactions play, however, an increasing role, raising the question of our ability to predict community composition at fine scales. Focusing on two lineages of land plants, bryophytes and tracheophytes, we determine the relative contribution of climatic, non-climatic environmental drivers, spatial effects, community architecture and composition of one lineage to predict community composition of the other lineage, and how our ability to predict community composition varies along an elevation gradient. The relationship between community composition of one lineage and 68 environmental variables at 2-25 m spatial resolution, architecture and composition of the other lineage, and spatial factors, was investigated by hierarchical and variance partitioning across 413 2x2m plots in the Swiss Alps. Climatic data, although significant, contributed less to the model than any other variable considered. Community composition of one lineage, reflecting both direct interactions and unmeasured (hidden) abiotic factors, was the best predictor of community composition of the other lineage. Total explained variance substantially varied with elevation, underlining the fact that the strength of the species composition-environment relationship varies depending on environmental conditions. Total variance explained increased towards high elevation up to 50 %, with an increasing importance of spatial effects and vegetation architecture, pointing to increasing positive interactions and aggregated species distribution patterns in alpine environments. In tracheophytes, an increase of the contribution of non-climatic environmental factors was also observed at high elevation, in line with the hypothesis of a stronger environmental control under harsher conditions. Further improvements of our ability to predict changes in plant community composition may involve the implementation of historical variables and higher-resolution climatic data to better describe the microhabitat conditions actually experienced by organisms.

Elevational patterns of genetic variation in the cosmopolitan moss Bryum argenteum (Bryaceae).

PREMISE OF THE STUDY: The Baas Becking tenet posits that 'everything is everywhere, but the environment selects' to explain cosmopolitan distributions in highly vagile taxa. Bryophyte species show wider distributions than vascular plants and include examples of truly cosmopolitan ranges, which have been interpreted as a result of high dispersal capacities and ecological plasticity. In the current study, we documented patterns of genetic structure and diversity in the cosmopolitan moss Bryum argenteum along an elevational gradient to determine if genetic diversity and structure is homogenized by intense migrations in the lack of ecological differentiation. • METHODS: 60 specimens were collected in the Sierra Nevada Mountains (Spain) between 100 and 2870 m and sequenced for ITS and rps4. Comparative analyses, genetic diversity estimators, and Mantel's tests were employed to determine the relationship between genetic variation, elevation, and geographic distance and to look for signs of demographic shifts. • KEY RESULTS: Genetic diversity peaked above 1900 m and no signs of demographic shifts were detected at any elevation. There was a strong phylogenetic component in elevational variation. Genetic variation was significantly correlated with elevation, but not with geographic distance. • CONCLUSIONS: The results point to the long-term persistence of Bryum argenteum in a range that was glaciated during the Late Pleistocene. Evidence for an environmentally driven pattern of genetic differentiation suggests adaptive divergence. This supports the Baas Becking tenet and indicates that ecological specialization might play a key role in explaining patterns of genetic structure in cosmopolitan mosses.

Morphology informed by phylogeny reveals unexpected patterns of species differentiation in the aquatic moss Rhynchostegium riparioides s.l.

Bryophyte floras typically exhibit extremely low levels of endemism. The interpretation, that this might reflect taxonomic shortcomings, is tested here for the Macaronesian flora, using the moss species complex of Rhynchostegium riparioides as a model. The deep polyphyly of R. riparioides across its distribution range reveals active differentiation that better corresponds to geographic than morphological differences. Morphometric analyses are, in fact, blurred by a size gradient that accounts for 80% of the variation observed among gametophytic traits. The lack of endemic diversification observed in R. riparioides in Macaronesia weakens the idea that the low rates of endemism observed in the Macaronesian bryophyte flora might solely be explained by taxonomic shortcomings. To the reverse, the striking polyphyly of North American and European lineages of R. riparioides suggests that the similarity between the floras of these continents has been over-emphasized. Discriminant analyses point to the existence of morphological discontinuities among the lineages resolved by the molecular phylogeny. The global rate of error associated to species identification based on morphology (0.23) indicates, however, that intergradation of shape and size characters among species in the group challenges their identification.

How do temperate bryophytes face the challenge of a changing environment? Lessons from the past and predictions for the future.

Bryophytes are a group of early land plants, whose specific ecophysiological and biological features, including poikilohydry, sensitivity to moderately high temperature and high dispersal ability, make them ideal candidates for investigating the impact of climate changes. Employing a combined approach of species distribution modelling (SDM) and molecular phylogeography in the temperate moss Homalothecium sericeum, we explore the significance of the Mediterranean refugia, contrasting the southern and northern refugia hypotheses, determine the extent to which recolonization of previously glaciated areas has been facilitated by the high dispersal ability of the species and make predictions on the extent to which it will be impacted by ongoing climate change. The Mediterranean areas exhibit the highest nucleotidic diversities and host a mixture of ancestral, endemic and more recently derived haplotypes. Extra-Mediterranean areas exhibit low genetic diversities and Euro-Siberian populations display a significant signal of expansion that is identified to be of Euro-Siberian origin, pointing to the northern refugia hypothesis. The SDMs predict a global net increase in range size owing to ongoing climate change, but substantial range reductions in southern areas. Presence of a significant phylogeographical signal at different spatial scales suggests, however, that dispersal limitations might constitute, as opposed to the traditional view of spore-producing plants as efficient dispersers, a constraint for migration. This casts doubts about the ability of the species to face the massive extinctions predicted in the southern areas, threatening their status of reservoir of genetic diversity.

Evolution of sexual systems, dispersal strategies and habitat selection in the liverwort genus Radula.

• Shifts in sexual systems are among the most common and important transitions in plants and are correlated with a suite of life-history traits. The evolution of sexual systems and their relationships to gametophyte size, sexual and asexual reproduction, and epiphytism are examined here in the liverwort genus Radula. • The sequence of trait acquisition and the phylogenetic correlations between those traits was investigated using comparative methods. • Shifts in sexual systems recurrently occurred from dioecy to monoecy within facultative epiphyte lineages. Production of specialized asexual gemmae was correlated to neither dioecy nor strict epiphytism. • The significant correlations among life-history traits related to sexual systems and habitat conditions suggest the existence of evolutionary trade-offs. Obligate epiphytes do not produce gemmae more frequently than facultative epiphytes and disperse by whole gametophyte fragments, presumably to avoid the sensitive protonemal stage in a habitat prone to rapid changes in moisture availability. As dispersal ranges correlate with diaspore size, this reinforces the notion that epiphytes experience strong dispersal limitations. Our results thus provide the evolutionary complement to metapopulation, metacommunity and experimental studies demonstrating trade-offs between dispersal distance, establishment ability, and life-history strategy, which may be central to the evolution of reproductive strategies in bryophytes.