Illuminating the detection chain of bacterial bioreporters.

Engineering bacteria for measuring chemicals of environmental or toxicological concern (bioreporter bacteria) has grown slowly into a mature research area. Despite many potential advantages, current bioreporters do not perform well enough to comply with environmental detection standards. Basically, the reasons for this are the lack of engineering principles in the detection chain in the bioreporters. Here, we dissect critical steps in the detection chain and illustrate how bioreporter design could be improved by mutagenizing specificity and selectivity of the sensing and regulatory proteins, by newer expression strategies and application of different signalling networks. Furthermore, we describe how redesigning bioreporter assays with respect to pollutant transport into the cells and application of other detection devices can decrease detection limits and increase the speed of detection.

Complex organization of CTF/NF-I, C/EBP, and HNF3 binding sites within the promoter of the liver-specific vitellogenin gene.

Vitellogenin genes are expressed specifically in the liver of female oviparous vertebrates under the strict control of estrogen. To explain this tissue-specific expression, we performed a detailed analysis of the Xenopus laevis vitellogenin gene B1 promoter by DNase I footprinting and gel mobility-shift assays. We characterized five binding sites for the ubiquitous factor CTF/NF-I. Two of these sites are close to the TATA-box, whereas the others are located on both sides of the estrogen responsive unit formed by two imperfect estrogen response elements. Moreover two liver-enriched factors, C/EBP and HNF3, were found to interact with multiple closely spaced proximal promoter elements in the first 100 base pairs upstream of the TATA-box. To confirm the physiological significance of this in vitro analysis, in vivo DNase I footprinting experiments were carried out using the ligation-mediated polymerase chain reaction technique. The various cis-elements characterized in vitro as binding sites for known transcription factors and more particularly for liver-enriched transcription factors are efficiently recognized in vivo as well, suggesting that they play an important role in the control of the liver-specific vitellogenin gene B1 expression.

The transcriptional repressor REST determines the cell-specific expression of the human MAPK8IP1 gene encoding IB1 (JIP-1).

Islet-brain 1 (IB1) is the human and rat homologue of JIP-1, a scaffold protein interacting with the c-Jun amino-terminal kinase (JNK). IB1 expression is mostly restricted to the endocrine pancreas and to the central nervous system. Herein, we explored the transcriptional mechanism responsible for this preferential islet and neuronal expression of IB1. A 731-bp fragment of the 5' regulatory region of the human MAPK8IP1 gene was isolated from a human BAC library and cloned upstream of a luciferase reporter gene. This construct drove high transcriptional activity in both insulin-secreting and neuron-like cells but not in unrelated cell lines. Sequence analysis of this promoter region revealed the presence of a neuron-restrictive silencer element (NRSE) known to bind repressor zinc finger protein REST. This factor is not expressed in insulin-secreting and neuron-like cells. By mobility shift assay, we confirmed that REST binds to the NRSE present in the IB1 promoter. Once transiently transfected in beta-cell lines, the expression vector encoding REST repressed IB1 transcriptional activity. The introduction of a mutated NRSE in the 5' regulating region of the IB1 gene abolished the repression activity driven by REST in insulin-secreting beta cells and relieved the low transcriptional activity of IB1 observed in unrelated cells. Moreover, transfection in non-beta and nonneuronal cell lines of an expression vector encoding REST lacking its transcriptional repression domain relieved IB1 promoter activity. Last, the REST-mediated repression of IB1 could be abolished by trichostatin A, indicating that deacetylase activity is required to allow REST repression. Taken together, these data establish a critical role for REST in the control of the tissue-specific expression of the human IB1 gene.

Genetics of cadmium hyperaccumulation variance in Thlaspi caerulescens natural populations from Switzerland

