Parallel evolution in Ugandan crater lakes: repeated evolution of limnetic body shapes in haplochromine cichlid fish

Evolutionary dynamics of pre- and postzygotic reproductive isolation in cichlid fishes

Cichlid fishes are exceptionally species-rich, speciated at explosive rates and, hence, are a model system in speciation research. Yet, their reproductive isolating barriers have, so far, not been comprehensively studied. Here, we review current knowledge on pre- and postzygotic mechanisms in cichlids. While premating isolation is the norm in cichlids, its strength varies across lineages and with the geographical setting. Moreover, manipulations of ambient conditions tended to reduce assortative mating among closely related species, suggesting that premating isolation in cichlids is often fragile and context dependent. The observed lack of complete reproductive isolation is supported by past and present hybridization events that have contributed to diversity by creating novel allelic combinations. On the other hand, our meta-analysis highlights that intrinsic postzygotic isolation might accumulate faster than assumed. Mild forms of genetic incompatibilities, such as sex ratio distortion, can already be observed among closely related species. Therefore, cessation of gene flow by strong reproductive isolation in cichlids requires a combination of premating prezygotic isolation supplemented with intrinsic and extrinsic postzygotic barriers. Further, we suggest crucial next steps to improve our knowledge about reproductive barriers in cichlids to understand the evolutionary dynamics of pre- and postzygotic isolation mechanisms during adaptive radiations. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Tracking parallel adaptation of shell morphology through geological times in the land snail genus Pupilla (Gastropoda: Stylommatophora: Pupillidae)

Changing environmental conditions force species either to disperse or to adapt locally either genetically or via phenotypic plasticity. Although limits of plasticity can be experimentally tested, the predictability of genetic adaptation is restricted due to its stochastic nature. Nevertheless, our understanding of evolutionary adaptation has been improving in particular through studies of parallel adaptation. Based on molecular phylogenetic inferences and morphological investigations of both recent and fossil shells we tracked the morphological changes in three land snails, Pupilla alpicola, Pupilla loessica and Pupilla muscorum. These species differ in habitat requirements as well as historical and extant distributions with P. alpicola and P. loessica being more similar to each other than to P. muscorum. Therefore, we hypothesized, that the three species reacted independently and individually to the conditions changing throughout the Pleistocene, but expected that changes within P. alpicola and P. loessica would be more similar compared to P. muscorum. Indeed, intraspecific shell shape differences across time were similar in P. alpicola and P. loessica, suggesting that similar niche shifts have led to similar transformations in parallel. In contrast, extant P. muscorum populations were practically identical in shape to their ancestors. They have probably tracked their ecological niches through time.

Molecular identification of Bulinus spp. intermediate host snails of Schistosoma spp. in crater lakes of western Uganda with implications for the transmission of the Schistosoma haematobium group parasites

Background

Human schistosomiasis is the second most important tropical disease and occurs in two forms in Africa (intestinal and urogenital) caused by the digenetic trematodes Schistosoma mansoni and Schistosoma haematobium, respectively. A proposed recent shift of schistosomiasis above a previously established altitudinal threshold of 1400 m above sea level in western Ugandan crater lakes has triggered more research interest there.

Methods

Based on extensive field sampling in western Uganda and beyond and employing an approach using sequences of the mitochondrial barcoding gene cytochrome c oxidase subunit 1 (cox1) this study aims were: (i) identification and establishment of the phylogenetic affinities of Bulinus species as potential hosts for Schistosoma spp.; (ii) determining diversity, frequency and distribution patterns of Bulinus spp.; and (iii) establishing genetic variability and phylogeographical patterns using Bayesian inference and parsimony network analyses.

Results

Out of the 58 crater lakes surveyed, three species of Bulinus snails were found in 34 crater lakes. Bulinus tropicus was dominating, Bulinus forskalii was found in two lakes and Bulinus truncatus in one. The latter two species are unconfirmed potential hosts for S. haematobium in this region. However, Bulinus tropicus is an important species for schistosomiasis transmission in ruminants. Bulinus tropicus comprised 31 haplotypes while both B. forskalii and B. truncatus exhibited only a single haplotype in the crater lakes. All species clustered with most of the haplotypes from surrounding lake systems forming source regions for the colonization of the crater lakes.

Conclusions

This first detailed malacological study of the crater lakes systems in western Uganda revealed presence of Bulinus species that are either not known or not regionally known to be hosts for S. haematobium, the causing agent of human urogenital schistosomiasis. Though this disease risk is almost negligible, the observed dominance of B. tropicus in the crater lakes shows that there is a likelihood of a high risk of infections with Schistosoma bovis. Thus, extra attention should be accorded to safeguard wild and domestic ruminants in this region as the population benefits from these animals.

