Ancient and Recent Hybridization in the Oreochromis Cichlid Fishes.

Cichlid fishes of the genus Oreochromis (tilapia) are among the most important fish for inland capture fisheries and global aquaculture. Deliberate introductions of non-native species for fisheries improvement and accidental escapees from farms have resulted in admixture with indigenous species. Such hybridization may be detrimental to native biodiversity, potentially leading to genomic homogenization of populations and the loss of important genetic material associated with local adaptation. By contrast, introgression may fuel diversification when combined with ecological opportunity, by supplying novel genetic combinations. To date, the role of introgression in the evolutionary history of tilapia has not been explored. Here we studied both ancient and recent hybridization in tilapia, using whole genome resequencing of 575 individuals from 23 species. We focused on Tanzania, a natural hotspot of tilapia diversity, and a country where hybridization between exotic and native species in the natural environment has been previously reported. We reconstruct the first genome-scale phylogeny of the genus and reveal prevalent ancient gene flow across the Oreochromis phylogeny. This has likely resulted in the hybrid speciation of one species, O. chungruruensis. We identify multiple cases of recent hybridization between native and introduced species in the wild, linked to the use of non-native species in both capture fisheries improvement and aquaculture. This has potential implications for both conservation of wild populations and the development of the global tilapia aquaculture industry.

Sympatry and parapatry among rocky reef cichlids of Lake Victoria explained by female mating preferences.

Work on the Lake Victoria cichlids Pundamilia nyererei (red dorsum males, deeper water), Pundamilia pundamilia (blue males, shallower water) and related species pairs has provided insights into processes of speciation. Here, we investigate the female mating behaviour of 5 Pundamilia species and 4 of their F1 hybrids through mate choice trials and paternity testing. Complete assortative mating was observed among all sympatric species. Parapatric species with similar depth habitat distributions interbred whereas other parapatric and allopatric species showed complete assortative mating. F1 hybrids mated exclusively with species accepted by females of the parental species. The existence of complete assortative mating among some currently allopatric species suggests that pre-existing mating barriers could be sufficient to explain current patterns of co-existence, although, of course, many other factors may be involved. Regardless of the mechanism, mating preferences may influence species distribution in potentially hybridizing taxa, such as in the adaptive radiation of cichlid fish. We suggest that this at least partly explains why some species fail to establish breeding populations in locations where they are occasionally recorded. Our results support the notion that the mating preferences of potentially cross-breeding species ought to be included in coexistence theory.

Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes.

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

Whole genome resequencing data enables a targeted SNP panel for conservation and aquaculture of Oreochromis cichlid fishes.

Cichlid fish of the genus Oreochromis form the basis of the global tilapia aquaculture and fisheries industries. Broodstocks for aquaculture are often collected from wild populations, which in Africa may be from locations containing multiple Oreochromis species. However, many species are difficult to distinguish morphologically, hampering efforts to maintain good quality farmed strains. Additionally, non-native farmed tilapia populations are known to be widely distributed across Africa and to hybridize with native Oreochromis species, which themselves are important for capture fisheries. The morphological identification of these hybrids is particularly unreliable. Here, we describe the development of a single nucleotide polymorphism (SNP) genotyping panel from whole-genome resequencing data that enables targeted species identification in Tanzania. We demonstrate that an optimized panel of 96 genome-wide SNPs based on FST outliers performs comparably to whole genome resequencing in distinguishing species and identifying hybrids. We also show this panel outperforms microsatellite-based and phenotype-based classification methods. Case studies indicate several locations where introduced aquaculture species have become established in the wild, threatening native Oreochromis species. The novel SNP markers identified here represent an important resource for assessing broodstock purity in hatcheries and helping to conserve unique endemic biodiversity.

A new species of deep-water Lethrinops (Cichlidae) from Lake Malawi.

