Genetic variation and demographic history of the Haplochromis laparogramma group of Lake Victoria-An analysis based on SINEs and mitochondrial DNA

Strong species structure but weak geographical structure in demersal Lake Victoria cichlids

Studying phenotypic and genetic differentiation between very young species can be very informative with regard to learning about processes of speciation. Identifying and characterizing genetic species structure and distinguishing it from spatial genetic structure within a species is a prerequisite for this and is often not given sufficient attention. Young radiations of cichlid fish are classical speciation study systems. However, it is only during the past decade that population genomics based on next-generation sequencing has begun to provide the power to resolve species and distinguish speciation from spatial population structure for the youngest of these radiations. The Lake Victoria haplochromine cichlids constitute the youngest large cichlid fish radiation, probably <20,000 years old. Earlier work showed that communities of rocky reef cichlids are composed of many reciprocally monophyletic species despite their very recent origins. Here, we build on this work by studying assemblages of offshore demersal cichlids, adding analyses of within-species spatial structure to the sympatric species structure. We sampled seven multispecies communities along a 6-km-long transect from one side of the Mwanza Gulf to the other side. We investigated whether phenotypically diagnosed putative species are reciprocally monophyletic and whether such monophyly is stable across species geographic ranges. We show that all species are genetically strongly differentiated in sympatry, that they are reciprocally monophyletic, and that monophyly is stable across distribution ranges. We found significant differentiation between geographically distinct populations in two species, but no or weak isolation by distance. We further found subtle but significant morphological differences between all species and a linear relationship between genomic and morphological distance which suggests that differences in morphology begin to accumulate after speciation has already affected genome-wide restrictions of gene flow.

Patterns of genomic differentiation between two Lake Victoria cichlid species, Haplochromis pyrrhocephalus and H. sp. 'macula'

BACKGROUND

The molecular basis of the incipient stage of speciation is still poorly understood. Cichlid fish species in Lake Victoria are a prime example of recent speciation events and a suitable system to study the adaptation and reproductive isolation of species.

RESULTS

Here, we report the pattern of genomic differentiation between two Lake Victoria cichlid species collected in sympatry, Haplochromis pyrrhocephalus and H. sp. 'macula,' based on the pooled genome sequences of 20 individuals of each species. Despite their ecological differences, population genomics analyses demonstrate that the two species are very close to a single panmictic population due to extensive gene flow. However, we identified 21 highly differentiated short genomic regions with fixed nucleotide differences. At least 15 of these regions contained genes with predicted roles in adaptation and reproductive isolation, such as visual adaptation, circadian clock, developmental processes, adaptation to hypoxia, and sexual selection. The nonsynonymous fixed differences in one of these genes, LWS, were reported as substitutions causing shift in absorption spectra of LWS pigments. Fixed differences were found in the promoter regions of four other differentially expressed genes, indicating that these substitutions may alter gene expression levels.

CONCLUSIONS

These diverged short genomic regions may have contributed to the differentiation of two ecologically different species. Moreover, the origins of adaptive variants within the differentiated regions predate the geological formation of Lake Victoria; thus Lake Victoria cichlid species diversified via selection on standing genetic variation.

Visual adaptation in Lake Victoria cichlid fishes: depth-related variation of color and scotopic opsins in species from sand/mud bottoms

Background

For Lake Victoria cichlid species inhabiting rocky substrates with differing light regimes, it has been proposed that adaptation of the long-wavelength-sensitive (LWS) opsin gene triggered speciation by sensory drive through color signal divergence. The extensive and continuous sand/mud substrates are also species-rich, and a correlation between male nuptial coloration and the absorption of LWS pigments has been reported. However, the factors driving genetic and functional diversity of LWS pigments in sand/mud habitats are still unresolved.

Results

