Patterns in interspecies similarity correlate with nucleotide composition in mammalian 3'UTRs

A massively parallel 3' UTR reporter assay reveals relationships between nucleotide content, sequence conservation, and mRNA destabilization

Compared to coding sequences, untranslated regions of the transcriptome are not well conserved, and functional annotation of these sequences is challenging. Global relationships between nucleotide composition of 3′ UTR sequences and their sequence conservation have been appreciated since mammalian genomes were first sequenced, but the functional relevance of these patterns remain unknown. We systematically measured the effect on gene expression of the sequences of more than 25,000 RNA-binding protein (RBP) binding sites in primary mouse T cells using a massively parallel reporter assay. GC-rich sequences were destabilizing of reporter mRNAs and come from more rapidly evolving regions of the genome. These sequences were more likely to be folded in vivo and contain a number of structural motifs that reduced accumulation of a heterologous reporter protein. Comparison of full-length 3′ UTR sequences across vertebrate phylogeny revealed that strictly conserved 3′ UTRs were GC-poor and enriched in genes associated with organismal development. In contrast, rapidly evolving 3′ UTRs tended to be GC-rich and derived from genes involved in metabolism and immune responses. Cell-essential genes had lower GC content in their 3′ UTRs, suggesting a connection between unstructured mRNA noncoding sequences and optimal protein production. By reducing gene expression, GC-rich RBP-occupied sequences act as a rapidly evolving substrate for gene regulatory interactions.

Comparative in silico analysis of EST-SSRs in angiosperm and gymnosperm tree genera

BACKGROUND

Simple Sequence Repeats (SSRs) derived from Expressed Sequence Tags (ESTs) belong to the expressed fraction of the genome and are important for gene regulation, recombination, DNA replication, cell cycle and mismatch repair. Here, we present a comparative analysis of the SSR motif distribution in the 5'UTR, ORF and 3'UTR fractions of ESTs across selected genera of woody trees representing gymnosperms (17 species from seven genera) and angiosperms (40 species from eight genera).

RESULTS

Our analysis supports a modest contribution of EST-SSR length to genome size in gymnosperms, while EST-SSR density was not associated with genome size in neither angiosperms nor gymnosperms. Multiple factors seem to have contributed to the lower abundance of EST-SSRs in gymnosperms that has resulted in a non-linear relationship with genome size diversity. The AG/CT motif was found to be the most abundant in SSRs of both angiosperms and gymnosperms, with a relative increase in AT/AT in the latter. Our data also reveals a higher abundance of hexamers across the gymnosperm genera.

CONCLUSIONS

Our analysis provides the foundation for future comparative studies at the species level to unravel the evolutionary processes that control the SSR genesis and divergence between angiosperm and gymnosperm tree species.

Hyper conserved elements in vertebrate mRNA 3'-UTRs reveal a translational network of RNA-binding proteins controlled by HuR

Little is known regarding the post-transcriptional networks that control gene expression in eukaryotes. Additionally, we still need to understand how these networks evolve, and the relative role played in them by their sequence-dependent regulatory factors, non-coding RNAs (ncRNAs) and RNA-binding proteins (RBPs). Here, we used an approach that relied on both phylogenetic sequence sharing and conservation in the whole mapped 3'-untranslated regions (3'-UTRs) of vertebrate species to gain knowledge on core post-transcriptional networks. The identified human hyper conserved elements (HCEs) were predicted to be preferred binding sites for RBPs and not for ncRNAs, namely microRNAs and long ncRNAs. We found that the HCE map identified a well-known network that post-transcriptionally regulates histone mRNAs. We were then able to discover and experimentally confirm a translational network composed of RNA Recognition Motif (RRM)-type RBP mRNAs that are positively controlled by HuR, another RRM-type RBP. HuR shows a preference for these RBP mRNAs bound in stem-loop motifs, confirming its role as a 'regulator of regulators'. Analysis of the transcriptome-wide HCE distribution revealed a profile of prevalently small clusters separated by unconserved intercluster RNA stretches, which predicts the formation of discrete small ribonucleoprotein complexes in the 3'-UTRs.

