Population genetic structure of parthenogenetic flatworm populations with occasional sex

1. Sexual populations are expected to perform better in fluctuating environments than asexuals because recombination provides the potential to adapt to changing environments due to increased genetic variation. Nevertheless, some asexual species show comparably high levels of genotypic diversity. Such diversity might be achieved through gene flow between coexisting sexual and asexual populations or through sexual events within asexual populations. 2. Evidence for occasional sex in the flatworm Schmidtea polychroa was previously found at one specific site that is inhabited by parthenogenetic forms. There, varying rates of sex between subpopulations, reaching up to 12%, were observed. Past recurrent sexual processes left a significant genetic signature in the population genetic structure of this population. In the present study, we examined the population genetic structure of six independent metapopulations (lakes) of the freshwater planarian flatworm S. polychroa, to confirm the presence of occasional sex and that its population genetic consequences can be generalised. 3. Using microsatellites, we found varying rates of occasional sex among subpopulations. Metapopulations showed medium to high levels of genotypic diversity that correlated with the rate of sex. 4. We conclude that occasional sex has considerable consequences for population genetic structure of parthenogenetic species and promotes diversity that might allow response to the particular type of selection that is usually predicted to favour sexual reproduction. This reproductive strategy provides genetic characteristics required for selection to act on, and might, therefore, explain the success of this parthenogenetic species.     

Sperm length and quality in sperm-dependent parthenogens

Is the male function of parthenogenetic hermaphrodites a historic remnant or an adaptive investment? In the first case, one would expect the male function to be ‘degraded’ and, in the second case, fully functional, despite being reduced. To answer this question, we investigated sperm morphology and quality in the hermaphroditic planarian Schmidtea polychroa. Parthenogenetic forms of this species are sperm‐dependent and require allosperm to trigger embryogenesis. Although there is generally no genetic paternal contribution, occasional fertilization of parthenogenetic eggs has been reported, leading to genetic changes in an otherwise apomictic lineage (‘occasional sex’). Comparing two locations, one with pure clonality and one with occasional sexuality, revealed substantial differences in both sperm length and fertility: Occasional sex appears to be coupled with longer sperm and higher sperm quality. We discuss the implications of these results for the maintenance of the male function in sperm‐dependent parthenogens. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 93 , 81–87.     

Geographic variation in sex‐chromosome differentiation in the common frog (Rana temporaria)

In sharp contrast with birds and mammals, sex‐determination systems in ectothermic vertebrates are often highly dynamic and sometimes multifactorial. Both environmental and genetic effects have been documented in common frogs (Rana temporaria). One genetic linkage group, mapping to the largest pair of chromosomes and harbouring the candidate sex‐determining gene Dmrt1, associates with sex in several populations throughout Europe, but association varies both within and among populations. Here, we show that sex association at this linkage group differs among populations along a 1500‐km transect across Sweden. Genetic differentiation between sexes is strongest (F_ST = 0.152) in a northern‐boreal population, where male‐specific alleles and heterozygote excesses (F_IS = ?0.418 in males, +0.025 in females) testify to a male‐heterogametic system and lack of X‐Y recombination. In the southernmost population (nemoral climate), in contrast, sexes share the same alleles at the same frequencies (F_ST = 0.007 between sexes), suggesting unrestricted recombination. Other populations show intermediate levels of sex differentiation, with males falling in two categories: some cluster with females, while others display male‐specific Y haplotypes. This polymorphism may result from differences between populations in the patterns of X‐Y recombination, co‐option of an alternative sex‐chromosome pair, or a mixed sex‐determination system where maleness is controlled either by genes or by environment depending on populations or families. We propose approaches to test among these alternative models, to disentangle the effects of climate and phylogeography on the latitudinal trend, and to sort out how this polymorphism relates to the ‘sexual races’ described in common frogs in the 1930s.     

Occasional recombination of a selfish X‐chromosome may permit its persistence at high frequencies in the wild