To address this issue, nucleotide sequences of eight opsin genes were compared in ten Lake Victoria cichlid species collected from sand/mud bottoms. Among eight opsins, the LWS and rod-opsin (RH1) alleles were diversified and one particular allele was dominant or fixed in each species. Natural selection has acted on and fixed LWS alleles in each species. The functions of LWS and RH1 alleles were measured by absorption of reconstituted A1- and A2-derived visual pigments. The absorption of pigments from RH1 alleles most common in deep water were largely shifted toward red, whereas those of LWS alleles were largely shifted toward blue in both A1 and A2 pigments. In both RH1 and LWS pigments, A2-derived pigments were closer to the dominant light in deep water, suggesting the possibility of the adaptation of A2-derived pigments to depth-dependent light regimes.

Conclusions

The RH1 and LWS sequences may be diversified for adaptation of A2-derived pigments to different light environments in sand/mud substrates. Diversification of the LWS alleles may have originally taken place in riverine environments, with a new mutation occurring subsequently in Lake Victoria.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1040-x) contains supplementary material, which is available to authorized users.

A microsatellite-based genetic linkage map and putative sex-determining genomic regions in Lake Victoria cichlids

Cichlid fishes in East Africa have undergone extensive adaptive radiation, which has led to spectacular diversity in their morphology and ecology. To date, genetic linkage maps have been constructed for several tilapias (riverine), Astatotilapia burtoni (Lake Tanganyika), and hybrid lines of Lake Malawi cichlids to facilitate genome-wide comparative analyses. In the present study, we constructed a genetic linkage map of the hybrid line of Lake Victoria cichlids, so that maps of cichlids from all the major areas of East Africa will be available. The genetic linkage map shown here is derived from the F2 progeny of an interspecific cross between Haplochromis chilotes and Haplochromis sauvagei and is based on 184 microsatellite and two single-nucleotide polymorphism (SNP) markers. Most of the microsatellite markers used in the present study were originally designed for other genetic linkage maps, allowing us to directly compare each linkage group (LG) among different cichlid groups. We found 25 LGs, the total length of which was 1133.2cM with an average marker spacing of about 6.09cM. Our subsequent linkage mapping analysis identified two putative sex-determining loci in cichlids. Interestingly, one of these two loci is located on cichlid LG5, on which the female heterogametic ZW locus and several quantitative trait loci (QTLs) related to adaptive evolution have been reported in Lake Malawi cichlids. We also found that V1R1 and V1R2, candidate genes for the fish pheromone receptor, are located very close to the recently detected sex-determining locus on cichlid LG5. The genetic linkage map study presented here may provide a valuable foundation for studying the chromosomal evolution of East African cichlids and the possible role of sex chromosomes in generating their genomic diversity.

Key processes for Cheirolophus (Asteraceae) diversification on oceanic islands inferred from AFLP data

The radiation of the genus Cheirolophus (Asteraceae) in Macaronesia constitutes a spectacular case of rapid diversification on oceanic islands. Twenty species - nine of them included in the IUCN Red List of Threatened Species - have been described to date inhabiting the Madeiran and Canarian archipelagos. A previous phylogenetic study revealed that the diversification of Cheirolophus in Macaronesia started less than 2 Ma. As a result of such an explosive speciation process, limited phylogenetic resolution was reported, mainly due to the low variability of the employed molecular markers. In the present study, we used highly polymorphic AFLP markers to i) evaluate species' boundaries, ii) infer their evolutionary relationships and iii) investigate the patterns of genetic diversity in relation to the potential processes likely involved in the radiation of Cheirolophus. One hundred and seventy-two individuals representing all Macaronesian Cheirolophus species were analysed using 249 AFLP loci. Our results suggest that geographic isolation played an important role in this radiation process. This was likely driven by the combination of poor gene flow capacity and a good ability for sporadic long-distance colonisations. In addition, we also found some traces of introgression and incipient ecological adaptation, which could have further enhanced the extraordinary diversification of Cheirolophus in Macaronesia. Last, we hypothesize that current threat categories assigned to Macaronesian Cheirolophus species do not reflect their respective evolutionary relevance, so future evaluations of their conservation status should take into account the results presented here.

