Natural Selection of the Rhodopsin Gene During the Adaptive Radiation of East African Great Lakes Cichlid Fishes

Tohru Sugawara, Yohey Terai and Norihiro Okada

Graduate School of Bioscience and Biotechnology, Tokyo Institute of Technology, Yokohama, Japan

The Great Lakes of the East African Rift Valley, Lakes Victoria, Malawi, and Tanganyika, harbor approximately 200, 400, and 170 endemic species of cichlid fishes, respectively (Fryer and Iles 1972Citation , pp. 6–590; Greenwood 1991Citation ). These fishes have fascinated evolutionary biologists as spectacular examples of explosive adaptive radiation among living vertebrates because the fishes have adapted to a variety of niches in the lakes (Fryer and Iles 1972Citation ; Greenwood 1984Citation ). Accordingly, they are extremely diverse, both ecologically and morphologically, despite having evolved over a very short period (Meyer et al. 1990Citation ; Sturmbauer and Meyer 1992Citation ; Johnson et al. 1996Citation ; Takahashi et al. 2001Citation ). The visual systems of these fishes are of particular interest because they are important for feeding (Fryer and Iles 1972Citation ) as well as for mate choice (Seehausen, van Alphen, and Witte 1997Citation ; Seehausen and van Alphen 1998Citation ), and cichlids have adapted to different photic conditions that vary with depth and time of day (Fryer and Iles 1972Citation ; Coulter 1991Citation ). Furthermore, differences in visual capabilities among coexisting cichlid species may reduce the competition for resources (van der Meer and Bowmaker 1995Citation ). Therefore, it is reasonable to postulate that genes that are essential for visual acuity must have evolved adaptively during adaptive radiation of the cichlids in these African Great Lakes. We attempted to identify such genes by focusing on cichlid opsin genes.

Visual pigments are composed of a light-absorbing chromophore, typically 11-cis-retinal, that is bound to the protein opsin (Wald 1968Citation ). The absorption spectrum of the chromophore can be altered by amino acid substitutions within opsin, and key substitution sites are located in the seven transmembrane {alpha}-helices that lie close to the retinal-binding pocket (Yokoyama R. and Yokoyama S. 1990Citation ; Yokoyama S. and Yokoyama R. 1996Citation ; Kochendoerfer et al. 1999Citation ). Cichlids have six known opsin genes whose products are sensitive to different wavelengths of visible-ultraviolet light, namely SWS-1 (ultraviolet), SWS-2A and SWS-2B (short wavelength), RH2 (midwavelength), LWS (long wavelength) (Carleton, Harosi, and Kocher 2000Citation ; Carleton and Kocher 2001Citation ), and rhodopsin (GenBank accession number AF315354). Rhodopsin is localized to rod cells that enable black and white images to be seen in dim light (Yokoyama R. and Yokoyama S. 1990Citation ; Yokoyama S. and Yokoyama R. 1996Citation ). The others are localized within cone cells that mediate color vision in bright light (Yokoyama S. and Yokoyama R. 1996Citation ; Yokoyama 2000Citation ). The evolution of visual pigments is a prime example of molecular adaptation in vertebrates (Yokoyama S. and Yokoyama R. 1996Citation ; Yokoyama 2000Citation ), and several such cases have been reported in diverse species such as primates (Shyue et al. 1995Citation ), birds (Yokoyama, Radlwimmer, and Blow 2000Citation ), and fish (Hunt et al. 1996Citation ; Yokoyama et al. 1999Citation ). As a first step toward understanding the molecular basis for the adaptation of opsins in cichlids, we cloned and characterized a cichlid rhodopsin gene.