The sex‐ratio X‐chromosome (SR) is a selfish chromosome that promotes its own transmission to the next generation by destroying Y‐bearing sperm in the testes of carrier males. In some natural populations of the fly Drosophila neotestacea, up to 30% of the X‐chromosomes are SR chromosomes. To investigate the molecular evolutionary history and consequences of SR, we sequenced SR and standard (ST) males at 11 X‐linked loci that span the ST X‐chromosome and at seven arbitrarily chosen autosomal loci from a sample of D. neotestacea males from throughout the species range. We found that the evolutionary relationship between ST and SR varies among individual markers, but genetic differentiation between SR and ST is chromosome‐wide and likely due to large chromosomal inversions that suppress recombination. However, SR does not consist of a single multilocus haplotype: we find evidence for gene flow between ST and SR at every locus assayed. Furthermore, we do not find long‐distance linkage disequilibrium within SR chromosomes, suggesting that recombination occurs in females homozygous for SR. Finally, polymorphism on SR is reduced compared to that on ST, and loci displaying signatures of selection on ST do not show similar patterns on SR. Thus, even if selection is less effective on SR, our results suggest that gene flow with ST and recombination between SR chromosomes may prevent the accumulation of deleterious mutations and allow its long‐term persistence at relatively high frequencies.     

Low rates of X‐Y recombination,not turnovers,account for homomorphic sex chromosomes in several diploid species of Palearctic green toads (Bufo viridis subgroup)

Contrasting with birds and mammals, most ectothermic vertebrates present homomorphic sex chromosomes, which might be due either to a high turnover rate or to occasional X‐Y recombination. We tested these two hypotheses in a group of Palearctic green toads that diverged some 3.3 million years ago. Using sibship analyses of sex‐linked markers, we show that all four species investigated share the same pair of sex chromosomes and a pattern of male heterogamety with drastically reduced X‐Y recombination in males. Phylogenetic analyses of sex‐linked sequences show that X and Y alleles cluster by species, not by gametolog. We conclude that X‐Y homomorphy and fine‐scale sequence similarity in these species do not stem from recent sex‐chromosome turnovers, but from occasional X‐Y recombination.     

A cryptic heterogametic transition revealed by sex-linked DNA markers in Palearctic green toads

In sharp contrast to birds and mammals, most cold‐blooded vertebrates have homomorphic (morphologically undifferentiated) sex chromosomes. This might result either from recurrent X‐Y recombination (occurring e.g. during occasional events of sex reversal) or from frequent turnovers (during which sex‐determining genes are overthrown by new autosomal mutations). Evidence for turnovers is indeed mounting in fish, but very few have so far been documented in amphibians, possibly because of practical difficulties in identifying sex chromosomes. Female heterogamety (ZW) has long been established in Bufo bufo, based on sex reversal and crossing experiments. Here, we investigate a sex‐linked marker identified from a laboratory cross between Palearctic green toads (Bufo viridis subgroup). The F₁ offspring produced by a female Bufo balearicus and a male Bufo siculus were phenotypically sexed, displaying an even sex ratio. A sex‐specific marker detected in highly reproducible AFLP genotypes was cloned. Sequencing revealed a noncoding, microsatellite‐containing fragment. Reamplification and genotyping of families of this and a reciprocal cross showed B. siculus to be male heterogametic (XY) and suggested the same system for B. balearicus. Our results thus reveal a cryptic heterogametic transition within bufonid frogs and help explain patterns of hybrid fitness within the B. viridis subgroup. Turnovers of genetic sex‐determination systems may be more frequent in amphibians than previously thought and thus contribute to the prevalence of homomorphic sex chromosomes in this group.     

Clonal reproduction shapes evolution in the lizard malaria parasite Plasmodium floridense

