Philosophical Transactions of the Royal Society B: Biological Sciences
You have accessReview articles

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

Sina J. Rometsch

Sina J. Rometsch

Department of Biology, University of Konstanz, Universitätsstraße 10, 78457 Konstanz, Germany

Google Scholar

Find this author on PubMed

,
Julián Torres-Dowdall

Julián Torres-Dowdall

Department of Biology, University of Konstanz, Universitätsstraße 10, 78457 Konstanz, Germany

Google Scholar

Find this author on PubMed

and
Axel Meyer

Axel Meyer

Department of Biology, University of Konstanz, Universitätsstraße 10, 78457 Konstanz, Germany

[email protected]

Google Scholar

Find this author on PubMed

    Abstract

    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'.

    1. Introduction

    Cichlid fishes are one of the most species-rich and phenotypically diverse families of vertebrates (e.g. reviewed in [1]). Approximately 10% of all fish species are cichlids and astonishing inter- and intraspecific variation can be observed in almost every trait including morphology and behaviour, but particularly in coloration, which is often sexually dimorphic [25]. Cichlids are not only exceptionally species-rich but also speciated at explosive rates. For example, around 250 species evolved in Lake Tanganyika in 10–12 Myr, more than 800 species are endemic to Lake Malawi originating in less than 4 Myr and around 500 species evolved in Lake Victoria in only 15 000–100 000 years [4,6,7]. In the crater lakes of Nicaragua, several species originated even in only a few hundred generations and are less than 2000 years old [8]. Considering these extremely rapid rates at which some radiations evolved in allopatry, and most remarkably speciation also occurred in sympatry, as mainly demonstrated for smaller lake systems (e.g. [810]), it is clear why cichlid fishes constitute a model system for the study of speciation and formation of adaptive radiations in different geographical settings and temporal scales.

    Hybridization in cichlids has been observed under laboratory conditions and inferred in nature [8,1113]. The ease with which often fertile hybrids can be generated is one of the reasons why cichlids are increasingly used to study the genomics of adaptation through forward genetic approaches such as quantitative trait loci-analyses (e.g. [1416]). More relevant, hybrids between nominal cichlid species are also found in nature [8,1720], suggesting that reproductive isolation barriers can be ‘leaky’. Evidence has accumulated pointing to an important role of hybridization in catalysing cichlid diversification in some adaptive radiations [8,11,2127]. However, incomplete reproductive isolation and ongoing hybridization in cichlids raise the question of what prevents divergent genetic clusters from collapsing and thereby disrupting or reversing the process of speciation in this lineage?

    Strong reproductive isolation is not necessarily the result of one specific barrier completely ceasing gene flow between divergent genetic clusters, but potentially the combination of multiple imperfect isolation barriers [28]. To identify the individual contributions of these barriers to the speciation of cichlid fishes, we review the literature on reproductive isolation and address the following questions: (i) is strong premating isolation the norm in cichlids, does its strength vary between species complexes belonging to different radiations, what is the role of the geographical setting of speciation and what are the cues affecting it? (ii) what is the evidence that prezygotic postmating barriers contribute to reproductive isolation in cichlids? and (iii) given the evolutionary youth of many cichlids, particularly in crater lake radiations, but also in young East African Rift Valley lakes, what is the contribution of intrinsic and extrinsic postzygotic isolation mechanisms to speciation?

    2. The importance of geographical settings for the evolution of reproductive isolation

    The geographical setting of speciation (i.e. sympatry or allopatry) has important implications for the establishment of reproductive isolation as it can contribute to pre- and postzygotic isolation in different ways [29]. Temporary spatial isolation, as expected under the model of traditional allopatric speciation, was suggested repeatedly to contribute to the astonishing species richness in African cichlids [3032]. Periods of allopatry have been common in all African Great Lakes owing to fluctuations in the water level, resulting in habitat fragmentation and limiting gene flow [33,34]. Most researches agree on the importance of such isolation events and/or limited dispersal, owing to strong philopatry that cichlids seem to exhibit, for building up genetic divergence between allopatric populations. Evidence for this is provided by phylogeographic studies that tend to find strong isolation by distance within radiations of African cichlids (e.g. [3538], but see [39]). Consequently, as a by-product of evolutionary divergence, incompatibilities between species accumulate (i.e. Dobzhansky–Muller incompatibilities, DMIs), that can be manifested in pre- and postzygotic isolation barriers [40]. However, most commonly DMIs have been implicated in postzygotic intrinsic isolation owing to the potentially detrimental effects of incompatible alleles being brought together in interspecific hybrids [40]. Theoretically, assuming that mate recognition systems are not under stabilizing selection, allopatric divergence can also favour the evolution of prezygotic isolation [41]. On the other hand, in sympatrically evolving species or in allopatric species coming into secondary contact, reinforcement through natural selection against hybrids with decreased fitness can be a powerful mechanism strengthening prezygotic barriers to gene flow [41], while the evolution of intrinsic postzygotic barriers has been suggested to be rather unlikely under conditions of continuous gene flow ([42], but see [4345]). As cichlid fishes represent a study system comprising closely related species evolving in allopatry and sympatry, investigating reproductive isolation in cichlids allows us to consider the importance of biogeographic settings on speciation mechanisms, to determine their influences on the relative prevalence and strength of different isolation barriers.

    3. Prezygotic isolation

    (a) Premating isolation

    Premating isolation encompasses all reproductive barriers that diminish heterospecific interactions between sexes to avoid hybridization [46]. In phylogenetically recent lineages, such as cichlids, premating barriers would be expected to be of particular importance for reproductive isolation [41,46]. Disruptive sexual selection based on intraspecific variation is argued to be one of the main contributors to the rapid rate of speciation seen in some cichlid lineages and to play an important role in maintaining cichlid species diversity in sympatry [6,4749].

    Here, we review the literature asking how commonly premating isolation is observed in cichlids, if it varies across lineages and among sympatric and allopatric species, and what influences its relative strength. A search in ‘Web of Science’ using keywords related to premating isolation in cichlids returned a total of 802 studies (electronic supplementary material, text and table S1). After removing theoretical papers, review articles and studies only citing cichlid literature but using focal species outside of the family Cichlidae, 497 relevant studies remained (electronic supplementary material, text and table S1). Additionally, we added 21 studies that investigated assortative mating in cichlids but were not included in the results of any of these searches. We filtered those studies for direct tests of assortative mating by species or populations. This reduced the number of relevant studies to 39 (electronic supplementary material, table S2), which we subsequently used to conduct a meta-analysis using OpenMEE [50] to investigate premating isolation (for more details, see the electronic supplementary material). Most studies (29 out of 39) reported count data on actual mating events (e.g. spawning events, egg laying, paternity analysis) which we used to conduct a phylogenetically controlled meta-analysis assuming a Brownian motion model (lambda fitted with λ = 1; figure 1). Overall, strong premating isolation seems to be the norm in cichlids (figure 1a). However, heterogeneity (I2, percentage of variation owing to study heterogeneity, [57]) across the investigated studies was high (I2 = 83.98, p < 0.001). There was significant phylogenetic variance (τ2λ = 1.3712) and including phylogenetic correlations as random effects significantly improved the model, as it reduced the amount of unexplained between study variance (τ2 = 0.796 versus τ2 = 0.982 when not controlling for phylogeny). Next, we explored the drivers of this variation. While the strength of premating barriers differed among species complexes, evidence at the meta-analysis level for assortative mate choice was observed in Neotropical lineages (Amphilophus cf. citrinellus and Apistogramma spp.), as well as in the African lineages Maylandia cf. zebra and Pundamilia cf. nyererei, with Tropheus cf. moorii from Lake Tanganyika being the only exception, where mating did not differ from random (figure 1b). All included studies that investigated mate choice in Tropheus cf. moorii were conducted in allopatric populations within species. As it is expected that premating isolation increases with genetic divergence and is favoured in sympatric settings through reinforcement, low levels thereof might potentially explain the absence of assortative mating in the Tropheus species complex. Interestingly, even upon human-induced secondary contact of multiple Tropheus populations differing in colour morphs, there was no evidence of reinforcement but rather signs of introgressive hybridization [58]. Accordingly, we determined if the geographical setting of speciation generally predicts the degree of premating isolation and if the strength of assortative mating differs between nominal species and populations within species (figure 1c). Sympatric and allopatric species differed only slightly in the degree of assortative mating. However, species that occur in allopatry showed more variation in the probability of mating with con- versus heterospecific individuals compared to sympatric ones. By contrast, while sympatric populations (i.e. polymorphisms among populations within species in traits thought to affect mate choice, such as coloration) mated as assortatively as sympatric nominal species, the strength of assortative mating was the weakest among allopatric populations. This might indicate that reinforcement could be especially important for the establishment of premating isolation in populations that have not yet evolved substantial genetic divergence [59]. Some studies reported preference scores as a behavioural proxy for mate choice, which we evaluated in a separate meta-analysis (electronic supplementary material, figure S1). Preference scores usually consider differences in female responsiveness to con- versus heterospecific male behaviour (e.g. quivers, lateral displays). Overall, meta-analyses based on preference scores and count data provided the same qualitative trend: females preferred conspecific over heterospecific males, but there was more variation across studies when considering preference scores. The only lineage differing in the pattern of premating isolation between the analyses was Tropheus cf. moorii. While there was no evidence of assortative mating when considering count data in the Tropheus species complex (figure 1b), the results of the meta-analysis using preference scores were more variable and, in some populations, females did show a preference to interact with conspecific males (electronic supplementary material, figure S1). These differences between count data and preference scores in Tropheus cf. moorii might reflect real differences between populations. However, the large variation across studies might also indicate that comparisons of preference scores across studies can be difficult owing to the differences in measuring female preference.

    Figure 1.

    Figure 1. Meta-analysis investigating premating isolation in cichlids. Premating isolation in cichlids (a) seems to be the norm as overall mating is assortative, (b) varies among lineages (all lineages comprise both, allopatric and sympatric comparisons except Tropheus cf. moori, which only includes allopatric comparisons) and (c) its strength is influenced by the geographical setting and genetic divergence. (d) Manipulative studies show that mate choice is context dependent. Log odds ratios depicted with 95% confidence intervals (CI) indicate the likelihood of premating isolation. (ac) Ratios of 1 indicate mating with con- versus heterospecifics is equally likely, values higher than 1 express the fold increase in the likelihood of mating assortatively and values lower than 1 express the fold increase in mating disassortatively. (d) Ratios of 1 indicate mating with con- versus heterospecifics is equally likely as without manipulation, values higher than 1 express the fold increase in the likelihood of mating assortatively after manipulation and values lower than 1 express the fold increase in mating disassortatively after manipulation; n, number of contrasts (this can deviate from the number of studies included, as studies can include more than one contrast); I2, percentage of variability that is owing to heterogeneity across studies rather than sampling variance, identities correspond to studies from which we extracted corresponding data: 1, [51]; 2, [52]; 3, [53]; 4, [54]; 5, [55]; 6, [56]. For more detailed representations of meta-analyses depicted in (b,c), see the electronic supplementary material, figures S2 and S3, respectively.

    Given the substantial diversity of the cichlid clade, large variation in the overall degree of premating isolation (figure 1a) is not surprising. Accordingly, it is unlikely that factors contributing to premating isolation can be generalized across all radiations. Below, we discuss in more detail the different cues contributing to the establishment of premating isolation and variation in their importance for assortative mating across lineages.

    (b) Mate choice: the importance of different cues

    (i) Visual cues: coloration and pattern

    Cichlids are renowned for their astonishing diversity in body coloration, with impressive variation in hue and pattern, and it has been repeatedly determined that visual cues are pivotal for their mating behaviour [51,6063]. Illustrating this diversity, nuptial coloration is positively associated with species richness at the phylogenetic level [64], often strongly varies between closely related species [65], and sometimes even considerable intraspecific variation among different populations can be observed [49,66,67]. Moreover, striking nuptial coloration in dimorphic species conveys information about the reproductive state, quality and status of the male and, therefore, plays a fundamental role in mate choice and sexual selection [6870]. Therefore, it is not surprising that attention in cichlid fish research has been focused predominantly on the role of visual cues in mate choice [65,67,71].

    Intraspecific sexual selection based on male coloration has been clearly established for some species that are polymorphic for coloration [62,67,72]. Several studies have shown that this extends to interspecific mate choice, with females tending to prefer males exhibiting the coloration of their conspecific males [49,72,73], as indicated by our meta-analysis (figure 1a). This is further supported as females prefer heterospecific males of the most similar species in the absence of conspecific males [65], or when exposed to hybrids segregating for nuptial coloration, females choose males that resemble their conspecifics [74,75]. Moreover, species or populations divergent in coloration have repeatedly been found to have high degrees of assortative mating (e.g. [7579]). However, these results may not apply universally (figure 1b). Sympatric and allopatric populations of several African cichlid species that differ in their body coloration show reduced degrees of assortative mating (e.g. [52,80,81]).

    Taken together, these studies highlight the importance of nuptial coloration but also indicate that their contribution to assortative mating might vary across different cichlid lineages. Given that the majority of studies investigating the effect of coloration on premating isolation have been conducted in only a small subset of species (M. cf. zebra from Lake Malawi, P. cf. nyererei from Lake Victoria, T. cf. moorii from Lake Tanganyika and Amphilophus. cf. citrinellus from the Nicaraguan crater lakes, figure 1b; electronic supplementary material, table S2), it is not clear to which degree those results can be generalized to other lineages.

    The importance of visual cues during mate choice and the magnificent variation in nuptial coloration suggest that visual perception might be important for reproductive isolation ([5]; figure 1d). Indeed, there is an impressive amount of diversity in visual sensitivities across cichlid fishes [82]. One of the best studied system in this regard is the Pundamilia complex, where light environment, nuptial coloration and visual sensitivity are strongly associated [83,84]. These associations have been interpreted in the context of sensory drive speciation, where variation in light conditions selected for different alleles of the lws gene that encodes the protein component of the red sensitive photoreceptors of the retina. In turn, this affects female preference for different male nuptial coloration. While female preference is correlated with the genotype at the lws locus [84], no association between opsin gene expression and female preference for male nuptial coloration has been found [85]. Moreover, introgression at this locus following hybridization events has been linked to the divergence among Lake Victoria cichlids [86,87]. Whether these associations hold in other cichlid lineages remains unclear. While in many radiations, variation in visual sensitivity is based on differences in opsin expression patterns [88], a general link between opsin expression and nuptial coloration has not been observed [89] and in some lineages, ecological factors rather than nuptial coloration seem to drive opsin expression [21,90].

    (ii) Olfactory cues

    Compared to the well-investigated nature of visual cues, our understanding of olfactory signals in cichlids is far less complete, partially because compounds relevant to species isolating mate choice have not yet been identified [91]. Nonetheless, it has been argued that the olfactory system plays an important role in various cichlid behaviours including mate choice [91]. The composition of chemical cues in cichlids is complex and only anecdotal evidence suggests differences that could confer species identity (e.g. major histocompatibility complex, [55,9294]). Males and females use variation in their chemical compounds to communicate physical condition, reproductive state and initiate courtship behaviour [9597]. Additionally, the composition of odorants provides a unique signature that is important for individual recognition [98,99] and during sexual imprinting [78,100].

    Although chemical cues alone appear not to be sufficient to drive assortative mating in cichlids [51,61,101], it has been suggested that visual signals often have to be supplemented by other cues, such as chemical ones, to guide mate choice ([52,55]; figure 1d). In fact, conspecific odour in the presence of random visual stimuli (i.e. heterospecific individuals) can trigger courtship. Male Astatotilapia burtoni will even court distantly related juvenile Tilapia mariae when their tank is supplemented with water conditioned by a gravid A. burtoni female, but will ignore them in the absence of such chemical cues [102]. Even in cases where visual signals are sufficient to initiate mating, olfactory cues may in some species be required to trigger behavioural and hormonal responses crucial for reproduction [95,96]. However, the role of odour in premating isolation rests on the assumption that cichlids can distinguish between conspecific and heterospecific cues, which to our knowledge has not been tested, at least not systematically. Experiments using a set-up as reported by Giaquinto et al. [97] could help to address this question and provide crucial information needed to further evaluate the importance of olfactory cues for premating isolation in cichlids.

    (iii) Acoustic cues

    Cichlids are capable of producing a variety of different sounds involving the pharyngeal jaw, stridulation or body movements [103105]. Variation in acoustic parameters is observed across species, among different populations within species and sometimes even high levels of inter- and intraindividual variation are reported [105108]. Several studies show that cichlids are sensitive to the range of frequencies of the sounds they produce, but compared to goldfish (Carassius auratus), they appear to have poor sensitivity [109111]. However, some cichlid species possess specialized connections between the swim bladder and the inner ear that improve their auditory abilities [112]. The role of auditory cues in mate choice is implied by the fact that male vocalization is especially prominent during courtship and female sensitivity to low-frequency sounds, as emitted by males during displays, is increased when they are ready to spawn [113,114].

    While cichlids are known to show pronounced interspecific differences in their acoustic parameters, most studies assume, but do not demonstrate, that females are actually capable of discriminating between con- or heterospecifics acoustic signals [105,106,114]. Playback experiments show that females perceive acoustic cues and preferably engage in courtship behaviour with males that are associated with sound playback over mere visual stimuli ([115]; figure 1d). Further, both sexes can discriminate between conspecific sounds and bursts of white noise [56]. However, if the specificity of sound perception allows conspecific to be distinguished from heterospecific vocalizations is yet unknown. Thus, predictions about the relative importance of acoustic cues, compared to visual and olfactory ones, in mating are difficult to make and might differ between species. As with chemical signals, sounds alone appear not to be sufficient for mate attraction and only constitute effective signals in combination with visual/olfactory cues [56,107]. Further evaluation of the contribution of acoustic cues to premating isolation will depend on future research investigating sound discrimination capability of cichlids. Playback experiments, as those conducted by Verzijden et al. [115], but testing con- versus heterospecific sounds, will be helpful to elucidate this issue.

    (iv) Multimodality: the interplay of different cues

    Integrating the knowledge on which cues govern premating isolation in cichlids makes it evident that species assortative mating is rarely based on a single cue, but rather relies on the combination of multiple cues ([52,55], but see [54]). However, depending on the focal species, the importance of these cues might differ and a hierarchical structure is likely [116]. Studies that manipulated various conditions allow us to disentangle the contribution of the different cues for premating isolation (figure 1d). A large percentage of the tests we evaluated, which manipulated or removed specific cues, resulted in reduced premating isolation. However, the effect varied across cues and species. For example, the complete impairment of visual cues led to the reduction in assortative mating without exception in all tested species. On the other hand, effects of manipulations masking colour differences between species by using monochromatic light strongly depended on the tested lineage. While females of P. cf. nyererei could not distinguish between conspecific and heterospecific males any longer ([54]; figure 1d, [117]), it had no effect on assortative mating in M. cf. zebra ([51,52]; figure 1d). This indicates that the importance of cues might differ, as in some cases not coloration per se but rather colour pattern might be important ([118]; figure 1d). By contrast, the effects of impeding olfactory and auditory cues were more variable and less pronounced compared to visual ones, suggestive of an overall smaller relevance of non-visual cues.

    Importantly, some of these studies exemplify that premating barriers are context dependent. Under the conditions in which the respective cues evolved, they seem to be a powerful driver for premating isolation. However, drastic environmental changes that constrain the perception of some cues may exacerbate mate choice, increasing interspecific gene flow. Eutrophication and associated turbidity in some regions of Lake Victoria, for example, can result in hybridization between the sympatric species pair P. cf. pundamilia and P. cf. nyererei, while these species mate strongly assortatively in clear water where visual signals are not affected by turbidity [18].

    Determining if the importance of different cues for premating isolation varies among cichlid radiations will be highly interesting and might be important to study the genetic basis of premating isolation. However, comprehensive analyses elucidating this issue will depend on further manipulation studies, disentangling the relative importance of different sensory modalities as conducted by Blais et al. [52] and Selz et al. [54] and including more species complexes.

    (c) Postmating prezygotic isolation

    Postmating prezygotic isolation (or gametic isolation) comprises all reproductive barriers acting between spawning and fertilization and is comparatively less studied in most organisms than premating isolation mechanisms [46], with cichlids being no exception. This is particularly unfortunate because the rapid evolution of reproductive proteins, mediating different processes after copulation including fertilization [119], has been attributed an important role in speciation in the last decades [46,120,121]. Gametic proteins evolve as by-product of sexual selection and an arms race between male and female proteins could result in particularly fast establishment of gametic isolation [46]. One particular good study, investigating interspecific fertilization rates in 26 heterospecific pairs of African cichlids, showed that fertilization failure linearly increases with divergence time, with complete fertilization failure observed after around 4 Myr (divergence time estimate depends on molecular clock used and can have wide confidence intervals; here, divergence time was calculated from a linear, internally calibrated clock using only recent biogeographic events, [12]). One of the few other examples testing specifically for postmating isolation between the closely related Lake Malawi species M. zebra and M. benetos, however, found no differences between interspecific fertilization rates and conspecific controls [122]. This indicates that gametic isolation in cichlids might not be as common and accumulating not as rapidly as proposed for other taxa [121,123,124].

    4. Postzygotic isolation

    (a) Intrinsic postzygotic isolation

    Intrinsic postzygotic isolation constitutes one of the most effective reproductive barriers owing to its irreversibility [125]. It can affect various life stages and is marked by flawed hybrid development that, at worst, results in hybrid inviability or physiological or behavioural hybrid sterility [46]. While the relative contribution of intrinsic barriers to reproductive isolation differs among species, there is agreement upon the fact that such barriers accumulate exponentially with genetic distance [12,46,126,127]. It has been implied that intrinsic incompatibilities are unlikely to evolve in the face of gene flow [42], however, this view has been challenged recently (e.g. [4345]). In cichlids, intrinsic barriers have been suggested to be neglectable among closely related species (e.g. [128,129]) and only cause complete hybrid inviability after around 4 Myr of divergence time (divergence time estimate depends on molecular clock used and can have wide confidence intervals; here, divergence time was calculated from a linear, internally calibrated clock using only recent biogeographic events; [12]). Considering that many cichlid radiations are much younger [21], it does not come as a surprise that only a minor role in conferring reproductive isolation in cichlids has been attributed to intrinsic barriers and some studies even consider their contribution to maintain species boundaries unlikely [49,128,130].

    We reviewed the literature to evaluate signatures of intrinsic incompatibilities in the form of inviability and sex ratio distortion among F1 hybrids of interspecific cichlid crosses (electronic supplementary material, table S3) with respect to genetic distance. D-loop sequences obtained from NCBI GenBank were used to compute pairwise genetic distances for interspecific crosses. The D-loop region is commonly sequenced and was, therefore, available for the majority of included species (for species without available D-loop sequences in GenBank, we used sequences of closely related species, electronic supplementary material, table S3). While we acknowledge that D-loop sequences might underestimate the actual level of genetic differentiation compared to more accurate estimates provided by larger datasets [131,132], those were not readily available for all the species used in this analysis. Intrinsic hybrid incompatibilities were analysed using a generalized linear model with binomial error distribution in R considering F1 viability or sex ratio distortion as response variable and genetic distance as explanatory variable (for more details on the methods, see the electronic supplementary material). While Stelkens et al. [12] provide a quantitative measure of hybrid inviability, the vast majority of studies only scored hybrids as viable or inviable. The same was done in studies reporting on sex ratio distortion. Therefore, we only consider binary data. While we acknowledge that this might not quantify the absolute strength of intrinsic incompatibilities, it provides comprehensive insights into the relative timing of the accumulation of hybrid inviability and hybrid sex ratio distortion, contributing to overall intrinsic postzygotic isolation.

    Based on our analysis, hybrid inviability is typically not observed among closely related species (figure 2a). These results are in line with previous studies [12], suggesting that complete hybrid inviability is reached only among more divergent cichlid species. The non-significant effect of genetic distance for F1 viability (estimate = 2.801, z = 0.581, p = 0.561) can probably be accounted for by the small number of studies reporting inviable hybrid crosses (electronic supplementary material, table S3). However, hybrid inviability appears not to be a necessary consequence of divergence in cichlids, because viable F1 offspring could be observed even in crosses with substantial genetic distance (figure 2a).

    Figure 2.

    Figure 2. Intrinsic incompatibilities in cichlids have been evaluated using (a) F1 hybrid inviability (n = 52 interspecific comparisons) and (b) F1 hybrid sex ratio distortion (n = 39 interspecific comparisons) as a function of pairwise genetic distances computed using D-loop sequences. Individual data points represent single interspecific crosses. Inviability of crosses and reports of sex ratio distortion were scored 1, while viable crosses and no reported skew in sex ratio were scored 0. Logistic regressions are depicted as black lines for F1 hybrid inviability (a; estimate = 2.801, z = 0.581, p = 0.561) and F1 sex ratio distortion (b; estimate = 16.330, z = 2.764, p = 0.006).

    Sex ratio distortion was also frequently detected in interspecific cichlid crosses (figure 2b; electronic supplementary material, table S3). Although its frequency increases significantly with genetic distance (estimate = 16.330, z = 2.764, p = 0.006), it can also be observed among closely related species (figure 2b), indicating that less severe intrinsic incompatibilities are manifested earlier than assumed. Skewed sex distributions are intriguing as they potentially point towards Haldane's rule, which states ‘when in offspring of two animal races one sex is absent, rare or sterile, that sex is the heterozygous sex’ [133, p. 1]. As the incompatibility described by this rule depicts an early stage in the evolution of postzygotic reproductive isolation, it is especially meaningful for young lineages where signs of intrinsic isolation seem to follow Haldane's rule in many species [46]. Because sex determination in cichlids is variable and complex [134], adherence to Haldane's rule as an intrinsic barrier remains only speculative [128,135]. Owing to the scant information on hybrid viability and other factors like hybrid sex ratio, our conclusions are only tentative. Yet, they might encourage future studies to report fitness consequences for hybrid crosses like growth rate, sex ratio and mortality, as this information will be required to improve our understanding of the importance of intrinsic isolation in cichlids.

    Whereas intrinsic incompatibilities between the parental genomes may be shielded in F1 hybrids, incompatible alleles causing severe developmental problems are often unmasked in the second generation by hybrid breakdown [136138]. Only few studies in cichlids have focused on the fitness consequences of these postzygotic barriers, i.e. DMIs [139,140]. Second-generation hybrids of seven interspecific crosses of African haplochromine cichlids (mostly comprising crosses among sympatric species but also some among allopatric ones, for details, see [139]) showed significantly lower fitness compared to the F1 hybrids and the pure grandparental lineages, and F2 inviability increased with divergence time [139]. Interestingly, segregation distortion in F2 hybrids between closely related sympatric species of Nicaraguan Midas cichlids indicate that hybrid incompatibilities in some cases might emerge earlier than assumed [140].

    Extreme cases of intrinsic incompatibilities are posed by major genomic alterations, such as changes in karyotypes or chromosomal rearrangements. Such forms of intrinsic incompatibilities have been reported in multiple different cichlid lineages (reviewed in [141]) and probably played a role in allopatric (e.g. Apistogramma spp., [142]), as well as in sympatric settings (e.g. Laetacara cf. dorsigera, [143]). While this form of intrinsic incompatibilities can generally cause strong and irreversible reproductive isolation, it is not clear if these incompatibilities were instrumental in the speciation process of the respective lineages or whether they only evolved after strong reproductive isolation was already established. Clearly, more comparative studies investigating such major genomic rearrangements (e.g. [144]) are needed to determine the generality of those findings and their contribution to reproductive isolation in different cichlid radiations. Taken together, all these studies stress that even though intrinsic barriers require divergence time to confer strong reproductive isolation, they could still play a role in young lineages, such as most cichlid radiations, and should not be excluded based solely on evolutionary youth.

    (b) Extrinsic postzygotic isolation

    Extrinsic ecological reproductive isolation, often referred to as ecological inviability, is characterized by decreased performance of hybrids in the respective environment of the parental species [46]. Studies investigating hybrid performance in the wild are rare (but see [145147]). Nonetheless, extrinsic isolation has been attributed an important role in speciation owing to the potential to diminish the number of hybrids, and Coyne & Orr [46, p. 255] even argued that ‘it has become fashionable to suggest that extrinsic, and especially ecological, postzygotic isolation is more common or more important than intrinsic in nature. This might well be true. But at present, such assertations rest more on intuition than data’. This lack of knowledge is partially owing to the difficulty of obtaining large-scale ecological information required to correctly assess the effect of extrinsic postzygotic isolation [46]. Unfortunately, more than 15 years later, this still reflects our rather poor understanding of extrinsic barriers compared to our knowledge of premating barriers in cichlids.

    Ecologically relevant phenotypes are frequently found to be intermediate in cichlid hybrids compared to the parental species [15,16]. Hence, it is commonly assumed that hybrids might have poor ecological performance and suffer from reduced fitness. Given that many cichlid species are highly specialized, it is possible that hybrids with intermediate phenotypes are often not as fit in nature as the parental phenotypes and some studies support this assumption [129,148150]. The importance of extrinsic ecological isolation rests on the strong assumption that there are no free niches and no ecological opportunities for hybrids with intermediate phenotypes in nature. In case those niches do exist, they might be occupied by other species that potentially outcompete hybrids [46,151], as the ecological space inhabited by cichlids tends to be heavily packed [47,152].

    Besides ecological inviability, extrinsic isolation also comprises behavioural hybrid sterility, as intermediate states of phenotypes important for mate choice can lower attractiveness and render it difficult for hybrids to find mates [46]. The importance of this barrier has been demonstrated for various taxa (e.g. [153,154]) and it probably also contributes to reproductive isolation in cichlids as hybrids often exhibit intermediate phenotypes in traits relevant for mate choice [23,74,75,148,155]. In most cases, females of the parental species show a clear preference for conspecific males over heterospecifics and hybrids [75,155]. Taken together, these studies suggest reduced mating success of hybrids consistent with behavioural sterility.

    5. Hybridization

    Traditionally, hybridization has been considered to disrupt the process of speciation [156158]. Owing to incomplete reproductive isolation, species can hybridize, which potentially causes divergent genetic clusters to collapse, negatively affecting species richness [18]. However, recent studies of East African cichlids have challenged this view by providing evidence that hybridization is not necessarily a destructive force to cichlid species diversity but actually fuelled explosive speciation bursts at early stages of these adaptive radiations [21,25,86]. Hybridization has the potential to drastically increase genetic and phenotypic diversity and it has been shown that entire species flocks are of hybrid origin for Lake Tanganyika [21], Lake Mweru [25] and the entire Lake Victoria region [86]. Ongoing hybridization is still frequently observed among cichlids and can create new phenotypes with extreme trait values [22,159]. Given ecological opportunity, such novelty can give hybrids the potential to outperform parental species outside of their respective niches [149] and if paired with non-random mating preferences, as demonstrated for some interspecific hybrids [75], this might set the stage for potential hybrid speciation. Therefore, strong but leaky reproductive isolation allowing for rare hybridization events in the presence of ecological opportunity has been demonstrated to be an important source of genetic variation and can catalyse radiation.

    6. Conclusion

    The importance of cichlid fishes as an evolutionary model system can be, among others, attributed to their incredible species richness and explosive diversification rates [1,2,4]. While cichlid diversification has been extensively studied (e.g. [1,3,4]), much less is known about what causes cessation of gene flow between diverging populations. However, in order to understand the origin and maintenance of the magnificent cichlid radiations, we urgently need to improve our knowledge concerning the barriers that confer reproductive isolation.

    Based on theory, premating isolation has been predicted to contribute most to overall reproductive isolation, as it acts early in the sequence of reproductive barriers and, therefore limits gene flow substantially [46,160,161]. Accordingly, premating isolation plays an important role in cichlids, as suggested by strong evidence at the meta-analysis level for assortative mating (figure 1a), which is in line with previously published work (e.g. [12,47,48]). However, its strength is variable. While fishes within certain cichlid lineages mate strongly assortatively, this does not necessarily apply for other lineages (figure 1b). Variation in the degree of premating isolation is probably influenced by a combination of genetic divergence and the geographical context, as expected under reinforcement (figure 1c). Studies investigating premating isolation are currently restricted to relatively few species complexes (i.e. Maylandia cf. zebra, Pundamilia cf. nyererei, Tropheus cf. moori, Amphilophus cf. citrinellus and Apistogramma spp.). Considering this lack of phylogenetic coverage and given the great diversity of the cichlid clade (e.g. reviewed in [1]), it is not surprising that strength and mechanisms of premating isolation differ across lineages and no general pattern has emerged, yet. Moreover, premating isolation in cichlids seems to be context dependent (figure 1d). While premating barriers can be very effective in limiting gene flow under the conditions they evolved in, they might not be highly resilient when environments are subject to drastic changes, as it is demonstrated for some lineages where species diversity is reduced after environmental change [18,162]. This indicates that strong reproductive isolation in cichlids most likely depends on the combination of multiple reproductive barriers. Our meta-analysis suggests that mild forms of intrinsic incompatibilities, such as sex ratio distortion, can already be observed among closely related species (figure 2b) and potentially supplement premating barriers. Moreover, a contribution of extrinsic postzygotic isolation is likely [129,148,149]. However, the importance of these different barriers for conferring reproductive isolation can vary along the speciation continuum [163,164].

    The speciation continuum has been suggested to consist of different stages based on divergence and strength of reproductive isolation among populations [44,165167]. Adaptive radiations are very informative to elucidate the signatures of these different stages because lineages within radiations may differ in the strength of reproductive isolation, as observed in stickleback fishes (Gasterosteidae) and Heliconius butterflies (Nymphalidae) [166,167]. Our meta-analyses suggest that members of different radiations of cichlid fishes also seem to vary in their stages along the speciation continuum (figures 1 and 2). Some cichlid lineages, such as Tropheus cf. moori, appear to be still at an early stage, where populations have undergone differentiation in some phenotypic traits (e.g. coloration), but levels of reproductive isolation are still low and genetic divergence might not be sufficient for reinforcement to act upon secondary contact ([58]; figure 1b). By contrast, other lineages, for example Pundamilia cf. nyererei, represent more advanced stages along the speciation continuum, marked by strong reproductive isolation (figure 1b). However, speciation is still incomplete, as changes in environmental conditions can disrupt species assortative mating [18]. As proposed for late stages in sticklebacks [167], completion of the speciation process in cichlids might require the establishment of major intrinsic incompatibilities, such as changes in karyotypes or chromosomal rearrangements, which can already be observed in some cichlid lineages [141144]. Because of great variation in the strength of reproductive isolation among different lineages, cichlid fishes might provide powerful contrasts to elucidate further characteristics distinguishing different stages along the speciation continuum associated with different biogeographic settings. However, this will require an increase in the number of comparative studies on different aspects of reproductive isolation across a broad taxonomic range of cichlid lineages, especially among radiations that differ in their evolutionary time.

    We suggest two lines of research that will substantially further our understanding of the barriers that might be essential to identify stages along the speciation continuum in cichlid fishes. Firstly, it will be crucial to improve the knowledge on linkage of genomic loci important for reproductive isolation. Early stages of differentiation are usually marked by signatures of divergent selection on few genomic loci and their surrounding regions (e.g. divergence hitchhiking, [165,168]). During later stages along the speciation continuum, linkage among multiple such loci under divergent selection might lead to an overall reduction in gene flow throughout the genome that results in genome-wide differentiation (e.g. genome hitchhiking, [44]). Given the increase in genomic resources available for cichlid fishes in the last decade (e.g. [169]), a necessary next step will be to look for linkage among genes known to contribute to reproductive isolation. Interesting candidates could be genes affecting visual sensitivity, female preference and nuptial coloration (see §3b(i), Visual cues: coloration and pattern). Knowledge about the genomic regions controlling visual sensitivity (e.g. [170]) and coloration (e.g. [15]) in cichlids is emerging. However, much less is known about the genetic basis of female preferences for nuptial coloration, although experimental studies suggest that it indeed has a genetic component [75,171]. Progress in determining the genomic regions affecting preference, coloration and visual sensitivity, as well as exploring the link among them will have a great impact on our understanding of premating isolation. Secondly, at the moment literature on reproductive isolation in cichlid fishes is highly skewed towards premating isolation, while the conclusions about the commonality and strength of postzygotic barriers rest on few studies. Therefore, it will be important to increase research directed towards determining the contribution of postzygotic reproductive isolation. It will be especially important to understand the role of divergence time and geographical settings for the evolution of intrinsic postzygotic barriers. The Lake Tanganyika radiations might be key in addressing the question of how genetic distance relates to the evolution of intrinsic incompatibilities. Lake Tanganyika is not only the oldest of the East African Rift Lakes, but its cichlid species flock is also the most diverse in ecological, morphological and genetic terms [172]. Different Tanganyikan lineages, like Ectodini, Lamprologini and Tropheini, have radiated to different degrees and vary significantly in their divergence time [21]. Thus, Lake Tanganyika offers the possibility for comparisons among radiations within the same lake, that might be essential to identify DMIs. Further, it also allows for meaningful contrasts to determine the strength of reproductive isolation during different stages along the speciation continuum. Following these aforementioned lines of research will not only improve our knowledge of the barriers that triggered and maintain the current cichlid diversity, but can also contribute to the understanding of the evolutionary dynamics that characterize different stages along the speciation continuum.

    Data accessibility

    All the data used in this Review are included in the electronic supplementary material

    Authors' contributions

    S.J.R., J.T.-D. and A.M. conceptualized the project and S.J.R. and J.T.-D. collected and analysed the data. S.J.R. wrote the manuscript with revisions from all authors.

    Competing interests

    We declare we have no competing interests.

    Funding

    This work was supported by the University of Konstanz. S.J.R. was funded by a fellowship from the Hector Fellow Academy and J.T.-D. was supported by grants of the DFG (grant no. TO914/3-1) and (grant no. ME1725/22-1) and the Zukunftskolleg of the University of Konstanz.

    Acknowledgements

    We thank Andreas F. Kautt for fruitful discussion and valuable comments.

    Footnotes

    One contribution of 19 to a theme issue ‘Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers’.

    Electronic supplementary material is available online at https://doi.org/10.6084/m9.figshare.c.5007992.

    Published by the Royal Society. All rights reserved.

    References