Phylogenomics reveals extensive introgression and a case of mito-nuclear discordance in the killifish genus Kryptolebias.

Introgression is a widespread evolutionary process leading to phylogenetic inconsistencies among distinct parts of the genomes, particularly between mitochondrial and nuclear-based phylogenetic reconstructions (e.g., mito-nuclear discordances). Here, we used mtDNA and genome-wide nuclear sites to provide the first phylogenomic-based hypothesis on the evolutionary relationships within the killifish genus Kryptolebias. In addition, we tested for evidence of past introgression in the genus given the multiple reports of undergoing hybridization between its members. Our mtDNA phylogeny generally agreed with the relationships previously proposed for the genus. However, our reconstruction based on nuclear DNA revealed an unknown lineage - Kryptolebias sp. 'ESP' - as the sister group of the self-fertilizing mangrove killifishes, K. marmoratus and K. hermaphroditus. All individuals sequenced of Kryptolebias sp. 'ESP' had the same mtDNA haplotype commonly observed in K. hermaphroditus, demonstrating a clear case of mito-nuclear discordance. Our analysis further confirmed extensive history of introgression between Kryptolebias sp. 'ESP' and K. hermaphroditus. Population genomics analyses indicate no current gene flow between the two lineages, despite their current sympatry and history of introgression. We also confirmed introgression between other species pairs in the genus that have been recently reported to form hybrid zones. Overall, our study provides a phylogenomic reconstruction covering most of the Kryptolebias species, reveals a new lineage hidden in a case of mito-nuclear discordance, and provides evidence of multiple events of ancestral introgression in the genus. These findings underscore the importance of investigating different genomic information in a phylogenetic framework, particularly in taxa where introgression is common as in the sexually diverse mangrove killifishes.

Extensive hybridization and past introgression between divergent lineages in a quasi-clonal hermaphroditic fish: Ramifications for species concepts and taxonomy.

Extreme inbreeding is expected to reduce the incidence of hybridization, serving as a prezygotic barrier. Mangrove rivulus is a small killifish that reproduces predominantly by self-fertilization, producing highly homozygous lines throughout its geographic range. The Bahamas and Caribbean are inhabited by two highly diverged phylogeographic lineages of mangrove rivulus, Kryptolebias marmoratus and a 'Central clade' closely related to K. hermaphroditus from Brazil. The two lineages are largely allopatric, but recently were found in syntopy on San Salvador, Bahamas, where a single hybrid was reported. To better characterize the degree of hybridization and the possibility of secondary introgression, here we conducted a detailed genetic analysis of the contact zone on San Salvador. Two mixed populations were identified, one of which contained sexually mature hybrids. The distribution of heterozygosity at diagnostic microsatellite loci in hybrids showed that one of these hybrids was an immediate offspring from the K. marmoratus x Central clade cross, whereas the remaining five hybrids were products of reproduction by self-fertilization for 1-3 generations following the initial cross. Two hybrids had mitochondrial haplotypes of K. marmoratus and the remaining four hybrids had a haplotype of the Central clade, indicating that crosses go in both directions. In hybrids, alleles of parental lineages were represented in equal proportions suggesting lack of recent backcrossing to either of the parental lineages. However, sympatric populations of two lineages were less diverged than allopatric populations, consistent with introgression. Results are discussed in terms of applicability of the biological species concept for isogenic, effectively clonal, organisms.

More than meets the eye: syntopic and morphologically similar mangrove killifish species show different mating systems and patterns of genetic structure along the Brazilian coast.

Different mating systems can strongly affect the extent of genetic diversity and population structure among species. Given the increased effects of genetic drift on reduced population size, theory predicts that species undergoing self-fertilisation should have greater population structure than outcrossed species; however, demographic dynamics may affect this scenario. The mangrove killifish clade is composed of the two only known examples of self-fertilising species among vertebrates (Kryptolebias marmoratus and Kryptolebias hermaphroditus). A third species in this clade, Kryptolebias ocellatus, inhabits mangrove forests in southeast Brazil; however, its mating system and patterns of genetic structure have been rarely explored. Here, we examined the genetic structure and phylogeographic patterns of K. ocellatus along its distribution, using mitochondrial DNA and microsatellites to compare its patterns of genetic structure with the predominantly selfing and often-syntopic, K. hermaphroditus. Our results indicate that K. ocellatus reproduces mainly by outcrossing, with no current evidence of selfing, despite being an androdioecious species. Our results also reveal a stronger population subdivision in K. ocellatus compared to K. hermaphroditus, contrary to the theoretical predictions based on reproductive biology of the two species. Our findings indicate that, although morphologically similar, K. ocellatus and K. hermaphroditus had remarkably different evolutionary histories when colonising the same mangrove areas in southeastern Brazil, with other factors (e.g., time of colonisation, dispersal/establishment capacity) having more profound effects on the current population structuring of those species than differences in mating systems.