Multilevel regulation of HIF-1 signaling by TTP

Phosphorylation of the RNA-binding protein tristetraprolin (TTP) by p38 MAPK/MK2 does not prevent its RNA interaction and switches the mode of TTP action from destabilization to stabilization of the HIF-1α mRNA and subsequent activation of HIF-1 signaling.

Hypoxia-inducible factor-1 (HIF-1) is a well-studied transcription factor mediating cellular adaptation to hypoxia. It also plays a crucial role under normoxic conditions, such as in inflammation, where its regulation is less well understood. The 3′-untranslated region (UTR) of HIF-1α mRNA is among the most conserved UTRs in the genome, hinting toward posttranscriptional regulation. To identify potential trans factors, we analyzed a large compilation of expression data. In contrast to its known function of being a negative regulator, we found that tristetraprolin (TTP) positively correlates with HIF-1 target genes. Mathematical modeling predicts that an additional level of posttranslational regulation of TTP can explain the observed positive correlation between TTP and HIF-1 signaling. Mechanistic studies revealed that TTP indeed changes its mode of regulation from destabilizing to stabilizing HIF-1α mRNA upon phosphorylation by p38 mitogen-activated protein kinase (MAPK)/MAPK-activated protein kinase 2. Using a model of monocyte-to-macrophage differentiation, we show that TTP-driven HIF-1α mRNA stabilization is crucial for cell migration. This demonstrates the physiological importance of a hitherto-unknown mechanism for multilevel regulation of HIF-1α in normoxia.

Evolution of alternative and constitutive regions of mammalian 5'UTRs

BACKGROUND

Alternative splicing (AS) in protein-coding sequences has emerged as an important mechanism of regulation and diversification of animal gene function. By contrast, the extent and roles of alternative events including AS and alternative transcription initiation (ATI) within the 5'-untranslated regions (5'UTRs) of mammalian genes are not well characterized.

RESULTS

We evaluated the abundance, conservation and evolution of putative regulatory control elements, namely, upstream start codons (uAUGs) and open reading frames (uORFs), in the 5'UTRs of human and mouse genes impacted by alternative events. For genes with alternative 5'UTRs, the fraction of alternative sequences (those present in a subset of the transcripts) is much greater than that in the corresponding coding sequence, conceivably, because 5'UTRs are not bound by constraints on protein structure that limit AS in coding regions. Alternative regions of mammalian 5'UTRs evolve faster and are subject to a weaker purifying selection than constitutive portions. This relatively weak selection results in over-abundance of uAUGs and uORFs in the alternative regions of 5'UTRs compared to constitutive regions. Nevertheless, even in alternative regions, uORFs evolve under a stronger selection than the rest of the sequences, indicating that some of the uORFs are conserved regulatory elements; some of the non-conserved uORFs could be involved in species-specific regulation.

CONCLUSION

The findings on the evolution and selection in alternative and constitutive regions presented here are consistent with the hypothesis that alternative events, namely, AS and ATI, in 5'UTRs of mammalian genes are likely to contribute to the regulation of translation.

Expansion of the human mu-opioid receptor gene architecture: novel functional variants

The mu-opioid receptor (OPRM1) is the principal receptor target for both endogenous and exogenous opioid analgesics. There are substantial individual differences in human responses to painful stimuli and to opiate drugs that are attributed to genetic variations in OPRM1. In searching for new functional variants, we employed comparative genome analysis and obtained evidence for the existence of an expanded human OPRM1 gene locus with new promoters, alternative exons and regulatory elements. Examination of polymorphisms within the human OPRM1 gene locus identified strong association between single nucleotide polymorphism (SNP) rs563649 and individual variations in pain perception. SNP rs563649 is located within a structurally conserved internal ribosome entry site (IRES) in the 5'-UTR of a novel exon 13-containing OPRM1 isoforms (MOR-1K) and affects both mRNA levels and translation efficiency of these variants. Furthermore, rs563649 exhibits very strong linkage disequilibrium throughout the entire OPRM1 gene locus and thus affects the functional contribution of the corresponding haplotype that includes other functional OPRM1 SNPs. Our results provide evidence for an essential role for MOR-1K isoforms in nociceptive signaling and suggest that genetic variations in alternative OPRM1 isoforms may contribute to individual differences in opiate responses.

