Strong patterns in homooligomer tracts occurrences in non-coding and in potential regulatory sites in eukaryotic genomes

Previous studies of the dinucleotides flanking both the 5' and 3' ends of homooligomer tracts have shown that some flanks are consistently preferred over others (1,2). In the first preferred group, the homooligomer tracts are flanked by the same nucleotide and/or the complementary nucleotides, e.g.,ATAn,TTAn,CCGn, where n = 2-5. Runs flanked by nucleotides with which they cannot base pair are distinctly disfavored. (In this group An/Tn are flanked by C and/or G; Gn/Cn are flanked by A/T, e.g.,CGAn,TnGG,GnAT). The frequencies of runs flanked by A or T, and G or C ("mixed"group) are as expected. Here we seek the origin of this effect and its relevance to protein-DNA interactions. Surprisingly, within the first group, runs flanked by their complements with a pyrimidine-purine junction (e.g.,TTAn,CnGG) are greatly preferred. The frequencies of their purine-pyrimidine junction mirror-images is just as expected. This effect, as well as additional ones enumerated below, is seen universally in eukaryotes and in prokaryotes, although it is stronger in the former. Detailed analysis of regulatory regions shows these strong trends, particularly in GC sequences. The potential relationship to DNA conformation and DNA-protein interaction is discussed.     

SNPs occur in regions with less genomic sequence conservation

Rates of SNPs (single nucleotide polymorphisms) and cross-species genomic sequence conservation reflect intra- and inter-species variation, respectively. Here, I report SNP rates and genomic sequence conservation adjacent to mRNA processing regions and show that, as expected, more SNPs occur in less conserved regions and that functional regions have fewer SNPs. Results are confirmed using both mouse and human data. Regions include protein start codons, 3' splice sites, 5' splice sites, protein stop codons, predicted miRNA binding sites, and polyadenylation sites. Throughout, SNP rates are lower and conservation is higher at regulatory sites. Within coding regions, SNP rates are highest and conservation is lowest at codon position three and the fewest SNPs are found at codon position two, reflecting codon degeneracy for amino acid encoding. Exon splice sites show high conservation and very low SNP rates, reflecting both splicing signals and protein coding. Relaxed constraint on the codon third position is dramatically seen when separating exonic SNP rates based on intron phase. At polyadenylation sites, a peak of conservation and low SNP rate occurs from 30 to 17 nt preceding the site. This region is highly enriched for the sequence AAUAAA, reflecting the location of the conserved polyA signal. miRNA 3' UTR target sites are predicted incorporating interspecies genomic sequence conservation; SNP rates are low in these sites, again showing fewer SNPs in conserved regions. Together, these results confirm that SNPs, reflecting recent genetic variation, occur more frequently in regions with less evolutionarily conservation.     

A PUF family portrait: 3'UTR regulation as a way of life

In eukaryotic cells, mRNAs are exquisitely controlled, often through regulatory elements in their 3' untranslated regions (3'UTRs). Proteins that bind to those sites are key players in controlling mRNA stability, translation and localization. One family of regulatory proteins--the PUF proteins--are not only structurally related, but also bind to 3'UTRs and modulate mRNA expression in a wide variety of eukaryotic species. They do so either by enhancing turnover or repressing translation, and act combinatorially with other regulatory proteins. Here, we discuss the evolution, biological function and mechanisms of action of the PUF protein family, and suggest that a primordial function of PUF proteins is to sustain mitotic proliferation of stem cells.     

A uniquely conserved regulatory signal is found around the translation initiation site in three different collagen genes

