Variable rates of hybridization among contact zones between a pair of topminnow species, Fundulus notatus and F. olivaceus.

Pairs of species that exhibit broadly overlapping distributions, and multiple geographically isolated contact zones, provide opportunities to investigate the mechanisms of reproductive isolation. Such naturally replicated systems have demonstrated that hybridization rates can vary substantially among populations, raising important questions about the genetic basis of reproductive isolation. The topminnows, Fundulus notatus and F. olivaceus, are reciprocally monophyletic, and co-occur in drainages throughout much of the central and southern United States. Hybridization rates vary substantially among populations in isolated drainage systems. We employed genome-wide sampling to investigate geographic variation in hybridization, and to assess the possible importance of chromosome fusions to reproductive isolation among nine separate contact zones. The species differ by chromosomal rearrangements resulting from Robertsonian (Rb) fusions, so we hypothesized that Rb fusion chromosomes would serve as reproductive barriers, exhibiting steeper genomic clines than the rest of the genome. We observed variation in hybridization dynamics among drainages that ranged from nearly random mating to complete absence of hybridization. Contrary to predictions, our use of genomic cline analyses on mapped species-diagnostic SNP markers did not indicate consistent patterns of variable introgression across linkage groups, or an association between Rb fusions and genomic clines that would be indicative of reproductive isolation. We did observe a relationship between hybridization rates and population phylogeography, with the lowest rates of hybridization tending to be found in populations inferred to have had the longest histories of drainage sympatry. Our results, combined with previous studies of contact zones between the species, support population history as an important factor in explaining variation in hybridization rates.

Invasion genetics of the mummichog (Fundulus heteroclitus): recent anthropogenic introduction in Iberia.

Human activities such as trade and transport have increased considerably in the last decades, greatly facilitating the introduction and spread of non-native species at a global level. In the Iberian Peninsula, Fundulus heteroclitus, a small euryhaline coastal fish with short dispersal, was found for the first time in the mid-1970s. Since then, F. heteroclitus has undergone range expansions, colonizing the southern region of Portugal, southwestern coast of Spain and the Ebro Delta in the Mediterranean Sea. Cytochrome b sequences were used to elucidate the species invasion pathway in Iberia. Three Iberian locations (Faro, Cádiz and Ebro Delta) and 13 other locations along the native range of F. heteroclitus in North America were sampled. Results revealed a single haplotype, common to all invasive populations, which can be traced to the northern region of the species' native range. We posit that the origin of the founder individuals is between New York and Nova Scotia. Additionally, the lack of genetic structure within Iberia is consistent with a recent invasion scenario and a strong founder effect. We suggest the most probable introduction vector is associated with the aquarium trade. We further discuss the hypothesis of a second human-mediated introduction responsible for the establishment of individuals in the Ebro Delta supported by the absence of adequate muddy habitats linking Cádiz and the Ebro Delta. Although the species has a high tolerance to salinity and temperature, ecological niche modelling indicates that benthic habitat constraints prevent along-shore colonisation suggesting that such expansions would need to be aided by human release.

Phylogenomic analysis of Fundulidae (Teleostei: Cyprinodotiformes) using RNA-sequencing data.

Fishes of the New World cyprinodontiform family Fundulidae display a wide variety of tolerance to environmental conditions, making them a valuable model system for comparative, evolutionary, and environmental studies. Despite numerous attempts to resolve the phylogenetic relationships of family Fundulidae, the basal structure of the phylogeny remains unresolved. The lack of a robust and fully resolved phylogeny for family Fundulidae and its most speciose genus Fundulus is an impediment to future research. This study utilized novel RNA-sequencing data for phylogenetic inference among16 members of Fundulidae to better refine the basal nodes of the family and confront long-standing questions regarding (1) the monophyletic status of genus Fundulus, and validity of the Lucania and recently synonymized Adinia genera; (2) the relationship of the west coast endemic Fundulus parvipinnis and F. lima to other Fundulus species; and (3) the validity of subgeneric classifications. In addition, previously published nuclear gene sequences for 32 Fundulidae species were re-analyzed in combination with novel RNA-sequencing data. Maximum likelihood and Bayesian analyses generated identical phylogenies with strong statistical support at nearly all nodes, demonstrating the utility of RNA-sequencing data in constructing robust phylogenies not achievable by previous methods. While many past hypothesized evolutionary relationships for Fundulidae were reinforced, several alternative relationships are hypothesized at basal nodes resulting in a re-analysis of the deeper structure of family Fundulidae. These results reveal family Fundulidae as a paraphyletic grouping of members of genus Fundulus and Lucania and supports the previous synonymy of genus Adinia with genus Fundulus.