The preponderant clonal evolution hypothesis (PCE) predicts that frequent clonal reproduction (sex between two clones) in many pathogens capable of sexual recombination results in strong linkage disequilibrium and the presence of discrete genetic subdivisions characterized by occasional gene flow. We expand on the PCE and predict that higher rates of clonal reproduction will result in: (1) morphologically cryptic species that exhibit (2) low within‐species variation and (3) recent between‐species divergence. We tested these predictions in the Caribbean lizard malaria parasite Plasmodium floridense using 63 single‐infection samples in lizards collected from across the parasite's range, and sequenced them at two mitochondrial, one apicoplast, and five nuclear genes. We identified 11 provisionally cryptic species within P. floridense, each of which exhibits low intraspecific variation and recent divergence times between species (some diverged approximately 110,000 years ago). Our results are consistent with the hypothesis that clonal reproduction can profoundly affect diversification of species capable of sexual recombination, and suggest that clonal reproduction may have led to a large number of unrecognized pathogen species. The factors that may influence the rates of clonal reproduction among pathogens are unclear, and we discuss how prevalence and virulence may relate to clonal reproduction.     

Population structure and geographical subdivision of the Leishmania major vector Phlebotomus papatasi as revealed by microsatellite variation

Abstract Multi‐locus microsatellite typing (MLMT) has been employed to infer the population structure of Phlebotomus papatasi (Scopoli) (Diptera: Psychodidae) sandflies and assign individuals to populations. Phlebotomus papatasi sandflies were collected from 35 sites in 15 countries. A total of 188 P. papatasi individuals were typed using five microsatellite loci, resulting in 113 different genotypes. Unique microsatellite signatures were observed for some of the populations analysed. Comparable results were obtained when the data were analysed with Bayesian model and distance‐based methods. Bayesian statistic‐based analyses split the dataset into two distinct genetic clusters, A and B, with further substructuring within each. Population A consisted of five subpopulations representing large numbers of alleles that were correlated with the geographical origins of the sandflies. Cluster B comprised individuals collected in the Middle East and the northern Mediterranean area. The subpopulations B1 and B2 did not, however, show any further correlation to geographical origin. The genetic differentiation between subpopulations was supported by F statistics showing statistically significant (Bonferroni‐corrected P < 0.005) values of 0.221 between B2 and B1 and 0.816 between A5 and A4. Identification of the genetic structure of P. papatasi populations is important for understanding the patterns of dispersal of this species and to developing strategies for sandfly control.     

Cryptic Sex in Symbiodinium (Alveolata,Dinoflagellata) is Supported by an Inventory of Meiotic Genes

Symbiodinium encompasses a diverse clade of dinoflagellates that are ecologically important as symbionts of corals and other marine organisms. Despite decades of study, cytological evidence of sex (karyogamy and meiosis) has not been demonstrated in Symbiodinium, although molecular population genetic patterns support the occurrence of sexual recombination. Here, we provide additional support for sex in Symbiodinium by uncovering six meiosis‐specific and 25 meiosis‐related genes in three published genomes. Cryptic sex may be occurring in Symbiodinium's seldom‐seen free‐living state while being inactive in the symbiotic state.     

Cross‐sex genetic covariances limit the evolvability of wing‐shape within and among species of Drosophila

The independent evolution of males and females is potentially constrained by both sexes inheriting the same alleles from their parents. This genetic constraint can limit the evolvability of complex traits; however, there are few studies of multivariate evolution that incorporate cross‐sex genetic covariances in their predictions. Drosophila wing‐shape has emerged as a model high‐dimensional phenotype; wing‐shape is highly evolvable in contemporary populations, and yet perplexingly stable across phylogenetic timescales. Here, we show that cross‐sex covariances in Drosophila melanogaster, given by the B ‐matrix, may considerably bias wing‐shape evolution. Using random skewers, we show that B would constrain the response to antagonistic selection by 90%, on average, but would double the response to concordant selection. Both cross‐sex within‐trait and cross‐sex cross‐trait covariances determined the predicted response to antagonistic selection, but only cross‐sex within‐trait covariances facilitated the predicted response to concordant selection. Similar patterns were observed in the direction of extant sexual dimorphism in D. melanogaster, and in directions of most and least dimorphic variation across the Drosophila phylogeny. Our results highlight the importance of considering between‐sex genetic covariances when making predictions about evolution on both macro‐ and microevolutionary timescales, and may provide one more explanatory piece in the puzzle of stasis.     

