Positive selection and gene conversion in SPP120, a fertilization-related gene, during the East African cichlid fish radiation

The ability to infer historical natural selection from sequence data aides in finding genes that might be important in adaptation and the formation of new species. As the fastest evolving and largest known vertebrate radiation, the cichlid fish of the African Great Lakes exhibit a wide range of recent morphological diversification. We used DNA databases, mostly of expressed sequence tags, to find candidate orthologous coding sequences from 2 tribes of cichlids and, using an automated procedure, scanned these sequence pairs for high dN/dS, the signal of positive selection and protein adaptation. The results included vertebrate genes commonly found to be under selection (e.g., major histocompatibility complex [MHC] loci) as well as genes known to be important specifically in the cichlid radiation (e.g., long-wave-sensitive opsins). Further investigation focused on a gene encoding a fertilization-related protein, SPP120, which was previously known only from cichlids. Using maximum likelihood analysis on novel SPP120 cDNA sequences from a range of African cichlids, we demonstrate the influence of positive selection in a specific subregion of the protein. We also show that SPP120 is a tandemly arranged, multicopy gene evolving with occasional interlocus gene conversion. A search of the Medaka genome database also revealed a tandem arrangement of multiple SPP120 copies and evolutionary rate differences between Medaka gene subregions mirroring those found for cichlids. Combined, these results suggest that SPP120 has been under repeated diversifying selection for over 100 Myr.     

Multilocus phylogeny of cichlid fishes (Pisces: Perciformes): evolutionary comparison of microsatellite and single-copy nuclear loci

Among vertebrates, cichlid fishes are the paradigmatic example of adaptive radiation and ecological specialization. In turn, molecular genetic studies of cichlids have focused primarily on more recently diverged groups. Here, we present an evolutionary hypothesis of the major lineages of cichlid fishes based on DNA sequence data from two nuclear loci. One marker, Tmo-4C4, is a single-copy locus containing a region of amino acid similarity to the muscle protein TITIN. Flanking sequence from a second, microsatellite, locus Tmo-M27, shows similarity to mammalian RAS guanine nucleotide-releasing factor. We compare and combine data from these loci to evaluate phylogenetic performance. In separate and combined analyses, the sequence data support and clarify previous morphological hypotheses of cichlid major-group relationships. Indian and Malagasy cichlids form a basal, paraphyletic group. Neotropical cichlids are the sister clade to an African assemblage composed of the paraphyletic west and Pan-African lineages and a group of east African rift lake taxa. We use a consensus phylogeny of the Cichlidae to trace evolutionary changes in the microsatellite repeat motif at Tmo-M27. Analysis reveals that the repeat region was nearly lost in the ancestor to cichlids and then amplified extensively in African taxa. Results demonstrate that the two new DNA markers could be widely applied in perciform systematics. Furthermore, the comparative approach can unveil mutational dynamics of simple-sequence repeat loci over long periods of fish evolution. Simple-sequence repeat regions are increasingly being found in introns of important regulatory genes. We address issues involving their function and suggest caution in making assumptions of strict neutrality.     

Cone opsin genes of african cichlid fishes: tuning spectral sensitivity by differential gene expression

Spectral tuning of visual pigments is typically accomplished through changes in opsin amino acid sequence. Within a given opsin class, changes at a few key sites control wavelength specificity. To investigate known differences in the visual pigment spectral sensitivity of the Lake Malawi cichlids, Metriaclima zebra (368, 488, and 533 nm) and Dimidiochromis compressiceps (447, 536, and 569 nm), we sequenced cone opsin genes from these species as well as Labeotropheus fuelleborni and Oreochromis niloticus. These cichlids have five distinct classes of cone opsin genes, including two unique SWS-2 genes. Comparisons of the inferred amino acid sequences from the five cone opsin genes of M. zebra, D. compressiceps, and L. fuelleborni show the sequences to be nearly identical. Therefore, evolution of key opsin sites cannot explain the differences in visual pigment sensitivities. Real-time PCR demonstrates that different cichlid species express different subsets of the available cone opsin genes. Metriaclima zebra and L. fuelleborni express a complement of genes which give them UV-shifted visual pigments, while D. compressiceps expresses a different set to produce a red-shifted visual system. Thus, variations in cichlid spectral sensitivity have arisen through evolution of gene regulation, rather than through changes in opsin amino acid sequence.     