Searching for Sympatric Speciation in the Genomic Era

Sympatric speciation illustrates how natural and sexual selection may create new species in isolation without geographic barriers. However, recent genomic reanalyses of classic examples of sympatric speciation reveal complex histories of secondary gene flow from outgroups into the radiation. In contrast, the rich theoretical literature on this process distinguishes among a diverse range of models based on simple genetic histories and different types of reproductive isolating barriers. Thus, there is a need to revisit how to connect theoretical models of sympatric speciation and their predictions to empirical case studies in the face of widespread gene flow. Here, theoretical differences among different types of sympatric speciation and speciation-with-gene-flow models are reviewed and summarized, and genomic analyses are proposed for distinguishing which models apply to case studies based on the timing and function of adaptive introgression. Investigating whether secondary gene flow contributed to reproductive isolation is necessary to test whether predictions of theory are ultimately borne out in nature.

Island- and lake-like parallel adaptive radiations replicated in rivers

Parallel adaptive radiations have arisen following the colonization of islands by lizards and lakes by fishes. In these classic examples, parallel adaptive radiation is a response to the ecological opportunities afforded by the colonization of novel ecosystems and similar adaptive landscapes that favour the evolution of similar suites of ecomorphs, despite independent evolutionary histories. Here, we demonstrate that parallel adaptive radiations of cichlid fishes arose in South American rivers. Speciation-assembled communities of pike cichlids (Crenicichla) have independently diversified into similar suites of novel ecomorphs in the Uruguay and Paraná Rivers, including crevice feeders, periphyton grazers and molluscivores. There were bursts in phenotypic evolution associated with the colonization of each river and the subsequent expansion of morphospace following the evolution of the ecomorphs. These riverine clades demonstrate that characteristics emblematic of textbook parallel adaptive radiations of island- and lake-dwelling assemblages are feasible evolutionary outcomes even in labile ecosystems such as rivers.

Lessons from a natural experiment: Allopatric morphological divergence and sympatric diversification in the Midas cichlid species complex are largely influenced by ecology in a deterministic way

Explaining why some lineages diversify while others do not and how are key objectives in evolutionary biology. Young radiations of closely related species derived from the same source population provide an excellent opportunity to disentangle the relative contributions of possible drivers of diversification. In these settings, lineage-specific effects are shared and can be ruled out. Moreover, the relevant demographic and ecological parameters can be estimated accurately. Midas cichlid fish in Nicaragua have repeatedly colonized several crater lakes, diverged from the same source populations, and, interestingly, diversified in some of them but not others. Here, using the most comprehensive molecular and geometric morphometric data set on Midas cichlids to date (∼20,000 SNPs, 12 landmarks, ∼700 individuals), we aim to understand why and how crater lake populations diverge and why some of them are more prone to diversify in sympatry than others. Taking ancestor-descendant relationships into account, we find that Midas cichlids diverged in parallel from their source population mostly-but not exclusively-by evolving more slender body shapes in all six investigated crater lakes. Admixture among crater lakes has possibly facilitated this process in one case, but overall, admixture and secondary waves of colonization cannot predict morphological divergence and intralacustrine diversification. Instead, morphological divergence is larger the more dissimilar a crater lake is compared to the source lake and happens rapidly after colonization followed by a slow-down with time. Our data also provide some evidence that founder effects may positively contribute to divergence. The depth of a crater lake is positively associated with variation in body shapes (and number of species), presumably by providing more ecological opportunities. In conclusion, we find that parallel morphological divergence in allopatry and the propensity for diversification in sympatry across the entire Midas cichlid fish radiation is partly predictable and mostly driven by ecology.

Ecological opportunity alters the timing and shape of adaptive radiation

The uneven distribution of diversity is a conspicuous phenomenon across the tree of life. Ecological opportunity is a prominent catalyst of adaptive radiation and therefore may alter patterns of diversification. We evaluated the distribution of shifts in diversification rates across the cichlid phylogeny and the distribution of major clades across phylogenetic space. We also tested if ecological opportunity influenced these patterns. Colonization-associated ecological opportunity altered the tempo and mode of diversification during the adaptive radiation of cichlid fishes. Clades that arose following colonization events diversified faster than other clades. Speciation rate shifts were nonrandomly distributed across the phylogeny such that they were disproportionally concentrated around nodes that corresponded with colonization events (i.e., of continents, river basins, or lakes). Young clades tend to expand faster than older clades; however, colonization-associated ecological opportunity accentuated this pattern. There was an interaction between clade age and ecological opportunity that explained the trajectory of clades through phylogenetic space over time. Our results indicate that ecological opportunities afforded by continental and ecosystem-scale colonization events explain the dramatic speciation rate heterogeneity and phylogenetic imbalance that arose during the evolutionary history of cichlid fishes.