A new species of cichlid fish, Lethrinops atrilabris is described from specimens collected by trawling at a depth of around 90 m off Monkey Bay, southern Lake Malawi. It is assigned to the genus Lethrinops on the basis of its vertical flank barring, lack of enlarged cephalic lateral line canal pores and the form of the lower jaw dental arcade. It can be distinguished from congeneric species by its male breeding dress of contrasting flank barring and dark ventral surface, most strikingly on the lips, throat and chest, its relatively small known maximum size [<75 mm standard length (SL)], large eyes (38%-41% head length), laterally compressed body (depth 2.5-2.7 times max head width) and lower gill raker count (13-14).

Ancestral Hybridization Facilitated Species Diversification in the Lake Malawi Cichlid Fish Adaptive Radiation.

The adaptive radiation of cichlid fishes in East African Lake Malawi encompasses over 500 species that are believed to have evolved within the last 800,000 years from a common founder population. It has been proposed that hybridization between ancestral lineages can provide the genetic raw material to fuel such exceptionally high diversification rates, and evidence for this has recently been presented for the Lake Victoria region cichlid superflock. Here, we report that Lake Malawi cichlid genomes also show evidence of hybridization between two lineages that split 3-4 Ma, today represented by Lake Victoria cichlids and the riverine Astatotilapia sp. "ruaha blue." The two ancestries in Malawi cichlid genomes are present in large blocks of several kilobases, but there is little variation in this pattern between Malawi cichlid species, suggesting that the large-scale mosaic structure of the genomes was largely established prior to the radiation. Nevertheless, tens of thousands of polymorphic variants apparently derived from the hybridization are interspersed in the genomes. These loci show a striking excess of differentiation across ecological subgroups in the Lake Malawi cichlid assemblage, and parental alleles sort differentially into benthic and pelagic Malawi cichlid lineages, consistent with strong differential selection on these loci during species divergence. Furthermore, these loci are enriched for genes involved in immune response and vision, including opsin genes previously identified as important for speciation. Our results reinforce the role of ancestral hybridization in explosive diversification by demonstrating its significance in one of the largest recent vertebrate adaptive radiations.

Newly discovered cichlid fish biodiversity threatened by hybridization with non-native species.

Invasive freshwater fishes are known to readily hybridize with indigenous congeneric species, driving loss of unique and irreplaceable genetic resources. Here we reveal that newly discovered (2013-2016) evolutionarily significant populations of Korogwe tilapia (Oreochromis korogwe) from southern Tanzania are threatened by hybridization with the larger invasive Nile tilapia (Oreochromis niloticus). We use a combination of morphology, microsatellite allele frequencies and whole genome sequences to show that O. korogwe from southern lakes (Nambawala, Rutamba and Mitupa) are distinct from geographically disjunct populations in northern Tanzania (Zigi River and Mlingano Dam). We also provide genetic evidence of O. korogwe × niloticus hybrids in three southern lakes and demonstrate heterogeneity in the extent of admixture across the genome. Finally, using the least admixed genomic regions we estimate that the northern and southern O. korogwe populations most plausibly diverged ~140,000 years ago, suggesting that the geographical separation of the northern and southern groups is not a result of a recent translocation, and instead these populations represent independent evolutionarily significant units. We conclude that these newly discovered and phenotypically unique cichlid populations are already threatened by hybridization with an invasive species, and propose that these irreplaceable genetic resources would benefit from conservation interventions.

Rare morph Lake Malawi mbuna cichlids benefit from reduced aggression from con- and hetero-specifics.