Contrasting phylogeographic histories between broadly sympatric topminnows in the Fundulus notatus species complex.

Sympatrically distributed closely related species provide opportunities for studying evolutionary patterns of diversification. Such studies must account for historical contingencies in interpreting contemporary patterns of variation. Topminnows in the Fundulus notatus species complex are distributed sympatrically across much of the southern and Midwestern United States. Throughout most of their ranges F. olivaceus is often found in headwater stream habitats, and F. notatus is more typically distributed along the margins of larger river habitats. However, in some drainages, ecological associations of the respective species are reversed, with F. notatus populations isolated in headwater streams and F. olivaceus in downstream river habitats. Phylogeographic analyses of AFLP marker and multi-locus sequence data detected historical isolation in F. notatus consistent with pre-Pleistocene drainage patterns. Four F. notatus clades corresponded to (i) the Western Gulf Slope, (ii) the southwestern Ouachita Highlands, (iii) the Mobile Basin, and (iv) central Coastal Plain and Mississippi River Basin. In contrast, a relative lack of range-wide geographic structure in F. olivaceus is consistent with recent range expansion over much of the same geographic area. The southwestern Ouachita Highlands and Mobile Basin F. notatus clades corresponded to regions where ecological associations between the two species are reversed, providing evidence of the independent evolution of variation in contemporary habitat associations. Fundulus olivaceus from several drainages demonstrated introgression of mitochondrial DNA from F. notatus, but none of the sites in this study included individuals with hybrid ancestry in their nuclear genome. Phylogenetic analyses that included only nuclear loci supported the reciprocal monophyly of F. notatus, F. olivaceus and a third narrowly endemic species, Fundulus euryzonus, and supported a sister relationship between F. olivaceus and F. euryzonus.

Evolution of a sexually dimorphic trait in a broadly distributed topminnow (Fundulus olivaceus).

Understanding the interaction between sexual and natural selection within variable environments is crucial to our understanding of evolutionary processes. The handicap principle predicts females will prefer males with exaggerated traits provided those traits are indicators of male quality to ensure direct or indirect female benefits. Spatial variability in ecological factors is expected to alter the balance between sexual and natural selection that defines the evolution of such traits. Male and female blackspotted topminnows (Fundulidae: Fundulus olivaceus) display prominent black dorsolateral spots that are variable in number across its broad range. We investigated variability in spot phenotypes at 117 sites across 13 river systems and asked if the trait was sexually dimorphic and positively correlated with measures of fitness (condition and gonadosomatic index [GSI]). Laboratory and mesocosm experiments assessed female mate choice and predation pressure on spot phenotypes. Environmental and community data collected at sampling locations were used to assess predictive models of spot density at the individual, site, and river system level. Greater number of spots was positively correlated with measures of fitness in males. Males with more spots were preferred by females and suffered greater mortality due to predation. Water clarity (turbidity) was the best predictor of spot density on the drainage scale, indicating that sexual and natural selection for the trait may be mediated by local light environments.

Relative influences of historical and contemporary forces shaping the distribution of genetic variation in the Atlantic killifish, Fundulus heteroclitus.

