Escherichia coli 6-pyruvoyltetrahydropterin synthase ortholog encoded by ygcM has a new catalytic activity for conversion of sepiapterin to 7,8-dihydropterin

Genetic Basis and Evolutionary Forces of Sexually Dimorphic Color Variation in a Toad-Headed Agamid Lizard

In the animal kingdom, sexually dimorphic color variation is a widespread phenomenon that significantly influences survival and reproductive success. However, the genetic underpinnings of this variation remain inadequately understood. Our investigation into sexually dimorphic color variation in the desert-dwelling Guinan population of the toad-headed agamid lizard (Phrynocephalus putjatai) utilized a multidisciplinary approach, encompassing phenotypic, ultrastructural, biochemical, genomic analyses, and behavioral experiments. Our findings unveil the association between distinct skin colorations and varying levels of carotenoid and pteridine pigments. The red coloration in males is determined by a genomic region on chromosome 14, housing four pigmentation genes: BCO2 and three 6-pyruvoyltetrahydropterin synthases. A Guinan population-specific nonsynonymous single nucleotide polymorphism in BCO2 is predicted to alter the electrostatic potential within the binding domain of the BCO2-β-carotene complex, influencing their interaction. Additionally, the gene MAP7 on chromosome 2 emerges as a potential contributor to the blue coloration in subadults and adult females. Sex-specific expression patterns point to steroid hormone-associated genes (SULT2B1 and SRD5A2) as potential upstream regulators influencing sexually dimorphic coloration. Visual modeling and field experiments support the potential selective advantages of vibrant coloration in desert environments. This implies that natural selection, potentially coupled with assortative mating, might have played a role in fixing color alleles, contributing to prevalence in the local desert habitat. This study provides novel insights into the genetic basis of carotenoid and pteridine-based color variation, shedding light on the evolution of sexually dimorphic coloration in animals. Moreover, it advances our understanding of the driving forces behind such intricate coloration patterns.

Red carotenoids and associated gene expression explain colour variation in frillneck lizards

A long-standing hypothesis in evolutionary ecology is that red-orange ornamental colours reliably signal individual quality owing to limited dietary availability of carotenoids and metabolic costs associated with their production, such as the bioconversion of dietary yellow carotenoids to red ketocarotenoids. However, in ectothermic vertebrates, these colours can also be produced by self-synthesized pteridine pigments. As a consequence, the relative ratio of pigment types and their biochemical and genetic basis have implications for the costs and information content of colour signals; yet they remain poorly known in most taxonomic groups. We tested whether red- and yellow-frilled populations of the frillneck lizard, Chlamydosaurus kingii, differ in the ratio of different biochemical classes of carotenoid and pteridine pigments, and examined associated differences in gene expression. We found that, unlike other squamate reptiles, red hues derive from a higher proportion of ketocarotenoids relative to both dietary yellow carotenoids and to pteridines. Whereas red frill skin showed higher expression of several genes associated with carotenoid metabolism, yellow frill skin showed higher expression of genes associated with steroid hormones. Based on the different mechanisms underlying red and yellow signals, we hypothesize that frill colour conveys different information in the two populations. More generally, the data expand our knowledge of the genetic and biochemical basis of colour signals in vertebrates.

Comparative Transcriptome and DNA methylation analyses of the molecular mechanisms underlying skin color variations in Crucian carp (Carassius carassius L.)

Background

Crucian carp is a popular ornamental strain in Asia with variants in body color. To further explore the genetic mechanisms underlying gray and red body color formation in crucian carp, the skin transcriptomes and partial DNA methylation sites were obtained from red crucian carp (RCC) and white crucian carp (WCC). Here, we show significant differences in mRNA expression and DNA methylation sites between skin tissues of RCC and WCC.

Results

Totals of 3434 and 3683 unigenes had significantly lower and higher expression in WCC, respectively, compared with unigenes expressed in RCC. Some potential genes for body color development were further identified by quantitative polymerase chain reaction, such as mitfa, tyr, tyrp1, and dct, which were down-regulated, and foxd3, hpda, ptps, and gch1, which were up-regulated. A KEGG pathway analysis indicated that the differentially expressed genes were mainly related to mitogen activated protein kinase (MAPK), Wnt, cell cycle, and endocytosis signaling pathways, as well as variations in melanogenesis in crucian carp. In addition, some differentially expressed DNA methylation site genes were related to pigmentation, including mitfa, tyr, dct, foxd3, and hpda. The differentially expressed DNA methylation sites were mainly involved in signaling pathways, including MAPK, cAMP, endocytosis, melanogenesis, and Hippo.

