Persistence of Mhc heterozygosity in homozygous clonal killifish,Rivulus marmoratus: implications for the origin of hermaphroditism

The mangrove killifish Rivulus marmoratus, a neotropical fish in the order Cyprinodontiformes, is the only known obligatorily selfing, synchronous hermaphroditic vertebrate. To shed light on its population structure and the origin of hermaphroditism, major histocompatibility complex (Mhc) class I genes of the killifish from seven different localities in Florida, Belize, and the Bahamas were cloned and sequenced. Thirteen loci and their alleles were identified and classified into eight groups. The loci apparently arose approximately 20 million years ago (MYA) by gene duplications from a single common progenitor in the ancestors of R. marmoratus and its closest relatives. Distinct loci were found to be restricted to different populations and different individuals in the same population. Up to 44% of the fish were heterozygotes at Mhc loci, as compared to near homozygosity at non-Mhc loci. Large genetic distances between some of the Mhc alleles revealed the presence of ancestral allelic lineages. Computer simulation designed to explain these findings indicated that selfing is incomplete in R. marmoratus populations, that Mhc allelic lineages must have diverged before the onset of selfing, and that the hermaphroditism arose in a population containing multiple ancestral Mhc lineages. A model is proposed in which hermaphroditism arose stage-wise by mutations, each of which spread through the entire population and was fixed independently in the emerging clones.     

A mixed-mating strategy in a hermaphroditic vertebrate

Mixed-mating systems, in which hermaphrodites can either self-fertilize or outcross, are common in many species of plants and invertebrates and have been informative models for studying the selective forces that can maintain both inbreeding and outbreeding in populations. Here, we document a remarkable instance of evolutionary convergence to an analogous mixed mating system by a vertebrate, the mangrove killifish (Kryptolebias marmoratus). In this androdioecious species, most individuals are simultaneous hermaphrodites that characteristically self-fertilize, resulting in local populations that consist of (nearly) homozygous lines. Most demes are also genetically diverse, an observation traditionally attributed to de novo mutation coupled with high levels of inter-site migration. However, data presented here, from a survey of 35 microsatellite loci in Floridian populations, show that genotypic diversity also stems proximally from occasional outcross events that release 'explosions' of transient recombinant variation. The result is a local population genetic pattern (of extensive genotypic variety despite low but highly heterogeneous intra-individual heterozygosities) that differs qualitatively from the genetic architectures known in any other vertebrate species. Advantages of a mixed-mating strategy in K. marmoratus probably relate to this fish's solitary lifestyle and its ability to colonize new habitats.     

Abundant males in populations of a selfing hermaphrodite fish, Rivulus marmoratus, from some Belize cays

The killifish Rivulus marmoratus is the only known selfing hermaphroditic vertebrate, and males of the species are usually exceedingly rare or non-existent in nature. Collections on several Belize cays in 1988 and 1989 yielded 13.5–24% males. Factors responsible for this unusually high proportion of males are currently not understood. Likewise, the biological significance of males in populations which otherwise consist of selfing hermaphrodites (with internal fertilization) is problematic and awaits further study.     

Association behaviour of the self-fertilizing Kryptolebias marmoratus (Poey): the influence of microhabitat use on the potential for a complex mating system

Preliminary observations were conducted to identify conspicuous body postures and movements of males and hermaphrodites in the mangrove killifish Kryptolebias marmoratus . These behaviours were used quantitatively to examine the social interactions for experimental pairings of K. marmoratus of different sexual states ( i.e. simultaneous hermaphrodite and male) in an aquarium with an open-water area and simulated crab burrows. This allowed observation of behaviours that could not be observed in the field. Kryptolebias marmoratus , regardless of sexual state and experimental treatment, spent 40% of the time in the burrow. Hermaphrodites exhibited a preference for associating with males rather than other hermaphrodites. The observed complexity of displayed behaviours and interactions between paired conspecifics indicate that K. marmoratus has a rich repertoire of social behaviour not predicted for a strictly selfing species. Also, land crab burrows play an important role in their social interactions.     

Environmental diel variation, parasite loads, and local population structuring of a mixed-mating mangrove fish

