Skip to main content
  • Research article
  • Open access
  • Published:

Evolution of the holozoan ribosome biogenesis regulon

Abstract

Background

The ribosome biogenesis (RiBi) genes encode a highly-conserved eukaryotic set of nucleolar proteins involved in rRNA transcription, assembly, processing, and export from the nucleus. While the mode of regulation of this suite of genes has been studied in the yeast, Saccharomyces cerevisiae, how this gene set is coordinately regulated in the larger and more complex metazoan genomes is not understood.

Results

Here we present genome-wide analyses indicating that a distinct mode of RiBi regulation co-evolved with the E(CG)-binding, Myc:Max bHLH heterodimer complex in a stem-holozoan, the ancestor of both Metazoa and Choanoflagellata, the protozoan group most closely related to animals. These results show that this mode of regulation, characterized by an E(CG)-bearing core-promoter, is specific to almost all of the known genes involved in ribosome biogenesis in these genomes. Interestingly, this holozoan RiBi promoter signature is absent in nematode genomes, which have not only secondarily lost Myc but are marked by invariant cell lineages typically producing small body plans of 1000 somatic cells. Furthermore, a detailed analysis of 10 fungal genomes shows that this holozoan signature in RiBi genes is not found in hemiascomycete fungi, which evolved their own unique regulatory signature for the RiBi regulon.

Conclusion

These results indicate that a Myc regulon, which is activated in proliferating cells during normal development as well as during tumor progression, has primordial roots in the evolution of an inducible growth regime in a protozoan ancestor of animals. Furthermore, by comparing divergent bHLH repertoires, we conclude that regulation by Myc but not by other bHLH genes is responsible for the evolutionary maintenance of E(CG) sites across the RiBi suite of genes.

Background

Ribosome biogenesis (RiBi) is a primary function of the nucleolus [1–3]. In the nucleolus, rRNA molecules are synthesized as precursors by DNA-directed RNA Pol I and Pol III. Nascent rRNAs then undergo extensive chemical modifications and RNA cleavage reactions. Numerous RiBi proteins are involved both in this enzymatic processing as well in assisting with proper rRNA folding to produce functional ribosomal subunits [1]. Synthesizing functional ribosomes requires immense coordination because gene products from all three DNA-dependent RNA polymerases are required to ensure proper stoichiometry of ribosomal components. This level of co-regulation is likely to be controlled through a highly-specific DNA signature as seen in other gene regulatory systems [4]. Such a signature would allow the appropriate factors to co-ordinately regulate the RiBi genes as a distinct regulon. For example, in the yeast Saccharomyces cerevisiae two important regulatory motifs consisting of the PAC (polymerase A and C) and RRPE (Ribosomal RNA Processing Element) motifs have been identified in RiBi genes [5–7]. In animals no known factor or motif is known to coordinate the entire RiBi set. However, the metazoan transcription factor Myc has been found to target at least some RiBi genes.

Myc is a bHLH DNA-binding protein present in animals where it plays an important role in cell growth and proliferation by regulating gene expression during development and tumorigenesis [8–10]. Myc functions by heterodimerizing with its obligate partner Max to bind the sequence 5'-CACGTG [a CG-core E-box, or E(CG)] and transactivate target genes [11–14]. Previous studies have implicated the Myc transcription factor in rRNA transcription [15–17], but the possibility that Myc could directly regulate the hundreds of gene products involved in RiBi, as opposed to a few early genes or key steps [18–20], has not been investigated. Also unknown is when such an RiBi circuit or other non-RiBi targets of Myc may have evolved and whether it was in the most recent common ancestor of Bilateria (protostomes and deuterostomes), Metazoa (animals), Holozoa (Metazoa + choanoflagellates), Opisthokonta (Holozoa + fungi), or Eukaryota.

Here, we use a multi-genomic approach to show that the vast majority of genes implicated in ribosome biogenesis are associated with E(CG)-bearing core promoters in all holozoan genomes containing Myc, and thus constitutes a uniquely holozoan RiBi regulon. Max, Mad, and Mnt, all members of the Myc bHLH superfamily, are all either insufficient or dispensable in explaining the correlation of E(CG) with RiBi as revealed by a comparison of multiple eukaryotic genomes, which differ in their bHLH repertoires. Thus, in addition to known RiBi targets of Myc [18–20] and the similar growth defects of both RiBi and myc mutant alleles [8–10, 18, 20, 21], our comparative genomic results suggest that the characteristic RiBi E(CG) core promoter architecture co-evolved with a proto-Myc:Max complex in a unicellular holozoan ancestor. This is consistent with the metabolic evolution of a unique, Myc-Max regulated, RiBi growth signalling pathway in an ancient unicellular heterotroph.

Results and discussion

Overall approach

To gain insight into the regulatory control of Ribosome Biogenesis (RiBi) in animals, we first wanted to characterize a co-regulated RiBi gene set definable by a shared regulatory signature, a component of which would be the Myc-Max binding site E(CG). A priori, we had no reason to expect whether a common regulatory signature would be found across the large number of known ribosome biogenesis genes or whether such a signature would be limited to only a subset (i.e. the sensitivity of the postulated signature for RiBi genes). Furthermore, we set no expectation on whether this signature would necessarily be present or absent in other non-RiBi genes devoted to cell growth or other functions (i.e. the specificity of the postulated signature for RiBi genes). For this purpose, we chose the Drosophila melanogaster genome because of its relatively compact size and absence of genome wide duplications characteristic of vertebrates. Both of these properties facilitate the identification of a sequence signature and the characterization of its sensitivity and specificity for entire biological functions because they simplify the ability to conduct and interpret computational queries of genome sequence.

Having identified the list of genes in this fly regulon we would then utilize the many available eukaryotic genome sequences to identify shared sequence signatures associated with this regulon in each genome. Because of the highly conserved nature of the protein functions associated with ribosome biogenesis, it is likely that this list of genes would remain co-regulated in other genomes despite evolution of the regulatory signatures associated with this regulon.

E(CG) is highly specific to the fly RiBi regulon

An overwhelming majority of confirmed Myc targets contain E(CG) near the transcriptional start site (TSS) [19]. We therefore examined all fly genes with promoter proximal E(CG) sites to determine whether they were related by a common cellular function. We searched the Drosophila melanogaster genome and identified only 390 promoter sequences containing an E(CG) site in the core promoter region (160 bp centered around +1) out of 20,468 annotated transcripts (14,752 genes). Remarkably, these genes include most proteins known to be involved throughout ribosome biogenesis [see Additional file 1]. To examine the relationship between the E(CG) motif and ribosome biogenesis in more detail, we used Gene Ontology (GO) classification of the yeast genome to map 121 fly orthologs of yeast nucleolar genes. Of these, ~75% possess E(CG), indicating that the majority are under control of a common regulatory motif in Drosophila [see Additional file 1]. As explained below, many additional fly E(CG)-bearing genes, which are not conserved as orthologs in yeast, are likely to be related to ribosome biogenesis. Furthermore, as detailed by multiple statistical tests conducted in this study, the rate of E(CG) motifs in the Drosophila RiBi gene promoters is significantly elevated relative to promoters of genes not known to be involved in nucleolar functions and/or ribosome biogenesis.

E(CG)-bearing, RiBi-type promoters are unique

Any examination restricted to gene orthologs precludes the identification of genes not conforming to a 1-to-1 orthology. To this end, it would be ideal to search the entire genome to provide a comprehensive analysis of the entire E(CG) fly RiBi regulon. This method could identify novel fly RiBi genes harboring E(CG) sites independent of a 1-to-1 orthology and ambiguous annotations of transcriptional initiation sites. To carry out this whole genome query, we first searched for additional motifs that comprised the full core promoter context of RiBi genes. With additional motifs in hand, highly specific whole genome queries for promoter-linked E(CG) sites could be achieved, thereby identifying a more complete RiBi regulon.