We have determined the nucleotide sequence of the 5' untranslated region and the sequence encoding the signal peptide for mRNAs of the chick alpha 1 type I and alpha 1 type III collagen. These sequences were obtained by synthesizing the corresponding cDNAs using as primers either a synthetic oligonucleotide to prime alpha 1 type I cDNA or a DNA fragment isolated from a genomic clone coding for alpha 1 type III collagen to prime the cognate cDNA. Both primers were selected so that the resulting cDNAs would be short and would contain sequence information for the 5' untranslated region and the signal peptide of the proteins. The nucleotide sequences of these cDNAs were compared with the corresponding sequence of alpha 2 type I collagen. In each mRNA the 5' untranslated segment is approximately 130 nucleotides and contains two or more AUG triplets preceding the AUG which serves as a translation initiation codon. A sequence of about 50 nucleotides surrounding the translation initiation codon is remarkably conserved in all three mRNAs, whereas the sequences preceding and following this segment diverge markedly. This homologous sequence contains an almost identical inverted repeat sequence which could form a stable stem-loop structure. The initiation codon and the AUG which precedes it are found at the same place within this symmetrical sequence and the distance between them is invariant. The rest of the conserved sequence shows a less perfect symmetry. This conserved sequence has not been found in other genes. Our data suggest that these three and perhaps other collagen genes contain an identical regulatory signal that may play a role in determining the level of expression of these genes by modulating translational efficiency.     

Genome-wide functional analysis of human 5' untranslated region introns

Background  

Approximately 35% of human genes contain introns within the 5' untranslated region (UTR). Introns in 5'UTRs differ from those in coding regions and 3'UTRs with respect to nucleotide composition, length distribution and density. Despite their presumed impact on gene regulation, the evolution and possible functions of 5'UTR introns remain largely unexplored.     

A catalog of polymorphisms falling in microRNA-binding regions of cancer genes

Recent evidence indicates that small, nonprotein-coding RNA molecules, called microRNAs (miRNAs), control cell growth, differentiation, and apoptosis, and are also involved in tumorigenesis. miRNAs can bind to the 3' untranslated regions (3'UTRs) of messenger RNAs and interfere with their translation. We hypothesized that common polymorphisms within their genes or within their targets could have an important impact for an individual's risk to develop complex diseases. In this study, we selected the 3'UTRs of 129 genes involved in pathways commonly acknowledged as important for cancer, and we identified putative miRNA-binding sites by means of specialized algorithms (PicTar, DIANA-MicroT, miRBase, miRanda, TargetScan, and MicroInspector). Then we investigated 79 single-nucleotide polymorphisms (SNPs) within the putative binding sites for their ability to affect or impair the binding with the miRNA by assessing the DeltaDeltaG, the variation of DeltaG (Gibbs free energy), through comparing the wild-type and their corresponding variant alleles. Moreover, we reported seven identified SNPs in seven pre-miRNA hairpin regions and one SNP in the mature sequence of miR-608. Considering the validation status of the SNPs and their frequencies, we found at least 23 candidate polymorphisms of biological relevance that we propose for further investigation in case-control association studies.     

Conserved sequences in both coding and 5'' flanking regions of mammalian opal suppressor tRNA genes.

The rabbit genome encodes an opal suppressor tRNA gene. The coding region is strictly conserved between the rabbit gene and the corresponding gene in the human genome. The rabbit opal suppressor gene contains the consensus sequence in the 3' internal control region but like the human and chicken genes, the rabbit 5' internal control region contains two additional nucleotides. The 5' flanking sequences of the rabbit and the human opal suppressor genes contain extensive regions of homology. A subset of these homologies is also present 5' to the chicken opal suppressor gene. Both the rabbit and the human genomes also encode a pseudogene. That of the rabbit lacks the 3' half of the coding region. Neither pseudogene has homologous regions to the 5' flanking regions of the genes. The presence of 5' homologies flanking only the transcribed genes and not the pseudogenes suggests that these regions may be regulatory control elements specifically involved in the expression of the eukaryotic opal suppressor gene. Moreover the strict conservation of coding sequences indicates functional importance for the opal suppressor tRNA genes.     

Efficient translation in chloroplasts requires element(s) upstream of the putative ribosome binding site from atpI