Whole-genome sequencing reveals the extent of heterozygosity in a preferentially self-fertilizing hermaphroditic vertebrate.

The mangrove rivulus, Kryptolebias marmoratus, is one of only two self-fertilizing hermaphroditic fish species and inhabits mangrove forests. While selfing can be advantageous, it reduces heterozygosity and decreases genetic diversity. Studies using microsatellites found that there are variable levels of selfing among populations of K. marmoratus, but overall, there is a low rate of outcrossing and, therefore, low heterozygosity. In this study, we used whole-genome data to assess the levels of heterozygosity in different lineages of the mangrove rivulus and infer the phylogenetic relationships among those lineages. We sequenced whole genomes from 15 lineages that were completely homozygous at microsatellite loci and used single nucleotide polymorphisms (SNPs) to determine heterozygosity levels. More variation was uncovered than in studies using microsatellite data because of the resolution of full genome sequencing data. Moreover, missense polymorphisms were found most often in genes associated with immune function and reproduction. Inferred phylogenetic relationships suggest that lineages largely group by their geographic distribution. The use of whole-genome data provided further insight into genetic diversity in this unique species. Although this study was limited by the number of lineages that were available, these data suggest that there is previously undescribed variation within lineages of K. marmoratus that could have functional consequences and (or) inform us about the limits to selfing (e.g., genetic load, accumulation of deleterious mutations) and selection that might favor the maintenance of heterozygosity. These results highlight the need to sequence additional individuals within and among lineages.

Natural hybridization between divergent lineages in a selfing hermaphroditic fish.

By definition, mating between individuals is infrequent in highly selfing organisms, and so too, therefore, hybridization should be rare between genetically divergent lineages in predominantly self-fertilizing species. Notwithstanding these expectations, here we report a remarkable case of natural hybridization between highly diverged phylogeographic lineages of the mangrove rivulus, a small killifish that reproduces predominantly by self-fertilization and typically is found as highly homozygous lines in most parts of its extensive geographical range. Two distinctive genetic lineages (Kryptolebias marmoratus and a 'Central clade' closely related to K. hermaphroditus) previously were not known in sympatry, but were found by us to co-occur on San Salvador, Bahamas. Genetic analyses of a mitochondrial and multiple nuclear markers determined the direction of a cross producing a hybrid fish. Furthermore, we show that this hybrid individual was viable, as it successfully reproduced by self-fertilization for two generations. Additional sampling of this population will be necessary to determine if backcrossing of hybrids to the parental lineages occurs in nature and to analyse whether such backcross progeny are viable. Application of the biological species concept (BSC) is traditionally difficult in clonally reproducing organisms. Our results show that although mangrove rivulus fish are mostly highly selfing in nature (resulting in isogenic, effectively clonal and homozygous progeny), classification within this taxonomic complex need not be incompatible with the BSC.

Transcriptomics of diapause in an isogenic self-fertilizing vertebrate.