Conclusions

Our study provides the results of comparative transcriptome and DNA methylation analyses between RCC and WCC skin tissues and reveals that the molecular mechanism of body color variation in crucian carp is strongly related to disruptions in gene expression and DNA methylation during pigmentation.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0564-9) contains supplementary material, which is available to authorized users.

Structural basis of a novel activity of bacterial 6-pyruvoyltetrahydropterin synthase homologues distinct from mammalian 6-pyruvoyltetrahydropterin synthase activity

Escherichia coli 6-carboxytetrahydropterin synthase (eCTPS), a homologue of 6-pyruvoyltetrahydropterin synthase (PTPS), possesses a much stronger catalytic activity to cleave the side chain of sepiapterin in vitro compared with genuine PTPS activity and catalyzes the conversion of dihydroneopterin triphosphate to 6-carboxy-5,6,7,8-tetrahydropterin in vivo. Crystal structures of wild-type apo eCTPS and of a Cys27Ala mutant eCTPS complexed with sepiapterin have been determined to 2.3 and 2.5 Å resolution, respectively. The structures are highly conserved at the active site and the Zn(2+) binding site. However, comparison of the eCTPS structures with those of mammalian PTPS homologues revealed that two specific residues, Trp51 and Phe55, that are not found in mammalian PTPS keep the substrate bound by stacking it with their side chains. Replacement of these two residues by site-directed mutagenesis to the residues Met and Leu, which are only found in mammalian PTPS, converted eCTPS to the mammalian PTPS activity. These studies confirm that these two aromatic residues in eCTPS play an essential role in stabilizing the substrate and in the specific enzyme activity that differs from the original PTPS activity. These aromatic residues Trp51 and Phe55 are a key signature of bacterial PTPS enzymes that distinguish them from mammalian PTPS homologues.

Functional promiscuity of the COG0720 family

The biosynthesis of GTP derived metabolites such as tetrahydrofolate (THF), biopterin (BH(4)), and the modified tRNA nucleosides queuosine (Q) and archaeosine (G(+)) relies on several enzymes of the Tunnel-fold superfamily. A subset of these proteins includes the 6-pyruvoyltetrahydropterin (PTPS-II), PTPS-III, and PTPS-I homologues, all members of the COG0720 family that have been previously shown to transform 7,8-dihydroneopterin triphosphate (H(2)NTP) into different products. PTPS-II catalyzes the formation of 6-pyruvoyltetrahydropterin in the BH(4) pathway, PTPS-III catalyzes the formation of 6-hydroxylmethyl-7,8-dihydropterin in the THF pathway, and PTPS-I catalyzes the formation of 6-carboxy-5,6,7,8-tetrahydropterin in the Q pathway. Genes of these three enzyme families are often misannotated as they are difficult to differentiate by sequence similarity alone. Using a combination of physical clustering, signature motif, phylogenetic codistribution analyses, in vivo complementation studies, and in vitro enzymatic assays, a complete reannotation of the COG0720 family was performed in prokaryotes. Notably, this work identified and experimentally validated dual function PTPS-I/III enzymes involved in both THF and Q biosynthesis. Both in vivo and in vitro analyses showed that the PTPS-I family could tolerate a translation of the active site cysteine and was inherently promiscuous, catalyzing different reactions on the same substrate or the same reaction on different substrates. Finally, the analysis and experimental validation of several archaeal COG0720 members confirmed the role of PTPS-I in archaeosine biosynthesis and resulted in the identification of PTPS-III enzymes with variant signature sequences in Sulfolobus species. This study reveals an expanded versatility of the COG0720 family members and illustrates that for certain protein families extensive comparative genomic analysis beyond homology is required to correctly predict function.

Biochemical characterization of the tetrahydrobiopterin synthesis pathway in the oleaginous fungus Mortierella alpina

