SEX CHROMOSOME LINKED GENETIC VARIANCE AND THE EVOLUTION OF SEXUAL DIMORPHISM OF QUANTITATIVE TRAITS

Theory predicts that sex chromsome linkage should reduce intersexual genetic correlations thereby allowing the evolution of sexual dimorphism. Empirical evidence for sex linkage has come largely from crosses and few studies have examined how sexual dimorphism and sex linkage are related within outbred populations. Here, we use data on an array of different traits measured on over 10,000 individuals from two pedigreed populations of birds (collared flycatcher and zebra finch) to estimate the amount of sex‐linked genetic variance (h²_z). Of 17 traits examined, eight showed a nonzero h²_Z estimate but only four were significantly different from zero (wing patch size and tarsus length in collared flycatchers, wing length and beak color in zebra finches). We further tested how sexual dimorphism and the mode of selection operating on the trait relate to the proportion of sex‐linked genetic variance. Sexually selected traits did not show higher h²_Z than morphological traits and there was only a weak positive relationship between h²_Z and sexual dimorphism. However, given the relative scarcity of empirical studies, it is premature to make conclusions about the role of sex chromosome linkage in the evolution of sexual dimorphism.     

GENETIC,ENVIRONMENTAL, AND CONDITION-DEPENDENT EFFECTS ON FEMALE AND MALE ORNAMENTATION IN THE BARN OWL TYTO ALBA

Secondary sexual characters are thought to indicate individual quality. Expression of sex-limited traits in an extravagant state may require both the underlying genes and the available nutrient resources. The assessment of the relative contribution of genes, environment, and body condition is relevant for understanding to that extent the extravagant trait may signal genotypic or phenotypic quality of the individual. In birds, usually only the males are ornamented. In the barn owl, Tyto alba, both females and males display sex-limited plumage traits. Males are commonly lighter colored and females spottier. In an experiment with combined cross-fostering and brood size manipulation we determined the relative contribution of genes, environment, and body condition to the variation in plumage coloration and plumage spottiness. The partial cross-fostering experiment tested the relative importance of shared genes and a shared environment for the resemblance of related birds. Siblings raised in different nests converged toward similar trait values, offspring resembled the true but not the foster parents, and plumage traits of unrelated nestlings sharing the same nest were not correlated. Results were not inflated by maternal effects detectable in the mother's phenotype, because middaughter to mother resemblance was not higher than midson to father resemblance. This suggests that plumage coloration and spottiness are largely genetically inherited traits, and that the rearing environment does not have a strong impact on the expression of these traits. To further investigate whether the two sex-limited traits are condition dependent, brood sizes were manipulated. Enlargement or reduction of broods by two nestlings resulted in lower and higher body mass of nestlings, respectively. However, nestlings raised in enlarged or reduced broods did not show either a significantly darker or lighter or a more or less spotted plumage. We did not detect any genotype-by-environment interaction. In conclusion, simultaneous cross-fostering and brood size manipulation demonstrate that additive genetic variance for plumage coloration and spottiness is maintained and that both the rearing environment and body condition do not account for a large proportion of the phenotypic variance in female and male ornamentations.     

CHARACTERIZATION OF GENETIC MARKERS LINKED TO SEX DETERMINATION IN THE HAPLOID‐DIPLOID RED ALGA GRACILARIA CHILENSIS1

Bulk segregant analysis, random amplified polymorphic DNA (RAPD), and sequence characterized amplified region (SCAR) methods were used to identify sex‐linked molecular markers in the haploid‐diploid rhodophyte Gracilaria chilensis C. J. Bird, McLachlan et E. C. Oliveira. One hundred and eighty 10 bp primers were tested on three bulks of DNA: haploid males, haploid females, and diploid tetrasporophytes. Three RAPD primers (OPD15, OPG16, and OPN20) produced male‐specific bands; and one RAPD primer (OPD12), a female‐specific band. The sequences of the cloned putative sex‐specific PCR fragments were used to design specific primers for the female marker SCAR‐D12‐386 and the male marker SCAR‐G16‐486. Both SCAR markers gave unequivocal band patterns that allowed sex and phase to be determined in G. chilensis. Thus, all the females presented only the female band, and all the males only the male band, while all the tetrasporophytes amplified both male and female bands. Despite this sex‐specific association, we were able to amplify SCAR‐D12‐386 and SCAR‐G16‐486 in both sexes at low melting temperature. The differences between male and female sequences were of 8%–9% nucleotide divergence for SCAR‐D12‐386 and SCAR‐G16‐486, respectively. SCAR‐D12‐386 and SCAR‐G16‐486 could represent degenerated or diverged sequences located in the nonrecombining region of incipient sex chromosomes or heteromorphic sex chromosomes with sequence differences at the DNA level such that PCR primers amplify only one allele and not the other in highly specific PCR conditions. Seven gametic progenies composed of 19 males, 19 females, and the seven parental tetrasporophytes were analyzed. In all of them, the two SCAR markers segregated perfectly with sexual phenotypes.     