Genetic variation within populations depends on population size, spatial structuring, and environmental variation, but is also influenced by mating system. Mangroves are some of the most productive and threatened ecosystems on earth and harbor a large proportion of species with mixed-mating (self-fertilization and outcrossing). Understanding population structuring in mixed-mating species is critical for conserving and managing these complex ecosystems. Kryptolebias marmoratus is a unique mixed-mating vertebrate inhabiting mangrove swamps under highly variable tidal regimes and environmental conditions. We hypothesized that geographical isolation and ecological pressures influence outcrossing rates and genetic diversity, and ultimately determine the local population structuring of K. marmoratus. By comparing genetic variation at 32 microsatellites, diel fluctuations of environmental parameters, and parasite loads among four locations with different degrees of isolation, we found significant differences in genetic diversity and genotypic composition but little evidence of isolation by distance. Locations also differed in environmental diel fluctuation and parasite composition. Our results suggest that mating system, influenced by environmental instability and parasites, underpins local population structuring of K. marmoratus. More generally, we discuss how the conservation of selfing species inhabiting mangroves and other biodiversity hotspots may benefit from knowledge of mating strategies and population structuring at small spatial scales.     

Population genetics and evolution of the mangrove rivulus Kryptolebias marmoratus,the world's only self‐fertilizing hermaphroditic vertebrate

The mangrove rivulus, Kryptolebias marmoratus (Rivulidae, Cyprinodontiformes), is phylogenetically embedded within a large clade of oviparous (egg laying) and otherwise mostly gonochoristic (separate sex) killifish species in the circumtropical suborder Aplocheiloidei. It is unique in its reproductive mode: K. marmoratus is essentially the world's only vertebrate species known to engage routinely in self‐fertilization as part of a mixed‐mating strategy of selfing plus occasional outcrossing with gonochoristic males. This unique form of procreation has profound population‐genetic and evolutionary‐genetic consequences that are the subject of this review.     

BEST OF BOTH WORLDS? ASSOCIATION BETWEEN OUTCROSSING AND PARASITE LOADS IN A SELFING FISH

Mixed-mating strategies (i.e., intermediate levels of self-fertilization and outcrossing in hermaphrodites) are relatively common in plants and animals, but why self-fertilization (selfing) rates vary so much in nature has proved difficult to explain. We tested the hypothesis that parasites help maintain mixed-mating using a partially selfing fish (Kryptolebias marmoratus) as a model. We show that outcrossed progeny in the wild are genetically more diverse and less susceptible to multiple parasite infections than their selfed counterparts. Given that outcrossing in K. marmoratus can only be attained by male-hermaphrodite matings, our data provide an explanation for the coexistence of males and hermaphrodites in androdioecious species where hermaphrodites are unable to outcross among themselves. Moreover, our study provides evidence that parasites contribute to maintaining mixed-mating in a natural animal population.     

Epigenetic regulation of sex ratios may explain natural variation in self-fertilization rates

Self-fertilization (selfing) favours reproductive success when mate availability is low, but renders populations more vulnerable to environmental change by reducing genetic variability. A mixed-breeding strategy (alternating selfing and outcrossing) may allow species to balance these needs, but requires a system for regulating sexual identity. We explored the role of DNA methylation as a regulatory system for sex-ratio modulation in the mixed-mating fish Kryptolebias marmoratus. We found a significant interaction between sexual identity (male or hermaphrodite), temperature and methylation patterns when two selfing lines were exposed to different temperatures during development. We also identified several genes differentially methylated in males and hermaphrodites that represent candidates for the temperature-mediated sex regulation in K. marmoratus. We conclude that an epigenetic mechanism regulated by temperature modulates sexual identity in this selfing species, providing a potentially widespread mechanism by which environmental change may influence selfing rates. We also suggest that K. marmoratus, with naturally inbred populations, represents a good vertebrate model for epigenetic studies.     

Microsatellite documentation of male-mediated outcrossing between inbred laboratory strains of the self-fertilizing mangrove killifish (Kryptolebias Marmoratus)

Primers for 36 microsatellite loci were developed and employed to characterize genetic stocks and detect possible outcrossing between highly inbred laboratory strains of the self-fertilizing mangrove killifish, Kryptolebias marmoratus. From attempted crosses involving hermaphrodites from particular geographic strains and gonochoristic males from others, 2 among a total of 32 surveyed progenies (6.2%) displayed multilocus heterozygosity clearly indicative of interstrain gametic syngamy. One of these outcross hybrids was allowed to resume self-fertilization, and microsatellite assays of progeny showed that heterozygosity decreased by approximately 50% after one generation, as expected. Although populations of K. marmoratus consist mostly of synchronous hermaphrodites with efficient mechanisms of internal self-fertilization, these laboratory findings experimentally confirm that conspecific males can mediate occasional outcross events and that this process can release extensive genic heterozygosity.     