Balancing selection is important for the maintenance of polymorphism as it can prevent either fixation of one morph through directional selection or genetic drift, or speciation by disruptive selection. Polychromatism, the presence of multiple genetically determined colour phenotypes, can be maintained if the fitness of alternative morphs depends on the relative frequency in a population. In aggressive species, negative frequency-dependent antagonism can prevent an increase in the frequency of rare morphs as they would only benefit from increased fitness while they are rare. Heterospecific aggression is common in nature and has the potential to contribute to rare morph advantage. Here we carry out field observations and laboratory aggression experiments with mbuna cichlids from Lake Malawi, to investigate the role of con- and heterospecific aggression in the maintenance of polychromatism and identify benefits to rare morphs which are likely to result from reduced aggression. We hypothesize that rare morph individuals receive less aggression than common morph individuals and therefore have an ecological advantage. Within species we found that males and females bias aggression towards their own morph, adding to the evidence that inherent own-morph aggression biases can contribute to balancing selection. Over-representation of rare morph territory owners may be influenced by two factors; higher tolerance of different morph individuals as neighbours, and the ability of rare morphs to spend more time feeding. Reduced aggression to rare morph individuals by heterospecifics may also contribute to rare morph advantage.

Molecular phylogeny of Oreochromis (Cichlidae: Oreochromini) reveals mito-nuclear discordance and multiple colonisation of adverse aquatic environments.

Although the majority of cichlid diversity occurs in the African Great Lakes, these fish have also diversified across the African continent. Such continental radiations, occurring in both rivers and lakes have received far less attention than lacustrine radiations despite some members, such as the oreochromine cichlids (commonly referred to as 'tilapia'), having significant scientific and socio-economic importance both within and beyond their native range. Unique among cichlids, several species of the genus Oreochromis exhibit adaptation to soda conditions (including tolerance to elevated temperatures and salinity), which are of interest from evolutionary biology research and aquaculture perspectives. Questions remain regarding the factors facilitating the diversification of this group, which to date have not been addressed within a phylogenetic framework. Here we present the first comprehensive (32/37 described species) multi-marker molecular phylogeny of Oreochromis and closely related Alcolapia, based on mitochondrial (1583 bp) and nuclear (3092 bp) sequence data. We show widespread discordance between nuclear DNA and mitochondrial DNA trees. This could be the result of incomplete lineage sorting and/or introgression in mitochondrial loci, although we did not find a strong signal for the latter. Based on our nuclear phylogeny we demonstrate that adaptation to adverse conditions (elevated salinity, temperature, or alkalinity) has occurred multiple times within Oreochromis, but that adaptation to extreme (soda) conditions (high salinity, temperature, and alkalinity) has likely arisen once in the lineage leading to O. amphimelas and Alcolapia. We also show Alcolapia is nested within Oreochromis, which is in agreement with previous studies, and here revise the taxonomy to synonymise the genus in Oreochromis, retaining the designation as subgenus Oreochromis (Alcolapia).

Losing cichlid fish biodiversity: genetic and morphological homogenization of tilapia following colonization by introduced species.

Among the many negative impacts of invasive species, hybridization with indigenous species has increasingly become recognized as a major issue. However, relatively few studies have characterized the phenotypic outcomes of hybridization following biological invasions. Here we investigate the genetic and morphological consequences of stocking invasive tilapia species in two water bodies in central Tanzania. We sampled individuals from the Mindu Reservoir on the Ruvu river system, and at Kidatu on the Great Ruaha-Rufiji river system. We screened individuals at 16 microsatellite loci, and quantified morphology using geometric morphometrics and linear measurements. In both the Mindu and Kidatu systems, we identified evidence of hybridization between indigenous Wami tilapia (Oreochromis urolepis) and the introduced Nile tilapia (Oreochromis niloticus) or blue-spotted tilapia (Oreochromis leucostictus). At both sites, purebred individuals could largely be separated using geometric morphometric variables, with hybrids occupying a broad morphospace among the parental species. Our data demonstrate that the gene pools and phenotypic identity of the indigenous O. urolepis have been severely impacted by the stocking of the invasive species. Given the lack of evidence for clear commercial benefits from stocking invasive tilapia species in waters already populated by indigenous congenerics, we suggest further spread of introduced species should be undertaken with considerable caution.