Posttranscriptional regulation of the breast cancer susceptibility gene BRCA1 by the RNA binding protein HuR

BRCA1 is a breast cancer susceptibility gene that is down-regulated in a significant proportion of sporadic breast cancers. BRCA1 is posttranscriptionally regulated by RNA-binding proteins, the identities of which are unknown. HuR is an RNA binding protein implicated in posttranscriptional regulation of many genes and is overexpressed in sporadic breast cancer. To investigate the possibility that these two molecules are functionally linked in breast cancer, we performed bioinformatic analysis of the BRCA1 3' untranslated region (UTR), RNA-protein assays with the HuR protein and the BRCA1 3'UTR, and immunohistochemical analysis of a cohort of breast tumors using antibodies against BRCA1 and HuR. Here, we describe the identification of two predicted HuR-binding sites in the BRCA1 3'UTR, one of which binds specifically to HuR. We also show that this interaction is disrupted by single nucleotide substitutions in the BRCA1 3'UTR and that endogenous HuR protein associates with BRCA1 transcripts in T47D and MCF7 breast cancer cells. Expression of ectopic HuR results in a significant decrease in BRCA1 protein expression and also BRCA1 3'UTR activity. Immunohistochemical analysis revealed that although BRCA1 and HuR expression were associated with some clinicopathologic features of the tumors, there was no statistically significant correlation between BRCA1 and HuR protein expression. These results identify the first posttranscriptional protein regulator of BRCA1 and have implications for understanding BRCA1 regulation in human breast cancer.

Expression patterns of protein kinases correlate with gene architecture and evolutionary rates

BACKGROUND

Protein kinase (PK) genes comprise the third largest superfamily that occupy approximately 2% of the human genome. They encode regulatory enzymes that control a vast variety of cellular processes through phosphorylation of their protein substrates. Expression of PK genes is subject to complex transcriptional regulation which is not fully understood.

PRINCIPAL FINDINGS

Our comparative analysis demonstrates that genomic organization of regulatory PK genes differs from organization of other protein coding genes. PK genes occupy larger genomic loci, have longer introns, spacer regions, and encode larger proteins. The primary transcript length of PK genes, similar to other protein coding genes, inversely correlates with gene expression level and expression breadth, which is likely due to the necessity to reduce metabolic costs of transcription for abundant messages. On average, PK genes evolve slower than other protein coding genes. Breadth of PK expression negatively correlates with rate of non-synonymous substitutions in protein coding regions. This rate is lower for high expression and ubiquitous PKs, relative to low expression PKs, and correlates with divergence in untranslated regions. Conversely, rate of silent mutations is uniform in different PK groups, indicating that differing rates of non-synonymous substitutions reflect variations in selective pressure. Brain and testis employ a considerable number of tissue-specific PKs, indicating high complexity of phosphorylation-dependent regulatory network in these organs. There are considerable differences in genomic organization between PKs up-regulated in the testis and brain. PK genes up-regulated in the highly proliferative testicular tissue are fast evolving and small, with short introns and transcribed regions. In contrast, genes up-regulated in the minimally proliferative nervous tissue carry long introns, extended transcribed regions, and evolve slowly.

CONCLUSIONS/SIGNIFICANCE

PK genomic architecture, the size of gene functional domains and evolutionary rates correlate with the pattern of gene expression. Structure and evolutionary divergence of tissue-specific PK genes is related to the proliferative activity of the tissue where these genes are predominantly expressed. Our data provide evidence that physiological requirements for transcription intensity, ubiquitous expression, and tissue-specific regulation shape gene structure and affect rates of evolution.

Identification of candidate regulatory sequences in mammalian 3' UTRs by statistical analysis of oligonucleotide distributions

Background

3' untranslated regions (3' UTRs) contain binding sites for many regulatory elements, and in particular for microRNAs (miRNAs). The importance of miRNA-mediated post-transcriptional regulation has become increasingly clear in the last few years.

Results

