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.

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.

Exceptional levels of species discovery ameliorate inferences of the biogeography and diversification of an Afrotropical catfish family.

Endeavours in species discovery, particularly the characterisation of cryptic species, have been greatly aided by the application of DNA molecular sequence data to phylogenetic reconstruction and inference of evolutionary and biogeographic processes. However, the extent of cryptic and undescribed diversity remains unclear in tropical freshwaters, where biodiversity is declining at alarming rates. To investigate how data on previously undiscovered biodiversity impacts inferences of biogeography and diversification dynamics, we generated a densely sampled species-level family tree of Afrotropical Mochokidae catfishes (220 valid species) that was ca. 70 % complete. This was achieved through extensive continental sampling specifically targeting the genus Chiloglanis a specialist of the relatively unexplored fast-flowing lotic habitat. Applying multiple species-delimitation methods, we report exceptional levels of species discovery for a vertebrate genus, conservatively delimiting a staggering ca. 50 putative new Chiloglanis species, resulting in a near 80 % increase in species richness for the genus. Biogeographic reconstructions of the family identified the Congo Basin as a critical region in the generation of mochokid diversity, and further revealed complex scenarios for the build-up of continental assemblages of the two most species rich mochokid genera, Synodontis and Chiloglanis. While Syndontis showed most divergence events within freshwater ecoregions consistent with largely in situ diversification, Chiloglanis showed much less aggregation of freshwater ecoregions, suggesting dispersal as a key diversification process in this older group. Despite the significant increase in mochokid diversity identified here, diversification rates were best supported by a constant rate model consistent with patterns in many other tropical continental radiations. While our findings highlight fast-flowing lotic freshwaters as potential hotspots for undescribed and cryptic species diversity, a third of all freshwater fishes are currently threatened with extinction, signifying an urgent need to increase exploration of tropical freshwaters to better characterise and conserve its biodiversity.

Sea temperature is the primary driver of recent and predicted fish community structure across Northeast Atlantic shelf seas.

Climate change has strongly influenced the distribution and abundance of marine fish species, leading to concern about effects of future climate on commercially harvested stocks. Understanding the key drivers of large-scale spatial variation across present-day marine assemblages enables predictions of future change. Here we present a unique analysis of standardised abundance data for 198 marine fish species from across the Northeast Atlantic collected by 23 surveys and 31,502 sampling events between 2005 and 2018. Our analyses of the spatially comprehensive standardised data identified temperature as the key driver of fish community structure across the region, followed by salinity and depth. We employed these key environmental variables to model how climate change will affect both the distributions of individual species and local community structure for the years 2050 and 2100 under multiple emissions scenarios. Our results consistently indicate that projected climate change will lead to shifts in species communities across the entire region. Overall, the greatest community-level changes are predicted at locations with greater warming, with the most pronounced effects at higher latitudes. Based on these results, we suggest that future climate-driven warming will lead to widespread changes in opportunities for commercial fisheries across the region.

Large mesopelagic fish biomass in the Southern Ocean resolved by acoustic properties.

The oceanic mesopelagic zone, 200-1000 m below sea level, holds abundant small fishes that play central roles in ecosystem function. Global mesopelagic fish biomass estimates are increasingly derived using active acoustics, where echosounder-generated signals are emitted, reflected by pelagic organisms and detected by transducers on vessels. Previous studies have interpreted a ubiquitous decline in acoustic reflectance towards the Antarctic continent as a reduction in mesopelagic fish biomass. Here, we use empirical data to estimate species-specific acoustic target strength for the dominant mesopelagic fish of the Scotia Sea in the Southern Ocean. We use these data, alongside estimates of fish relative abundance from net surveys, to interpret signals received in acoustic surveys and calculate mesopelagic biomass of the broader Southern Ocean. We estimate the Southern Ocean mesopelagic fish biomass to be approximately 274 million tonnes if Antarctic krill contribute to the acoustic signal, or 570 million tonnes if mesopelagic fish alone are responsible. These quantities are approximately 1.8 and 3.8 times greater than previous net-based biomass estimates. We also show a peak in fish biomass towards the seasonal ice-edge, corresponding to the preferred feeding grounds of penguins and seals, which may be at risk under future climate change scenarios. Our study provides new insights into the abundance and distributions of ecologically significant mesopelagic fish stocks across the Southern Ocean ecosystem.

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.

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.

Behavior-dependent cis regulation reveals genes and pathways associated with bower building in cichlid fishes.

