Genomic Signatures for Species-Specific Adaptation in Lake Victoria Cichlids Derived from Large-Scale Standing Genetic Variation.

The cichlids of Lake Victoria are a textbook example of adaptive radiation, as >500 endemic species arose in just 14,600 years. The degree of genetic differentiation among species is very low due to the short period of time after the radiation, which allows us to ascertain highly differentiated genes that are strong candidates for driving speciation and adaptation. Previous studies have revealed the critical contribution of vision to speciation by showing the existence of highly differentiated alleles in the visual opsin gene among species with different habitat depths. In contrast, the processes of species-specific adaptation to different ecological backgrounds remain to be investigated. Here, we used genome-wide comparative analyses of three species of Lake Victoria cichlids that inhabit different environments-Haplochromis chilotes, H. sauvagei, and Lithochromis rufus-to elucidate the processes of adaptation by estimating population history and by searching for candidate genes that contribute to adaptation. The patterns of changes in population size were quite distinct among the species according to their habitats. We identified many novel adaptive candidate genes, some of which had surprisingly long divergent haplotypes between species, thus showing the footprint of selective sweep events. Molecular phylogenetic analyses revealed that a large fraction of the allelic diversity among Lake Victoria cichlids was derived from standing genetic variation that originated before the adaptive radiation. Our analyses uncovered the processes of species-specific adaptation of Lake Victoria cichlids and the complexity of the genomic substrate that facilitated this adaptation.

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.

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.

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.

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.

Multiple episodic evolution events in V1R receptor genes of East-African cichlids.

Fish use olfaction to detect a variety of nonvolatile chemical signals, and thus, this sense is key to survival and communication. However, the contribution of the olfactory sense to social-especially reproductive-interactions in cichlids is still controversial. To obtain insights into this issue, we investigated the genes encoding V1Rs-possible candidates for reproductive pheromone receptors-among East-African cichlids. Interestingly, we found an excess of nonsynonymous over synonymous substitutions in four of six V1R genes in multiple cichlid lineages. First, we found that highly dimorphic V1R2 allele groups were shared among the cichlids inhabiting all East-African Great Lakes emerged through the episodic accumulation of the nonsynonymous substitutions prior to the radiation of the Lake Tanganyika species flock. We further detected such episodic events in V1R1 of the tribe Tropheini, and in V1R3 and V1R6 of the tribe Trematocarini. The excess of nonsynonymous substitutions in these examples were indicated as dN/dS > 1, which were all statistically significant by Fisher's exact test. Furthermore, we speculate that the amino acid changes in these episodic events are likely functional switch because they occurred in the putative ligand-binding pocket. Our finding of the occurrence of multiple episodic events and the unexpected gene diversity in one unique gene family is suggestive of the contribution of the V1R to the species diversification and the social interaction in cichlids.

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.

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.

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.

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

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

More than 500 endemic haplochromine cichlid species inhabit Lake Victoria. This striking species diversity is a classical example of recent explosive adaptive radiation thought to have happened within the last approximately 15,000 years. In this study, we examined the population structure and historical demography of 3 pelagic haplochromine cichlid species that resemble in morphology and have similar niche, Haplochromis (Yssichromis) laparogramma, Haplochromis (Y.) pyrrhocephalus, and Haplochromis (Y.) sp. "glaucocephalus". We investigated the sequences of the mitochondrial DNA control region and the insertion patterns of short interspersed elements (SINEs) of 759 individuals. We show that sympatric forms are genetically differentiated in 4 of 6 cases, but we also found apparent weakening of the genetic differentiation in areas with turbid water. We estimated the timings of population expansion and species divergence to coincide with the refilling of the lake at the Pleistocene/Holocene boundary. We also found that estimates can be altered significantly by the choice of the shape of the molecular clock. If we employ the nonlinear clock model of evolutionary rates in which the rates are higher towards the recent, the population expansion was dated at around the event of desiccation of the lake ca. 17,000 YBP. Thus, we succeeded in clarifying the species and population structure of closely related Lake Victoria cichlids and in showing the importance of applying appropriate clock calibrations in elucidating recent evolutionary events.

Population structure of two closely related pelagic cichlids in Lake Victoria, Haplochromis pyrrhocephalus and H. laparogramma.

Cichlid fishes in Lake Victoria show spectacular diversification that is thought to be recent. Therefore, by investigating those fishes, we may be able to elucidate recently completed or ongoing speciation processes. We studied the population structures of two closely related pelagic cichlid species, Haplochromis pyrrhocephalus and H. laparogramma, using a mitochondrial DNA locus and 12 nuclear microsatellite loci as putative neutral markers. Ten and two populations of H. pyrrhocephalus and H. laparogramma, respectively, were sampled from the southern part of Lake Victoria. We grouped those 12 populations into four mutually differentiated regional populations, one of which consisted of the two H. laparogramma populations. The levels of differentiation were substantial at the mitochondrial locus (F(ST) = 0.03-0.54), but very low at microsatellite loci (R(ST) = 0.008-0.116). The data from both types of loci indicated that the regional population of H. laparogramma was first separated from those of H. pyrrhocephalus if we set aside one erratic population of H. pyrrhocephalus. The data also suggested recent population expansions of the two species, the time scales for which were estimated to be on the order of 10(4)-10(5) years. These data suggested that dynamic speciation processes accompanied occasional spawning of new species and population size changes in this lake.