Different molecular changes underlie the same phenotypic transition: Origins and consequences of independent shifts to homostyly within species

The repeated transition from outcrossing to selfing is a key topic in evolutionary biology. However, the molecular basis of such shifts has been rarely examined due to lack of knowledge of the genes controlling these transitions. A classic example of mating system transition is the repeated shift from heterostyly to homostyly. Occurring in 28 angiosperm families, heterostyly is characterized by the reciprocal position of male and female sexual organs in two (or three) distinct, usually self-incompatible floral morphs. Conversely, homostyly is characterized by a single, self-compatible floral morph with reduced separation of male and female organs, facilitating selfing. Here, we investigate the origins of homostyly in Primula vulgaris and its microevolutionary consequences by integrating surveys of the frequency of homostyles in natural populations, DNA sequence analyses of the gene controlling the position of female sexual organs (CYPᵀ), and microsatellite genotyping of both progeny arrays and natural populations characterized by varying frequencies of homostyles. As expected, we found that homostyles displace short-styled individuals, but long-style morphs are maintained at low frequencies within populations. We also demonstrated that homostyles repeatedly evolved from short-styled individuals in association with different types of loss-of- function mutations in CYPᵀ. Additionally, homostyly triggers a shift to selfing, promoting increased inbreeding within and genetic differentiation among populations. Our results elucidate the causes and consequences of repeated transitions to homostyly within species, and the putative mechanisms precluding its fixation in P. vulgaris. This study represents a benchmark for future analyses of losses of heterostyly in other angiosperms

A new classification of Cyperaceae (Poales) supported by phylogenomic data

Cyperaceae (sedges) are the third largest monocot family and are of considerable economic and ecological importance. Sedges represent an ideal model family to study evolutionary biology due to their species richness, global distribution, large discrepancies in lineage diversity, broad range of ecological preferences, and adaptations including multiple origins of C4 photosynthesis and holocentric chromosomes. Goetghebeur′s seminal work on Cyperaceae published in 1998 provided the most recent complete classification at tribal and generic level, based on a morphological study of Cyperaceae inflorescence, spikelet, flower, and embryo characters, plus anatomical and other information. Since then, several family-level molecular phylogenetic studies using Sanger sequence data have been published. Here, more than 20 years after the last comprehensive classification of the family, we present the first family-wide phylogenomic study of Cyperaceae based on targeted sequencing using the Angiosperms353 probe kit sampling 311 accessions. In addition, 62 accessions available from GenBank were mined for overlapping reads and included in the phylogenomic analyses. Informed by this backbone phylogeny, a new classification for the family at the tribal, subtribal, and generic levels is proposed. The majority of previously recognized suprageneric groups are supported, and for the first time, we establish support for tribe Cryptangieae as a clade including the genus Koyamaea. We provide a taxonomic treatment including identification keys and diagnoses for the 2 subfamilies, 24 tribes, and 10 subtribes, and basic information on the 95 genera. The classification includes five new subtribes in tribe Schoeneae: Anthelepidinae, Caustiinae, Gymnoschoeninae, Lepidospermatinae, and Oreobolinae

Whole-genome analyses disentangle reticulate evolution of primroses in a biodiversity hotspot

Biodiversity hotspots, such as the Caucasus mountains, provide unprecedented opportunities for understanding the evolutionary processes that shape species diversity and richness. Therefore, we investigated the evolution of Primula sect. Primula, a clade with a high degree of endemism in the Caucasus. We performed phylogenetic and network analyses of whole-genome resequencing data from the entire nuclear genome, the entire chloroplast genome, and the entire heterostyly supergene. The different characteristics of the genomic partitions and the resulting phylogenetic incongruences enabled us to disentangle evolutionary histories resulting from tokogenetic vs cladogenetic processes. We provide the first phylogeny inferred from the heterostyly supergene that includes all species of Primula sect. Primula. Our results identified recurrent admixture at deep nodes between lineages in the Caucasus as the cause of non-monophyly in Primula. Biogeographic analyses support the 'out-of-the-Caucasus' hypothesis, emphasizing the importance of this hotspot as a cradle for biodiversity. Our findings provide novel insights into causal processes of phylogenetic discordance, demonstrating that genome-wide analyses from partitions with contrasting genetic characteristics and broad geographic sampling are crucial for disentangling the diversification of species-rich clades in biodiversity hotspots

A tale of worldwide success: Behind the scenes of Carex (Cyperaceae) biogeography and diversification