We characterized the de novo biosynthetic pathway of tetrahydrobiopterin (BH₄) in the lipid-producing fungus Mortierella alpina. The BH₄ cofactor is essential for various cell processes, and is probably present in every cell or tissue of higher organisms. Genes encoding two copies of GTP cyclohydrolase I (GTPCH-1 and GTPCH-2) for the conversion of GTP to dihydroneopterin triphosphate (H₂-NTP), 6-pyruvoyltetrahydropterin synthase (PTPS) for the conversion of H₂-NTP to 6-pyruvoyltetrahydropterin (PPH₄), and sepiapterin reductase (SR) for the conversion of PPH₄ to BH₄, were expressed heterologously in Escherichia coli. The recombinant enzymes were produced as His-tagged fusion proteins and were purified to homogeneity to investigate their enzymic activities. Enzyme products were analysed by HPLC and electrospray ionization-MS. Kinetic parameters and other properties of GTPCH, PTPS and SR were investigated. Physiological roles of BH₄ in M. alpina are discussed, and comparative analyses between GTPCH, PTPS and SR proteins and other homologous proteins were performed. The presence of two functional GTPCH enzymes has, as far as we are aware, not been reported previously, reflecting the unique ability of this fungus to synthesize both BH₄ and folate, using the GTPCH product as a common substrate. To our knowledge, this study is the first to report the comprehensive characterization of a BH₄ biosynthesis pathway in a fungus.

The evolution of metabolic networks of E. coli

Background

Despite the availability of numerous complete genome sequences from E. coli strains, published genome-scale metabolic models exist only for two commensal E. coli strains. These models have proven useful for many applications, such as engineering strains for desired product formation, and we sought to explore how constructing and evaluating additional metabolic models for E. coli strains could enhance these efforts.

Results

We used the genomic information from 16 E. coli strains to generate an E. coli pangenome metabolic network by evaluating their collective 76,990 ORFs. Each of these ORFs was assigned to one of 17,647 ortholog groups including ORFs associated with reactions in the most recent metabolic model for E. coli K-12. For orthologous groups that contain an ORF already represented in the MG1655 model, the gene to protein to reaction associations represented in this model could then be easily propagated to other E. coli strain models. All remaining orthologous groups were evaluated to see if new metabolic reactions could be added to generate a pangenome-scale metabolic model (iEco1712_pan). The pangenome model included reactions from a metabolic model update for E. coli K-12 MG1655 (iEco1339_MG1655) and enabled development of five additional strain-specific genome-scale metabolic models. These additional models include a second K-12 strain (iEco1335_W3110) and four pathogenic strains (two enterohemorrhagic E. coli O157:H7 and two uropathogens). When compared to the E. coli K-12 models, the metabolic models for the enterohemorrhagic (iEco1344_EDL933 and iEco1345_Sakai) and uropathogenic strains (iEco1288_CFT073 and iEco1301_UTI89) contained numerous lineage-specific gene and reaction differences. All six E. coli models were evaluated by comparing model predictions to carbon source utilization measurements under aerobic and anaerobic conditions, and to batch growth profiles in minimal media with 0.2% (w/v) glucose. An ancestral genome-scale metabolic model based on conserved ortholog groups in all 16 E. coli genomes was also constructed, reflecting the conserved ancestral core of E. coli metabolism (iEco1053_core). Comparative analysis of all six strain-specific E. coli models revealed that some of the pathogenic E. coli strains possess reactions in their metabolic networks enabling higher biomass yields on glucose. Finally the lineage-specific metabolic traits were compared to the ancestral core model predictions to derive new insight into the evolution of metabolism within this species.

Conclusion

Our findings demonstrate that a pangenome-scale metabolic model can be used to rapidly construct additional E. coli strain-specific models, and that quantitative models of different strains of E. coli can accurately predict strain-specific phenotypes. Such pangenome and strain-specific models can be further used to engineer metabolic phenotypes of interest, such as designing new industrial E. coli strains.

Quantitative mass spectrometry catalogues Salmonella pathogenicity island-2 effectors and identifies their cognate host binding partners

Gram-negative bacterial pathogens have developed specialized secretion systems to transfer bacterial proteins directly into host cells. These bacterial effectors are central to virulence and reprogram host cell processes to favor bacterial survival, colonization, and proliferation. Knowing the complete set of effectors encoded by a particular pathogen is the key to understanding bacterial disease. In addition, the identification of the molecular assemblies that these effectors engage once inside the host cell is critical to determining the mechanism of action of each effector. In this work we used stable isotope labeling of amino acids in cell culture (SILAC), a powerful quantitative proteomics technique, to identify the proteins secreted by the Salmonella pathogenicity island-2 type three secretion system (SPI-2 T3SS) and to characterize the host interaction partners of SPI-2 effectors. We confirmed many of the known SPI-2 effectors and were able to identify several novel substrate candidates of this secretion system. We verified previously published host protein-effector binding pairs and obtained 11 novel interactions, three of which were investigated further and confirmed by reciprocal co-immunoprecipitation. The host cell interaction partners identified here suggest that Salmonella SPI-2 effectors target, in a concerted fashion, cellular processes such as cell attachment and cell cycle control that are underappreciated in the context of infection. The technology outlined in this study is specific and sensitive and serves as a robust tool for the identification of effectors and their host targets that is readily amenable to the study of other bacterial pathogens.

