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.

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.

Microgeographic population structure of green swordail fish: genetic differentiation despite abundant migration.

Swordtails (Xiphophorus; Poeciliidae) have figured prominently in research on fish mating behaviours, sexual selection, and carcinogenesis, but their population structures and dispersal patterns have been relatively neglected. Using nine microsatellite loci, we estimated genetic differentiation in Xiphophorus helleri within and between adjacent streams in Belize. The genetic data were complemented by a tagging study of movement within one stream. In the absence of physical dispersal barriers (waterfalls), population structure followed an isolation by distance (IBD) pattern. Genetic differentiation (F(ST) up to 0.07) was significant between and within creeks, despite high dispersal in the latter as judged by the tagging data. Such heterogeneity apparently was a result of genetic drift in local demes, due to small population sizes and highly skewed paternity. The IBD pattern was interrupted by waterfalls, boosting F(ST) above 0.30 between adjacent samples across these barriers. Overall, our results are helpful in understanding the interplay of evolutionary forces and population dynamics in a small fish living in a changeable habitat.

Evolving genomic metaphors: a new look at the language of DNA.

Recent genome-sequencing efforts have confirmed that traditional "good-citizen" genes (those that encode functional RNA and protein molecules of obvious benefit to the organism) constitute only a small fraction of the genomic populace in humans and other multicellular creatures. The rest of the DNA sequence includes an astonishing collection of noncoding regions, regulatory modules, deadbeat pseudogenes, legions of repetitive elements, and hosts of oft-shifty, self-interested nomads, renegades, and immigrants. To help visualize functional operations in such intracellular genomic societies and to better encapsulate the evolutionary origins of complex genomes, new and evocative metaphors may be both entertaining and research-stimulating.

Molecular genetic relationships of the extinct dusky seaside sparrow.

Mitochondrial DNA from the extinct dusky seaside sparrow (Ammodramus maritimus nigrescens) was compared in terms of nucleotide sequence divergence to mitochondrial DNAs from extant populations of seaside sparrows. Analyses of restriction sites revealed a close phylogenetic affinity of A. m. nigrescens to other sparrow populations along the Atlantic coast of the United States but considerable genetic distance from Gulf coast birds. Concerns and applied management strategies for the seaside sparrow have been based on a morphological taxonomy that does not adequately reflect evolutionary relationships within the complex.

Balancing selection at allozyme loci in oysters: implications from nuclear RFLPs.

Population genetic analyses that depend on the assumption of neutrality for allozyme markers are used widely. Restriction fragment length polymorphisms in nuclear DNA of the American oyster evidence a pronounced population subdivision concordant with mitochondrial DNA. This finding contrasts with a geographic uniformity in allozyme frequencies previously thought to reflect high gene flow mediated by the pelagic gametes and larvae. The discordance likely is due to selection on protein electrophoretic characters that balances allozyme frequencies in the face of severe constraints to gene flow. These results raise a cautionary note for studies that rely on assumptions of neutrality for allozyme markers.

A genetic test of the natal homing versus social facilitation models for green turtle migration.

Female green turtles exhibit strong nest-site fidelity as adults, but whether the nesting beach is the natal site is not known. Under the natal homing hypothesis, females return to their natal beach to nest, whereas under the social facilitation model, virgin females follow experienced breeders to nesting beaches and after a "favorable" nesting experience, fix on that site for future nestings. Differences shown in mitochondrial DNA genotype frequency among green turtle colonies in the Caribbean Sea and Atlantic Ocean are consistent with natal homing expectations and indicate that social facilitation to nonnatal sites is rare.

Investigating sea turtle migration using DNA markers.

The past year has seen a further marshaling of genetic evidence for 'natal homing' in several species of marine turtles, a phenomenon wherein females, upon reaching sexual maturity, exhibit a propensity to return to nest at or near the respective beaches upon which they hatched some two or more decades earlier. This genetics-based inference stems from the strong spatial patterning observed in mitochondrial DNA lineages among nesting sites. Rookery-specific mitochondrial DNA markers are now being employed to monitor the natal sources of individuals captured at other phases of the life cycle, and the genetic findings have important conservation ramifications.

