Pollen success, functional gender and assortative mating in an experimental plant population

Antagonistic Effects of Assortative Mating on the Evolution of Phenotypic Plasticity along Environmental Gradients

AbstractPrevious theory has shown that assortative mating for plastic traits can maintain genetic divergence across environmental gradients despite high gene flow. Yet these models did not examine how assortative mating affects the evolution of plasticity. We here describe patterns of genetic variation across elevation for plasticity in a trait under assortative mating, using multiple-year observations of budburst date in a common garden of sessile oaks. Despite high gene flow, we found significant spatial genetic divergence for the intercept, but not for the slope, of reaction norms to temperature. We then used individual-based simulations, where both the slope and the intercept of the reaction norm evolve, to examine how assortative mating affects the evolution of plasticity, varying the intensity and distance of gene flow. Our model predicts the evolution of either suboptimal plasticity (reaction norms with a slope shallower than optimal) or hyperplasticity (slopes steeper than optimal) in the presence of assortative mating when optimal plasticity would evolve under random mating. Furthermore, a cogradient pattern of genetic divergence for the intercept of the reaction norm (where plastic and genetic effects are in the same direction) always evolves in simulations with assortative mating, consistent with our observations in the studied oak populations.

Sequenced-based paternity analysis to improve breeding and identify self-incompatibility loci in intermediate wheatgrass (Thinopyrum intermedium)

KEY MESSAGE

Paternity assignment and genome-wide association analyses for fertility were applied to a Thinopyrum intermedium breeding program. A lack of progeny between combinations of parents was associated with loci near self-incompatibility genes. In outcrossing species such as intermediate wheatgrass (IWG, Thinopyrum intermedium), polycrossing is often used to generate novel recombinants through each cycle of selection, but it cannot track pollen-parent pedigrees and it is unknown how self-incompatibility (SI) genes may limit the number of unique crosses obtained. This study investigated the potential of using next-generation sequencing to assign paternity and identify putative SI loci in IWG. Using a reference population of 380 individuals made from controlled crosses of 64 parents, paternity was assigned with 92% agreement using Cervus software. Using this approach, 80% of 4158 progeny (n = 3342) from a polycross of 89 parents were assigned paternity. Of the 89 pollen parents, 82 (92%) were represented with 1633 unique full-sib families representing 42% of all potential crosses. The number of progeny per successful pollen parent ranged from 1 to 123, with number of inflorescences per pollen parent significantly correlated to the number of progeny (r = 0.54, p < 0.001). Shannon's diversity index, assessing the total number and representation of families, was 7.33 compared to a theoretical maximum of 8.98. To test our hypothesis on the impact of SI genes, a genome-wide association study of the number of progeny observed from the 89 parents identified genetic effects related to non-random mating, including marker loci located near putative SI genes. Paternity testing of polycross progeny can impact future breeding gains by being incorporated in breeding programs to optimize polycross methodology, maintain genetic diversity, and reveal genetic architecture of mating patterns.

An analysis of mating biases in trees

Assortative mating is a deviation from random mating based on phenotypic similarity. As it is much better studied in animals than in plants, we investigate for trees whether kinship of realized mating pairs deviates from what is expected from the set of potential mates and use this information to infer mating biases that may result from kin recognition and/or assortative mating. Our analysis covers 20 species of trees for which microsatellite data is available for adult populations (potential mates) as well as seed arrays. We test whether mean relatedness of observed mating pairs deviates from null expectations that only take pollen dispersal distances into account (estimated from the same data set). This allows the identification of elevated as well as reduced kinship among realized mating pairs, indicative of positive and negative assortative mating, respectively. The test is also able to distinguish elevated biparental inbreeding that occurs solely as a result of related pairs growing closer to each other from further assortativeness. Assortative mating in trees appears potentially common but not ubiquitous: nine data sets show mating bias with elevated inbreeding, nine do not deviate significantly from the null expectation, and two show mating bias with reduced inbreeding. While our data sets lack direct information on phenology, our investigation of the phenological literature for each species identifies flowering phenology as a potential driver of positive assortative mating (leading to elevated inbreeding) in trees. Since active kin recognition provides an alternative hypothesis for these patterns, we encourage further investigations on the processes and traits that influence mating patterns in trees.