Cryptic recombination in the ever-young sex chromosomes of Hylid frogs

Sex chromosomes are expected to evolve suppressed recombination, which leads to degeneration of the Y and heteromorphism between the X and Y. Some sex chromosomes remain homomorphic, however, and the factors that prevent degeneration of the Y in these cases are not well understood. The homomorphic sex chromosomes of the European tree frogs (Hyla spp.) present an interesting paradox. Recombination in males has never been observed in crossing experiments, but molecular data are suggestive of occasional recombination between the X and Y. The hypothesis that these sex chromosomes recombine has not been tested statistically, however, nor has the X‐Y recombination rate been estimated. Here, we use approximate Bayesian computation coupled with coalescent simulations of sex chromosomes to quantify X‐Y recombination rate from existent data. We find that microsatellite data from H. arborea, H. intermedia and H. molleri support a recombination rate between X and Y that is significantly different from zero. We estimate that rate to be approximately 10⁵ times smaller than that between X chromosomes. Our findings support the notion that very low recombination rate may be sufficient to maintain homomorphism in sex chromosomes.     

Local genes for local bacteria: Evidence of allopatry in the genomes of transatlantic Campylobacter populations

The genetic structure of bacterial populations can be related to geographical locations of isolation. In some species, there is a strong correlation between geographical distance and genetic distance, which can be caused by different evolutionary mechanisms. Patterns of ancient admixture in Helicobacter pylori can be reconstructed in concordance with past human migration, whereas in Mycobacterium tuberculosis it is the lack of recombination that causes allopatric clusters. In Campylobacter, analyses of genomic data and molecular typing have been successful in determining the reservoir host species, but not geographical origin. We investigated biogeographical variation in highly recombining genes to determine the extent of clustering between genomes from geographically distinct Campylobacter populations. Whole‐genome sequences from 294 Campylobacter isolates from North America and the UK were analysed. Isolates from within the same country shared more recently recombined DNA than isolates from different countries. Using 15 UK/American closely matched pairs of isolates that shared ancestors, we identify regions that have frequently and recently recombined to test their correlation with geographical origin. The seven genes that demonstrated the greatest clustering by geography were used in an attribution model to infer geographical origin which was tested using a further 383 UK clinical isolates to detect signatures of recent foreign travel. Patient records indicated that in 46 cases, travel abroad had occurred <2 weeks prior to sampling, and genomic analysis identified that 34 (74%) of these isolates were of a non‐UK origin. Identification of biogeographical markers in Campylobacter genomes will contribute to improved source attribution of clinical Campylobacter infection and inform intervention strategies to reduce campylobacteriosis.     

RECENT GENE‐CAPTURE ON THE UV SEX CHROMOSOMES OF THE MOSS CERATODON PURPUREUS

Sex chromosomes evolve from ordinary autosomes through the expansion and subsequent degeneration of a region of suppressed recombination that is inherited through one sex. Here we investigate the relative timing of these processes in the UV sex chromosomes of the moss Ceratodon purpureus using molecular population genetic analyses of eight newly discovered sex‐linked loci. In this system, recombination is suppressed on both the female‐transmitted (U) sex chromosome and the male‐transmitted (V) chromosome. Genes on both chromosomes therefore should show the deleterious effects of suppressed recombination and sex‐limited transmission, while purifying selection should maintain homologs of genes essential for both sexes on both sex chromosomes. Based on analyses of eight sex‐linked loci, we show that the nonrecombining portions of the U and V chromosomes expanded in at least two events (～0.6–1.3 MYA and ～2.8–3.5 MYA), after the divergence of C. purpureus from its dioecious sister species, Trichodon cylindricus and Cheilothela chloropus. Both U‐ and V‐linked copies showed reduced nucleotide diversity and limited population structure, compared to autosomal loci, suggesting that the sex chromosomes experienced more recent selective sweeps that the autosomes. Collectively these results highlight the dynamic nature of gene composition and molecular evolution on nonrecombining portions of the U and V sex chromosomes.     