The Integrated Genomic Architecture and Evolution of Dental Divergence in East African Cichlid Fishes (Haplochromis chilotes x H. nyererei)

The independent evolution of the two toothed jaws of cichlid fishes is thought to have promoted their unparalleled ecological divergence and species richness. However, dental divergence in cichlids could exhibit substantial genetic covariance and this could dictate how traits like tooth numbers evolve in different African Lakes and on their two jaws. To test this hypothesis, we used a hybrid mapping cross of two trophically divergent Lake Victoria species (Haplochromis chilotes × Haplochromis nyererei) to examine genomic regions associated with cichlid tooth diversity. Surprisingly, a similar genomic region was found to be associated with oral jaw tooth numbers in cichlids from both Lake Malawi and Lake Victoria. Likewise, this same genomic location was associated with variation in pharyngeal jaw tooth numbers. Similar relationships between tooth numbers on the two jaws in both our Victoria hybrid population and across the phylogenetic diversity of Malawi cichlids additionally suggests that tooth numbers on the two jaws of haplochromine cichlids might generally coevolve owing to shared genetic underpinnings. Integrated, rather than independent, genomic architectures could be key to the incomparable evolutionary divergence and convergence in cichlid tooth numbers.

Genome-wide analysis of UDP-glycosyltransferase super family in Brassica rapa and Brassica oleracea reveals its evolutionary history and functional characterization

Background

Glycosyltransferases comprise a highly divergent and polyphyletic multigene family that is involved in widespread modification of plant secondary metabolites in a process called glycosylation. According to conserved domains identified in their amino acid sequences, these glycosyltransferases can be classified into a single UDP-glycosyltransferase (UGT) 1 superfamily.

Results

We performed genome-wide comparative analysis of UGT genes to trace evolutionary history in algae, bryophytes, pteridophytes, and angiosperms; then, we further investigated the expansion mechanisms and function characterization of UGT gene families in Brassica rapa and Brassica oleracea. Using Hidden Markov Model search, we identified 3, 21, 140, 200, 115, 147, and 147 UGTs in Chlamydomonas reinhardtii, Physcomitrella patens, Selaginella moellendorffii, Oryza sativa, Arabidopsis thaliana, B. rapa, and B. oleracea, respectively. Phylogenetic analysis revealed that UGT80 gene family is an ancient gene family, which is shared by all plants and UGT74 gene family is shared by ferns and angiosperms, but the remaining UGT gene families were shared by angiosperms. In dicot lineage, UGTs among three species were classified into three subgroups containing 3, 6, and 12 UGT gene families. Analysis of chromosomal distribution indicates that 98.6 and 71.4% of UGTs were located on B. rapa and B. oleracea pseudo-molecules, respectively. Expansion mechanism analyses uncovered that whole genome duplication event exerted larger influence than tandem duplication on expansion of UGT gene families in B. rapa, and B. oleracea. Analysis of selection forces of UGT orthologous gene pairs in B. rapa, and B. oleracea compared to A. thaliana suggested that orthologous genes in B. rapa, and B. oleracea have undergone negative selection, but there were no significant differences between A. thaliana –B. rapa and A. thaliana –B. oleracea lineages. Our comparisons of expression profiling illustrated that UGTs in B. rapa performed more discrete expression patterns than these in B. oleracea indicating stronger function divergence. Combing with phylogeny and expression analysis, the UGTs in B. rapa and B. oleracea experienced parallel evolution after they diverged from a common ancestor.

Conclusion

We first traced the evolutionary history of UGT gene families in plants and revealed its evolutionary and functional characterization of UGTs in B. rapa, and B. oleracea. This study provides novel insights into the evolutionary history and functional divergence of important traits or phenotype-related gene families in plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3844-x) contains supplementary material, which is available to authorized users.

On the maintenance of genetic variation and adaptation to environmental change: considerations from population genomics in fishes