We propose two complementary approaches to the statistical analysis of oligonucleotide frequencies in mammalian 3' UTRs aimed at the identification of candidate binding sites for regulatory elements. The first method is based on the identification of sets of genes characterized by evolutionarily conserved overrepresentation of an oligonucleotide. The second method is based on the identification of oligonucleotides showing statistically significant strand asymmetry in their distribution in 3' UTRs.

Conclusion

Both methods are able to identify many previously known binding sites located in 3'UTRs, and in particular seed regions of known miRNAs. Many new candidates are proposed for experimental verification.

Substrate dependence of conformational changes in the RNA-binding domain of tristetraprolin assessed by fluorescence spectroscopy of tryptophan mutants

Association of tristetraprolin (TTP) with mRNAs containing selected AU-rich mRNA-destabilizing elements (AREs) initiates rapid cytoplasmic degradation of these transcripts. The RNA-binding activity of TTP is mediated by an internal tandem zinc finger domain that preferentially recognizes U-rich RNA ligands containing adjacent UUAU half-sites and is accompanied by conformational changes within the peptide. Here, we have used analogues of the TTP RNA-binding domain containing specific tryptophan substitutions to probe the Zn2+ and RNA substrate dependence of conformational events within individual zinc fingers. Fluorescence methods demonstrate that the N-terminal, but not C-terminal, zinc finger domain adopts a stably folded conformation in the presence of Zn2+. Denaturant titrations suggest that both the N- and C-terminal zinc fingers exhibit limited structural heterogeneity in the absence of RNA substrates, although this is more pronounced for the C-terminal finger. Binding to a cognate ARE substrate induced significant conformational changes within each zinc finger, which also included increased resistance to chemical denaturation. Studies with mutant ARE ligands revealed that a single UUAU half-site was sufficient to induce structural modulation of the N-terminal finger. However, RNA-dependent folding of the C-terminal zinc finger was only observed in the presence of tandem UUAU half-sites, suggesting that the conformation of this domain is linked not only to RNA substrate recognition but also to the ligand occupancy and/or conformational status of the N-terminal finger. Coupled with previous structural and thermodynamic analyses, these data provide a mechanistic framework for discrimination of RNA substrates involving ligand-dependent conformational adaptation of both zinc fingers within the TTP RNA-binding domain.

RNA expression in a cartilaginous fish cell line reveals ancient 3' noncoding regions highly conserved in vertebrates

We have established a cartilaginous fish cell line [Squalus acanthias embryo cell line (SAE)], a mesenchymal stem cell line derived from the embryo of an elasmobranch, the spiny dogfish shark S. acanthias. Elasmobranchs (sharks and rays) first appeared >400 million years ago, and existing species provide useful models for comparative vertebrate cell biology, physiology, and genomics. Comparative vertebrate genomics among evolutionarily distant organisms can provide sequence conservation information that facilitates identification of critical coding and noncoding regions. Although these genomic analyses are informative, experimental verification of functions of genomic sequences depends heavily on cell culture approaches. Using ESTs defining mRNAs derived from the SAE cell line, we identified lengthy and highly conserved gene-specific nucleotide sequences in the noncoding 3' UTRs of eight genes involved in the regulation of cell growth and proliferation. Conserved noncoding 3' mRNA regions detected by using the shark nucleotide sequences as a starting point were found in a range of other vertebrate orders, including bony fish, birds, amphibians, and mammals. Nucleotide identity of shark and human in these regions was remarkably well conserved. Our results indicate that highly conserved gene sequences dating from the appearance of jawed vertebrates and representing potential cis-regulatory elements can be identified through the use of cartilaginous fish as a baseline. Because the expression of genes in the SAE cell line was prerequisite for their identification, this cartilaginous fish culture system also provides a physiologically valid tool to test functional hypotheses on the role of these ancient conserved sequences in comparative cell biology.

A periodic pattern of mRNA secondary structure created by the genetic code