Application of a random walk model to geographic distributions of animal mitochondrial DNA variation.

In rapidly evolving molecules, such as animal mitochondrial DNA, mutations that delineate specific lineages may not be dispersed at sufficient rates to attain an equilibrium between genetic drift and gene flow. Here we predict conditions that lead to nonequilibrium geographic distributions of mtDNA lineages, test the robustness of these predictions and examine mtDNA data sets for consistency with our model. Under a simple isolation by distance model, the variance of an mtDNA lineage's geographic distribution is expected be proportional to its age. Simulation results indicated that this relationship is fairly robust. Analysis of mtDNA data from natural populations revealed three qualitative distributional patterns: (1) significant departure of lineage structure from equilibrium geographic distributions, a pattern exhibited in three rodent species with limited dispersal; (2) nonsignificant departure from equilibrium expectations, exhibited by two avian and two marine fish species with potentials for relatively long-distance dispersal; and (3) a progression from nonequilibrium distributions for younger lineages to equilibrium distributions for older lineages, a condition displayed by one surveyed avian species. These results demonstrate the advantages of considering mutation and genealogy in the interpretation of mtDNA geographic variation.

Molecular zoogeography of freshwater fishes in the southeastern United States.

Restriction fragment length polymorphisms in mitochondrial DNA (mtDNA) were used to reconstruct evolutionary relationships of conspecific populations in four species of freshwater fish-Amia calva, Lepomis punctatus, L. gulosus, and L. microlophus. A suite of 14-17 endonucleases was employed to assay mtDNAs from 305 specimens collected from 14 river drainages extending from South Carolina to Louisiana. Extensive mtDNA polymorphism was observed within each assayed species. In both phenograms and Wagner parsimony networks, mtDNA clones that were closely related genetically were usually geographically contiguous. Within each species, major mtDNA phylogenetic breaks also distinguished populations from separate geographic regions, demonstrating that dispersal and gene flow have not been sufficient to override geographic influences on population subdivision.-Importantly, there were strong patterns of congruence across species in the geographic placements of the mtDNA phylogenetic breaks. Three major boundary regions were characterized by concentrations of phylogenetic discontinuities, and these zones agree well with previously described zoogeographic boundaries identified by a different kind of data base-distributional limits of species-suggesting that a common set of historical factors may account for both phenomena. Repeated episodes of eustatic sea level change along a relatively static continental morphology are the likely causes of several patterns of drainage isolation and coalescence, and these are discussed in relation to the genetic data.-Overall, results exemplify the positive role that intraspecific genetic analyses may play in historical zoogeographic reconstruction. They also point out the potential inadequacies of any interpretations of population genetic structure that fail to consider the influences of history in shaping that structure.

A genetic discontinuity in a continuously distributed species: mitochondrial DNA in the American oyster, Crassostrea virginica.

Restriction site variation in mitochondrial DNA (mtDNA) of the American oyster (Crassostrea virginica) was surveyed in continuously distributed populations sampled from the Gulf of St. Lawrence, Canada, to Brownsville, Texas. mtDNA clonal diversity was high, with 82 different haplotypes revealed among 212 oysters with 13 endonucleases. The mtDNA clones grouped into two distinct genetic arrays (estimated to differ by about 2.6% in nucleotide sequence) that characterized oysters collected north vs. south of a region on the Atlantic mid-coast of Florida. The population genetic "break" in mtDNA contrasts with previous reports of near uniformity of nuclear (allozyme) allele frequencies throughout the range of the species, but agrees closely with the magnitude and pattern of mtDNA differentiation reported in other estuarine species in the southeastern United States. This concordance of mtDNA phylogenetic pattern across independently evolving species provides strong evidence for vicariant biogeographic processes in initiating intraspecific population structure. The post-Miocene ecological history of the region suggests that reduced precipitation levels in an enlarged Floridian peninsula may have created discontinuities in suitable estuarine habitat for oysters during glacial periods, and that today such population separations are maintained by the combined influence of ecological gradients and oceanic currents on larval dispersal. The results are consistent with the hypothesis that historical vicariant events, in conjunction with contemporary environmental influences on gene flow, can result in genetic discontinuities in continuously distributed species with high dispersal capability.