Thousands of proteins make up a chloroplast, but fewer than 100 are encoded by the chloroplast genome. Despite this low number, expression of chloroplast-encoded genes is essential for plant survival. Every chloroplast has its own gene expression system with a major regulatory point at the initiation of protein synthesis (translation). In chloroplasts, most protein-encoding genes contain elements resembling the ribosome binding sites (RBS) found in prokaryotes. In vitro, these putative chloroplast ribosome binding sequences vary in their ability to support translation. Here we report results from an investigation into effects of the predicted RBS for the tobacco chloroplast atpI gene on translation in vivo. Two reporter constructs, differing only in their 5'-untranslated regions (5'UTRs) were stably incorporated into tobacco chloroplast genomes and their expression analyzed. One 5'UTR was derived from the wild-type (WT) atpI gene. The second, Holo-substitution (Holo-sub), had nonchloroplast sequence replacing all wild-type nucleotides, except for the putative RBS. The abundance of reporter RNA was the same for both 5'UTRs. However, translation controlled by Holo-sub was less than 4% that controlled by WT. These in vivo experiments support the idea that translation initiation in land plant chloroplasts depends on 5'UTR elements outside the putative RBS.     

Cooperative enhancement at the Drosophila Sgs-3 locus

The Drosophila glue gene Sgs-3 is specifically expressed in the secretory cells of the salivary glands of third instar larvae. We have assayed the expression of gene fusions to determine the role of cis-acting Sgs-3 sequences in conferring this pattern of expression. These experiments define two regulatory regions required for expression of reporter genes from the Sgs-3 promoter. One region, between 106 and 56 bp upstream of the Sgs-3 mRNA 5' end is sufficient for low but correct tissue- and stage-specific expression. A second region, lying between 629 and 130 bp 5' of the RNA start site is functionally equivalent; that is, it alone will also direct low level, specific expression. These two regions act synergistically to give high level expression. More distant upstream regions function to further increase levels of expression. These two regulatory elements can confer a salivary gland-specific pattern of expression on a heterologous promoter and are also sufficient to drive gene expression in other Drosophila species, implying conservation of regulators.     

Evidence that the respiratory syncytial virus polymerase is recruited to nucleotides 1 to 11 at the 3' end of the nucleocapsid and can scan to access internal signals

The 3'-terminal end of the respiratory syncytial virus genomic RNA contains a 44-nucleotide leader (Le) region adjoining the gene start signal of the first gene. Previous mapping studies demonstrated that there is a promoter located at the 3' end of Le, which can signal initiation of antigenome synthesis. The aim of this study was to investigate the role of the 3' terminus of the RNA template in (i) promoter recognition and (ii) determining the initiation site for antigenome synthesis. A panel of minigenomes containing additional sequence at the 3' end of the Le were analyzed for their ability to direct antigenome and mRNA synthesis. Minigenomes containing heterologous extensions of 6 nucleotides or more were unable to support efficient RNA synthesis. However, the activity of a minigenome with a 56-nucleotide extension could be restored by insertion of Le nucleotides 1 to 11 or 1 to 13 at the 3' end, indicating that these nucleotides, in conjunction with the 3' terminus, are sufficient to recruit polymerase to the template. Northern blot and 5' rapid amplification of cDNA ends analysis of antigenome RNA indicated that antigenome initiation occurred at the first position of Le, irrespective of the terminal extension. This finding demonstrates that the 3' terminus of the RNA is not necessary for determining the antigenome initiation site. Data are presented which suggest that following recruitment to a promoter at the 3' end of Le, the polymerase is able to scan and respond to a promoter signal embedded within the RNA template.     

RNA-binding protein Dnd1 inhibits microRNA access to target mRNA

MicroRNAs (miRNAs) are inhibitors of gene expression capable of controlling processes in normal development and cancer. In mammals, miRNAs use a seed sequence of 6-8 nucleotides (nt) to associate with 3' untranslated regions (3'UTRs) of mRNAs and inhibit their expression. Intriguingly, occasionally not only the miRNA-targeting site but also sequences in its vicinity are highly conserved throughout evolution. We therefore hypothesized that conserved regions in mRNAs may serve as docking platforms for modulators of miRNA activity. Here we demonstrate that the expression of dead end 1 (Dnd1), an evolutionary conserved RNA-binding protein (RBP), counteracts the function of several miRNAs in human cells and in primordial germ cells of zebrafish by binding mRNAs and prohibiting miRNAs from associating with their target sites. These effects of Dnd1 are mediated through uridine-rich regions present in the miRNA-targeted mRNAs. Thus, our data unravel a novel role of Dnd1 in protecting certain mRNAs from miRNA-mediated repression.