Diverse mating consequences of the evolutionary breakdown of the sexual polymorphism heterostyly.

Reproductive systems of flowering plants are evolutionarily fluid, with mating patterns changing in response to shifts in abiotic conditions, pollination systems, and population characteristics. Changes in mating should be particularly evident in species with sexual polymorphisms that become ecologically destabilized, promoting transitions to alternative reproductive systems. Here, we decompose female mating portfolios (incidence of selfing, outcross mate number, and intermorph mating) in eight populations of Primula oreodoxa, a self-compatible insect-pollinated herb. This species is ancestrally distylous, with populations subdivided into two floral morphs that usually mate with each other (disassortative mating). Stages in the breakdown of polymorphism also occur, including "mixed" populations of distylous and homostylous (self-pollinating) morphs and purely homostylous populations. Population morph ratios vary with elevation in association with differences in pollinator availability, providing an unusual opportunity to investigate changes in mating patterns accompanying transitions in reproductive systems. Unexpectedly, individuals mostly outcrossed randomly, with substantial disassortative mating in at most two distylous populations. As predicted, mixed populations had higher selfing rates than distylous populations, within mixed populations, homostyles selfed almost twice as much as the distylous morphs, and homostylous populations exhibited the highest selfing rates. Populations with homostyles outcrossed with fewer mates and mate number varied negatively with population selfing rates. These differences indicate maintenance of distyly at low elevation, transition to monomorphic selfing at high elevation, and uncertain, possibly variable fates at intermediate elevation. By quantifying the earliest changes in mating that initiate reproductive transitions, our study highlights the key role of mating in promoting evolutionary divergence.

Phased Assembly of Neo-Sex Chromosomes Reveals Extensive Y Degeneration and Rapid Genome Evolution in Rumex hastatulus.

Y chromosomes are thought to undergo progressive degeneration due to stepwise loss of recombination and subsequent reduction in selection efficiency. However, the timescales and evolutionary forces driving degeneration remain unclear. To investigate the evolution of sex chromosomes on multiple timescales, we generated a high-quality phased genome assembly of the massive older (<10 MYA) and neo (<200,000 yr) sex chromosomes in the XYY cytotype of the dioecious plant Rumex hastatulus and a hermaphroditic outgroup Rumex salicifolius. Our assemblies, supported by fluorescence in situ hybridization, confirmed that the neo-sex chromosomes were formed by two key events: an X-autosome fusion and a reciprocal translocation between the homologous autosome and the Y chromosome. The enormous sex-linked regions of the X (296 Mb) and two Y chromosomes (503 Mb) both evolved from large repeat-rich genomic regions with low recombination; however, the complete loss of recombination on the Y still led to over 30% gene loss and major rearrangements. In the older sex-linked region, there has been a significant increase in transposable element abundance, even into and near genes. In the neo-sex-linked regions, we observed evidence of extensive rearrangements without gene degeneration and loss. Overall, we inferred significant degeneration during the first 10 million years of Y chromosome evolution but not on very short timescales. Our results indicate that even when sex chromosomes emerge from repetitive regions of already-low recombination, the complete loss of recombination on the Y chromosome still leads to a substantial increase in repetitive element content and gene degeneration.

Parallel evolution of morphological and genomic selfing syndromes accompany the breakdown of heterostyly.

Evolutionary transitions from outcrossing to selfing in flowering plants have convergent morphological and genomic signatures and can involve parallel evolution within related lineages. Adaptive evolution of morphological traits is often assumed to evolve faster than nonadaptive features of the genomic selfing syndrome. We investigated phenotypic and genomic changes associated with transitions from distyly to homostyly in the Primula oreodoxa complex. We determined whether the transition to selfing occurred more than once and investigated stages in the evolution of morphological and genomic selfing syndromes using 22 floral traits and both nuclear and plastid genomic data from 25 populations. Two independent transitions were detected representing an earlier and a more recently derived selfing lineage. The older lineage exhibited classic features of the morphological and genomic selfing syndrome. Although features of both selfing syndromes were less developed in the younger selfing lineage, they exhibited parallel development with the older selfing lineage. This finding contrasts with the prediction that some genomic changes should lag behind adaptive changes to morphological traits. Our findings highlight the value of comparative studies on the timing and extent of transitions from outcrossing to selfing between related lineages for investigating the tempo of morphological and molecular evolution.