Many behaviors are associated with heritable genetic variation [Kendler and Greenspan (2006) Am J Psychiatry 163:1683-1694]. Genetic mapping has revealed genomic regions or, in a few cases, specific genes explaining part of this variation [Bendesky and Bargmann (2011) Nat Rev Gen 12:809-820]. However, the genetic basis of behavioral evolution remains unclear. Here we investigate the evolution of an innate extended phenotype, bower building, among cichlid fishes of Lake Malawi. Males build bowers of two types, pits or castles, to attract females for mating. We performed comparative genome-wide analyses of 20 bower-building species and found that these phenotypes have evolved multiple times with thousands of genetic variants strongly associated with this behavior, suggesting a polygenic architecture. Remarkably, F1 hybrids of a pit-digging and a castle-building species perform sequential construction of first a pit and then a castle bower. Analysis of brain gene expression in these hybrids showed that genes near behavior-associated variants display behavior-dependent allele-specific expression with preferential expression of the pit-digging species allele during pit digging and of the castle-building species allele during castle building. These genes are highly enriched for functions related to neurodevelopment and neural plasticity. Our results suggest that natural behaviors are associated with complex genetic architectures that alter behavior via cis-regulatory differences whose effects on gene expression are specific to the behavior itself.

Meta-Fish-Lib: A generalised, dynamic DNA reference library pipeline for metabarcoding of fishes.

The accuracy and reliability of DNA metabarcoding analyses depend on the breadth and quality of the reference libraries that underpin them. However, there are limited options available to obtain and curate the huge volumes of sequence data that are available on public repositories such as NCBI and BOLD. Here, we provide a pipeline to download, clean and annotate mitochondrial DNA sequence data for a given list of fish species. Features of this pipeline include (a) support for multiple metabarcode markers; (b) searches on species synonyms and taxonomic name validation; (c) phylogeny assisted quality control for identification and removal of misannotated sequences; (d) automatically generated coverage reports for each new GenBank release update; and (e) citable, versioned DOIs. As an example we provide a ready-to-use curated reference library for the marine and freshwater fishes of the U.K. To augment this reference library for environmental DNA metabarcoding specifically, we generated 241 new MiFish-12S sequences for 88 U.K. marine species, and make available new primer sets useful for sequencing these. This brings the coverage of common U.K. species for the MiFish-12S fragment to 93%, opening new avenues for scaling up fish metabarcoding across wide spatial gradients. The Meta-Fish-Lib reference library and pipeline is hosted at https://github.com/genner-lab/meta-fish-lib.

Swimbladder morphology masks Southern Ocean mesopelagic fish biomass.

Within the twilight of the oceanic mesopelagic realm, 200-1000 m below sea level, are potentially vast resources of fish. Collectively, these mesopelagic fishes are the most abundant vertebrates on Earth, and this global fish community plays a vital role in the function of oceanic ecosystems. The biomass of these fishes has recently been estimated using acoustic survey methods, which rely on echosounder-generated signals being reflected from gas-filled swimbladders and detected by transducers on vessels. Here, we use X-ray computed tomography scans to demonstrate that several of the most abundant species of mesopelagic fish in the Southern Ocean lack gas-filled swimbladders. We also show using catch data from survey trawls that the fish community switches from fish possessing gas-filled swimbladders to those lacking swimbladders as latitude increases towards the Antarctic continent. Thus, the acoustic surveys that repeatedly show a decrease in mesopelagic fish biomass towards polar environments systematically overlook a large proportion of fish species that dominate polar seas. Importantly, this includes lanternfish species that are key prey items for top predators in the region, including king penguins and elephant seals. This latitudinal community switch, from gas to non-gas dominance, has considerable implications for acoustic biomass estimation, ecosystem modelling and long-term monitoring of species at risk from climate change and potential exploitation.

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

Minimalist barcodes for sponges: a case study classifying African freshwater Spongillida.

African sponges, particularly freshwater sponges, are understudied relative to demosponges in most other geographical regions. Freshwater sponges (Spongillida) likely share a common ancestor; however, their evolutionary history, particularly during their radiation into endemic and allegedly cosmopolitan groups, is unclear. Freshwater sponges of at least 58 species of 17 genera and four families are described from Central and Eastern Africa, but the diversity is underestimated due to limited distinguishable morphological features. The discovery of additional cryptic species is very likely with the use of molecular techniques such as DNA barcoding. The Royal Museum of Central Africa (MRAC, Tervuren, Belgium) hosts one of the largest collections of (Central) African freshwater sponge type material. Type specimens in theory constitute ideal targets for molecular taxonomy; however, the success is frequently hampered by DNA degradation and deamination, which are a consequence of suboptimal preservation techniques. Therefore, we genotyped African demosponge holotype material of the MRAC with specific short primers suitable for degenerated tissue and compare the results with the current, morphology-based classification. Our results demonstrate the utility of minimalistic barcodes for identification of sponges, potentially enabling efficient identification of individuals in taxonomic or metabarcoding studies, and highlight inconsistencies in the current freshwater sponge classification.

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.