Genetic structure of pelagic and littoral cichlid fishes from Lake Victoria

The approximately 700 species of cichlids found in Lake Victoria in East Africa are thought to have evolved over a short period of time, and they represent one of the largest known examples of adaptive radiation. To understand the processes that are driving this spectacular radiation, we must determine the present genetic structure of these species and elucidate how this structure relates to the ecological conditions that caused their adaptation. We analyzed the genetic structure of two pelagic and seven littoral species sampled from the southeast area of Lake Victoria using sequences from the mtDNA control region and 12 microsatellite loci as markers. Using a Bayesian model-based clustering method to analyze the microsatellite data, we separated these nine species into four groups: one group composed of pelagic species and another three groups composed mainly of rocky-shore species. Furthermore, we found significant levels of genetic variation between species within each group at both marker loci using analysis of molecular variance (AMOVA), although the nine species often shared mtDNA haplotypes. We also found significant levels of genetic variation between populations within species. These results suggest that initial groupings, some of which appear to have been related to habitat differences, as well as divergence between species within groups took place among the cichlid species of Lake Victoria.

Two new species of zooplanktivorous haplochromine cichlids from Lake Victoria, Tanzania

Two new species of zooplanktivorous haplochromine cichlids from Lake Victoria, Tanzania, are described and illustrated. These species closely resemble each other. Their affinities to other zooplanktivorous haplochromines from Lake Victoria are discussed. Haplochromis argens sp. n., which featured under nicknames (mainly Haplochromis "argens") in more than 50 papers, was caught both in the Mwanza Gulf and the Emin Pasha Gulf, whereas Haplochromis goldschmidti sp. n. was only found in the Emin Pasha Gulf. Of the latter species only males are available, but it seems unlikely that it represents a case of male colour polymorphism as several presumably unrelated characters differ in sympatry between the two species, suggesting that there is no gene flow. Statistical analysis revealed that the overall difference between the two species is greater than that between the populations from the two locations. Body depth of the two species in sympatry in the Emin Pasha Gulf was more similar than that of Haplochromis goldschmidti sp. n. and the allopatric population of Haplochromis argens sp. n. from the Mwanza Gulf,which mayindicate an overall environmental effect. However, several measurements related to the width of snout and mouth differed more between the populations of the two species in sympatry than between the allopatric populations. In contrast to a group of zooplanktivorous species that recovered successfully after environmental changes in the lake, Haplochromis argens sp. n. is among a group that became extremely rare and probably is in danger of extinction; the conservation status of Haplochromis goldschmidti sp. n. is currently unknown.

Genome-wide RAD sequence data provide unprecedented resolution of species boundaries and relationships in the Lake Victoria cichlid adaptive radiation

Although population genomic studies using next generation sequencing (NGS) data are becoming increasingly common, studies focusing on phylogenetic inference using these data are in their infancy. Here, we use NGS data generated from reduced representation genomic libraries of restriction-site-associated DNA (RAD) markers to infer phylogenetic relationships among 16 species of cichlid fishes from a single rocky island community within Lake Victoria's cichlid adaptive radiation. Previous attempts at sequence-based phylogenetic analyses in Victoria cichlids have shown extensive sharing of genetic variation among species and no resolution of species or higher-level relationships. These patterns have generally been attributed to the very recent origin (<15,000 years) of the radiation, and ongoing hybridization between species. We show that as we increase the amount of sequence data used in phylogenetic analyses, we produce phylogenetic trees with unprecedented resolution for this group. In trees derived from our largest data supermatrices (3 to >5.8 million base pairs in width), species are reciprocally monophyletic with high bootstrap support, and the majority of internal branches on the tree have high support. Given the difficulty of the phylogenetic problem that the Lake Victoria cichlid adaptive radiation represents, these results are striking. The strict interpretation of the topologies we present here warrants caution because many questions remain about phylogenetic inference with very large genomic data set and because we can with the current analysis not distinguish between effects of shared ancestry and post-speciation gene flow. However, these results provide the first conclusive evidence for the monophyly of species in the Lake Victoria cichlid radiation and demonstrate the power that NGS data sets hold to resolve even the most difficult of phylogenetic challenges.