Myc's known co-localization to core promoters [16] suggests it may define or associate with a distinct core promoter architecture. We therefore looked for novel elements that may be specific to RiBi promoters but infrequent across other promoters. We also investigated whether promoters bearing E(CG) are associated with specific core promoter elements such as TATA-boxes, Initiator sequences, downstream promoter elements (DPEs), and other common motifs [22]. As detailed below, these results support the idea that the fly RiBi-type promoter, uniquely characterized by the E(CG) site, defines a distinctive and highly specific promoter architecture, which is useful in identifying this gene set in Drosophila.

By searching for novel motifs, we first noted that the flanks of the E(CG) site often match an extended consensus MAACACGTGYG (M = A/C, Y = C/T). Three out of every four core promoters that contain this extended E(CG) consensus map to RiBi genes (Fig. 1A). As a negative control for the specificity of the flanking sequence, we also searched the entire genome with an E(CG) motif in which the flanking pattern was maximally divergent from the observed, RiBi-specific, CG-core E-box flanking pattern. This "anti-flank" E(CG) motif, KKYCACGTGRMK (K = G/T, R = A/G), maps to almost 3-fold more sites than the extended E(CG) consensus, but nonetheless is absent from the core promoters of known nucleolar or RiBi orthologs (data not shown).

Figure 1
figure 1

E(CG) is a core promoter element of Drosophila ribosome biogenesis (RiBi) genes. (A) Fly RiBi genes (5 examples shown) generally possess three common features around the transcriptional start site (rightward pointing arrow) and upstream of the translational start (ATG). This distinct promoter architecture is characterized by a CG-core E-box (blue box), a specific E(CG) flanking motif (green box) and a coordinating cluster of sites matching the DNA Replication Element, DRE, (red boxes) spanning a distance less than 100 bp. The distance of E(CG) to the TSS for each gene is indicated above E(CG). G6375 corresponds to the pit gene, which is a known Myc target and an RiBi gene [18]. (B) This core promoter architecture identifies several functional groups of genes associated with RiBi (green circles). The number of genes is indicated for functional groups with more than 3 members. The sum of 151 genes (large circle) is the sum of all of the individual subfunctions with specific roles in Ribosome Biogenesis. The RiBi genes encode a variety of domains and protein folds including RNA-binding regions (RNP-1), C-terminal helicases, DEAD/DEAH box helicases, WD-40 repeats, ARM repeats, Histone-folds, AAA ATPases and many others [see Additional file 1]. (C) The results of genome queries in Drosophila for E(CG) type core promoters results in a highly significant enrichment of GO terms directly related to ribosome biogenesis (nucleolar, rRNA metabolism, rRNA binding, snoRNA complex, pseudouridine synthesis, ribosomal subunits, etc.)

We also identified sequences corresponding to the DRE core promoter element, 5'-CTATCGATA, as previously reported [23]. DRE (DNA-replication element) binds DREF (DRE factor) and is associated with promoters involved in DNA replication [24]. We observe a cluster of DRE sites in many RiBi gene promoters, where it occurs immediately upstream of E(CG) (Fig. 1A). When the fly genome is queried for regions containing at least 2 DRE motifs and an E(CG) within an 80 bp window, this highly specific signature identifies 126 loci in the genome. Over half of these are known to be associated with RiBi [see Additional file 1]. Tellingly, DRE motifs are not associated with E(CG) sites containing the anti-flank sequence (i.e. KKYCACGTGRMK).