To identify possible adaptive changes in the rhodopsin gene during cichlid adaptive radiation, we sequenced the coding region (924 bp) from 19 species of the major cichlid lineages from the African Great Lakes and from five species of riverine cichlids (fig. 1a ). Southern blot analysis showed that the cichlid rhodopsin gene is encoded in a single gene at this level of analysis (data not shown). The gene was amplified using a pair of polymerase chain reaction (PCR) primers designated RhodF1 (5'-AGCCAGAAGAAACACCTCTGAAG-3') and RhodR1 (5'-TTGGAGGCAGTAGAAGATGCT-3') on the basis of the rhodopsin sequence of Astatotilapia burtoni (GenBank accession number AF315354). Amplifications were carried out in a PTC-100 Programmable Thermal Controller (MJ Research). Genomic DNA (50 ng) was added to the reaction mixture containing 1x ExTaq PCR buffer, 0.2 mM each of the four deoxynucleoside triphosphates, 0.1 µM primers, and 1 U ExTaq polymerase (Takara, Shiga, Japan). The PCR program consisted of a denaturation step for 3 min at 94°C, followed by 30 cycles consisting of 30 s denaturation at 94°C, 30 s annealing at 55°C, and 2 min extension at 72°C. Direct sequencing of the amplified gene was performed using a DNA sequencing kit (Applied Biosystems) and the Applied Biosystems Automated 3100 Sequencer. The nucleotide sequences were deposited in GenBank under accession numbers AB084924AB084947.



View larger version (49K):
[in this window]
[in a new window]
 
Fig. 1.—(a) Phylogeny of the cichlid species used in this study, and the average values of Dn/Ds within each group. The consensus phylogenetic tree of cichlids is based on molecular data from several sources (Nishida 1991Citation ; Kocher et al. 1993Citation ; Sturmbauer and Meyer 1993Citation ; Mayer, Tichy, and Klein 1998Citation ; Takahashi et al. 2001Citation ). The A292S above certain branches indicates that this amino acid substitution occurred in the lineage (if C. frontosa and B. tricoti form a sister-group, then alanine was substituted by serine in a common ancestor of these species). (b) The informative positions of amino acid substitutions among rhodopsin genes from African cichlids. Read vertically; the top three lines specify the amino acid positions in the consensus sequence of cichlid rhodopsin. A dot indicates identity with the consensus sequence. Shaded columns indicate that identical amino acid substitutions occurred independently in separate lineages

 
All pairwise values of nonsynonymous substitutions per nonsynonymous site (Dn) and synonymous substitutions per synonymous site (Ds) were calculated using version 2.1 of the program MEGA (Kumar et al. 2001Citation ). The ratio Dn/Ds provides an estimate of the evolutionary rate of amino acid substitutions as well as standard errors by bootstrap resampling (Miyata and Yasunaga 1980Citation ; Felsenstein 1985Citation ).

To elucidate the evolution of a particular gene, in general, it is necessary to have prior knowledge of the phylogeny of the species in question. On the basis of the currently accepted phylogeny of cichlid species (fig. 1a, see legend), African cichlid fishes are divided into two groups (fig. 1a ). The riverine group includes a small number of species and thus may serve as an out-group for comparison purposes (Greenwood 1991Citation ; Ribbink 1991Citation ; Mayer, Tichy, and Klein 1998Citation ). The Great Lakes group includes members of the tribes in Lake Tanganyika and the East African riverine Haplochromine, as well as those of the Lake Victoria and Lake Malawi flocks (Mayer, Tichy, and Klein 1998Citation ). For comparison, average Dn/Ds values were calculated for both riverine and Great Lakes cichlid rhodopsin genes. We postulated that the rhodopsin gene has evolved along with the adaptation of each of the cichlid species in the lakes. If so, then amino acid substitutions should have occurred at an accelerated rate. Such a hypothesis would imply that average Dn/Ds values for the Great Lakes group should be higher than those for the riverine group because the adaptation for each niche occurred more frequently in the Great Lakes lineage. But if the rhodopsin gene did not play a role in adaptation in the Great Lakes cichlids, then the Dn/Ds values should be similar between the two groups. Figure 1a presents the cichlid phylogeny used in this study, as well as average Dn/Ds values within each group. Average values for the rhodopsin gene were estimated as 0.176 ± 0.028 for the riverine group and 0.915 ± 0.058 for the Great Lakes group, indicating that the average value for the Great Lakes group is fivefold higher than that for the riverine group. The confidence intervals of the estimates for the two groups did not overlap, and the difference between the estimated Dn/Ds values for the two groups was statistically significant, as estimated by bootstrap resampling (Felsenstein 1985Citation ). The higher average Dn/Ds value for the Great Lakes group suggests that the amino acids encoded by the rhodopsin gene underwent substitutions at a relatively accelerated rate in this group. Thus, there appears to be a correlation between the adaptation to each niche in the Great Lakes and the high rate of change in amino acids encoded by the rhodopsin gene. We also calculated average Dn/Ds values for the rhodopsin gene within a subset of the Great Lakes group, rapid speciated group, consisting of the tribes of Tropheini and Haplochromini (Poll 1986Citation ) in Lake Tanganyika as well as flocks in Lake Victoria and Lake Malawi (designated the THMV group; see fig. 1a ). The average Dn/Ds value for the THMV group was only slightly higher than that for the Great Lakes group, and the difference was not significant (fig. 1a ). The average Dn/Ds value for the THMV group was expected to be significantly higher; however, there were no synonymous substitutions in the four sequence pairs within the THMV group and thus the Dn/Ds values for these sequence pairs could not be included in the analysis.

