Evolution of CpG islands within the myc gene family

CpG islands are discrete regions of DNA with significantly greater frequencies of CpG doublets than bulk genomic DNA. They are most frequently associated with the 5'-ends of housekeeping genes and are involved in the regulation of their expression. In this study, the structure and evolution of CpG islands within genes of the myc family were evaluated with the protein-coding sequences of animals and their transducing viruses. These evaluations relied on a gene tree for the entire myc family to test the origins of CpG islands within their two protein-coding exons. Overall, CG-very rich and CG-rich islands are associated with exon 2 of the different myc genes of warm-blooded vertebrates and with exon 3 of the N-myc and s-myc sequences of mammals, but not birds. These overall distributions of well-developed islands can be related to the major transitions of the CG-rich genomes of warm-blooded vertebrates from the CG-poor ones of other animals. In turn, the greater variability of well-developed islands within exon 3 of the N-myc gene and among the different retrogenes of the myc family can be attributed to their reduced functional constraints, as evidenced by their limited and very restricted patterns of expression, respectively.     

CpG islands: features and distribution in the genomes of vertebrates

We have investigated the distribution of unmethylated CpG islands in vertebrate genomes fractionated according to their base composition. Genomes from warm-blooded vertebrates (man, mouse and chicken) are characterized by abundant CpG islands, whose frequency increases in DNA fractions of increasing % of guanine + cytosine; % G + C (GC), in parallel with the distribution of genes and CpG doublets. Small, yet significant, differences in the distribution of CpG islands were found in the three genomes. In contrast, genomes from cold-blooded vertebrates (two reptiles, one amphibian, and two fishes) were characterized by an extreme scarcity or absence of CpG islands (detected in these experiments as HpaII tiny fragments or HTF). CpG islands associated with homologous genes from cold- and warm-blooded vertebrates were then compared by analyzing CpG frequencies, GC levels, HpaII sites, rare-cutter sites and G/C boxes (GGGGCGGGGC and closely related motifs) in sequences available in gene banks. Small, yet significant, differences were again detected among the CpG islands associated with homologous genes from warm-blooded vertebrates, in that CpG islands associated with mouse or rat genes often showed low CpG and/or GC levels, as well as low numbers of HpaII sites, rare-cutter sites and G/C boxes, compared to homologous human genes; more rarely, CpG islands were just absent. As far as cold-blooded vertebrates were concerned, a number of genes showed CpG islands, which exhibited a much lower frequency of CpG doublets than that found in CpG islands of warm-blooded vertebrates, but still approached the statistically expected frequency; none of the other features of CpG islands associated with genes from warm-blooded vertebrates were present. Other genes did not show any associated CpG islands, unlike their homologues from warm-blooded vertebrates.     

The amplification of oligonucleotide themes in the evolution of themyc protooncogene family

The evolutionary past of intragenic repeats in protein-coding exons of c-, N-, L-, and s-myc-protooncogene subfamilies was elucidated. Apparently these genes evolved by succession of distinct unit events rather than by a steady flow of random point mutations. An evolutionary event probably involved a duplication of the whole gene, which was followed by amplification of progressively shorter oligonucleotide themes and motifs. The repeats were either joined in tandem or one of the copies was transposed and integrated elsewhere within the same exon. In some instances multiple fragments of an amplified theme were integrated at several sites. Direct repeats were found to prevail over inverted ones. By reconstructing the fate of repeats in the course of evolution of vertebrates, the origins of some functional domains could be traced to the initial amplification event. For example, an N-myc-specific domain was created by tandem duplication of a single-copy theme of L-myc exon; at the time of divergence of the c-myc and N-myc, the tandem duplex underwent a new round of duplication followed by transposition of the new copy, thus accounting for the formation of a new domain specific for c-myc. The model proposed here may be regarded as a molecular-level equivalent of the theory of punctuated equilibria. This is the author's last paper; it was submitted shortly before his death. The proofs were corrected by Michal Dvorak Correspondence to: M. Dvorak     