This analysis provides enough information to distinguish most core-promoter linked E(CG) sites from the majority of ~15,000 E(CG) sites that lie outside of promoter regions and likely represent random background occurrences. We queried the entire Drosophila genome for the different E(CG)-type core promoter signatures [E(CG) + flanks, or E(CG) + DRE, or E(CG) + mapped 5'-end]. The largest functional group of genes identified across the genome is the RiBi gene set (151/321 genes; Fig. 1). These genes are involved at all steps of ribosomal processing including factors involved in rRNA transcription (11 genes), snoRNA processing (3 genes), snoRNPs (12 genes), 90S particles (30 genes), pre-60S particles (48 genes), 40S particles (15 genes), ribosome structure (19 genes), unknown steps in ribosome biogenesis (10 genes), and unknown functions of the nucleolus (3 genes). Some of these RiBi genes have previously been implicated as potential direct or indirect Myc targets [16, 18, 19, 23]. Some genes encode products that are only known to be localized to the nucleolus (e.g. spermidine synthase) [16]. Others participate in tRNA modification in addition to rRNA modification [25, 26]. Altogether these results support the hypothesis that the many RiBi genes are co-regulated and identifiable through a specific E(CG)-bearing promoter architecture.

Characterization of the RiBi regulon across Eukaryota

To determine the evolutionary origins of the RiBi regulon, we first measured the conservation of E(CG) motifs in the RiBi genes of other bilaterian genomes, including humans and nematodes. Of the human RiBi genes, 77%, possess an E(CG) in a window ± 600 bp from the 5' annotated end, which is a significantly elevated compared to the ~20% background level in control promoters (Fig. 2; Table 1). In contrast with humans, there is no elevated level of E(CG) in RiBi genes in the nematode genome of C. elegans despite the conservation of such genes (Fig. 2). The presence of an E(CG)-RiBi signature in both a deuterostome (humans) and a protostome (flies) suggests that the absence in another protostome (nematodes) is a secondary loss.

Table 1 Highly conserved genes with an E(CG)-bearing promoter.
Figure 2
figure 2

Holozoan RiBi promoters are enriched in E(CG) sites. The percentage of E(CG)-bearing promoters in RiBi genes is two to four fold higher in D. melanogaster (Dm), H. sapiens (Hs), N. vectensis (Nv), and M. brevicollis (Mb) relative to negative control sequences composed of promoter regions of downstream conserved genes (C1), the 3' regions of RiBi genes (C2), or the promoters of genes with GO mitochondrial classification (C M ). This difference between RiBi and C1, C2, or C M is lacking in outgroup genomes such as S. cerevisiae (Sc), which lack Myc, as well as in the nematode genome of C. elegans (Ce), which has secondarily lost Myc (Fig. 3). Inset depicts phylogenetic relationships among these organisms.

To test whether components of the fly RiBi regulon are conserved outside of Bilateria, we analyzed the presence of the E(CG) signature in the RiBi orthologous cohort in the genomes of more distantly-related organisms (Figs. 2 and 3). For a metazoan outgroup to Bilateria, we used the cnidarian genome of Nematostella vectensis [27]. For a holozoan outgroup to Metazoa, we used the choanoflagellate genome of Monosiga brevicollis [28]. For an opisthokont out-group to Holozoa, we used the baker's yeast genome of Saccharomyces cerevisiae.

Figure 3
figure 3

The RiBi-E(CG) regulon occurs only in Myc-bearing holozoan genomes. (A) Specific amino acid residues in holozoan MAX (red), MAD (blue), MNT (pink), and MYC (green) allow identification among the MYC/MAX superfamily bHLH genes (common superfamily residues in yellow and underlined). Only three bHLH genes were found in the choanoflagellate genome of Monosiga brevicollis: Mb-MYC, MbMAX and MbMUSH, corresponding to Myc and Max orthologs, and a distant MITF/USF/SREBP homolog (not shown). No Myc and Max orthologs were found outside of Holozoa. The predicted amino acid sequences of the bHLH regions of the M. brevicollis Myc/Max family of genes are shown aligned to Drosophila, Caenorhabditis, and Nematostella orthologs. (B) The presence of Myc (green filled boxes) is correlated with multiple genomes possessing the E(CG)-RiBi signature (green filled boxes). Other bHLH genes in the Myc superfamily (Max, Mad, Mnt; gray filled boxes) are either not necessary (Mad or Mnt) or insufficient (Max) to explain the occurrence of E(CG) sites in the RiBi regulon. "X" boxes indicate absence of a gene or E(CG) signature as indicated.

We identified all of the RiBi orthologs between the fly, human, and yeast genomes and their apparent orthologs in the Nematostella and Monosiga genomes and were able to find at least 100 such genes in each genome. We found that each holozoan genome (except C. elegans) possesses E(CG) sites across ~50% to 90% of its identifiable RiBi genes in the region ± 600 bp from the 5'-most end (Fig. 2). This corresponds to a statistically significant (p < 0.001) two-fold to four-fold elevated level of E(CG) relative to the core promoter regions of adjacent control genes (C1 in Fig. 2). As additional negative controls, we analyzed the frequency of E(CG) in the region ± 600 bp from the 3' end of the same test genes (C2), or in the promoters of genes with mitochondrial GO terms (C M ). We again found no enrichment over background (Fig. 2).

Importantly, we also failed to find elevated levels of E(CG) in yeast RiBi versus control groups (C1 and C M in Fig. 2). Of relevance, the yeast RiBi genes are known to be regulated by motifs that are distinct from E(CG) [29–31]. We were also unable to find the E(CG)-RiBi signature in other fungal and more distantly related eukaryotic genomes (see Figs. 4, 5 and Methods).

Figure 4
figure 4

Frequency of E(CG) in opisthokont RiBi promoters. The frequency of E(CG) in 25 opisthokont RiBi orthologs was investigated. These 25 RiBi orthologs were selected based on the presence of conserved E(CG) sites in the human, fly, sea anemone, and choanoflagellate orthologs (Table 1B). DNA sequences 500 bp upstream from the translational start sites in the RiBi orthologs of S. cerevisiae (Sc), C. glabrata (Cg), K. lactis (Kl), A. gossypii (Ag), P. stipitis (Ps), D. hansenii (Dh), C. albicans (Ca), Y. lipolytica (Yl), N. crassa (Nc), S. pombe (Sp), and M. brevicollis (Mb) were collected. For D. melanogaster (Dm), 500 bp of DNA sequence (± 250 bp) from the 5' annotated end was collected. The sequences of each of these opisthokont promoters was analyzed for the presence of E(CG) motifs. The percentage of (ECG) in each species' RiBi orthologs is depicted on the Y-axis with the number of orthologs containing E(CG) over the total ortholog number of orthologs displayed above each genome. Key nodes for latest common ancestors (LCAs) are depicted in the phylogenetic tree [59–61].

Figure 5
figure 5

Evolution of shared motifs in opisthokont RiBi core promoters. The promoter gene sets for each species depicted in Figure 4 weres analyzed by MEME to identify common cis-regulatory motifs for each lineage [57, 58]. Motifs found in greater than two-thirds of RiBi genes are depicted from left to right (highest scoring motifs to left). The frequency of each motif is expressed as a percentage in the upper left corner. The known fungal motifs RRPE [5, 6] and PAC [5–7] are shown in shown in black and blue, respectfully. The holozoan E(CG) motif identified in this work is shown in orange.

Interestingly, there are a few fly E(CG)-bearing genes whose orthologs are E(CG)-bearing across humans, cnidarians, and choanoflagellates, but are not yet known to be related to ribosome biogenesis (Table 1). One of these is the vasa (DDX4) locus, which is an RNA-binding protein and metazoan germline determinant. This gene is maintained as an E(CG)-bearing promoter in humans, ascidians, flies, cnidarians, and choanoflagellates, but not in the nematode C. elegans and could conceivably be another DDX/DHX-containing RiBi gene, which was secondarily co-opted as a metazoan germline determinant. Another corresponds to the Dph5 gene, which is involved in the diphthamide modification of a histidine residue on elongation factor 2 (EF2) [32, 33]. Dph5 is possibly required for frame-shift suppression during translation [34], but could conceivably play a role in unknown ribosomal modifications.

In contrast to the highly conserved nature of the RiBi regulon across Holozoa, there is little evolutionary conservation of E(CG) sites for genes not known to be involved in RiBi [see Additional file 1]. This indicates that these fly E(CG)-bearing genes, including some known Myc targets in flies and/or humans [e.g. CAD [35] and TIMM10 [36]], acquired the E(CG) site either in a stem bilaterian or in subsequent independent occurrences. These sites could also represent background noise due to sequence drift in a particular lineage. Many of these genes have fly E(CG) sites that are ± 600 bp of the initiation site, but are not as tightly linked to the core promoter (± 80 bp) as RiBi genes.

The E(CG)-bearing RiBi regulon is found only in Myc-bearing genomes

Like the bilaterian genomes of flies and humans, the cnidarian Nematostella, a non-bilaterian metazoan, possesses Myc and Max homologs [27, 37, 38]. Nematostella Myc homologs can also bind to E(CG) in vitro as well as rescue the proliferative defects of myc null mammalian cells (Janice Ascano, S.J.B, and M.D.C.; in preparation). Intriguingly, the non-metazoan choanoflagellate genome of Monosiga also possesses clear orthologs of both Myc and Max (Fig. 3). However, outside of Holozoa, Myc and Max genes appear to be absent. Yeast do not possess Myc or Max orthologs even though E(CG)-binding bHLH dimers are present [39]. Moreover, the yeast RiBi genes, which lack E(CG) sites, are known to be regulated by non-bHLH factors [29–31, 39]. We were also unable to find members of the Myc clade of genes in other fungal and more distantly related eukaryotic genomes (see Methods).

Secondary loss of Myc is linked to secondary loss E(CG) RiBi signature in nematodes

Unlike flies and humans, the nematode represents a bilaterian that appears to have secondarily lost Myc but not Max (Fig. 3) [38]. This loss of a nematode Myc together with the loss of the E(CG) signature in the RiBi regulon is intriguing and suggests a few hypotheses and predictions. We explore these here because the loss of both an important metazoan transcription factor and a statistically significant association of E(CG) motifs with the holozoan RiBi gene battery is noteworthy and instructive of Myc function.

One major developmental hypothesis of nematode loss of both Myc and E(CG) sites in RiBi core promoters might stem from the relatively small number of 1000 somatic cells in the adult nematode [40]. This extremely low quantity of adult cells and limited cell-proliferation may render Myc's induction of the RiBi regulon unnecessary. Thus, under this hypothesis, the RiBi regulon is under two modes of regulation in Holozoa. First, there is a basal rate of low-level expression, and second there is a Myc-dependent induced rate of elevated expression that is associated with proliferating cells. Thus, nematodes with their particular small body sizes, and limited cell numbers, might have evolved to forgo Myc-RiBi induction. If this hypothesis is correct, we might expect to find further pseudogenization of the Myc locus, in animals with similar character. Other phyla known to include animals of such type include gastrotrichs, kinorhynchs, loriciferans, nematomorphs, and some priapulids.

An alternative molecular hypothesis of nematode-specific loss of both Myc and E(CG) sites in RiBi core promoters involves the nematode phenomenon of trans-splicing [41, 42]. Nematodes employ a distinctive 5'-mRNA capping process that depends on trans-splicing of pre-capped 5' leader sequences expressed at independent loci. Because Myc has recently been found to upregulate mRNA capping of target genes [43, 44], trans-splicing might have alleviated the need for keeping the Myc gene.

Yet a third hypothesis for loss of both Myc and E(CG) sites in RiBi core promoters involves novel shared motifs in C. elegans RiBi genes. Interestingly, as described further below our investigation analysis across opisthokont genomes, reveals several additional motifs resembling the RRPE motif that is present across Opisthokonts (Fig. 5). This could indicate a greater reliance on factors binding the RRPE motif or an augmented basal level of expression. Several studies have documented the evolution of cis-regulatory signatures via substitution of transcription factors that regulate even highly conserved gene sets [45, 46].

All three hypotheses for loss of nematode myc along with loss of the E(CG)-RiBi signature are mutually non-exclusive. For instance, loss of proliferative cells and the underlying Myc genetic circuitry might have been permissive for the development of trans-splicing. Alternatively, the evolution of nematode trans-splicing might have occured first and been permissive for loss of Myc regulation.

Mad and Mnt are dispensable for the E(CG)-RiBi core promoter signature

The loss of specific bHLH genes, such as the loss of Myc in nematodes, results in different bHLH repertoires present across holozoan genomes. We therefore next examined genomic bHLH repertoires among different organisms in order to consider potential trans-factors other than Myc that might correlate with E(CG) core elements in RiBi promoters (Fig. 3).

Among the Myc superfamily, Max:Max, Mad:Max, Mnt:Max, and Myc:Max dimers, can all bind to E(CG) [12, 13, 23]. The continued presence of both Mad and Max in nematodes [38] suggests that Mad:Max or Max:Max complexes do not target the RiBi regulon via the E(CG) target site, which is absent in nematode genomes (Fig. 3A). Furthermore, a Mad ortholog is not present in Drosophila, whereas both Mnt and Mad are apparently absent in Monosiga suggesting that both are dispensable for the function of the E(CG)-RiBi signature (Fig. 3A). [Both Mad and Mnt are present in other eumetazoan genomes such as cnidiarians (N. vectensis) and sea urchins (S. purpuratus), but are reciprocally lost in flies (no Mad) and nematodes (no Mnt) (Fig. 3A).]

All of these genomic configurations suggest that E(CG) promoter signatures are direct targets of Myc:Max complexes and are consistent with both the known biochemically confirmed RiBi targets of Myc as well as the growth-related phenotypes of myc mutant alleles. Nonetheless, these results do not exclude the possibility that other E(CG)-binding proteins potentially modulate the RiBi regulon when they coexist with Myc in the genome [47]. Thus, in conjunction with our RiBi-E(CG) genomic analyses, a consideration of the bHLH repertoire of holozoan genomes supports the hypothesis that the need for regulation by Myc, but not by other bHLH genes, is responsible for the evolutionary maintenance of E(CG) sites across the RiBi suite of genes.

Evolution of RiBi core promoter sequences across Opisthokonts

While we were not able to find either Myc or Max orthologs in any available genome outside of Holozoa, it is still possible that the E(CG) motif is still associated with RiBi genes and controlled by other factors. For instance, Saccharomyces cerevisiae possesses two E(CG)-binding bHLH complexes from its small repertoire of bHLH genes [39]. These correspond to Pho2p/Pho4p heterodimers and Cbf1p homodimers. The Pho2p/Pho4p complex is involved in the regulation of phosphate biogenesis in responses to phosphate starvation, whereas Cbf in involved in centromeric function, methionine biosynthesis, sulfur metabolism, and regulating ribosomal structural proteins, [39, 46]. We therefore looked for the presence of E(CG) in the core promoters of RiBi genes of 10 different fungal genomes (Fig. 4). We specifically looked at the 25 RiBi orthologs that are E(CG)-bearing across Holozoa (Table 1B). We find that the average rate of E(CG) sites across all identifiable RiBi orthologs in the 500 bp window immediately upstream of the start ATG (this window size corresponds to the typical intergenic distances in these fungal genomes) is 16%. Three separate fungal genomes have 0% E(CG)-bearing RiBi core promoters (0/25 RiBi promoters), while Ashbya gossipyi has the highest at 37.5% (9/24 RiBi promoters). By comparison, in the 500 bp core promoter window of the 25 RiBi genes in choanoflagellates the rate is 72% (18/25), while in flies it is 84% (21/25); the few missing promoters have E(CG) motifs in the adjacent upstream 100 bp. Thus, this core set of RiBi genes is not likely to be regulated by E(CG) motifs in fungi as the majority of their promoters lack this site.

We next conducted a MEME analysis of 500 bp core promoter fragments across multiple holozoan and fungal core promoter sequences for RiBi genes to identify all potential motifs that might serve as common binding sites in each system. We find that the previously identified PAC (Polymerase A and C) RiBi motif [5–7] can be readily identified in almost all Hemiascomycete fungi except the distantly related Yarrowia lipolytica. This motif, when present, is usually found in 100% of all RiBi genes analyzed in each species. Thus, the same genes that are likely co-regulated by Myc:Max via E(CG) in Holozoa are instead co-regulated by PAC binding factors in the Hemiascomycota sub-phylum. Interestingly, among the Myc targets identified in flies, we have found subunits of DNA-dependent RNA polymerases [see Additional file 1].

Unexpectedly, the RRPE motif, which was identified as a co-occuring motif in yeast RiBi genes along with PAC motifs, appears in both fungal and holozoan RiBi genes promoters. A yeast-specific factor, Stb3, has been proposed to promote cell growth by binding to at least some RRPEs in target genes in a glucose-dependent manner [48]. However, not all RRPE-bearing yeast promoters appear to require Stb3 for induction [48]. Nonetheless, this motif may represent a more ancient and possibly conserved RiBi binding motif than previously appreciated. In this case, Holozoan-specific and Hemiascomycete-specific modes of RiBi regulation might be relatively more recent additions since their divergence. Altogether, these analyses support the conclusion that E(CG) is uniquely associated with Myc-bearing Holozoan genomes and that different signatures and factors control the same regulon in distant taxa.

Conclusion

We successfully identified a specific core-promoter signature, partially composed of the Myc:Max binding site E(CG), which is highly specific to the entire suite of genes devoted to ribosome biogenesis in Holozoa but not in fungal genomes. Based on these whole-genome analyses, and the confirmation of individual RiBi genes as Myc targets, we conclude that the entire RiBi gene set constitutes a bona fide Myc-targeted regulon. By analyzing a wide diversity of eukaryotic genomes, we show that this specific core-promoter signature is present only in holozoan genomes that still contain Myc. Furthermore, gene loss in other Myc:Max superfamily members, such as Mad or Mnt, while retaining Myc, is apparently not sufficient for loss of the E(CG) signature in the RiBi gene set. Thus, nematode genomes, which lack Myc, but not Max or Mad, do not possess the E(CG) signature across the well-conserved RiBi gene set.

A toolkit of animal-specific genes, including the bHLH family of DNA-binding factors, is thought to have been assembled early in pre-animal evolution. The bHLH family has been intriguing because of its diversification into cell-type specific functions modulating proliferation, differentiation, and metabolic programs across eukaryotes [37–39, 49, 50]. Consequently, the evolution of animals is likely to have involved the establishment of a canonical set of bHLH transcription factors regulating these downstream genomic programs. Here, we describe a large RiBi regulon that co-evolved with Myc:Max in a stem-holozoan. Significantly, Myc and Max mark the beginning of an animal-like bHLH repertoire in a pre-metazoan ancestor.

Choanoflagellates represent the sister-group to animals. Their ability to form colonies is indicative of a possible precondition in the evolution of multicellularity in metazoa [51]. Furthermore, at least one Receptor Tyrosine Kinase (RTK) has been identified in choanoflagellates [52]. RTKs were once thought to be exclusive to animals, which use them in cell-cell communication and as signaling inputs into growth factor-mediated pathways of gene activation such as Myc [53]. These results suggest a model for the origins of a Myc-induced RiBi regulon, which is commonly misregulated in diverse human cancers. Around 750 to 1000 million years ago [54], a protozoan, heterotrophic ancestor of Holozoa either adapted, or co-opted through duplication and divergence, a proto-Myc:Max bHLH heterodimer complex and evolved the capability to induce the primordial Myc moiety in response to RTK-mediated growth signals. Core promoter-bound Myc:Max complexes would then co-ordinately up-regulate the ribosome biogenesis regulon and thereby commit to cell growth and/or proliferation.

An alternative hypothesis would be that the RiBi genes have always been inducible in diverse taxa, but that this mode of regulation has diverged and/or been supplanted by distinct mechanisms in different taxa. Thus, the Holozoan E(CG) RiBi signature would not represent a new ability to up-regulate this gene set but rather a unique mechanism for inducing the RiBi gene cohort. Similarly, in the hemiascomycete fungi (except for the distantly related Yarrowia) the RiBi gene set is co-regulated by non-bHLH factors via the PAC motif.

A recent study has also identified E(CG) as present in a subset of core promoters of yeast ribosomal structural protein-encoding genes driven by Cbf1p [46]. Yeast are well known to have at least two E(CG)-binding bHLH systems in Pho2/Pho4p and Cbf1p complexes, which are involved in phosphate biogenesis and methionine biosynthesis/translation, respectively. We also see a slight elevated level of E(CG) in some fungal genomes such as Pichia and Ashbya. Thus, the Holozoan Myc and Max system might have evolved out of an Opisthokont ancestor in which E(CG) motifs might have been loosely tied to generic growth programs controlling both RiBi and RP gene sets. Under this scenario, a general cell growth pathway culminating in a bHLH induction of an undefined regulon might have existed in an ancient opisthokont ancestor. This pathway then diverged separately in fungi and Holozoa. In Holozoa, this bHLH gene evolved into a bHLH-ZIP encoding gene and subsequently duplicated and diverged to produce Max, and a growth-inducible activating form with Myc. The Myc-Max then specialized in induction of RiBi genes. In fungi, this system perhaps retained its ancestral homodimeric form in Cbf1p and controlled a more generic cell growth pathway that included ribosomal proteins, a few RiBi genes, as well as other cell growth functions. Furthermore, the RiBi genes in hemiascomycetes evolved to be largely controlled specifically by PAC-binding factors. Future research will have to be conducted in both fungal and animal genomes to explore these ideas.

Methods

Statistical analysis

Statistical tests were conducted to test for the significance of the difference between RiBi, mitochondrial (C M ), and other control (C1 or C2) gene sets. A Mann-Whitney test for the significance of the difference between RiBi genes and control orthologs (C1 or C2, see Fig. 2) was statistically significant (p < 0.001) for Hs, Dm, Nv, and Mb data. There was no statistical significance between RiBi orthologs and control data for Ce or Sc. A Mann-Whitney test for the significance of the difference between non-RiBi genes with RiBi-type promoters, and mitochondrial or control orthologs (C1 or C2, see Fig. 2) was statistically significant (p < 0.001) for Hs and Dm data. However, there was no statistical significance between RiBi-type promoter bearing non-RiBi gene orthologs and control data or mitochondrial orthologs for Nv, Mb, Ce, or Sc data. Statistical validation of over-represented GO terms shared by genes that we identified in Fig. 1 was carried out as previously described [55]. The GO term "Ribosome Biogenesis and Assembly" was found to be statistical significant (p < 1.2e - 27) based on Fisher's Exact Test (one-sided P-value) of the association between attribute and query.

Human and fly orthology

Lists of fly and human genes together with their corresponding DNA sequences (± 600 bp from 5'-annotated end) meeting orthology gene tree tests between genomes were retrieved from Ensembl data mining tool BioMart. The annotation data corresponded to Ensembl Release 44, 2007 using genome builds Drosophila melanogaster BDGP 4.3 and Homo sapiens NCBI 36. This resulted in 3066 human/fly pairs although some groups of pairs correspond to multiple isoforms in one genome. DNA sequences were searched for CACGTG E-boxes [E(CG)] using the UNIX grep and perl tools. Orthologous gene pairs with E(CG) in each genome were compared with a list of verified ribosome biogenesis (RiBi) genes [56], the list of yeast orthologs with "nucleolar" GO classification, and current results in the literature (Pubmed).

Genome assemblies and orthology identification

We used the following genome builds: Drosophila melanogaster BDGP 4.3, Homo sapiens NCBI 36, Saccharomyces cerevisiae SGD1.01, and Caenorhabditis elegans WS170. We used Ensembl Release 44 (2007) for orthology calls between these four genomes. ENSEMBL orthology calls use best reciprocal hits between genomes to cluster proteins followed by construction of maximum likelihood phylogenetic gene trees (NJTREE) and distinguish orthologs from paralogs. For identifying orthologous loci in the Nematostella vectensis and Monosiga brevicollis we used BLASTP to identify the best matches in the respective genomes to the Drosophila amino acid sequence. EST coverage in each genome allowed independent confirmation of the majority of homologous sequences. We identified all hits with E <e - 10. We then performed a reciprocal BLASTP to weed out 1-to-many hits.

Promoter sequence analysis

In searching for core promoter linked E(CG) sites we looked used the BDGP assembly release 4, flybase annotation rel.4.3-20060130. To identify potential E(CG) flanking patterns, the sequences adjacent to E(CG) sites in the conserved human/fly RiBi genes as well as other known confirmed Myc targets were aligned to identify the reported information content in the immediate flanking sequences. Additionally, these sequences were compared to a variety of control sequences composed of anti-signature 2 motifs, promoter sequences of adjacent genes, or unrelated developmental genes (anterior/posterior Drosophila developmental loci) to identify over-represented motifs spanning 6 to 8 bp with 0, 1 or 2 wild cards. This resulted in the identification of three classes of motifs: 1) E(CG) with flanking sequence, 2) DRE sequences, 3) A-rich sequences (Fig. 1). Genome queries were conducted by direct searches of the most recent Drosophila genome (BDGP 4.3) using UNIX grep and perl. Genomic queries for signatures were defined as follows. Signature 1 identifies 15,434 in the Drosophila BDGP 4.3 genome searches. Signature 2 was any sequence matching one of the following CAACACGTGCG, AAACACGTGTG, and AAACACGTGCG. Signature 3 was defined as a window of 80 bp containing CACGTG and 2 of the following sequences: CTATCG or TATCGA. Loci that mapped within 600 bp, 240 bp or 1 kb of these three signatures, respectively, were identified.