There are two reasonable explanations for the accelerated evolution of the rhodopsin gene, namely, positive selection or relaxation of negative selection. To discern these two alternatives, we performed a one-tailed Z-test (n = 1,000 replicates) and determined the level of significance of positive selection (Dn/Ds > 1 for each pairwise value in the Great Lakes group) using version 2.1 of the program MEGA (Kumar et al. 2001Citation ). Among the pairwise Dn/Ds values from the Great Lakes group, 25% were >1 (data not shown). Table 1 shows that 16 of the pairwise Dn/Ds values were >2, and seven of these were statistically significant (P < 0.05). These results suggest that the amino acids encoded by the cichlid rhodopsin gene have evolved at an accelerated rate by positive selection, at least in several lineages, during the adaptive radiation of cichlids in the Great Lakes.


View this table:
[in this window]
[in a new window]
 
Table 1 Positive Selection of Rhodopsin

 
The fact that there is evolutionary convergence of amino acids encoded by the rhodopsin gene represents further evidence that positive selection occurred. At 19 amino acid positions (fig. 1b, shaded), substitutions of the same amino acid occurred in separate lineages, independently and repeatedly (for example, at position 292 in fig. 1a; see below). These substitutions were assumed to result from functional convergence.

These findings raise the question of how variation relates to adaptation among Great Lakes cichlid fishes. Molecular adaptation of the rhodopsin gene has been reported in certain fish (Hunt et al. 1996Citation ; Yokoyama et al. 1999Citation ). The Comoran coelacanth (living fossil) lives at depths of ~200 m and perceives only a narrow range of color around 480 nm. Two amino acid substitutions in the Comoran coelacanth proteins RH1 (rhodopsin) and RH2 (RH1-like opsin) elicit a shift in peak wavelength absorbance ({lambda}max) of ~20 nm toward the blue portion of the spectrum, relative to corresponding orthologous pigments ({lambda}max ~ 500 nm; Yokoyama et al. 1999Citation ). Another example is the Cottoid fishes that live at various depths in Lake Baikal, the deepest lake in eastern Siberia. The {lambda}max of rhodopsin in Cottoid fish is shifted in a stepwise manner from 516 nm in species that live in shallow water to about 484 nm in species living in deep water (Hunt et al. 1996Citation ). This shifting of {lambda}max is mediated by amino acid substitutions at four positions (Hunt et al. 1996Citation ).

Among the amino acid substitutions in cichlid rhodopsins (fig. 1b ), it is presumed that position 292 might shift the absorption spectra. This position corresponds to position 308 in human red and mouse green pigments, and the amino acid change from alanine to serine (A308S) shifts the {lambda}max of each pigment toward the blue portion of the spectrum by ~18 nm (Sun, Macke, and Nathans 1997Citation ; Yokoyama 2000Citation ). Therefore, it is reasonable to expect that the cichlid rhodopsin A292S substitution would produce a similar shift (Yokoyama and Radlwimmer 1998Citation ). This substitution occurs in Trematocara macrostoma, Cyphotilapia frontosa, and Benthochromis tricoti (table 1 and fig. 1a ). The Trematocara species live in clear deepwater habitats during the daytime, and the fact that they typically have large eyes suggests adaptation under low illumination so as to improve visual acuity (Coulter 1991Citation ). Benthochromis tricoti has been caught more commonly near the deeper end of its range, and C. frontosa prefers a rocky bottom habitat in deep water, being caught most abundantly between 60 and 120 m in Lake Tanganyika (Coulter 1991Citation ). Therefore, the evolutionary convergence of the A292S substitution in cichlid rhodopsin might be due to adaptation for clear deepwater habitat (not muddy) via a blueshift of the {lambda}max of the rod pigment.