How multilocus genotypic pattern helps to understand the history of selfing populations: a case study in Medicago truncatula

The occurrence of populations exhibiting high genetic diversity in predominantly selfing species remains a puzzling question, since under regular selfing genetic diversity is expected to be depleted at a faster rate than under outcrossing. Fine-scale population genetics approaches may help to answer this question. Here we study a natural population of the legume Medicago truncatula in which both the fine-scale spatial structure and the selfing rate are characterized using three different methods. Selfing rate estimates were very high ( approximately 99%) irrespective of the method used. A clear pattern of isolation by distance reflecting small seed dispersal distances was detected. Combining genotypic data over loci, we could define 34 multilocus genotypes. Among those, six highly inbred genotypes (lines) represented more than 75% of the individuals studied and harboured all the allelic variation present in the population. We also detected a large set of multilocus genotypes resembling recombinant inbred lines between the most frequent lines occurring in the population. This finding illustrates the importance of rare recombination in redistributing available allelic diversity into new genotypic combinations. This study shows how multilocus and fine-scale spatial analyses may help to understand the population history of self-fertilizing species, especially to make inferences about the relative role of foundation/migration and recombination events in such populations.     

Spatial effects and rare outcrossing events in Medicago truncatula (Fabaceae)

In order to elucidate the mechanisms underlying the large amount of RAPD polymorphism found in 1990 in a population of the selfing annual Medicago truncatula GAERTN. (Fabaceae), we have analysed most of the individuals (n = 363) from the same population 6 years later using microsatellite loci. We confirm the result of the earlier study, namely that this population is very polymorphic and highly subdivided, with approximately 37% of the variance distributed among subpopulations, only 50 m apart one from another. We use standard F-statistics analyses, linkage disequilibria, minimum spanning network, multilocus assignment tests and spatial autocorrelation analyses to test the hypotheses that spatial structure and outcrossing events are involved in maintaining the large amount of genetic diversity at the level of each subpopulation. Interestingly, fine-scale spatial structure could be observed in only one subpopulation suggesting that other mechanisms are acting elsewhere. To the best of our knowledge, this is the first study of fine spatial genetic structure in a predominantly selfing species.     

A Markov chain Monte Carlo approach for joint inference of population structure and inbreeding rates from multilocus genotype data

Nonrandom mating induces correlations in allelic states within and among loci that can be exploited to understand the genetic structure of natural populations (Wright 1965). For many species, it is of considerable interest to quantify the contribution of two forms of nonrandom mating to patterns of standing genetic variation: inbreeding (mating among relatives) and population substructure (limited dispersal of gametes). Here, we extend the popular Bayesian clustering approach STRUCTURE (Pritchard et al. 2000) for simultaneous inference of inbreeding or selfing rates and population-of-origin classification using multilocus genetic markers. This is accomplished by eliminating the assumption of Hardy-Weinberg equilibrium within clusters and, instead, calculating expected genotype frequencies on the basis of inbreeding or selfing rates. We demonstrate the need for such an extension by showing that selfing leads to spurious signals of population substructure using the standard STRUCTURE algorithm with a bias toward spurious signals of admixture. We gauge the performance of our method using extensive coalescent simulations and demonstrate that our approach can correct for this bias. We also apply our approach to understanding the population structure of the wild relative of domesticated rice, Oryza rufipogon, an important partially selfing grass species. Using a sample of n = 16 individuals sequenced at 111 random loci, we find strong evidence for existence of two subpopulations, which correlates well with geographic location of sampling, and estimate selfing rates for both groups that are consistent with estimates from experimental data (s approximately 0.48-0.70).     

Extreme homogeneity and low genetic diversity in Kryptolebias ocellatus from south-eastern Brazil suggest a recent foundation for this androdioecious fish population

This study documents unexpectedly low levels of intra and interpopulation genetic diversity in Kryptolebias ocellatus, an androdioecious and predominantly self-fertilizing killifish from south-eastern Brazil. This finding generally is inconsistent with the established opinion that the K. ocellatus and K. marmoratus clade originated in this geographic region and later dispersed northward into the Caribbean.     