bHLH phylogenetic analyses

Over 150 bHLH amino acid sequences from plant (Arabidopsis thaliana), ciliate (Tetrahymena thermophilia), yeast (Saccharomyces cerevisiae), choanoflagellate (Monosiga brevicollis), sponge (Amphimedon queenslandica), cnidarian (Nematostella vectensis), protostomes (Drosophila melanogaster and Caenorhabditis elegans) and deuterostome (Strongylocentrotus purpuratus) organisms were aligned using CLUSTALW and used to make a primary alignment and phylogenetic guide tree. Alignments were adjusted by hand to reduce the number of insertions, and were subsequently used to generate random samples using PHYLIP Seqboot. Phylogenetic trees were generated using neighbor-joining (PHYLIP Neighbor), parsimony (PHYLIP Protpars), or maximum liklihood (PHYLIP Proml). The analyses were conducted by using all, or diverse subsets of bHLH sequences, all of which supported the Myc and Max clades. The Mnt/Mad clades alternatively group with a Myc/Max clade or to the Max clade with a Myc sister clade to that. Choanoflagellate bHLH sequences were identified by BLASTP to bHLH sequences from yeast and Drosophila. This process identified 3 choanoflagellate bHLH sequences (Fig. 3).

Genome conservation

Human, Drosophila, and Saccharomyces genomic data were retrieved for RiBi orthologs from Ensembl using their 1-to-1 orthology classification. Nematostella, and Monosiga single best BLAST matches (p < 10 - 5) were identified from the Joint Genome Institute data sets annotated by EST sequences and their genomic sequences retrieved. To ascertain that the absence of the E(CG) regulon in the yeast genome was not a secondarily derived trait akin to nematodes, we investigated other fungal genomes as well as more distantly related eukaryotic genomes, including at 10 other sequenced fungal genomes (JGI), as well as the ciliated protist Tetrahymena thermophilia, the single celled green alga Chlamydomonas reinhardtii (JGI), and the poplar tree Populus trichocarpa (JGI). Each genome was searched for BLASTP alignments using the fly Myc and Max bHLH amino acid sequences. Only the Tetrahymena genome produced a match, which upon various CLUSTALW alignments falls outside of the Group B superclade from which the Myc and Max clades group.