Single-stranded mRNA molecules form secondary structures through complementary self-interactions. Several hypotheses have been proposed on the relationship between the nucleotide sequence, encoded amino acid sequence and mRNA secondary structure. We performed the first transcriptome-wide in silico analysis of the human and mouse mRNA foldings and found a pronounced periodic pattern of nucleotide involvement in mRNA secondary structure. We show that this pattern is created by the structure of the genetic code, and the dinucleotide relative abundances are important for the maintenance of mRNA secondary structure. Although synonymous codon usage contributes to this pattern, it is intrinsic to the structure of the genetic code and manifests itself even in the absence of synonymous codon usage bias at the 4-fold degenerate sites. While all codon sites are important for the maintenance of mRNA secondary structure, degeneracy of the code allows regulation of stability and periodicity of mRNA secondary structure. We demonstrate that the third degenerate codon sites contribute most strongly to mRNA stability. These results convincingly support the hypothesis that redundancies in the genetic code allow transcripts to satisfy requirements for both protein structure and RNA structure. Our data show that selection may be operating on synonymous codons to maintain a more stable and ordered mRNA secondary structure, which is likely to be important for transcript stability and translation. We also demonstrate that functional domains of the mRNA [5′-untranslated region (5′-UTR), CDS and 3′-UTR] preferentially fold onto themselves, while the start codon and stop codon regions are characterized by relaxed secondary structures, which may facilitate initiation and termination of translation.

Evidence for conservation and selection of upstream open reading frames suggests probable encoding of bioactive peptides

Background

Approximately 40% of mammalian mRNA sequences contain AUG trinucleotides upstream of the main coding sequence, with a quarter of these AUGs demarcating open reading frames of 20 or more codons. In order to investigate whether these open reading frames may encode functional peptides, we have carried out a comparative genomic analysis of human and mouse mRNA 'untranslated regions' using sequences from the RefSeq mRNA sequence database.

Results

We have identified over 200 upstream open reading frames which are strongly conserved between the human and mouse genomes. Consensus sequences associated with efficient initiation of translation are overrepresented at the AUG trinucleotides of these upstream open reading frames, while comparative analysis of their DNA and putative peptide sequences shows evidence of purifying selection.

Conclusion

The occurrence of a large number of conserved upstream open reading frames, in association with features consistent with protein translation, strongly suggests evolutionary maintenance of the coding sequence and indicates probable functional expression of the peptides encoded within these upstream open reading frames.

Identification and characterization of cis-acting elements in the human and bovine PTH mRNA 3'-untranslated region

The human PTH mRNA 3'-UTR has a cis element homologous to the rat cis-acting instability element and a more proximal element identical to the single binding element identified in bovine PTH mRNA 3'-UTR. The function of the elements was shown in vitro.

INTRODUCTION

In the rat, Ca(2+) and phosphate regulate PTH mRNA stability by the interaction of trans-acting proteins with a defined cis-acting instability element in the distal region of the PTH mRNA 3'-untranslated region (UTR). This element has been characterized in the rat and is conserved in human, canine, feline, and murine 3'-UTRs but not in bovine and porcine 3'-UTRs.

MATERIALS AND METHODS

Parathyroid protein-binding assays to the PTH mRNA transcripts were performed. Functionality was studied in reporter genes that were transiently transfected into HEK293 cells.

RESULTS

Protein-RNA binding experiments identified an element in bovine PTH mRNA at the proximal end of the 3'-UTR that is different from the rat protein-binding element. The human 3'-UTR contains both elements, but only the distal element binds proteins. Functional studies with HEK293 cells transiently transfected with reporter genes containing the different elements and flanking nucleotides (nt) showed that the human distal element destabilized a reporter mRNA similar to the effect of this element in the rat. A reporter mRNA containing the single bovine PTH mRNA protein-binding element was also destabilized, and this was prevented by coexpression of AU-rich element binding factor 1 (AUF1).

CONCLUSION

Our results identify a new protein-binding element in the PTH mRNA 3'-UTR. In bovine PTH mRNA, it is the only element, and it is functional in destabilizing a reporter gene. It is also present in other species, including human PTH mRNA, where it is not functional, possibly because of differences in flanking sequences. The human PTH mRNA 3'-UTR distal element is highly homologous to the rat cis-acting instability element and destabilized a reporter gene, indicating its functionality. Therefore, different species have alternative cis-acting protein-binding elements that may determine the regulation of PTH mRNA stability in response to changes in serum calcium and phosphate.

