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.

Genetic analysis of Octopus cyanea reveals high gene flow in the South-West Indian Ocean.

Octopus cyanea (Gray, 1849), abundant in the South-West Indian Ocean (SWIO), constitutes a vital resource for both subsistence and commercial fisheries. However, despite this socioeconomic importance, and recent indications of overfishing, little is known about the population structure of O. cyanea in the region. To inform sustainable management strategies, this study assessed the spatio-temporal population structure and genetic variability of O. cyanea at 20 sites in the SWIO (Kenya, Tanzania, Mozambique, Madagascar, Mauritius, Rodrigues, and the Seychelle Islands) by complementary analysis of mitochondrial DNA (mtDNA) noncoding region (NCR) sequences and microsatellite markers. MtDNA analysis revealed a shallow phylogeny across the region, with demographic tests suggesting historic population fluctuations that could be linked to glacial cycles. Contrary to expectations, NCR variation was comparable to other mtDNA regions, indicating that the NCR is not a hypervariable region. Both nuclear and mtDNA marker types revealed a lack of genetic structure compatible with high gene flow throughout the region. As adults are sedentary, this gene flow likely reflects connectivity by paralarval dispersal. All samples reported heterozygote deficits, which, given the overall absence of structure, likely reflect ephemeral larval recruitment variability. Levels of mtDNA and nuclear variability were similar at all locations and congruent with those previously reported for harvested Octopodidae, implying resilience to genetic erosion by drift, providing current stock sizes are maintained. However, as O. cyanea stocks in the SWIO represent a single, highly connected population, fisheries may benefit from additional management measures, such as rotational closures aligned with paralarval ecology and spanning geopolitical boundaries.

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.

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.

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.

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.

Environmental DNA-based xenomonitoring for determining Schistosoma presence in tropical freshwaters.

BACKGROUND: Schistosomiasis is a neglected tropical disease that infects over 200 million people worldwide. Control measures can benefit from improved surveillance methods in freshwaters, with environmental DNA (eDNA) surveys having the potential to offer effective and rapid detection of schistosomes. However, sampling eDNA directly from natural water bodies can lead to inaccurate estimation of infection risk if schistosome eDNA is rare in the environment. Here we report a xenomonitoring method that allows schistosome infections of host snail species to be determined from eDNA in water used to house those snails. METHODS: Host snail species were collected and placed in containers of water and allowed to shed cercariae, and then water samples were filtered and tested using qPCR assays specific to the African species Schistosoma mansoni and Schistosoma haematobium. We evaluated this "eDNA-based xenomonitoring" approach by experimentally comparing the results to those obtained from direct qPCR screening of tissue sourced from the snails in the experiment. RESULTS: We found that our method accurately diagnosed the presence of S. mansoni-infected snails in all tests, and S. haematobium-infected snails in 92% of tests. Moreover, we found that the abundance of Schistosoma eDNA in experiments was directly dependent on the number and biomass of infected snails. CONCLUSIONS: These results provide a strong indication that this surveillance method combining the utility of eDNA-based monitoring with the reliability of traditional xenomonitoring approaches could be used to accurately assay the presence of Schistosoma species in natural habitats. This approach may be well-suited for epidemiological studies and monitoring in endemic areas, where it can assist schistosomiasis control by indicating infection risk from freshwaters and guiding necessary interventions to eliminate the disease.

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.

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

Widespread colonisation of Tanzanian catchments by introduced Oreochromis tilapia fishes: the legacy from decades of deliberate introduction.

From the 1950s onwards, programmes to promote aquaculture and improve capture fisheries in East Africa have relied heavily on the promise held by introduced species. In Tanzania these introductions have been poorly documented. Here we report the findings of surveys of inland water bodies across Tanzania between 2011 and 2017 that clarify distributions of tilapiine cichlids of the genus Oreochromis. We identified Oreochromis from 123 sampling locations, including 14 taxa restricted to their native range and three species that have established populations beyond their native range. Of these three species, the only exotic species found was blue-spotted tilapia (Oreochromis leucostictus), while Nile tilapia (Oreochromis niloticus) and Singida tilapia (Oreochromis esculentus), which are both naturally found within the country of Tanzania, have been translocated beyond their native range. Using our records, we developed models of suitable habitat for the introduced species based on recent (1960-1990) and projected (2050, 2070) East African climate. These models indicated that presence of suitable habitat for these introduced species will persist and potentially expand across the region. The clarification of distributions provided here can help inform the monitoring and management of biodiversity, and inform policy related to the future role of introduced species in fisheries and aquaculture.