The megadiverse genus Carex (c. 2000 species, Cyperaceae) has a nearly cosmopolitan distribution, displaying an inverted latitudinal richness gradient with higher species diversity in cold-temperate areas of the Northern Hemisphere. Despite great expansion in our knowledge of the phylogenetic history of the genus and many molecular studies focusing on the biogeography of particular groups during the last few decades, a global analysis of Carex biogeography and diversification is still lacking. For this purpose, we built the hitherto most comprehensive Carex-dated phylogeny based on three markers (ETS–ITS–matK), using a previous phylogenomic Hyb-Seq framework, and a sampling of two-thirds of its species and all recognized sections. Ancestral area reconstruction, biogeographic stochastic mapping, and diversification rate analyses were conducted to elucidate macroevolutionary biogeographic and diversification patterns. Our results reveal that Carex originated in the late Eocene in E Asia, where it probably remained until the synchronous diversification of its main subgeneric lineages during the late Oligocene. E Asia is supported as the cradle of Carex diversification, as well as a “museum” of extant species diversity. Subsequent “out-of-Asia” colonization patterns feature multiple asymmetric dispersals clustered toward present times among the Northern Hemisphere regions, with major regions acting both as source and sink (especially Asia and North America), as well as several independent colonization events of the Southern Hemisphere. We detected 13 notable diversification rate shifts during the last 10 My, including remarkable radiations in North America and New Zealand, which occurred concurrently with the late Neogene global cooling, which suggests that diversification involved the colonization of new areas and expansion into novel areas of niche space.This work was carried out with financial support by the National Science Foundation (Award #1255901 to ALH and Award #1256033 to EHR), the Spanish Ministry of Economy and Competitiveness (project CGL2016–77401‐P to SM-B and ML), the USDA National Institute of Food and Agriculture (McIntire Stennis project 1018692 to DS) as well as postdoctoral fellowships towards SM‐B (Universidad Pablo de Olavide, PP16/12‐APP), and PJ‐M (National Science Foundation, Award #1256033, and the Smithsonian Postdoctoral Fellowship program)

A new classification of Cyperaceae (Poales) supported by phylogenomic data

Cyperaceae (sedges) are the third largest monocot family and are of considerable economic and ecological importance. Sedges represent an ideal model family to study evolutionary biology because of their species richness, global distribution, large discrepancies in lineage diversity, broad range of ecological preferences, and adaptations including multiple origins of C4 photosynthesis and holocentric chromosomes. Goetghebeur’s seminal work on Cyperaceae published in 1998 provided the most recent complete classification at tribal and generic level, based on a morphological study of Cyperaceae inflorescence, spikelet, flower and embryo characters plus anatomical and other information. Since then, several family‐level molecular phylogenetic studies using Sanger sequence data have been published. Here, more than 20 years after the last comprehensive classification of the family, we present the first family‐wide phylogenomic study of Cyperaceae based on targeted sequencing using the Angiosperms353 probe kit sampling 311 accessions. Additionally, 62 accessions available from GenBank were mined for overlapping reads and included in the phylogenomic analyses. Informed by this backbone phylogeny, a new classification for the family at the tribal, subtribal and generic levels is proposed. The majority of previously recognized suprageneric groups are supported, and for the first time we establish support for tribe Cryptangieae as a clade including the genus Koyamaea. We provide a taxonomic treatment including identification keys and diagnoses for the 2 subfamilies, 24 tribes and 10 subtribes and basic information on the 95 genera. The classification includes five new subtribes in tribe Schoeneae: Anthelepidinae, Caustiinae, Gymnoschoeninae, Lepidospermatinae and Oreobolinae. This article is protected by copyright. All rights reserved

Suppl. Fig. 1. matK tree

Maximum likelihood topology based on matK sequence data alone, with parsimony/likelihood bootstrap percentages for every branch with >50% parsimony or ML bootstrap. Branches with >85% parsimony bootstrap support are emphasized with bold lines. A dash (-) indicates <50% bootstrap

Suppl. Fig. 8. ETS indels tree

Maximum likelihood topology based on ETS-1f indels alone, with parsimony/likelihood bootstrap percentages for every branch with >50% parsimony or ML bootstrap. Branches with >85% parsimony bootstrap support are emphasized with bold lines. A dash (-) indicates <50% bootstrap

Suppl. Fig. 11. morphology tree

Maximum likelihood topology based on morphology alone, with parsimony/likelihood bootstrap percentages for every branch with >50% parsimony or ML bootstrap. Branches with >85% parsimony bootstrap support are emphasized with bold lines. A dash (-) indicates <50% bootstrap

Suppl. data: alignments

NEXUS concatenated data matrix including DNA sequence, indel and morphological partitions, and relaxed PHYLIP alignments of each individual partition

Suppl. Fig. 3. ndhF tree

Maximum likelihood topology based on ndhF sequence data alone, with parsimony/likelihood bootstrap percentages for every branch with >50% parsimony or ML bootstrap. Branches with >85% parsimony bootstrap support are emphasized with bold lines. A dash (-) indicates <50% bootstrap