A major goal of population genetics research is to identify the relative influences of historical and contemporary processes that serve to structure genetic variation. Most population genetic models assume that populations exist in a state of migration-drift equilibrium. However, in the past this assumption has rarely been verified, and is likely rarely achieved in natural populations. We assessed the equilibrium status at both local and regional scales of the Atlantic killifish, Fundulus heteroclitus. This species is a model organism for the study of adaptive clinal variation, but has also experienced a complicated history of range expansion and secondary contact following allopatric divergence, potentially obscuring the influence of contemporary evolutionary processes. Presumptively neutral genetic markers (microsatellites) demonstrated zones of secondary intergradation among coastal populations centred around northern New Jersey and the Chesapeake Bay region. Analysis of genetic variation indicated isolation by distance among some populations and provided supporting evidence that the Delaware Bay, but not the Chesapeake Bay, has acted as a barrier to dispersal among coastal populations. Bayesian estimates indicated large effective population sizes and low migration rates, and were in good agreement with empirically derived estimates of population and neighbourhood size from mark-recapture studies. These data indicate that populations are not in migration-drift equilibrium at a regional scale, and suggest that contributing factors include large population size combined with relatively low migration rates. These conditions should be considered when interpreting the evolutionary significance of the distribution of genetic variation among F. heteroclitus populations.

Acquisition of endonuclease specificity during evolution of L1 retrotransposon.

L1 is the most proliferative autonomous retroelement that comprises about 20% of mammalian genomes. Why L1s have proliferated so extensively in mammalian genomes is an important yet unsolved question. L1 copies are amplified via retrotransposition, in which the DNA cleavage specificity by the L1-encoded endonuclease (EN) primarily dictates sites of insertion. Whereas mammalian L1s show target preference for 5'-TTAAAA-3', other L1-like elements exhibit various degrees of target specificity. To gain insights on diversification of the EN specificity during L1 evolution, ENs of zebrafish L1 elements were analyzed here. We revealed that they form 3 discrete clades, M, F, and Tx1, which is in stark contrast to a single L1 clade in mammalian species. Interestingly, zebrafish clade M elements cluster as a sister group of mammalian L1s and show target-site preference for 5'-TTAAAA-3'. In contrast, elements of the clade F, the immediate outgroup of the clade M, show little specificity. We identified certain clade-specific amino acid residues in EN, many of which are located in the cleft that recognizes the substrate, suggesting that these amino acid alterations have generated 2 types of ENs with different substrate specificities. The distribution pattern of the 3 clades suggests a possibility that the acquisition of target specificity by the L1 ENs improved the L1 fitness under the circumstances in mammalian hosts.

Microsatellite analysis of the phylogeography, Pleistocene history and secondary contact hypotheses for the killifish, Fundulus heteroclitus.

The mummichog, Fundulus heteroclitus, exhibits extensive latitudinal clinal variation in a number of physiological and biochemical traits, coupled with phylogeographical patterns at mitochondrial and nuclear DNA loci that suggest a complicated history of spatially variable selection and secondary intergradation. This species continues to serve as a model for understanding local and regional adaptation to variable environments. Resolving the influences of historical processes on the distribution of genetic variation within and among extant populations of F. heteroclitus is crucial to a better understanding of how populations evolve in the context of contemporary environments. In this study, we analysed geographical patterns of genetic variation at eight microsatellite loci among 15 populations of F. heteroclitus distributed throughout the North American range of the species from Nova Scotia to Georgia. Genetic variation in Northern populations was lower than in Southern populations and was strongly correlated with latitude throughout the species range. The most common Northern alleles at all eight loci exhibited concordant latitudinal clinal patterns, and the existence of an abrupt transition zone in allele frequencies between Northern and Southern populations was similar to that observed for mitochondrial DNA and allozyme loci. A significant pattern of isolation by distance was observed both within and between northern and southern regions. This pattern was unexpected, particularly for northern populations, given the recent colonization history of post-Pleistocene habitats, and was inconsistent with either a recent northward population expansion or a geographically restricted northern Pleistocene refugium. The data provided no evidence for recent population bottlenecks, and estimates of historical effective population sizes suggest that post-Pleistocene populations have been large throughout the species distribution. These results suggest that F. heteroclitus was broadly distributed throughout most of its current range during the last glacial event and that the abrupt transition in allele frequencies that separate Northern and Southern populations may reflect regional disequilibrium conditions associated with the post-Pleistocene colonization history of habitats in that region.

Teleost fish genomes contain a diverse array of L1 retrotransposon lineages that exhibit a low copy number and high rate of turnover.