Positive Darwinian Selection Drives the Evolution of the Morphology-Related Gene, EPCAM, in Particularly Species-Rich Lineages of African Cichlid Fishes

The study of genetic evolution within the context of adaptive radiations offers insights to genes and selection pressures that result in rapid morphological change. Cichlid fishes are very species-rich and variable in coloration, behavior, and morphology, and so provide a classical model system for studying the genetics of adaptive radiation. In this study, we researched the evolution of the epithelial cell adhesion molecule (EPCAM), a candidate gene for the adaptive evolution of morphology broadly, and skin development specifically, in fishes. We compared EPCAM gene sequences from a rapidly speciating African cichlid lineage (the haplochromines), a species-poor African lineage (Nile tilapia Oreochromis niloticus), and a very young adaptive radiation in the Neotropics (sympatric crater lake Midas cichlids, Amphilophus sp.). Our results, based on a hierarchy of evolutionary analyses of nucleotide substitution, demonstrate that there are different selection pressures on the EPCAM gene among the cichlid lineages. Several waves of positive natural selection were identified not only on the terminal branches, but also on ancestral branches. Interestingly, significant positive or directional selection was found in the haplochromine cichlids only but not the comparatively species-poor tilapia lineage. We hypothesize that the strong signal of selection in the ancestral African cichlid lineage coincided with the transition from riverine to lacustrine habitat. The two neotropical species for which we collected new sequence data were invariant in the EPCAM locus. Our results suggest that functional changes promoted by positive Darwinian selection are widespread in the EPCAM gene during African cichlid evolution.     

池蝶蚌Sox2 基因在不同组织及不同月龄精巢中的表达

Sox 基因家族在胚胎发育过程和性别分化中起重要作用, 为研究池蝶蚌中Sox 基因的功能, 以人SRY基因HMG-box 保守区的序列设计简并引物, 以雌、雄池蝶蚌基因组DNA 和精巢cDNA 为模板进行扩增, 获得了2 个不完全相同的序列, 分别为DNA-HMG1、DNA-HMG2 和cDNA-HMG, 长度均为220 bp, 编码73个氨基酸。与人等物种Sox1、Sox2、Sox3 及Sox14 有很高的同源性, 雌雄个体之间没有序列差异性。采用RACE-PCR 扩增获得了池蝶蚌性腺Sox2 部分cDNA 片段, 长度为1774 bp, 该序列核苷酸与欧洲帽贝的SoxB和人类的Sox2 的同源性最高; 在部分开放阅读框249 个氨基酸残基中, 具有Sox 家族典型的HMG-box 结构域, 与人类、小鼠、原鸡和斑马鱼等Sox2 的HMG-box 同源性为98%。为了解该基因在各组织中的表达情况,采用实时荧光定量PCR 方法分析了外套膜、闭壳肌、鳃、肠、肝、肾、精巢和卵巢在内的8 种组织hs-Sox2的表达情况, 结果显示, hs-Sox2 基因在8 种组织中均有表达, 其中在肾脏中的表达量最高, 其次是肠与闭壳肌, 在雄性性腺中的表达量明显高于雌性性腺, 在肝脏中的表达量最低; 为了解hs-Sox2 在不同性腺发育时期的表达情况, 采用实时荧光定量PCR 方法分析了5 个不同月龄的精巢组织中hs-Sox2 的表达情况, 结果显示在39 月龄性腺的表达量最高, 其次是16 月龄性腺, 63 月龄蚌中的表达量最少。以上结果表明, hs-Sox2 基因可能参与了池蝶蚌精巢的发育及功能的维持。       

Identification of chick rax/rx genes with overlapping patterns of expression during early eye and brain development.

We have isolated chick rax/rx cDNAs, cRaxL (chick Rax/Rx-like) and cRax, (chick Rax) and examined their expression patterns during early eye and brain development. The cRaxL cDNA encodes a 228 amino acid protein that is most closely related to the zebrafish Rx1 and Rx2. The cRax cDNA encodes a 317 amino acid protein, which shares higher homology with the Xenopus Rx. In addition to the homeodomain, the octapeptide and paired tail domains are conserved between the cRax and other vertebrate Rax/Rx, while cRaxL lacks the octapeptide containing N-terminal region which is conserved among all other members of the rax/rx gene family identified so far. The chick rax/rx genes are expressed in overlapping domains in the anterior neural ectoderm which corresponds to the forebrain and retina field, and later in the optic vesicle. cRax mRNA can be detected earlier than cRaxL prior to the formation of the notochord and its expression domain appears broader than that of cRaxL.     