Genomic evidence supports the genetic convergence of a supergene controlling the distylous floral syndrome.

Heterostyly, a plant sexual polymorphism controlled by the S-locus supergene, has evolved numerous times among angiosperm lineages and represents a classic example of convergent evolution in form and function. Determining whether underlying molecular convergence occurs could provide insights on constraints to floral evolution. Here, we investigated S-locus genes in distylous Gelsemium (Gelsemiaceae) to determine whether there is evidence of molecular convergence with unrelated distylous species. We used several approaches, including anatomical measurements of sex-organ development and transcriptome and whole-genome sequencing, to identify components of the S-locus supergene. We also performed evolutionary analysis with candidate S-locus genes and compared them with those reported in Primula and Turnera. The candidate S-locus supergene of Gelsemium contained four genes, of which three appear to have originated from gene duplication events within Gelsemiaceae. The style-length genes GeCYP in Gelsemium and CYP734A50 in Primula likely arose from duplication of the same gene, CYP734A1. Three out of four S-locus genes in Gelsemium elegans were hemizygous, as previously reported in Primula and Turnera. We provide genomic evidence on the genetic convergence of the supergene underlying distyly among distantly related angiosperm lineages and help to illuminate the genetic architecture involved in the evolution of heterostyly.

Haplotype-resolved genome assembly provides insights into the evolution of S-locus supergene in distylous Nymphoides indica.

Distyly has evolved independently in numerous animal-pollinated angiosperm lineages. Understanding of its molecular basis has been restricted to a few species, primarily Primula. Here, we investigate the genetic architecture of the single diallelic locus (S-locus) supergene, a linkage group of functionally associated genes, and explore how it may have evolved in distylous Nymphoides indica, a lineage of flowering plants not previously investigated. We assembled haplotype-resolved genomes, used read-coverage-based genome-wide association study (rb-GWAS) to locate the S-locus supergene, co-expression network analysis to explore gene networks underpinning the development of distyly, and comparative genomic analyses to investigate the origins of the S-locus supergene. We identified three linked candidate S-locus genes - NinBAS1, NinKHZ2, and NinS1 - that were only evident in the short-styled morph and were hemizygous. Co-expression network analysis suggested that brassinosteroids contribute to dimorphic sex organs in the short-styled morph. Comparative genomic analyses indicated that the S-locus supergene likely evolved via stepwise duplications and has been affected by transposable element activities. Our study provides novel insight into the structure, regulation, and evolution of the supergene governing distyly in N. indica. It also provides high-quality genomic resources for future research on the molecular mechanisms underlying the striking evolutionary convergence in form and function across heterostylous taxa.

The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly.

The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous-homostylous species and identified patterns of nonsynonymous to synonymous divergence (dN/dS) and polymorphism (πN/πS) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.

Widespread Recombination Suppression Facilitates Plant Sex Chromosome Evolution.

Classical models suggest that recombination rates on sex chromosomes evolve in a stepwise manner to localize sexually antagonistic variants in the sex in which they are beneficial, thereby lowering rates of recombination between X and Y chromosomes. However, it is also possible that sex chromosome formation occurs in regions with preexisting recombination suppression. To evaluate these possibilities, we constructed linkage maps and a chromosome-scale genome assembly for the dioecious plant Rumex hastatulus. This species has a polymorphic karyotype with a young neo-sex chromosome, resulting from a Robertsonian fusion between the X chromosome and an autosome, in part of its geographic range. We identified the shared and neo-sex chromosomes using comparative genetic maps of the two cytotypes. We found that sex-linked regions of both the ancestral and the neo-sex chromosomes are embedded in large regions of low recombination. Furthermore, our comparison of the recombination landscape of the neo-sex chromosome to its autosomal homolog indicates that low recombination rates mainly preceded sex linkage. These patterns are not unique to the sex chromosomes; all chromosomes were characterized by massive regions of suppressed recombination spanning most of each chromosome. This represents an extreme case of the periphery-biased recombination seen in other systems with large chromosomes. Across all chromosomes, gene and repetitive sequence density correlated with recombination rate, with patterns of variation differing by repetitive element type. Our findings suggest that ancestrally low rates of recombination may facilitate the formation and subsequent evolution of heteromorphic sex chromosomes.

