<title>Introduction</title> <p>Cholesterol and other sterols such as ergosterol and phytosterols are universal components of eukaryotic plasma membranes and are absent from the membranes of prokaryotes. These sterols (cholesterol, ergosterol and phytosterol) modulate membrane fluidity ##REF##15726825##[1]##, ##REF##11412894##[2]##. In addition to this structural role, cholesterol is essential for signaling processes. In fact, it is a precursor of steroid hormones, oxysterols, ecdysones (in insects) and vitamin D. Moreover, it may influence intercellular signaling through its covalent attachment to proteins such as the protein Hedgehog (Hh) in <italic>Drosophila</italic>
##REF##12894170##[3]##–##UREF##0##[5]##. The nematode <italic>Caenorhabditis elegans</italic> possesses several genes encoding proteins with regions similar to Hh and potentially undergoing cholesterylation ##REF##10523520##[6]##, ##REF##16873581##[7]##.</p> <p>Yeast, plants and mammals synthesize sterols through a series of complex reactions that occur in the endoplasmic reticulum (ER) and, therefore, most of the involved enzymes have transmembrane domains (##FIG##0##Figure 1##). We outline below the series of reactions to produce ergosterol in the yeast <italic>Saccharomyces cerevisiae</italic>, which involve <italic>ERG</italic> (i.e. from ERGosterol) genes and the respective ERG proteins.</p> <p>The first three steps of sterol biosynthesis are catalized by ERG9p (squalene synthase), ERG1p (squalene epoxidase) and ERG7p (lanosterol synthase), respectively. These proteins are essential for aerobic viability and their absence results in an inability to synthesize ergosterol. The third enzyme, ERG7p, converts squalene epoxide into lanosterol, the first cyclic component of the cholesterol biosynthesis cascade (i.e. the first sterol). The remaining enzymes of the pathway metabolize lanosterol into ergosterol. The sequential action of ERG11p (lanosterol demethylase) and ERG24p (C-14 reductase) leads to 4,4-dimethylzymosterol. They are also essential in aerobic conditions ##REF##7791529##[8]##–##REF##15995173##[10]##. Removal of the two C-4 methyl groups of 4,4-dimethylzymosterol is a complex reaction involving the products of genes <italic>ERG25</italic>, <italic>ERG26</italic>, and <italic>ERG27</italic>, in cooperation with ERG28p. Finally, ERG6p (C-24 methylase) converts zymosterol to fecosterol, which is further transformed into ergosterol by ERG2p, ERG3p, ERG5p and ERG4p ##REF##7791529##[8]##–##REF##15995173##[10]## (##FIG##0##Figure 1##).</p> <p>Recent results show that, in yeast, ergosterol biosynthetic enzymes display specific protein-protein interactions and form a functional complex called ergosome. Proteins ERG11p, ERG25p, ERG27p and ERG28p, appear to form a core center of the complex and would interact with other enzymes of the pathway ##UREF##1##[9]##, ##REF##15995173##[10]##. Indeed, the small transmembrane protein ERG28p functions as a scaffold to tether the C-4 demethylation complex (ERG25p-ERG26p-ERG27p) to the ER ##UREF##2##[11]## and also interacts with the downstream enzyme ERG6p. More recent results have indicated that ERG28p also interacts with ERG11p and ERG1p. Moreover, ERG27p is required for oxidosqualene cyclase (ERG7p) activity. The physical interaction of ERG27p with ERG7p might indeed contribute to yeast sterol biosynthesis regulation ##UREF##3##[12]##. These results altogether suggest that many sterol biosynthetic proteins, if not all, may be tethered to the ER as a large complex ##UREF##1##[9]##, ##REF##15995173##[10]##.</p> <p>Most animals synthesize cholesterol. However, some animals such as <italic>D. melanogaster</italic> and <italic>C. elegans</italic> cannot synthesize sterols <italic>de novo</italic>
##REF##4869093##[13]##, ##REF##1779708##[14]##. They are auxotrophic for sterols because they do not possess the enzymatic activities necessary to complete this process ##REF##4869093##[13]##, ##REF##1779708##[14]##. <italic>C. elegans</italic> takes sterols from animal feces or yeast and plant remnants ##REF##12894170##[3]##. <italic>Drosophila</italic> obtains sterols from the diet: ergosterol from yeast and phytosterols from plants ##REF##9632664##[15]##. However, both animals express the homologs of the enzymes that produce the initial intermediates of sterol biosynthesis up to the very early precursor farnesyl pyrophosphate. However, they cannot synthesize either squalene or lanosterol, key intermediates of sterol biosynthesis ##REF##4869093##[13]##, ##REF##1779708##[14]##.</p> <p>As already noted, sterols in these animals are required as hormone precursors and/or developmental effectors ##REF##12518036##[16]##. Not surprisingly, <italic>dauer larva</italic> formation and molting depend on sterols. Cholesterol, or more likely its derivatives, seem to act as hormones. Indeed, recent papers report the identification of natural steroid ligands for the DAF-12 nuclear receptor ##REF##16564008##[17]##.</p> <p>In <italic>C. elegans</italic>, the distribution and transport of cholesterol <italic>in vivo</italic> has been studied by using dehydroergosterol (DHE), a fluorescent analog which mimics many cholesterol properties ##UREF##4##[18]##, ##UREF##5##[19]##. DHE accumulates primarily in the pharynx, nerve ring, excretory gland cell, and gut of L1–L3 larvae ##REF##11408580##[20]##. Interestingly, sterols present in the pharynx and in the intestine of feeding animals are distributed in a proximal-to-distal gradient. Cholesterol in <italic>C. elegans</italic> might be involved in the structural and functional organization of the plasma membrane cell types that are richer in this lipid ##REF##10559886##[21]## and in modulating the activity of signaling molecules (such as Hh-like proteins).</p> <p>A previous comparative genome analysis of <italic>D. melanogaster</italic> with <italic>Anopheles gambiae</italic> and prototrophs has suggested that these insects have lost most of the genes involved in sterol synthesis ##REF##12364792##[22]##. This makes sense, knowing that <italic>Drosophila</italic> is unable to synthesize cholesterol. However, in a previous work we have shown that <italic>Drosophila</italic> contains an <italic>ERG28</italic> ortholog that has undergone a process of acceleration in its evolution, and is undetectable using the current techniques for ortholog detection by sequence homology ##REF##11910120##[23]##. Thus, here we have revisited this question for both <italic>C. elegans</italic> and <italic>D. melanogaster</italic>, not only for <italic>ERG28</italic>, but most of the genes/enzymes involved in the sterol synthesis pathway, in the light of new genomic and functional data.</p>
<title>Materials and Methods</title> <title>Ortholog search</title> <p>We used BLASTp ##REF##2231712##[24]## and considered as orthologs the best reciprocal hits. In short, a BLAST search is performed starting with a protein sequence from organism A to detect the matching sequences in organism B. Then, the sequence from B displaying the highest score is compared (reverse BLAST) against the all sequences of A and the highest scoring hit must be the initial sequence. This is a widely accepted criterion of orthology ##REF##16948844##[25]##. We also exploited the significance of matches using Psi-BLAST ##REF##9254694##[35]## or with the conserved domain database (CDD), where the scoring of the protein alignment uses very sensitive “tailor made” scoring matrices.</p> <title>Computation of hydrophobic profiles as further evidence of structural relationship between potentially orthologous proteins</title> <p>For some distant homologues we computed the hydrophobic profiles and calculated the Pearson correlation coefficient “R” for various pairs of profiles. The values of R range between 1 (perfect match between the profiles) and −1 (profiles are mirror images). Strong positive correlation is not a proof of orthology but strengthens the idea of structural relationship at the protein level.</p> <title>Calculation of Ka/Ks</title> <p>The protein-coding nucleotide sequences were first translated into amino acid sequences and aligned. Then, this alignment was used to define the alignment of the corresponding nucleotide sequences to avoid frame-shifting indels as alignment artifacts ##UREF##11##[46]##. Ka and Ks were estimated in two ways: i) by using the PBL method ##REF##8433381##[47]##, ##REF##8487630##[48]## implemented in DAMBE ##UREF##11##[46]##, ##REF##11535656##[49]##, and ii) by the likelihood-based method implemented in the YN00 program in the PAML package with the resulting Ka and Ks designated by dN and dS, respectively in ##TAB##0##table 1##
##REF##10666704##[50]##.</p> <title>pI estimation and outlier detection</title> <p>We have investigated whether the pI of presumably orthologous proteins were similar. In order to test this, we gathered the protein sequences of the orthologs of ERG2, 6, 24, 25, 26 and 28 listed in the HomoloGene division of the NCBI database and estimated their pI (<ext-link ext-link-type="uri" xlink:href="http://www.expasy.ch/tools/pi_tool.html">http://www.expasy.ch/tools/pi_tool.html</ext-link>). Then, we asked whether the corresponding <italic>Drosophila</italic> and <italic>Caenorhabditis</italic> sequences displayed outlier pI values. In short, we performed an extreme studentized deviate test to determine whether one of the values in the pI list was a significant outlier from the rest.</p>
<title>Results and Discussion</title> <title>Looking for <italic>ERG</italic> orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic>
</title> <p>In this section we present our search for <italic>ERG</italic> gene orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic> following their order in the sterol synthesis cascade (as shown in ##FIG##0##Figure 1##). We have taken advantage of the fact that the full genomic sequences of these two animals are available. We have used BLASTp ##REF##2231712##[24]## and considered as orthologs the best reciprocal hits ##REF##16948844##[25]##. The results are summarized in ##FIG##0##Figures 1## and ##FIG##1##2##.</p> <p>We failed to detect the squalene synthase (ERG9p) homolog, which catalyzes the first committed step in cholesterogenesis. In the search for squalene epoxidase (ERG1p) in a BLAST using ERG1p from yeast, we detected the gene CBG19254 in <italic>C. briggsae</italic>, with a marginal Score (S = 48 bits, E-value = 10<sup>−04</sup>). This protein, which also exists in <italic>C. elegans</italic>, contains two functional domains: monooxygenase and UbiH (2-polyprenyl-6-methoxyphenol hydroxylase and FAD-dependent oxidase). This is also the case in yeast ERG1p. However, in a reverse BLAST, CBG19254 recognized with very high score yeast Coq6 and only marginally ERG1p (S = 48, E = 10<sup>−06</sup>). A similar behavior was observed for GA20231-RA from <italic>D. pseudoobscura</italic>, so we did not proceed to a further analysis of these sequences.</p> <p>For lanosterol synthase (ERG7p), which follows in the classical pathway, we could not gather any convincing evidence for the existence of orthologs in the fruitfly and the nematode.</p> <p>The following enzyme in the pathway is a sterol 14α-demethylase, ERG11p/Cyp51, involved in the biosynthesis of cholesterol, phytosterols and ergosterol. Thus, it is the only cytochrome P450 having an ortholog common to animals, plant and fungi ##REF##11179683##[26]##. Similarity BLAST hits with yeast ERG11p/Cyp51, were obtained in <italic>Drosophila</italic> (CG2397, CG10247 and other Cyps). However, the <italic>Drosophila</italic> genes are likely to belong to other Cyp subfamilies (not Cyp51). Cyp51 is probably missing, which is in agreement with the results of Tijet, Helvig &amp; Feyereisen ##REF##11179683##[26]## who analyzed 90 sequences of the cytochrome <italic>P450</italic> gene superfamily. Cyp51 is also absent in <italic>C. elegans</italic>
##REF##11179683##[26]##.</p> <p>Potential <italic>ERG24</italic> orthologs in <italic>D. melanogaster</italic> (CG17952) and in <italic>C. elegans</italic> (B0250.9) were easily found. The corresponding proteins contain the ERG4-24 domain. <italic>D. melanogaster</italic> produces three isoforms that are longer than the yeast ortholog, a peculiarity that they share with the human ortholog. The ortholog of <italic>ERG24</italic> in mammals encodes the Lamin B receptor (LBR), a nuclear envelope protein first described in vertebrates. LBR bears extensive structural similarities with the members of the sterol reductase family (ERG24p and ERG4p). Human LBR (hLBR) cannot restore ergosterol biosynthesis in an <italic>ERG4</italic> yeast mutant, whereas it is able to restore ergosterol prototrophy in an <italic>ERG24</italic> mutant. This strongly suggests that hLBR is a sterol C14-reductase ##UREF##6##[27]##. Not surprisingly, a mutation in the <italic>hLBR</italic> gene causes an autosomal recessive disease called hydrops-ectopic calcification-‘moth-eaten’ (HEM). This mutation leads to high levels of cholesta-8,14-dien-3-beta-ol in cultured skin fibroblasts, which is compatible with a deficiency of the cholesterol biosynthetic enzyme 3-beta-hydroxysterol delta(14)-reductase ##REF##12618959##[28]##.</p> <p>The hLBR contains two major domains: a ∼220-amino-acid N-terminal segment highly charged, and a hydrophobic C-terminal half with eight putative transmembrane segments ##UREF##7##[29]##, ##REF##9169472##[30]##. Interestingly, it has been hypothesized that the region encoding the N-terminal domain of the <italic>LBR</italic> gene arose from an ancestral gene coding for a soluble nuclear protein (which provides a nuclear localization signal) and that the rest of the protein evolved from another gene, similar to yeast <italic>ERG24</italic>. Indeed, the C-terminal hydrophobic domain of LBR can be retained in the endoplasmic reticulum when expressed in transfected cells, as expected for the ortholog of <italic>ERG24</italic> in mammals. In turn, the N-terminal domain is transported to the nucleus ##REF##7679672##[31]##. This domain might be responsible for the targeting of the hydrophobic domain to the nucleus and for the interaction with lamin B ##REF##8157663##[32]##. So far, cholesterol synthesis is supposed to occur in the smooth ER. Since the N-terminal domain of LBR is responsible for its nuclear localization, it would be interesting to investigate whether the LBR transcript, or protein, undergoes some processing leading to the production of a C-terminal domain sorted to the ER.</p> <p>Recent functional studies show that the <italic>Drosophila</italic> CG17952 gene is the ortholog of vertebrate <italic>LBR</italic>
##REF##15054108##[33]##. The protein encoded by CG17952 shares some properties with hLBR. The <italic>Drosophila</italic> LBR (dLBR) possesses a highly charged N-terminal domain of 307 amino acids followed by eight transmembrane segments. Transmembrane segments 1–6 are similar in length and position to the transmembrane domains 1–6 of hLBR. However, the putative membrane domains 7 and 8 of dLBR are shorter than those of hLBR. Thus, dLBR is expected to have a topological organization similar to that of its vertebrate orthologs ##REF##15054108##[33]##. dLBR is able to bind to the <italic>Drosophila</italic> lamin B, a function residing in the N-terminal domain. Not unexpectedly, dLBR does not display sterol C14 reductase activity when expressed in the yeast <italic>ERG24</italic> mutant. This shows that, during insect evolution, although the enzymatic activity of this protein has been lost, its capacity to bind lamin B has not. However, depletion of dLBR by RNA interference does not lead to any obvious effect on nuclear architecture, or viability, of treated cells and embryos. Thus, although dLBR might be important, it is not a limiting component of the nuclear architecture in <italic>Drosophila</italic> cells, at least during the first days of development ##REF##15054108##[33]##. Our BLAST search shows that sequence B0250.9 is the potential <italic>ERG24</italic> ortholog in <italic>C. elegans</italic>. It would be interesting to experimentally assess if it has kept LBR activity.</p> <p>Sequences homologous to ERG25p (C-4 methyloxidase) were also easily found in both <italic>D. melanogaster</italic> (CG1998/<italic>dERG25A</italic> and CG11162/<italic>dERG25B</italic>) and <italic>C. elegans</italic> (F49E12.9/<italic>ERG25A</italic> and F49E12.10<italic>/ERG25B</italic>) using the sequence of yeast ERG25p as a starting point. In both organisms the duplicated copies of <italic>ERG25</italic> are located in the same chromosome (chromosome II for <italic>C. elegans</italic> and chromosome X for <italic>D. melanogaster</italic>). The two paralogs in <italic>D. melanogaster</italic> are separated by 0.25 Mb and contain a different number of predicted exons. Namely, CG1998 contain 6 exons, while CG11162 contains only 2. However, the last intron of both genes interrupts the coding sequence at very similar positions (i.e. between the second and the third positions of a Lys codon, ##FIG##2##Figure 3##). In <italic>C. elegans</italic> the paralogs are at less than 3 kb away from each other. F49E12.9 contains 8 exons while F49E12.10 has 5 exons that may have been produced through exon fusion/splitting events.</p> <p>ERG26p (C-3 dehydrogenase) belongs to the 3β-hydroxysteroid dehydrogenase family and convincing evidence for the presence of orthologs in <italic>C. elegans</italic> (ZC8.1) and <italic>D. melanogaster</italic> (CG7724) was obtained in BLAST searches taking the sequence from yeast as the starting point (See details in ##FIG##1##Figure 2##). In the case of <italic>Drosophila</italic> CG7724, similarity with yeast ERG26p extends over the first 250 amino acids while the remaining amino acids are more divergent. Interestingly, when performing a BLAST with the <italic>Drosophila</italic> sequence, the first hit in <italic>C. elegans</italic> was C32D5.12, but it was not the best hit either with yeast, or with <italic>A. thaliana</italic>, or even with the human (NSDHL) orthologs. However, a BLAST with C32D5.12 detected as the first hits NSDHL in man and ERG26 in yeast (##FIG##1##Figure 2##). Thus, it is tempting to invoke some kind of sequence convergence as an explanation for this behavior.</p> <p>In agreement with Breitling <italic>et al.</italic>
##UREF##8##[34]## we failed to find any clear homologue of ERG27p (C-3 ketoreductase) in both <italic>Drosophila</italic> and <italic>C. elegans</italic>, although several oxidoreductases were detected.</p> <p>As outlined above, ERG28p might tether many other ERG proteins to the ER. The ERG28p ortholog of <italic>C. elegans</italic> (C14C10.6) was hardly detectable by BLAST starting with yeast sequences. This precludes the use of standard phylogenetics methods to show orthology. However, further evidence of sequence relatedness was gathered using Psi-BLAST with the yeast sequence against the Metazoa division of Genbank ##REF##9254694##[35]##. We included in the iterations the very divergent sequence CG17270 of <italic>D. melanogaster</italic>, which is the ortholog of ERG28p according to our previous results ##REF##11910120##[23]##. This allowed us to detect C14C10.6 as a signinficant hit (S = 120 bits and E = 10<sup>−27</sup>). Reverse Psi-BLAST also suggested orthology (i.e. significant scores). The sequence of C14C10.6 is so divergent that it had no match in the conserved domain database (CDD). Moreover, we also computed the hydrophobic profiles for some of the orthologs ##REF##7704669##[36]## and calculated the Pearson correlation coefficient “R” for various pairs of profiles. Significant R-values were obtained for the pair-wise comparisons (##FIG##3##Figure 4##). This is not a proof of orthology but strengthens the idea of structural relationship at the protein level. Finally, we found that all proteins had similar lengths and were basic, with isoelectric points (pI) &gt;8.5.</p> <p>The ortholog of ERG6p in <italic>C. elegans</italic> (H14E04.1) was easily detected by BLAST with protein sequences from either yeast or <italic>A. thaliana</italic> (SMT1). We used the plant sequence since no clear <italic>ERG6</italic> homolog could be detected in human. The issue with <italic>Drosophila</italic> turned out to be more complex because, when starting the search with the yeast sequence, we detected CG8067 marginally (E&gt;1). This was even worse when starting with the sequence of <italic>A. thaliana</italic>. However, when using the <italic>C. elegans</italic> sequence as a starting point, the first significant hit in <italic>Drosophila</italic> (E&lt;10<sup>−6</sup>) was CG2453, which proved to be the ortholog of yeast Coq5, but not of ERG6p. Finally, considering i) the similar lengths of the previously marginally detected CG8067, of H14E04.1 and of SMT1 proteins and ii) their similar pI, we have preferred CG8067 as the most likely ortholog. Indeed, the hydrophobic profiles of SMT1 and the protein encoded by CG8067 displayed a strong correlation. Namely, we obtained an R = 0.43 with a p-value 10<sup>−16</sup>.</p> <p>The ortholog of ERG2p in <italic>C. elegans</italic> (W08F4.3) was marginally detected by BLAST with the yeast and the human protein (opioid sigma-1 receptor, OPRS1). However, W08F4.3 was found to contain a sigma1-receptor domain when compared with the CDD. This strengthens the idea that this gene is the ortholog of <italic>OPRS1</italic> and <italic>ERG2</italic>. The situation in <italic>Drosophila</italic> was more complex. When starting our BLAST with either yeast or human sequences, we detected the sequence HDC14735 (DAA04220) very marginally (E≫1) and no conserved domain was found. However, when starting with the <italic>C. elegans</italic> sequence, it came as the best hit with E = 0.002 (the reverse was also true, with E = 0.003). Although not significant, this result was taken as suggestive of similarity. The results were improved using Psi-BLAST. Again, we computed the hydrophobic profiles for the various potential orthologs and we found strong correlations (##FIG##3##Figure 4##). Finally, the pI of the protein encoded by HDC14735 (pI = 6.14) was comparable to those of OPRS1 (5.61) and ERG2 (5.54) (see below).</p> <p>In BLAST searches with yeast ERG3p, we detected again CG1998 and CG11162 in <italic>D. melanogaster</italic> and F49E12.9 and F49E12.10 in <italic>C. elegans</italic>, but with worse scores than in BLASTs with ERG25 (S∼55 bits versus 85 bits respectively). Therefore, we propose that the orthologs of ERG3 are potentially missing in both organisms.</p> <p>For ERG5p, which belongs to the big Cyp protein family, no clear orthologs could be established. However, three potential candidates were found: CG4321-PA (Cyp4d8), CG3540-PA (Cyp4d14) and CG8859-PA (Cyp6g2). In the reverse BLAST, they all matched ERG5p as the best scoring hit in yeast.</p> <p>Finally, the search for ERG4p orthologs led to the same ERG24 orthologs in the nematode and the fruitfly. Moreover, with ERG4p the BLAST scores were worse than with ERG24p. Thus, either ERG4 orthologs are missing or they have been replaced by ERG24.</p> <title>Evolution of the <italic>ERG</italic> orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic>
</title> <p>To assess whether some of the <italic>ERG</italic> orthologs mentioned above undergo a selective pressure, we have examined the ratio of the number of non-synonymous substitutions per non-synonymous site (Ka) and the number of synonymous substitutions per synonymous site (Ks). The ratio Ka/Ks is indicative of the mode of evolution operating on the sequences. If selection is dominantly purifying, then we expect few non-synonymous substitutions per background synonymous changes and hence, a low ratio. If selection is absent, then a ratio of unity is expected (##TAB##0##Table 1##). To obtain the Ka/Ks values for the fly we compared the sequences from <italic>D. melanogaster</italic> and <italic>D. pseudoobscura</italic> while for the nematode we compared the sequences from <italic>C. elegans</italic> and <italic>C. briggsae</italic>. In spite of the divergence of the various <italic>ERG</italic> orthologs in the lineages leading to insects and to nematodes, the low value of Ka/Ks (dN/dS) ratios show that these genes are under selective pressure (##TAB##0##Table 1##).</p> <p>It is natural to ask what protein properties have been conserved by the purifying selection. We have investigated a protein-level property, namely, the isoelectric point (pI). This property is important for enzymes because protein-protein and enzyme-substrate interactions are often electrostatic in nature. Thus, we should expect the pI of an enzyme to be similar in different organisms if protein-protein and/or enzyme-substrate interactions are conserved. In order to test this idea we gathered the protein sequences of the orthologs of ERG2, 6, 24, 25, 26 and 28. Next, we asked whether the corresponding <italic>Drosophila</italic> and <italic>Caenorhabditis</italic> sequences displayed outlier pI values. None of the orthologous proteins in the auxotrophs was detected as outlier, suggesting conservation of this physico-chemical property within the nematode and the fly lineages in spite of primary sequence divergence (##TAB##1##Table 2##). Similarity in pI is expected, when the orthologs are detected as the best reciprocal significant hits in BLAST searches. However, the existence of similarity in the hydrophobic profiles and the pI between the potential orthologs, for which the BLAST comparisons failed to be significant, is more surprising. It means that, in spite of important sequence divergence, which renders BLAST unreliable in many cases, the proteins have maintained substantial structural and physico-chemical similarity.</p> <title>Analyzing expression-profiling data to explore the potential function to the divergent <italic>ERG</italic> orthologs</title> <p>Co-expression is indicative of i) physical interaction between proteins and ii) of membership to the same complex or molecular process ##REF##11779829##[37]##. This is called the paradigm of “guilt by association” ##REF##14551426##[38]##. Thus, we have used published microarray expression data, downloadable from the Gene Expression Omnibus of the NCBI, to investigate potential co-expression patterns of the <italic>ERG</italic> orthologs that we have described above for both <italic>D. melanogaster</italic> and <italic>C. elegans</italic>. Co-expression can be assessed by determining the correlation coefficient. The correlation coefficient can be artificially inflated by flat profiles (no changes in the expression of the relevant genes). To avoid this, we focused on experiments where the genes of interest display strong variation (see ##FIG##4##Figure 5A and B##). Thus, we gathered data concerning 51 different microarray experiments for <italic>D. melanogaster</italic> respecting the criterion outlined previously. We performed a similar analysis for <italic>C. elegans</italic> but, unfortunately, the most interesting genes showed flat profiles (close to 0 in all experiments) and it was not possible to proceed further with the analysis.</p> <p>First, we asked whether the expression profiles of the <italic>D. melanogaster ERG</italic> orthologs were correlated (##TAB##2##Table 3##). The strongest correlation was found between the <italic>dLBR</italic> (<italic>ERG24</italic>) and CG1998 (<italic>ERG25</italic>) with an associated p-value of 10<sup>−16</sup> (after a Bonferroni correction). Such a p-value means that only one correlation coefficient out of 10<sup>16</sup> is expected to be as high as 0.89 just by chance (for n = 51 experimental points). Considering the maximum number of possible correlations for the 14000 transcripts in the microarrays (representing the <italic>Drosophila</italic> genome), such a high R cannot be found by chance. The behavior of dLBR and CG1998 might be reminiscent of the situation in yeast because ERG24p and ERG25p are supposed to interact, according to Epistatic MiniArray Profiling experiments ##REF##16269340##[39]##. CG7724 and CG11162 displayed the second highest R (R = 0.4, p = 0.03), but this R is not relevant in genomic terms.</p> <p>At first, we were expecting good expressional correlation among the <italic>ERG</italic> orthologs in <italic>Drosophila</italic>. However, it seems that only dLBR and CG1998 still “remember” their ancestral belonging to the sterol biosynthesis pathway. Poor expressional correlation among the rest of the <italic>ERG</italic> orthologs also suggests that the corresponding proteins either have lost their ability to physically interact in order to form stable complexes, or they do so in conditions/moments not covered by the microarray experiments explored here. Then, we focused our attention on <italic>dLBR</italic> and CG1998 by determining which other genes were expressionally correlated with them. For this, we gathered 84 genes displaying R≥0.875 with respect to both genes. For a correlation involving 51 data points, this R cut-off is associated with a safe p-value of 10<sup>−13</sup> after correction (##FIG##4##Figure 5C##). In order to get insights about these 84 genes, we used the functional classification tool of the DAVID database (<ext-link ext-link-type="uri" xlink:href="http://david.abcc.ncifcrf.gov/">http://david.abcc.ncifcrf.gov/</ext-link>
##REF##12734009##[40]##). This software provides a rapid means to organize large lists of genes into functionally related groups and to unravel biological relationships.</p> <p>In the analysis using DAVID, the most overrepresented class included genes encoding membrane proteins, often targeted to the ER, where sterol biosynthesis takes place in prototrophs (Group 1, in ##TAB##3##Table 4##). Interestingly, several of these proteins are supposed to be involved in co-translational protein targeting to membranes, signal peptide recognition, heat shock protein-binding, as well as unfolded protein binding, or to be elements of the translocon (a complex of proteins associated with the translocation of nascent polypeptides into the ER ##REF##15130556##[41]##). The following functional category (Group 2) contained four transporter proteins while the last functional group involved chaperones (i.e. peptidyl-prolyl <italic>cis-trans</italic> isomerase), chaperone cofactors or unfolded protein-binding factors. A similar analysis conducted using g∶profiler ##REF##17478515##[42]## confirmed that genes encoding protein folding actors were overrepresented among the genes displaying strong expressional correlation with dLBR and CG1998 (p&lt;10<sup>−5</sup>). The existence of expressional correlation does not imply any causality. In fact, from this exploration it is not possible to determine whether dLBR and CG1998 somehow interact with other partners to participate in intracellular protein trafficking or folding, or on the contrary, they undergo the action of the latter. Since dLBR has been shown to be a nuclear protein ##REF##7679672##[31]##, ##REF##8157663##[32]##, it would be interesting to investigate whether the dLBR transcript or protein are somehow processed to produce a C-terminal polypeptide that might be sorted to the ER. That would explain the expressional correlation between dLBR and ER proteins observed above. On the other hand, we have also explored the annotation of the genes expressionally correlated with the rest of the <italic>ERG</italic> orthologs. However, no unifying theme emerged from this analysis (data not shown, available upon request). All in all, the strong expressional correlation between dLBR and CG1998 with proteins involved in intracellular protein trafficking or folding, and the absence of such correlation with other Erg orthologs (that also require chaperons) suggest that the involvement of dLBR and CG1998 in both processes is worth exploring.</p> <p>In a previous paper, given the structural similarity between cholesterol and ecdysteroids, we had proposed that divergent ERGp orthologs might somehow participate in the synthesis of the latter ##REF##11910120##[23]##. We have therefore assessed the expressional correlation between candidate genes involved in this process: <italic>Dare1</italic>/CG12390 ##REF##10498693##[43]##, <italic>Jhamt</italic>
##UREF##9##[44]## and <italic>Start1</italic>
##UREF##10##[45]##, with <italic>dLBR</italic> and <italic>CG1998</italic>. While for <italic>Jahmt</italic> and <italic>Dare1</italic> the values of R are below 0.6, a very strong correlation (R = 0.8) was found for <italic>Start1</italic>. Interestingly, <italic>Start1</italic>, which is involved in intra-mitochondrial sterol transport, is expressed ubiquitously. However, <italic>in situ</italic> hybridization demonstrates a stronger expression in the prothoracic gland, where ecdysteroids are synthesized from cholesterol. These and other observations are consistent with the idea that Start1 plays a key role in the regulation of ecdysteroid synthesis ##UREF##10##[45]##. The potential functional link between dLBR, CG1998 and Start1 is also worth exploring.</p> <p>In conclusion, we detected a preservation of <italic>ERG</italic> genes in <italic>Drosophila melanogaster</italic> and <italic>Caenorhabditis elegans</italic>. In spite of their sequence divergence with respect to the corresponding orthologs in sterol prototrophs, they still are under selective pressure. Since insects are unable to synthesize cholesterol <italic>de novo</italic>, an appealing way to explain this evolutionary acceleration is that ERGp orthologs have other biological functions in addition to sterol synthesis. This is clearly the case of the LBR, which is also a reductase in sterol prototrophs. Shut-down of cholesterogenesis in insects and nematodes would have allowed these proteins to evolve as much as their other functions were not compromised ##REF##11910120##[23]##. Another, less parsimonious, explanation would be the evolution of different novel functions. Our microarray meta-analysis shows strong expressional correlation between the orthologs of <italic>ERG24</italic> and <italic>ERG25</italic> in <italic>D. melanogaster</italic> and genes encoding factors involved in intracellular protein trafficking and folding. This is compatible with our idea that ERGp might be involved in other biological roles in addition to sterol synthesis. The potential link between ERG proteins and intracellular protein trafficking and folding deserves experimental exploration not only in <italic>Drosophila</italic> and <italic>Caenorhabditis</italic> but also in sterol prototrophs. Moreover, the potential link between dLBR, CG1998 and Start1 is to be explored in <italic>D. melanogaster</italic>. This is compatible with our previous idea of a potential implication of these proteins in the synthesis of ecdysteroids. We hope that this genomic exploration and the hypotheses prompted here might open new avenues of experimental research.</p>
<title>Results and Discussion</title> <title>Looking for <italic>ERG</italic> orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic>
</title> <p>In this section we present our search for <italic>ERG</italic> gene orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic> following their order in the sterol synthesis cascade (as shown in ##FIG##0##Figure 1##). We have taken advantage of the fact that the full genomic sequences of these two animals are available. We have used BLASTp ##REF##2231712##[24]## and considered as orthologs the best reciprocal hits ##REF##16948844##[25]##. The results are summarized in ##FIG##0##Figures 1## and ##FIG##1##2##.</p> <p>We failed to detect the squalene synthase (ERG9p) homolog, which catalyzes the first committed step in cholesterogenesis. In the search for squalene epoxidase (ERG1p) in a BLAST using ERG1p from yeast, we detected the gene CBG19254 in <italic>C. briggsae</italic>, with a marginal Score (S = 48 bits, E-value = 10<sup>−04</sup>). This protein, which also exists in <italic>C. elegans</italic>, contains two functional domains: monooxygenase and UbiH (2-polyprenyl-6-methoxyphenol hydroxylase and FAD-dependent oxidase). This is also the case in yeast ERG1p. However, in a reverse BLAST, CBG19254 recognized with very high score yeast Coq6 and only marginally ERG1p (S = 48, E = 10<sup>−06</sup>). A similar behavior was observed for GA20231-RA from <italic>D. pseudoobscura</italic>, so we did not proceed to a further analysis of these sequences.</p> <p>For lanosterol synthase (ERG7p), which follows in the classical pathway, we could not gather any convincing evidence for the existence of orthologs in the fruitfly and the nematode.</p> <p>The following enzyme in the pathway is a sterol 14α-demethylase, ERG11p/Cyp51, involved in the biosynthesis of cholesterol, phytosterols and ergosterol. Thus, it is the only cytochrome P450 having an ortholog common to animals, plant and fungi ##REF##11179683##[26]##. Similarity BLAST hits with yeast ERG11p/Cyp51, were obtained in <italic>Drosophila</italic> (CG2397, CG10247 and other Cyps). However, the <italic>Drosophila</italic> genes are likely to belong to other Cyp subfamilies (not Cyp51). Cyp51 is probably missing, which is in agreement with the results of Tijet, Helvig &amp; Feyereisen ##REF##11179683##[26]## who analyzed 90 sequences of the cytochrome <italic>P450</italic> gene superfamily. Cyp51 is also absent in <italic>C. elegans</italic>
##REF##11179683##[26]##.</p> <p>Potential <italic>ERG24</italic> orthologs in <italic>D. melanogaster</italic> (CG17952) and in <italic>C. elegans</italic> (B0250.9) were easily found. The corresponding proteins contain the ERG4-24 domain. <italic>D. melanogaster</italic> produces three isoforms that are longer than the yeast ortholog, a peculiarity that they share with the human ortholog. The ortholog of <italic>ERG24</italic> in mammals encodes the Lamin B receptor (LBR), a nuclear envelope protein first described in vertebrates. LBR bears extensive structural similarities with the members of the sterol reductase family (ERG24p and ERG4p). Human LBR (hLBR) cannot restore ergosterol biosynthesis in an <italic>ERG4</italic> yeast mutant, whereas it is able to restore ergosterol prototrophy in an <italic>ERG24</italic> mutant. This strongly suggests that hLBR is a sterol C14-reductase ##UREF##6##[27]##. Not surprisingly, a mutation in the <italic>hLBR</italic> gene causes an autosomal recessive disease called hydrops-ectopic calcification-‘moth-eaten’ (HEM). This mutation leads to high levels of cholesta-8,14-dien-3-beta-ol in cultured skin fibroblasts, which is compatible with a deficiency of the cholesterol biosynthetic enzyme 3-beta-hydroxysterol delta(14)-reductase ##REF##12618959##[28]##.</p> <p>The hLBR contains two major domains: a ∼220-amino-acid N-terminal segment highly charged, and a hydrophobic C-terminal half with eight putative transmembrane segments ##UREF##7##[29]##, ##REF##9169472##[30]##. Interestingly, it has been hypothesized that the region encoding the N-terminal domain of the <italic>LBR</italic> gene arose from an ancestral gene coding for a soluble nuclear protein (which provides a nuclear localization signal) and that the rest of the protein evolved from another gene, similar to yeast <italic>ERG24</italic>. Indeed, the C-terminal hydrophobic domain of LBR can be retained in the endoplasmic reticulum when expressed in transfected cells, as expected for the ortholog of <italic>ERG24</italic> in mammals. In turn, the N-terminal domain is transported to the nucleus ##REF##7679672##[31]##. This domain might be responsible for the targeting of the hydrophobic domain to the nucleus and for the interaction with lamin B ##REF##8157663##[32]##. So far, cholesterol synthesis is supposed to occur in the smooth ER. Since the N-terminal domain of LBR is responsible for its nuclear localization, it would be interesting to investigate whether the LBR transcript, or protein, undergoes some processing leading to the production of a C-terminal domain sorted to the ER.</p> <p>Recent functional studies show that the <italic>Drosophila</italic> CG17952 gene is the ortholog of vertebrate <italic>LBR</italic>
##REF##15054108##[33]##. The protein encoded by CG17952 shares some properties with hLBR. The <italic>Drosophila</italic> LBR (dLBR) possesses a highly charged N-terminal domain of 307 amino acids followed by eight transmembrane segments. Transmembrane segments 1–6 are similar in length and position to the transmembrane domains 1–6 of hLBR. However, the putative membrane domains 7 and 8 of dLBR are shorter than those of hLBR. Thus, dLBR is expected to have a topological organization similar to that of its vertebrate orthologs ##REF##15054108##[33]##. dLBR is able to bind to the <italic>Drosophila</italic> lamin B, a function residing in the N-terminal domain. Not unexpectedly, dLBR does not display sterol C14 reductase activity when expressed in the yeast <italic>ERG24</italic> mutant. This shows that, during insect evolution, although the enzymatic activity of this protein has been lost, its capacity to bind lamin B has not. However, depletion of dLBR by RNA interference does not lead to any obvious effect on nuclear architecture, or viability, of treated cells and embryos. Thus, although dLBR might be important, it is not a limiting component of the nuclear architecture in <italic>Drosophila</italic> cells, at least during the first days of development ##REF##15054108##[33]##. Our BLAST search shows that sequence B0250.9 is the potential <italic>ERG24</italic> ortholog in <italic>C. elegans</italic>. It would be interesting to experimentally assess if it has kept LBR activity.</p> <p>Sequences homologous to ERG25p (C-4 methyloxidase) were also easily found in both <italic>D. melanogaster</italic> (CG1998/<italic>dERG25A</italic> and CG11162/<italic>dERG25B</italic>) and <italic>C. elegans</italic> (F49E12.9/<italic>ERG25A</italic> and F49E12.10<italic>/ERG25B</italic>) using the sequence of yeast ERG25p as a starting point. In both organisms the duplicated copies of <italic>ERG25</italic> are located in the same chromosome (chromosome II for <italic>C. elegans</italic> and chromosome X for <italic>D. melanogaster</italic>). The two paralogs in <italic>D. melanogaster</italic> are separated by 0.25 Mb and contain a different number of predicted exons. Namely, CG1998 contain 6 exons, while CG11162 contains only 2. However, the last intron of both genes interrupts the coding sequence at very similar positions (i.e. between the second and the third positions of a Lys codon, ##FIG##2##Figure 3##). In <italic>C. elegans</italic> the paralogs are at less than 3 kb away from each other. F49E12.9 contains 8 exons while F49E12.10 has 5 exons that may have been produced through exon fusion/splitting events.</p> <p>ERG26p (C-3 dehydrogenase) belongs to the 3β-hydroxysteroid dehydrogenase family and convincing evidence for the presence of orthologs in <italic>C. elegans</italic> (ZC8.1) and <italic>D. melanogaster</italic> (CG7724) was obtained in BLAST searches taking the sequence from yeast as the starting point (See details in ##FIG##1##Figure 2##). In the case of <italic>Drosophila</italic> CG7724, similarity with yeast ERG26p extends over the first 250 amino acids while the remaining amino acids are more divergent. Interestingly, when performing a BLAST with the <italic>Drosophila</italic> sequence, the first hit in <italic>C. elegans</italic> was C32D5.12, but it was not the best hit either with yeast, or with <italic>A. thaliana</italic>, or even with the human (NSDHL) orthologs. However, a BLAST with C32D5.12 detected as the first hits NSDHL in man and ERG26 in yeast (##FIG##1##Figure 2##). Thus, it is tempting to invoke some kind of sequence convergence as an explanation for this behavior.</p> <p>In agreement with Breitling <italic>et al.</italic>
##UREF##8##[34]## we failed to find any clear homologue of ERG27p (C-3 ketoreductase) in both <italic>Drosophila</italic> and <italic>C. elegans</italic>, although several oxidoreductases were detected.</p> <p>As outlined above, ERG28p might tether many other ERG proteins to the ER. The ERG28p ortholog of <italic>C. elegans</italic> (C14C10.6) was hardly detectable by BLAST starting with yeast sequences. This precludes the use of standard phylogenetics methods to show orthology. However, further evidence of sequence relatedness was gathered using Psi-BLAST with the yeast sequence against the Metazoa division of Genbank ##REF##9254694##[35]##. We included in the iterations the very divergent sequence CG17270 of <italic>D. melanogaster</italic>, which is the ortholog of ERG28p according to our previous results ##REF##11910120##[23]##. This allowed us to detect C14C10.6 as a signinficant hit (S = 120 bits and E = 10<sup>−27</sup>). Reverse Psi-BLAST also suggested orthology (i.e. significant scores). The sequence of C14C10.6 is so divergent that it had no match in the conserved domain database (CDD). Moreover, we also computed the hydrophobic profiles for some of the orthologs ##REF##7704669##[36]## and calculated the Pearson correlation coefficient “R” for various pairs of profiles. Significant R-values were obtained for the pair-wise comparisons (##FIG##3##Figure 4##). This is not a proof of orthology but strengthens the idea of structural relationship at the protein level. Finally, we found that all proteins had similar lengths and were basic, with isoelectric points (pI) &gt;8.5.</p> <p>The ortholog of ERG6p in <italic>C. elegans</italic> (H14E04.1) was easily detected by BLAST with protein sequences from either yeast or <italic>A. thaliana</italic> (SMT1). We used the plant sequence since no clear <italic>ERG6</italic> homolog could be detected in human. The issue with <italic>Drosophila</italic> turned out to be more complex because, when starting the search with the yeast sequence, we detected CG8067 marginally (E&gt;1). This was even worse when starting with the sequence of <italic>A. thaliana</italic>. However, when using the <italic>C. elegans</italic> sequence as a starting point, the first significant hit in <italic>Drosophila</italic> (E&lt;10<sup>−6</sup>) was CG2453, which proved to be the ortholog of yeast Coq5, but not of ERG6p. Finally, considering i) the similar lengths of the previously marginally detected CG8067, of H14E04.1 and of SMT1 proteins and ii) their similar pI, we have preferred CG8067 as the most likely ortholog. Indeed, the hydrophobic profiles of SMT1 and the protein encoded by CG8067 displayed a strong correlation. Namely, we obtained an R = 0.43 with a p-value 10<sup>−16</sup>.</p> <p>The ortholog of ERG2p in <italic>C. elegans</italic> (W08F4.3) was marginally detected by BLAST with the yeast and the human protein (opioid sigma-1 receptor, OPRS1). However, W08F4.3 was found to contain a sigma1-receptor domain when compared with the CDD. This strengthens the idea that this gene is the ortholog of <italic>OPRS1</italic> and <italic>ERG2</italic>. The situation in <italic>Drosophila</italic> was more complex. When starting our BLAST with either yeast or human sequences, we detected the sequence HDC14735 (DAA04220) very marginally (E≫1) and no conserved domain was found. However, when starting with the <italic>C. elegans</italic> sequence, it came as the best hit with E = 0.002 (the reverse was also true, with E = 0.003). Although not significant, this result was taken as suggestive of similarity. The results were improved using Psi-BLAST. Again, we computed the hydrophobic profiles for the various potential orthologs and we found strong correlations (##FIG##3##Figure 4##). Finally, the pI of the protein encoded by HDC14735 (pI = 6.14) was comparable to those of OPRS1 (5.61) and ERG2 (5.54) (see below).</p> <p>In BLAST searches with yeast ERG3p, we detected again CG1998 and CG11162 in <italic>D. melanogaster</italic> and F49E12.9 and F49E12.10 in <italic>C. elegans</italic>, but with worse scores than in BLASTs with ERG25 (S∼55 bits versus 85 bits respectively). Therefore, we propose that the orthologs of ERG3 are potentially missing in both organisms.</p> <p>For ERG5p, which belongs to the big Cyp protein family, no clear orthologs could be established. However, three potential candidates were found: CG4321-PA (Cyp4d8), CG3540-PA (Cyp4d14) and CG8859-PA (Cyp6g2). In the reverse BLAST, they all matched ERG5p as the best scoring hit in yeast.</p> <p>Finally, the search for ERG4p orthologs led to the same ERG24 orthologs in the nematode and the fruitfly. Moreover, with ERG4p the BLAST scores were worse than with ERG24p. Thus, either ERG4 orthologs are missing or they have been replaced by ERG24.</p> <title>Evolution of the <italic>ERG</italic> orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic>
</title> <p>To assess whether some of the <italic>ERG</italic> orthologs mentioned above undergo a selective pressure, we have examined the ratio of the number of non-synonymous substitutions per non-synonymous site (Ka) and the number of synonymous substitutions per synonymous site (Ks). The ratio Ka/Ks is indicative of the mode of evolution operating on the sequences. If selection is dominantly purifying, then we expect few non-synonymous substitutions per background synonymous changes and hence, a low ratio. If selection is absent, then a ratio of unity is expected (##TAB##0##Table 1##). To obtain the Ka/Ks values for the fly we compared the sequences from <italic>D. melanogaster</italic> and <italic>D. pseudoobscura</italic> while for the nematode we compared the sequences from <italic>C. elegans</italic> and <italic>C. briggsae</italic>. In spite of the divergence of the various <italic>ERG</italic> orthologs in the lineages leading to insects and to nematodes, the low value of Ka/Ks (dN/dS) ratios show that these genes are under selective pressure (##TAB##0##Table 1##).</p> <p>It is natural to ask what protein properties have been conserved by the purifying selection. We have investigated a protein-level property, namely, the isoelectric point (pI). This property is important for enzymes because protein-protein and enzyme-substrate interactions are often electrostatic in nature. Thus, we should expect the pI of an enzyme to be similar in different organisms if protein-protein and/or enzyme-substrate interactions are conserved. In order to test this idea we gathered the protein sequences of the orthologs of ERG2, 6, 24, 25, 26 and 28. Next, we asked whether the corresponding <italic>Drosophila</italic> and <italic>Caenorhabditis</italic> sequences displayed outlier pI values. None of the orthologous proteins in the auxotrophs was detected as outlier, suggesting conservation of this physico-chemical property within the nematode and the fly lineages in spite of primary sequence divergence (##TAB##1##Table 2##). Similarity in pI is expected, when the orthologs are detected as the best reciprocal significant hits in BLAST searches. However, the existence of similarity in the hydrophobic profiles and the pI between the potential orthologs, for which the BLAST comparisons failed to be significant, is more surprising. It means that, in spite of important sequence divergence, which renders BLAST unreliable in many cases, the proteins have maintained substantial structural and physico-chemical similarity.</p> <title>Analyzing expression-profiling data to explore the potential function to the divergent <italic>ERG</italic> orthologs</title> <p>Co-expression is indicative of i) physical interaction between proteins and ii) of membership to the same complex or molecular process ##REF##11779829##[37]##. This is called the paradigm of “guilt by association” ##REF##14551426##[38]##. Thus, we have used published microarray expression data, downloadable from the Gene Expression Omnibus of the NCBI, to investigate potential co-expression patterns of the <italic>ERG</italic> orthologs that we have described above for both <italic>D. melanogaster</italic> and <italic>C. elegans</italic>. Co-expression can be assessed by determining the correlation coefficient. The correlation coefficient can be artificially inflated by flat profiles (no changes in the expression of the relevant genes). To avoid this, we focused on experiments where the genes of interest display strong variation (see ##FIG##4##Figure 5A and B##). Thus, we gathered data concerning 51 different microarray experiments for <italic>D. melanogaster</italic> respecting the criterion outlined previously. We performed a similar analysis for <italic>C. elegans</italic> but, unfortunately, the most interesting genes showed flat profiles (close to 0 in all experiments) and it was not possible to proceed further with the analysis.</p> <p>First, we asked whether the expression profiles of the <italic>D. melanogaster ERG</italic> orthologs were correlated (##TAB##2##Table 3##). The strongest correlation was found between the <italic>dLBR</italic> (<italic>ERG24</italic>) and CG1998 (<italic>ERG25</italic>) with an associated p-value of 10<sup>−16</sup> (after a Bonferroni correction). Such a p-value means that only one correlation coefficient out of 10<sup>16</sup> is expected to be as high as 0.89 just by chance (for n = 51 experimental points). Considering the maximum number of possible correlations for the 14000 transcripts in the microarrays (representing the <italic>Drosophila</italic> genome), such a high R cannot be found by chance. The behavior of dLBR and CG1998 might be reminiscent of the situation in yeast because ERG24p and ERG25p are supposed to interact, according to Epistatic MiniArray Profiling experiments ##REF##16269340##[39]##. CG7724 and CG11162 displayed the second highest R (R = 0.4, p = 0.03), but this R is not relevant in genomic terms.</p> <p>At first, we were expecting good expressional correlation among the <italic>ERG</italic> orthologs in <italic>Drosophila</italic>. However, it seems that only dLBR and CG1998 still “remember” their ancestral belonging to the sterol biosynthesis pathway. Poor expressional correlation among the rest of the <italic>ERG</italic> orthologs also suggests that the corresponding proteins either have lost their ability to physically interact in order to form stable complexes, or they do so in conditions/moments not covered by the microarray experiments explored here. Then, we focused our attention on <italic>dLBR</italic> and CG1998 by determining which other genes were expressionally correlated with them. For this, we gathered 84 genes displaying R≥0.875 with respect to both genes. For a correlation involving 51 data points, this R cut-off is associated with a safe p-value of 10<sup>−13</sup> after correction (##FIG##4##Figure 5C##). In order to get insights about these 84 genes, we used the functional classification tool of the DAVID database (<ext-link ext-link-type="uri" xlink:href="http://david.abcc.ncifcrf.gov/">http://david.abcc.ncifcrf.gov/</ext-link>
##REF##12734009##[40]##). This software provides a rapid means to organize large lists of genes into functionally related groups and to unravel biological relationships.</p> <p>In the analysis using DAVID, the most overrepresented class included genes encoding membrane proteins, often targeted to the ER, where sterol biosynthesis takes place in prototrophs (Group 1, in ##TAB##3##Table 4##). Interestingly, several of these proteins are supposed to be involved in co-translational protein targeting to membranes, signal peptide recognition, heat shock protein-binding, as well as unfolded protein binding, or to be elements of the translocon (a complex of proteins associated with the translocation of nascent polypeptides into the ER ##REF##15130556##[41]##). The following functional category (Group 2) contained four transporter proteins while the last functional group involved chaperones (i.e. peptidyl-prolyl <italic>cis-trans</italic> isomerase), chaperone cofactors or unfolded protein-binding factors. A similar analysis conducted using g∶profiler ##REF##17478515##[42]## confirmed that genes encoding protein folding actors were overrepresented among the genes displaying strong expressional correlation with dLBR and CG1998 (p&lt;10<sup>−5</sup>). The existence of expressional correlation does not imply any causality. In fact, from this exploration it is not possible to determine whether dLBR and CG1998 somehow interact with other partners to participate in intracellular protein trafficking or folding, or on the contrary, they undergo the action of the latter. Since dLBR has been shown to be a nuclear protein ##REF##7679672##[31]##, ##REF##8157663##[32]##, it would be interesting to investigate whether the dLBR transcript or protein are somehow processed to produce a C-terminal polypeptide that might be sorted to the ER. That would explain the expressional correlation between dLBR and ER proteins observed above. On the other hand, we have also explored the annotation of the genes expressionally correlated with the rest of the <italic>ERG</italic> orthologs. However, no unifying theme emerged from this analysis (data not shown, available upon request). All in all, the strong expressional correlation between dLBR and CG1998 with proteins involved in intracellular protein trafficking or folding, and the absence of such correlation with other Erg orthologs (that also require chaperons) suggest that the involvement of dLBR and CG1998 in both processes is worth exploring.</p> <p>In a previous paper, given the structural similarity between cholesterol and ecdysteroids, we had proposed that divergent ERGp orthologs might somehow participate in the synthesis of the latter ##REF##11910120##[23]##. We have therefore assessed the expressional correlation between candidate genes involved in this process: <italic>Dare1</italic>/CG12390 ##REF##10498693##[43]##, <italic>Jhamt</italic>
##UREF##9##[44]## and <italic>Start1</italic>
##UREF##10##[45]##, with <italic>dLBR</italic> and <italic>CG1998</italic>. While for <italic>Jahmt</italic> and <italic>Dare1</italic> the values of R are below 0.6, a very strong correlation (R = 0.8) was found for <italic>Start1</italic>. Interestingly, <italic>Start1</italic>, which is involved in intra-mitochondrial sterol transport, is expressed ubiquitously. However, <italic>in situ</italic> hybridization demonstrates a stronger expression in the prothoracic gland, where ecdysteroids are synthesized from cholesterol. These and other observations are consistent with the idea that Start1 plays a key role in the regulation of ecdysteroid synthesis ##UREF##10##[45]##. The potential functional link between dLBR, CG1998 and Start1 is also worth exploring.</p> <p>In conclusion, we detected a preservation of <italic>ERG</italic> genes in <italic>Drosophila melanogaster</italic> and <italic>Caenorhabditis elegans</italic>. In spite of their sequence divergence with respect to the corresponding orthologs in sterol prototrophs, they still are under selective pressure. Since insects are unable to synthesize cholesterol <italic>de novo</italic>, an appealing way to explain this evolutionary acceleration is that ERGp orthologs have other biological functions in addition to sterol synthesis. This is clearly the case of the LBR, which is also a reductase in sterol prototrophs. Shut-down of cholesterogenesis in insects and nematodes would have allowed these proteins to evolve as much as their other functions were not compromised ##REF##11910120##[23]##. Another, less parsimonious, explanation would be the evolution of different novel functions. Our microarray meta-analysis shows strong expressional correlation between the orthologs of <italic>ERG24</italic> and <italic>ERG25</italic> in <italic>D. melanogaster</italic> and genes encoding factors involved in intracellular protein trafficking and folding. This is compatible with our idea that ERGp might be involved in other biological roles in addition to sterol synthesis. The potential link between ERG proteins and intracellular protein trafficking and folding deserves experimental exploration not only in <italic>Drosophila</italic> and <italic>Caenorhabditis</italic> but also in sterol prototrophs. Moreover, the potential link between dLBR, CG1998 and Start1 is to be explored in <italic>D. melanogaster</italic>. This is compatible with our previous idea of a potential implication of these proteins in the synthesis of ecdysteroids. We hope that this genomic exploration and the hypotheses prompted here might open new avenues of experimental research.</p>
<p>Conceived and designed the experiments: GV XX RAV. Performed the experiments: RAV. Analyzed the data: GV XX RAV. Wrote the paper: GV XX RAV.</p> <title>Background</title> <p>It is known that primary sequences of enzymes involved in sterol biosynthesis are well conserved in organisms that produce sterols <italic>de novo</italic>. However, we provide evidence for a preservation of the corresponding genes in two animals unable to synthesize cholesterol (auxotrophs): <italic>Drosophila melanogaster</italic> and <italic>Caenorhabditis elegans</italic>.</p> <title>Principal Findings</title> <p>We have been able to detect <italic>bona fide</italic> orthologs of several <italic>ERG</italic> genes in both organisms using a series of complementary approaches. We have detected strong sequence divergence between the orthologs of the nematode and of the fruitfly; they are also very divergent with respect to the orthologs in organisms able to synthesize sterols <italic>de novo</italic> (prototrophs). Interestingly, the orthologs in both the nematode and the fruitfly are still under selective pressure. It is possible that these genes, which are not involved in cholesterol synthesis anymore, have been recruited to perform different new functions. We propose a more parsimonious way to explain their accelerated evolution and subsequent stabilization. The products of <italic>ERG</italic> genes in prototrophs might be involved in several biological roles, in addition to sterol synthesis. In the case of the nematode and the fruitfly, the relevant genes would have lost their ancestral function in cholesterogenesis but would have retained the other function(s), which keep them under pressure.</p> <title>Conclusions</title> <p>By exploiting microarray data we have noticed a strong expressional correlation between the orthologs of <italic>ERG24</italic> and <italic>ERG25</italic> in <italic>D. melanogaster</italic> and genes encoding factors involved in intracellular protein trafficking and folding and with <italic>Start1</italic> involved in ecdysteroid synthesis. These potential functional connections are worth being explored not only in <italic>Drosophila</italic>, but also in <italic>Caenorhabditis</italic> as well as in sterol prototrophs.</p>
<fig id="pone-0002883-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002883.g001</object-id><label>Figure 1</label><caption><title>Outline of the ergosterol synthesis pathway in yeast.</title><p>(+) the corresponding gene is present in <italic>C. elegans</italic> and <italic>D. melanogaster</italic>, according to our exploration. (−) the corresponding gene is absent. (?) not convincing evidence for the presence of the ortholog.</p></caption></fig> <fig id="pone-0002883-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002883.g002</object-id><label>Figure 2</label><caption><title>Details of the BLAST analysis that allowed the detection of <italic>ERG</italic> orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic>.</title><p>The comparisons were performed in the directions indicated by the arrows (i.e. Yeast→Human represents a BLASTp search with an ERG protein from yeast in the division <italic>Homo sapiens</italic> of Genbank). Low Scores S and high E-values (&gt;10<sup>−6</sup>) are classically considered as non-significant (unrelated or divergent sequences). The small tables display the S and E-values for the comparisons using the species-specific divisions of Genbank.</p></caption></fig> <fig id="pone-0002883-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002883.g003</object-id><label>Figure 3</label><caption><title>Segments of the <italic>Drosophila</italic> paralogs CG1998 and CG11162 (homologs of Erg25) and the corresponding conceptual translations.</title><p>The interruption of the open reading frames of both genes, by their last intron, is shown by vertical lines.</p></caption></fig> <fig id="pone-0002883-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002883.g004</object-id><label>Figure 4</label><caption><title>Hydrophobic profiles of several potential ERG orthologous proteins.</title><p>A way to show a structural relationship is to predict the hydrophobic profiles of the relevant proteins. Here we have used the TopPred program ##REF##7704669##[36]## as implemented in the server of Pasteur Institute (<ext-link ext-link-type="uri" xlink:href="http://www.pasteur.fr">http://www.pasteur.fr</ext-link>). Left panels show the results for the potential homologs of ERG2p while the right panels display the profiles for ERG28p homologs (using the Kyte-Doolitle scale, with the default parameters). Negative (positive) values represent hydrophobic (hydrophilic) segments. A way to statistically assess the similarity of two profiles is to calculate their correlation coefficient R. R-values for pairwise comparisons with the human sequence are reported. We tried to maximize the R-value by slightly sliding one profile over the other (that is why the frames of the profiles are not perfectly aligned).</p></caption></fig> <fig id="pone-0002883-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002883.g005</object-id><label>Figure 5</label><caption><title>Expression profiles of several genes expressionally correlated with dLBR and CG1998.</title><p>Panels A and B show schematic profiles displaying very strong R. However, in panel A “co-variation” is not very informative, in this case either the red and blue genes are not expressed or do not change their expression in the conditions analyzed. To avoid this artificial inflation of R, we have focused on experiments where the genes of interest display strong variation (as in panel B). Panel C: Expression profiles of several genes involved in protein trafficking and folding, which co-vary with the expression of dLBR and CG1998. The profile of Start1 is also represented.</p></caption></fig>
<table-wrap id="pone-0002883-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002883.t001</object-id><label>Table 1</label><caption><title>Ka/Ks and dN/dS values for several orthologs of <italic>ERG</italic> genes in the fruitfly and the worm.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Orthologs</td><td colspan="2" align="left" rowspan="1">
<italic>C. elegans</italic> vs <italic>C. briggsae</italic>
</td><td colspan="2" align="left" rowspan="1">
<italic>D. melanogaster</italic> vs. <italic>D. pseudobscura</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ERG</italic>
</td><td align="left" rowspan="1" colspan="1">Ka/Ks</td><td align="left" rowspan="1" colspan="1">dN/dS</td><td align="left" rowspan="1" colspan="1">Ka/Ks</td><td align="left" rowspan="1" colspan="1">dN/dS</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.127</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.080</td><td align="left" rowspan="1" colspan="1">0.029</td></tr><tr><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1">0.044</td><td align="left" rowspan="1" colspan="1">0.029</td><td align="left" rowspan="1" colspan="1">0.041</td><td align="left" rowspan="1" colspan="1">0.016</td></tr><tr><td align="left" rowspan="1" colspan="1">24</td><td align="left" rowspan="1" colspan="1">0.055</td><td align="left" rowspan="1" colspan="1">0.055</td><td align="left" rowspan="1" colspan="1">0.122</td><td align="left" rowspan="1" colspan="1">0.085</td></tr><tr><td align="left" rowspan="1" colspan="1">25A</td><td align="left" rowspan="1" colspan="1">0.041</td><td align="left" rowspan="1" colspan="1">0.039</td><td align="left" rowspan="1" colspan="1">0.064</td><td align="left" rowspan="1" colspan="1">0.041</td></tr><tr><td align="left" rowspan="1" colspan="1">25B</td><td align="left" rowspan="1" colspan="1">0.055</td><td align="left" rowspan="1" colspan="1">0.031</td><td align="left" rowspan="1" colspan="1">0.138</td><td align="left" rowspan="1" colspan="1">0.087</td></tr><tr><td align="left" rowspan="1" colspan="1">28</td><td align="left" rowspan="1" colspan="1">0.080</td><td align="left" rowspan="1" colspan="1">0.064</td><td align="left" rowspan="1" colspan="1">0.033</td><td align="left" rowspan="1" colspan="1">0.013</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002883-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002883.t002</object-id><label>Table 2</label><caption><title>Isoelectric points (pI) of the ERGp orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic>.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Proteins (n)</td><td align="left" rowspan="1" colspan="1">
<italic>C. elegans</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>D. melanogaster</italic>
</td><td align="left" rowspan="1" colspan="1">pI (w/o) mean+/−std</td><td align="left" rowspan="1" colspan="1">pI (w) mean+/−std</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<bold>Erg6 (7)</bold>
</td><td align="left" rowspan="1" colspan="1">5.95 (H14E04.1)</td><td align="left" rowspan="1" colspan="1">5.79 (CG8067)</td><td align="left" rowspan="1" colspan="1">5.98+/−0.38</td><td align="left" rowspan="1" colspan="1">6.19+/−0.83</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>Erg26 (16)</bold>
</td><td align="left" rowspan="1" colspan="1">6.52 (ZC8.1)</td><td align="left" rowspan="1" colspan="1">8.99 (CG7724)</td><td align="left" rowspan="1" colspan="1">7.93+/−1.23</td><td align="left" rowspan="1" colspan="1">7.91+/−1.22</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>Erg24 (13)</bold>
</td><td align="left" rowspan="1" colspan="1">9.38 (B0250.9)</td><td align="left" rowspan="1" colspan="1">9.89 (CG17952)</td><td align="left" rowspan="1" colspan="1">9.19+/−0.67</td><td align="left" rowspan="1" colspan="1">9.25+/−0.64</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>Erg28 (4)</bold>
</td><td align="left" rowspan="1" colspan="1">9.02 (C14C10.6)</td><td align="left" rowspan="1" colspan="1">9.62 (CG17270)</td><td align="left" rowspan="1" colspan="1">9.31+/−0.39</td><td align="left" rowspan="1" colspan="1">9.31+/−0.35</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>Erg2 (12)</bold>
</td><td align="left" rowspan="1" colspan="1">9.39 (W08F4.3)</td><td align="left" rowspan="1" colspan="1">5.49 (HDC14735)</td><td align="left" rowspan="1" colspan="1">6.28+/−1.02</td><td align="left" rowspan="1" colspan="1">6.44+/−1.28</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>Erg25 (15)</bold>
</td><td align="left" rowspan="1" colspan="1">7.57 (F49E12.9/25A)</td><td align="left" rowspan="1" colspan="1">8.17 (CG1998/25A)</td><td align="left" rowspan="1" colspan="1">7.83+/−0.63</td><td align="left" rowspan="1" colspan="1">7.96+/−0.73</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">8.4 (F49E12.10/25B)</td><td align="left" rowspan="1" colspan="1">9.76 (CG11162/25B)</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002883-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002883.t003</object-id><label>Table 3</label><caption><title>Expressional correlation among <italic>D. melanogaster ERG</italic> orthologs.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">CG1998 ERG25A</td><td align="left" rowspan="1" colspan="1">CG11162 ERG25B</td><td align="left" rowspan="1" colspan="1">LBR ERG24</td><td align="left" rowspan="1" colspan="1">CG17270 ERG28</td><td align="left" rowspan="1" colspan="1">CG8067 ERG6</td><td align="left" rowspan="1" colspan="1">CG2453 ERG6?</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<bold>CG7724</bold>
</td><td align="left" rowspan="1" colspan="1">0.22</td><td align="left" rowspan="1" colspan="1">0.40</td><td align="left" rowspan="1" colspan="1">0.18</td><td align="left" rowspan="1" colspan="1">0.14</td><td align="left" rowspan="1" colspan="1">0.32</td><td align="left" rowspan="1" colspan="1">0.27</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>CG1998</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.03</td><td align="left" rowspan="1" colspan="1">0.89</td><td align="left" rowspan="1" colspan="1">−0.08</td><td align="left" rowspan="1" colspan="1">0.01</td><td align="left" rowspan="1" colspan="1">0.14</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>CG11162</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.07</td><td align="left" rowspan="1" colspan="1">0.27</td><td align="left" rowspan="1" colspan="1">0.28</td><td align="left" rowspan="1" colspan="1">0.34</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>LBR</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">−0.06</td><td align="left" rowspan="1" colspan="1">−0.02</td><td align="left" rowspan="1" colspan="1">0.17</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>CG17270</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.31</td><td align="left" rowspan="1" colspan="1">0.41</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>CG8067</bold>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">0.55</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002883-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002883.t004</object-id><label>Table 4</label><caption><title>Functional clustering of genes whose expression profiles stronlgy correlate with those of dLBR and CG1998 (using the DAVID classification tool at Medium stringency).</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gene Group 1</td><td align="left" rowspan="1" colspan="1">Enrichment Score: 4.16</td><td align="left" rowspan="1" colspan="1">Gene Name</td><td align="left" rowspan="1" colspan="1">Key words</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG5885</bold>
</td><td align="left" rowspan="1" colspan="1">cotranslational protein targeting to membrane, integral to ER membrane, SRP, translocon complex</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>SURF4</bold>
</td><td align="left" rowspan="1" colspan="1">surfeit 4, receptor signaling protein activity, asymmetric protein localization, ER membrane</td></tr><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG8583</bold>
</td><td align="left" rowspan="1" colspan="1">HSP binding, SRP binding, unfolded protein binding, SRP receptor complex</td></tr><tr><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>SEC61ALPHA</bold>
</td><td align="left" rowspan="1" colspan="1">protein translocase activity, SRP-dependent cotranslational protein targeting to membrane, translocon complex</td></tr><tr><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG32700</bold>
</td><td align="left" rowspan="1" colspan="1">oxidoreductase activity</td></tr><tr><td align="left" rowspan="1" colspan="1">6</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG33162</bold>
</td><td align="left" rowspan="1" colspan="1">SRP receptor</td></tr><tr><td align="left" rowspan="1" colspan="1">7</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>GP210</bold>
</td><td align="left" rowspan="1" colspan="1">transporter activity, protein targeting, integral to membrane</td></tr><tr><td align="left" rowspan="1" colspan="1">8</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>SSRBETA</bold>
</td><td align="left" rowspan="1" colspan="1">signal sequence receptor</td></tr><tr><td align="left" rowspan="1" colspan="1">9</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>PROMININ-LIKE</bold>
</td><td align="left" rowspan="1" colspan="1">prominin-like protein, intracellular protein transport, integral to membrane,</td></tr><tr><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG33105</bold>
</td><td align="left" rowspan="1" colspan="1">intracellular protein transport, Golgi apparatus, integral to membrane</td></tr><tr><td align="left" rowspan="1" colspan="1">11</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG1967</bold>
</td><td align="left" rowspan="1" colspan="1">intracellular transporter activity, post-Golgi vesicle-mediated transport, coated vesicle,integral to membrane</td></tr><tr><td align="left" rowspan="1" colspan="1">12</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG11857</bold>
</td><td align="left" rowspan="1" colspan="1">vesicle-mediated transport</td></tr><tr><td align="left" rowspan="1" colspan="1">13</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>GFR</bold>
</td><td align="left" rowspan="1" colspan="1">GDP-fucose transport, Golgi-associated vesicle</td></tr><tr><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG1751</bold>
</td><td align="left" rowspan="1" colspan="1">Probable signal peptidase complex subunit 2</td></tr><tr><td align="left" rowspan="1" colspan="1">15</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG33214</bold>
</td><td align="left" rowspan="1" colspan="1">FGF binding, receptor binding, cell adhesion, intracellular protein transport</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>Gene Group 2</bold>
</td><td align="left" rowspan="1" colspan="1">Enrichment Score: 1.84</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG5594</bold>
</td><td align="left" rowspan="1" colspan="1">amino acid-polyamine transporter activity, cation transporter activity</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG15094</bold>
</td><td align="left" rowspan="1" colspan="1">high affinity inorganic phosphate∶sodium symporter activity,serine-type endopeptidase inhibitor activity</td></tr><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG8291</bold>
</td><td align="left" rowspan="1" colspan="1">neurotransmitter transporter activity</td></tr><tr><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG6293</bold>
</td><td align="left" rowspan="1" colspan="1">L-ascorbate∶sodium symporter activity, nucleotide and nucleic acid transport</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold>Gene Group 3</bold>
</td><td align="left" rowspan="1" colspan="1">Enrichment Score: 1.64</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>TORP4A</bold>
</td><td align="left" rowspan="1" colspan="1">torsin-like protein precursor, unfolded protein binding, protein folding</td></tr><tr><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG7872</bold>
</td><td align="left" rowspan="1" colspan="1">HSP binding</td></tr><tr><td align="left" rowspan="1" colspan="1">3</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>FKBP13</bold>
</td><td align="left" rowspan="1" colspan="1">peptidyl-prolyl cis-trans isomerase activity, protein folding, ER</td></tr><tr><td align="left" rowspan="1" colspan="1">4</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>CG14715</bold>
</td><td align="left" rowspan="1" colspan="1">peptidyl-prolyl cis-trans isomerase activity, protein folding, ER</td></tr><tr><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<bold>GP93</bold>
</td><td align="left" rowspan="1" colspan="1">glycoprotein 93, unfolded protein binding, response to stress, ER</td></tr></tbody></table></alternatives></table-wrap>
<table-wrap-foot><fn id="nt101"><label/><p>The mean pI values (and standard deviations –std) of the orthologous proteins in cholesterol prototrophs are shown for comparison. “pI (w/o)” stands for means +/− std calculated <bold>w</bold>ith<bold>o</bold>ut taking into account the orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic> while “pI (w)” stands for means +/− std calculated including the orthologs in <italic>C. elegans</italic> and <italic>D. melanogaster</italic>. “n” is the number of orthologs in cholesterol prototrophs used to calculate the means.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt102"><label/><p>The analyses were performed using data downloaded from the Gene expression-Omnibus database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgiCMDsearchDBgeo">http://www.ncbi.nlm.nih.gov/entrez/query.fcgiCMDsearchDBgeo</ext-link>). We considered the datasets GDS192 (wing imaginal disc spatial gene expression), GDS653 (neurotransmitter-specific neuronal gene expression), GDS664 (splicing factor mutant at permissive and restrictive temperatures) and GDS667 (mRNA splicing factor knock-down) which contain 51 data points for 14 000 transcripts.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt103"><label/><p>ER: endoplasmic reticulum, HSP: heat shock protein, SRP: signal recognition particle. At high stringency, only the first functional cluster is obtained.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>Supported by the INSERM and Paris VII university.</p></fn></fn-group>
<graphic xlink:href="pone.0002883.g001"/> <graphic xlink:href="pone.0002883.g002"/> <graphic xlink:href="pone.0002883.g003"/> <graphic xlink:href="pone.0002883.g004"/> <graphic id="pone-0002883-t001-1" xlink:href="pone.0002883.t001"/> <graphic id="pone-0002883-t002-2" xlink:href="pone.0002883.t002"/> <graphic xlink:href="pone.0002883.g005"/> <graphic id="pone-0002883-t003-3" xlink:href="pone.0002883.t003"/> <graphic id="pone-0002883-t004-4" xlink:href="pone.0002883.t004"/>
[{"label": ["5"], "element-citation": ["\n"], "surname": ["Mann", "Beachy"], "given-names": ["RK", "PA"], "year": ["2000"], "article-title": ["Cholesterol modification of proteins."], "source": ["Bioch Bioph Acta"], "volume": ["1529"], "fpage": ["188"], "lpage": ["202"]}, {"label": ["9"], "element-citation": ["\n"], "surname": ["Mo", "Bard"], "given-names": ["C", "M"], "year": ["2005a"], "article-title": ["A systematic study of yeast sterol biosynthetic protein-protein interactions using the split-ubiquitin system."], "source": ["Bioch Bioph Acta"], "volume": ["1737"], "fpage": ["152"], "lpage": ["160"]}, {"label": ["11"], "element-citation": ["\n"], "surname": ["Mo", "Valachovic", "Bard"], "given-names": ["C", "M", "M"], "year": ["2004"], "article-title": ["The ERG28-encoded protein, Erg28p, interacts with both the sterol C-4 demethylation enzyme complex as well as the late biosynthetic protein, the C-24 sterol methyltransferase (Erg6p)."], "source": ["Bioch Bioph Acta"], "volume": ["1686"], "fpage": ["30"], "lpage": ["36"]}, {"label": ["12"], "element-citation": ["\n"], "surname": ["Mo", "Milla", "Athenstaedt", "Ott", "Balliano", "Daum", "Bard"], "given-names": ["C", "P", "K", "R", "G", "G", "M"], "year": ["2003"], "article-title": ["In yeast sterol biosynthesis the 3-keto reductase protein (Erg27p) is required for oxidosqualene cyclase (Erg7p) activity."], "source": ["Bioch Bioph Acta"], "volume": ["1633"], "fpage": ["68"], "lpage": ["74"]}, {"label": ["18"], "element-citation": ["\n"], "surname": ["Schroeder", "Jefferson", "Kier", "Knittel", "Scallen"], "given-names": ["F", "JR", "AB", "J", "TJ"], "year": ["1991"], "article-title": ["Membrane cholesterol dynamics: cholesterol domains and kinetic pools."], "source": ["Proc Socr Exp Biol Med"], "volume": ["196"], "fpage": ["235"], "lpage": ["252"]}, {"label": ["19"], "element-citation": ["\n"], "surname": ["Schroeder", "Woodford", "Kavecansky", "Wood", "Joiner"], "given-names": ["F", "JK", "J", "WG", "C"], "year": ["1995"], "article-title": ["Cholesterol domains in biological membranes."], "source": ["Mol Membrane Biol"], "volume": ["12"], "fpage": ["113"], "lpage": ["119"]}, {"label": ["27"], "element-citation": ["\n"], "surname": ["Silve", "Dupuy", "Ferrara", "Loison"], "given-names": ["S", "PH", "P", "G"], "year": ["1998"], "article-title": ["Human lamin B receptor exhibits sterol C14-reductase activity in Saccharomyces cerevisiae."], "source": ["Bioch Biophys Acta"], "volume": ["1392"], "fpage": ["233"], "lpage": ["44"]}, {"label": ["29"], "element-citation": ["\n"], "surname": ["Ye", "Barton", "Worman", "Hermann", "Robin Harris"], "given-names": ["Q", "RM", "HJ", "H", "J"], "year": ["1998"], "article-title": ["Nuclear lamin binding proteins."], "source": ["Subcellular Biochemistry: Intermediate Filaments."], "publisher-loc": ["London, UK"], "publisher-name": ["Plenum Publishing"], "fpage": ["587"], "lpage": ["610"]}, {"label": ["34"], "element-citation": ["\n"], "surname": ["Breitling", "Krazeisen", "Moller", "Adamski"], "given-names": ["R", "A", "G", "J"], "year": ["2001"], "article-title": ["17beta-hydroxysteroid dehydrogenase type 7\u2013an ancient 3-ketosteroid reductase of cholesterogenesis."], "source": ["Mol Cell Endo"], "volume": ["171"], "fpage": ["199"], "lpage": ["204"]}, {"label": ["44"], "element-citation": ["\n"], "surname": ["Shinoda", "Itoyama"], "given-names": ["T", "K"], "year": ["2003"], "article-title": ["Juvenile hormone acid methyltransferase: a key regulatory enzyme for insect metamorphosis."], "source": ["Proc of the Natl Acad Sci U S A"], "volume": ["100"], "fpage": ["11986"], "lpage": ["11991"]}, {"label": ["45"], "element-citation": ["\n"], "surname": ["Roth", "Gierl", "Vollborn", "Meise", "Lintermann"], "given-names": ["GE", "MS", "L", "M", "R"], "year": ["2004"], "article-title": ["The Drosophila gene Start1: a putative cholesterol transporter and key regulator of ecdysteroid synthesis."], "source": ["Proc of the Natl Acad Sci U S A"], "volume": ["101"], "fpage": ["1601"], "lpage": ["1606"]}, {"label": ["46"], "element-citation": ["\n"], "surname": ["Xia"], "given-names": ["X"], "year": ["2001"], "source": ["Data analysis in molecular biology and evolution"], "publisher-loc": ["Boston"], "publisher-name": ["Kluwer Academic Publishers"], "fpage": ["277, 38"], "lpage": ["39"]}]
PLoS One. 2008 Aug 6; 3(8):e2883