Correlation between nuptial colors and visual sensitivities tuned by opsins leads to species richness in sympatric Lake Victoria cichlid fishes

Reproductive isolation that prevents interspecific hybridization between closely related coexisting species maintains sympatric species diversity. One of the reproductive isolations is mate choice based on color signals (breeding color perceived by color vision). This is well known in several animal taxa, yet little is known about its genetic and molecular mechanism. Lake Victoria cichlid fishes are thought to be an example of sympatric species diversity. In the species inhabiting different light environments in rocky shore, speciation by sensory drive through color signals has been proposed by analyses of the long wavelength-sensitive (LWS) opsin gene and the male nuptial coloration. However, the genetic and molecular mechanism of how diversity of sympatric species occurring in the same habitat is maintained remains unknown. To address this issue, we determined nucleotide sequences of eight opsins of six sympatric species collected from a sandy-muddy shore--an ideal model system for studying sympatric species. Among eight opsins, the LWS and RH1 alleles were diversified and one particular allele is dominant or fixed in each species, and we propose that this is due to natural selection. The functions of their LWS alleles were also diversified as shown by absorption measurements of reconstituted visual pigments. To analyze the relationship between nuptial coloration and the absorption of LWS pigments, we systematically evaluated and defined nuptial coloration. We showed that the coloration was species specific with respect to hue and significantly differentiated by the index values of hue (dominant wavelength: λ(d)). The λ(d) value of the male nuptial coloration correlated with the absorption of LWS pigments from all the species, suggesting that reproductive isolation through mate choice using color signals may prevent sympatric interspecific hybridization, thereby maintaining the species diversity in sympatric species in Lake Victoria.

Community Genetics Reveal Elevated Levels of Sympatric Gene Flow among Morphologically Similar but Not among Morphologically Dissimilar Species of Lake Victoria Cichlid Fish

We examined genetic structure among five species of Lake Victoria haplochromine cichlids in four island communities, using a full factorial sampling design that compared genetic differentiation between pairs of species and populations of varying morphological similarity and geographical proximity. We found that allopatric conspecific populations were on average significantly more strongly differentiated than sympatric heterospecific populations of morphologically similar species. Allopatric heterospecific populations of morphologically dissimilar species were most differentiated. Our work demonstrates that phenotypic divergence can be maintained and perhaps even evolve in sympatry despite considerable gene flow between species. Conversely, phenotypic resemblance among conspecific populations can be maintained despite geographical isolation. Additionally we show that anthropogenically increased hybridization does not affect all sympatric species evenly but predominantly affects morphologically similar and closely related species. This has important implications for the evolution of reproductive isolation between species These findings are also consistent with the hypothesis of speciation reversal due to weakening of divergent selection and reproductive isolation as a consequence of habitat homogenization and offers an evolutionary mechanistic explanation for the observation that species poor assemblages in turbid areas of the lake are characterized by just one or two species in each of a few morphologically distinct genera.

Ancient origin and recent divergence of a haplochromine cichlid lineage from isolated water bodies in the East African Rift system

Phylogenetic analysis identified haplochromine cichlids from isolated water bodies in the eastern branch of the East African Rift system as an ancient lineage separated from their western sister group in the course of the South Kenyan-North Tanzanian rift system formation. Within this lineage, the close phylogenetic relatedness among taxa indicates a recent common ancestry and historical connections between now separated water bodies. In connection with a total lack of local genetic diversity attributable to population bottlenecks, the data suggest cycles of extinction and colonization in the unstable habitat provided by the small lakes and rivers in this geologically highly active area.