Global analysis of floral longevity reveals latitudinal gradients and biotic and abiotic correlates.

The length of time a flower remains open and functional - floral longevity - governs important reproductive processes influencing pollination and mating and varies considerably among angiosperm species. However, little is known about large-scale biogeographic patterns and the correlates of floral longevity. Using published data on floral longevity from 818 angiosperm species in 134 families and 472 locations world-wide, we present the first global quantification of the latitudinal pattern of floral longevity and the relationships between floral longevity and a range of biotic and abiotic factors. Floral longevity exhibited a significant phylogenetic signal and was longer at higher latitudes in both northern and southern hemispheres, even after accounting for elevation. This latitudinal variation was associated with several biotic and abiotic variables. The mean temperature of the flowering season had the highest predictive power for floral longevity, followed by pollen number per flower. Surprisingly, compatibility status, flower size, pollination mode, and growth form had no significant effects on flower longevity. Our results suggest that physiological processes associated with floral maintenance play a key role in explaining latitudinal variation in floral longevity across global ecosystems, with potential implications for floral longevity under global climate change and species distributions.

Genetic variation and clonal diversity in floating aquatic plants: Comparative genomic analysis of water hyacinth species in their native range.

Many eukaryotic organisms reproduce by sexual and asexual reproduction. Genetic diversity in populations can be strongly dependent on the relative importance of these two reproductive modes. Here, we compare the amounts and patterns of genetic diversity in related water hyacinths that differ in their propensity for clonal propagation - highly clonal Eichhornia crassipes and moderately clonal E. azurea (Pontederiaceae). Our comparisons involved genotype-by-sequencing (GBS) of 137 E. crassipes ramets from 60 locations (193,495 nucleotide sites) and 118 E. azurea ramets from 53 locations (198,343 nucleotide sites) among six hydrological basins in central South America, the native range of both species. We predicted that because of more prolific clonal propagation, E. crassipes would exhibit lower clonal diversity than E. azurea. This prediction was supported by all measures of clonal diversity that we examined. Eichhornia crassipes also had a larger excess of heterozygotes at variant sites, another signature of clonality. However, genome-wide heterozygosity was not significantly different between the species. Eichhornia crassipes had weaker spatial genetic structure and lower levels of differentiation among hydrological basins than E. azurea, probably because of higher clonality and more extensive dispersal of its free-floating life form. Our findings for E. crassipes contrast with earlier studies from the invasive range which have reported very low levels of clonal diversity and extensive geographic areas of genetic uniformity.

Effects of the neo-X chromosome on genomic signatures of hybridization in Rumex hastatulus.

Natural hybrid zones provide opportunities for studies of the evolution of reproductive isolation in wild populations. Although recent investigations have found that the formation of neo-sex chromosomes is associated with reproductive isolation, the mechanisms remain unclear in most cases. Here, we assess the contemporary structure of gene flow in the contact zone between largely allopatric cytotypes of the dioecious plant Rumex hastatulus, a species with evidence of sex chromosome turn-over. Males to the west of the Mississippi river, USA, have an X and a single Y chromosome, whereas populations to the east of the river have undergone a chromosomal rearrangement giving rise to a larger X and two Y chromosomes. Using reduced-representation sequencing, we provide evidence that hybrids form readily and survive multiple backcross generations in the field, demonstrating the potential for ongoing gene flow between the cytotypes. Cline analysis of each chromosome separately captured no signals of difference in cline shape between chromosomes. However, principal component regression revealed a significant increase in the contribution of individual SNPs to inter-cytotype differentiation on the neo-X chromosome, but no correlation with recombination rate. Cline analysis revealed that the only SNPs with significantly steeper clines than the genome average were located on the neo-X. Our data are consistent with a role for neo-sex chromosomes in reproductive isolation between R. hastatulus cytotypes. Our investigation highlights the importance of studying plant hybrid zones for understanding the evolution of sex chromosomes.

Contrasting patterns of genetic diversity and differentiation across the continental disjunct range of a sexually polymorphic aquatic plant.

