Sexual Lineages Research Articles

Niche similarity accompanying COI genetic differentiations in asexual and sexual Corbicula lineages within East and Southeast Asia

Niche similarity accompanying COI genetic differentiations in asexual and sexual Corbicula lineages within East and Southeast Asia

Contrasting Evolutionary Patterns Between Sexual and Asexual Lineages in a Genomic Region Linked to Reproductive Mode Variation in the pea aphid.

Although asexual lineages evolved from sexual lineages in many different taxa, the genetics of sex loss remains poorly understood. We addressed this issue in the pea aphid Acyrthosiphon pisum, whose natural populations encompass lineages performing cyclical parthenogenesis (CP) and producing one sexual generation per year, as well as obligate parthenogenetic (OP) lineages that can no longer produce sexual females but can still produce males. An SNP-based, whole-genome scan of CP and OP populations sequenced in pools (103 individuals from 6 populations) revealed that an X-linked region is associated with the variation in reproductive mode. This 840-kb region is highly divergent between CP and OP populations (FST = 34.9%), with >2,000 SNPs or short Indels showing a high degree of association with the phenotypic trait. In OP populations specifically, this region also shows reduced diversity and Tajima's D, consistent with the OP phenotype being a derived trait in aphids. Interestingly, the low genetic differentiation between CP and OP populations at the rest of the genome (FST = 2.5%) suggests gene flow between them. Males from OP lineages thus likely transmit their op allele to new genomic backgrounds. These genetic exchanges, combined with the selection of the OP and CP reproductive modes under different climates, probably contribute to the long-term persistence of the cp and op alleles.

Clones on the run: The genomics of a recently expanded partially clonal species.

Why species that in their core areas mainly reproduce sexually become enriched with clones in marginal populations ("geographic parthenogenesis") remains unclear. Earlier hypotheses have emphasized that selection might promote clonality because it protects locally adapted genotypes. On the other hand, it also hampers recombination and adaptation to changing conditions. The aim of the present study was to investigate the early stages of range expansion in a partially clonal species and what drives an increase in cloning during such expansion. We used genome-wide sequencing to investigate the origin and evolution of large clones formed in a macroalgal species (Fucus vesiculosus) during a recent expansion into the postglacial Baltic Sea. We found low but persistent clonality in core populations, while at range margins, large dominant clonal lineages had evolved repeatedly from different sexual populations. A range expansion model showed that even when asexual recruitment is less favourable than sexual recruitment in core populations, repeated bottlenecks at the expansion front can establish a genetically eroded clonal wave that spreads ahead of a sexual wave into the new area. Genetic variation decreases by drift following repeated bottlenecks at the expansion front. This results in the emerging clones having low expected heterozygosity, which corroborated our empirical observations. We conclude that Baker's Law (clones being favoured by uniparental reproductive assurance in new areas) can play an important role during range expansion in partially clonal species, resulting in a complex spatiotemporal mosaic of clonal and sexual lineages that might persist during thousands of generations.

Evidence for stronger discrimination between conspecific and heterospecific mating partners in sexual vs. asexual female freshwater snails.

Once-useful traits that no longer contribute to fitness tend to decay over time. Here, we address whether the expression of mating-related traits that increase the fitness of sexually reproducing individuals but are likely less useful or even costly to asexual counterparts seems to exhibit decay in the latter. Potamopyrgus antipodarum is a New Zealand freshwater snail characterized by repeated transitions from sexual to asexual reproduction. The frequent coexistence of sexual and asexual lineages makes P. antipodarum an excellent model for the study of mating-related trait loss. Under the presumption (inherent in the Biological Species Concept) that failure to discriminate between conspecific and heterospecific mating partners represents a poor mate choice, we used a mating choice assay including sexual and asexual P. antipodarum females and conspecific (presumed better choice) vs. heterospecific (presumed worse choice) males to evaluate the loss of behavioral traits related to sexual reproduction. We found that sexual females engaged in mating behaviors with conspecific mating partners more frequently and for a greater duration than with heterospecific mating partners. By contrast, asexual females mated at similar frequency and duration as sexual females, but did not mate more often or for longer duration with conspecific vs. heterospecific males. While further confirmation will require inclusion of a more diverse array of sexual and asexual lineages, these results are consistent with a scenario where selection acting to maintain effective mate discrimination in asexual P. antipodarum is weak or ineffective relative to sexual females and, thus, where asexual reproduction is associated with the evolutionary decay of mating-related traits in this system.