SEXUAL SELECTION CAN INCREASE THE EFFECT OF RANDOM GENETIC DRIFT—A QUANTITATIVE GENETIC MODEL OF POLYMORPHISM IN OOPHAGA PUMILIO,THE STRAWBERRY POISON‐DART FROG

The variation in color pattern between populations of the poison‐dart frog Oophaga pumilio across the Bocas del Toro archipelago in Panama is suggested to be due to sexual selection, as two other nonsexually selecting Dendrobatid species found in the same habitat and range do not exhibit this variation. We theoretically test this assertion using a quantitative genetic sexual selection model incorporating aposematic coloration and random drift. We find that sexual selection could cause the observed variation via a novel process we call “coupled drift.” Within our model, for certain parameter values, sexual selection forces frog color to closely follow the evolution of female preference. Any between‐population variation in preference due to genetic drift is passed on to color. If female preference in O. pumilio is strongly affected by drift, whereas color in the nonsexually selecting Dendrobatid species is not, coupled drift will cause increased between‐population phenotypic variation. However, with different parameter values, coupled drift will result in between‐population variation in color being suppressed compared to its neutral value, or in little or no effect. We suggest that coupled drift is a novel theoretical process that could have a role linking sexual selection with speciation both in O. pumilio, and perhaps more generally.     

NEUTRAL THEORY OF QUANTITATIVE GENETIC VARIANCE IN AN ISLAND MODEL WITH LOCAL EXTINCTION AND COLONIZATION

Additive genetic variance maintained by mutation in a selectively neutral quantitative character is analyzed for an ideal population distributed on n islands, each with local effective size N, that exchange migrants at a small rate, m. In a stable population structure, the expected genetic variance maintained within islands is identical to that in a panmictic population of the same total size, regardless of the migration rate (m > 0). This result contrasts with Wright's classical conclusion, based on inbreeding coefficients, that at least one immigrant per island every other generation (Nm > ½) is necessary for the genetic variance within local populations to approach that under panmixia. The expected genetic variance maintained among islands is inversely proportional to m and increases with the number of islands, but is independent of N. Local extinction and colonization diminish the genetic variance maintained within islands by reducing the effective size of island populations through the founder effect, although the expected genetic variance within islands is nearly as large as that in a panmictic population of the same total effective size. If the founders of new colonies originate from more than one island, rates of local extinction and colonization larger than about twice the migration rate will substantially reduce the genetic variance maintained among islands. These results indicate the importance of mutation and migration in maintaining quantitative genetic variance within small local populations.     

THE EVOLUTION OF THRESHOLD TRAITS: A QUANTITATIVE GENETIC ANALYSIS OF THE PHYSIOLOGICAL AND LIFE-HISTORY CORRELATES OF WING DIMORPHISM IN THE SAND CRICKET

Many traits are phenotypically discrete but polygenically determined. Such traits can be understood using the threshold model of quantitative genetics that posits a continuously distributed underlying trait, called the liability, and a threshold of response, individuals above the threshold displaying one morph and individuals below the threshold displaying the alternate morph. For many threshold traits the liability probably consists of a hormone or a suite of hormones. Previous experiments have implicated juvenile hormone esterase (JHE), a degratory enzyme of juvenile hormone, as a physiological determinant of wing dimorphism in the crickets Gryllus rubens and G. firmus. The present study uses a half-sib experiment to measure the heritability of JHE in the last nymphal stadium of G. firmus and its genetic correlation with fecundity, a trait that is itself genetically correlated with wing morph. The phenotypic and genetic parameters are consistent with the hypothesis that JHE is a significant component of the liability. Comparison of sire and dam estimates suggest that nonadditive effects may be important. Two models have been proposed to account for the fitness differences between morphs: the dichotomy model, which assumes that each morph can be characterized by a particular suite of traits, and the continuous model, which assumes that the associated fitness traits are correlated with the liability rather than the morphs themselves. The latter model predicts that the fitness differences will not be constant but change with the morph frequencies. Variation in fecundity and flight muscle histolysis are shown to be more consistent with the continuous model. Data from the present experiment on JHE are inconclusive, but results from a previous selection experiment also suggest that variation in JHE is consistent only with the continuous model.     