Temporal population genetic structure in the pollen pool for flowering time: A field experiment with Brassica rapa (Brassicaceae)

PREMISE OF THE STUDY

Assortative mating by flowering time can cause temporal genetic structure in species with heritable flowering times. A strong temporal structure, when coupled with a seasonal shift in selection, may lead to adaptive temporal clines. We implemented a prospective and retrospective method to estimate the temporal genetic structure in the pollen pool of Brassica rapa.

METHODS

The prospective method uses flowering schedules to estimate the seasonal shift in the pollen donors' phenotype. By examining the offspring generation, we can get a direct estimate of temporal genetic structure, i.e., a retrospective estimate. However, this estimate is problematic because of the phenotypic correlation of the trait of interest, flowering time, between dam and sire. We developed a novel retrospective method that isolates flowering time by holding the maternal contribution constant and sampled the pollen pool in eight open-pollinated field plots throughout the flowering season.

KEY RESULTS

We found temporal genetic structure for flowering time in seven of the eight field plots. Interestingly, the direct (retrospective) temporal structure estimate was 35% larger than the prospective estimate based on flowering schedules. Spatial clumping of pollen donors did not affect temporal structure, but structure intensified when heritability was experimentally enhanced.

CONCLUSIONS

Temporal genetic structure, especially for flowering time, likely occurs in many plant populations and may be underestimated using a prospective method. We discuss the genome-wide consequences of temporal genetic structure and the potential for adaptive temporal clines in plant populations.

THE EFFECT OF INFLORESCENCE SIZE ON MALE FITNESS: EXPERIMENTAL TESTS IN THE ANDROMONOECIOUS LILY, ZIGADENUS PANICULATUS

We studied the relationship between inflorescence size and male fitness in the andromonoecious lily Zigadenus paniculatus, using experimentally manipulated inflorescences to eliminate possible correlations between flower number, resource availability, and other floral traits. Allozyme markers were used to determine the siring success of large versus small plants in 14 arrays of plants, each array containing five large and five small plants. The inflorescence size of small plants was held constant both within and among arrays; the size of large plants was held constant within an array but was varied among arrays. Large plants sired more than half the seeds in 12 of the 14 arrays, and significantly more than half in six of these 12. However, in eight of the arrays, large plants sired significantly fewer seeds than expected on the basis of their size advantage. Furthermore, there was no significant relationship between relative size and relative siring success in comparisons among arrays. A maximum-likelihood model estimated that 28% of seeds were sired by imported pollen, with 95% confidence limits of 13% and 50%. Within these limits, high import rates tended to mask the relative success of large plants in several arrays. These results suggest that the evolution of inflorescence size in Z. paniculatus is at least partly driven by selection for increased male success, assuming genetic variation for flower number. However, the data also support a growing body of evidence that estimates of male fitness in plants can be highly variable. We discuss the sources of this variability and the possible effects of inflorescence design on the relationship between inflorescence size and fitness.

EXPERIMENTAL STUDIES ON THE FUNCTIONAL SIGNIFICANCE OF HETEROSTYLY

Heterostyly has been viewed as both an antiselfing device and a mechanism that increases the proficiency of pollen transfer between plants. We used experimental manipulation of the morph structure of garden populations of self-compatible, tristylous Eichhornia paniculata to investigate the function of floral polymorphism. Outcrossing rates (t), levels of intermorph mating (d), and morph-specific male and female reproductive success were compared in replicate trimorphic and monomorphic populations. In trimorphic populations, t and d averaged 0.81 (2 SE = 0.03) and 0.77 (2 SE = 0.03) respectively, with no difference in either parameter among morphs. Ninety-five percent of outcrossed seeds were therefore the result of intermorph fertilizations. Male reproductive success of the long-styled morph was low, especially in comparison with plants of the short-styled morph. Outcrossing rates for each morph were higher in trimorphic than monomorphic populations where t averaged 0.71 (2 SE = 0.01), 0.30 (2 SE = 0.04) and 0.43 (2 SE = 0.1) for the long-, mid-, and short-styled morphs, respectively. Seed set was lower in monomorphic populations, particularly those composed of the L morph, reflecting reduced pollen deposition. Floral polymorphism therefore increased both outcrossing rate and fecundity but the magnitude of the differences varied among morphs. If the ancestral condition in heterostylous groups resembled the L morph, as has been suggested, data from this study suggests that the selective basis for the establishment of floral polymorphism could have been increased pollen transfer rather than higher levels of outcrossing.

