Kathryn R. Elmer

Rapid evolution and selection inferred from the transcriptomes of sympatric crater lake cichlid fishes

Crater lakes provide a natural laboratory to study speciation of cichlid fishes by ecological divergence. Up to now, there has been a dearth of transcriptomic and genomic information that would aid in understanding the molecular basis of the phenotypic differentiation between young species. We used next‐generation sequencing (Roche 454 massively parallel pyrosequencing) to characterize the diversity of expressed sequence tags between ecologically divergent, endemic and sympatric species of cichlid fishes from crater lake Apoyo, Nicaragua: benthic Amphilophus astorquii and limnetic Amphilophus zaliosus. We obtained 24 174 A. astorquii and 21 382 A. zaliosus high‐quality expressed sequence tag contigs, of which 13 106 pairs are orthologous between species. Based on the ratio of nonsynonymous to synonymous substitutions, we identified six sequences exhibiting signals of strong diversifying selection (Ka/Ks > 1). These included genes involved in biosynthesis, metabolic processes and development. This transcriptome sequence variation may be reflective of natural selection acting on the genomes of these young, sympatric sister species. Based on Ks ratios and p‐distances between 3′‐untranslated regions (UTRs) calibrated to previously published species divergence times, we estimated a neutral transcriptome‐wide substitutional mutation rate of ∼1.25 × 10−6 per site per year. We conclude that next‐generation sequencing technologies allow us to infer natural selection acting to diversify the genomes of young species, such as crater lake cichlids, with much greater scope than previously possible.

Local variation and parallel evolution: morphological and genetic diversity across a species complex of neotropical crater lake cichlid fishes

The polychromatic and trophically polymorphic Midas cichlid fish species complex (Amphilophus cf. citrinellus) is an excellent model system for studying the mechanisms of speciation and patterns of phenotypic diversification in allopatry and in sympatry. Here, we first review research to date on the species complex and the geological history of its habitat. We analyse body shape variation from all currently described species in the complex, sampled from six crater lakes (maximally 1.2–23.9 kyr old) and both great lakes in Nicaragua. We find that Midas cichlid populations in each lake have their own characteristic body shape. In lakes with multiple sympatric species of Midas cichlid, each species has a distinct body shape. Across the species complex, most body shape change relates to body depth, head, snout and mouth shape and caudal peduncle length. There is independent parallel evolution of an elongate limnetic species in at least two crater lakes. Mitochondrial genetic diversity is higher in crater lakes with multiple species. Midas cichlid species richness increases with the size and age of the crater lakes, though no such relationship exists for the other syntopic fishes. We suggest that crater lake Midas cichlids follow the predicted pattern of an adaptive radiation, with early divergence of each crater lake colonization, followed by intralacustrine diversification and speciation by ecological adaptation and sexual selection.

Color Assortative Mating Contributes to Sympatric Divergence of Neotropical Cichlid Fish

It is still debated vigorously whether sexual selection can result in speciation without physical barriers to gene flow. In this study, we used field data and molecular methods to investigate the gold-normal color polymorphism in two endemic cichlid fish species of crater lake Xiloá, Nicaragua. We found significant assortative mating by color in both Amphilophus xiloaensis and A. sagittae. Focusing on A. xiloaensis, microsatellite allele frequencies, an assignment test, and model-based cluster analysis demonstrates significant and clear genetic differentiation (F ST = 0.03) between gold and normal individuals in sympatry. In addition, we find genetic differentiation between all three sympatric and ecologically distinct Midas cichlid species of Lake Xiloá, A. amarillo, A. sagittae, and A. xiloaensis (F ST = 0.03 – 0.19), and clear genetic isolation of these species from their closest relative (A. citrinellus) in the neighboring great lake Managua. The A. xiloaensis gold morph is genetically more distinct from the lakes other two Midas cichlid species than is A. xiloaensis-normal. Thus, we have identified sexual isolation based on color that is evident in population genetics and mate choice. Our results suggest that sexual selection through color assortative mating may play an important role in incipient sympatric speciation in Midas cichlids of Nicaragua.