CpG islands, genes and isochores in the genomes of vertebrates

We have shown that human genes associated with CpG islands increase in number as they increase in % of guanine + cytosine (GC) levels, and that most genes associated with CpG islands are located in the GC-richest compartment of the human genome. This is an independent confirmation of the concentration gradient of CpG islands (detected as HpaII tiny fragments, or HTF) which was demonstrated in the genome of warm-blooded vertebrates [A?ssani and Bernardi, Gene 106 (1991) 173-183]. We then reassessed the location of CpG islands using the data currently available and confirmed that CpG islands are most frequently located in the 5'-flanking sequences of genes and that they overlap genes to variable extents. We have shown that such extents increase with the increasing GC levels of genes, the GC-richest genes being completely included in CpG islands. Under such circumstances, we have investigated the properties of the 'extragenic' CpG islands located in the 5'-flanking segments of homologous genes from both warm- and cold-blooded vertebrates. We have confirmed that, in cold-blooded vertebrates, CpG islands are often absent; when present, they have lower GC and CpG levels; the latter attain, however, statistically expected values. Finally, we have shown that CpG doublets increase with the increasing GC of exons, introns and intergenic sequences (including 'extragenic' CpG islands) in the genomes from both warm- and cold-blooded vertebrates. The correlations found are the same for both classes of vertebrates, and are similar for exons, introns and intergenic sequences (including 'extragenic' CpG islands). The findings just outlined indicate that the origin and evolution of CpG islands in the vertebrate genome are associated with compositional transitions (GC increases) in genes and isochores.     

The growth-inhibitory Ndrg1 gene is a Myc negative target in human neuroblastomas and other cell types with overexpressed N- or c-myc

A major prognostic marker for neuroblastoma (Nb) is N-myc gene amplification, which predicts a poor clinical outcome. We sought genes differentially expressed on a consistent basis between multiple human Nb cell lines bearing normal versus amplified N-myc, in hopes of finding target genes that might clarify how N-myc overexpression translates into poor clinical prognosis. Using differential display, we find the previously described growth-inhibitory gene Ndrg1 is strongly repressed in all tested Nb cell lines bearing N-myc amplification, as well as in a neuroepithelioma line with amplified c-myc. Overexpression of N-myc in non-amplified Nb cells leads to repression of Ndrg1, as does activation of an inducible c-myc transgene in fibroblasts. Conversely, N-myc downregulation in N-myc-amplified Nb cells results in re-expression of the Ndrg1, and stimuli known to induce Ndrg1 do so in Nb cells while simultaneously down-regulating N-myc. Relevant to these results, we demonstrate an in vitro interaction of Myc protein with the Ndrg1 core promoter. We also find that Ndrg1 levels increase dramatically during in vitro differentiation of two cell lines modeling neural and glial development, while c- and N-myc levels decline. Our results combined with previous information on the Ndrg1 gene product suggest that downregulation of this gene is an important component of N-Myc effects in neuroblastomas with poor clinical outcome. In support of this notion, we find that re-expression of Ndrg1 in high-Myc Nb cells results in smaller cells with reduced colony size in soft-agar assays, further underscoring the functional significance of this gene in human neuroblastoma cells.     

N-Myc Regulates Expression of Pluripotency Genes in Neuroblastoma Including lif,klf2, klf4, and lin28b