Genetic change and rates of cladogenesis.

Models are introduced which predict ratios of mean levels of genetic divergence in species-rich versus species-poor phylads under two competing assumptions: (1) genetic differentiation is a function of time, unrelated to the number of cladogenetic events and (2) genetic differentiation is proportional to the number of speciation events in the group. The models are simple, general, and biologically real, but not precise. They lead to qualitatively distinct predictions about levels of genetic divergence depending upon the relationship between rates of speciation and amount of genetic change. When genetic distance between species is a function of time, mean genetic distances in speciose and depauperate phylads of equal evolutionary age are very similar. On the contrary, when genetic distance is a function of the number of speciations in the history of a phylad, the ratio of mean genetic distances separating species in speciose versus depauperate phylads is greater than one, and increases rapidly as the frequency of speciations in one group relative to the other increases. The models may be tested with data from natural populations to assess (1) possible correlations between rates of anagenesis and cladogenesis and (2) the amount of genetic differentiation accompanying the speciation process. The data collected in electrophoretic surveys and other kinds of studies can be used to test the predictions of the models. For this purpose genetic distances need to be measured in speciose and depauperate phylads of equal evolutionary age. The limited information presently available agrees better with the model predicting that genetic change is primarily a function of time, and is not correlated with rates of speciation. Further testing of the models is, however, required before firm conclusions can be drawn.

An assessment of "hidden" heterogeneity within electromorphs at three enzyme loci in deer mice.

Allelic heterogeneity within protein electromorphs at three loci was examined in populations of deer mice (Peromyscus maniculatus) collected from five localities across North America. We used a variety of electrophoretic techniques (including several starch and acrylamide conditions, gel-sieving, and isoelectric focusing) plus heat denaturation. Of particular interest was the supernatant glutamate oxalate transaminase system (GOT-1; aspartate aminotransferase-1 of some authors), which under standard electrophoretic conditions had been shown to exhibit basically a two-allele polymorphism throughout the range of maniculatus. The use of all of the above techniques failed to uncover any additional variation for GOT-1 in these populations. Similarly, no new scorable variation was resolved at the essentially monomorphic malate dehydrogenase-1 locus by additional conditions of electrophoresis. In marked contrast to the results for the above two enzymes, the use of multiple conditions of electrophoresis resolved the 8 standard-condition electromorphs of esterase-1 into a total of 23 variants showing strong geographic differentiation in frequency. These 23 electromorphs were further divided into a total of 35 variants by thermal stability studies. However, the allelic nature of all of the thermal stability esterase variants remains to be documented. The results of this study, taken together with the remarkable geographic heterogeneity for this species in ecology, morphology, karyotype and mitochondrial DNA sequence, suggest that some form of balancing selection may be acting to maintain the GOT-1 polymorphism.

Hybridization and introgression among species of sunfish (Lepomis): analysis by mitochondrial DNA and allozyme markers.

We explore the potential of mitochondrial DNA (mtDNA) analysis, alone and in conjunction with allozymes, to study low-frequency hybridization and introgression phenomena in natural populations. MtDNAs from small samples of nine species of sunfish (Lepomis, Centrarchidae) were purified and digested with each of 13 informative restriction enzymes. Digestion profiles for all species were highly distinct: estimates of overall fragment homology between pairs of species ranged from 0-36%. Allozymes encoded by nine nuclear genes also showed large freqency differences among species and together with mtDNA provided many genetic markers for hybrid identification. A genetic analysis of 277 sunfish from two locations in north Georgia revealed the following: (1) a low frequency of interspecific hybrids, all of which appeared to be F1's; (2) the involvement of five sympatric Lepomis species in the production of these hybrids; (3) no evidence for introgression between species in our study locales (although for rare hybridization, most later-generation backcrosses would not be reliably distinguished from parentals); (4) a tendency for hybridizations to take place preferentially between parental species differing greatly in abundance; (5) a tendency for the rare species in a hybrid cross to provide the female parent. Our data suggest that absence of conspecific pairing partners and mating stimuli for females of rarer species may be important factors in increasing the likelihood of interspecific hybridization. The maternal inheritance of mtDNA offers at least two novel advantages for hybridization analysis: (1) an opportunity to determine direction in hybrid crosses; and (2) due to the linkage among mtDNA markers, an increased potential to distinguish effects of introgression from symplesiomorphy or character convergence.