Substantial genetic mixing among sexual and androgenetic lineages within the clam genus Corbicula

“Occasional” sexuality occurs when a species combines clonal reproduction and genetic mixing. This strategy is predicted to combine the advantages of both asexuality and sexuality, but its actual consequences on the genetic diversity and species longevity are poorly understood. Androgenesis, a reproductive mode in which the offspring inherits its entire nuclear genome from the father, is often reported as a strictly clonal reproductive mode. Androgenesis is the predominant reproductive mode within the hermaphroditic, invasive lineages of the mollusk genus Corbicula. Their ability to reproduce clonally through androgenesis has been determinant in their invasive success, having colonized during the 20th century American and European freshwater systems, where they became notorious invaders with a widespread, global distribution. However, in androgenetic Corbicula clams, occasional genetic mixing between distinct lineages has also been observed when the sperm of one lineage fertilizes the oocyte of another one. Because of these occasional introgressions, the genetic relationships between Corbicula species remained unclear, and the biogeographic origins of the invasive androgenetic lineages have been challenging to identify. To address these issues, we analyzed the patterns of allele sharing for several nuclear and mitochondrial molecular markers among Corbicula individuals collected across both the native and invasive range. Our results show the occurrence of an allelic pool encompassing all Corbicula freshwater species worldwide, including sexual and androgenetic ones, which highlights the substantial genetic mixing within this genus. However, the differences in allele sharing patterns between invasive lineages, and the low diversity within each lineage, suggest recent, distinct biogeographic origins of invasive Corbicula androgenetic lineages. Finally, the polyploidy, high heterozygosity, and hybrid phenotypes and genotypes found in our study probably originated from hybridization events following egg parasitism between distinct Corbicula lineages. This extensive cross-lineage mixing found in Corbicula may generate nuclear diversity in an otherwise asexually reproducing species.

Seed size, endosperm and germination variation in sexual and apomictic Boechera.

Asexual reproduction results in offspring that are genetically identical to the mother. Among apomictic plants (reproducing asexually through seeds) many require paternal genetic contribution for proper endosperm development (pseudogamous endosperm). We examined phenotypic diversity in seed traits using a diverse panel of sexual and apomictic accessions from the genus Boechera. While genetic uniformity resulting from asexual reproduction is expected to reduce phenotypic diversity in seeds produced by apomictic individuals, pseudogamous endosperm, variable endosperm ploidy, and the deviations from 2:1 maternal:paternal genome ratio in endosperm can all contribute to increased phenotypic diversity among apomictic offspring. We characterized seed size variation in 64 diploid sexual and apomictic (diploid and triploid) Boechera lineages. In order to find out whether individual seed size was related to endosperm ploidy we performed individual seed measurements (projected area and mass) using the phenoSeeder robot system and flow cytometric seed screen. In order to test whether individual seed size had an effect on resulting fitness we performed a controlled growth experiment and recorded seedling life history traits (germination success, germination timing, and root growth rate). Seeds with triploid embryos were 33% larger than those with diploid embryos, but no average size difference was found between sexual and apomictic groups. We identified a maternal effect whereby chloroplast lineage 2 had 30% larger seeds than lineage 3, despite having broad and mostly overlapping geographic ranges. Apomictic seeds were not more uniform in size than sexual seeds, despite genetic uniformity of the maternal gametophyte in the former. Among specific embryo/endosperm ploidy combinations, seeds with tetraploid (automomous) endosperm were on average smaller, and the proportion of such seeds was highest in apomicts. Larger seeds germinated more quickly than small seeds, and lead to higher rates of root growth in young seedlings. Seed mass is under balancing selection in Boechera, and it is an important predictor of several traits, including germination probability and timing, root growth rates, and developmental abnormalities in apomictic accessions.

The protective effect of L-carnitine on testosterone synthesis pathway, and spermatogenesis in monosodium glutamate-induced rats