Control Sequences

Control genes (C1) were obtained by finding neighboring downstream genes preserving orthology [see Additional file 2]. First, yeast control genes were identified by finding for each gene with a nucleolar GO term, the next downstream gene with an ortholog in the human genome. Orthology with a human gene ensured that control sequences were derived from genes that are as conserved as RiBi genes. Fly and worm control genes were generated by finding for each human RiBi ortholog, the nearest downstream gene a human ortholog. Monosiga and Nematostella downstream controls (C1) were generated by finding the nearest downstream EST to the corresponding RiBi gene. The presence of one or more E(CG) sites in the region ± 600 bp from the 5'-most annotated end was ascertained for each such control gene. Additional control sequences (C2) were obtained for the Nematostella and Monosiga genomes by taking ± 600 bp from the 3' annotated end of each test gene [see Additional file 1].

MEME analysis and sequence logos

A set of RiBi orthologs meeting orthology gene tree tests between genomes were retrieved from the Ensembl data mining tool BioMart for Hs, Dr, Dm, Ce, and Sc genomes respectively. The annotation data used genome builds Hs (NCBI 36), Dr (Zv7), Dm (BDGP 5.4), CE (WS180) and Sc (SDG1.01). For identifying orthologous loci in Nv, Mb, and Ps, we used best reciprocal hits to identify orthologs in the respective genomes of the closest related species. A set of RiBi orthologs for Ag, Nc, and Sp were retrieved from the Ashbya Genome Database project base on Ensembl release 40. For identifying orthologous loci in Kl, Dh, Cg, Ca, and Yl we used ortholgy calls from the Génolevures Genomic Project. The DNA sequences for each species RiBi gene promoters were then collected from either Ensembl, JGI, or Génolevures for each respective genome. For fungal and choanoflagellate RiBi promoter sequences, -500 bp from the translational start site was collected for each ortholog. For Dm, Ce, Nv, Dr, and Hs, ± 250 bp from the 5' annotated end was collected for each ortholog. Each organism's RiBi sequence group was analyzed my MEME to determine overrepresented motifs [57, 58]. Motifs 6–10 bp in length that occured either zero or one time in each sequence of each set per species such that at least 10–75 motifs per set were identified. Furthermore, a secondary cut-off stipulating that all motifs be found in at least 2/3 of all sequences per set was applied. Sequence logos from the matrices of over-represented motifs derived from MEME were then created using WebLogo.