When can stress facilitate divergence by altering time to flowering?

Stressors and heterogeneity are ubiquitous features of natural environments, and theory suggests that when environmental qualities alter flowering schedules through phenotypic plasticity, assortative mating can result that promotes evolutionary divergence. Therefore, it is important to determine whether common ecological stressors induce similar changes in flowering time. We review previous studies to determine whether two important stressors, water restriction and herbivory, induce consistent flowering time responses among species; for example, how often do water restriction and herbivory both delay flowering? We focus on the direction of change in flowering time, which affects the potential for divergence in heterogeneous environments. We also tested whether these stressors influenced time to flowering and nonphenology traits using Mimulus guttatus. The literature review suggests that water restriction has variable effects on flowering time, whereas herbivory delays flowering with exceptional consistency. In the Mimulus experiment, low water and herbivory advanced and delayed flowering, respectively. Overall, our results temper theoretical predictions for evolutionary divergence due to habitat‐induced changes in flowering time; in particular, we discuss how accounting for variation in the direction of change in flowering time can either increase or decrease the potential for divergence. In addition, we caution against adaptive interpretations of stress‐induced phenology shifts.

Sexually antagonistic polymorphism in simultaneous hermaphrodites

In hermaphrodites, pleiotropic genetic trade-offs between female and male reproductive functions can lead to sexually antagonistic (SA) selection, where individual alleles have conflicting fitness effects on each sex function. Although an extensive theory of SA selection exists for dioecious species, these results have not been generalized to hermaphrodites. We develop population genetic models of SA selection in simultaneous hermaphrodites, and evaluate effects of dominance, selection on each sex function, self-fertilization, and population size on the maintenance of polymorphism. Under obligate outcrossing, hermaphrodite model predictions converge exactly with those of dioecious populations. Self-fertilization in hermaphrodites generates three points of divergence with dioecious theory. First, opportunities for stable polymorphism decline sharply and become less sensitive to dominance with increased selfing. Second, selfing introduces an asymmetry in the relative importance of selection through male versus female reproductive functions, expands the parameter space favorable for the evolutionary invasion of female-beneficial alleles, and restricts invasion criteria for male-beneficial alleles. Finally, contrary to models of unconditionally beneficial alleles, selfing decreases genetic hitchhiking effects of invading SA alleles, and should therefore decrease these population genetic signals of SA polymorphisms. We discuss implications of SA selection in hermaphrodites, including its potential role in the evolution of "selfing syndromes."

Phenological patterns in a natural population of a tropical timber tree species, Milicia excelsa (Moraceae): Evidence of isolation by time and its interaction with feeding strategies of dispersers

PREMISE OF THE STUDY

Population genetic structuring over limited timescales is commonly viewed as a consequence of spatial constraints. Indirect approaches have recently revealed reproductive isolation resulting from flowering time (so-called isolation by time, IBT). Since phenological processes can be subject to selection, the persistence of flowering asynchrony may be due to opposing selective pressures during mating, dispersal, and regeneration phases. Our study aimed to investigate phenology, fruit handling by animals, and their interaction in a timber tree species, Milicia excelsa.

METHODS

We analyzed phenological data collected over 6 years on 69 genotyped trees in a Cameroonian natural rainforest complemented by data from germination trials and field observations of dispersers.

KEY RESULTS