QUANTITATIVE AND QUALITATIVE ANALYSES OF PROTEIN SYNTHESIS DURING HEAT SHOCK IN THE MARINE DIATOM NITZSCHIA ALBA (BACILLARIOPHYCEAE)1

Protein synthesis in the diatom Nitzschia alba Lewin and Lewin was drastically altered when the cells were incubated at a supraoptimal temperaeture. Quantitatively, the overall protein synthesis was greatly suppressed as indicated by teh rate of [³⁵S] methionine incorporation. The extent of suppression of protein synthesis was proportional to the severity of the heat-shock treatment which was a combination of elevated temperature and treatment duration. The in viro synthesized proteins were also qualitativelty anlayzed by two-dimensional gel electrophoresis. Dependeing on the treatment condition, a set of heat-shock proteins (HSPs) were induced. They were best detected when the cells were subjected to shocks of 35°C for 60 min or 40°C for 10 min followed by a 60 min labelling at 30°C. The results revealed 16 HSps which had moluecular weights ranging from 24–94 kD and isoelectric points ranging from 5.50–7.10. Some of the HSps were identical in molelcular weights but with differeentr isoelectric points. The induction and accumulation of HSPs in Nitzschia alba were transitory. Prologned heat-shock treatments resulted in a complete cessation of protein syntehsis and no further induction of HSPs. In all aspects, the heat shock response of diatoms was similar to that in higher plants such as soybean, maize and tobacco but differenet from most animal systems.     

蛙下丘听觉神经元对双耳强度差检测功能的建模和仿真

建立了蛙下丘听觉神经元对双耳刺激强度差检测功能的一个数学模型。按此模型所作的计算机仿真和相应实验结果比较的一致性支持了下列假设 :下丘中的EO神经元对同侧刺激不产生反应可能是由于接受了来自同侧的强烈抑制性输入 ,从而掩盖了它同时接受到的来自同侧耳的兴奋性输入。而来自同侧的抑制性输入 ,与来自对侧的兴奋性输入可能通过突触前抑制的相互作用 ,则导致了EE神经元的双耳抑制现象。     

用PFU法研究微型生物群集过程中数据的处理

根据MacArthur-Wilson的岛屿区系平衡模型S_t=S_(eq)(1-e~(GT)),可以从野外生态效应试验和室内毒性试验中,提出3个功能参数(S_(eq)、G、t_(90％))进行比较。本文提出两种计算方法:复合梯形法和最小二乘法,后者已在计算机上实现了BASIC计算程序。从数学理论上论证,最小二乘法误差较小,但如果实验布局合理,两种计算方法能得到十分一致的结果。实验模型是否符合理论模型,可以用统计学上的拟合差异度检验法来检验。     

SEX AND SPECIATION: GENETIC ARCHITECTURE AND EVOLUTIONARY POTENTIAL OF SEXUAL VERSUS NONSEXUAL TRAITS IN THE SIBLING SPECIES OF THE DROSOPHILA MELANOGASTER COMPLEX

Phenotypic divergence in the male reproductive system (genitalia and gonads) between species of the Drosophila melanogaster complex and their hybrids was quantified to decipher the role of these traits in species differentiation and speciation. Internal as well as external, sexual and nonsexual traits were analyzed with respect to genetic variation and trait asymmetry between strains within species, genetic divergence between species, and dominance and asymmetry in species and hybrids. The variation between strains within species was significant among sexual traits, and only external traits were less asymmetric than internal ones, which suggests that sexual traits are not strongly constrained within species. Three main findings show that sexual traits are most divergent between species: (1) testis length and area, and the area of the posterior lobe of the genital arch (sexual traits) showed the highest proportion of variation between species; (2) linear discriminant functions with the highest components associated to sexual traits were better predictors of species membership; and (3) testis length and area revealed a departure from a linear relationship between members of the species group. Examination of interspecific hybrids showed that sexual traits had higher asymmetry in species hybrids than in the parental species and that sexual traits showed additivity or dominance whereas nonsexual traits showed overdominance (with the exception of malpighian tubules length). These results suggest that sexual traits have undergone more genetic changes and, as a result, tend to show higher divergence and stronger hybrid breakdown between species than nonsexual traits. We propose that sexual selection in the broad sense, affecting all aspects of sexuality, may be responsible for the diversified appearance of sexual traits among closely related species and that the genetic architecture underlying sexual traits may be more prone to disruption during the early stages of speciation.     