Contrasting demographic history and population structure in Capsella rubella and Capsella grandiflora, two closely related species with different mating systems

Both mating system and population history can have large impacts on genetic diversity and population structure. Here, we use multilocus sequence data to investigate how these factors impact two closely related Brassicaceae species: the selfing Capsella rubella and the outcrossing C. grandiflora. To do this, we have sequenced 16 loci in approximately 70 individuals from 7 populations of each species. Patterns of population structure differ strongly between the two species. In C. grandiflora, we observe an isolation-by-distance pattern and identify three clearly delineated genetic groups. In C. rubella, where we estimate the selfing rate to be 0.90-0.94, the pattern is less clear with some sampling populations forming separate genetic clusters while others are highly mixed. The two species also have divergent histories. Our analysis gives support for a bottleneck approximately 73 kya (20-139 kya) in C. rubella, which most likely represents speciation from C. grandiflora. In C. grandiflora, there is moderate support for the standard neutral model in 2 of 3 genetic clusters, while the third cluster and the total data set show evidence of expansion. It is clear that mating system has an impact on these two species, for example affecting the level of genetic variation and the genetic structure. However, our results also clearly show that a combination of past and present processes, some of which are not affected by mating system, is needed to explain the differences between C. rubella and C. grandiflora.     

Consequences of multiple mating‐system shifts for population and range‐wide genetic structure in a coastal dune plant

Evolutionary transitions from outcrossing to selfing can strongly affect the genetic diversity and structure of species at multiple spatial scales. We investigated the genetic consequences of mating‐system shifts in the North American, Pacific coast dune endemic plant Camissoniopsis cheiranthifolia (Onagraceae) by assaying variation at 13 nuclear (n) and six chloroplast (cp) microsatellite (SSR) loci for 38 populations across the species range. As predicted from the expected reduction in effective population size (N_e) caused by selfing, small‐flowered, predominantly selfing (SF) populations had much lower nSSR diversity (but not cpSSR) than large‐flowered, predominantly outcrossing (LF) populations. The reduction in nSSR diversity was greater than expected from the effects of selfing on N_e alone, but could not be accounted for by indirect effects of selfing on population density. Although selfing should reduce gene flow, SF populations were not more genetically differentiated than LF populations. We detected five clusters of nSSR genotypes and three groups of cpSSR haplotypes across the species range consisting of parapatric groups of populations that usually (but not always) differed in mating system, suggesting that selfing may often initiate ecogeographic isolation. However, lineage‐wide genetic variation was not lower for selfing clusters, failing to support the hypothesis that selection for reproductive assurance spurred the evolution of selfing in this species. Within three populations where LF and SF plants coexist, we detected genetic differentiation among diverged floral phenotypes suggesting that reproductive isolation (probably postzygotic) may help maintain the striking mating‐system differentiation observed across the range of this species.     

Characterization of 24 microsatellite markers in Primula chungensis (Primulaceae), a distylous homostylous species,using MiSeq sequencing

Primula chungensis is a species with considerable floral and mating-system variation, including distylous (outcrossing), homostylous (selfing) and mixed populations that contain both outcrossing and selfing forms. We isolated 24 microsatellite markers from P.chungensis using Illumina MiSeq sequencing. Polymorphism and genetic diversity were then measured based on a sample of 24 individuals from a natural population in southern Tibet. All loci were polymorphic with the number of alleles per locus ranging from 2 to4. The observed and expected heterozygosity ranged from 0 to 1 and 0.219 to 0.708, respectively. The microsatellite markers we have identified will serve as valuable tools for the investigation of the population genetic structure and phylogeography of P.chungensis and will inform models of the evolutionary history of mating systems in the species.     

Characterization of 24 microsatellite markers in Primula chungensis(Primulaceae),a distylous-homostylous species,using MiSeq sequencing

Primula chungensis is a species with considerable floral and mating-system variation,including distylous(outcrossing),homostylous(selfing) and mixed populations that contain both outcrossing and selfing forms.We isolated 24 microsatellite markers from P.chungensis using Illumina Mi Seq sequencing.Polymorphism and genetic diversity were then measured based on a sample of 24 individuals from a natural population in southern Tibet.All loci were polymorphic with the number of alleles per locus ranging from 2 to 4.The observed and expected heterozygosity ranged from 0 to 1 and 0.219 to 0.708,respectively.The microsatellite markers we have identified will serve as valuable tools for the investigation of the population genetic structure and phylogeography of P.chungensis and will inform models of the evolutionary history of mating systems in the species.     