Fine‐scale genetic structure in an eastern Alpine black grouse Tetrao tetrix metapopulation

Understanding genetic consequences of habitat fragmentation is crucial for the management and conservation of wildlife populations, especially in case of species sensitive to environmental changes and landscape alteration. In central Europe, the Alps are the core area of black grouse Tetrao tetrix distribution. There, black grouse dispersal is limited by high altitude mountain ridges and recent black grouse habitats are known to show some degree of natural fragmentation. Additionally, substantial anthropogenic fragmentation has occurred within the past ninety years. Facing losses of peripheral subpopulations and ongoing range contractions, we explored genetic variability and the fine‐scale genetic structure of the Alpine black grouse metapopulation at the easternmost fringe of the species’ Alpine range. Two hundred and fifty tissue samples and non‐invasive faecal and feather samples of eleven a priori defined subpopulations were used for genetic analysis based on nine microsatellite loci. Overall, eastern Alpine black grouse show similar amounts of genetic variation (H_O = 0.65, H_E = 0.66) to those found in more continuous populations like in Scandinavia. Despite of naturally and anthropogenically fragmented landscapes, genetic structuring was weak (global F_ST < 0.05), suggesting that the actual intensity of habitat fragmentation does not completely hamper dispersal, but probably restricts it to some extent. The most peripheral subpopulations at the edge of the species range show signs of genetic differentiation. The present study gives new insights into the population genetic structure of black grouse in the eastern Alps and provides a more fine‐scale view of genetic structure than previously available. Our findings will contribute to monitor the current and future status of the population under human pressures and to support supra‐regional land use planning as well as decision making processes in responsibilities of public administration.     

Isolated in an ocean of grass: low levels of gene flow between termite subpopulations

Habitat fragmentation is one of the most important causes of biodiversity loss, but many species are distributed in naturally patchy habitats. Such species are often organized in highly dynamic metapopulations or in patchy populations with high gene flow between subpopulations. Yet, there are also species that exist in stable patchy habitats with small subpopulations and presumably low dispersal rates. Here, we present population genetic data for the ‘magnetic’ termite Amitermes meridionalis, which show that short distances between subpopulations do not hinder exceptionally strong genetic differentiation (F_ST: 0.339; R_ST: 0.636). Despite the strong genetic differentiation between subpopulations, we did not find evidence for genetic impoverishment. We propose that loss of genetic diversity might be counteracted by a long colony life with low colony turnover. Indeed, we found evidence for the inheritance of colonies by so‐called ‘replacement reproductives’. Inhabiting a mound for several generations might result in loss of gene diversity within a colony but maintenance of gene diversity at the subpopulation level.     

Trans‐species variation in Dmrt1 is associated with sex determination in four European tree‐frog species

Empirical studies on the relative roles of occasional XY recombination versus sex‐chromosome turnover in preventing sex‐chromosome differentiation may shed light on the evolutionary forces acting on sex‐determination systems. Signatures of XY recombination are difficult to distinguish from those of homologous transitions (i.e., transitions in sex‐determination systems that keep sex‐chromosome identity): both models predict X and Y alleles at sex‐linked genes to cluster by species. However, the XY‐recombination model specifically predicts the reverse pattern (clustering by gametologs) for those genes that are directly involved in sex determination. Hence, the latter model can only be validated by identification of an ancestral sex‐determining region (SDR) with trans‐species polymorphism associated to sex. Here we combine a candidate‐gene approach with a genome scan to identify a small SDR shared by four species of a monophyletic clade of European tree frogs. This SDR encompasses at least the N‐terminal part of Dmrt1 and immediate upstream sequences. Our findings provide definitive evidence that sex‐chromosome homomorphy in this clade results only from XY recombination, and take an important step toward the identification of the sex‐determining locus. Moreover, the sex‐diagnostic markers we identify will enable research on environmental sex reversal in a wider range of frog species.     