myc genes are best known for causing tumors when overexpressed, but recent studies suggest endogenous myc regulates pluripotency and self-renewal of stem cells. For example, N-myc is associated with a number of tumors including neuroblastoma, but also plays a central role in the function of normal neural stem and precursor cells (NSC). Both c- and N-myc also enhance the production of induced pluripotent stem cells (iPSC) and are linked to neural tumor stem cells. The mechanisms by which myc regulates normal and neoplastic stem-related functions remain largely open questions. Here from a global, unbiased search for N-Myc bound genes using ChIP-chip assays in neuroblastoma, we found lif as a putative N-Myc bound gene with a number of strong N-Myc binding peaks in the promoter region enriched for E-boxes. Amongst putative N-Myc target genes in expression microarray studies in neuroblastoma we also found lif and three additional important embryonic stem cell (ESC)-related factors that are linked to production of iPSC: klf2, klf4, and lin28b. To examine the regulation of these genes by N-Myc, we measured their expression using neuroblastoma cells that contain a Tet-regulatable N-myc transgene (TET21N) as well as NSC with a nestin-cre driven N-myc knockout. N-myc levels closely correlated with the expression of all of these genes in neuroblastoma and all but lif in NSC. Direct ChIP assays also indicate that N-Myc directly binds the lif promoter. N-Myc regulates trimethylation of lysine 4 of histone H3 in the promoter of lif and possibly in the promoters of several other stem-related genes. Together these findings indicate that N-Myc regulates overlapping stem-related gene expression programs in neuroblastoma and NSC, supporting a novel model by which amplification of the N-myc gene may drive formation of neuroblastoma. They also suggest mechanisms by which Myc proteins more generally contribute to maintenance of pluripotency and self-renewal of ESC as well as to iPSC formation.     

Non-methylated islands in fish genomes are GC-poor.

In the vertebrate genomes studied to date the 5' end of many genes are associated with distinctive sequences known as CpG islands. CpG islands have three properties: they are non-methylated; the dinucleotide CpG occurs at the frequency predicted by base composition; and they are GC-rich. Unexpectedly we have found that CpG islands in certain fish only have the first two properties; that is, their GC-content is not elevated compared to bulk genomic DNA. Based on this finding, we speculate that the GC-richness of CpG islands in vertebrates other than fish is a passive consequence of a higher mutation rate in regions of open chromatin under conditions where the nucleotide precursor pools are biased.     

Characterization of N-myc amplification in a human neuroblastoma cell line by clones isolated following the phenol emulsion reassociation technique and by hexagonal field gel electrophoresis

The N-myc amplification of human neuroblastomas was characterized by the amplified DNA cloned from the cell line MC-NB-1 using the phenol emulsion reassociation technique (PERT). A number of PERT clones exhibiting amplification in this cell line were tested for amplification in other neuroblastoma cell lines. In almost all cell lines examined, only a few clones were co-amplified with N-myc and most of the others were exclusively amplified in a subset of the cell lines. The total aggregate size of theHind III fragment identified by the PERT clones was approximately 350 kb. Most of the PERT clones were mapped to human chromosome (chr) 2p23-2pter, where the N-myc gene is located. Four types of amplicons, the 100, 420, 480 and 520 kb fragments, shown to beNot I fragments, were identified by hexagonal field gel electrophoresis. Three fragments are ordered in a head-to-tail array, and the remaining fragment is either ordered in a tail-to-head array or something else. Despite the extremely unusual construction of the amplified sequences in this cell line as compared with others, there was a low degree of sequence heterogeneity among the amplicons within this cell line. These observations lead to the idea that the complex rearrangements that give rise to the heterogeneous organization of the amplified sequences among the different cell lines precede the amplification of these sequences.     

Compositional compartmentalization and gene composition in the genome of vertebrates

Summary The compositional distribution of coding sequences from five vertebrates (Xenopus, chicken, mouse, rat, and human) is shifted toward higher GC values compared to that of the DNA molecules (in the 35–85-kb size range) isolated from the corresponding genomes. This shift is due to the lower GC levels of intergenic sequences compared to coding sequences. In the cold-blooded vertebrate, the two distributions are similar in that GC-poor genes and GC-poor DNA molecules are largely predominant. In contrast, in the warm-blooded vertebrates, GC-rich genes are largely predominant over GC-poor genes, whereas GC-poor DNA molecules are largely predominant over GC-rich DNA molecules. As a consequence, the genomes of warm-blooded vertebrates show a compositional gradient of gene concentration. The compositional distributions of coding sequences (as well as of DNA molecules) showed remarkable differences between chicken and mammals, and between mouse (or rat) and human. Differences were also detected in the compositional distribution of housekeeping and tissue-specific genes, the former being more abundant among GC-rich genes.     

