- Research article
- Open Access
Functional diversity in the color vision of cichlid fishes
BMC Biologyvolume 8, Article number: 133 (2010)
Color vision plays a critical role in visual behavior. An animal's capacity for color vision rests on the presence of differentially sensitive cone photoreceptors. Spectral sensitivity is a measure of the visual responsiveness of these cones at different light wavelengths. Four classes of cone pigments have been identified in vertebrates, but in teleost fishes, opsin genes have undergone gene duplication events and thus can produce a larger number of spectrally distinct cone pigments. In this study, we examine the question of large-scale variation in color vision with respect to individual, sex and species that may result from differential expression of cone pigments. Cichlid fishes are an excellent model system for examining variation in spectral sensitivity because they have seven distinct cone opsin genes that are differentially expressed.
To examine the variation in the number of cones that participate in cichlid spectral sensitivity, we used whole organism electrophysiology, opsin gene expression and empirical modeling. Examination of over 100 spectral sensitivity curves from 34 individuals of three species revealed that (1) spectral sensitivity of individual cichlids was based on different subsets of four or five cone pigments, (2) spectral sensitivity was shaped by multiple cone interactions and (3) spectral sensitivity differed between species and correlated with foraging mode and the spectral reflectance of conspecifics. Our data also suggest that there may be significant differences in opsin gene expression between the sexes.
Our study describes complex opponent and nonopponent cone interactions that represent the requisite neural processing for color vision. We present the first comprehensive evidence for pentachromatic color vision in vertebrates, which offers the potential for extraordinary spectral discrimination capabilities. We show that opsin gene expression in cichlids, and possibly also spectral sensitivity, may be sex-dependent. We argue that females and males sample their visual environment differently, providing a neural basis for sexually dimorphic visual behaviour. The diversification of spectral sensitivity likely contributes to sensory adaptations that enhance the contrast of transparent prey and the detection of optical signals from conspecifics, suggesting a role for both natural and sexual selection in tuning color vision.
Vision is central to the survival of animals. Visual cues are used for orientation, detecting prey, avoiding predators and communication. The visual process starts with visual pigments absorbing light and initiating a photochemical cascade that leads to neural signaling, perception and ultimately visually mediated behavior. A common method for studying vision is the measurement of spectral sensitivity. Spectral sensitivity is the relative efficiency of detection of light as a function of wavelength. Spectral sensitivity is used to describe the characteristics of visual pigments found in cone photoreceptors in the retina, and it is particularly useful in describing the mechanisms of color vision [1, 2].
Comparative studies in vertebrates established the presence of four spectrally distinct classes of cone visual pigments produced by several different opsin genes: SWS1 being the ultraviolet/violet-sensitive class (355-440 nm), SWS2 being the short wavelength-sensitive class (410-490 nm), RH2 being the mid-wavelength-sensitive class (470-530 nm) and LWS being the long-wavelength-sensitive class (495-570 nm) [3, 4]. These opsin genes have been identified in the earliest vertebrate lineage, the jawless fishes or Agnathans (lamprey) , and arose through duplications of a single ancestral opsin gene. Mutations in duplicated genes can lead to the evolution of additional spectral classes of cone visual pigments within a class of opsin gene. In fact, in many teleost fishes, cone opsin genes have undergone gene duplication to produce a wide range of opsin genes [6–8]. The expression of different subsets of these genes can potentially generate large-scale variation in spectral sensitivity and the mechanisms of color vision in fish [9–11].
The cichlid fishes of Lake Malawi are an excellent model system for examining large-scale variation in color vision. Lake Malawi has 700-800 species of cichlids [12–14] that have evolved from a common ancestor in a brief period of evolutionary time (2-4 million years) . Malawi cichlids are notable for their diversity in male nuptial color patterns, sexual dimorphism in color patterns and visual communication processes governing mate choice [16–18]. Cichlids have undergone multiple opsin gene duplication events, producing seven classes of cone opsin genes [10, 19]. These seven opsin genes include an ultraviolet-sensitive (SWS1), two short-wavelength-sensitive (SWS2a and SWS2b), three mid-wavelength-sensitive (RH2aα, RH2aβ, RH2b) and long-wavelength-sensitive (LWS) cone opsins. Of particular interest is that Malawi cichlids show differential expression of primarily three of the seven available cone opsin genes [9, 20, 21], with some evidence of differential expression through ontogeny . This raises an intriguing question of how seven opsin genes are maintained in Malawi cichlids. Our central research focus concerns the large-scale diversity in cichlid visual systems and differences in the number of cone classes that participate in color vision. Here we examine variation in cone classes contributing to spectral sensitivity in cichlids between individuals, sexes and species and discuss the adaptive significance of this variation in color vision.
To quantify the number of physiologically functional cone classes, we recorded electroretinograms (ERG) from whole fish preparations of three Malawi cichlids: Metriaclima zebra, Melanochromis auratus and Protomelas taeniolatus. We measured the b-wave amplitude of ERGs, representative of ON (response to light stimulus onset) bipolar cell activity, to assess the integrative response of functional cones. We were primarily interested in evaluating the number of cone classes contributing to spectral sensitivity and their corresponding visual pigments, but we also used multiple-mechanism modeling to understand the opponent and nonopponent cone interactions at play that can shape the spectral tuning of cone mechanisms. Our analysis was focused on describing variation in spectral sensitivity with respect to individual, sex and species differences.
Individual differences in spectral sensitivity
Over 100 spectral sensitivity curves were recorded from a total of 15 females and 19 males of the three species studied. The spectral sensitivity of each individual was evaluated under three different color background conditions, each aimed at isolating the sensitivity of specific cones (Figure 1; see Additional file 1 for the background conditions and the quantum catches of the various cone mechanisms). For each individual, we identified the number of sensitivity peaks and their spectral locations for the three background conditions. We then fitted cone pigment absorption templates  to the sensitivity peaks using a least-squares method.
Individuals of all three species of Malawi cichlids examined in this study possessed 11 different cone subsets (Table 1; see Table 2 for specific λmax values and goodness of fit of visual pigment templates). A maximum of five different cone classes were present in 41% of fish, four different cone classes were present in 56% of fish, and three different cone classes were present in only 3% of the fish (one individual) (Table 1). All cone subsets (excluding no. 5) consisted of at least two pigments corresponding to single cones and two pigments corresponding to double cones.
Examining the individuals of each species separately revealed that all individuals in M. zebra possessed four cone pigments, while 55% of individuals in M. auratus and 67% of individuals in P. taeniolatus possessed five cone pigments (Table 1).
The collection of all cone subsets exhibited by either M. auratus or P. taeniolatus included the complete opsin gene set (Table 1). That is, each of these species utilizes the complete set of six cone pigments (RH2aα and RH2aβ pooled because of genetic and functional similarity [10, 19]). However, the collection of all subsets exhibited by M. zebra encompassed only five cone pigments and excluded the LWS pigment.
Sex differences in spectral sensitivity
To evaluate how sex differences contribute to the variation in cone subsets, we calculated the frequency of cone subsets in each species and sex. In all species, the individual variation in pigment subsets used was high, and a strong trend for differences between females and males in the identity and frequency of the pigment subsets they utilized emerged (Table 3). The degree to which pigment subsets were shared between sexes varied across species. Half of M. zebra females shared subset 2 with all conspecific males, 28% of M. auratus males shared subset 3 with half the females and 60% of P. taeniolatus females shared subsets 6 and 7 with 71% of males.
The frequency of each cone pigment present in females and males was calculated (Figure 2). Sex differences in the frequency of cone pigments were evident in all species, with M. auratus and M. zebra showing the largest differences. The frequency of M. auratus females showing SWS1 pigments was 75% larger than in males, while the frequency of males showing SWS2a and LWS pigments was 71% and 100% larger than in females, respectively. The frequency of M. zebra males showing SWS2a cone pigment was 100% larger than in females, while SWS2b cone pigment was found only in females. In contrast, the frequency of cone pigments in P. taeniolatus showed less variation across sexes, with pigment frequencies differing 12%-42% between sexes.
We also examined the effect of sex on the pattern of cone opsin gene expression. We focused our analysis on M. auratus and P. taeniolatus, the species that exhibited the lowest and highest degree of common cone pigment subsets between the sexes, respectively. Variation in opsin gene expression between females and males in M. auratus was larger than in P. taeniolatus (Figure 3). Specifically, in M. auratus, the expression of RH2b and RH2a opsin genes differed across sexes (t-test, df = 8, P < 0.0002 and P < 0.005, respectively). In contrast, in P. taeniolatus, no significant differences in the relative opsin gene expression were detected between sexes (t-test, P > 0.2 for all genes). See Additional file 2 for detailed t-test results and Additional file 3 for primer specifications for cone opsin genes. The frequency of cone pigments and opsin gene expression are not quantitatively comparable. In the calculation of pigment frequency, each individual is scored as either 1 (pigment present) or 0 (pigment absent), whereas gene expression can assume any value between 0 and 1. Moreover, the expression level of an opsin gene does not necessarily dictate its contribution to sensitivity, since patterns of convergence of cones onto retinal interneurons and the network processing of these neurons  are critical to governing visual responsiveness. Our results show that females and males differ in the cone opsin gene expression profiles. Our results also suggest that females and males differ in cone pigment subsets and in the frequency of cone pigments. However, owing to sample size limitation imposed by the complexity of spectral sensitivity measurements, the statistical significance of the sex differences observed in cone pigment subsets and in the frequency of cone pigments cannot be evaluated at this time. It is important to note, however, that the variability of all three properties of the visual system of cichlids was consistently larger in M. auratus than in the other two species.
Species differences in spectral sensitivity
In general, different cichlid species used different cone subsets (Table 3). Closer examination reveals that females of the three species did not use the same cone subsets, but 43% of male M. auratus shared three cone subsets with P. taeniolatus males.
The frequency of cone pigments varied across species (Figure 4a). The frequency of RH2a (SD = 14%) and RH2b (SD = 5%) showed the smallest variation and did not differ significantly between species (Fisher's exact test; P = 0.091 and P = 1, respectively). In contrast, the frequency of LWS (SD = 53%) and SWS1 (SD = 38%) showed the largest variation and differed significantly between species (Fisher's exact test; P < 0.001 for both pigments). See Figure 4 caption for detailed statistics. The frequency of SWS1 was highest in M. zebra (100%), lower in M. auratus (54%) and lowest in P. taeniolatus (25%). The frequency of LWS in P. taeniolatus (100%) and M. auratus (81%) was higher than in M. zebra (0%).
The frequency of the SWS1 cone pigment across species qualitatively correlated with planktivory. M. zebra may be highly planktivorous, feeding on both phytoplankton (mainly diatoms) and zooplankton (Cladocerans and copepods) . M. auratus is omnivorous, consuming algae and a wide range of nonalgal dietary items such as cyclopoid copepods [24, 25]. P. taeniolatus is almost strictly herbivorous and feeds on the biocover detached from rocks, mainly comprising algae and diatoms [21, 26].
The frequency of the LWS cone pigment across species was quantitatively correlated with the proportion of long-wavelength reflectance in the color pattern of conspecific males (Figure 4b). As a first-order approximation, the quantum catch of an LWS cone pigment was calculated for the spectral reflectance of the species used in this study. The normalized quantum catch was highest in P. taeniolatus (100%), lower in M. auratus (80%) and lowest in M. zebra (17%).
Color vision in Lake Malawi cichlids
Spectral sensitivity curves determined using different color background conditions were used to evaluate mechanisms of color vision. Color vision requires the possession of at least two differentially sensitive cones that interact through opponent and nonopponent processes to enable wavelength discrimination . Our results show that 97% of individuals possessed four or five cone classes. To study cone interactions, we used a multiple-cone mechanism (MCM) model, employing strict criteria (see Methods), to reconstruct spectral sensitivity under different background conditions. To do this, we selected three individual fish that exhibited the most frequent cone subset for a given species. These were subsets 2, 3 and 7, occurring in 73%, 36% and 42% of the M. zebra, M. auratus and P. taeniolatus individuals, respectively.
Spectral sensitivity was successfully reconstructed using the MCM model (R2 = 0.83-0.99 for all species and background conditions) (Figure 5). Cone interaction was expressed as an additive (nonopponent) or subtractive (opponent) contribution to the modeled spectral sensitivity, and we summarize these cone interactions in Table 4 (see Additional file 4 for more comprehensive information that details the specific coefficients and goodness of fit of the model). Typically, the cone class exhibiting the highest sensitivity for a spectral range contributed positively to the modeled spectral sensitivity, whereas other cone classes contributed either positively (nonopponent) or negatively (opponent) to the modeled spectral sensitivity. The successful reconstruction of spectral sensitivity curves further validated our analysis and the distribution of the cone pigment subsets. Most importantly, we identified and characterized cone interactions representative of retinal neural network processing, essential to color vision.
Pentachromatic color vision in Lake Malawi cichlids
Ninety-seven percent of fish examined in this study employed cone subsets of four or five pigments, with most M. auratus and P. taeniolatus exhibiting five functional cone classes. The spectral sensitivity of all three species of cichlids was shaped by both opponent and nonopponent cone interactions, and thus our data provide the first comprehensive evidence for pentachromatic color vision in vertebrates. Although the presence of five (or more) cone classes was previously reported in microspectrophotometry studies [10, 28], this is the first demonstration that five cone mechanisms were used by a single individual. In this regard, designating the dimensionality of color vision requires the essential evidence that shows not only that the different cone mechanisms participate in spectral sensitivity but also that the different cone mechanisms show multiple neuronal interactions. P. taeniolatus and M. auratus showed five distinct cone classes spectrally compressed in a range of 146 nm (subset 7) and 187 nm (subset 3), respectively. The spectral compression of this sensitivity provides the basis for extraordinary spectral discrimination, which could be important in visual communication. To better understand how pentachromatic color vision would be beneficial to cichlids, we need to consider factors at play in the visual environment of cichlids where visual communication takes place: (1) the high species diversity of cichlids in Lake Malawi (700-800 species) [14, 25, 29], (2) the large within-genera variation in the color pattern of fish and the small within-genera variation in their size and shape [30–32] and (3) the prevalence of visual cues driving mate choice [33, 34] and intrasexual competition . Visual communication would depend on exquisite sensory performance in this visual environment, since cichlids are continually challenged to make fine-scale discriminations of complex color patterns. This visual system capability would be particularly important in guiding visual behavior related to mate choice, where the detection and recognition of conspecific optical signals is critical. These visual adaptations could be fundamental to mate choice fidelity and contribute to the maintenance of biodiversity in the cichlid communities of Lake Malawi.
It is important to point out that our findings differ from the notion that Lake Malawi cichlids primarily exhibit three cone pigments [10, 36]. Recent reports show a significant variation in the number of opsin genes expressed among Malawi species , with several species (one of which is P. taeniolatus) expressing more than three cone opsin genes . Furthermore, microspectrophotometry studies have revealed classes of cones in low abundance in the retina of various species that are of unknown significance [10, 36]. Our study suggests that these rare cone classes contribute to the spectral sensitivity of Lake Malawi cichlids.
Sex differences in the color vision of Lake Malawi cichlids
Our study provides some evidence for sex differences in opsin gene expression, and possibly in spectral sensitivity, for cichlids. While the sample size in electrophysiology studies tends to be low, a possible pattern of differentiation between males and females emerged. The large diversity of cone pigment subsets found in each species was in part qualitatively related to sex differences. This, together with qualitative sex differences in pigment frequency and significant differences in cone opsin gene expression, suggests the possibility for sex differences in the visual system of some Lake Malawi cichlids. Our results suggest that the most prominent sex differences occurred in M. auratus. We argue here that males and females may have quite different visual requirements in visual communication scenarios [38–40]. M. auratus differs from the two other species in two aspects. (1) M. zebra and P. taeniolatus males hold a permanent territory while females choose among potential mates ; M. auratus males, on the other hand, assume a territory only temporarily around breeding events . (2) M. zebra and P. taeniolatus males are conspicuous while the color pattern of females is dull, whereas both female and male M. auratus are conspicuous and exhibit colorful, but very different, body patterns .
We propose a possible explanation for the observed distribution of sex differences across species. Considering the significant role of visual cues in female mate choice [33, 34] and male-male competition for territory , both females and territorial males of M. zebra and P. taeniolatus would likely benefit from possessing color vision attuned to provide the best recognition ability of conspecific males. In contrast, M. auratus males assume a territory only temporarily and thus are less dependent on recognition of conspecific males to increase their mating success. Additionally, not holding a permanent territory potentially allows M. auratus males to choose between females. Indeed, in the cichlid Astatotilapia flaviijosephi, males defend territories only during the breeding season and actively choose between females . Furthermore, the investment of M. auratus females in their color pattern adds support to the hypothesis that M. auratus males choose between females. Thus, the distribution of sex differences in the visual system across species might be associated with territorial behavior. However, the adaptive role of sexual dimorphism in color vision of Lake Malawi cichlids requires further examination.
Adaptive significance of species differences in the color vision of Lake Malawi cichlids
The variation in cone pigment frequency across species reflects visual adaptations that both enhance the contrast of transparent prey and detect signals important for mate choice, suggesting that both natural selection and sexual selection played a role in shaping the spectral sensitivity of cichlids.
The frequency of the SWS1 cone pigment across species qualitatively correlated with planktivory. UV photoreception mediated by the SWS1 cone pigment enhances the contrast of transparent zooplankton against the water background and thus aids in their detection [43, 44]. However, UV sensitivity does not confer any advantage over longer-wavelength sensitivity in the detection of opaque food items such as loose algae. Thus, UV sensitivity is expected to be highly advantageous for the planktivorous M. zebra, less advantageous for the omnivorous M. auratus and least advantageous for the herbivorous P. taeniolatus, which is in agreement with our findings. Thus, the presence of the SWS1 pigment in the retina of cichlids has likely significant adaptive value. In this regard, the expression of the SWS1 opsin gene was recently reported to be highest among species foraging on zooplankton, phytoplankton and algae and lowest among species foraging on fish or benthic invertebrates ; however, the relationship between opsin gene expression and spectral sensitivity in cichlids is still largely unknown. UV photoreception is also tightly linked to polarization sensitivity, and this may provide yet greater contrast enhancement of transparent prey [45–47]. Researchers in our lab are currently examining polarization sensitivity in cichlid fishes and its role in visual behavior.
The frequency of the LWS cone pigment across species was quantitatively correlated with the proportion of long-wavelength reflectance in the color pattern of conspecific males. While the possession of the LWS cone pigment by M. auratus and P. taeniolatus may facilitate the detection of long-wavelength signals of conspecifics, this would not be the case for M. zebra that do not have prominent long-wavelength signals, which is in agreement with our findings. Therefore, the presence of the LWS pigment in M. auratus and P. taeniolatus has a significant adaptive value. Interestingly, the detection of UV reflections from the color pattern of these fish species may be mediated through the β-band of the LWS cone pigment (Figure 4b), thus potentially eliminating the need for the possession of the SWS1 pigment to detect UV signals from conspecifics.
Both natural and sexual selection work sequentially to contribute to the divergence of Lake Malawi rock-dwelling cichlids . In the first episode of cladogenesis, competition for trophic resources resulted in the differentiation of trophic morphology and the diversification of the visual system to allow the utilization of different foraging styles. The correlation of the frequency of the SWS1 cone pigment across species to planktivory (this study), and the recent report that the expression of the SWS1 opsin gene correlates with planktivory , support the idea of diversification of the visual system based on competition for trophic resources. In the second episode of cladogenesis, sexual selection contributed to the differentiation of male nuptial coloration and the accompanying diversification of the visual system to allow high recognition ability of conspecifics. The correlation of the frequency of the LWS cone pigment across species to the proportion of long-wavelength reflectance from male color patterns supports the idea of visual system diversification based on the diverse male nuptial coloration.
Different cichlid species used different cone subsets, with the spectral sensitivity of females being more constrained than that of males. This is suggestive of specific tuning of color vision in females to allow for the best recognition ability of the distinctly colored conspecific males. Previous studies have shown that different cichlid species express diverse opsin gene subsets [9, 21], supporting our results. However, a close examination of the different constraints imposed on the visual systems of females and males needs to be undertaken.
The retention of cone opsin genes in Lake Malawi cichlids
None of the individuals examined possessed all available six cone pigments (RH2aα and RH2aβ pooled). However, the large number of pigments utilized and the high diversity in pigment subsets within a species add support for the retention of the complete set of cone opsin genes in Malawi cichlids.
Several factors may contribute to the observed pigment diversity within a species. First, a significant genetic component of opsin expression was reported  where individuals with different genetic histories may utilize different pigment subsets. Second, thyroid hormones (TH) are important in vertebrate development and can modulate the opsin gene expression and induce cone loss in fish [49, 50]. The level of TH in fish has been associated with environmental stress and subordinance via the effect of cortisol [51, 52]. Thus, individuals differing in their stress level or social status may utilize diverse pigment subsets through TH level modulation. Third, social visual cues in fish modulate the activity of GnRH neurons . Since GnRH has been shown to affect retinal neurons , it is possible that it may be responsible for changes in opsin gene expression and spectral sensitivity in fish.
The cichlid model system illustrates that the visual systems of fish may differ across individuals, sexes and species. The large number of available cone opsin genes facilitates this variation in the spectral sensitivity of fish. We show that even rare cone populations and opsin gene expression at low levels contribute to the spectral sensitivity of fish. Our results suggest that the diversification of color vision across species contributes to sensory adaptations that both enhance the contrast of transparent prey and the detection of optical signals of conspecifics. Therefore, both natural and sexual selection may work in concert to shape spectral sensitivity in fish. Taken together, our findings have important implications for understanding the variable nature of fish color vision and the selective forces shaping detection and recognition capabilities.
Animals and holding conditions
Three species of Lake Malawi cichlids were used in this study: Metriaclima zebra (blue top manda), Melanochromis auratus and Protomelas taeniolatus (Old World Exotic Fish, Homestead, FL, USA). M. zebra and M. auratus belong to a rock-dwelling evolutionary lineage, the mbuna, whereas P. taeniolatus is a member of the nonmbuna lineage. These species are sexually dimorphic and occur sympatrically in the rocky habitat in Lake Malawi. The males have distinctive nuptial color patterns, and they use the rocky habitat differently with respect to reproductive behavior [25, 29]. Adult fish were held in our aquatic facility tanks under a 12 h:12 h light-dark photoperiod at 25 ± 1°C. Facility lighting featured enhanced full spectrum fluorescent lamps (UV-Blue actinic and BlueMax lamps; Full Spectrum Solutions, Jackson, MI, USA). All experimental and animal care procedures were approved by Queen's University Animal Care Committee under the auspices of the Canadian Council for Animal Care.
Electroretinogram (ERG) experimental apparatus
The general design of the optical system and recording apparatus has been described previously [2, 54]. Two background channels each with a 250-W quartz-halogen lamp (Osram, Danvers, MA, USA) were used to provide constant background, to which test fish were light-adapted. Long- and short-pass interference filters (Fused Silica, OD 2; Edmund Optics, Barrington, NJ, USA), band-pass interference filters (Edmund Optics), broadband color filters (Schott, Elmsford, NY, USA) and reflective neutral density filters (Edmund Optics), were used to manipulate the spectral irradiance of each background channel. Light from the two background channels was guided to the electrophysiology rig using a bifurcated optical fiber (fused silica, numerical aperture (NA) = 0.22; Fiberoptic Systems, Simi Valley, CA, USA).
The stimulus channel used a 300-W xenon arc lamp system (Thermo Oriel, Stratford, CT, USA). The optical path consisted of a monochrometer (Instruments SA, Metuchen, New Jersey, USA), Inconel quartz neutral density wedge (0-4.0 neutral density; CVI Melles Griot, Albuquerque, NM, USA), shutter (Uniblitz, Rochester, NY, USA), optical filters to block spectral sidebands and UV lenses to match the numerical aperture of the liquid light pipe (fused silica; NA = 0.55; Fiberoptic Systems). The background and stimulus optical fibers were fitted to a beam splitter to produce a stimulus beam (diameter 0.5 cm at the plane of the fish eye) contained within the background beam (diameter 1 cm). This setup ensured that the chromatically adapted portion of the fish retina was also the one stimulated.
Preparation of fish
Prior to ERG recordings, fish were immersed in a solution of 150 mg l-1 tricaine methanesulfonate (MS-222) until the fish reached stage IV anesthesia . Intramuscular injections of a general anaesthetic, Maranil (0.1 mg g-1 body mass) and an immobilizing agent, pancuronium bromide (0.04 mg g-1 body mass), were administered at several sites. Test fish were then placed in a holding cradle in a Faraday cage. Experimental fish were irrigated with aerated fresh water (20°C, flow rate ~3 ml s-1), and their body was covered with a moist cloth.
ERG recording procedure
ERG recordings started at least 1 hour following the onset of the light phase and concluded before the onset of the dark phase to avoid any effects related to circadian rhythm [2, 56, 57]. A glass electrode (1.5 mm outer diameter, 1 mm inner diameter, borosilicate glass; World Precision Instruments, Sarasota, FL, USA) that was pulled to a tip diameter of 80-125 μm (P-97 Flaming/Brown Micropipette puller; Sutter Instruments, Novato, CA, USA) was loaded with saline (0.684 M sodium chloride) and inserted into a saline-filled chlorided AgCl half-cell (A-M systems, Sequim, WA, USA). The electrode tip was positioned using a micromanipulator on the dorsal-nasal surface of the right eye. A ground electrode was attached to the caudal fin and a chlorided-silver reference electrode was placed on the head of the test fish. Fish were chromatically adapted for 1 hour prior to experiments. The stimulus duration was 500 ms with an interstimulus interval of 5 s. An isolated bioamplifier (ISO-80; World Precision Instruments, Sarasota, FL, USA) amplified the ERG signal and filtered the signal using band-pass filter settings (5-Hz low-pass, 100-Hz high-pass). The amplified signal was analyzed with a 16-bit A/D data acquisition system and Signal 4.0 software (Micro 1401; Cambridge Electronic Design Limited, Cambridge, UK). A custom-designed software analysis module determined the b-wave amplitude (corresponding to the ON response of bipolar cells) by measuring the potential change between the a-wave and b-wave peaks. Spectral sensitivity was measured in 20-nm increments, from 320 to 700 nm, using a staggered wavelength presentation to prevent adaptation to a spectral region.
Analysis of electroretinograms
A response versus intensity (RI) curve was generated for each wavelength examined. To determine the sensitivity at a given wavelength, the empirical Naka-Rushton function was least-squares fitted to the RI curve with a slope parameter of 1 [58, 59]. The log irradiance level corresponding to half response (LogK) was determined, from which the sensitivity was calculated by taking the reciprocal of this value. Typically, the upper asymptote of the Naka-Rushton function could be reached for all wavelengths. There were a few exceptions at wavelengths shorter than 360 nm and longer than 600 nm. At these wavelengths, the value of the maximum response parameter (Rmax) was set to that of the neighboring wavelength examined, prior to fitting. A log relative sensitivity curve was constructed by normalizing the log absolute sensitivity values to the maximum sensitivity. This procedure was repeated for each individual under each of the background conditions.
Background isolation conditions
Four background conditions were designed and used for isolating cones dominating certain spectral regions. A long-wavelength-isolation condition (LW) was used for isolating cone mechanisms most sensitive to long wavelengths, dim and intense short-wavelength-isolation conditions (Dim-SW and SW) aimed at isolating cone mechanisms most sensitive to short wavelengths and a relatively spectrally flat background condition was used as a control. To design and carefully control the level of light adaptation of the respective cone mechanisms, a quantum catch model was used:
where Q i (photons cm-2 s-1) denotes the quantum catch of cone mechanism i (i = 1, 2, ..., 6), A i (λ) represents the visual pigment absorption coefficient of cone mechanism i at a wavelength λ (nm), and E i (λ) denotes the photon irradiance of the background light field at a wavelength λ.
Characterization of background conditions
The irradiance provided under the various background conditions was characterized by measuring the spectral irradiance of the background beam. Spectral irradiance was measured using a spectroradiometer (QE65000; Ocean Optics, Dunedin, FL, USA) connected to a 2-m optical fiber (QP600-2-UV/VIS; Ocean Optics) that was fitted with a cosine corrector (CC-3-UV; Ocean Optics). The spectroradiometer utilized a 1024 × 58-element square silicon CCD (charge-coupled device) array and was configured with a 25-μm slit and a variable blaze wavelength grating (HC-1, groove density = 300 mm-1; Ocean Optics), resulting in an effective spectral resolution of 1.9 nm (FWHM) between 200 and 950 nm. The fiber end was held approximately 5 cm away from the emergence plane of the background beam, ensuring sampling of the entire beam diameter. The spectroradiometer setup was calibrated for absolute irradiance using a NIST (National Institute of Standards and Technology, Gaithersberg, MD, USA) calibrated Halogen-Deuterium dual light source (200-1000 nm, DH-2000-CAL; Ocean Optics). In cases where spectral variation exceeded the dynamic range of the spectroradiometer and thus a reliable measurement could not be obtained, the irradiance delivered under the background condition was calculated by multiplying the measured output of each light source by the spectral transmission of the filters used in producing the background condition.
Analysis of spectral sensitivity curves
1.Identification of functional cones
Two criteria were devised for the identification of cone mechanisms. A sensitivity peak was considered to be a cone mechanism if it satisfied the following conditions: (1) the sensitivity peak appeared under at least two background isolation conditions, and (2) the sensitivity peak exhibited changes that were directly related to changes in the spectral composition of the background conditions. For instance, a sensitivity peak in the ultraviolet range (UV; 340-400 nm) was identified as a UV cone when sensitivity increased under a short-wavelength isolation background or as the β-band of mid- and long-wavelength cones when sensitivity decreased under a short-wavelength isolation background.
2.Correspondence between visual pigments and cone mechanisms
To relate visual pigments with cone mechanisms, absorption templates of visual pigments were fitted to spectral sensitivity curves using a least-squares fit. This technique was previously used in numerous studies and repeatedly has been shown to allow for the identification of the cone classes in the retinas of animals [60–62]. In this regard, interactions between different cone classes produce sensitivity peaks that are narrower than the absorption templates of visual pigments, resulting in a somewhat reduced goodness of fit [60, 61]. Visual pigment absorbance templates  were constructed for the opsin genes previously reported in cichlids: Note that the λmax of A1-reconstituted visual pigment is provided in parentheses SWS1 (368 nm), SWS2b (423 nm), SWS2a (456 nm), Rh2b (484 nm), Rh2aα (519 nm), Rh2aβ (528 nm) and LWS (560 nm) [10, 19]. In the case of Rh2aα and Rh2aβ, the spectral overlap of the two visual pigment curves necessitated calculating an average λmax of 523 nm. The 532-nm λmax absorption spectrum was used for subsequent analysis and is hereafter referred to as RH2a. To generate visual pigment templates we combined absorption spectra for the A1 and A2 chromophores (A1, equations 1, 2, 4, 5a, 5b; A2, equations 1, 4, 6a, 6b, 8a, 8b) . The proportion of the A2 state was presented using a fraction parameter, a (0 ≤ a ≤ 1), and therefore the absorption spectra of a given cone type exhibiting an A2 proportion of a was calculated as A2 (a) = A1·(1 - a) + A2·a. The λmax of each cone type exhibits a defined wavelength shift as the A2 proportion changes . This shift was also taken into account when generating the visual pigment templates for varying A2 proportions.
Typically, the long-wavelength peak in spectral sensitivity was located between 560 and 600 nm, and thus it could correspond to either the RH2a (523-560 nm; λmax A1-λmax A2) or the LWS (560-626 nm; λmax A1-λmax A2) visual pigments. Therefore, to correlate the long-wavelength peak with a visual pigment absorption curve, it was necessary to determine the A2 proportion in the retina of each individual. Visual pigment absorption templates for the RH2a and LWS pigments were fitted to the long-wavelength limb of the spectral sensitivity curve. The least-squares fit was performed while leaving the A2 proportion and a magnitude coefficient unrestricted, allowing the software to find a visual pigment template with an A2 proportion that best describes the long-wavelength limb. Estimating the A2 proportion on the basis of the long-wavelength limb is the best possible approach because (1) the λmax shift between the A1 and A2 states is largest for long-wavelength-sensitive cones, and (2) the spectral sensitivity is based on one cone class and thus is not affected by opponent interactions. For each individual, the A2 proportion estimate was used for fitting all visual pigment templates to the spectral sensitivity curve. All absorption templates were corrected for the spectral transmission of the fish lens (P. taeniolatus ; M. zebra and M. auratus measured ).
3.Modeling cone interactions
A multiple-cone mechanism (MCM) model was used to determine the relative contribution of each cone class to each spectral sensitivity curve by assigning weights, which can be positive (excitatory) or negative (inhibitory) for each cone mechanism. This "upper envelope" model that was previously used in primates [65, 66] and fish  assumes that the spectral sensitivity of the eye is determined by the cones that are most sensitive over the spectral region of concern. This linear model takes the general form Sλ1, λ2 (λ) = ∑ k i ·A i (λ), where S λ1, λ2 (λ) denotes the modeled spectral sensitivity in the spectral range enclosed by λ1 and λ2 while accounting for cone interactions; A i (λ) represents the absorption of the visual pigment template, corresponding to cone i (i = 1, 2, ..., 6) at a wavelength λ (nm); and k i denotes the weight representing the contribution of cone i. The MCM model was least-squares fitted to the spectral sensitivity curves determined under the various background conditions for each individual separately. The spectral range between two sensitivity notches was set as the λ1 - λ2 range [61, 65, 66]. For modeling spectral sensitivity, we selected the cone interaction that (1) included the minimum possible number of cones, (2) allowed for the best fit under all background conditions, and (3) ensured the same type of cone interaction (opponent and nonopponent) under all background conditions. To prevent overparameterization of the model, the number of possible cone interactions was restricted to three. That is, the maximum number of free parameters in the model was three, while the number of data points used to fit the model was always larger than three. The weights assigned to each cone, in addition to the R2 values (total amount of variance accounted for by the model across the spectrum), were determined. Log relative sensitivity was transformed to percent relative sensitivity and normalized to range between 0 and 1 prior to fitting .
Preparation of retinal samples
Upon completion of ERG recordings, fish were dark-adapted for 1 hour and then killed by cervical transection. Under deep red illumination (>650 nm), both eyes were enucleated and hemisected along an anterior-posterior axis. The neural retina was then dissected free of pigmented epithelium. Each isolated retina was preserved in 0.5 ml of RNAlater (Ambion), stored at 4°C for the first 24 hours, and then kept at -80°C until further processing. The sex of the fish was determined using the acetocarmine technique [68, 69]. Fish were dissected and their gonads were removed, stained with acetocarmine and viewed under a dissecting microscope to determine the fish sex. All fish were sexually mature adults showing well-developed eggs or sperm and were at least 9 months old. Cichlids typically reach sexual maturity at 6 months from release [70, 71].
Relative gene expression by qPCR
Quantitative real-time RT-PCR (qPCR) was used to quantify the relative levels of mRNA expression corresponding to the various cone opsin genes. Unless specified otherwise, all procedures described below were performed following the manufacturer's protocols. Total RNA was extracted from retinas (Absolutely RNA Miniprep Kit; Agilent Technologies, Santa Clara, CA, USA), and its amount and quality were determined by reading the absorbance at 260 nm and calculating the absolute 260/280 absorbance ratio using a benchtop spectrophotometer (Cary 300 Bio; Varian, Santa Clara, CA, USA). Total cDNA was synthesized using 250 ng of total retinal RNA and a Superscript III first-strand synthesis SuperMix (Invitrogen, Carlsbad, CA, USA) at 50°C for 30 minutes. See Additional file 3 for primer sequences for the amplification of the cDNA opsin genes in cichlids (RH2aα and RH2aβ combined). All primers were analyzed using the Primer-Blast primer design tool (National Center for Biotechnology Information, Bethesda, MD, USA) for product size, melting temperature, guanine-cytosine content (%) and sequence specificity. The specificity of all primer pairs was tested by amplifying target sequences present within the cDNA synthesized from retinal total RNA. Amplification was carried out in a Mastercycler gradient 5331 (Eppendorf) using the following concentrations: 2.5 mM MgCl2, 200 μM deoxynucleotide triphosphate (dNTP), 100 μM forward and reverse primers, 0.05 U GotTaq Flexi DNA polymerase (Promega, Madison, WI, USA), 1× PCR buffer and 1 μl of cDNA template, in a final volume of 25 μl. The PCR cycling profile consisted of an initial hold at 92°C for 2 minutes, 40 cycles (92°C for 25 s, 55°C for 25 s, 72°C for 25 s), and a final extension at 72°C for 5 min. To verify that the amplified product obtained from each primer pair consisted of a single band and was of the correct size, amplified DNA products were separated in a 2.0% agarose gel 1× TBE (Tris/Borate/EDTA) buffer, visualized by GelRed on a gel imager (AlphaImager, Cell Biosciences, Santa Clara, CA, USA) and compared with a benchtop 100-bp DNA ladder (Promega).
qPCR analysis of individual retinal cDNA samples was carried out for each opsin gene using Brilliant SYBR green qPCR Master Mix (Agilent Technologies) in a real-time quantitative system MX3000P™ (Agilent Technologies). Each 25-μl reaction contained 1× Brilliant SYBR Green Master Mix, 200 mM quantities of both forward and reverse primers, 0.05 μl of ROX passive reference dye and 1 μl of four times diluted cDNA. The reaction mix was placed in 96-well nonskirted polypropylene PCR plate and capped with optical strip caps (Agilent Technologies). The plate was briefly centrifuged and eye-inspected for the presence of bubbles. The thermocycle program was 95°C for 10 minutes, followed by 40 cycles of 95°C for 25 seconds, 55°C for 25 seconds and 72°C for 25 seconds. Controls included a reaction lacking cDNA template) and a nontranscribed reaction for genomic DNA contamination (No-RT). All samples were run in triplicate, and fractional cycle values (Cq) were averaged. Amplification efficiencies per sample per target were calculated from the slope of the amplification curve in the exponential phase using LinRegPCR 11.0 software . Relative gene expression was determined for the six opsin genes as a fraction of the total cone opsin genes expressed for an individual :
where T g,i /T all,i is the relative gene expression ratio for gene g normalized by the total cone opsin genes expressed for individual i, E g,i is the amplification efficiency for gene g and individual i, and Cq g,i is the fractional cycle value for gene g and individual i. Finally, relative gene expression was averaged for each species and sex.
Underwater irradiance measurements
Underwater spectral irradiance at a horizontal line of sight was measured in July 2008 at a near-shore site at Cape Maclear on the northwestern part of the Nankumba Peninsula, Lake Malawi (14° 1" 0" S 34° 51" 0" E; local time: 11:00-11:30, GMT +02:00, average solar zenith angle within the water = 25.6°). The sampling site, Otter Point, is exposed to wind and wave action and exhibits a rocky bottom that subsides with depth and becomes a sandy bottom at ~15 m depth . Irradiance was measured using a thermoelectrically cooled spectroradiometer (QE65000; Ocean Optics) connected to a 30-m optical fiber (ZPK600-30-UV/VIS; Ocean Optics) that was fitted with a cosine corrector (CC-3-UV; Ocean Optics). See the specifications of the spectroradiometer in the Characterization of background conditions section. Holding the irradiance probe, a SCUBA diver attained position at depths of 5, 10 and 15 m, and readings were saved on a computer placed on a boat (positioned as far as possible and never between the diver and the sun to prevent shading). The spectroradiometer setup was calibrated for absolute irradiance following the same procedure described in the Characterization of background conditions section.
Spectral reflectance measurements
Spectral reflectance of the color pattern of fish was measured using a spectroradiometer (USB2000; Ocean Optics) connected to one arm of a 2-m bifurcated optical fiber (BIF600-2-UV/VIS; Ocean Optics). The other arm of the fiber was connected to a light source (DH-2000-BAL; Ocean Optics). The spectroradiometer utilized a 2048-element linear silicon CCD array and was configured with a 50-μm slit and a grating (groove density, 600 mm-1; blaze wavelength, 400 nm; grating 2, Ocean Optics), resulting in an effective spectral resolution of 2.06 nm (FWHM) between 200 and 950 nm. The light source integrated two lamps, tungsten-halogen and deuterium, providing a high and spectrally balanced output between 200 and 1000 nm. The common end of the bifurcated fiber was fitted with a flat black reflectance probe exhibiting a 3-mm diameter tip cut at an angle of 45°. Prior to each measurement session, the light source was allowed to warm up for at least 40 minutes. Then a measurement of a Spectralon diffuse reflectance standard (WS-1-SL; Ocean Optics) was taken as 100% reflectance, and a dark measurement was taken as zero reflectance. Fish were immersed in a 1:10 clove oil:ethanol solution immediately after capture until the fish reached stage IV anesthesia . Fish were held submerged in lake water, and spectral reflectance was measured at 10 points along the trunk of the fish (one individual per species). Care was given to sealing the reflectance probe against the fish skin to reduce stray light. Readings were acquired and saved on a laptop computer.
To calculate the quantum catch of the LWS cone pigment, we used the spectral radiance reflected from the color pattern of the species studied (equation 1). In this case, however, Ei(λ) was substituted by the product of the sidewelling irradiance at a depth of 5 m, Eh(λ) times the spectral reflectance of the color pattern of the fish, R(λ). Absorption template for the LWS cone pigment was generated while assuming A2 proportion of 0.5 .
To study the effect of sex on the relative opsin gene expression, we performed multiple t-tests where an experiment-wise error rate of 5% was corrected to 0.83% (α = 0.05/6 = 0.0083) following Bonferroni correction for six hypothesis tests  that correspond to six opsin genes (RH2aα and RH2aβ were pooled because of genetic and functional similarity [10, 19]). Prior to performing statistical analyses, the normality of all data was confirmed using the Kolmogorov-Smirnov test, and the homogeneity of variance across treatment groups was confirmed using Cochran's C test . To study the effect of species on the frequency of cone pigments, we performed multiple Fisher's exact tests (two-tailed) with a corrected experimental error (α = 0.0083). Statistical analysis was performed using the Statistica (t-test) and R version 2.11.1 (Fisher's exact test) software.
dim short-wavelength-isolation background condition
full width half max
short-wavelength-isolation background condition
long-wavelength-isolation background condition
long wavelength sensitive
National Institute of standards and technology
response versus intensity
short wavelength sensitive
Lisney TJ, Studd E, Hawryshyn CW: Electrophysiological assessment of spectral sensitivity in adult Nile tilapia Oreochromis niloticus: evidence for violet sensitivity. J Exp Biol. 2010, 213: 1453-1463. 10.1242/jeb.036897.
Parkyn DC, Hawryshyn CW: Spectral and ultraviolet-polarisation sensitivity in juvenile salmonids: A comparative analysis using electrophysiology. J Exp Biol. 2000, 203: 1173-1191.
Yokoyama S, Yokoyama R: Adaptive evolution of photoreceptors and visual pigments in vertebrates. Annu Rev Ecol Syst. 1996, 27: 543-567. 10.1146/annurev.ecolsys.27.1.543.
Bowmaker JK: Evolution of vertebrate visual pigments. Vision Res. 2008, 48: 2022-2041. 10.1016/j.visres.2008.03.025.
Collin SP, Knight MA, Davies WL, Potter IC, Hunt DM, Trezise AEO: Ancient colour vision: Multiple opsin genes in the ancestral vertebrates. Curr Biol. 2003, 13: R864-R865. 10.1016/j.cub.2003.10.044.
Christoffels A, Koh EGL, Chia JM, Brenner S, Aparicio S, Venkatesh B: Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes. Mol Biol Evol. 2004, 21: 1146-1151. 10.1093/molbev/msh114.
Furutani-Seiki M, Wittbrodt J: Medaka and zebrafish, an evolutionary twin study. Mech Dev. 2004, 121: 629-637. 10.1016/j.mod.2004.05.010.
Meyer A, Schartl M: Gene and genome duplications in vertebrates: the one-to-four (-to-eight in fish) rule and the evolution of novel gene functions. Curr Opin Cell Biol. 1999, 11: 699-704. 10.1016/S0955-0674(99)00039-3.
Carleton KL, Kocher TD: Cone opsin genes of African cichlid fishes: Tuning spectral sensitivity by differential gene expression. Mol Biol Evol. 2001, 18: 1540-1550.
Parry JWL, Carleton KL, Spady T, Carboo A, Hunt DM, Bowmaker JK: Mix and match color vision: Tuning spectral sensitivity by differential opsin gene expression in Lake Malawi cichlids. Curr Biol. 2005, 15: 1734-1739. 10.1016/j.cub.2005.08.010.
Shand J, Davies WL, Thomas N, Balmer L, Cowing JA, Pointer M, Carvalho LS, Trezise AEO, Collin SP, Beazley LD, Hunt DM: The influence of ontogeny and light environment on the expression of visual pigment opsins in the retina of the black bream, Acanthopagrus butcheri. J Exp Biol. 2008, 211: 1495-1503. 10.1242/jeb.012047.
Fryer G, Iles TD: The Cichlid Fishes of the Great Lakes of Africa: Their Biology and Evolution. 1972, Neptune City, NJ: T. F. H. Publications
Salzburger W: The interaction of sexually and naturally selected traits in the adaptive radiations of cichlid fishes. Mol Ecol. 2009, 18: 169-185. 10.1111/j.1365-294X.2008.03981.x.
Turner GF, Seehausen O, Knight ME, Allender CJ, Robinson RL: How many species of cichlid fishes are there in African lakes?. Mol Ecol. 2001, 10: 793-806. 10.1046/j.1365-294x.2001.01200.x.
Genner MJ, Seehausen O, Lunt DH, Joyce DA, Shaw PW, Carvalho GR, Turner GF: Age of cichlids: new dates for ancient lake fish radiations. Mol Biol Evol. 2007, 24: 1269-1282. 10.1093/molbev/msm050.
Danley PD, Kocher TD: Speciation in rapidly diverging systems: lessons from Lake Malawi. Mol Ecol. 2001, 10: 1075-1086. 10.1046/j.1365-294X.2001.01283.x.
Kocher TD: Adaptive evolution and explosive speciation: the cichlid fish model. Nature Reviews Genetics. 2004, 5: 288-298. 10.1038/nrg1316.
Jordan R, Kellogg K, Juanes F, Stauffer J: Evaluation of female mate choice cues in a group of Lake Malawi mbuna (Cichlidae). Copeia. 2003, 181-186. 10.1643/0045-8511(2003)003[0181:EOFMCC]2.0.CO;2.
Spady TC, Parry JWL, Robinson PR, Hunt DM, Bowmaker JK, Carleton KL: Evolution of the cichlid visual palette through ontogenetic subfunctionalization of the opsin gene arrays. Mol Biol Evol. 2006, 23: 1538-1547. 10.1093/molbev/msl014.
Carleton KL, Spady TC, Streelman JT, Kidd MR, McFarland WN, Loew ER: Visual sensitivities tuned by heterochronic shifts in opsin gene expression. BMC Biol. 2008, 6: 22-10.1186/1741-7007-6-22.
Hofmann CM, O'Quin KE, Marshall NJ, Cronin TW, Seehausen O, Carleton KL: The eyes have it: regulatory and structural changes both underlie cichlid visual pigment diversity. PLoS Biol. 2009, 7: e1000266-10.1371/journal.pbio.1000266.
Govardovskii VI, Fyhrquist N, Reuter T, Kuzmin DG, Donner K: In search of the visual pigment template. Vis Neurosci. 2000, 17: 509-528. 10.1017/S0952523800174036.
Djamgoz MBA, Yamada M: Electrophysiological characteristics of retinal neurons; interactions and functional output. The Visual System of Fish. Edited by: Douglas R, Djamgoz MBA. 1990, London: Chapman and Hall, 161-210.
Reinthal PN: The feeding-habits of group of herbivorous rock-dwelling cichlid fishes (Cichlidae, Perciformes) from Lake Malawi, Africa. Environ Biol Fishes. 1990, 27: 215-233. 10.1007/BF00001674.
Ribbink AJ, Marsh BA, Marsh AC, Ribbink AC, Sharp BJ: A preliminary survey of the cichlid fishes of rocky habitats in Lake Malawi. S Afr J Zool. 1983, 18: 149-310.
Konings A: Konings's Book of Cichlids and All the Other Fishes of Lake Malawi. 1990, Neptune, NJ: T.F.H publications
Hurvich LM: Color Vision. 1981, Sunderland, MA: Sinauer Associates
Fuller RC, Fleishman LJ, Leal M, Travis J, Loew E: Intraspecific variation in retinal cone distribution in the bluefin killifish, Lucania goodei. J Comp Physiol A. 2003, 189: 609-616. 10.1007/s00359-003-0435-x.
Konings A: Lake Malawi Cichlids in their Natural Habitat. 1995, Germany: Cichlid Press
Deutsch JC: Colour diversification in Malawi cichlids: Evidence for adaptation, reinforcement or sexual selection?. Biol J Linn Soc Lond. 1997, 62: 1-14. 10.1111/j.1095-8312.1997.tb01614.x.
McElroy DM, Kornfield I, Everett J: Coloration in African cichlids: diversity and constraints in Lake Malawi endemics. Neth J Zool. 1991, 41: 250-268. 10.1163/156854291X00180.
Stauffer JR, Bowers NJ, Kellogg KA, McKaye KR: A revision of the blue-black Pseudotropheus zebra (Teleostei: Cichlidae) complex from Lake Malawi, Africa, with a description of a new genus and ten new species. Proc Acad Nat Sci Phila. 1997, 148: 189-230.
Couldridge VCK, Alexander GJ: Color patterns and species recognition in four closely related species of Lake Malawi cichlid. Behav Ecol. 2002, 13: 59-64. 10.1093/beheco/13.1.59.
Kidd MR, Danley PD, Kocher TD: A direct assay of female choice in cichlids: all the eggs in one basket. J Fish Biol. 2006, 68: 373-384. 10.1111/j.0022-1112.2006.00896.x.
Dijkstra PD, Seehausen O, Groothuis TGG: Direct male-male competition can facilitate invasion of new colour types in Lake Victoria cichlids. Behav Ecol Sociobiol. 2005, 58: 136-143. 10.1007/s00265-005-0919-5.
Jordan R, Kellogg K, Howe D, Juanes F, Stauffer J, Loew E: Photopigment spectral absorbance of Lake Malawi cichlids. J Fish Biol. 2006, 68: 1291-1299. 10.1111/j.0022-1112.2006.00992.x.
Carleton KL: Cichlid fish visual systems: mechanisms of spectral tuning. Integr Zool. 2009, 4: 75-86. 10.1111/j.1749-4877.2008.00137.x.
Darwin C: The Descent of Man and Selection in Relation to Sex. 1871, London: Murray
Shine R: Ecological causes for the evolution of sexual dimorphism: a review of the evidence. Q Rev Biol. 1989, 64: 419-461. 10.1086/416458.
Slatkin M: Ecological causes of sexual dimorphism. Evolution. 1984, 38: 622-630. 10.2307/2408711.
Barlow GW: Mating systems among cichlid fishes. Cichlid Fishes: Behaviour, Ecology and Evolution. Edited by: Keenleyside MHA. 1991, London: Chapman & Hall, 173-190.
Werner NY, Lotem A: Choosy males in a haplochromine cichlid: First experimental evidence for male mate choice in a lekking species. Anim Behav. 2003, 66: 293-298. 10.1006/anbe.2003.2208.
Browman HI, Novales Flamarique I, Hawryshyn CW: Ultraviolet photoreception contributes to prey search behavior in 2 species of zooplanktivorous fishes. J Exp Biol. 1994, 186: 187-198.
Johnsen S, Widder EA: Ultraviolet absorption in transparent zooplankton and its implications for depth distribution and visual predation. Mar Biol. 2001, 138: 717-730. 10.1007/s002270000499.
Novales Flamarique I, Browman HI: Foraging and prey-search behaviour of small juvenile rainbow trout (Oncorhynchus mykiss) under polarized light. J Exp Biol. 2001, 204: 2415-2422.
Sabbah S, Shashar N: Polarization contrast of zooplankton: a model for polarization-based sighting distance. Vision Res. 2006, 46: 444-456. 10.1016/j.visres.2005.05.017.
Shashar N, Hanlon RT, Petz AD: Polarization vision helps detect transparent prey. Nature. 1998, 393: 222-223. 10.1038/30380.
Carleton KL, Hofmann CM, Klisz C, Patel Z, Chircus LM, Simenauer LH, Soodoo N, Albertson RC, Ser JR: Genetic basis of differential opsin gene expression in cichlid fishes. J Evol Biol. 2010, 23: 840-853. 10.1111/j.1420-9101.2010.01954.x.
Temple SE, Ramsden SD, Haimberger TJ, Veldhoen KM, Veldhoen NJ, Carter NL, Roth WM, Hawryshyn CW: Effects of exogenous thyroid hormones on visual pigment composition in Coho salmon (Oncorhynchus kisutch). J Exp Biol. 2008, 211: 2134-2143. 10.1242/jeb.009365.
Raine JC, Coffin AB, Hawryshyn CW: Systemic thyroid hormone is necessary and sufficient to induce ultraviolet-sensitive cone loss in the juvenile rainbow trout retina. J Exp Biol. 2010, 213: 493-501. 10.1242/jeb.036301.
Ejike C, Schreck CB: Stress and social hierarchy rank in Coho salmon. Trans Am Fish Soc. 1980, 109: 423-426. 10.1577/1548-8659(1980)109<423:SASHRI>2.0.CO;2.
Pickering AD: Growth and stress in fish production. Aquaculture. 1993, 111: 51-63. 10.1016/0044-8486(93)90024-S.
Ramakrishnan S, Wayne NL: Social cues from conspecifics alter electrical activity of gonadotropin-releasing hormone neurons in the terminal nerve via visual signals. Am J Physiol Regul Integr Comp Physiol. 2009, 297: R135-R141.
Hawryshyn CW, Moyer HD, Allison WT, Haimberger TJ, McFarland WN: Multidimensional polarization sensitivity in damselfishes. J Comp Physiol A. 2003, 189: 213-220.
Jolly DW, Mawdesley-Thomas LE, Bucke D: Anaesthesia of fish. Vet Rec. 1972, 91: 424-426. 10.1136/vr.91.18.424.
Li P, Temple S, Gao Y, Haimberger TJ, Hawryshyn CW, Li L: Circadian rhythms of behavioral cone sensitivity and long wavelength opsin mRNA expression: a correlation study in zebrafish. J Exp Biol. 2005, 208: 497-504. 10.1242/jeb.01424.
Parkyn DC, Hawryshyn CW: Polarized-light sensitivity in rainbow-trout (Oncorhynchus-Mykiss): characterization from multiunit responses in the optic-nerve. J Comp Physiol A. 1993, 172: 493-500. 10.1007/BF00213531.
Dowling JE, Ripps H: Adaptation in skate photoreceptors. J Gen Physiol. 1972, 60: 698-719. 10.1085/jgp.60.6.698.
Naka KI, Rushton WAH: S-potentials from colour units in retina of fish (Cyprinidae). J Physiol (London). 1966, 185: 536-555.
Hawryshyn CW, Ramsden SD, Betke KM, Sabbah S: Spectral and polarization sensitivity of juvenile Atlantic salmon (Salmo salar): phylogenetic considerations. J Exp Biol. 2010, 213: 3187-3197. 10.1242/jeb.038760.
Hughes A, Saszik S, Bilotta J, DeMarco PJ, Patterson WF: Cone contributions to the photopic spectral sensitivity of the zebrafish ERG. Vis Neurosci. 1998, 15: 1029-1037. 10.1017/S095252389815602X.
Patterson WF, McDowell AL, Hughes A, Bilotta J: Opponent and nonopponent contributions to the zebrafish electroretinogram using heterochromatic flicker photometry. J Comp Physiol A. 2002, 188: 283-293. 10.1007/s00359-002-0302-1.
Harosi FI: An analysis of two spectral properties of vertebrate visual pigments. Vision Res. 1994, 34: 1359-1367. 10.1016/0042-6989(94)90134-1.
Hofmann CM, O'Quin KE, Marshall NJ, Carleton KL: The relationship between lens transmission and opsin gene expression in cichlids from Lake Malawi. Vision Res. 2010, 50: 357-363. 10.1016/j.visres.2009.12.004.
Kalloniatis M, Harwerth RS: Spectral sensitivity and adaptation characteristics of cone mechanisms under white-light adaptation. J Opt Soc Am A. 1990, 7: 1912-1928. 10.1364/JOSAA.7.001912.
Sperling HG, Harwerth RS: Red-green cone interactions in increment-threshold spectral sensitivity of primates. Science. 1971, 172: 180-184. 10.1126/science.172.3979.180.
Demarco PJ, Powers MK: Spectral sensitivity of On and Off responses from the optic-nerve of goldfish. Visual Neurosci. 1991, 6: 207-217. 10.1017/S0952523800006222.
Guerrero RD, Shelton WL: An aceto-carmine squash technique for sexing juvenile fishes. The Progressive Fish Culturist. 1974, 36: 56-10.1577/1548-8659(1974)36[56:AASMFS]2.0.CO;2.
Wassermann GJ, Afonso LOB: Validation of the aceto-caramin technique for evaluating phenotypic sex in Nile tilapia (Oreochromis niloticus) fry. Ciência Rural. 2002, 32: 133-139. 10.1590/S0103-84782002000100023.
Balon EK: Early ontogeny of Labeotropheus Ahl 1927 (Mbuna, Cichlidae, Lake Malawi) with a discussion on advanced protective styles in fish reproduction and development. Environ Biol Fishes. 1977, 2: 147-176. 10.1007/BF00005370.
Morrison CM, Miyake T, Wright JR: Histological study of the development of the embryo and early larva of Oreochromis niloticus (Pisces: Cichlidae). J Morphol. 2001, 247: 172-195. 10.1002/1097-4687(200102)247:2<172::AID-JMOR1011>3.0.CO;2-H.
Ruijter JM, Ramakers C, Hoogaars WMH, Karlen Y, Bakker O, van den Hoff MJB, Moorman AFM: Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009, 37: e45-10.1093/nar/gkp045.
Quinn GP, Keough MJ: Experimental Design and Data Analysis for Biologists. 2002, New York: Cambridge University Press
Underwood AJ: Experiments in Ecology. 1997, Cambridge, UK: Cambridge University Press
We thank Rotem Goldshmid, Prof. Barrie Frost, Dr. Nicholas Roberts, Dr. Shelby Temple, Shyh-Chi Chen, Ilse Dedden, and Dr. Taro Matsumoto for valuable comments. We thank Dr. James Fraser and Dr. Tom Lisney for help in optimizing the ERG experimental apparatus and procedures, Dr. William Nelson and Dr. Nikolaus Troje for suggestions regarding modeling procedures, Dr. Chris Hofmann for providing lens transmission data for cichlids, and Dr. Karen Carleton for qPCR primer sequences and helpful suggestions for gene expression analysis. This research was supported by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant, NSERC Research Tools and Instruments Grant, the Canada Research Chair Program, the Canada Foundation for Innovation, and the Ontario Innovation Trust to CWH. SS was supported by an NSERC Vanier Canada Graduate Scholarship. SMG was supported by a NSERC Postdoctoral Fellowship.
The authors declare that they have no competing interests.
SS carried out the design of the experiments, executed the spectral sensitivity measurements, performed all data analysis and participated in writing the manuscript; CWH participated in the design of the study and writing the manuscript; RLL performed the qPCR gene expression profiling; SMG performed the spectral reflectance measurements. All authors read and approved the final manuscript.