The genetics of mate preferences in hybrids between two young and sympatric Lake Victoria cichlid species.

The genetic architecture of mate preferences is likely to affect significant evolutionary processes, including speciation and hybridization. Here, we investigate laboratory hybrids between a pair of sympatric Lake Victoria cichlid fish species that appear to have recently evolved from a hybrid population between similar predecessor species. The species demonstrate strong assortative mating in the laboratory, associated with divergent male breeding coloration (red dorsum versus blue). We show in a common garden experiment, using DNA-based paternity testing, that the strong female mate preferences among males of the two species are fully recovered in a large fraction of their F2 hybrid generation. Individual hybrid females often demonstrated consistent preferences in multiple mate choice trials (more than or equal to five) across a year or more. This result suggests that female mate preference is influenced by relatively few major genes or genomic regions. These preferences were not changed by experience of a successful spawning event with a male of the non-preferred species in a no-choice single-male trial. We found no evidence for imprinting in the F2 hybrids, although the F1 hybrid females may have been imprinted on their mothers. We discuss this nearly Mendelian inheritance of consistent innate mate preferences in the context of speciation theory.

Geographical ancestry of Lake Malawi's cichlid fish diversity.

The Lake Malawi haplochromine cichlid flock is one of the largest vertebrate adaptive radiations. The geographical source of the radiation has been assumed to be rivers to the south and east of Lake Malawi, where extant representatives of the flock are now present. Here, we provide mitochondrial DNA evidence suggesting the sister taxon to the Lake Malawi radiation is within the Great Ruaha river in Tanzania, north of Lake Malawi. Estimates of the time of divergence between the Lake Malawi flock and this riverine sister taxon range from 2.13 to 6.76 Ma, prior to origins of the current radiation 1.20-4.06 Ma. These results are congruent with evaluations of 2-3.75 Ma fossil material that suggest past faunal connections between Lake Malawi and the Ruaha. We propose that ancestors of the Malawi radiation became isolated within the catchment during Pliocene rifting that formed both Lake Malawi and the Kipengere/Livingstone mountain range, before colonizing rivers to the south and east of the lake region and radiating within the lake basin. Identification of this sister taxon allows tests of whether standing genetic diversity has predisposed Lake Malawi cichlids to rapid speciation and adaptive radiation.

Nuclear environmental DNA resolves fine-scale population genetic structure in an aquatic habitat.

There is considerable potential for nuclear genomic material in environmental DNA (eDNA) to inform us of population genetic structure within aquatic species. We tested if nuclear allelic composition data sourced from eDNA can resolve fine scale spatial genetic structure of the cichlid fish Astatotilapia calliptera in Lake Masoko, Tanzania. In this ∼35 m deep crater lake the species is diverging into two genetically distinguishable ecomorphs, separated by a thermo-oxycline at ∼15 m that divides biologically distinct water masses. We quantified population genetic structure along a depth transect using single nucleotide polymorphisms (SNPs) derived from genome sequencing of 530 individuals. This population genetic structure was reflected in a focal set of SNPs that were also reliably amplified from eDNA - with allele frequencies derived from eDNA reflecting those of fish within each depth zone. Thus, by targeting known genetic variation between populations within aquatic eDNA, we measured genetic structure within the focal species.

Lateral line system diversification during the early stages of ecological speciation in cichlid fish.