BACKGROUND AND AIMS: Reproductive systems enabling opportunities for self-fertilization influence population genetic structure and play a key role in colonization and genetic differentiation during range expansion. Because of their well-developed powers of dispersal, aquatic plants often have widespread disjunct geographical distributions, providing opportunities to investigate the role of reproductive systems in structuring genetic variation between parts of the range that differ in migration history and ecology. METHODS: We compared reproductive systems and spatial genetic structure of the freshwater aquatic macrophyte Sagittaria latifolia between disjunct western and eastern ranges of North America (NA). Populations of this species are most commonly either monoecious or dioecious. We examined chloroplast DNA haplotype variation to test the hypothesis that the western range of this species represents a secondary colonization from the east, and evaluated the roles of reproductive system differences and geography in structuring contemporary patterns of genetic variation at 11 polymorphic SSR (simple sequence repeat) loci. KEY RESULTS: Chloroplast haplotyping revealed a single haplotype in western NA compared to numerous haplotypes in eastern NA, consistent with a genetic bottleneck during westward migration. Estimates of genetic diversity in eastern NA populations differed significantly between reproductive systems, but this pattern was not evident in the western range. Eastern populations could be reliably assigned to genetic clusters based on their reproductive systems, whereas western populations clustered primarily by geographical location. CONCLUSIONS: The sparser distribution of aquatic habitats in the drier western range of S. latifolia, combined with secondary colonization of this region, probably cause the lower genetic diversity and increased differentiation among populations, possibly overriding the effects of reproductive system evident in the eastern portion of the range. Our findings demonstrate that the complex interplay between migratory history, reproductive systems and habitat availability plays an important role in structuring spatial patterns of genetic variation in disjunct plant populations.

Do annual and perennial populations of an insect-pollinated plant species differ in mating system?

BACKGROUND AND AIMS: Theory predicts that outcrossing should be more prevalent among perennials than annuals, a pattern confirmed by comparative evidence from diverse angiosperm families. However, intraspecific comparisons between annual and perennial populations are few because such variation is uncommon among flowering plants. Here, we test the hypothesis that perennial populations outcross more than annual populations by investigating Incarvillea sinensis, a wide-ranging insect-pollinated herb native to China. The occurrence of both allopatric and sympatric populations allows us to examine the stability of mating system differences between life histories under varying ecological conditions. METHODS: We estimated outcrossing rates and biparental inbreeding in 16 allopatric and five sympatric populations in which both life histories coexisted using 20 microsatellite loci. In each population we measured height, branch number, corolla size, tube length and herkogamy for ~30 individuals. In a sympatric population, we recorded daily flower number, pollinator visitation and the fruit and seed set of annual and perennial plants. KEY RESULTS: As predicted, outcrossing rates (t) were considerably higher in perennial (mean = 0.76) than annual (mean = 0.09) populations. This difference in mating system was also maintained at sympatric sites where plants grew intermixed. In both allopatric and sympatric populations the degree of herkogamy was consistently larger in outcrossing than selfing plants. Perennials were more branched, with more and larger flowers than in annuals. In a sympatric population, annuals had a significantly higher fruit and seed set than perennials. CONCLUSIONS: Genetically based differences in herkogamy between annuals and perennials appear to play a key role in governing outcrossing rates in populations, regardless of variation in local ecological conditions. The maintenance of mating system and life history trait differentiation between perennial and annual populations of I. sinensis probably results from correlated evolution in response to local environmental conditions.

Influence of local density and sex ratio on pollination in an ambophilous flowering plant.