Molecular Characterization of Two Endothelin Pathways in East African Cichlid Fishes

The adaptive radiations of cichlid fishes in East Africa have been associated with the acquisition of evolutionary novelties as well as the ecological opportunities existing in the East African Great lakes. Two remarkable evolutionary innovations are the pharyngeal jaw apparatus, found in all cichlid species, and the anal fin egg-spots of mouthbrooding cichlids. Based on their conserved functions during the development of both the jaw apparatus and pigmentation, the endothelin ligands and receptors form a putative link between these naturally and sexually selected traits. Here we study the evolutionary history of four members of two endothelin pathways (Edn1/EdnrAa and Edn3b/EdnrB1a) to elucidate their possible roles during the evolution and development of key innovations in East African cichlids species. The analyses performed on partial sequences (ca. 6,000 bp per taxon) show that all four endothelin family members evolved under purifying selection, although both ligands are characterized by an accelerated rate of protein evolution in comparison to the receptors. In accordance with earlier findings, we show that the mature protein sequence of Edn1 and Edn3 are highly conserved, also in cichlids, whereas the preproendothelin parts are variable indicating relaxed selective constraints. In the receptors, nonsynonymous substitutions were mainly found in the ligand-binding domains suggesting functional divergence. Gene expression assays with Real-Time PCR indeed reveal that the two studied endothelin pathways are expressed in the cichlid pharyngeal jaw and in the haplochromine egg-spot (among other pigment-cell containing tissues), suggesting their involvement during morphogenesis of naturally and sexually selected traits in cichlids.     

Identification of Cichlid Fishes from Lake Malawi Using Computer Vision

Background

The explosively radiating evolution of cichlid fishes of Lake Malawi has yielded an amazing number of haplochromine species estimated as many as 500 to 800 with a surprising degree of diversity not only in color and stripe pattern but also in the shape of jaw and body among them. As these morphological diversities have been a central subject of adaptive speciation and taxonomic classification, such high diversity could serve as a foundation for automation of species identification of cichlids.

Methodology/Principal Finding

Here we demonstrate a method for automatic classification of the Lake Malawi cichlids based on computer vision and geometric morphometrics. For this end we developed a pipeline that integrates multiple image processing tools to automatically extract informative features of color and stripe patterns from a large set of photographic images of wild cichlids. The extracted information was evaluated by statistical classifiers Support Vector Machine and Random Forests. Both classifiers performed better when body shape information was added to the feature of color and stripe. Besides the coloration and stripe pattern, body shape variables boosted the accuracy of classification by about 10%. The programs were able to classify 594 live cichlid individuals belonging to 12 different classes (species and sexes) with an average accuracy of 78%, contrasting to a mere 42% success rate by human eyes. The variables that contributed most to the accuracy were body height and the hue of the most frequent color.

Conclusions

Computer vision showed a notable performance in extracting information from the color and stripe patterns of Lake Malawi cichlids although the information was not enough for errorless species identification. Our results indicate that there appears an unavoidable difficulty in automatic species identification of cichlid fishes, which may arise from short divergence times and gene flow between closely related species.     

High prevalence of non-synonymous substitutions in mtDNA of cichlid fishes from Lake Victoria

When a population size is reduced, genetic drift may fix slightly deleterious mutations, and an increase in nonsynonymous substitution is expected. It has been suggested that past aridity has seriously affected and decreased the populations of cichlid fishes in Lake Victoria, while geographical studies have shown that the water levels in Lake Tanganyika and Lake Malawi have remained fairly constant. The comparably stable environments in the latter two lakes might have kept the populations of cichlid fishes large enough to remove slightly deleterious mutations. The difference in the stability of cichlid fish population sizes between Lake Victoria and the Lakes Tanganyika and Malawi is expected to have caused differences in the nonsynonymous/synonymous ratio, ω (= d_N/d_S), of the evolutionary rate. Here, we estimated ω and compared it between the cichlids of the three lakes for 13 mitochondrial protein-coding genes using maximum likelihood methods. We found that the lineages of the cichlids in Lake Victoria had a significantly higher ω for several mitochondrial loci. Moreover, positive selection was indicated for several codons in the mtDNA of the Lake Victoria cichlid lineage. Our results indicate that both adaptive and slightly deleterious molecular evolution has taken place in the Lake Victoria cichlids' mtDNA genes, whose nonsynonymous sites are generally conserved.     