Cryptic diversity and deep divergence in an upper Amazonian leaflitter frog, Eleutherodactylus ockendeni

BackgroundThe forests of the upper Amazon basin harbour some of the worlds highest anuran species richness, but to date we have only the sparsest understanding of the distribution of genetic diversity within and among species in this region. To quantify region-wide genealogical patterns and to test for the presence of deep intraspecific divergences that have been documented in some other neotropical anurans, we developed a molecular phylogeny of the wide-spread terrestrial leaflitter frog Eleutherodactylus ockendeni (Leptodactylidae) from 13 localities throughout its range in Ecuador using data from two mitochondrial genes (16S and cyt b; 1246 base pairs). We examined the relation between divergence of mtDNA and the nuclear genome, as sampled by five species-specific microsatellite loci, to evaluate indirectly whether lineages are reproductively isolated where they co-occur. Our extensive phylogeographic survey thus assesses the spatial distribution of E. ockendeni genetic diversity across eastern Ecuador.ResultsWe identified three distinct and well-supported clades within the Ecuadorean range of E. ockendeni: an uplands clade spanning north to south, a northeastern and central lowlands clade, and a central and southeastern clade, which is basal. Clades are separated by 12% to 15% net corrected p-distance for cytochrome b, with comparatively low sequence divergence within clades. Clades marginally overlap in some geographic areas (e.g., Napo River basin) but are reproductively isolated, evidenced by diagnostic differences in microsatellite PCR amplification profiles or DNA repeat number and coalescent analyses (in MDIV) best modelled without migration. Using Bayesian (BEAST) and net phylogenetic estimates, the Southeastern Clade diverged from the Upland/Lowland clades in the mid-Miocene or late Oligocene. Lowland and Upland clades speciated more recently, in the early or late Miocene.ConclusionOur findings uncover previously unsuspected cryptic species diversity within the common leaflitter frog E. ockendeni, with at least three different species in Ecuador. While these clades are clearly geographically circumscribed, they do not coincide with any existing landscape barriers. Divergences are ancient, from the Miocene, before the most dramatic mountain building in the Ecuadorean Andes. Therefore, this diversity is not a product of Pleistocene refuges. Our research coupled with other studies suggests that species richness in the upper Amazon is drastically underestimated by current inventories based on morphospecies.

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits — including body shape, colour, lips and jaws — is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high‐bodied benthic species (Amphilophus astorquii). A total of 453 genome‐wide informative SNPs were identified in 240 F2 hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy‐two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL‐linked genomic regions, genes that might contribute to divergence in body shape along the benthic–limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.

Pleistocene desiccation in East Africa bottlenecked but did not extirpate the adaptive radiation of Lake Victoria haplochromine cichlid fishes

The Great Lakes region of East Africa, including Lake Victoria, is the center of diversity of the mega-diverse cichlid fishes (Perciformes: Teleostei). Paleolimnological evidence indicates dramatic desiccation of this lake ca. 18,000–15,000 years ago. Consequently, the hundreds of extant endemic haplochromine species in the lake must have either evolved since then or refugia must have existed, within that lake basin or elsewhere, from which Lake Victoria was recolonized. We studied the population history of the Lake Victoria region superflock (LVRS) of haplochromine cichlids based on nuclear genetic analysis (12 microsatellite loci from 400 haplochomines) of populations from Lake Kivu, Lake Victoria, and the connected and surrounding rivers and lakes. Population genetic analyses confirmed that Lake Kivu haplochromines colonized Lake Victoria. Coalescent analyses show a 30- to 50-fold decline in the haplochromine populations of Lake Victoria, Lake Kivu, and the region ca. 18,000–15,000 years ago. We suggest that this coincides with drastic climatic and geological changes in the late Pleistocene. The most recent common ancestor of the Lake Victoria region haplochromines was estimated to have existed about 4.5 million years ago, which corresponds to the first radiation of cichlids in Lake Tanganyika and the origin of the tribe Haplochrominii. This relatively old evolutionary origin may explain the high levels of polymorphism still found in modern haplochromines. This degree of polymorphism might have acted as a “genetic reservoir” that permitted the explosive radiation of hundreds of haplochromines and their array of contemporary adaptive morphologies.