BackgroundMonosodium glutamate (MSG) is a food ingredient that is increasingly used commercially. MSG leads to oxidative stress, consequently suppressing steroid hormone production that causes defects in male reproductive system. This study aimed to evaluate the effect of L-carnitine as an antioxidant on testicular damage in MSG-induced male rats.MethodsSixty adult male Spargue-Dawley rats were randomly divided into six groups of ten as follows: control (water), sham (normal saline), L-carnitine (200 mg/kg b.w), MSG (3 g/kg b.w), MSG + L-carnitine 100 (3 g/kg b.w of MSG and 100 mg/kg b.w of L-carnitine), and MSG + L-carnitine 200 (3 g/kg b.w of MSG and 200 mg/kg b.w of L-carnitine). The treatment was administered by oral gavage for six months. Serum levels of Malondialdehyde (MDA), Total Anti-oxidant Capacity (TAC), LH, FSH, testosterone, and mRNA expressions of Star, Cyp11a1, and Hsd17b3 genes, and histological and stereological changes were assessed.ResultsL-carnitine led to a significant decrease in the level of MDA and a significant rise in the serum levels of TAC, LH, FSH, and mRNA expression of Star and Cyp11a1 compared to the MSG group (p < 0.05). Furthermore, stereological results indicated a significant increment in the number of sexual lineage cells, the total volume of the testis, length, diameter, and volume of seminiferous tubules, the height of the germinal epithelium, sperm count, and sperm motility (p < 0.05) in MSG + L-carnitine 200 compare to MSG group.ConclusionThe study’s findings demonstrated that L-carnitine due to its anti-oxidant properties, ameliorated the reproductive abnormalities in the male rats exposed to MSG.

ZW sex-chromosome evolution and contagious parthenogenesis in Artemia brine shrimp.

Eurasian brine shrimp (genus Artemia) have closely related sexual and asexual lineages of parthenogenetic females, which produce rare males at low frequencies. Although they are known to have ZW chromosomes, these are not well characterized, and it is unclear whether they are shared across the clade. Furthermore, the underlying genetic architecture of the transmission of asexuality, which can occur when rare males mate with closely related sexual females, is not well understood. We produced a chromosome-level assembly for the sexual Eurasian species Artemia sinica and characterized in detail the pair of sex chromosomes of this species. We combined this new assembly with short-read genomic data for the sexual species Artemia sp. Kazakhstan and several asexual lineages of Artemia parthenogenetica, allowing us to perform an in-depth characterization of sex-chromosome evolution across the genus. We identified a small differentiated region of the ZW pair that is shared by all sexual and asexual lineages, supporting the shared ancestry of the sex chromosomes. We also inferred that recombination suppression has spread to larger sections of the chromosome independently in the American and Eurasian lineages. Finally, we took advantage of a rare male, which we backcrossed to sexual females, to explore the genetic basis of asexuality. Our results suggest that parthenogenesis is likely partly controlled by a locus on the Z chromosome, highlighting the interplay between sex determination and asexuality.

Patterns of gene expression in ovaries of sexual vs. asexual lineages of a freshwater snail

Why sexual reproduction is so common when asexual reproduction should be much more efficient and less costly remains an open question in evolutionary biology. Comparisons between otherwise similar sexual and asexual taxa allow us to characterize the genetic architecture underlying asexuality, which can, in turn, illuminate how this reproductive mode transition occurred and the mechanisms by which it is maintained or disrupted. Here, we used transcriptome sequencing to compare patterns of ovarian gene expression between actively reproducing obligately sexual and obligately asexual females from multiple lineages of Potamopyrgus antipodarum, a freshwater New Zealand snail characterized by frequent separate transitions to asexuality and coexistence of otherwise similar sexual and asexual lineages. We also used these sequence data to evaluate whether population history accounts for variation in patterns of gene expression. We found that source population was a major source of gene expression variation, and likely more influential than reproductive mode. This outcome for these common garden-raised snails is strikingly similar to earlier results from field-collected snails. While we did not identify a likely set of candidate genes from expression profiles that could plausibly explain how transitions to asexuality occurred, we identified around 1,000 genes with evidence of differential expression between sexual and asexual reproductive modes, and 21 genes that appear to exhibit consistent expression differences between sexuals and asexuals across genetic backgrounds. This second smaller set of genes provides a good starting point for further exploration regarding a potential role in the transition to asexual reproduction. These results mark the first effort to characterize the causes of asexuality in P. antipodarum, demonstrate the apparently high heritability of gene expression patterns in this species, and hint that for P. antipodarum, transitions to asexuality might not necessarily be strongly associated with broad changes in gene expression.

Substantial genetic mixing among sexual and androgenetic lineages within the clam genus Corbicula

Substantial genetic mixing among sexual and androgenetic lineages within the clam genus Corbicula

Enhanced heterozygosity from male meiotic chromosome chains is superseded by hybrid female asexuality in termites