An Fe2+-dependent cyclic phosphodiesterase catalyzes the hydrolysis of 7,8-dihydro-D-neopterin 2',3'-cyclic phosphate in methanopterin biosynthesis

7,8-Dihydro-D-neopterin 2',3'-cyclic phosphate (H(2)N-cP) is the first intermediate in biosynthesis of the pterin portion of tetrahydromethanopterin (H(4)MPT), a C(1) carrier coenzyme first identified in the methanogenic archaea. This intermediate is produced from GTP by MptA (MJ0775 gene product), a new class of GTP cyclohydrolase I [Grochowski, L. L., Xu, H., Leung, K., and White, R. H. (2007) Biochemistry 46, 6658-6667]. Here we report the identification of a cyclic phosphodiesterase that hydrolyzes the cyclic phosphate of H(2)N-cP and converts it to a mixture of 7,8-dihydro-D-neopterin 2'-monophosphate and 7,8-dihydro-d-neopterin 3'-monophosphate. The enzyme from Methanocaldococcus jannachii is designated MptB (MJ0837 gene product) to indicate that it catalyzes the second step of the biosynthesis of methanopterin. MptB is a member of the HD domain superfamily of enzymes, which require divalent metals for activity. Direct metal analysis of the recombinant enzyme demonstrated that MptB contained 1.0 mol of zinc and 0.8 mol of iron per protomer. MptB requires Fe(2+) for activity, the same as observed for MptA. Thus the first two enzymes involved in H(4)MPT biosynthesis in the archaea are Fe(2+) dependent.

6-pyruvoyltetrahydropterin synthase paralogs replace the folate synthesis enzyme dihydroneopterin aldolase in diverse bacteria

Dihydroneopterin aldolase (FolB) catalyzes conversion of dihydroneopterin to 6-hydroxymethyldihydropterin (HMDHP) in the classical folate biosynthesis pathway. However, folB genes are missing from the genomes of certain bacteria from the phyla Chloroflexi, Acidobacteria, Firmicutes, Planctomycetes, and Spirochaetes. Almost all of these folB-deficient genomes contain an unusual paralog of the tetrahydrobiopterin synthesis enzyme 6-pyruvoyltetrahydropterin synthase (PTPS) in which a glutamate residue replaces or accompanies the catalytic cysteine. A similar PTPS paralog from the malaria parasite Plasmodium falciparum is known to form HMDHP from dihydroneopterin triphosphate in vitro and has been proposed to provide a bypass to the FolB step in vivo. Bacterial genes encoding PTPS-like proteins with active-site glutamate, cysteine, or both residues were accordingly tested together with the P. falciparum gene for complementation of the Escherichia coli folB mutation. The P. falciparum sequence and bacterial sequences with glutamate or glutamate plus cysteine were active; those with cysteine alone were not. These results demonstrate that PTPS paralogs with an active-site glutamate (designated PTPS-III proteins) can functionally replace FolB in vivo. Recombinant bacterial PTPS-III proteins, like the P. falciparum enzyme, mediated conversion of dihydroneopterin triphosphate to HMDHP, but other PTPS proteins did not. Neither PTPS-III nor other PTPS proteins exhibited significant dihydroneopterin aldolase activity. Phylogenetic analysis indicated that PTPS-III proteins may have arisen independently in various PTPS lineages. Consistent with this possibility, merely introducing a glutamate residue into the active site of a PTPS protein conferred incipient activity in the growth complementation assay, and replacing glutamate with alanine in a PTPS-III protein abolished complementation.

Escherichia coli QueD is a 6-carboxy-5,6,7,8-tetrahydropterin synthase

To elucidate the early steps required during biosynthesis of a broad class of 7-deazapurine-containing natural products, we have studied the reaction catalyzed by Escherichia coli QueD, a 6-pyruvoyl-5,6,7,8-tetrahydropterin synthase (PTPS) homologue possibly involved in queuosine biosynthesis. While mammalian PTPS homologues convert 7,8-dihydroneopterin triphosphate (H(2)NTP) to 6-pyruvoyltetrahydropterin (PPH(4)) in biopterin biosynthesis, E. coli QueD catalyzes the conversion of H(2)NTP to 6-carboxy-5,6,7,8-tetrahydropterin (CPH(4)). E. coli QueD can also convert PPH(4) and sepiapterin to CPH(4), allowing a mechanism to be proposed.