Differences in photic conditions with respect to water color may also have influenced the extent of adaptive evolution of cichlid rhodopsin. In contrast to the waters of Lake Malawi and Lake Tanganyika, those of Lake Victoria are generally quite turbid and, therefore, largely impenetrable to shorter-wavelength light. Color discrimination under these conditions therefore shift to longer wavelengths of the visible spectrum (van der Meer and Bowmaker 1995Citation ). As shown in figure 1b, amino acid substitutions at five positions were specific to Lake Victoria. These substitutions may be a consequence of adaptation to water color. The significance of these five substitutions as well as most of the others observed in cichlid rhodopsins remains unknown at present, and an analysis of potential substitution-dependent shifts in absorption in these rhodopsins is now in progress.

The adaptive changes that are apparent in cichlid rhodopsin may help to explain the adaptation of cichlid fishes to the various niches in the East African Great Lakes. Analysis of the genes that play an important role in the cichlid visual system provides insight into the mechanisms of species diversification at the molecular level.

Acknowledgements

We thank Drs. Y. Shichida and H. Imai (Graduate School of Science, Kyoto University, Kyoto, Japan) for their expertise with opsins and for helpful discussions and Dr. T. Sato (Enviromental Education WWF Coral Reef Conservation and Research Centre) for helpful discussions. This work was supported by the Ministry of Education, Culture, Sports, Science, and Technology of Japan.

Footnotes

Dan Graur, Reviewing Editor

Keywords: cichlid adaptive radiation visual pigment gene rhodopsin positive selection Back

Address for correspondence and reprints: Norihiro Okada, Graduate School of Bioscience and Biotechnology, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8501, Japan. nokada{at}bio.titech.ac.jp Back