BACKGROUND: Many vertebrate species have the ability to undergo weeks or even months of diapause (a temporary arrest of development during early ontogeny). Identification of diapause genes has been challenging due in part to the genetic heterogeneity of most vertebrate animals. RESULTS: Here we take the advantage of the mangrove rivulus fish (Kryptolebias marmoratus or Kmar)-the only vertebrate that is extremely inbred due to consistent self-fertilization-to generate isogenic lineages for transcriptomic dissection. Because the Kmar genome is not publicly available, we built de novo genomic (642 Mb) and transcriptomic assemblies to serve as references for global genetic profiling of diapause in Kmar, via RNA-Seq. Transcripts unique to diapause in Kmar proved to constitute only a miniscule fraction (0.1 %) of the total pool of transcribed products. Most genes displayed lower expression in diapause than in post-diapause. However, some genes (notably dusp27, klhl38 and sqstm1) were significantly up-regulated during diapause, whereas others (col9a1, dspp and fmnl1) were substantially down-regulated, compared to both pre-diapause and post-diapause. CONCLUSION: Kmar offers a strong model for understanding patterns of gene expression during diapause. Our study highlights the importance of using a combination of genome and transcriptome assemblies as references for NGS-based RNA-Seq analyses. As for all identified diapause genes, in future studies it will be critical to link various levels of RNA expression with the functional roles of the coded products.

A Genetic Test for Whether Pairs of Hermaphrodites Can Cross-Fertilize in a Selfing Killifish.

Kryptolebias marmoratus, a small killifish that lives in mangrove habitat from southern Florida to Brazil, is one of the planet's only known self-fertilizing hermaphroditic vertebrates. Generation after generation, hermaphroditic individuals simultaneously produce sperm and eggs and internally self-fertilize to produce what are, in effect, highly inbred clones of themselves. Although populations are composed primarily of hermaphrodites, they also contain some true males. The frequency of males in a population varies geographically, from <2% in Florida to as high as 25% in Belize. Males are known to mate occasionally with hermaphrodites, thereby releasing genetic variation that has profound consequences for population genetic structure. However, it is unknown whether hermaphrodites can or do sporadically mate with each other also. Here, we test whether hermaphroditic individuals of the killifish Kryptolebias marmoratus are capable of crossing with one another, in addition to their much more common habits of self-fertilization and occasional outcrossing with pure males. We employ an experimental design in which replicate hermaphrodite pairs were housed together and allowed to reproduce naturally. Among 173 embryos screened at diagnostic microsatellite loci, all were found to result from selfing (i.e., no embryos were the product of a hermaphrodite cross). We thus conclude that hermaphrodite pairs are unlikely to cross, or do so exceedingly rarely.

Genetic Subdivision and Variation in Selfing Rates Among Central American Populations of the Mangrove Rivulus, Kryptolebias marmoratus.

We used 32 polymorphic microsatellite loci to investigate how a mixed-mating system affects population genetic structure in Central American populations (N = 243 individuals) of the killifish Kryptolebias marmoratus (mangrove rivulus), 1 of 2 of the world's only known self-fertilizing vertebrates. Results were also compared with previous microsatellite surveys of Floridian populations of this species. For several populations in Belize and Honduras, population structure and genetic differentiation were pronounced and higher than in Florida, even though the opposite trend was expected because populations in the latter region were presumably smaller and highly selfing. The deduced frequency of selfing (s) ranged from s = 0.39-0.99 across geographic locales in Central America. This heterogeneity in selfing rates was in stark contrast to Florida, where s > 0.9. The frequency of outcrossing in a population (t = 1 - s) was tenuously correlated with local frequencies of males, suggesting that males are one of many factors influencing outcrossing. Observed distributions of individual heterozygosity showed good agreement with expected distributions under an equilibrium mixed-mating model, indicating that rates of selfing remained relatively constant over many generations. Overall, our results demonstrate the profound consequences of a mixed-mating system for the genetic architecture of a hermaphroditic vertebrate.

Allard's argument versus Baker's contention for the adaptive significance of selfing in a hermaphroditic fish.

Fertilization assurance (Baker's contention) and multilocus coadaptation (Allard's argument) are two distinct hypotheses for the adaptive significance of self-fertilization in hermaphroditic taxa, and both scenarios have been invoked to rationalize isogenicity via incest in various plants and invertebrate animals with predominant selfing. Here we contrast Allard's argument and Baker's contention as applied to the world's only known vertebrate that routinely self-fertilizes. We pay special attention to frequencies of locally most common multilocus genotypes in Floridian populations of the Mangrove Rivulus (Kryptolebias marmoratus). Isogenicity patterns in this fish appear inconsistent with Allard's argument, thus leaving Baker's contention as the more plausible scenario (a result also supported by natural history information for this species). These results contrast with the isogenicity patterns and conclusions previously drawn from several self-fertilizing plants and invertebrate animal species. Thus, the adaptive significance of selfing apparently varies across hermaphroditic taxa.