PREMISE: Variation in local density and sex ratio in dioecious plants can affect mating success through the actions of pollen vectors, principally generalist insects or wind. Increased density and male-biased sex ratios should promote pollen transfer and seed production, but their combined effects have not been investigated for ambophilous species, which exhibit both insect and wind pollination. METHODS: We manipulated density (low vs. high) and sex ratio (1:1 vs. 3:1 male-biased) in arrays of dioecious ambophilous Thalictrum pubescens. We quantified visitation rates and foraging times to examine whether pollinators exhibited sex-specific preferences and determined the seed set of arrays. RESULTS: Pollinators visited more plants per foraging bout at high than low density. Visitation rates and foraging times of visitors were greater for male than for female plants but did not depend on the density or sex ratio of arrays. However, whereas solitary bees displayed a strong preference for males, hover flies were indifferent to plant sex phenotype. Solitary bees also visited significantly more plants per foraging bout than hover flies. There was a significant interaction between density and sex ratio on seed set. At low density, seed set was greater for 3:1 than for 1:1 arrays, but at high density the opposite pattern occurred. CONCLUSIONS: The demographic factors we investigated had complex influences on pollinator foraging behavior and patterns of seed set. Several factors may explain our results, including the influence of density and sex ratio on pollen export from arrays, grooming by pollinators, and the contribution of wind pollination.

'A most complex marriage arrangement': recent advances on heterostyly and unresolved questions.

Heterostylous genetic polymorphisms provide paradigmatic systems for investigating adaptation and natural selection. Populations are usually comprised of two (distyly) or three (tristyly) mating types, maintained by negative frequency-dependent selection resulting from disassortative mating. Theory predicts this mating system should result in equal style-morph ratios (isoplethy) at equilibrium. Here, I review recent advances on heterostyly, focusing on examples challenging stereotypical depictions of the polymorphism and unresolved questions. Comparative analyses indicate multiple origins of heterostyly, often within lineages. Ecological studies demonstrate that structural components of heterostyly are adaptations improving the proficiency of animal-mediated cross-pollination and reducing pollen wastage. Both neutral and selective processes cause deviations from isoplethy in heterostylous populations, and, under some ecological and demographic conditions, cause breakdown of the polymorphism, resulting in either the evolution of autogamy and mixed mating, or transitions to alternative outcrossing systems, including dioecy. Earlier ideas on the genetic architecture of the S-locus supergene governing distyly have recently been overturned by discovery that the dominant S-haplotype is a hemizygous region absent from the s-haplotype. Ecological, phylogenetic and molecular genetic data have validated some features of theoretical models on the selection of the polymorphism. Although heterostyly is the best-understood floral polymorphism in angiosperms, many unanswered questions remain.

Comparative analysis of pollen release biomechanics in Thalictrum: implications for evolutionary transitions between animal and wind pollination.

Transitions from animal to wind pollination have occurred repeatedly in flowering plants, driven by structural and biomechanical modifications to flowers. But the initial changes promoting wind pollination are poorly understood, especially those required to release pollen into airflows - the critical first stage of wind pollination. Using a wind tunnel, we performed a comparative study of pollen release biomechanics in 36 species of animal- and wind-pollinated Thalictrum. We quantified pollination syndromes and stamen natural frequency (fn ), a key vibration parameter, to determine if floral traits reliably predicted pollen release probability. We then investigated if pollen release was caused by wind-induced resonance vibration of stamens. We detected wind-induced stamen resonance in 91% of species and a strong effect of stamen acceleration on pollen release, inversely driven by fn . However, unlike fn , pollination syndromes did not reliably predict the probability of pollen release among species. Our results directly link fn to the capacity of stamens to release pollen by wind and suggest that structural mechanisms reducing fn are likely to be important for initiating transitions from animal to wind pollination. Our inability to predict the probability of pollen release based on pollination syndromes suggests diverse phenotypic trajectories from animal to wind pollination.

Variation in sexual dimorphism in a wind-pollinated plant: the influence of geographical context and life-cycle dynamics.

Understanding the mechanisms causing phenotypic differences between females and males has long fascinated evolutionary biologists. An extensive literature exists on animal sexual dimorphism but less information is known about sex differences in plants, particularly the extent of geographical variation in sexual dimorphism and its life-cycle dynamics. Here, we investigated patterns of genetically based sexual dimorphism in vegetative and reproductive traits of a wind-pollinated dioecious plant, Rumex hastatulus, across three life-cycle stages using open-pollinated families from 30 populations spanning the geographic range and chromosomal variation (XY and XY1 Y2 ) of the species. The direction and degree of sexual dimorphism was highly variable among populations and life-cycle stages. Sex-specific differences in reproductive function explained a significant amount of temporal change in sexual dimorphism. For several traits, geographical variation in sexual dimorphism was associated with bioclimatic parameters, likely due to the differential responses of the sexes to climate. We found no systematic differences in sexual dimorphism between chromosome races. Sex-specific trait differences in dioecious plants largely result from a balance between sexual and natural selection on resource allocation. Our results indicate that abiotic factors associated with geographical context also play a role in modifying sexual dimorphism during the plant life-cycle.