BACKGROUND: The mechanosensory lateral line system is an important sensory modality in fishes, informing multiple behaviours related to survival including finding food and navigating in dark environments. Given its ecological importance, we may expect lateral line morphology to be under disruptive selection early in the ecological speciation process. Here we quantify the lateral line system morphology of two ecomorphs of the cichlid fish Astatotilapia calliptera in crater Lake Masoko that have diverged from common ancestry within the past 1,000 years. RESULTS: Based on geometric morphometric analyses of CT scans, we show that the zooplanktivorous benthic ecomorph that dominates the deeper waters of the lake has large cranial lateral line canal pores, relative to those of the nearshore invertebrate-feeding littoral ecomorph found in the shallower waters. In contrast, fluorescence imaging revealed no evidence for divergence between ecomorphs in the number of either superficial or canal neuromasts. We illustrate the magnitude of the variation we observe in Lake Masoko A. calliptera in the context of the neighbouring Lake Malawi mega-radiation that comprises over 700 species. CONCLUSIONS: These results provide the first evidence of divergence in this often-overlooked sensory modality in the early stages of ecological speciation, suggesting that it may have a role in the broader adaptive radiation process.

Ancient hybridization and phenotypic novelty within Lake Malawi's cichlid fish radiation.

Does hybridization play a broad innovative role in evolution? Many studies have shown hybrid origins of individual species, particularly in major adaptive radiations, but this may be a consequence, rather than a cause, of the existence of many closely related species. Cases of hybridization in the early stages of major adaptive radiations are comparatively rare. Here, we report phylogenetic evidence for ancient introgression between distinct lineages of the species-rich Lake Malawi haplochromine cichlid fishes. Mitochondrial DNA (mtDNA) sequences indicated surprisingly close relationships between the shallow-water rocky habitat "Mbuna" species and a group of dark-adapted "Deep-Benthic" species specialized for feeding in low-light conditions (dawn/dusk, under overhangs, and deep water). By contrast, analyses of nuclear amplified fragment length polymorphism data demonstrated that these Deep-Benthic cichlids were more closely related to shallow water "Shallow-Benthic" soft-sediment feeders, a group that shares similar head and body morphology. A coalescent-based computer simulation indicated that the mtDNA similarity of rocky habitat Mbuna species and dark-adapted Deep-Benthic species was due to hybridization rather than incomplete lineage sorting. Comparisons of morphology indicated that some Deep-Benthic species possessed novel morphology not present in other Lake Malawi species groups. Thus, these analyses support the hypothesis that ancient hybridization occurred within the Lake Malawi cichlid radiation, that the event occurred before the radiation of a species group adapted to low-light benthic habitats, and that this group went on to dominate the deep-water regions of Lake Malawi. The results of this study contribute to a growing literature consistent with a creative role of hybridization in the evolution of species diversity and adaptive radiations.

A new species of Lethrinops (Cichliformes: Cichlidae) from a Lake Malawi satellite lake, believed to be extinct in the wild.

A new species of cichlid fish, Lethrinops chilingali is described from specimens collected from Lake Chilingali, near Nkhotakota, Malawi. It is assigned to the genus Lethrinops based on the form of the lower jaw dental arcade and by the absence of traits diagnostic of the phenotypically similar Ctenopharynx, Taeniolethrinops and Tramitichromis. It also lacks the enlarged cephalic lateral line canal pores found in species of Alticorpus and Aulonocara. The presence of a broken horizontal stripe on the flanks of females and immature/non-territorial males of Lethrinops chilingali distinguishes them from all congeners, including Lethrinops lethrinus, in which the stripe is typically continuous. Lethrinops chilingali also has a relatively shorter snout, shorter lachrymal bone and less ventrally positioned mouth than Lethrinops lethrinus. It appears likely that Lethrinops chilingali is now extinct in the wild, as this narrow endemic species has not been positively recorded in the natural environment since 2009. Breeding populations remain in captivity.

The influence of oceanographic fronts and early-life-history traits on connectivity among littoral fish species.

