Comprehensive analysis of the base composition around the transcription start site in Metazoa

Cytosine methylation and DNA repair

Cytosine methylation is a common form of post-replicative DNA modification seen in both bacteria and eukaryotes. Modified cytosines have long been known to act as hotspots for mutations due to the high rate of spontaneous deamination of this base to thymine, resulting in a G/T mismatch. This will be fixed as a C-->T transition after replication if not repaired by the base excision repair (BER) pathway or specific repair enzymes dedicated to this purpose. This hypermutability has led to depletion of the target dinucleotide CpG outside of special CpG islands in mammals, which are normally unmethylated. We review the importance of C-->T transitions at non-island CpGs in human disease: When these occur in the germline, they are a common cause of inherited diseases such as epidermolysis bullosa and mucopolysaccharidosis, while in the soma they are frequently found in the genes for tumor suppressors such as p53 and the retinoblastoma protein, causing cancer. We also examine the specific repair enzymes involved, namely the endonuclease Vsr in Escherichia coli and two members of the uracil DNA glycosylase (UDG) superfamily in mammals, TDG and MBD4. Repair brings its own problems, since it will require remethylation of the replacement cytosine, presumably coupling repair to methylation by either the maintenance methylase Dnmt1 or a de novo enzyme such as Dnmt3a. Uncoupling of methylation from repair may be one way to remove methylation from DNA. We also look at the possible role of specific cytosine deaminases such as Aid and Apobec in accelerating deamination of methylcytosine and consequent DNA demethylation.

Developmentally regulated DNA methylation in Dictyostelium discoideum

Methylation of cytosine residues in DNA plays a critical role in the silencing of gene expression, organization of chromatin structure, and cellular differentiation of eukaryotes. Previous studies failed to detect 5-methylcytosine in Dictyostelium genomic DNA, but the recent sequencing of the Dictyostelium genome revealed a candidate DNA methyltransferase gene (dnmA). The genome sequence also uncovered an unusual distribution of potential methylation sites, CpG islands, throughout the genome. DnmA belongs to the Dnmt2 subfamily and contains all the catalytic motifs necessary for cytosine methyltransferases. Dnmt2 activity is typically weak in Drosophila melanogaster, mouse, and human cells and the gene function in these systems is unknown. We have investigated the methylation status of Dictyostelium genomic DNA with antibodies raised against 5-methylcytosine and detected low levels of the modified nucleotide. We also found that DNA methylation increased during development. We searched the genome for potential methylation sites and found them in retrotransposable elements and in several other genes. Using Southern blot analysis with methylation-sensitive and -insensitive restriction endonucleases, we found that the DIRS retrotransposon and the guaB gene were indeed methylated. We then mutated the dnmA gene and found that DNA methylation was reduced to about 50% of the wild-type level. The mutant cells exhibited morphological defects in late development, indicating that DNA methylation has a regulatory role in Dictyostelium development. Our findings establish a role for a Dnmt2 methyltransferase in eukaryotic development.

Large-scale structure of genomic methylation patterns

The mammalian genome depends on patterns of methylated cytosines for normal function, but the relationship between genomic methylation patterns and the underlying sequence is unclear. We have characterized the methylation landscape of the human genome by global analysis of patterns of CpG depletion and by direct sequencing of 3073 unmethylated domains and 2565 methylated domains from human brain DNA. The genome was found to consist of short (<4 kb) unmethylated domains embedded in a matrix of long methylated domains. Unmethylated domains were enriched in promoters, CpG islands, and first exons, while methylated domains comprised interspersed and tandem-repeated sequences, exons other than first exons, and non-annotated single-copy sequences that are depleted in the CpG dinucleotide. The enrichment of regulatory sequences in the relatively small unmethylated compartment suggests that cytosine methylation constrains the effective size of the genome through the selective exposure of regulatory sequences. This buffers regulatory networks against changes in total genome size and provides an explanation for the C value paradox, which concerns the wide variations in genome size that scale independently of gene number. This suggestion is compatible with the finding that cytosine methylation is universal among large-genome eukaryotes, while many eukaryotes with genome sizes <5 x 10(8) bp do not methylate their DNA.

Identification of highly specific localized sequence motifs in human ribosomal protein gene promoters