Multiple mating and clutch size in invertebrate brooders versus pregnant vertebrates.

We summarize the genetic literature on polygamy rates and sire numbers per clutch in invertebrate animals that brood their offspring and then compare findings with analogous data previously compiled for vertebrate species displaying viviparity or other pregnancy-like syndromes. As deduced from molecular parentage analyses of several thousand broods from more than 100 "pregnant" species, invertebrate brooders had significantly higher mean incidences of multiple mating than pregnant vertebrates, a finding generally consistent with the postulate that clutch size constrains successful mate numbers in species with extended parental care. However, we uncovered no significant correlation in invertebrates between brood size and genetically deduced rates of multiple mating by the incubating sex. Instead, in embryo-gestating animals otherwise as different as mammals and mollusks, polygamy rates and histograms of successful mates per brooder proved to be strikingly similar. Most previous studies have sought to understand why gestating parents have so many mates and such high incidences of successful multiple mating; an alternative perspective based on logistical constraints turns the issue on its head by asking why mate numbers and polygamy rates are much lower than they theoretically could be, given the parentage-resolving power of molecular markers and the huge sizes of many invertebrate broods.

Molecular evidence for multiple paternity in a population of the Viviparous Tule Perch Hysterocarpus traski.

Population density might be an important variable in determining the degree of multiple paternity. In a previous study, a high level of multiple paternity was detected in the shiner perch Cymatogaster aggregata, a species with high population density and a high mate encounter rate. The tule perch Hysterocarpus traski is phylogenetically closely related to C. aggregata, but it has relatively lower population density, which may result in distinct patterns of multiple paternity in these 2 species. To test the hypothesis that mate encounter rate may affect the rate of successful mating, we used polymorphic microsatellite markers to identify multiple paternity in the progeny arrays of 12 pregnant females from a natural population of tule perch. Multiple paternity was detected in 11 (92%) of the 12 broods. The number of sires per brood ranged from 1 to 4 (mean 2.5) but with no correlation between sire number and brood size. Although the brood size of tule perch is considerably larger than that of shiner perch (40.7 vs. 12.9, respectively), the average number of sires per brood in tule perch is much lower than that in shiner perch (2.5 vs. 4.6, respectively). These results are consistent with the hypothesis that mate encounter rate is an important factor affecting multiple mating.

Out to sea: ocean currents and patterns of asymmetric gene flow in an intertidal fish species.

Passive dispersal via wind or ocean currents can drive asymmetric gene flow, which influences patterns of genetic variation and the capacity of populations to evolve in response to environmental change. The mangrove rivulus fish (Kryptolebias marmoratus), hereafter "rivulus," is an intertidal fish species restricted to the highly fragmented New World mangrove forests of Central America, the Caribbean, the Bahamas, and Florida. Mangrove patches are biological islands with dramatic differences in both abiotic and biotic conditions compared to adjacent habitat. Over 1,000 individual rivulus across 17 populations throughout its range were genotyped at 32 highly polymorphic microsatellites. Range-wide population genetic structure was evaluated with five complementary approaches that found eight distinct population clusters. However, an analysis of molecular variance indicated significant population genetic structure among regions, populations within regions, sampling locations within populations, and individuals within sampling locations, indicating that rivulus has both broad- and fine-scale genetic differentiation. Integrating range-wide genetic data with biophysical modeling based on 10 years of ocean current data showed that ocean currents and the distance between populations over water drive gene flow patterns on broad scales. Directional migration estimates suggested some significant asymmetries in gene flow that also were mediated by ocean currents and distance. Specifically, populations in the center of the range (Florida Keys) were identified as sinks that received migrants (and alleles) from other populations but failed to export individuals. These populations thus harbor genetic variation, perhaps even from extirpated populations across the range, but ocean currents and complex arrangements of landmasses might prevent the distribution of that genetic variation elsewhere. Hence, the inherent asymmetry of ocean currents shown to impact both genetic differentiation and directional migration rates may be responsible for the complex distribution of genetic variation across the range and observed patterns of metapopulation structure.