An atypical orthologue of 6-pyruvoyltetrahydropterin synthase can provide the missing link in the folate biosynthesis pathway of malaria parasites

Folate metabolism in malaria parasites is a long-standing, clinical target for chemotherapy and prophylaxis. However, despite determination of the complete genome sequence of the lethal species Plasmodium falciparum, the pathway of de novo folate biosynthesis remains incomplete, as no candidate gene for dihydroneopterin aldolase (DHNA) could be identified. This enzyme catalyses the third step in the well-characterized pathway of plants, bacteria, and those eukaryotic microorganisms capable of synthesizing their own folate. Utilizing bioinformatics searches based on both primary and higher protein structures, together with biochemical assays, we demonstrate that P. falciparum cell extracts lack detectable DHNA activity, but that the parasite possesses an unusual orthologue of 6-pyruvoyltetrahydropterin synthase (PTPS), which simultaneously gives rise to two products in comparable amounts, the predominant of which is 6-hydroxymethyl-7,8-dihydropterin, the substrate for the fourth step in folate biosynthesis (catalysed by 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase; PPPK). This can provide a bypass for the missing DHNA activity and thus a means of completing the biosynthetic pathway from GTP to dihydrofolate. Supported by site-directed mutagenesis experiments, we ascribe the novel catalytic activity of the malarial PTPS to a Cys to Glu change at its active site relative to all previously characterized PTPS molecules, including that of the human host.

Evolution of new function in the GTP cyclohydrolase II proteins of Streptomyces coelicolor

The genome sequence of Streptomyces coelicolor contains three open reading frames (sco1441, sco2687, and sco6655) that encode proteins with significant (>40%) amino acid identity to GTP cyclohydrolase II (GCH II), which catalyzes the committed step in the biosynthesis of riboflavin. The physiological significance of the redundancy of these proteins in S. coelicolor is not known. However, the gene contexts of the three proteins are different, suggesting that they may serve alternate biological niches. Each of the three proteins was overexpressed in Escherichia coli and characterized to determine if their functions are biologically overlapping. As purified, each protein contains 1 molar equiv of zinc/mol of protein and utilizes guanosine 5'-triphosphate (GTP) as substrate. Two of these proteins (SCO 1441 and SCO 2687) produce the canonical product of GCH II, 2,5-diamino-6-ribosylamino-4(3H)-pyrimidinone 5'-phosphate (APy). Remarkably, however, one of the three proteins (SCO 6655) converts GTP to 2-amino-5-formylamino-6-ribosylamino-4(3H)-pyrimidinone 5'-phosphate (FAPy), as shown by UV-visible spectrophotometry, mass spectrometry, and NMR. This activity has been reported for a GTP cyclohydrolase III protein from Methanocaldococcus jannaschii [Graham, D. E., Xu, H., and White, R. H. (2002) Biochemistry 41, 15074-15084], which has no amino acid sequence homology to SCO 6655. Comparison of the sequences of these proteins and mapping onto the structure of the E. coli GCH II protein [Ren, J., Kotaka, M., Lockyer, M., Lamb, H. K., Hawkins, A. R., and Stammers, D. K. (2005) J. Biol. Chem. 280, 36912-36919] allowed identification of a switch residue, Met120, which appears to be responsible for the altered fate of GTP observed with SCO 6655; a Tyr is found in the analogous position of all proteins that have been shown to catalyze the conversion of GTP to APy. The Met120Tyr variant of SCO 6655 acquires the ability to catalyze the conversion of GTP to APy, suggesting a role for Tyr120 in the late phase of the reaction. Our data are consistent with duplication of GCH II in S. coelicolor promoting evolution of a new function. The physiological role(s) of the gene clusters that house GCH II homologues will be discussed.

6-Pyruvoyltetrahydropterin synthase orthologs of either a single or dual domain structure are responsible for tetrahydrobiopterin synthesis in bacteria