Evolutionary history of the buildup and breakdown of the heterostylous syndrome in Plumbaginaceae.

The evolutionary pathways leading to the heterostylous syndrome are not well understood, and models concerning the origins of distyly differ in the order in which reciprocal herkogamy and self-incompatibility evolve. We investigated the evolution and breakdown of distyly in Plumbaginaceae, a family with considerable diversity of floral traits and reproductive systems. Using Bayesian Markov chain Monte Carlo analyses and stochastic character mapping, we examined the evolutionary assembly and breakdown of the heterostylous syndrome based on a well-resolved phylogeny of 121 species of Plumbaginaceae and six outgroup taxa using five nuclear and plastid gene regions. We used the distribution of reproductive traits and reconstructed ancestral characters across phylogenies to evaluate competing models for the evolution of distyly. The most likely common ancestor of Plumbaginaceae was self-incompatible and monomorphic for sex-organ arrangement and pollen-stigma characters. Character state reconstructions indicated that reciprocal herkogamy evolved at least three times and that shifts to selfing and apomixis occurred on multiple occasions. Our results provide comparative support for the early ideas of H. G. Baker on evolutionary pathways in Plumbaginaceae, and the more recent selfing avoidance model by D. & B. Charlesworth in which distyly evolves from self-incompatible ancestors.

Water mediates fertilization in a terrestrial flowering plant.

Water-mediated fertilization is ubiquitous in early land plants. This ancestral mode of fertilization has, however, generally been considered to have been lost during the evolutionary history of terrestrial flowering plants. We investigated reproductive mechanisms in the subtropical ginger Cautleya gracilis (Zingiberaceae), which has two pollen conditions - granular and filiform masses - depending on external conditions. We tested whether rain transformed granular pollen into filiform masses and whether this then promoted pollen-tube growth and fertilization of ovules. Using experimental manipulations in the field we investigated the contribution of water-mediated fertilization to seed production. Rain caused granular pollen to form filiform masses of germinating pollen tubes, which transported sperm to ovules, resulting in fertilization and seed set. Flowers exposed to rain produced significantly more seeds than those protected from the rain, which retained granular pollen. Insect pollination made only a limited contribution to seed set because rainy conditions limited pollinator service. Our results reveal a previously undescribed fertilization mechanism in flowering plants involving water-mediated fertilization stimulated by rain. Water-mediated fertilization is likely to be adaptive in the subtropical monsoon environments in which C. gracilis occurs by ensuring reproductive assurance when persistent rain prevents insect-mediated pollination.

Phylogenomic analysis reveals multiple evolutionary origins of selfing from outcrossing in a lineage of heterostylous plants.

Evolutionary transitions from outcrossing to selfing often occur in heterostylous plants. Selfing homostyles originate within distylous populations and frequently evolve to become reproductively isolated species. We investigated this process in 10 species of Primula section Obconicolisteri using phylogenomic approaches and inferred how often homostyly originated from distyly and its consequences for population genetic diversity and floral trait evolution. We estimated phylogenetic relationships and reconstructed character evolution using the whole plastome comprised of 76 protein-coding genes. To investigate mating patterns and genetic diversity we screened 15 microsatellite loci in 40 populations. We compared floral traits among distylous and homostylous populations to determine how phenotypically differentiated homostyles were from their distylous ancestors. Section Obconicolisteri was monophyletic and we estimated multiple independent transitions from distyly to homostyly. High selfing rates characterised homostylous populations and this was associated with reduced genetic diversity. Flower size and pollen production were reduced in homostylous populations, but pollen size was significantly larger in some homostyles than in distylous morphs. Repeated transitions to selfing in section Obconicolisteri are likely to have been fostered by the complex montane environments that species occupy. Unsatisfactory pollinator service is likely to have promoted reproductive assurance in homostyles leading to subsequent population divergence through isolation.