Detection of Horizontal Gene Transfer Mediated by Natural Conjugative Plasmids in E. coli.

Conjugation represents one of the main mechanisms facilitating horizontal gene transfer in Gram-negative bacteria. This work describes methods for the study of the mobilization of naturally occurring conjugative plasmids, using two naturally-occurring plasmids as an example. These protocols rely on the differential presence of selectable markers in donor, recipient, and conjugative plasmid. Specifically, the methods described include 1) the identification of natural conjugative plasmids, 2) the quantification of conjugation rates in solid culture, and 3) the diagnostic detection of the antibiotic resistance genes and plasmid replicon types in transconjugant recipients by polymerase chain reaction (PCR). The protocols described here have been developed in the context of studying the evolutionary ecology of horizontal gene transfer, to screen for the presence of conjugative plasmids carrying antibiotic-resistance genes in bacteria found in the environment. The efficient transfer of conjugative plasmids observed in these experiments in culture highlights the biological relevance of conjugation as a mechanism promoting horizontal gene transfer in general and the spread of antibiotic resistance in particular.

Long-term retention of self-fertilization in a fish clade.

Among vertebrate animals, only the mangrove rivulus (Kryptolebias marmoratus) was known to self-fertilize. Here, we use microsatellite analyses to document a high selfing rate (97%) in a related nominal species, Kryptolebias ocellatus, which likewise is androdioecious (populations consist of males and hermaphrodites). In contrast, we find no evidence of self-fertilization in Kryptolebias caudomarginatus (an androdioecious species closely related to the marmoratus-ocellatus clade) or in Kryptolebias brasiliensis (a dioecious outgroup). These findings indicate that the initiation of self-fertilization predated the origin of the marmoratus-ocellatus clade. From mitochondrial DNA sequences and microsatellite data, we document a substantial genetic distance between Kryptolebias marmoratus and K. ocellatus, implying that the selfing capacity has persisted in these fishes for at least several hundred thousand years.

Multiple Copies of a qnrB19 Gene Are Carried by Tandem Repeats of an IS26 Composite Transposon in an Escherichia coli Plasmid.

The fully assembled genome of Escherichia coli strain BR1220 shows a triple tandemly arrayed IS26 composite transposon carrying a qnrB19 quinolone resistance gene in a 100-kb multidrug resistance plasmid (1.6 copies per chromosome [CPC]) and a 2.6-kb Col(pHAD28) plasmid (27.8 CPC) with a nearly identical qnrB19 gene region.

Effects of Pleistocene climatic fluctuations on the phylogeographic and demographic histories of Pacific herring (Clupea pallasii).

We gathered mitochondrial DNA sequences (557 bp from the control region in 935 specimens and 668 bp of the cytochrome b gene in 139 specimens) of Pacific herring collected from 20 nearshore localities spanning the species' extensive range along the North Pacific coastlines of Asia and North America. Haplotype diversity and nucleotide diversity were high, and three major phylogeographic lineages (sequence divergences ca. 1.5%) were detected. Using a variety of phylogenetic methods, coalescent reasoning, and molecular dating interpreted in conjunction with paleoclimatic and physiographic evidence, we infer that the genetic make-up of extant populations of C. pallasii was shaped by Pleistocene environmental impacts on the historical demography of this species. A deep genealogical split that cleanly distinguishes populations in the western vs. eastern North Pacific probably originated as a vicariant separation associated with a glacial cycle that drove the species southward and isolated two ancestral populations in Asia and North America. Another deep genealogical split may have involved either a vicariant isolation of a third herring lineage (perhaps originally in the Gulf of California) or it may have resulted simply from the long coalescent times that are possible in large populations. Coalescent analyses showed that all the three evolutionary lineages of C. pallasii experienced major expansions in their most recent histories after having remained more stable in the preceding periods. Independent of the molecular calibration chosen, populations of C. pallasii appear to have remained stable or grown throughout the periods that covered at least two major glaciations, and probably more.

Against the Odds: Hybrid Zones between Mangrove Killifish Species with Different Mating Systems.