Initiation of flowering was correlated with variation in temperature and relative humidity, but was also affected by genetic factors: pairwise differences in flowering time between nearby individuals correlated with kinship coefficient, and earliness of flowering remained stable over time. A decrease in mean seed production per fruit with increasing flowering time suggests selection against late bloomers. However, germination rate was not affected by seed collection date, and the main seed disperser, the bat Eidolon helvum, seemed to increase in abundance at the end of the reproductive season and preferred trees in open habitats where early and late bloomers are expected.

CONCLUSIONS

The pairwise approach performs well in detecting IBT. The persistence of different mating pools in such a case may result from a trade off between selective forces during the mating and seed dispersal processes.

Population structure attributable to reproductive time: isolation by time and adaptation by time

Many populations are composed of a mixture of individuals that reproduce at different times, and these times are often heritable. Under these conditions, gene flow should be limited between early and late reproducers, even within populations having a unimodal temporal distribution of reproductive activity. This temporal restriction on gene flow might be called "isolation by time" (IBT) to acknowledge its analogy with isolation by distance (IBD). IBD and IBT are not exactly equivalent, however, owing to differences between dispersal in space and dispersal in time. We review empirical studies of natural populations that provide evidence for IBT based on heritabilities of reproductive time and on molecular genetic differences associated with reproductive time. When IBT is present, variation in selection through the reproductive season may lead to adaptive temporal variation in phenotypic traits [adaptation by time (ABT)]. We introduce a novel theoretical model that shows how ABT increases as (i) selection on the trait increases; (ii) environmental influences on reproductive time decrease; (iii) the heritability of reproductive time increases; and (iv) the temporal distribution of reproductive activity becomes increasingly uniform. We then review empirical studies of natural populations that provide evidence for ABT by documenting adaptive temporal clines in phenotypic traits. The best evidence for IBT and ABT currently comes from salmonid fishes and flowering plants, but we expect that future work will show these processes are more widespread.

Direct and indirect assortative mating: a multivariate approach to plant flowering schedules

This paper develops methods to partition the phenotypic correlation between mates for a focal trait--the standard measure for assortative mating--into a direct component and additional indirect components. Indirect assortative mating occurs when a nonassorting trait is correlated within individuals to a directly assorting trait. Direct and indirect assortative mating is assessed for flowering phenology in Brassica rapa. The flowering time of pollen recipients (mothers) was strongly correlated (rho=0.67) to that of potential pollen donors (fathers). Similarly, recipients and donors were correlated for duration of their flowering periods (rho=0.32) and stem diameters (rho=0.52). A partitioning of between-mate correlations revealed direct assortative mating for flowering time and period duration. However, assortment for stem diameter is explained solely through its correlation to flowering time. Examination of standard quantitative genetic theory shows that indirect assortative mating inflates genetic variance in a focal trait and the genetic covariance between focal and phenotypically correlated traits.

Genetic variation in flowering time induces phenological assortative mating: quantitative genetic methods applied to Brassica rapa

It has been argued from first principles that plants mate assortatively by flowering time. However, there have been very few studies of phenological assortative mating, perhaps because current methods to infer paternal phenotype are difficult to apply to natural populations. Two methods are presented to estimate the phenotypic correlation between mates-the quantitative genetic metric for assortative mating-for phenological traits. The first method uses individual flowering schedules to estimate mating probabilities for every potential pairing in a sample. These probabilities are then incorporated into a weighted phenotypic correlation between all potential mates and thus yield a prospective estimate based on mating opportunities. The correlation between mates can also be estimated retrospectively by comparing the regression of offspring phenotype over one parent, which is inflated by assortative mating, to the regression over mid-parent, which is not. In a demonstration experiment with Brassica rapa, the prospective correlation between flowering times (days from germination to anthesis) of pollen recipients and their potential donors was 0.58. The retrospective estimate of this correlation strongly agreed with the prospective estimate. The prospective method is easily employed in field studies that explore the effect of phenological assortative mating on selection response and population differentiation.

Current knowledge of gene flow in plants: implications for transgene flow