Although males are a ubiquitous feature of animals, they have been lost repeatedly in diverse lineages. The tendency for obligate asexuality to evolve is thought to be reduced in animals whose males play a critical role beyond the contribution of gametes, for example, via care of offspring or provision of nuptial gifts. To our knowledge, the evolution of obligate asexuality in such species is unknown. In some species that undergo frequent inbreeding, males are hypothesized to play a key role in maintaining genetic heterozygosity through the possession of neo-sex chromosomes, although empirical evidence for this is lacking. Because inbreeding is a key feature of the life cycle of termites, we investigated the potential role of males in promoting heterozygosity within populations through karyotyping and genome-wide single-nucleotide polymorphism analyses of the drywood termite Glyptotermes nakajimai We showed that males possess up to 15 out of 17 of their chromosomes as sex-linked (sex and neo-sex) chromosomes and that they maintain significantly higher levels of heterozygosity than do females. Furthermore, we showed that two obligately asexual lineages of this species-representing the only known all-female termite populations-arose independently via intraspecific hybridization between sexual lineages with differing diploid chromosome numbers. Importantly, these asexual females have markedly higher heterozygosity than their conspecific males and appear to have replaced the sexual lineages in some populations. Our results indicate that asexuality has enabled females to supplant a key role of males.

Condition-dependent sexual reproduction is driven by benefits, not costs of sex

AbstractFacultative sexual organisms must allocate resources to both asexual and sexual reproduction. Optimal patterns of investment in sex depend on the relative costs and benefits of each reproductive mode, and may consequently be context- and condition-dependent. Two proposed explanations for the observed variation in investment in sex among facultative sexual lineages invoke alternative condition-dependent scenarios. Under the fitness-associated sex hypothesis, sex is predicted when individuals are in poor condition or experience stressful environments. Under the resource-demanding sex hypothesis, sex is only affordable to individuals in good condition experiencing favourable environments. Direct tests of these contrasting hypotheses are rare; moreover, investment in different components of sexual reproduction responds differently to cues promoting sex, and may be subject to different energetic constraints. Using genotypes of facultative sexual Daphnia carinata that differ in their level of investment in sex, we manipulated resource availability while accounting for day length (a seasonal cue for sex) to evaluate these hypotheses. The sexual response to day length depended on resource availability: increased day lengths and reduced food availability increased the production of sexual eggs, and relative investment in males, in a manner consistent with the fitness-associated sex hypothesis. The pattern of condition-dependence was specific to each component of reproductive investment – while male production covaried with asexual fecundity across genotypes, increased sexual egg production was associated with reduced asexual reproduction. Our results suggest that investment in sex is determined largely by its context-dependent advantages, and that this investment is not moderated by immediate costs to asexual reproduction.

Transitions to asexuality and evolution of gene expression in Artemia brine shrimp.

While sexual reproduction is widespread among many taxa, asexual lineages have repeatedly evolved from sexual ancestors. Despite extensive research on the evolution of sex, it is still unclear whether this switch represents a major transition requiring major molecular reorganization, and how convergent the changes involved are. In this study, we investigated the phylogenetic relationship and patterns of gene expression of sexual and asexual lineages of Eurasian Artemia brine shrimp, to assess how gene expression patterns are affected by the transition to asexuality. We find only a few genes that are consistently associated with the evolution of asexuality, suggesting that this shift may not require an extensive overhauling of the meiotic machinery. While genes with sex-biased expression have high rates of expression divergence within Eurasian Artemia, neither female- nor male-biased genes appear to show unusual evolutionary patterns after sexuality is lost, contrary to theoretical expectations.

Calculating Orthologous Protein-Coding Sequence Set Probability Using the Poisson Process.

We extend the popular Jukes-Cantor evolution model and calculate the probability of an orthologous nucleotide sequence set [a reference sequence (B1) stays with the other sequences (B-1)], where the sequence evolution [from a last common ancestral sequence (ɑ)] follows the (prospective) Poisson process with the overall event rate λ prorated among mutation types (nucleotide/codon substitution, insertion, and deletion) and sites along each sequence. The corresponding retrospective process (reversing the prospective process) facilitates developing algorithms to calculate the marginal probability [Pr(B1)] (Monte Carlo integration) and sample ɑ (given B1). We calculate probability Pr(B-1|ɑ) based on the identified events (during "ɑ→B-1") from pairwise sequence alignment to implement Pr(B-1|B1) calculation (Monte Carlo integration). Event queue sampling and probability magnifiers are used to improve the computational efficiency when the number of events is large. We finally test our procedure on both simulated and recently studied hexapod transcriptome data (Brandt et al.), where each asexual lineage pairs with its closest related sexual lineage. Rate estimates (for Phasmatodea and Zygentoma) and model comparison indicate that the asexual lineages likely mutate several times faster than their sexual relatives.