The first goal of this paper was to overview modern approaches to local adaptation, with a focus on the use of population genomics data to detect signals of natural selection in fishes. Several mechanisms are discussed that may enhance the maintenance of genetic variation and evolutionary potential, which have been overlooked and should be considered in future theoretical development and predictive models: the prevalence of soft sweeps, polygenic basis of adaptation, balancing selection and transient polymorphisms, parallel evolution, as well as epigenetic variation. Research on fish population genomics has provided ample evidence for local adaptation at the genome level. Pervasive adaptive evolution, however, seems to almost never involve the fixation of beneficial alleles. Instead, adaptation apparently proceeds most commonly by soft sweeps entailing shifts in frequencies of alleles being shared between differentially adapted populations. One obvious factor contributing to the maintenance of standing genetic variation in the face of selective pressures is that adaptive phenotypic traits are most often highly polygenic, and consequently the response to selection should derive mostly from allelic co-variances among causative loci rather than pronounced allele frequency changes. Balancing selection in its various forms may also play an important role in maintaining adaptive genetic variation and the evolutionary potential of species to cope with environmental change. A large body of literature on fishes also shows that repeated evolution of adaptive phenotypes is a ubiquitous evolutionary phenomenon that seems to occur most often via different genetic solutions, further adding to the potential options of species to cope with a changing environment. Moreover, a paradox is emerging from recent fish studies whereby populations of highly reduced effective population sizes and impoverished genetic diversity can apparently retain their adaptive potential in some circumstances. Although more empirical support is needed, several recent studies suggest that epigenetic variation could account for this apparent paradox. Therefore, epigenetic variation should be fully integrated with considerations pertaining to role of soft sweeps, polygenic and balancing selection, as well as repeated adaptation involving different genetic basis towards improving models predicting the evolutionary potential of species to cope with a changing world.

Incipient sympatric speciation in Midas cichlid fish from the youngest and one of the smallest crater lakes in Nicaragua due to differential use of the benthic and limnetic habitats?

Understanding how speciation can occur without geographic isolation remains a central objective in evolutionary biology. Generally, some form of disruptive selection and assortative mating are necessary for sympatric speciation to occur. Disruptive selection can arise from intraspecific competition for resources. If this competition leads to the differential use of habitats and variation in relevant traits is genetically determined, then assortative mating can be an automatic consequence (i.e., habitat isolation). In this study, we caught Midas cichlid fish from the limnetic (middle of the lake) and benthic (shore) habitats of Crater Lake Asososca Managua to test whether some of the necessary conditions for sympatric speciation due to intraspecific competition and habitat isolation are given. Lake As. Managua is very small (<900 m in diameter), extremely young (maximally 1245 years of age), and completely isolated. It is inhabited by, probably, only a single endemic species of Midas cichlids, Amphilophus tolteca. We found that fish from the limnetic habitat were more elongated than fish collected from the benthic habitat, as would be predicted from ecomorphological considerations. Stable isotope analyses confirmed that the former also exhibit a more limnetic lifestyle than the latter. Furthermore, split‐brood design experiments in the laboratory suggest that phenotypic plasticity is unlikely to explain much of the observed differences in body elongation that we observed in the field. Yet, neutral markers (microsatellites) did not reveal any genetic clustering in the population. Interestingly, demographic inferences based on RAD‐seq data suggest that the apparent lack of genetic differentiation at neutral markers could simply be due to a lack of time, as intraspecific competition may only have begun a few hundred generations ago.

Do relaxed selection and habitat temperature facilitate biased mitogenomic introgression in a narrowly endemic fish?

Introgression might be exceptionally common during the evolution of narrowly endemic species. For instance, in the springs of the small and isolated Cuatro Ciénegas Valley, the mitogenome of the cichlid fish Herichthys cyanoguttatus could be rapidly introgressing into populations of the trophically polymorphic H. minckleyi. We used a combination of genetic and environmental data to examine the factors associated with this mitochondrial introgression. A reduced representation library of over 6220 single nucleotide polymorphisms (SNPs) from the nuclear genome showed that mitochondrial introgression into H. minckleyi is biased relative to the amount of nuclear introgression. SNP assignment probabilities also indicated that cichlids with more hybrid ancestry are not more commonly female providing no support for asymmetric backcrossing or hybrid‐induced sex‐ratio distortion in generating the bias in mitochondrial introgression. Smaller effective population size in H. minckleyi inferred from the SNPs coupled with sequences of all 13 mitochondrial proteins suggests that relaxed selection on the mitogenome could be facilitating the introgression of “H. cyanoguttatus” haplotypes. Additionally, we showed that springs with colder temperatures had greater amounts of mitochondrial introgression from H. cyanoguttatus. Relaxed selection in H. minckleyi coupled with temperature‐related molecular adaptation could be facilitating mitogenomic introgression into H. minckleyi.