Staying out of the heat: how habitat use is determined by local temperature.

Atlantic cod (Gadus morhua) and the Tvedestrand fjord on the Norwegian Skagerrak coast studied by Freitas et al. (). Photographs by Øystein Paulsen (left), and Institute of Marine Research, Norway (right). In Focus: Freitas, C., Olsen, E. M., Knutsen, H., Albretsen, J. & Moland, E. (2016) Temperature-associated habitat selection in a cold-water marine fish. Journal of Animal Ecology, 85, 611-613. In the marine environment, species distributions are closely linked to temperature gradients, but how individual behaviour is affected by local temperatures is less well understood. Freitas et al. () tracked Atlantic cod within a Norwegian fjord using electronic acoustic tags. They showed that when surface waters were warm, cod occupied the cold deep non-vegetated habitats. However, when surface waters cooled, fish moved into shallow seagrass and macroalgae beds that were previously out-of-bounds. The study provides a clear example of how thermal regimes determine habitat use over fine spatial and temporal scales, with potential implications for population dynamics under climate warming.

Patterns of species range evolution in Indo-Pacific reef assemblages reveal the Coral Triangle as a net source of transoceanic diversity.

The Coral Triangle in the Indo-Pacific is a region renowned for exceptional marine biodiversity. The area could have acted as a 'centre of origin' where speciation has been prolific or a 'centre of survival' by providing refuge during major environmental shifts such as sea-level changes. The region could also have acted as a 'centre of accumulation' for species with origins outside of the Coral Triangle, owing to it being at a central position between the Indian and Pacific oceans. Here, we investigated support for these hypotheses using population-level DNA sequence-based reconstructions of the range evolution of 45 species (314 populations) of Indo-Pacific reef-associated organisms. Our results show that populations undergoing the most ancient establishment were significantly more likely to be closer to the centre of the Coral Triangle than to peripheral locations. The data are consistent with the Coral Triangle being a net source of coral-reef biodiversity for the Indo-Pacific region, suggesting that the region has acted primarily as a centre of survival, a centre of origin or both. These results provide evidence of how a key location can influence the large-scale distributions of biodiversity over evolutionary timescales.

Secondary contact seeds phenotypic novelty in cichlid fishes.

Theory proposes that genomic admixture between formerly reproductively isolated populations can generate phenotypic novelty for selection to act upon. Secondary contact may therefore be a significant promoter of phenotypic novelty that allows species to overcome environmental challenges and adapt to novel environments, including during adaptive radiation. To date, this has largely been considered from the perspective of interspecific hybridization at contact zones. However, it is also possible that this process occurs more commonly between natural populations of a single species, and thus its importance in adaptive evolution may have been underestimated. In this study, we tested the consequences of genomic introgression during apparent secondary contact between phenotypically similar lineages of the riverine cichlid fish Astatotilapia calliptera. We provide population genetic evidence of a secondary contact zone in the wild, and then demonstrate using mate-choice experiments that both lineages can reproduce together successfully in laboratory conditions. Finally, we show that genomically admixed individuals display extreme phenotypes not observed in the parental lineages. Collectively, the evidence shows that secondary contact can drive the evolution of phenotypic novelty, suggesting that pulses of secondary contact may repeatedly seed genetic novelty, which when coupled with ecological opportunity could promote rapid adaptive evolution in natural circumstances.

Warming shelf seas drive the subtropicalization of European pelagic fish communities.

Pelagic fishes are among the most ecologically and economically important fish species in European seas. In principle, these pelagic fishes have potential to demonstrate rapid abundance and distribution shifts in response to climatic variability due to their high adult motility, planktonic larval stages, and low dependence on benthic habitat for food or shelter during their life histories. Here, we provide evidence of substantial climate-driven changes to the structure of pelagic fish communities in European shelf seas. We investigated the patterns of species-level change using catch records from 57,870 fisheries-independent survey trawls from across European continental shelf region between 1965 and 2012. We analysed changes in the distribution and rate of occurrence of the six most common species, and observed a strong subtropicalization of the North Sea and Baltic Sea assemblages. These areas have shifted away from cold-water assemblages typically characterized by Atlantic herring and European sprat from the 1960s to 1980s, to warmer-water assemblages including Atlantic mackerel, Atlantic horse mackerel, European pilchard and European anchovy from the 1990s onwards. We next investigated if warming sea temperatures have forced these changes using temporally comprehensive data from the North Sea region. Our models indicated the primary driver of change in these species has been sea surface temperatures in all cases. Together, these analyses highlight how individual species responses have combined to result in a dramatic subtropicalization of the pelagic fish assemblage of the European continental shelf.