Different mating systems are expected to affect the extent and direction of hybridization. Due to the different levels of sexual conflict, the weak inbreeder/strong outbreeder (WISO) hypothesis predicts that gametes from self-incompatible (SI) species should outcompete gametes from self-compatible (SC) ones. However, other factors such as timing of selfing and unilateral incompatibilities may also play a role on the direction of hybridization. In addition, differential mating opportunities provided by different mating systems are also expected to affect the direction of introgression in hybrid zones involving outcrossers and selfers. Here, we explored these hypotheses with a unique case of recent hybridization between two mangrove killifish species with different mating systems, Kryptolebias ocellatus (obligately outcrossing) and K. hermaphroditus (predominantly self-fertilizing) in two hybrid zones in southeast Brazil. Hybridization rates were relatively high (~20%), representing the first example of natural hybridization between species with different mating systems in vertebrates. All F1 individuals were sired by the selfing species. Backcrossing was small, but mostly asymmetrical with the SI parental species, suggesting pattern commonly observed in plant hybrid zones with different mating systems. Our findings shed light on how contrasting mating systems may affect the direction and extent of gene flow between sympatric species, ultimately affecting the evolution and maintenance of hybrid zones.

Pronounced reproductive skew in a natural population of green swordtails, Xiphophorus helleri.

For many species in nature, a sire's progeny may be distributed among a few or many dams. This poses logistical challenges--typically much greater across males than across females--for assessing means and variances in mating success (number of mates) and reproductive success (number of progeny). Here we overcome these difficulties by exhaustively analyzing a population of green swordtail fish (Xiphophorus helleri) for genetic paternity (and maternity) using a suite of highly polymorphic microsatellite loci. Genetic analyses of 1476 progeny from 69 pregnant females and 158 candidate sires revealed pronounced skews in male reproductive success both within and among broods. These skews were statistically significant, greater than in females, and correlated in males but not in females with mating success. We also compare the standardized variances in swordtail reproductive success to the few such available estimates for other taxa, notably several mammal species with varied mating systems and degrees of sexual dimorphism. The comparison showed that the opportunity for selection on male X. helleri is among the highest yet reported in fishes, and it is intermediate compared to estimates available for mammals. This study is one of a few exhaustive genetic assessments of joint-sex parentage in a natural fish population, and results are relevant to the operation of sexual selection in this sexually dimorphic, high-fecundity species.

Rapid concerted evolution in animal mitochondrial DNA.

Recombinational genetic processes are thought to be rare in the uniparentally inherited mitochondrial (mt) DNA molecules of vertebrates and other animals. Here, however, we document extremely rapid concerted microevolution, probably mediated by frequent gene conversion events, of duplicated sequences in the mtDNA control region of mangrove killifishes (Kryptolebias marmoratus). In local populations, genetic distances between paralogous loci within an individual were typically smaller (and often zero) than those between orthologous loci in different specimens. These findings call for the recognition of concerted evolution as a microevolutionary process and gene conversion as a likely recombinational force in animal mtDNA. The previously unsuspected power of these molecular phenomena could greatly impact mtDNA dynamics within germ cell lineages and in local animal populations.

Nucleotide variation at the dopa decarboxylase (Ddc) gene in natural populations of Drosophila melanogaster.

We studied nucleotide sequence variation at the gene coding for dopa decarboxylase (Ddc) in seven populations of Drosophila melanogaster. Strength and pattern of linkage disequilibrium are somewhat distinct in the extensively sampled Spanish and Raleigh populations. In the Spanish population, a few sites are in strong positive association, whereas a large number of sites in the Raleigh population are associated nonrandomly but the association is not strong. Linkage disequilibrium analysis shows presence of two groups of haplotypes in the populations, each of which is fairly diverged, suggesting epistasis or inversion polymorphism. There is evidence of two forms of natural selection acting on Ddc. The McDonald-Kreitman test indicates a deficit of fixed amino acid differences between D. melanogaster and D. simulans, which may be due to negative selection. An excess of derived alleles at high frequency, significant according to the H-test, is consistent with the effect of hitchhiking. The hitchhiking may have been caused by directional selection downstream of the locus studied, as suggested by a gradual decrease of the polymorphism-to-divergence ratio. Altogether, the Ddc locus exhibits a complicated pattern of variation apparently due to several evolutionary forces. Such a complex pattern may be a result of an unusually high density of functionally important genes.