References

    Carleton K. L., F. I. Harosi, T. D. Kocher, 2000 Visual pigments of African cichlid fishes: evidence for ultraviolet vision from microspectrophotometry and DNA sequences Vision Res 40:879-890[ISI][Medline]

    Carleton K. L., T. D. Kocher, 2001 Cone opsin genes of African cichlid fishes: tuning spectral sensitivity by differential gene expression Mol. Biol. Evol 18:1540-1550[Abstract/Free Full Text]

    Coulter G. W., 1991 The benthic fish community Pp. 151–199 in G. W. Coulter, ed. Lake Tanganyika and its life. Oxford University Press, London

    Felsenstein J., 1985 Confidence limits on phylogenies: an approach using the bootstrap Evolution 39:783-791[ISI]

    Fryer G., T. D. Iles, 1972 The cichlid fishes of the Great Lakes of Africa Oliver & Boyd, Edinburgh

    Greenwood P. H., 1984 African cichlids and evolutionary theories Pp. 141–154 in A. A. Echelle and I. Kornfield, eds. Evolution of fish species flock. University of Maine at Orono Press, Orono, Me

    ———. 1991 Speciation Pp. 86–102 in M. H. A. Keenleyside, ed. Cichlid fishes: behavior, ecology and evolution. Chapman & Hall, London

    Hunt D. M., J. Fitzgibbon, J. S. Sergey, J. K. Bowmaker, 1996 Spectral tuning and molecular evolution of rod visual pigments in the species flock of cottoid fish in Lake Baikal Vision Res 36:1217-1224[Medline]

    Johnson T. C., C. A. Scholz, M. R. Talbot, K. Kelts, R. D. Ricketts, G. Ngobi, K. Beuning, I. Ssemmanda, J. W. McGill, 1996 Late Pleistocene desiccation of Lake Victoria and rapid evolution of cichlid fishes Science 273:1091-1093[Abstract]

    Kochendoerfer G. G., S. W. Lin, T. P. Sakmar, R. A. Mathies, 1999 How color visual pigments are tuned Trends Biochem. Sci 24:300-305[ISI][Medline]

    Kocher T. D., J. A. Conroy, K. R. McKaye, J. R. Stauffer, 1993 Similar morphologies of cichlid fish in Lakes Tanganyika and Malawi are due to convergence Mol. Phylogenet. Evol 2:158-165[Medline]

    Kumar S., K. Tamura, I. B. Jakobsen, M. Nei, 2001 MEGA 2: molecular evolutionary genetic analysis software Arizona State University, Tempe, Ariz

    Mayer W. E., H. Tichy, J. Klein, 1998 Phylogeny of African cichlid fishes as revealed by molecular markers Heredity 80:702-714[ISI][Medline]

    Meyer A., T. D. Kocher, P. Basasibwaki, A. C. Wilson, 1990 Monophyletic origin of Lake Victoria cichlid fishes suggested by mitochondrial DNA sequences Nature 347:550-553[ISI][Medline]

    Miyata T., T. J. Yasunaga, 1980 Molecular evolution of mRNA: a method for estimating evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences and its application J. Mol. Evol 16:23-36[ISI][Medline]

    Nishida M., 1991 Lake Tanganyika as an evolutionary reservoir of old lineages of East African cichlid fishes: inference from allozyme data Experientia 47:974-979[ISI]

    Poll M., 1986 Classification des Cichlidae du lac Tanganyika: tribus, genres et especes Acad. R. Belg. (2E serie) XLV:1-163

    Ribbink A. J., 1991 Distribution and ecology of the cichlids of the African Great Lakes Pp. 36–59 in M. H. A. Keenleyside, ed. Cichlid fishes: behavior, ecology and evolution. Chapman & Hall, London

    Seehausen O., J. J. M. van Alphen, 1998 The effect of male coloration on female mate choice in closely related Lake Victoria cichlids (Haplochromis nyererei complex) Behav. Ecol. Sociobiol 42:1-8[ISI]

    Seehausen O., J. J. M. van Alphen, F. Witte, 1997 Cichlid fish diversity threatened by eutrophication that curbs sexual selection Science 277:1808-1811[Abstract/Free Full Text]

    Shyue S. K., D. Hewett-Emmett, H. G. Sperling, D. M. Hunt, J. K. Bowmaker, J. D. Mollon, W. H. Li, 1995 Adaptive evolution of color vision genes in higher primates Science 269:1265-1267[ISI][Medline]

    Sturmbauer C., A. Meyer, 1992 Genetic divergence, speciation and morphological status in a lineage of African cichlid fishes Nature 358:578-581[ISI][Medline]

    ———. 1993 Mitochondrial phylogeny of the endemic mouthbrooding lineages from Lake Tanganyika in East Africa Mol. Biol. Evol 10:751-768[Abstract]

    Sun H., J. P. Macke, J. Nathans, 1997 Mechanisms of spectral tuning in the mouse green cone pigment Proc. Natl. Acad. Sci. USA 94:8860-8865[Abstract/Free Full Text]

    Takahashi K., Y. Terai, M. Nishida, N. Okada, 2001 Phylogenetic relationships and ancient incomplete lineage sorting among cichlid fishes in Lake Tanganyika as revealed by analysis of the insertion of retroposons Mol. Biol. Evol 18:2057-2066[Abstract/Free Full Text]

    van der Meer H. J., J. K. Bowmaker, 1995 Interspecific variation of photoreceptors in four co-existing Haplochromine cichlid fishes Brain Behav. Evol 45:232-240[ISI][Medline]

    Wald G., 1968 The molecular basis of visual excitation Nature 219:800-807[ISI][Medline]

    Yokoyama R., S. Yokoyama, 1990 Convergent evolution of the red- and green-like visual pigment genes in fish, Astyanax fasciatus, and human Proc. Natl. Acad. Sci. USA 87:9315-9318[Abstract]

    Yokoyama S., 2000 Molecular evolution of vertebrate visual pigments Prog. Retin. Eye Res 19:385-419[ISI][Medline]

    Yokoyama S., F. B. Radlwimmer, 1998 The "five-sites" rule and the evolution of red and green color vision in mammals Mol. Biol. Evol 15:560-577[Abstract]

    Yokoyama S., F. B. Radlwimmer, N. S. Blow, 2000 Ultraviolet pigments in birds evolved from violet pigments by a single amino acid change Proc. Natl. Acad. Sci. USA 97:7366-7371[Abstract/Free Full Text]

    Yokoyama S., R. Yokoyama, 1996 Adaptive evolution of photoreceptors and visual pigments in vertebrates Annu. Rev. Ecol. Syst 27:543-567[ISI]

    Yokoyama S., H. Zhang, F. B. Radlwimmer, N. S. Blow, 1999 Adaptive evolution of color vision of the Comoran coelacanth (Latimeria chalumnae) Proc. Natl. Acad. Sci. USA 96:6279-6284[Abstract/Free Full Text]

Accepted for publication May 30, 2002.