Plant evolutionary biologists' view of gene flow and hybridization has undergone a revolution. Twenty-five years ago, both were considered rare and largely inconsequential. Now gene flow and hybridization are known to be idiosyncratic, varying with the specific populations involved. Gene flow typically occurs at evolutionarily significant rates and at significant distances. Spontaneous hybridization occasionally has important applied consequences, such as stimulating the evolution of more aggressive invasives and increasing the extinction risk for rare species. The same problems have occurred for spontaneous hybridization between crops and their wild relatives. These new data have implications for transgenic crops: (i) for most crops, gene flow can act to introduce engineered genes into wild populations; (ii) depending on the specific engineered gene(s) and populations involved, gene flow may have the same negative impacts as those observed for traditionally improved crops; (iii) gene flow's idiosyncratic nature may frustrate management and monitoring attempts; and (iv) intercrop transgene flow, although rarely discussed, is equally worthy of study.

Paternity and gene dispersal in limber pine (Pinus flexilis James)

This study provides empirical information on intrapopulation gene dispersal via pollen, the size of genetic 'neighbourhoods', and interpopulation gene flow in a long-lived conifer, limber pine (Pinus flexilis). We used allozyme loci for a paternity analysis of 518 seeds produced in an isolated population of limber pine located in north-eastern Colorado, U.S.A., separated by 2 km from the nearest conspecific trees and nearly 100 km from populations in the Rocky Mountains. We also used indirect techniques (FST analyses) to estimate gene flow rates among subdivisions of the study population and among five widely separated populations. Within the main study population limber pine exhibited a polymorphism level of 50%, observed heterozygosity of 0.159, and 2.36 alleles per polymorphic locus. Mountain populations were slightly more variable. The main study population showed significant differentiation in allozymes among neighbouring subpopulations. The mean FST was 0.031 and the gene flow rate among subpopulations was estimated as 7.8 migrants per generation. Among widely separated populations the mean FST was 0.035 and the gene flow rate was estimated as 6.9 migrants per generation. The paternity analysis indicated a best estimate of 6.5% pollen immigration (minimum 1.1%) from populations 2 km to 100+ km away. For 4% of the seeds examined, paternity could be ascribed to a single tree in the study population. Fractional paternity and likelihood methods were used to estimate pollen dispersal distances for the remainder of the seeds. Mean pollen dispersal distance was estimated at 140 m using the fractional method, similar to results from the other techniques. This compares with a mean distance of 172 m between potential mates. These results suggest near-panmictic pollen dispersal over this population, which covers about 15 ha. The observed allozyme differences and surprisingly low estimates of among-subpopulation gene flow are ascribed to a probable restriction of gene dispersal by seeds.

Fractional paternity assignment: theoretical development and comparison to other methods

There has recently been a burgeoning interest in the analysis of paternity patterns for natural populations because of its relevance to population genetic phenomena such as the distance between successful mates, relative male reproductive success and gene flow. In this paper we develop a method of analyzing populational patterns of paternity, the fractional paternity method, and compare its performance to two other commonly used methods of paternity analysis (simple exclusion and the most-likely methods). We show that the fractional method is the most accurate method for determining populational patterns of paternity because it assigns paternity to all progeny examined, and because it avoids biases inherent in the other paternity analysis methods when model assumptions are met. In particular, it avoids a systematic bias of the most-likely paternity assignment method, which has a tendency to over-assign paternity of progeny to certain male parents with a greater than average number of homozygous marker loci. We also demonstrate the effect of linkage of some of the marker loci on paternity assignment, showing how the knowledge of the linkage phase of male and female parents in the population can significantly improve the accuracy of the estimates of populational patterns of paternity. Knowledge of the linkage phase of individuals in a population is usually unknown and difficult to assess without progeny testing, which involves considerable labor. However, we show how the linkage phase of hermaphroditic individuals in a population can be obtained in conjunction with the paternity analysis if progeny can be obtained from each hermaphroditic individual in the population, thereby avoiding the problem of traditional progeny testing. Applications of the fractional paternity approach developed herein should contribute significantly to our understanding of the mating patterns in, and hence the evolution of, natural populations.