Cis- and trans-acting genetic factors contribute to heterogeneity in the rate of crossing over between the Drosophila simulans clade species

In the genus Drosophila, variation in recombination rates has been found within and between species. Genetic variation for both cis‐ and trans‐acting factors has been shown to affect recombination rates within species, but little is known about the genetic factors that affect differences between species. Here, we estimate rates of crossing over for seven segments that tile across the euchromatic length of the X chromosome in the genetic backgrounds of three closely related Drosophila species. We first generated a set of Drosophila mauritiana lines each having two semidominant visible markers on the X chromosome and then introgressed these doubly marked segments into the genetic backgrounds of its sibling species, Drosophila simulans and Drosophila sechellia. Using these 21 lines (seven segments, three genetic backgrounds), we tested whether recombination rates within the doubly marked intervals differed depending on genetic background. We find significant heterogeneity among intervals and among species backgrounds. Our results suggest that a combination of both cis‐ and trans‐acting factors have evolved among the three D. simulans clade species and interact to affect recombination rate.     

Sex‐chromosome differentiation parallels postglacial range expansion in European tree frogs (Hyla arborea)

Occasional XY recombination is a proposed explanation for the sex‐chromosome homomorphy in European tree frogs. Numerous laboratory crosses, however, failed to detect any event of male recombination, and a detailed survey of NW‐European Hyla arborea populations identified male‐specific alleles at sex‐linked loci, pointing to the absence of XY recombination in their recent history. Here, we address this paradox in a phylogeographic framework by genotyping sex‐linked microsatellite markers in populations and sibships from the entire species range. Contrasting with postglacial populations of NW Europe, which display complete absence of XY recombination and strong sex‐chromosome differentiation, refugial populations of the southern Balkans and Adriatic coast show limited XY recombination and large overlaps in allele frequencies. Geographically and historically intermediate populations of the Pannonian Basin show intermediate patterns of XY differentiation. Even in populations where X and Y occasionally recombine, the genetic diversity of Y haplotypes is reduced below the levels expected from the fourfold drop in copy numbers. This study is the first in which X and Y haplotypes could be phased over the distribution range in a species with homomorphic sex chromosomes; it shows that XY‐recombination patterns may differ strikingly between conspecific populations, and that recombination arrest may evolve rapidly (<5000 generations).     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

Spatial genetic structure of a tropical understory shrub,PSYCHOTRIA OFFICINALIS (RuBIACEAE)

Analyses of fine-scale and macrogeographic genetic structure in plant populations provide an initial indication of how gene flow, natural selection, and genetic drift may collectively influence the distribution of genetic variation. The objective of our study is to evaluate the spatial dispersion of alleles within and among subpopulations of a tropical shrub, Psychotria officinalis (Rubiaceae), in a lowland wet forest in Costa Rica. This insect-pollinated, self-incompatible understory plant is dispersed primarily by birds, some species of which drop the seeds immediately while others transport seeds away from the parent plant. Thus, pollination should promote gene flow while at least one type of seed dispersal agent might restrict gene flow. Sampling from five subpopulations in undisturbed wet forest at Estación Biologíca La Selva, Costa Rica, we used electrophoretically detected isozyme markers to examine the spatial scale of genetic structure. Our goals are: 1) describe genetic diversity of each of the five subpopulations of Psychotria officinalis sampled within a contiguous wet tropical forest; 2) evaluate fine-scale genetic structure of adults of P. officinalis within a single 2.25-ha mapped plot; and 3) estimate genetic structure of P. officinalis using data from five subpopulations located up to 2 km apart. Using estimates of coancestry, statistical analyses reveal significant positive genetic correlations between individuals on a scale of 5 m but no significant genetic relatedness beyond that interplant distance within the studied subpopulation. Multilocus estimates of genetic differentiation among subpopulations were low, but significant (F_st = 0.095). Significant F_st estimates were largely attributable to a single locus (Lap-2). Thus, multilocus estimates of F_st may be influenced by microgeographic selection. If true, then the observed levels of IBD may be overestimates.