References

  1. Brown DD, Gurdon JB: Absence of ribosomal RNA synthesis in the anucleate mutant of Xenopus laevis. Proc Natl Acad Sci USA. 1964, 51: 139-146. 10.1073/pnas.51.1.139.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  2. Miller OLJ, Beatty BR: Visualization of nucleolar genes. Science. 1969, 164 (882): 955-957. 10.1126/science.164.3882.955.

    Article  PubMed  Google Scholar 

  3. Perry RP: On the nucleolar and nuclear dependence of cytoplasmic RNA synthesis in HeLa cells. Exp Cell Res. 1960, 20: 216-220. 10.1016/0014-4827(60)90240-8.

    Article  PubMed  CAS  Google Scholar 

  4. Erives A, Levine M: Coordinate enhancers share common organizational features in the Drosophila genome. Proc Natl Acad Sci USA. 2004, 101 (11): 3851-3856. 10.1073/pnas.0400611101.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Hughes JD, Estep PW, Tavazoie S, Church GM: Computational identification of cis-regulatory elements associated with groups of functionally related genes in Saccharomyces cerevisiae. J Mol Biol. 2000, 296 (5): 1205-1214. 10.1006/jmbi.2000.3519.

    Article  PubMed  CAS  Google Scholar 

  6. Tanay A, Regev A, Shamir R: Conservation and evolvability in regulatory networks: the evolution of ribosomal regulation in yeast. Proc Natl Acad Sci USA. 2005, 102 (20): 7203-7208. 10.1073/pnas.0502521102.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Dequard-Chablat M, Riva M, Carles C, Sentenac A: RPC19, the gene for a subunit common to yeast RNA polymerases A (I) and C (III). J Biol Chem. 1991, 266 (23): 15300-15307.

    PubMed  CAS  Google Scholar 

  8. Schreiber-Agus N, Stein D, Chen K, Goltz JS, Stevens L, DePinho RA: Drosophila Myc is oncogenic in mammalian cells and plays a role in the diminutive phenotype. Proc Natl Acad Sci USA. 1997, 94 (4): 1235-1240. 10.1073/pnas.94.4.1235.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  9. Johnston LA, Prober DA, Edgar BA, Eisenman RN, Gallant P: Drosophila myc regulates cellular growth during development. Cell. 1999, 98 (6): 779-790. 10.1016/S0092-8674(00)81512-3.

    Article  PubMed  CAS  Google Scholar 

  10. Gallant P, Shiio Y, Cheng PF, Parkhurst SM, Eisenman RN: Myc and Max homologs in Drosophila. Science. 1996, 274 (5292): 1523-1527. 10.1126/science.274.5292.1523.

    Article  PubMed  CAS  Google Scholar 

  11. Berberich S, Hyde-DeRuyscher N, Espenshade P, Cole M: max encodes a sequence-specific DNA-binding protein and is not regulated by serum growth factors. Oncogene. 1992, 7 (4): 775-779.

    PubMed  CAS  Google Scholar 

  12. Blackwell TK, Kretzner L, Blackwood EM, Eisenman RN, Weintraub H: Sequence-specific DNA binding by the c-Myc protein. Science. 1990, 250 (4984): 1149-1151. 10.1126/science.2251503.

    Article  PubMed  CAS  Google Scholar 

  13. Blackwood EM, Eisenman RN: Max: a helix-loop-helix zipper protein that forms a sequence-specific DNA-binding complex with Myc. Science. 1991, 251 (4998): 1211-1217. 10.1126/science.2006410.

    Article  PubMed  CAS  Google Scholar 

  14. Kerkhoff E, Bister K, Klempnauer KH: Sequence-specific DNA binding by Myc proteins. Proc Natl Acad Sci USA. 1991, 88 (10): 4323-4327. 10.1073/pnas.88.10.4323.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  15. Grandori C, Gomez-Roman N, Felton-Edkins ZA, Ngouenet C, Galloway DA, Eisenman RN, White RJ: c-Myc binds to human ribosomal DNA and stimulates transcription of rRNA genes by RNA polymerase I. Nat Cell Biol. 2005, 7 (3): 311-318. 10.1038/ncb1224.

    Article  PubMed  CAS  Google Scholar 

  16. Schlosser I, Holzel M, Murnseer M, Burtscher H, Weidle UH, Eick D: A role for c-Myc in the regulation of ribosomal RNA processing. Nucleic Acids Res. 2003, 31 (21): 6148-6156. 10.1093/nar/gkg794.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Greasley PJ, Bonnard C, Amati B: Myc induces the nucleolin and BN51 genes: possible implications in ribosome biogenesis. Nucleic Acids Res. 2000, 28 (2): 446-453. 10.1093/nar/28.2.446.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Zaffran S, Chartier A, Gallant P, Astier M, Arquier N, Doherty D, Gratecos D, Semeriva M: A Drosophila RNA helicase gene, pitchoune, is required for cell growth and proliferation and is a potential target of d-Myc. Development. 1998, 125 (18): 3571-3584.

    PubMed  CAS  Google Scholar 

  19. Hulf T, Bellosta P, Furrer M, Steiger D, Svensson D, Barbour A, Gallant P: Whole-genome analysis reveals a strong positional bias of conserved dMyc-dependent E-boxes. Mol Cell Biol. 2005, 25 (9): 3401-3410. 10.1128/MCB.25.9.3401-3410.2005.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  20. Grewal SS, Li L, Orian A, Eisenman RN, Edgar BA: Myc-dependent regulation of ribosomal RNA synthesis during Drosophila development. Nat Cell Biol. 2005, 7 (3): 295-302. 10.1038/ncb1223.

    Article  PubMed  CAS  Google Scholar 

  21. Giordano E, Peluso I, Senger S, Furia M: minifly, a Drosophila gene required for ribosome biogenesis. J Cell Biol. 1999, 144 (6): 1123-1133. 10.1083/jcb.144.6.1123.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Ohler U, Liao Gc, Niemann H, Rubin GM: Computational analysis of core promoters in the Drosophila genome. Genome Biol. 2002, 3 (12): RESEARCH0087-10.1186/gb-2002-3-12-research0087.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Orian A, van Steensel B, Delrow J, Bussemaker HJ, Li L, Sawado T, Williams E, Loo LWM, Cowley SM, Yost C, Pierce S, Edgar BA, Parkhurst SM, Eisenman RN: Genomic binding by the Drosophila Myc, Max, Mad/Mnt transcription factor network. Genes Dev. 2003, 17 (9): 1101-1114. 10.1101/gad.1066903.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  24. Hochheimer A, Zhou S, Zheng S, Holmes MC, Tjian R: TRF2 associates with DREF and directs promoter-selective gene expression in Drosophila. Nature. 2002, 420 (6914): 439-445. 10.1038/nature01167.

    Article  PubMed  CAS  Google Scholar 

  25. Ansmant I, Massenet S, Grosjean H, Motorin Y, Branlant C: Identification of the Saccharomyces cerevisiae RNA:pseudouridine synthase responsible for formation of psi(2819) in 21S mitochondrial ribosomal RNA. Nucleic Acids Res. 2000, 28 (9): 1941-1946. 10.1093/nar/28.9.1941.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Stolc V, Altman S: Rpp1, an essential protein subunit of nuclear RNase P required for processing of precursor tRNA and 35S precursor rRNA in Saccharomyces cerevisiae. Genes Dev. 1997, 11 (21): 2926-2937. 10.1101/gad.11.21.2926.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  27. Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, Terry A, Shapiro H, Lindquist E, Kapitonov VV, Jurka J, Genikhovich G, Grigoriev IV, Lucas SM, Steele RE, Finnerty JR, Technau U, Martindale MQ, Rokhsar DS: Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007, 317 (5834): 86-94. 10.1126/science.1139158.

    Article  PubMed  CAS  Google Scholar 

  28. King N, Westbrook MJ, Young SL, Kuo A, Abedin M, Chapman J, Fairclough S, Hellsten U, Isogai Y, Letunic I, Marr M, Pincus D, Putnam N, Rokas A, Wright KJ, Zuzow R, Dirks W, Good M, Goodstein D, Lemons D, Li W, Lyons JB, Morris A, Nichols S, Richter DJ, Salamov A, Sequencing JGI, Bork P, Lim WA, Manning G, Miller WT, McGinnis W, Shapiro H, Tjian R, Grigoriev IV, Rokhsar D: The genome of the choanoflagellate Monosiga brevicollis and the origin of metazoans. Nature. 2008, 451 (7180): 783-788. 10.1038/nature06617.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  29. Wade CH, Umbarger MA, McAlear MA: The budding yeast rRNA and ribosome biosynthesis (RRB) regulon contains over 200 genes. Yeast. 2006, 23 (4): 293-306. 10.1002/yea.1353.

    Article  PubMed  CAS  Google Scholar 

  30. Jorgensen P, Nishikawa JL, Breitkreutz BJ, Tyers M: Systematic identification of pathways that couple cell growth and division in yeast. Science. 2002, 297 (5580): 395-400. 10.1126/science.1070850.

    Article  PubMed  CAS  Google Scholar 

  31. Jorgensen P, Rupes I, Sharom JR, Schneper L, Broach JR, Tyers M: A dynamic transcriptional network communicates growth potential to ribosome synthesis and critical cell size. Genes Dev. 2004, 18 (20): 2491-2505. 10.1101/gad.1228804.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  32. Liu S, Milne GT, Kuremsky JG, Fink GR, Leppla SH: Identification of the proteins required for biosynthesis of diphthamide, the target of bacterial ADP-ribosylating toxins on translation elongation factor 2. Mol Cell Biol. 2004, 24 (21): 9487-9497. 10.1128/MCB.24.21.9487-9497.2004.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  33. Mattheakis LC, Shen WH, Collier RJ: DPH5, a methyltransferase gene required for diphthamide biosynthesis in Saccharomyces cerevisiae. Mol Cell Biol. 1992, 12 (9): 4026-4037.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Ortiz PA, Ulloque R, Kihara GK, Zheng H, Kinzy TG: Translation elongation factor 2 anticodon mimicry domain mutants affect fidelity and diphtheria toxin resistance. J Biol Chem. 2006, 281 (43): 32639-32648. 10.1074/jbc.M607076200.

    Article  PubMed  CAS  Google Scholar 

  35. Boyd KE, Farnham PJ: Myc versus USF: discrimination at the cad gene is determined by core promoter elements. Mol Cell Biol. 1997, 17 (5): 2529-2537.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  36. Marinkovic D, Marinkovic T, Kokai E, Barth T, Moller P, Wirth T: Identification of novel Myc target genes with a potential role in lymphomagenesis. Nucleic Acids Res. 2004, 32 (18): 5368-5378. 10.1093/nar/gkh877.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  37. Simionato E, Ledent V, Richards G, Thomas-Chollier M, Kerner P, Coornaert D, Degnan BM, Vervoort M: Origin and diversification of the basic helix-loop-helix gene family in metazoans: insights from comparative genomics. BMC Evol Biol. 2007, 7: 33-10.1186/1471-2148-7-33.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Ledent V, Paquet O, Vervoort M: Phylogenetic analysis of the human basic helix-loop-helix proteins. Genome Biol. 2002, 3 (6): RESEARCH0030-10.1186/gb-2002-3-6-research0030.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Robinson KA, Lopes JM: SURVEY AND SUMMARY: Saccharomyces cerevisiae basic helix-loop-helix proteins regulate diverse biological processes. Nucleic Acids Res. 2000, 28 (7): 1499-1505. 10.1093/nar/28.7.1499.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  40. Herman MA: Hermaphrodite cell-fate specification. WormBook. 2006, 1-16.

    Google Scholar 

  41. Blumenthal T, Thomas J: Cis and trans mRNA splicing in C. elegans. Trends Genet. 1988, 4 (11): 305-308. 10.1016/0168-9525(88)90107-2.

    Article  PubMed  CAS  Google Scholar 

  42. Blumenthal T: Trans-splicing and operons. WormBook. 2005, 1-9.

    Google Scholar 

  43. Cowling VH, Cole MD: HATs off to capping: a new mechanism for Myc. Cell Cycle. 2007, 6 (8): 907-909.

    Article  PubMed  CAS  Google Scholar 

  44. Cowling VH, Cole MD: The Myc transactivation domain promotes global phosphorylation of the RNA polymerase II carboxy-terminal domain independently of direct DNA binding. Mol Cell Biol. 2007, 27 (6): 2059-2073. 10.1128/MCB.01828-06.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  45. Gasch AP, Moses AM, Chiang DY, Fraser HB, Berardini M, Eisen MB: Conservation and evolution of cis-regulatory systems in ascomycete fungi. PLoS Biol. 2004, 2 (12): e398-10.1371/journal.pbio.0020398.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Hogues H, Lavoie H, Sellam A, Mangos M, Roemer T, Purisima E, Nantel A, Whiteway M: Transcription factor substitution during the evolution of fungal ribosome regulation. Mol Cell. 2008, 29 (5): 552-562. 10.1016/j.molcel.2008.02.006.

    Article  PubMed  CAS  Google Scholar 

  47. Pierce SB, Yost C, Anderson SAR, Flynn EM, Delrow J, Eisenman RN: Drosophila growth and development in the absence of dMyc and dMnt. Dev Biol. 2008, 315 (2): 303-316. 10.1016/j.ydbio.2007.12.026.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  48. Liko D, Slattery MG, Heideman W: Stb3 binds to ribosomal RNA processing element motifs that control transcriptional responses to growth in Saccharomyces cerevisiae. J Biol Chem. 2007, 282 (36): 26623-26628. 10.1074/jbc.M704762200.

    Article  PubMed  CAS  Google Scholar 

  49. Moore AW, Barbel S, Jan LY, Jan YN: A genomewide survey of basic helix-loop-helix factors in Drosophila. Proc Natl Acad Sci USA. 2000, 97 (19): 10436-10441. 10.1073/pnas.170301897.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  50. Ramsay NA, Glover BJ: MYB-bHLH-WD40 protein complex and the evolution of cellular diversity. Trends Plant Sci. 2005, 10 (2): 63-70. 10.1016/j.tplants.2004.12.011.

    Article  PubMed  CAS  Google Scholar 

  51. King N: The unicellular ancestry of animal development. Dev Cell. 2004, 7 (3): 313-325. 10.1016/j.devcel.2004.08.010.

    Article  PubMed  CAS  Google Scholar 

  52. King N, Carroll SB: A receptor tyrosine kinase from choanoflagellates: molecular insights into early animal evolution. Proc Natl Acad Sci USA. 2001, 98 (26): 15032-15037. 10.1073/pnas.261477698.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  53. Kelly K, Cochran BH, Stiles CD, Leder P: Cell-specific regulation of the c-myc gene by lymphocyte mitogens and platelet-derived growth factor. Cell. 1983, 35 (3 Pt 2): 603-610. 10.1016/0092-8674(83)90092-2.

    Article  PubMed  CAS  Google Scholar 

  54. Peterson KJ, Lyons JB, Nowak KS, Takacs CM, Wargo MJ, McPeek MA: Estimating metazoan divergence times with a molecular clock. Proc Natl Acad Sci USA. 2004, 101 (17): 6536-6541. 10.1073/pnas.0401670101.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  55. Berriz GF, King OD, Bryant B, Sander C, Roth FP: Characterizing gene sets with FuncAssociate. Bioinformatics. 2003, 19 (18): 2502-2504. 10.1093/bioinformatics/btg363.

    Article  PubMed  CAS  Google Scholar 

  56. Fromont-Racine M, Senger B, Saveanu C, Fasiolo F: Ribosome assembly in eukaryotes. Gene. 2003, 313: 17-42. 10.1016/S0378-1119(03)00629-2.

    Article  PubMed  CAS  Google Scholar 

  57. Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.

    PubMed  CAS  Google Scholar 

  58. Bailey TL, Williams N, Misleh C, Li WW: MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, W369-73. 10.1093/nar/gkl198. 34 Web Server

  59. Dujon B, Sherman D, Fischer G, Durrens P, Casaregola S, Lafontaine I, De Montigny J, Marck C, Neuveglise C, Talla E, Goffard N, Frangeul L, Aigle M, Anthouard V, Babour A, Barbe V, Barnay S, Blanchin S, Beckerich JM, Beyne E, Bleykasten C, Boisrame A, Boyer J, Cattolico L, Confanioleri F, De Daruvar A, Despons L, Fabre E, Fairhead C, Ferry-Dumazet H, Groppi A, Hantraye F, Hennequin C, Jauniaux N, Joyet P, Kachouri R, Kerrest A, Koszul R, Lemaire M, Lesur I, Ma L, Muller H, Nicaud JM, Nikolski M, Oztas S, Ozier-Kalogeropoulos O, Pellenz S, Potier S, Richard GF, Straub ML, Suleau A, Swennen D, Tekaia F, Wesolowski-Louvel M, Westhof E, Wirth B, Zeniou-Meyer M, Zivanovic I, Bolotin-Fukuhara M, Thierry A, Bouchier C, Caudron B, Scarpelli C, Gaillardin C, Weissenbach J, Wincker P, Souciet JL: Genome evolution in yeasts. Nature. 2004, 430 (6995): 35-44. 10.1038/nature02579.

    Article  PubMed  Google Scholar 

  60. Fitzpatrick DA, Logue ME, Stajich JE, Butler G: A fungal phylogeny based on 42 complete genomes derived from supertree and combined gene analysis. BMC Evol Biol. 2006, 6: 99-10.1186/1471-2148-6-99.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Jeffries TW, Grigoriev IV, Grimwood J, Laplaza JM, Aerts A, Salamov A, Schmutz J, Lindquist E, Dehal P, Shapiro H, Jin YS, Passoth V, Richardson PM: Genome sequence of the lignocellulose-bioconverting and xylose-fermenting yeast Pichia stipitis. Nat Biotechnol. 2007, 25 (3): 319-326. 10.1038/nbt1290.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We thank Jonathan Brown for help in statistical analyses. We thank John M. Wallace from the Dartmouth College Research Computing group for consulting on computer programming. We thank Mark McPeek and Kevin Peterson for comments on the manuscript. We also thank anonymous reviewers who have helped us with our manuscript. This work was supported by a Pre-doctoral Training Grant in Cancer Biology and Carcino-genesis (NCI/NIH) to S.J.B., a grant from the National Cancer Institute to M.D.C., and a start-up grant from Dartmouth College to A.J.E.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Michael D Cole or Albert J Erives.