6-Pyruvoyltetrahydropterin synthase (PTPS) catalyzes the second step of tetrahydrobiopterin (BH4) synthesis. We previously identified PTPS orthologs (bPTPS-Is) in bacteria which do not produce BH4. In this study we disrupted the gene encoding bPTPS-I in Synechococcus sp. PCC 7942, which produces BH4-glucoside. The mutant was normal in BH4-glucoside production, demonstrating that bPTPS-I does not participate in BH4 synthesis in vivo and bringing us a new PTPS ortholog (bPTPS-II) of a bimodular polypeptide. The recombinant Synechococcus bPTPS-II was assayed in vitro to show PTPS activity higher than human enzyme. Further computational analysis revealed the presence of mono and bimodular bPTPS-II orthologs mostly in green sulfur bacteria and cyanobacteria, respectively, which are well known for BH4-glycoside production. In summary we found new bacterial PTPS orthologs, having either a single or dual domain structure and being responsible for BH4 synthesis in vivo, thereby disclosing all the bacterial PTPS homologs.

Sepiapterin reductases from Chlorobium tepidum and Chlorobium limicola catalyze the synthesis of L-threo-tetrahydrobiopterin from 6-pyruvoyltetrahydropterin

The ORF sequences of the gene encoding sepiapterin reductase were cloned from the genomic DNAs of Chlorobium tepidum and Chlorobium limicola, which are known to produce L-threo- and L-erythro-tetrahydrobiopterin (BH4)-N-acetylglucosamine, respectively. The deduced amino acid sequence of C. limicola consists of 241 residues, while C. tepidum SR has three residues more at the C-terminal. The overall protein sequence identity was 87.7%. Both recombinant proteins generated from Escherichia coli were identified to catalyze reduction of diketo compound 6-pyruvoyltetrahydropterin to L-threo-BH4. This result suggests that C. limicola needs an additional enzyme for L-erythro-BH4 synthesis to yield its glycoside. The catalytic activity of Chlorobium SRs also supports the previously proposed mechanism of two consecutive reductions of C1' carbonyl group of 6-pyruvoyltetrahydropterin via isomerization reaction.

Identification of four genes necessary for biosynthesis of the modified nucleoside queuosine

Queuosine (Q) is a hypermodified 7-deazaguanosine nucleoside located in the anticodon wobble position of four amino acid-specific tRNAs. In bacteria, Q is produced de novo from GTP via the 7-deazaguanosine precursor preQ1 (7-aminoethyl 7-deazaguanine) by an uncharacterized pathway. PreQ1 is subsequently transferred to its specific tRNA by a tRNA-guanine transglycosylase (TGT) and then further modified in situ to produce Q. Here we use comparative genomics to implicate four gene families (best exemplified by the B. subtilis operon ykvJKLM) as candidates in the preQ1 biosynthetic pathway. Deletions were constructed in genes for each of the four orthologs in Acinetobacter. High pressure liquid chromatography analysis showed the Q nucleoside was absent from the tRNAs of each of four deletion strains. Electrospray ionization mass spectrometry confirmed the absence of Q in each mutant strain. Finally, introduction of the Bacillus subtilis ykvJKLM operon in trans complemented the Q deficiency of the two deletion mutants that were tested. Thus, the products of these four genes (named queC, -D, -E, and -F) are essential for the Q biosynthetic pathway.

The pteridine pathway in zebrafish: regulation and specification during the determination of neural crest cell-fate

This review describes pteridine biosynthesis and its relation to the differentiation of neural crest derivatives in zebrafish. During the embryonic development of these fish, neural crest precursor cells segregate into neural elements, ectomesenchymal cells and pigment cells; the latter then diversifying into melanophores, iridophores and xanthophores. The differentiation of neural cells, melanophores, and xanthophores is coupled closely with the onset of pteridine synthesis which starts from GTP and is regulated through the control of GTP cyclohydrolase I activity. De novo pteridine synthesis in embryos of this species increases during the first 72-h postfertilization, producing H4biopterin, which serves as a cofactor for neurotransmitter synthesis in neural cells and for tyrosine production in melanophores. Thereafter, sepiapterin (6-lactoyl-7,8-dihydropterin) accumulates as yellow pigment in xanthophores, together with 7-oxobiopterin, isoxanthopterin and 2,4,7-trioxopteridine. Sepiapterin is the key intermediate in the formation of 7-oxopteridines, which depends on the availability of enzymes belonging to the xanthine oxidoreductase family. Expression of the GTP cyclohydrolase I gene (gch) is found in neural cells, in melanoblasts and in early xanthophores (xanthoblasts) of early zebrafish embryos but steeply declines in xanthophores by 42-h postfertilization. The mechanism(s) whereby sepiapterin branches off from the GTP-H4biopterin pathway is currently unknown and will require further study. The surge of interest in zebrafish as a model for vertebrate development and its amenability to genetic manipulation provide powerful tools for analysing the functional commitment of neural crest-derived cells and the regulation of pteridine synthesis in mammals.