Parallel evolution of Nicaraguan crater lake cichlid fishes via non-parallel routes

Fundamental to understanding how biodiversity arises and adapts is whether evolution is predictable in the face of stochastic genetic and demographic factors. Here we show rapid parallel evolution across two closely related but geographically isolated radiations of Nicaraguan crater lake cichlid fishes. We find significant morphological, ecological and genetic differentiation between ecomorphs in sympatry, reflected primarily in elongated versus high-bodied shape, differential ecological niche use and genetic differentiation. These eco-morphological divergences are significantly parallel across radiations. Based on 442,644 genome-wide single nucleotide polymorphisms, we identify strong support for the monophyly of, and subsequent sympatric divergence within, each radiation. However, the order of speciation differs across radiations; in one lake the limnetic ecomorph diverged first while in the other a benthic ecomorph. Overall our results demonstrate that complex parallel phenotypes can evolve very rapidly and repeatedly in similar environments, probably due to natural selection, yet this evolution can proceed along different evolutionary genetic routes.

A hybrid genetic linkage map of two ecologically and morphologically divergent Midas cichlid fishes (Amphilophus spp.) obtained by massively parallel DNA sequencing (ddRADSeq)

Cichlid fishes are an excellent model system for studying speciation and the formation of adaptive radiations because of their tremendous species richness and astonishing phenotypic diversity. Most research has focused on African rift lake fishes, although Neotropical cichlid species display much variability as well. Almost one dozen species of the Midas cichlid species complex (Amphilophus spp.) have been described so far and have formed repeated adaptive radiations in several Nicaraguan crater lakes. Here we apply double-digest restriction-site associated DNA sequencing to obtain a high-density linkage map of an interspecific cross between the benthic Amphilophus astorquii and the limnetic Amphilophus zaliosus, which are sympatric species endemic to Crater Lake Apoyo, Nicaragua. A total of 755 RAD markers were genotyped in 343 F2 hybrids. The map resolved 25 linkage groups and spans a total distance of 1427 cM with an average marker spacing distance of 1.95 cM, almost matching the total number of chromosomes (n = 24) in these species. Regions of segregation distortion were identified in five linkage groups. Based on the pedigree of parents to F2 offspring, we calculated a genome-wide mutation rate of 6.6 × 10−8 mutations per nucleotide per generation. This genetic map will facilitate the mapping of ecomorphologically relevant adaptive traits in the repeated phenotypes that evolved within the Midas cichlid lineage and, as the first linkage map of a Neotropical cichlid, facilitate comparative genomic analyses between African cichlids, Neotropical cichlids and other teleost fishes.

Genomics of adaptation and speciation in cichlid fishes: recent advances and analyses in African and Neotropical lineages.

Cichlid fishes are remarkably phenotypically diverse and species-rich. Therefore, they provide an exciting opportunity for the study of the genetics of adaptation and speciation by natural and sexual selection. Here, we review advances in the genomics and transcriptomics of cichlids, particularly regarding ecologically relevant differences in body shape, trophic apparatus, coloration and patterning, and sex determination. Research conducted so far has focused almost exclusively on African cichlids. To analyse genomic diversity and selection in a Neotropical radiation, we conducted a comparative transcriptomic analysis between sympatric, ecologically divergent crater-lake Midas cichlids (Lake Xiloá Amphilophus amarillo and Amphilophus sagittae). We pyrosequenced (Roche 454) expressed sequence tag (EST) libraries and generated more than 178 000 000 ESTs and identified nine ESTs under positive selection between these sister species (Ka/Ks > 1). None of these ESTs were found to be under selection in African cichlids. Of 11 candidate genes for ecomorphological differentiation in African cichlids, none showed signs of selection between A. amarillo and A. sagittae. Although more population-level studies are now needed to thoroughly document patterns of divergence during speciation of cichlids, available information so far suggests that adaptive phenotypic diversification in Neotropical and African cichlids may be evolving through non-parallel genetic bases.

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ∼ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.