A whole-body micro-CT scan library that captures the skeletal diversity of Lake Malawi cichlid fishes.

Here we describe a dataset of freely available, readily processed, whole-body µCT-scans of 56 species (116 specimens) of Lake Malawi cichlid fishes that captures a considerable majority of the morphological variation present in this remarkable adaptive radiation. We contextualise the scanned specimens within a discussion of their respective ecomorphological groupings and suggest possible macroevolutionary studies that could be conducted with these data. In addition, we describe a methodology to efficiently µCT-scan (on average) 23 specimens per hour, limiting scanning time and alleviating the financial cost whilst maintaining high resolution. We demonstrate the utility of this method by reconstructing 3D models of multiple bones from multiple specimens within the dataset. We hope this dataset will enable further morphological study of this fascinating system and permit wider-scale comparisons with other cichlid adaptive radiations.

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.

Multiple environmental stressors affect predation pressure in a tropical freshwater system.

Environmental change can alter predator-prey dynamics. However, studying predators in the context of co-occurring environmental stressors remains rare, especially under field conditions. Using in situ filming, we examined how multiple stressors, including temperature and turbidity, impact the distribution and behaviour of wild fish predators of Trinidadian guppies (Poecilia reticulata). The measured environmental variables accounted for 17.6% of variance in predator species composition. While predator species differed in their associations with environmental variables, the overall prevalence of predators was greatest in slow flowing, deeper, warmer and less turbid habitats. Moreover, these warmer and less turbid habitats were associated with earlier visits to the prey stimulus by predators, and more frequent predator visits and attacks. Our findings highlight the need to consider ecological complexity, such as co-occurring stressors, to better understand how environmental change affects predator-prey interactions.

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.

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.

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.

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.

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.

Lateral line morphology, sensory perception and collective behaviour in African cichlid fish.

The lateral line system of fishes provides cues for collective behaviour, such as shoaling, but it remains unclear how anatomical lateral line variation leads to behavioural differences among species. Here we studied associations between lateral line morphology and collective behaviour using two morphologically divergent species and their second-generation hybrids. We identify collective behaviours associated with variation in canal and superficial lateral line morphology, with closer proximities to neighbouring fish associated with larger canal pore sizes and fewer superficial neuromasts. A mechanistic understanding of the observed associations was provided by hydrodynamic modelling of an artificial lateral line sensor, which showed that simulated canal-based neuromasts were less susceptible to saturation during unidirectional movement than simulated superficial neuromasts, while increasing the canal pore size of the simulated lateral line sensor elevated sensitivity to vortices shed by neighbouring fish. Our results propose a mechanism behind lateral line flow sensing during collective behaviour in fishes.

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.

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.

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.

Mapping epigenetic divergence in the massive radiation of Lake Malawi cichlid fishes.

Epigenetic variation modulates gene expression and can be heritable. However, knowledge of the contribution of epigenetic divergence to adaptive diversification in nature remains limited. The massive evolutionary radiation of Lake Malawi cichlid fishes displaying extensive phenotypic diversity despite extremely low sequence divergence is an excellent system to study the epigenomic contribution to adaptation. Here, we present a comparative genome-wide methylome and transcriptome study, focussing on liver and muscle tissues in phenotypically divergent cichlid species. In both tissues we find substantial methylome divergence among species. Differentially methylated regions (DMR), enriched in evolutionary young transposons, are associated with transcription changes of ecologically-relevant genes related to energy expenditure and lipid metabolism, pointing to a link between dietary ecology and methylome divergence. Unexpectedly, half of all species-specific DMRs are shared across tissues and are enriched in developmental genes, likely reflecting distinct epigenetic developmental programmes. Our study reveals substantial methylome divergence in closely-related cichlid fishes and represents a resource to study the role of epigenetics in species diversification.

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.

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.

Correction: Schistosoma species detection by environmental DNA assays in African freshwaters.

[This corrects the article DOI: 10.1371/journal.pntd.0008129.].

Schistosoma species detection by environmental DNA assays in African freshwaters.

BACKGROUND: Schistosomiasis is a neglected tropical parasitic disease associated with severe pathology, mortality and economic loss worldwide. Programs for disease control may benefit from specific and sensitive diagnostic methods to detect Schistosoma trematodes in aquatic environments. Here we report the development of novel environmental DNA (eDNA) qPCR assays for the presence of the human-infecting species Schistosoma mansoni, S. haematobium and S. japonicum. METHODOLOGY/PRINCIPAL FINDINGS: We first tested the specificity of the assays across the three species using genomic DNA preparations which showed successful amplification of target sequences with no cross amplification between the three focal species. In addition, we evaluated the specificity of the assays using synthetic DNA of multiple Schistosoma species, and demonstrated a high overall specificity; however, S. japonicum and S. haematobium assays showed cross-species amplification with very closely-related species. We next tested the effectiveness of the S. mansoni assay using eDNA samples from aquaria containing infected host gastropods, with the target species revealed as present in all infected aquaria. Finally, we evaluated the effectiveness of the S. mansoni and S. haematobium assays using eDNA samples from eight discrete natural freshwater sites in Tanzania, and demonstrated strong correspondence between infection status established using eDNA and conventional assays of parasite prevalence in host snails. CONCLUSIONS/SIGNIFICANCE: Collectively, our results suggest that eDNA monitoring is able to detect schistosomes in freshwater bodies, but refinement of the field sampling, storage and assay methods are likely to optimise its performance. We anticipate that environmental DNA-based approaches will help to inform epidemiological studies and contribute to efforts to control and eliminate schistosomiasis in endemic areas.

Climate Change Drives Poleward Increases and Equatorward Declines in Marine Species.

Marine environments have increased in temperature by an average of 1°C since pre-industrial (1850) times [1]. Given that species ranges are closely allied to physiological thermal tolerances in marine organisms [2], it may therefore be expected that ocean warming would lead to abundance increases at poleward side of ranges and abundance declines toward the equator [3]. Here, we report a global analysis of abundance trends of 304 widely distributed marine species over the last century, across a range of taxonomic groups from phytoplankton to fish and marine mammals. Specifically, using a literature database, we investigate the extent that the direction and strength of long-term species abundance changes depend on the sampled location within the latitudinal range of species. Our results show that abundance increases have been most prominent where sampling has taken place at the poleward side of species ranges, and abundance declines have been most prominent where sampling has taken place at the equatorward side of species ranges. These data provide evidence of omnipresent large-scale changes in abundance of marine species consistent with warming over the last century and suggest that adaptation has not provided a buffer against the negative effects of warmer conditions at the equatorward extent of species ranges. On the basis of these results, we suggest that projected sea temperature increases of up to 1.5°C over pre-industrial levels by 2050 [4] will continue to drive latitudinal abundance shifts in marine species, including those of importance for coastal livelihoods.