Limited hybridization between introduced and Critically Endangered indigenous tilapia fishes in northern Tanzania.

Hybridization between introduced and indigenous species can lead to loss of unique genetic resources and precipitate extinction. In Tanzania, the Nile tilapia (Oreochromis niloticus) and blue-spotted tilapia (Oreochromis leucostictus) have been widely introduced to non-native habitats for aquaculture and development of capture fisheries. Here, we aimed to quantify interspecific hybridization between these introduced species and the indigenous species Oreochromis esculentus, Oreochromis jipe and Oreochromis korogwe. In the Pangani basin, several hybrids were observed (O. niloticus × O. jipe, O. leucostictus × O. jipe, O. niloticus × O. korogwe), although hybrids were relatively uncommon within samples relative to purebreds. Hybrids between the native O. jipe × O. korogwe were also observed. In the Lake Victoria basin, no evidence of hybrids was found. Analysis of body shape using geometric morphometrics suggested that although purebreds could be discriminated from one another, hybrids could not be readily identified on body and head shape alone. These results provide the first evidence of hybridization between the introduced species and the Critically Endangered O. jipe in Tanzania. Given uncertainty regarding benefits of introduced species over large-bodied indigenous species in aquaculture and capture fisheries, we suggest that future introductions of hybridization-prone species should be carefully evaluated.

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.

Multiple colonisations of the Lake Malawi catchment by the genus Opsaridium (Teleostei: Cyprinidae).

It has been proposed that the fish faunas of African rivers assemble through multiple colonisation events, while lake faunas form additionally through intralacustine speciation. While this pattern has been established for many lineages, most notably cichlids, there are opportunities to further investigate the concept using phylogenies of congeneric endemic species within ancient lake catchments. The Lake Malawi catchment contains three river-spawning cyprinids of the genus Opsaridium, two of which are endemic. These species differ in body size, migratory behaviour and habitat use, but it has never previously been tested if these represent a monophyletic radiation, or have instead colonised the lake independently. We placed these species in a broader phylogeny of Opsaridium and the related genus Raiamas, including all known species from the river systems surrounding Lake Malawi. Our results suggest that each of the species has independently colonised the lake catchment, with all three taxa having well-defined sister taxa outside of the lake, and all sharing a common ancestor ∼14.9million years ago, before the Lake Malawi basin started to form ∼8.6million years ago. Additionally, the results strongly support previous observations that Opsaridium is not a monophyletic group, but instead contains Raiamas from the Congo drainage. Together these results are supportive of the concept that river fish faunas within African catchments are primarily assembled through a process of accumulation from independent origins, rather than within-catchment speciation and adaptive radiation. In light of these results we also suggest there is scope for a re-evaluation of systematics of both Opsaridium and Raiamas.

Genomic islands of speciation separate cichlid ecomorphs in an East African crater lake.

The genomic causes and effects of divergent ecological selection during speciation are still poorly understood. Here we report the discovery and detailed characterization of early-stage adaptive divergence of two cichlid fish ecomorphs in a small (700 meters in diameter) isolated crater lake in Tanzania. The ecomorphs differ in depth preference, male breeding color, body shape, diet, and trophic morphology. With whole-genome sequences of 146 fish, we identified 98 clearly demarcated genomic "islands" of high differentiation and demonstrated the association of genotypes across these islands with divergent mate preferences. The islands contain candidate adaptive genes enriched for functions in sensory perception (including rhodopsin and other twilight-vision-associated genes), hormone signaling, and morphogenesis. Our study suggests mechanisms and genomic regions that may play a role in the closely related mega-radiation of Lake Malawi.

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.

Single-male paternity in coelacanths.

Latimeria chalumnae, a 'living fossil,' is of great scientific interest, as it is closely related to the aquatic ancestors of land-living tetrapods. Latimeria show internal fertilization and bear live young, but their reproductive behaviour is poorly known. Here we present for the first time a paternity analysis of the only available material from gravid females and their offspring. We genotype two L. chalumnae females and their unborn brood for 14 microsatellite loci. We find that the embryos are closely related to each other and never show more than three different alleles per locus, providing evidence for a single father siring all of the offspring. We reconstruct the father's genotype but cannot identify it in the population. These data suggest that coelacanths have a monogamous mating system and that individual relatedness is not important for mate choice.

Coelacanth genomes reveal signatures for evolutionary transition from water to land.