Retrotransposable elements exhibit a wide range of variation in population dynamics, abundance, and lineage diversity among host genomes across taxa. This range of diversity is illustrated by a single well-defined constituent monophyletic clade of L1 non-LTR retrotransposons that is shared between mammalian and teleost fish genomes. Despite the clear phylogenetic relationships that exist between mammalian and teleost L1 sequences, these elements exhibit markedly different dynamics within their respective taxa. While mammalian genomes typically contain a single, abundant lineage of L1 elements that traces millions of years of evolution, the zebraflsh genome was recently shown to exhibit a high diversity of ancient lineages coexisting at a very low copy number and apparently exhibiting a high rate of turnover. In the present study, a combination of degenerate PCR, lineage-specific PCR, and genomic Southern blot analysis is utilized to demonstrate high L1 lineage diversity, low copy number, and a high proportion of polymorphic inserts in the genomes of the killifish species, Fundulus heteroclitus. Additional species surveyed by degenerate PCR include Cyprinodon variegatus, Rivulus marmoratus, and Menidia beryllina. These results further support the generality of the differences that exist in host-element dynamics between teleost fish and mammalian genomes with regard to L1 retrotransposons.

Single-locus latitudinal clines and their relationship to temperate adaptation in metabolic genes and derived alleles in Drosophila melanogaster.

We report a study in Drosophila melanogaster of latitudinal clines for 23 SNPs embedded in 13 genes (Pgi, Gapdh1, UGPase, Pglym78, Pglym87, Eno, Men, Gdh, Sod, Pgk, Mdh1, TreS, Treh) representing various metabolic enzymes. Our samples are from 10 populations spanning latitude from southern Florida to northern Vermont. Three new clines with latitude were detected. These are the amino acid polymorphisms in the NAD-dependent glutamate dehydrogenase (Gdh) and trehalase (Treh) genes, and a silent site polymorphism in the UDP-glucose pyrophosphorylase gene (UGPase). The result, when combined with the overall incidence and pattern of reports for six other genes (Adh, Gpdh, Pgm, G6pd, 6Pgd, Hex-C), presents a picture of latitudinal clines in metabolic genes prevalent around the branch point of competing pathways. For six of the seven amino acid polymorphisms showing significant latitudinal clines in North America, the derived allele is the one increasing with latitude, suggesting temperate adaptation. This is consistent with a model of an Afrotropical ancestral species adapting to temperate climates through selection favoring new mutations.

L1 (LINE-1) retrotransposon diversity differs dramatically between mammals and fish.

L1 retrotransposons replicate (amplify) by copying (reverse transcribing) their RNA transcript into genomic DNA. The evolutionary history of L1 in mammals has been unique. In mice and humans approximately 80 million years of L1 evolution and replication produced a single evolutionary lineage of L1 elements while generating approximately 20% of the genomic mass in each species. By contrast, zebrafish contain >30 distinct L1 lineages that have generated approximately one-tenth as much DNA. We contend that, by becoming far more permissive of interspersed repeated DNA than other organisms, mammals are conducive to competition between L1 families for replicative dominance, and that this competition, perhaps for the host factors required for L1 replication, results in a single L1 lineage.

Clines and adaptive evolution in the methuselah gene region in Drosophila melanogaster.

In an effort to characterize further the patterns of selection and adaptive evolution at the methuselah locus in Drosophila species, we extended an analysis of geographical variation to include single nucleotide polymorphisms (SNPs) in adjacent genes on either side of the mth locus, and examined the molecular variation in a neighbouring methuselah paralogue (mth2). An analysis of 13 SNPs spanning a region of nearly 19 kilobases surrounding the mth locus demonstrated that a clinal pattern associated with the most common mth haplotype does not extend to adjacent gene loci, providing compelling evidence that the clinal pattern results from selection on as yet unidentified sites associated with the functional mth locus. mth2 exhibited a significant pattern of adaptive divergence among D. melanogaster, D. simulans and D. yakuba similar to that seen at mth. However, Ka : Ks ratios indicate a difference in levels of functional constraint at the two methuselah, loci with mth2 exhibiting a five- to six-fold reduction in levels of amino acid divergence relative to mth.