The spatial distribution of neutral genetic diversity is mainly influenced by barriers to dispersal. The nature of such barriers varies according to the dispersal means and capabilities of the organisms concerned. Although these barriers are often obvious on land, in the ocean they can be more difficult to identify. Determining the relative influence of physical and biotic factors on genetic connectivity remains a major challenge for marine ecologists. Here, we compare gene flow patterns of 7 littoral fish species from 6 families with a range of early-life-history traits sampled at the same geographic locations across common environmental discontinuities in the form of oceanic fronts in the Western Mediterranean. We show that these fronts represent major barriers to gene flow and have a strong influence on the population genetic structure of some fish species. We also found no significant relation between the early-life-history traits most commonly investigated (egg type, pelagic larval duration, and inshore-offshore spawning) and gene flow patterns, suggesting that other life-history factors should deserve attention. The fronts analyzed and the underlying physical mechanisms are not site-specific but common among the oceans, suggesting the generality of our findings.

Oca2 targeting using CRISPR/Cas9 in the Malawi cichlid Astatotilapia calliptera.

Identifying genetic loci underlying trait variation provides insights into the mechanisms of diversification, but demonstrating causality and characterizing the role of genetic loci requires testing candidate gene function, often in non-model species. Here we establish CRISPR/Cas9 editing in Astatotilapia calliptera, a generalist cichlid of the remarkably diverse Lake Malawi radiation. By targeting the gene oca2 required for melanin synthesis in other vertebrate species, we show efficient editing and germline transmission. Gene edits include indels in the coding region, probably a result of non-homologous end joining, and a large deletion in the 3' untranslated region due to homology-directed repair. We find that oca2 knock-out A. calliptera lack melanin, which may be useful for developmental imaging in embryos and studying colour pattern formation in adults. As A. calliptera resembles the presumed generalist ancestor of the Lake Malawi cichlids radiation, establishing genome editing in this species will facilitate investigating speciation, adaptation and trait diversification in this textbook radiation.

Epigenetic divergence during early stages of speciation in an African crater lake cichlid fish.

Epigenetic variation can alter transcription and promote phenotypic divergence between populations facing different environmental challenges. Here, we assess the epigenetic basis of diversification during the early stages of speciation. Specifically, we focus on the extent and functional relevance of DNA methylome divergence in the very young radiation of Astatotilapia calliptera in crater Lake Masoko, southern Tanzania. Our study focuses on two lake ecomorphs that diverged approximately 1,000 years ago and a population in the nearby river from which they separated approximately 10,000 years ago. The two lake ecomorphs show no fixed genetic differentiation, yet are characterized by different morphologies, depth preferences and diets. We report extensive genome-wide methylome divergence between the two lake ecomorphs, and between the lake and river populations, linked to key biological processes and associated with altered transcriptional activity of ecologically relevant genes. Such genes differing between lake ecomorphs include those involved in steroid metabolism, hemoglobin composition and erythropoiesis, consistent with their divergent habitat occupancy. Using a common-garden experiment, we found that global methylation profiles are often rapidly remodeled across generations but ecomorph-specific differences can be inherited. Collectively, our study suggests an epigenetic contribution to the early stages of vertebrate speciation.

An extant cichlid fish radiation emerged in an extinct Pleistocene lake.

The haplochromine cichlid fish of the East African Great Lakes represent some of the fastest and most species-rich adaptive radiations known, but rivers in most of Africa accommodate only a few morphologically similar species of haplochromine cichlid fish. This has been explained by the wealth of ecological opportunity in large lakes compared with rivers. It is therefore surprising that the rivers of southern Africa harbour many, ecologically diverse haplochromines. Here we present genetic, morphological and biogeographical evidence suggesting that these riverine cichlids are products of a recent adaptive radiation in a large lake that dried up in the Holocene. Haplochromine species richness peaks steeply in an area for which geological data reveal the historical existence of Lake palaeo-Makgadikgadi. The centre of this extinct lake is now a saltpan north of the Kalahari Desert, but it once hosted a rapidly evolving fish species radiation, comparable in morphological diversity to that in the extant African Great Lakes. Importantly, this lake seeded all major river systems of southern Africa with ecologically diverse cichlids. This discovery reveals how local evolutionary processes operating during a short window of ecological opportunity can have a major and lasting effect on biodiversity on a continental scale.