Coelacanths are known as "living fossils," as they show remarkable morphological resemblance to the fossil record and belong to the most primitive lineage of living Sarcopterygii (lobe-finned fishes and tetrapods). Coelacanths may be key to elucidating the tempo and mode of evolution from fish to tetrapods. Here, we report the genome sequences of five coelacanths, including four Latimeria chalumnae individuals (three specimens from Tanzania and one from Comoros) and one L. menadoensis individual from Indonesia. These sequences cover two African breeding populations and two known extant coelacanth species. The genome is ∼2.74 Gbp and contains a high proportion (∼60%) of repetitive elements. The genetic diversity among the individuals was extremely low, suggesting a small population size and/or a slow rate of evolution. We found a substantial number of genes that encode olfactory and pheromone receptors with features characteristic of tetrapod receptors for the detection of airborne ligands. We also found that limb enhancers of bmp7 and gli3, both of which are essential for limb formation, are conserved between coelacanth and tetrapods, but not ray-finned fishes. We expect that some tetrapod-like genes may have existed early in the evolution of primitive Sarcopterygii and were later co-opted to adapt to terrestrial environments. These coelacanth genomes will provide a cornerstone for studies to elucidate how ancestral aquatic vertebrates evolved into terrestrial animals.

Extremely slow rate of evolution in the HOX cluster revealed by comparison between Tanzanian and Indonesian coelacanths.

Coelacanths are known as "living fossils" because their morphology has changed very little from that in the fossil record. To elucidate why coelacanths have evolved so slowly is thus of primary importance in evolutionary biology. In the present study, we determined the entire sequence of the HOX cluster of the Tanzanian coelacanth (Latimeria chalumnae) and compared it with that of the Indonesian coelacanth (L. menadoensis), which was available in the literature. The most intriguing result was the extremely small genetic divergence between the two coelacanths. The synonymous divergence of the HOX coding region between the two coelacanths was estimated to be 0.07%, which is ~11-fold smaller than that of human-chimp. When we applied the estimated divergence time of the two coelacanths of 6 million years ago (MYA) and 30 MYA, which were proposed in independent mitochondrial DNA analyses, the synonymous substitution rate of the coelacanth HOX cluster was estimated to be ~11-fold and 56-fold smaller than that of human-chimp, respectively. Thus, the present study implies that the reduction of the nucleotide substitution rate in coelacanth HOX genes may account for the conservation of coelacanth morphology during evolution.

Population divergence in East African coelacanths.

The coelacanth, Latimeria chalumnae, occurs at the Eastern coast of Africa from South Africa up to Kenya. It is often referred to as a living fossil mainly because of its nearly unchanged morphology since the Middle Devonian. As it is a close relative to the last common ancestor of fish and tetrapods, molecular studies mostly focussed on their phylogenetic relationships. We now present a population genetic study based on 71 adults from the whole known range of the species. Despite an overall low genetic diversity, there is evidence for divergence of local populations. We assume that originally the coelacanths at the East African Coast derived from the Comoros population, but have since then diversified into additional independent populations: one in South Africa and another in Tanzania. Unexpectedly, we find a split of the Comoran coelacanths into two sympatric subpopulations. Despite its undeniably slow evolutionary rate, the coelacanth still diversifies and is therefore able to adapt to new environmental conditions.

Genetically distinct coelacanth population off the northern Tanzanian coast.

Since the sensational discovery of a living coelacanth off the east coast of South Africa, the geographic distribution of viable coelacanth populations has been a subject of debate. In the past, the coelacanths off the African mainland were thought to be strays from the Comoros because most coelacanths captured were caught in the waters surrounding the Comoros archipelagos. However, in recent years, a large number of coelacanths were captured off the coast of Tanzania, including nine living specimens observed in a remotely operated vehicles survey. Thus, it is possible that there is a reproducing population inhabiting waters off the Tanzania coast. We have sequenced the complete mitochondrial genomes of 21 Tanzanian and 2 Comoran coelacanths and analyzed these sequences together with two additional full mitochondrial genomes and 47 d-loop sequences from the literature. We found that the coelacanth population off the northern Tanzanian coast is genetically differentiated from those of the southern Tanzania coast and the Comoros, whereas no significant genetic differentiation occurs between the latter two localities. The differentiation between the northern and southern Tanzanian coast populations is consistent with the hypothesis that the existence of northward-flowing ocean current along the Tanzanian coast may reduce or prevent gene flow from the northern to the southern population. Finally, we estimated that the population localized to the southern Tanzanian coast and the Comoros diverged from other coelacanths at least 200,000 y ago. These results indicate that the coelacanths off the northern Tanzania coast are not strays but a genetically distinct group. Our study provides important information for the conservation of this threatened "living fossil."