For ribosomal protein (RP) genes the start of transcription is rigidly controlled to maintain the 5'-TOP signal on the messenger RNA. The responsible regulatory mechanism is not yet fully understood. Careful comparative analysis of their proximal promoter sequences reveals common characteristics and thus provides clues to the underlying mechanism. We have extracted the proximal promoters of the 80 human cytosolic ribosomal protein genes together with the orthologous mouse sequences. After annotating the set with transcription factor binding sites based on the available literature, we searched for over-represented sequence motifs. We uncovered a novel motif that is localized at a fixed distance downstream to the transcription start. 31 out of the 80 promoters contain the motif in the same orientation around position +62 (standard deviation 6). A second evolutionary conserved and palindromic motif is found 13 times in the RP promoter set, 9 instances of which are located upstream around position -40. In addition, we see a characteristic profile of the GC-content and of the CpG dinucleotide frequencies. Our results support a model for the transcription of ribosomal protein genes in which the maintenance of the accurate start of transcription is provided by specific transcription factors. Such a factor binds the target DNA at a fixed location relative to the TSS, and possibly interacts directly with the basal transcription machinery.

Comparative genomics of trypanosomatid parasitic protozoa

A comparison of gene content and genome architecture of Trypanosoma brucei, Trypanosoma cruzi, and Leishmania major, three related pathogens with different life cycles and disease pathology, revealed a conserved core proteome of about 6200 genes in large syntenic polycistronic gene clusters. Many species-specific genes, especially large surface antigen families, occur at nonsyntenic chromosome-internal and subtelomeric regions. Retroelements, structural RNAs, and gene family expansion are often associated with syntenic discontinuities that-along with gene divergence, acquisition and loss, and rearrangement within the syntenic regions-have shaped the genomes of each parasite. Contrary to recent reports, our analyses reveal no evidence that these species are descended from an ancestor that contained a photosynthetic endosymbiont.

Begin at the beginning: predicting genes with 5' UTRs

The retrainable, comparative gene predictor N-SCAN integrates multigenome modeling and 5' untranslated region (5' UTR) modeling. In this article, we evaluate N-SCAN's transcription-start site (TSS) and first exon predictions both computationally and experimentally. The computational results indicate that N-SCAN is more accurate than any of the other tools we tested at predicting the TSS and the complete first exon. It is the only one of these tools that can predict complete gene structures together with 5' UTRs. Experimental evaluation shows that N-SCAN can be used to validate novel UTR introns in human gene predictions that do not overlap any RefSeq gene and even to correct RefSeq mRNAs by adding validated UTR exons that are missing from RefSeq.

GC-compositional strand bias around transcription start sites in plants and fungi

BACKGROUND

A GC-compositional strand bias or GC-skew (=(C-G)/(C+G)), where C and G denote the numbers of cytosine and guanine residues, was recently reported near the transcription start sites (TSS) of Arabidopsis genes. However, it is unclear whether other eukaryotic species have equally prominent GC-skews, and the biological meaning of this trait remains unknown.

RESULTS

Our study confirmed a significant GC-skew (C > G) in the TSS of Oryza sativa (rice) genes. The full-length cDNAs and genomic sequences from Arabidopsis and rice were compared using statistical analyses. Despite marked differences in the G+C content around the TSS in the two plants, the degrees of bias were almost identical. Although slight GC-skew peaks, including opposite skews (C < G), were detected around the TSS of genes in human and Drosophila, they were qualitatively and quantitatively different from those identified in plants. However, plant-like GC-skew in regions upstream of the translation initiation sites (TIS) in some fungi was identified following analyses of the expressed sequence tags and/or genomic sequences from other species. On the basis of our dataset, we estimated that > 70 and 68% of Arabidopsis and rice genes, respectively, had a strong GC-skew (> 0.33) in a 100-bp window (that is, the number of C residues was more than double the number of G residues in a +/-100-bp window around the TSS). The mean GC-skew value in the TSS of highly-expressed genes in Arabidopsis was significantly greater than that of genes with low expression levels. Many of the GC-skew peaks were preferentially located near the TSS, so we examined the potential value of GC-skew as an index for TSS identification. Our results confirm that the GC-skew can be used to assist the TSS prediction in plant genomes.

CONCLUSION

The GC-skew (C > G) around the TSS is strictly conserved between monocot and eudicot plants (ie. angiosperms in general), and a similar skew has been observed in some fungi. Highly-expressed Arabidopsis genes had overall a more marked GC-skew in the TSS compared to genes with low expression levels. We therefore propose that the GC-skew around the TSS in some plants and fungi is related to transcription. It might be caused by mutations during transcription initiation or the frequent use of transcription factor-biding sites having a strand preference. In addition, GC-skew is a good candidate index for TSS prediction in plant genomes, where there is a lack of correlation among CpG islands and genes.