Molecular characterization of two isoforms of ZFAND3 cDNA from the Japanese quail and the leopard gecko, and different expression patterns between testis and ovary

Zing finger AN1-type domain 3 (ZFAND3), also known as testis expressed sequence 27 (Tex27), is a gene found in the mouse testis, but its physiological function is unknown. We identified the full-length sequences of two isoforms (short and long) of ZFAND3 cDNA from Japanese quail and leopard gecko. This is the first cloning of avian and reptilian ZFAND3 cDNA. The two isoforms are generated by alternative polyadenylation in the 3'UTR and have the same ORF sequences encoding identical proteins. There were highly conserved regions in the 3'UTR of the long form near the polyadenylation sites from mammals to amphibians, suggesting that the features for determining the stability of mRNA or translation efficiency differ between isoforms. The deduced amino acid sequence of ZFAND3 has two putative zinc finger domains, an A20-like zinc finger domain at the N-terminal and an AN1-like zinc finger domain at the C-terminal. Sequence analysis revealed an additional exon in the genomic structures of the avian and reptilian ZFAND3 genes which is not present in mammals, amphibians, or fish, and this exon produces additional amino acid residues in the A20-like zinc finger domain. Expression analysis in Japanese quail revealed that the expression level of ZFAND3 mRNA was high in not only the testis but also the ovary, and ZFAND3 mRNA was expressed in both spermatides of the testis and oocytes of the ovary. While the short form mRNA was mainly expressed in the testis, the expression level of the long form mRNA was high in the ovary. These results suggest that ZFAND3 has physiological functions related to germ cell maturation and regulatory mechanisms that differ between the testis and ovary.     

Mitochondrial DNA Phylogeny of the Family Cichlidae: Monophyly and Fast Molecular Evolution of the Neotropical Assemblage

A mitochondrial DNA (mtDNA) phylogeny of cichlid fish is presented for the most taxonomically inclusive data set compiled to date (64 taxa). 16S rDNA data establish with confidence relationships among major lineages of cichlids, with a general pattern congruent with previous morphological studies and less inclusive molecular phylogenies based on nuclear genes. Cichlids from Madagascar and India are the most basal groups of the family Cichlidae and sister to African–Neotropical cichlids. The cichlid phylogeny suggests drift-vicariance events, consistent with the fragmentation of Gondwana, to explain current biogeographic distributions. Important phylogenetic findings include the placement of the controversial genus Heterochromis basal among African cichlids, the South American genus Retroculus as the most basal taxon of the Neotropical cichlid assemblage, and the close relationship of the Neotropical genera Cichla with Astronotus rather than with the crenicichlines. Based on a large number of South American genera, the Neotropical cichlids are defined as a monophyletic assemblage and shown to harbor significantly higher levels of genetic variation than their African counterparts. Relative rate tests suggest that Neotropical cichlids have experienced accelerated rates of molecular evolution. But these high evolutionary rates were significantly higher among geophagine cichlids. Received: 18 September 1998 / Accepted: 16 December 1998     

Age and spread of the haplochromine cichlid fishes in Africa

The Haplochromini are by far the most species-rich cichlid fish tribe that originated along with the so-called primary radiation of the Lake Tanganyika cichlid species flock, i.e. at the same time during which the majority of the endemic Lake Tanganyika cichlid tribes emerged. Unlike the other tribes, the haplochromines are not restricted to Lake Tanganyika but distributed throughout Africa, except for the northwestern part of the continent. Haplochromine cichlids seeded the adaptive radiation of cichlid fishes in Lakes Malawi, Kivu, Victoria, Turkana, as well as in the now extinct paleo-Lake Makgadikgadi. Here we present a comprehensive phylogenetic and phylogeographic analysis of haplochromine cichlids that is based upon DNA sequences of two mitochondrial gene segments of riverine taxa covering all major African biogeographic regions where haplochromines are found. Our analysis revealed that six lineages of haplochromines originated within a short period of time, about 5.3-4.4 MYA. These haplochromine lineages show a highly complex phylogeographic pattern, probably severely influenced by climate- and/or geology-induced changes of the environment, with river capture events most likely playing an important role for species dispersal.