双列杂交分析估算遗传参数两种方法的比较

本文就遗传参数的估算,对双列杂交分析中应用最广的Griffing法和Hayman法作了比较。结果表明,Hayman法估算非加性方差的结果比Griffing法小(p-1)/P~2_5~2, 而估算加性方差和遗传方差的结果则分别比Griffing法大2(p-1)/P~2_s~2和(p-1)/P~2_s~2。另外,本文还就造成两种方法分析结果不同的原因作了讨论。     

MATERNAL INHERITANCE AND ITS EFFECT ON ADAPTIVE EVOLUTION: A QUANTITATIVE GENETIC ANALYSIS OF MATERNAL EFFECTS IN A NATURAL PLANT POPULATION

A mother can influence a trait in her offspring both by the genes she transmits (Mendelian inheritance) and by maternal attributes that directly affect that trait in her offspring (maternal inheritance). Maternal inheritance can alter the direction, rate, and duration of adaptive evolution from standard Mendelian models and its impact on adaptive evolution is virtually unexplored in natural populations. In a hierarchical quantitative genetic analysis to determine the magnitude and structure of maternal inheritance in the winter annual plant, Collinsia verna, I consider three potential models of inheritance. These range from a standard Mendelian model estimating only direct (i.e., Mendelian) additive and environmental variance components to a maternal inheritance model estimating six additive and environmental variance components: direct additive and environmental variances; maternal additive and environmental variances; and the direct-maternal additive () and environmental covariances. The structure of maternal inheritance differs among the 10 traits considered at four stages in the life cycle. Early in the life cycle, seed weight and embryo weight display substantial , a negative , and a positive . Subsequently, cotyledon diameter displays and of roughly the same magnitude and negative . For fall rosettes, leaf number and length are best described by a Mendelian model. In the spring, leaf length displays maternal inheritance with significant and and a negative . All maternally inherited traits show significant negative . Predicted response to selection under maternal inheritance depends on and as well as . Negative results in predicted responses in the opposite direction to selection for seed weight and embryo weight and predicted responses near zero for all subsequent maternally inherited traits. Maternal inheritance persists through the life cycle of this annual plant for a number of size-related traits and will alter the direction and rate of evolutionary response in this population.     

SEX LINKAGE,SEX‐SPECIFIC SELECTION,AND THE ROLE OF RECOMBINATION IN THE EVOLUTION OF SEXUALLY DIMORPHIC GENE EXPRESSION

Sex‐biased genes—genes that are differentially expressed within males and females—are nonrandomly distributed across animal genomes, with sex chromosomes and autosomes often carrying markedly different concentrations of male‐ and female‐biased genes. These linkage patterns are often gene‐ and lineage‐dependent, differing between functional genetic categories and between species. Although sex‐specific selection is often hypothesized to shape the evolution of sex‐linked and autosomal gene content, population genetics theory has yet to account for many of the gene‐ and lineage‐specific idiosyncrasies emerging from the empirical literature. With the goal of improving the connection between evolutionary theory and a rapidly growing body of genome‐wide empirical studies, we extend previous population genetics theory of sex‐specific selection by developing and analyzing a biologically informed model that incorporates sex linkage, pleiotropy, recombination, and epistasis, factors that are likely to vary between genes and between species. Our results demonstrate that sex‐specific selection and sex‐specific recombination rates can generate, and are compatible with, the gene‐ and species‐specific linkage patterns reported in the genomics literature. The theory suggests that sexual selection may strongly influence the architectures of animal genomes, as well as the chromosomal distribution of fixed substitutions underlying sexually dimorphic traits.     

MORPHOLOGICAL AND GENETIC VARIATION IN THE POPULATIONS OF SARGASSUM HEMIPHYLLUM (PHAEOPHYCEAE) IN THE NORTHWESTERN PACIFIC1