Additional information

Authors' contributions

SJB, MDC, and AJE conceived and designed the bioinformatic and analytic strategies; AJE worked on the fly signature and whole fly genome queries; SJB worked on documenting the presence of the E(CG) signature in RiBi and control gene groups across eukaryotic genomes; both SJB and AJE made the final figures and tables; SJB, MDC, and AJE contributed to the writing and discussion of the paper.

Electronic supplementary material

12864_2008_1635_MOESM1_ESM.xls

Additional file 1: Catalog of E(CG)-bearing promoters by genome and gene function. All fly genes matching the fly RiBi promoter architecture (see Methods), or else constituting the fly ortholog to a yeast nucleolar or RiBi gene is listed. For each such fly gene, the corresponding ortholog (human, nematode, and yeast) or closest 1-to-1 homolog (cnidarian, choanoflagellate) is listed along with the presence (1, red) or absence (0) of E(CG) within 600 bp of the annotated start site. The E values scores for BLASTP matches in the Nematostella and Monosiga genomes are given in parentheses for all matches with E <e - 5. Matches below e-10 are highlighted in yellow. Boxes containing N/A indicate the presence of the cis-element could not be addressed because a distinct 1-to-1 ortholog could not be identified. Red highlighted yeast genes were listed as RiBi genes in a recent comprehensive review of eukaryotic ribosome biogenesis [56]. C2 control data sets are shown where applicable. (XLS 178 KB)

12864_2008_1635_MOESM2_ESM.xls

Additional file 2: C 1 control data sets for Figure 2. The presence (presence = "1", red) or absence (absence = "0") of promoter-proximal (± 600 bp) E(CG) sites within downstream control genomic loci (C1) is indicated. (XLS 73 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Brown, S.J., Cole, M.D. & Erives, A.J. Evolution of the holozoan ribosome biogenesis regulon. BMC Genomics 9, 442 (2008). https://doi.org/10.1186/1471-2164-9-442

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-9-442

Keywords