Male and female reproductive fitness costs of an immune response in natural populations.

Parasites can mediate host fitness both directly, via effects on survival and reproduction, or indirectly by inducing host immune defense with costly side-effects. The evolution of immune defense is determined by a complex interplay of costs and benefits of parasite infection and immune response, all of which may differ for male and female hosts in sexual lineages. Here, we examine fitness costs associated with an inducible immune defense in a fish-cestode host-parasite system. Cestode infection induces peritoneal fibrosis in threespine stickleback (Gasterosteus aculeatus), constraining cestode growth and sometimes encasing and killing the parasite. Surveying two wild populations of stickleback, we confirm that the presence of fibrosis scar tissue is associated with reduced parasite burden in both male and female fish. However, fibrotic fish had lower foraging success and reproductive fitness (reduced female egg production and male nesting success), indicating strong costs of the lingering immunopathology. Consistent with substantial sexually concordant fitness effects of immune response, we find alignment of multivariate selection across the sexes despite sexual antagonism over morphological shape. Although both sexes experienced costs of fibrosis, the net impacts are unequal because in the two study populations females had higher cestode exposure. To evaluate whether this difference in risk should drive sex-specific immune strategies, we analyze a quantitative genetic model of host immune response to a trophically transmitted parasite. The model and empirical data illustrate how shared costs and benefits of immune response lead to shared evolutionary interests of male and female hosts, despite unequal infection risks across the sexes.

Asexuality Associated with Marked Genomic Expansion of Tandemly Repeated rRNA and Histone Genes.

How does asexual reproduction influence genome evolution? Although is it clear that genomic structural variation is common and important in natural populations, we know very little about how one of the most fundamental of eukaryotic traits—mode of genomic inheritance—influences genome structure. We address this question with the New Zealand freshwater snail Potamopyrgus antipodarum, which features multiple separately derived obligately asexual lineages that coexist and compete with otherwise similar sexual lineages. We used whole-genome sequencing reads from a diverse set of sexual and asexual individuals to analyze genomic abundance of a critically important gene family, rDNA (the genes encoding rRNAs), that is notable for dynamic and variable copy number. Our genomic survey of rDNA in P. antipodarum revealed two striking results. First, the core histone and 5S rRNA genes occur between tandem copies of the 18S–5.8S–28S gene cluster, a unique architecture for these crucial gene families. Second, asexual P. antipodarum harbor dramatically more rDNA–histone copies than sexuals, which we validated through molecular and cytogenetic analysis. The repeated expansion of this genomic region in asexual P. antipodarum lineages following distinct transitions to asexuality represents a dramatic genome structural change associated with asexual reproduction—with potential functional consequences related to the loss of sexual reproduction.

Rapid local adaptation in both sexual and asexual invasive populations of monkeyflowers (Mimulus spp.).

Traditionally, local adaptation has been seen as the outcome of a long evolutionary history, particularly with regard to sexual lineages. By contrast, phenotypic plasticity has been thought to be most important during the initial stages of population establishment and in asexual species. We evaluated the roles of adaptive evolution and phenotypic plasticity in the invasive success of two closely related species of invasive monkeyflowers (Mimulus) in the UK that have contrasting reproductive strategies: M. guttatus combines sexual (seeds) and asexual (clonal growth) reproduction while M. × robertsii is entirely asexual. We compared the clonality (number of stolons), floral and vegetative phenotype, and phenotypic plasticity of native (M. guttatus) and invasive (M. guttatus and M. × robertsii) populations grown in controlled environment chambers under the environmental conditions at each latitudinal extreme of the UK. The goal was to discern the roles of temperature and photoperiod on the expression of phenotypic traits. Next, we tested the existence of local adaptation in the two species within the invasive range with a reciprocal transplant experiment at two field sites in the latitudinal extremes of the UK, and analysed which phenotypic traits underlie potential local fitness advantages in each species. Populations of M. guttatus in the UK showed local adaptation through sexual function (fruit production), while M. × robertsii showed local adaptation via asexual function (stolon production). Phenotypic selection analyses revealed that different traits are associated with fitness in each species. Invasive and native populations of M. guttatus had similar phenotypic plasticity and clonality. M. × robertsii presents greater plasticity and clonality than native M. guttatus, but most populations have restricted clonality under the warm conditions of the south of the UK. This study provides experimental evidence of local adaptation in a strictly asexual invasive species with high clonality and phenotypic plasticity. This indicates that even asexual taxa can rapidly (<200 years) adapt to novel environmental conditions in which alternative strategies may not ensure the persistence of populations.