Definition and properties of disequilibrium statistics for associations between nuclear and cytoplasmic genotypes.

We define and establish the interrelationships of four components of statistical association between a diploid nuclear gene and a uniparentally transmitted, haploid cytoplasmic gene: an allelic (gametic) disequilibrium (D), which measures associations between alleles at the two loci; and three genotypic disequilibria (D1, D2, D3), which measure associations between two cytotypes and the three respective nuclear backgrounds. We also consider an alternative set of measures, including D and the residual disequilibrium (d). The dynamics of these disequilibria are then examined under three conventional models of the mating system: (1) random mating; (2a) assortative mating without dominance (the "mixed-mating model"); and (2b) assortative mating with dominance ("O'Donald's model"). The trajectories of gametic disequilibria are similar to those for pairs of unlinked nuclear loci. The dynamics of genotypic disequilibria exhibit a variety of behaviors depending on the model and the initial conditions. Procedures for statistical estimation of cytonuclear disequilibria are developed and applied to several real and hypothetical data sets. Special attention is paid to the biological interpretations of various categories of allelic and genotypic disequilibria in hybrid zones. Genetic systems for which these statistics might be appropriate include nuclear genotype frequencies in conjunction with those for mitochondrial DNA, chloroplast DNA, or cytoplasmically inherited microorganisms.

The effects of assortative mating and migration on cytonuclear associations in hybrid zones.

We examine the influence of nonrandom mating and immigration on the evolutionary dynamics of cytonuclear associations in hybrid zones. Recursion equations for allelic and genotypic cytonuclear disequilibria were generated under models of (1) migration alone, assuming hybrid zone matings are random with respect to cytonuclear genotype; and (2) migration in conjunction with refined epistatic mating, in which females of the pure parental species preferentially mate with conspecific males. Major results are as follows: (a) even the slightest migration removes the dependency of the final outcome on initial conditions, producing a unique equilibrium in which both pure parental genotypes are maintained in the hybrid zone; (b) in contrast to nuclear genes, the dynamics of cytoplasmic allele frequencies appear robust to changes in the assumed mating system, yet are particularly sensitive to gene flow; (c) continued immigration can generate permanent cytonuclear disequilibria, whether mating is random or assortative; and (d) the order of population censusing (before versus after reproduction by immigrants) can have a dramatic effect on the magnitude but not the pattern of cytonuclear disequilibria. Using the maximum likelihood method, the parameter space of migration rates and assortative mating rates was examined for best fit to observed cytonuclear disequilibria data in a hybrid population of Hyla tree frogs. An epistatic mating model with a total immigration rate of about 32% per generation produces equilibrium gene frequencies and cytonuclear disequilibria consistent with the empirical observations.

Genetic parentage in large half-sib clutches: theoretical estimates and empirical appraisals.

Nearly all of the 906 embryos from a male-tended nest of the sand goby (Pomatoschistus minutus) were genotyped at two hypervariable microsatellite loci to document conclusively the number of mothers and their relative genetic contributions to the nest. The true number of mothers determined by this nearly exhaustive genetic appraisal was compared to computer simulation treatments based on allele frequencies in the population, assumptions about reproductive skew, and statistical sampling strategies of progeny subsets. The "ground-truthed" appraisal and the theoretical estimates showed good agreement, indicating that for this nest a random sample of approximately 20 offspring would have sufficed for assessing the true number of biological parents (but not necessarily their relative genetic contributions). Also, a general dilocus matrix procedure is suggested for organizing and interpreting otherwise cumbersome data sets when extremely large numbers of full-sib and half-sib embryos from a nest are genotyped at two or more hypervariable loci.