Difficulty in species identification of Sargassum (Sargassaceae, Fucales) is partly attributed to the high polymorphism among its individuals and populations. This study aimed at assessing morphological and genetic variations in two varieties, var. hemiphyllum J. Agardh and var. chinense J. Agardh, of Sargassum hemiphyllum (Turner) C. Agardh, a widely distributed species in the northwestern Pacific. We investigated 26 measurable, five numerical, and 33 categorical morphological parameters associated with different branching levels of specimens from each of six localities within its distribution range using cluster analysis (CA) and principal coordinate analysis (PCoA). Leaf size of the primary and secondary branching levels and the vesicle size of the secondary branches of the specimens examined were determined to be the most important morphological parameters that were significantly different among populations. Change in leaf and vesicle length of individuals among the six populations followed a latitudinal gradient, with smaller leaves and vesicles associated with northern populations and larger ones in the southern populations. The possible influence of the gradual change in sea surface temperatures (SSTs) along this gradient in the northwestern Pacific on leaf and vesicle morphologies of this species was suggested. PCR‐RFLP analysis of the RUBISCO spacer in the chloroplast genome revealed two distinct and highly homogenous clades, a China clade and a Japan‐Korea clade, which corresponded to var. chinense and var. hemiphyllum, respectively. The formation of refugia along the “Paleo‐coast” in the East China Sea during glacial periods is suggested to have led to the vicariance of ancestral populations of S. hemiphyllum and thus to have promoted genetic differentiation. The massive freshwater outflow of the Yellow and Yangtze rivers may continue to act as a barrier, prolonging the allopatric distribution of the two varieties.     

CADMIUM AND ZINC BIOAVAILABILITY TO SELENASTRUM CAPRICORNUTUM (CHLOROPHYCEAE): ACCIDENTAL METAL UPTAKE AND TOXICITY IN THE PRESENCE OF CITRATE

Toxicities of cadmium (Cd) and zinc (Zn) to the green alga Selenastrum capricornutum Printz were determined over 72 h in defined synthetic media buffered by citrate (FRAQ_CIT ; [citrate] = 100 μM or 5 μM) or nitrilotriacetate (FRAQ_NTA ; [NTA] = 5 μM). Algal sensitivity to free Cd²+ or free Zn²+ in FRAQ_CIT was much higher than in FRAQ_NTA. In parallel experiments, short-term intracellular uptake of radiolabeled ¹⁰⁹Cd was measured as a function of time (0–30 min) in FRAQ_CIT and FRAQ_NTA; for a given free Cd²+ concentration (8, 250, or 610 nM), intracellular accumulation of Cd was consistently higher in FRAQ_CIT than in FRAQ_NTA. Under the same conditions, the alga accumulated ¹⁴C-labeled citrate almost linearly over a 2-h period. Loss of ¹⁰⁹Cd from algal cells that had been prelabeled with the radionuclide occurred slowly, and the loss rate was insensitive to the presence or absence of citrate, indicating that the overall permeability of the algal membrane to Cd was unaffected by citrate. The enhanced bioavailability of Cd in the presence of citrate could be explained by membrane transport of a charged Cd–citrate complex, presumably by accidental transport.     

EFFECTS OF LIGHT AND MELATONIN TREATMENT ON BODY TEMPERATURE AND MELATONIN SECRETION DAILY RHYTHMS IN A DIURNAL RODENT,THE FAT SAND RAT