Animals mix it up too: the distribution of self-fertilization among hermaphroditic animals

Excluding insects, hermaphroditism occurs in about one-third of animal species, providing numerous opportunities for the evolution of selfing. Here we provide an overview of reproductive traits in hermaphroditic animal species, review the distribution of selfing rates in animals, and test for ecological correlates of selfing. Our dataset (1342 selfing-rate estimates for 142 species) is 97% based on estimates derived from the analysis of population structure (F(IS)-estimates) using genetic markers. The distribution of selfing is slightly U-shaped and differs significantly from the more strongly U-shaped plant distribution with 47% of animal t-estimates being intermediate (falling between 0.2 and 0.8) compared to 42% for plants. The influence of several factors on the distribution of selfing rates was explored (e.g., number of populations studied per species, habitat, coloniality, sessility, or fertilization type), none of which significantly affect the distribution. Our results suggest that genetic forces might contribute to the evolution of self-fertilization to the same extent in animals and plants, although the high proportion of intermediate outcrossing suggests a significant role of ecological factors (e.g., reproductive assurance) in animals. However, we caution that the distribution of selfing rates in animals is affected by various factors that might bias F(IS)-estimates, including phylogenetic underrepresentation of highly selfing and outcrossing species, various genotyping errors (e.g., null alleles) and inbreeding depression. This highlights the necessity of obtaining better estimates of selfing for hermaphroditic animals, such as genotyping progeny arrays, as in plants.     

A metapopulation perspective on genetic diversity and differentiation in partially self-fertilizing plants

Abstract.— Partial self-fertilization is common in higher plants. Mating system variation is known to have important consequences for how genetic variation is distributed within and among populations. Selfing is known to reduce effective population size, and inbreeding species are therefore expected to have lower levels of genetic variation than comparable out crossing taxa. However, several recent empirical studies have shown that reductions in genetic diversity within populations of inbreeding species are far greater than the expected reductions based on the reduced effective population size. Two different processes have been argued to cause these patterns, selective sweeps (or hitchhiking) and background selection. Both are expected to be most effective in reducing genetic variation in regions of low recombination rates. Selfing is known to reduce the effective recombination rate, and inbreeding taxa are thus thought to be particularly vulnerable to the effects of hitchhiking or background selection. Here I propose a third explanation for the lower-than-expected levels of genetic diversity within populations of selfing species; recurrent extinctions and recolonizations of local populations, also known as metapopulation dynamics. I show that selfing in a metapopulation setting can result in large reductions in genetic diversity within populations, far greater than expected based the lower effective population size inbreeding species is expected to have. The reason for this depends on an interaction between selfing and pollen migration.     

Migration without interbreeding: Evolutionary history of a highly selfing Mediterranean grass inferred from whole genomes

Wild plant populations show extensive genetic subdivision and are far from the ideal of panmixia which permeates population genetic theory. Understanding the spatial and temporal scale of population structure is therefore fundamental for empirical population genetics – and of interest in itself, as it yields insights into the history and biology of a species. In this study we extend the genomic resources for the wild Mediterranean grass Brachypodium distachyon to investigate the scale of population structure and its underlying history at whole-genome resolution. A total of 86 accessions were sampled at local and regional scales in Italy and France, which closes a conspicuous gap in the collection for this model organism. The analysis of 196 accessions, spanning the Mediterranean from Spain to Iraq, suggests that the interplay of high selfing and seed dispersal rates has shaped genetic structure in B. distachyon. At the continental scale, the evolution in B. distachyon is characterized by the independent expansion of three lineages during the Upper Pleistocene. Today, these lineages may occur on the same meadow yet do not interbreed. At the regional scale, dispersal and selfing interact and maintain high genotypic diversity, thus challenging the textbook notion that selfing in finite populations implies reduced diversity. Our study extends the population genomic resources for B. distachyon and suggests that an important use of this wild plant model is to investigate how selfing and dispersal, two processes typically studied separately, interact in colonizing plant species.