The use of restriction endonucleases to measure mitochondrial DNA sequence relatedness in natural populations. I. Population structure and evolution in the genus Peromyscus.

In this study we introduce to natural population analysis a molecular technique that involves the use of restriction endonucleases to compare mitochondrial DNA (mtDNA) sequences. We have examined the fragment patterns produced by six restriction endonucleases acting upon mtDNA isolated from 23 samples of three species of the rodent Peromyscus. Our observations confirm the following conclusions derived from previous experiments with laboratory animals: (1) mtDNA within an individual homogeneous; (2) at least the majority of mtDNA present in an individual is inherited from the female parent. Our experiments demonstrate for the first time that there is detectable heterogeneity in mtDNA sequences within and among natural geographic populations of a species and that this heterogeneity can readily be used to estimate relatedness between individuals and populations. Individuals collected within a single locale show less than 0.5% sequence divergence, while those collected from conspecific populations separated by 50 ti 500 miles differ by approximately 1.5%. The mtDNAs of the closely related sibling species P. polionotus and P. maniculatus differ from each other by 13 to 17%; nonsibling species differ by more than 20%. Qualitative and quantitative approaches to analysis of digestion patterns are suggested. The results indicate that restriction analysis of mtNDA may become the most sensitive and powerful technique yet available for reconstructing evolutionary relationships among conspecific organisms.

Molecular evidence for multiple origins of hybridogenetic fish clones (Poeciliidae:Poeciliopsis).

Hybrid matings between the sexual species Poeciliopsis monacha and Poeciliopsis lucida produced a series of diploid all-female lineages of P. monacha-lucida that inhabit the Río Fuerte of northwestern Mexico. Restriction site analyses of mitochondrial DNA (mtDNA) clearly revealed that P. monacha was the maternal ancestor of these hybrids. The high level of mtDNA diversity in P. monacha was mirrored by similarly high levels in P. monacha-lucida; thus hybridizations giving rise to unisexual lineages have occurred many times. However, mtDNA variability among P. monacha-lucida lineages revealed a geographical component. Apparently the opportunity for the establishment of unisexual lineages varies among tributaries of the Río Fuerte. We hypothesize that a dynamic complex of sexual and clonal fishes appear to participate in a feedback process that maintains genetic diversity in both the sexual and asexual components.

Genetic Variation and Geographic Differentiation in Mitochondrial DNA of the Horseshoe Crab, LIMULUS POLYPHEMUS.

Restriction site variation in mitochondrial DNA (mtDNA) of the horseshoe crab (Limulus polyphemus) was surveyed in populations ranging from New Hampshire to the Gulf Coast of Florida. MtDNA clonal diversity was moderately high, particularly in southern samples, and a major genetic "break" (nucleotide sequence divergence approximately 2%) distinguished all sampled individuals which were north vs. south of a region in northeastern Florida. The area of genotypic divergence in Limulus corresponds to a long-recognized zoogeographic boundary between warm-temperate and tropical marine faunas, and it suggests that selection pressures and/or gene flow barriers associated with water mass differences may also influence the evolution of species widely distributed across such transition zones. On the other hand, a comparison of the mtDNA divergence patterns in Limulus with computer models involving stochastic lineage extinction in species with limited gene flow demonstrates that deterministic explanations need not necessarily be invoked to account for the observations. Experiments to distinguish stochastic from deterministic possibilities are suggested. Overall, the pattern and magnitude of mtDNA differentiation in horseshoe crabs is very similar to that typically reported for freshwater and terrestrial species assayed over a comparable geographic range. Results demonstrate for the first time that, geographically, at least some continuously distributed marine organisms can show considerable mtDNA genetic differentiation.