Breeding system of diploid sexuals within the Ranunculus auricomus complex and its role in a geographical parthenogenesis scenario.

The larger distribution area of asexuals compared with their sexual relatives in geographical parthenogenesis (GP) scenarios has been widely attributed to the advantages of uniparental reproduction and polyploidy. However, potential disadvantages of sexuals due to their breeding system have received little attention so far. Here, we study the breeding system of five narrowly distributed sexual lineages of Ranunculus notabilis s.l. (R. auricomus complex) and its effects on outcrossing, inbreeding, female fitness, and heterozygosity. We performed selfing and intra‐ and interlineage crossings by bagging 481 flowers (59 garden individuals) followed by germination experiments. We compared seed set and germination rates, and related them to genetic distance and genome‐wide heterozygosity (thousands of RADseq loci). Selfings (2.5%) unveiled a significantly lower seed set compared with intra‐ (69.0%) and interlineage crossings (69.5%). Seed set of intra‐ (65%) compared to interpopulation crossings (78%) was significantly lower. In contrast, all treatments showed comparable germination rates (32%–43%). Generalized linear regressions between seed set and genetic distance revealed positive relationships in general and between lineages, and a negative one within lineages. Seed set was the main decisive factor for female fitness. Germination rates were not related to genetic distance at any level, but were positively associated with heterozygosity in interlineage crossings. Experiments confirmed full crossability and predominant outcrossing among sexual R. notabilis s.l. lineages. However, up to 5% (outliers 15%–31%) of seeds were formed by selfing, probably due to semi‐self‐compatibility in a multi‐locus gametophytic SI system. Less seed set in intrapopulation crossings, and higher seed set and germination rates from crossings of genetically more distant and heterozygous lineages (interlineage) indicate negative inbreeding and positive outbreeding effects. In GP scenarios, sexual species with small and/or isolated populations can suffer from decreased female fitness due to their breeding system. This factor, among others, probably limits range expansion of sexuals.

Identifying the fitness consequences of sex in complex natural environments.

In the natural world, sex prevails, despite its costs. Although much effort has been dedicated to identifying the intrinsic costs of sex (e.g., the cost of males), few studies have identified the ecological fitness consequences of sex. Furthermore, correlated biological traits that differ between sexuals and asexuals may alter these costs, or even render the typical costs of sex irrelevant. We conducted a large‐scale, multisite, reciprocal transplant using multiple sexual and asexual genotypes of a native North American wildflower to show that sexual genotypes have reduced lifetime fitness, despite lower herbivory. We separated the effects of sex from those of hybridity, finding that overwinter survival is elevated in asexuals regardless of hybridity, but herbivores target hybrid asexuals more than nonhybrid asexual or sexual genotypes. Survival is lowest in homozygous sexual lineages, implicating inbreeding depression as a cost of sex. Our results show that the consequences of sex are shaped not just by sex itself, but by complex natural environments, correlated traits, and the identity and availability of mates.

Toxic males: Density‐dependent male mating harassment can explain geographic parthenogenesis

AbstractUnderstanding evolutionary forces that maintain sexual reproduction has been a central question in evolutionary biology, because the well‐known “two‐fold” cost of sex demands an explanation, especially since asexual reproduction (vegetative reproduction and parthenogenesis) often arises by mutation and hybridization. Geographic parthenogenesis, where asexual reproduction by close relatives of sexual lineages tends to be found in “marginal” habitats, has been considered a key pattern for studying the maintenance of sex. Recent studies of sexually reproducing populations have revealed that mating attempts by males can be harmful to females due to sexual conflict. Where that occurs and where there is a mixture of both sexual and asexual females, if male mating harassment is more harmful to asexual females than to sexual females, this may prevent the invasion of rare asexuals. Here, we use a simple ecological model with density‐dependent male mating harassment to show how asexual reproduction can become prevalent in habitats when environmental conditions are poor (i.e., low carrying capacities). The maintenance of sex is possible in a bistable system where abundant males are so harmful that they prevent the invasion of rare asexuals. However, male density is not large enough to suppress asexuals and sex is lost when carrying capacities are low or when sex ratio is highly male‐biased. While previous theoretical studies employed a frequency‐dependent model to study male mating harassment, we demonstrate that density‐dependence may be important to explain the dynamics in low population densities where there is geographical parthenogenesis.