Many mammals display predictable daily rhythmicity in both neuroendocrine function and behavior. The basic rest-activity cycles are usually consistent for a given species and vary from night-active (nocturnal), those mostly active at dawn and dusk (i.e., crepuscular), and to day-active (diurnal) species. A number of daily rhythms are oppositely phased with respect to the light/dark (LD) cycle in diurnal compared with nocturnal mammals, whereas others are equally phased with respect to the LD cycle, regardless of diurnality/nocturnality. Pineal produced melatonin (MLT) perfectly matches this phase-locked feature in that its production and secretion always occurs during the night in both diurnal and nocturnal mammals. As most rodents studied to date in the field of chronobiology are nocturnal, the aim in this study was to evaluate the effect of light manipulations and different photoperiods on a diurnal rodent, the fat sand rat, Psammomys obesus. The authors studied its daily rhythms of body temperature (T_b) and 6-sulphatoxymelatonin (6-SMT) under various photoperiodic regimes and light manipulations (acute and chronic exposures) while maintaining a constant ambient temperature of 30°C?±?1°C. The following protocols were used: (A) Control (CON) conditions 12L:12D; (A1) exposure to one light interference (LI) of CON-acclimated individuals for 30?min, 5?h after lights-off; (A2) short photoperiod (SP) acclimation (8L:16D) for 3 wks; (A3) 3 wks of SP acclimation with chronic LI of 15?min, three times a night at 4-h intervals; (A4) chronic exposure to constant dim blue light (470nm, 30 lux) for 24?h for 3 wks (LL). (B) The response to exogenous MLT administration, provided in drinking water, was measured under the following protocols: (B1) After chronic exposure to SP with LI, MLT was provided once, starting 1?h before the end of photophase; (B2) after a continuous exposure to dim blue light, MLT was provided at 15:00?h for 2?h for 2 wks; (B3) to CON animals, MLT was given intraperitoneally (i.p.) at 14:00?h. The results demonstrate that under CON acclimation, Psammomys obesus has robust T_b and 6-SMT daily rhythms in which the acrophase (peak time) of T_b is during the photophase, whereas that of 6-SMT is during scotophase. LI resulted in an elevation of T_b and a reduction of 6-SMT levels. A significant difference in the response was noted between acute and chronic exposure to LI, particularly in 6-SMT levels, which were lower than CON after LI and higher after chronic LI, implying an acclimation process. Constant exposure to blue light abolished T_b and 6-SMT rhythms in all the animals. MLT administration resumed the T_b daily rhythm in these animals, and had a recovery effect on the chronic LI-exposed animals, resulting in a T_b decrease. Altogether, the authors show in this study the different modifications of T_b rhythms and MLT levels in response to environmental light manipulations. These series of experiments may serve as a basis for establishing P. obesus as an animal model for further studies in chronobiology. (Author correspondence: hagitsc@research.haifa.ac.il)     

EVOLVABILITY OF INDIVIDUAL TRAITS IN A MULTIVARIATE CONTEXT: PARTITIONING THE ADDITIVE GENETIC VARIANCE INTO COMMON AND SPECIFIC COMPONENTS

Genetic covariation among multiple traits will bias the direction of evolution. Although a trait's phenotypic context is crucial for understanding evolutionary constraints, the evolutionary potential of one (focal) trait, rather than the whole phenotype, is often of interest. The extent to which a focal trait can evolve independently depends on how much of the genetic variance in that trait is unique. Here, we present a hypothesis‐testing framework for estimating the genetic variance in a focal trait that is independent of variance in other traits. We illustrate our analytical approach using two Drosophila bunnanda trait sets: a contact pheromone system comprised of cuticular hydrocarbons (CHCs), and wing shape, characterized by relative warps of vein position coordinates. Only 9% of the additive genetic variation in CHCs was trait specific, suggesting individual traits are unlikely to evolve independently. In contrast, most (72%) of the additive genetic variance in wing shape was trait specific, suggesting relative warp representations of wing shape could evolve independently. The identification of genetic variance in focal traits that is independent of other traits provides a way of studying the evolvability of individual traits within the broader context of the multivariate phenotype.     

VARIATION IN THE STRENGTH OF MALE MATE CHOICE ALLOWS LONG‐TERM COEXISTENCE OF SPERM‐DEPENDENT ASEXUALS AND THEIR SEXUAL HOSTS

In several asexual taxa, reproduction requires mating with related sexual species to stimulate egg development, even though genetic material is not incorporated from the sexuals (gynogenesis). In cases in which gynogens do not invest in male function, they can potentially have a twofold competitive advantage over sexuals because the asexuals avoid the cost of producing males. If unmitigated, however, the competitive success of the asexuals would ultimately lead to their own demise, following the extinction of the sexual species that stimulate egg development. We have studied a model of mate choice among sexual individuals and asexual gynogens, where males of the sexual species preferentially mate with sexual females over gynogenetic females, to determine if such mating preferences can stably maintain both gynogenetic and sexual individuals within a community. Our model shows that stable coexistence of gynogens and their sexual hosts can occur when there is variation among males in the degree of preference for mating with sexual females and when pickier males pay a higher cost of preference.