Conservation of Bmp2 post-transcriptional regulatory mechanisms

Bone morphogenetic protein (BMP) orthologs from diverse species like flies and humans are functionally interchangeable and play key roles in fundamental processes such as dorso-ventral axis formation in metazoans. Because both transcriptional and post-transcriptional mechanisms play central roles in modulating developmental protein levels, we have analyzed the 3'-untranslated region (3'UTR) of the Bmp 2 gene. This 3'UTR is unusually long and is alternatively polyadenylated. Mouse, human, and dog mRNAs are 83-87% identical within this region. A 265-nucleotide sequence, conserved between mammals, birds, frogs, and fish, is present in Bmp2 but not Bmp4. The ability of AmphiBMP2/4, a chordate ortholog to Bmp2 and Bmp4, to align with this sequence suggests that its function may have been lost in Bmp4. Activation of reporter genes by the conserved region acts by a post-transcriptional mechanism. Mouse, human, chick, and zebrafish Bmp2 synthetic RNAs decay rapidly in extracts from cells not expressing Bmp2. In contrast, these RNAs are relatively stable in extracts from Bmp2-expressing cells. Thus, Bmp2 RNA half-lives in vitro correlate with natural Bmp2 mRNA levels. The fact that non-murine RNAs interact appropriately with the mouse decay machinery suggests that the function of these cis-regulatory regions has been conserved for 450 million years since the fish and tetrapod lineages diverged. Overall, our results suggest that the Bmp2 3'UTR contains essential regulatory elements that act post-transcriptionally.

The mammalian transcriptome and the function of non-coding DNA sequences

Many non-coding sequences transcribed from the mammalian genome are proving to have important regulatory roles, but the functions of the majority remain mysterious.

For decades, researchers have focused most of their attention on protein-coding genes and proteins. With the completion of the human and mouse genomes and the accumulation of data on the mammalian transcriptome, the focus now shifts to non-coding DNA sequences, RNA-coding genes and their transcripts. Many non-coding transcribed sequences are proving to have important regulatory roles, but the functions of the majority remain mysterious.

Comparative analysis of orthologous eukaryotic mRNAs: potential hidden functional signals

Sequencing of multiple, nearly complete eukaryotic genomes creates opportunities for detecting previously unnoticed, subtle functional signals in non-coding regions. A genome-wide comparative analysis of orthologous sets of mammalian and yeast mRNAs revealed distinct patterns of evolutionary conservation at the boundaries of the untranslated regions (UTRs) and the coding region (CDS). Elevated sequence conservation was detected in approximately 30 nt regions around the start codon. There seems to be a complementary relationship between sequence conservation in the approximately 30 nt regions of the 5'-UTR immediately upstream of the start codon and that in the synonymous positions of the 5'-terminal 30 nt of the CDS: in mammalian mRNAs, the 5'-UTR shows a greater conservation than the CDS, whereas the opposite trend holds for yeast mRNAs. Unexpectedly, a approximately 30 nt region downstream of the stop codon shows a substantially lower level of sequence conservation than the downstream portions of the 3'-UTRs. However, the sequence in this poorly conserved 30 nt portion of the 3'-UTR is non-random in that it has a higher GC content than the rest of the UTR. It is hypothesized that the elevated sequence conservation in the region immediately upstream of the start codon is related to the requirement for initiation factor binding during pre-initiation ribosomal scanning. In contrast, the poorly conserved region downstream of the stop codon could be involved in the post- termination scanning and dissociation of the ribosomes from the mRNA, which requires only the mRNA-ribosome interaction. Additionally, it was found that the choice of the stop codon in mammals, but not in yeasts, and the context in the immediate vicinity of the stop codons in both mammals and yeasts are subject to strong selection. Thus, genome-wide analysis of orthologous gene sets allows detection of previously unrecognized patterns of sequence conservation, which are likely to reflect hidden functional signals, such as ribosomal filters that could regulate translation by modulating the interaction between the mRNA and ribosomes.