Compositional properties of nuclear genes from cold-blooded vertebrates

Summary We have investigated the compositional properties of coding sequences from cold-blooded vertebrates and we have compared them with those from warm-blooded vertebrates. Moreover, we have studied the compositional correlations of coding sequences with the genomes in which they are contained, as well as the compositional correlations among the codon positions of the genes analyzed.The distribution of GC levels of the third codon positions of genes from cold-blooded vertebrates are distinctly different from those of warm-blooded vertebrates in that they do not reach the high values attained by the latter. Moreover, coding sequences from cold-blooded vertebrates are either equal, or, in most cases, lower in GC (not only in third, but also in first and second codon positions) than homologous coding sequences from warm-blooded vertebrates; higher values are exceptional. These results at the gene level are in agreement with the compositional differences between cold-blooded and warm-blooded vertebrates previously found at the whole genome (DNA) level (Bernardi and Bernardi 1990a,b).Two linear correlations were found: one between the GC levels of coding sequences (or of their third codon positions) and the GC levels of the genomes of cold-blooded vertebrates containing them; and another between the GC levels of third and first+ second codon positions of genes from cold-blooded vertebrates. The first correlation applies to the genomes (or genome compartments) of all vertebrates and the second to the genes of all living organisms. These correlations are tantamount to a genomic code.     

Gene-associated CpG islands in plants as revealed by analyses of genomic sequences

We screened plant genome sequences, primarily from rice and Arabidopsis thaliana, for CpG islands, and identified DNA segments rich in CpG dinucleotides within these sequences. These CpG-rich clusters appeared in the analysed sequences as discrete peaks and occurred at the frequencies of one per 4.7 kb in rice and one per 4.0 kb in A. thaliana. In rice and A. thaliana, most of the CpG-rich clusters were associated with genes, which suggests that these clusters are useful landmarks in genome sequences for identifying genes in plants with small genomes. In contrast, in plants with larger genomes, only a few of the clusters were associated with genes. These plant CpG-rich clusters satisfied the criteria used for identifying human CpG islands, which suggests that these CpG clusters may be regarded as plant CpG islands. The position of each island relative to the 5'-end of its associated gene varied considerably. Genes in the analysed sequences were grouped into five classes according to the position of the CpG islands within their associated genes. A large proportion of the genes belonged to one of two classes, in which a CpG island occurred near the 5'-end of the gene or covered the whole gene region. The position of a plant CpG island within its associated gene appeared to be related to the extent of tissue-specific expression of the gene; the CpG islands of most of the widely expressed rice genes occurred near the 5'-end of the genes.     

Thermal Adaptation of the Small Subunit Ribosomal RNA Gene: A Comparative Study

We carried out a comprehensive survey of small subunit ribosomal RNA sequences from archaeal, bacterial, and eukaryotic lineages in order to understand the general patterns of thermal adaptation in the rRNA genes. Within each lineage, we compared sequences from mesophilic, moderately thermophilic, and hyperthermophilic species. We carried out a more detailed study of the archaea, because of the wide range of growth temperatures within this group. Our results confirmed that there is a clear correlation between the GC content of the paired stem regions of the 16S rRNA genes and the optimal growth temperature, and we show that this correlation cannot be explained simply by phylogenetic relatedness among the thermophilic archaeal species. In addition, we found a significant, positive relationship between rRNA stem length and growth temperature. These correlations are found in both bacterial and archaeal rRNA genes. Finally, we compared rRNA sequences from warm-blooded and cold-blooded vertebrates. We found that, while rRNA sequences from the warm-blooded vertebrates have a higher overall GC content than those from the cold-blooded vertebrates, this difference is not concentrated in the paired regions of the molecule, suggesting that thermal adaptation is not the cause of the nucleotide differences between the vertebrate lineages. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Nicolas Galtier]