SUMMARY Heavy metal presence in the environment is a serious concern since some of them can be toxic to plants, animals and humans once accumulated along the food chain. Cadmium (Cd) is one of the most toxic heavy metal. It is naturally present in soils at various levels and its concentration can be increased by human activities. Several plants however have naturally developed strategies allowing them to grow on heavy metal enriched soils. One of them consists in the accumulation and sequestration of heavy metals in the above-ground biomass. Some plants present in addition an extreme strategy by which they accumulate a limited number of heavy metals in their shoots in amounts 100 times superior to those expected for a non-accumulating plant in the same conditions. Understanding the genetic basis of the hyperaccumulation trait - particularly for Cd - remains an important challenge which may lead to biotechnological applications in the soil phytoremediation. In this thesis, Thlaspi caerulescens J. & C. Presl (Brassicaceae) was used as a model plant to study the Cd hyperaccumulation trait, owing to its physiological and genetic characteristics. Twenty-four wild populations were sampled in different regions of Switzerland. They were characterized for environmental and soil parameters as well as intrinsic characteristics of plants (i.e. metal concentrations in shoots). They were as well genetically characterized by AFLPs, plastid DNA polymorphism and genes markers (CAPS and microsatellites) mainly developed in this thesis. Some of the investigated genes were putatively linked to the Cd hyperaccumulation trait. Since the study of the Cd hyperaccumulation in the field is important as it allows the identification of patterns of selection, the present work offered a methodology to define the Cd hyperaccumulation capacity of populations from different habitats permitting thus their comparison in the field. We showed that Cd, Zn, Fe and Cu accumulations were linked and that populations with higher Cd hyperaccumulation capacity had higher shoot and reproductive fitness. Using our genetic data, statistical methods (Beaumont & Nichols's procedure, partial Mantel tests) were applied to identify genomic signatures of natural selection related to the Cd hyperaccumulation capacity. A significant genetic difference between populations related to their Cd hyperaccumulation capacity was revealed based on somè specific markers (AFLP and candidate genes). Polymorphism at the gene encoding IRTl (Iron-transporter also participating to the transport of Zn) was suggested as explaining part of the variation in Cd hyperaccumulation capacity of populations supporting previous physiological investigations. RÉSUMÉ La présence de métaux lourds dans l'environnement est un phénomène préoccupant. En effet, certains métaux lourds - comme le cadmium (Cd) -sont toxiques pour les plantes, les animaux et enfin, accumulés le long de la chaîne alimentaire, pour les hommes. Le Cd est naturellement présent dans le sol et sa concentration peut être accrue par différentes activités humaines. Certaines plantes ont cependant développé des stratégies leur permettant de pousser sur des sols contaminés en métaux lourds. Parmi elles, certaines accumulent et séquestrent les métaux lourds dans leurs parties aériennes. D`autres présentent une stratégie encore plus extrême. Elles accumulent un nombre limité de métaux lourds en quantités 100 fois supérieures à celles attendues pour des espèces non-accumulatrices sous de mêmes conditions. La compréhension des bases génétiques de l'hyperaccumulation -particulièrement celle du Cd - représente un défi important avec des applications concrètes en biotechnologies, tout particulièrement dans le but appliqué de la phytoremediation des sols contaminés. Dans cette thèse, Thlaspi caerulescens J. & C. Presl (Brassicaceae) a été utilisé comme modèle pour l'étude de l'hyperaccumulation du Cd de par ses caractéristiques physiologiques et génétiques. Vingt-quatre populations naturelles ont été échantillonnées en Suisse et pour chacune d'elles les paramètres environnementaux, pédologique et les caractéristiques intrinsèques aux plantes (concentrations en métaux lourds) ont été déterminés. Les populations ont été caractérisées génétiquement par des AFLP, des marqueurs chloroplastiques et des marqueurs de gènes spécifiques, particulièrement ceux potentiellement liés à l'hyperaccumulation du Cd (CAPS et microsatellites). La plupart ont été développés au cours de cette thèse. L'étude de l'hyperaccumulation du Cd en conditions naturelles est importante car elle permet d'identifier la marque, éventuelle de sélection naturelle. Ce travail offre ainsi une méthodologie pour définir et comparer la capacité des populations à hyperaccumuler le Cd dans différents habitats. Nous avons montré que les accumulations du Cd, Zn, Fe et Cu sont liées et que les populations ayant une grande capacité d'hyperaccumuler le Cd ont également une meilleure fitness végétative et reproductive. Des méthodes statistiques (l'approche de Beaumont & Nichols, tests de Martel partiels) ont été utilisées sur les données génétiques pour identifier la signature génomique de la sélection naturelle liée à la capacité d'hyperaccumuler le Cd. Une différenciation génétique des populations liée à leur capacité d'hyperaccumuler le Cd a été mise en évidence sur certains marqueurs spécifiques. En accord avec les études physiologiques connues, le polymorphisme au gène codant IRT1 (un transporteur de Fe impliqué dans le transport du Zn) pourrait expliquer une partie de la variance de la capacité des populations à hyperaccumuler le Cd.

Adaptive divergence of ancient gene duplicates in the avian MHC class II beta.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

Infrequent expression of the MAGE gene family in uveal melanomas.

It has previously been reported that MAGE-1, -2, -3 and -4 genes are expressed in human cancers including cutaneous melanoma. MAGE-1 and MAGE-3 represent targets for specific immunotherapy because they encode peptide antigens which are recognised by cytotoxic T lymphocytes (CTL) when presented by HLA class I molecules, and pilot clinical trials with these peptides are currently in progress. It is likely that other members of the MAGE gene family may also encode antigens recognised by CTL. Uveal melanomas, like cutaneous melanomas, arise from melanocytes that are derived from the neural crest. To determine if uveal melanoma patients would be suitable for MAGE-peptide immunotherapy, the expression of MAGE-1, -2, -3 and -4 genes was assessed by reverse transcription followed by polymerase chain reaction (RT-PCR) amplification and ethidium bromide staining. Expression of MAGE genes was not detected in any of 27 primary tumours. Either MAGE-1 or MAGE-4 was expressed in only 2 of 26 metastatic samples, but expression of MAGE-2 or -3 was not detected. Our data suggest that, unlike cutaneous melanomas, uveal melanomas may not be suitable candidates for MAGE-peptide immunotherapy.

Repercussion of a deficiency in mitochondrial ß-oxidation on the carbon flux of short-chain fatty acids to the peroxisomal ß-oxidation cycle in Aspergillus nidulans.

The fungus Aspergillus nidulans contains both a mitochondrial and peroxisomal ß-oxidation pathway. This work was aimed at studying the influence of mutations in the foxA gene, encoding a peroxisomal multifunctional protein, or in the scdA/echA genes, encoding a mitochondrial short-chain dehydrogenase and an enoyl-CoA hydratase, respectively, on the carbon flux to the peroxisomal ß-oxidation pathway. A. nidulans transformed with a peroxisomal polyhydroxyalkanoate (PHA) synthase produced PHA from the polymerization of 3-hydroxyacyl-CoA intermediates derived from the peroxisomal ß-oxidation of external fatty acids. PHA produced from erucic acid or heptadecanoic acid contained a broad spectrum of monomers, ranging from 5 to 14 carbons, revealing that the peroxisomal ß-oxidation cycle can handle both long and short-chain intermediates. While the ∆foxA mutant grown on erucic acid or oleic acid synthesized 10-fold less PHA compared to wild type, the same mutant grown on octanoic acid or heptanoic acid produced 3- to 6-fold more PHA. Thus, while FoxA has an important contribution to the degradation of long-chain fatty acids, the flux of short-chain fatty acids to peroxisomal ß-oxidation is actually enhanced in its absence. While no change in PHA was observed in the ∆scdA∆echA mutant grown on erucic acid or oleic acid compared to wild type, there was a 2- to 4-fold increased synthesis of PHA in ∆scdA∆echA cells grown in octanoic acid or heptanoic acid. These results reveal that a compensatory mechanism exists in A. nidulans that increases the flux of short-chain fatty acids towards the peroxisomal ß-oxidation cycle when the mitochondrial ß-oxidation pathway is defective.

Use of the frc gene as a molecular marker to characterize oxalate-oxidizing bacterial abundance and diversity structure in soil.

Oxalate catabolism, which can have both medical and environmental implications, is performed by phylogenetically diverse bacteria. The formyl-CoA-transferase gene was chosen as a molecular marker of the oxalotrophic function. Degenerated primers were deduced from an alignment of frc gene sequences available in databases. The specificity of primers was tested on a variety of frc-containing and frc-lacking bacteria. The frc-primers were then used to develop PCR-DGGE and real-time SybrGreen PCR assays in soils containing various amounts of oxalate. Some PCR products from pure cultures and from soil samples were cloned and sequenced. Data were used to generate a phylogenetic tree showing that environmental PCR products belonged to the target physiological group. The extent of diversity visualised on DGGE pattern was higher for soil samples containing carbonate resulting from oxalate catabolism. Moreover, the amount of frc gene copies in the investigated soils was detected in the range of 1.64x10(7) to 1.75x10(8)/g of dry soil under oxalogenic tree (representing 0.5 to 1.2% of total 16S rRNA gene copies), whereas the number of frc gene copies in the reference soil was 6.4x10(6) (or 0.2% of 16S rRNA gene copies). This indicates that oxalotrophic bacteria are numerous and widespread in soils and that a relationship exists between the presence of the oxalogenic trees Milicia excelsa and Afzelia africana and the relative abundance of oxalotrophic guilds in the total bacterial communities. This is obviously related to the accomplishment of the oxalate-carbonate pathway, which explains the alkalinization and calcium carbonate accumulation occurring below these trees in an otherwise acidic soil. The molecular tools developed in this study will allow in-depth understanding of the functional implication of these bacteria on carbonate accumulation as a way of atmospheric CO(2) sequestration.

Increased endometrial placenta growth factor (PLGF) gene expression in women with successful implantation.

Objective: To analyze the vascularization of the endometrium via hysteroscopy and to assess its correlation with angiogenic factor gene expression and embryo implantation rate.Design: Cross-sectional study.Setting: Public university hospital.Patient(s): Patients undergoing hysteroscopy for supposed infertility.Intervention(s): Endometrial quality assessment according to Sakumoto-Masamoto, performed in the early secretory phase of the cycle. Collection of an endometrial tissue biopsy.Main Outcome Measure(s): RNA extraction, reverse transcription, and determination of gene expression of angiogenesis- and implantation-relevant factors using quantitative polymerase chain reaction. Retrieval of pregnancy information from the medical records.Result(s): Good quantity/quality RNA with infertility history was obtained from 63 participating women. Those with a "good" endometrium and subsequent pregnancy showed increased gene expression for placenta growth factor when compared with patients with a "bad" endometrium and who did not succeed with pregnancy to date. Nonpregnant women with a "good" endometrium presented an intermediate result. No significant differences were observed for several other genes tested, but trends in the same direction were observed.Conclusion(s): This study demonstrates for the first time that endometrial PLGF expression corresponds to the hysteroscopic appearance of the endometrium, and therefore has potential as a clinically relevant prognosticator for infertility treatment success. (Fertil Steril (R) 2011;96:663-8. (C)2011 by American Society for Reproductive Medicine.)

Differential Gene Expression between Fall- and Spring-Run Chinook Salmon Assessed by Long Serial Analysis of Gene Expression

Of all Pacific salmonids, Chinook salmon Oncorhynchus tshawytscha display the greatest variability in return times to freshwater. The molecular mechanisms of these differential return times have not been well described. Current methods, such as long serial analysis of gene expression (LongSAGE) and microarrays, allow gene expression to be analyzed for thousands of genes simultaneously. To investigate whether differential gene expression is observed between fall- and spring-run Chinook salmon from California's Central Valley, LongSAGE libraries were constructed. Three libraries containing between 25,512 and 29,372 sequenced tags (21 base pairs/tag) were generated using messenger RNA from the brains of adult Chinook salmon returning in fall and spring and from one ocean-caught Chinook salmon. Tags were annotated to genes using complementary DNA libraries from Atlantic salmon Salmo salar and rainbow trout O. mykiss. Differentially expressed genes, as estimated by differences in the number of sequence tags, were found in all pairwise comparisons of libraries (freshwater versus saltwater = 40 genes; fall versus spring = 11 genes: and spawning versus nonspawning = 51 genes). The gene for ependymin, an extracellular glycoprotein involved in behavioral plasticity in fish, exhibited the most differential expression among the three groupings. Reverse transcription polymerase chain reaction analysis verified the differential expression of ependymin between the fall- and spring-run samples. These LongSAGE libraries, the first reported for Chinook salmon, provide a window of the transcriptional changes during Chinook salmon return migration to freshwater and spawning and increase the amount of expressed sequence data.

A 338-bp proximal fragment of the glucose transporter type 2 (GLUT2) promoter drives reporter gene expression in the pancreatic islets of transgenic mice.

The high Km glucose transporter GLUT2 is a membrane protein expressed in tissues involved in maintaining glucose homeostasis, and in cells where glucose-sensing is necessary. In many experimental models of diabetes, GLUT2 gene expression is decreased in pancreatic beta-cells, which could lead to a loss of glucose-induced insulin secretion. In order to identify factors involved in pancreatic beta-cell specific expression of GLUT2, we have recently cloned the murine GLUT2 promoter and identified cis-elements within the 338-bp of the proximal promoter capable of binding islet-specific trans-acting factors. Furthermore, in transient transfection studies, this 338-bp fragment could efficiently drive the expression of the chloramphenicol acetyl transferase (CAT) gene in cell lines derived from the endocrine pancreas, but displayed no promoter activity in non-pancreatic cells. In this report, we tested the cell-specific expression of a CAT reporter gene driven by a short (338 bp) and a larger (1311 bp) fragment of the GLUT2 promoter in transgenic mice. We generated ten transgenic lines that integrated one of the constructs. CAT mRNA expression in transgenic tissues was assessed using the RNAse protection assay and the quantitative reverse transcribed polymerase chain reaction (RT-PCR). Overall CAT mRNA expression for both constructs was low compared to endogenous GLUT2 mRNA levels but the reporter transcript could be detected in all animals in the pancreatic islets and the liver, and in a few transgenic lines in the kidney and the small intestine. The CAT protein was also present in Langerhans islets and in the liver for both constructs by immunocytochemistry. These findings suggest that the proximal 338 bp of the murine GLUT2 promoter contain cis-elements required for the islet-specific expression of GLUT2.

Reverse transcription quantitative real-time polymerase chain reaction reference genes in the spared nerve injury model of neuropathic pain: validation and literature search.

BACKGROUND: The reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is a widely used, highly sensitive laboratory technique to rapidly and easily detect, identify and quantify gene expression. Reliable RT-qPCR data necessitates accurate normalization with validated control genes (reference genes) whose expression is constant in all studied conditions. This stability has to be demonstrated.We performed a literature search for studies using quantitative or semi-quantitative PCR in the rat spared nerve injury (SNI) model of neuropathic pain to verify whether any reference genes had previously been validated. We then analyzed the stability over time of 7 commonly used reference genes in the nervous system - specifically in the spinal cord dorsal horn and the dorsal root ganglion (DRG). These were: Actin beta (Actb), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal proteins 18S (18S), L13a (RPL13a) and L29 (RPL29), hypoxanthine phosphoribosyltransferase 1 (HPRT1) and hydroxymethylbilane synthase (HMBS). We compared the candidate genes and established a stability ranking using the geNorm algorithm. Finally, we assessed the number of reference genes necessary for accurate normalization in this neuropathic pain model. RESULTS: We found GAPDH, HMBS, Actb, HPRT1 and 18S cited as reference genes in literature on studies using the SNI model. Only HPRT1 and 18S had been once previously demonstrated as stable in RT-qPCR arrays. All the genes tested in this study, using the geNorm algorithm, presented gene stability values (M-value) acceptable enough for them to qualify as potential reference genes in both DRG and spinal cord. Using the coefficient of variation, 18S failed the 50% cut-off with a value of 61% in the DRG. The two most stable genes in the dorsal horn were RPL29 and RPL13a; in the DRG they were HPRT1 and Actb. Using a 0.15 cut-off for pairwise variations we found that any pair of stable reference gene was sufficient for the normalization process. CONCLUSIONS: In the rat SNI model, we validated and ranked Actb, RPL29, RPL13a, HMBS, GAPDH, HPRT1 and 18S as good reference genes in the spinal cord. In the DRG, 18S did not fulfill stability criteria. The combination of any two stable reference genes was sufficient to provide an accurate normalization.

Identification of a new cancer/testis gene family, CT47, among expressed multicopy genes on the human X chromosome.

Cancer/testis (CT) genes are normally expressed in germ cells only, yet are reactivated and expressed in some tumors. Of the approximately 40 CT genes or gene families identified to date, 20 are on the X chromosome and are present as multigene families, many with highly conserved members. This indicates that novel CT gene families may be identified by detecting duplicated expressed genes on chromosome X. By searching for transcript clusters that map to multiple locations on the chromosome, followed by in silico analysis of their gene expression profiles, we identified five novel gene families with testis-specific expression and >98% sequence identity among family members. The expression of these genes in normal tissues and various tumor cell lines and specimens was evaluated by qualitative and quantitative RT-PCR, and a novel CT gene family with at least 13 copies was identified on Xq24, designated as CT47. mRNA expression of CT47 was found mainly in the testes, with weak expression in the placenta. Brain tissue was the only positive somatic tissue tested, with an estimated CT47 transcript level 0.09% of that found in testis. Among the tumor specimens tested, CT47 expression was found in approximately 15% of lung cancer and esophageal cancer specimens, but not in colorectal cancer or breast cancer. The putative CT47 protein consists of 288 amino acid residues, with a C-terminus rich in alanine and glutamic acid. The only species other than human in which a gene homologous to CT47 has been detected is the chimpanzee, with the predicted protein showing approximately 80% identity in its carboxy terminal region.

Gene expression changes in chronic inflammatory demyelinating polyneuropathy skin biopsies.

Chronic-inflammatory demyelinating polyneuropathy (CIDP) is an immune-mediated disease with no known biomarkers for diagnosing the disease or assessing its prognosis. We performed transcriptional profiling microarray analysis on skin punch biopsies from 20 CIDP patients and 17 healthy controls to identify disease-associated gene expression changes. We demonstrate changes in expression of genes involved in immune and chemokine regulation, growth and repair. We also found a combination of two upregulated genes that can be proposed as a novel biomarker of the disorder.

Adaptive sequence evolution in a color gene involved in the formation of the characteristic egg-dummies of male haplochromine cichlid fishes.

BACKGROUND: The exceptionally diverse species flocks of cichlid fishes in East Africa are prime examples of parallel adaptive radiations. About 80% of East Africa's more than 1 800 endemic cichlid species, and all species of the flocks of Lakes Victoria and Malawi, belong to a particularly rapidly evolving lineage, the haplochromines. One characteristic feature of the haplochromines is their possession of egg-dummies on the males' anal fins. These egg-spots mimic real eggs and play an important role in the mating system of these maternal mouthbrooding fish. RESULTS: Here, we show that the egg-spots of haplochromines are made up of yellow pigment cells, xanthophores, and that a gene coding for a type III receptor tyrosine kinase, colony-stimulating factor 1 receptor a (csf1ra), is expressed in egg-spot tissue. Molecular evolutionary analyses reveal that the extracellular ligand-binding and receptor-interacting domain of csf1ra underwent adaptive sequence evolution in the ancestral lineage of the haplochromines, coinciding with the emergence of egg-dummies. We also find that csf1ra is expressed in the egg-dummies of a distantly related cichlid species, the ectodine cichlid Ophthalmotilapia ventralis, in which markings with similar functions evolved on the pelvic fin in convergence to those of the haplochromines. CONCLUSION: We conclude that modifications of existing signal transduction mechanisms might have evolved in the haplochromine lineage in association with the origination of anal fin egg-dummies. That positive selection has acted during the evolution of a color gene that seems to be involved in the morphogenesis of a sexually selected trait, the egg-dummies, highlights the importance of further investigations of the comparative genomic basis of the phenotypic diversification of cichlid fishes.