{"id":787,"date":"2023-03-16T13:57:19","date_gmt":"2023-03-16T06:57:19","guid":{"rendered":"https:\/\/conf.icgbio.ru\/bgrs98\/?page_id=787"},"modified":"2023-09-04T15:47:36","modified_gmt":"2023-09-04T08:47:36","slug":"060_eukaryotic-mrnas-encoding-abundant-and-scarce-proteins-are-dissimilar-in-many-structural-features-of-5-untranslated-leaders","status":"publish","type":"page","link":"https:\/\/conf.icgbio.ru\/bgrs98\/abstracts\/abstract-list\/060_eukaryotic-mrnas-encoding-abundant-and-scarce-proteins-are-dissimilar-in-many-structural-features-of-5-untranslated-leaders\/","title":{"rendered":"EUKARYOTIC MRNAS ENCODING ABUNDANT AND SCARCE PROTEINS ARE DISSIMILAR IN MANY STRUCTURAL FEATURES OF 5\u2019-UNTRANSLATED LEADERS"},"content":{"rendered":"

KOCHETOV A.V.<\/a>+<\/sup>,\u00a0PONOMARENKO M.P.<\/a>,\u00a0VOROBIEV D.G.<\/a>,\u00a0FROLOV A.S.<\/a>,\u00a0KISSELEV L.L.<\/a>1<\/sup>,\u00a0KOLCHANOV N.A.<\/a><\/p>\n

Institute of Cytology and Genetics, Siberian Branch of the Russian Academy of Sciences, 10 Lavrentiev Ave., Novosibirsk, 630090, Russia; e-mail:ak@bionet.nsc.ru<\/a><\/p>\n

1<\/sup>Engelhardt Institute of Molecular Biology, Moscow, 117984 Russia, e-mail:kissel@imb.imb.msk.ru<\/p>\n

+<\/sup>Corresponding author<\/p>\n

Keywords<\/a>: eukarioyic mRNA, structural features, 5′-untraslated leaders, low-expression mRNA, high-expression mRNA, computer program<\/p>\n

 <\/p>\n

Abstract<\/b><\/p>\n

Structural characteristics of 5\u2019-untranslated regions of eukaryotic mRNAs encoding highly and low abundant proteins were compared. It was found that the leader sequences of the low-expression mRNAs are longer, their guanine-plus-cytosine content is higher, they form more stable secondary structures, exhibit weaker context of the translation initiation codons and more frequent occurrence of AUGs compared with the 5\u2019UTRs of high-expression mRNAs. These structural features of the low- and high-mRNAs may contribute to their different translational efficiency and abundance of their protein products. Computer program capable of efficient discrimination between the high- and low-expression genes of dicot plants basing on analysis of their mRNA 5\u2019UTR features was designed.<\/p>\n

1. Introduction<\/b><\/p>\n

The level of gene expression, evaluated as the amount of a protein product, varies considerably from a few molecules per cell (e.g., some transcription factors) to 1-2% of the total cellular protein content (e.g., translation elongation factor I\"\"\u00a0and actin). Until recently, the major contribution to the regulation of gene expression was thought to occur at the transcriptional level. However, other processes, such as maturation of mRNA, rates of its translation and degradation, can also be regulated and thereby affect the level of gene expression. It is known that eukaryotic mRNAs considerably differ in their translation efficiency thought to be attributed to varying efficiency of translation initiation1<\/sup>, since the contextual and structural features of the 5\u2019-untranslated region (5\u2019-UTR) significantly influence translation initiation. According to the scanning model2<\/sup>, 40S ribosomal subunits recognize the cap at the 5\u2019 end of mRNA, bind to it, and linearly scan the mRNA in the 5\u2019=> 3\u2019 direction in search for the nearest AUG codon; once the initiation codon is reached, binding of the 40S ribosomal subunit follows, and in most cases translation proceeds to the elongation step. The productive recognition of the AUG codon as an initiator depends on its nucleotide context. Adenine at position -3 relative to AUG and, to less extent, guanine at position +4 provide the optimal context for translation initiation, while pyrimidine at position -3 decreases its efficiency.<\/p>\n

There are some 5\u2019UTR features influencing the mRNA translation efficiency. AUGs occurring within 5\u2019UTR can also interact with 40S ribosomal subunits and decrease the translation rate from the true start site2,3<\/sup>. The secondary structure of 5\u2019-UTRs negatively affects the migration of 40S ribosomal subunits along mRNA\u00a02,3<\/sup>, presumably because it takes time to overcome the hindrance, thereby decreasing the translation rate. Negative influence of a hairpin in the\u00a0in vivo<\/i>\u00a0translation of the mammalian mRNA depends on the hairpin stability and its location in the leader2<\/sup>.<\/p>\n

The mRNAs of the eukaryotic genes contrasting in their expression levels are likely to differ in their contextual features and structural organization. To test this assumption, we compared the mRNA features of several groups of the housekeeping genes, highly expressed in eukaryotic cells (H-mRNAs), and regulatory genes, whose expression is low and under stringent control (L-mRNAs).<\/p>\n

2. Materials and methods<\/b><\/p>\n

mRNA sequences were taken from EMBL database, release 49. The similarity scores for all the 5\u2019UTR sequence pairs were determined.\u00a0<\/b>If two sequences in the database shared a similarity greater than 90%, the shorter sequence was deleted from the set.<\/p>\n

3. Results<\/b><\/p>\n

3.1 Selection of mRNAs for proteins of high and low abundancy<\/i><\/b><\/p>\n

In this study, we analyzed sequences of two distinct taxa: mammals and higher plants. Mammalian high-expression mRNAs (H-mRNA set) was consisted of mRNAs encoding highly abundant eukaryotic proteins of the following families: translation elongation factor 1 alpha (EF1\"\") and ribosomal proteins, actins, tubulins, 70-kDa heat shock proteins, myosins, and histones. Plant H-mRNA set was consisted of mRNAs of photosynthesis-related genes (RbcS, Cab), actins, histones, and stress-responsive genes (HSPs, ADH, ALD, GST). Totally, 114 5\u2019UTRs of mammalian mRNAs and 146 5\u2019UTRs of plant mRNAs (only dicot plants were used) were extractedb (Table 1).<\/p>\n

The L-mRNAs, whose protein products are present in the cells in small amounts, was represented by 198 sequences of mammalian and 247 sequences of plant 5\u2019UTRs, including the genes of regulatory proteins (growth factors, receptors, transcription factors extracted from the TRANSFAC database, and regulatory protein molecules of other types). The expression of these genes is known to be under stringent control not only at the transcriptional level, but also through a decrease in stability of mRNA4<\/sup>\u00a0and proteins5<\/sup>.<\/p>\n

3.2. Length of the 5\u2019UTRs of H- and L-mRNAs<\/b><\/i><\/p>\n

The mean length of the 5\u2019UTR of L-mRNAs (157 nt for the plant and 203 nt for the mammalian 5\u2019UTR sets) considerably exceeded that of the H-mRNAs (85 nt for the plant and 80 nt for the mammalian 5\u2019UTR sets). 5\u2019UTRs of mammalian mRNAs were not longer than 100 nucleotides in 70% of H-mRNAs and in 30% of L-mRNAs. According to Kolmogorov-Smirnov’s test, the difference between the length distributions for these two 5\u2019UTR sets is significant (P<0.001) for both mammalian and plant samples. Consequently, the leader length clearly distinguishes the H- and L-mRNAs.<\/p>\n

3.3. Nucleotide composition of the H- and L-mRNAs<\/b><\/i><\/p>\n

G and C significantly contribute to the stability of the RNA secondary structure. For this reason, we analyzed the G+C content in the mRNA sets considering heterogeneous taxonomic composition of the sets.<\/p>\n

 <\/p>\n

Table 1.\u00a0<\/b>G+C content (%, mean \"\"s.d.) in 5\u2019UTRs of the H- and L-mRNAs from distinct taxonomic groups<\/p>\n\n\n\n\n\n\n
<\/td>\nG+C content in the 5\u2019UTRs<\/b><\/td>\n<\/tr>\n
Taxa<\/b><\/td>\n\n

H-mRNAs<\/p>\n<\/td>\n

\n

L-mRNAs<\/p>\n<\/td>\n<\/tr>\n

Mammalia<\/td>\n\n

114<\/p>\n<\/td>\n

\n

56.0\"\"12.0<\/p>\n<\/td>\n

\n

198<\/p>\n<\/td>\n

\n

62.8\"\"12.5<\/p>\n<\/td>\n<\/tr>\n

Plant (dicot)<\/td>\n\n

146<\/p>\n<\/td>\n

\n

36.8\"\"9.3<\/p>\n<\/td>\n

\n

247<\/p>\n<\/td>\n

\n

37.7\"\"7.1<\/p>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

 <\/p>\n

It was demonstrated that G+C content of the L-mRNA 5\u2019UTRs is higher than that of the H-mRNAs of the same taxon. The difference in the G+C content of the H- and L- mammalian mRNA 5\u2019UTRs is significant (P<0.001), as shown by K.-S. test.<\/p>\n

At the next step, we analyzed the mRNA 5\u2019UTRs for the ratios of the frequencies of the complementary nucleotides involved in formation of the secondary structure (G\/C and A\/U). The data are listed in Table 2.<\/p>\n

 <\/p>\n

Table 2.G\/C (A\/U) ratio (%, mean \"\"s.d.) in the 5\u2019UTRs of the H- and L-mRNAs from distinct taxonomic groups<\/p>\n\n\n\n\n\n\n
 <\/p>\n

Taxon<\/b><\/p>\n<\/td>\n

 <\/p>\n

Mammalian mRNA 5\u2019UTRs<\/b><\/p>\n<\/td>\n

 <\/p>\n

Dicot plant mRNA 5\u2019UTRs<\/b><\/p>\n<\/td>\n<\/tr>\n

 <\/p>\n

Expression level<\/i><\/p>\n<\/td>\n

 <\/p>\n

H-mRNAs<\/i><\/p>\n<\/td>\n

 <\/p>\n

L-mRNAs<\/i><\/p>\n<\/td>\n

 <\/p>\n

H-mRNAs<\/i><\/p>\n<\/td>\n

 <\/p>\n

L-mRNAs<\/i><\/p>\n<\/td>\n<\/tr>\n

0.8<G\/C<1.2<\/td>\n\n

26.9<\/p>\n<\/td>\n

\n

34.3<\/p>\n<\/td>\n

\n

17.8<\/p>\n<\/td>\n

\n

23.5<\/p>\n<\/td>\n<\/tr>\n

0.8<A\/U<1.2<\/td>\n\n

20.8<\/p>\n<\/td>\n

\n

38.9<\/p>\n<\/td>\n

\n

21.2<\/p>\n<\/td>\n

\n

26.3<\/p>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

 <\/p>\n

It was shown that the contents of the complementary nucleotides were considerably more asymmetric in the H-mRNA 5\u2019UTRs.<\/p>\n

3.4. Context of the initiator AUG codon<\/b><\/i><\/p>\n

To analyze the context of translation initiation site, the mRNA fragments (from -10 nt to -1 nt upstream of AUG) were aligned (data not shown). The contexts of the AUG codons of the mammalian H- and L- mRNAs differed considerably. Pyrimidines at position -3 prior to AUG occurred in 23.2% of the L-mRNAs and in 4.35% of the H-mRNAs. Adenine with the strongest positive effect on the context \u201cstrength\u201d occurred most frequently at position -3 in both the H- and L-mRNAs. However, A in this position occurred 1.5 times less frequently in the L-mRNAs (40.4%) than in the H-mRNAs (59.1%). Context of the translation initiation site in the dicot H-mRNAs is also more optimal: 45.9% of high- versus 30% of the low-expression mRNAs contained the optimal combination A-3<\/sup>\/G+4<\/sup>. This difference indicates that the context of the translation initiation codon is more often non-optimal in the L-mRNAs.<\/p>\n

3.5. AUG codons in 5\u2019UTRs of the H- and L-mRNAs<\/b><\/i><\/p>\n

The AUG codon frequencies in 5\u2019UTRs were compared for the sets of mammalian H- and L-mRNAs and found to differ significantly (P<0.001). AUG codons were recorded in 40 of 327 H-mRNA leaders and in 112 of 287 L-mRNA leaders. Hence, the proportion of the AUG-containing 5\u2019UTRs in the L-mRNAs was 2.5-fold higher than in the H-mRNAs. In 16 of 40 AUG-containing 5\u2019UTRs of the H-mRNAs, AUG codons are in a non-optimal context, i.e. with pyrimidine in position -3 and non-G in position +4. As to the L-mRNAs, 27 of 112 mRNA 5\u2019UTRs contained only non-optimal AUG codons. Similarly, AUG-containing plant 5\u2019UTRs constituted 45.3% of the L-mRNA and 8.8% of H-mRNA sets.<\/p>\n

It is possible that the increased content of AUGs in the 5\u2019UTRs of L-mRNAs is due to their greater length. To verify this, we calculated the observed AUG frequency in 5\u2019UTRs of the mammalian H- and L-mRNAs normalized to the corresponding 5\u2019UTR lengths. Expected AUG frequencies were calculated according to equation: Paug=Pa*Pu*Pg, where Px is the expected content of the nucleotide X in the 5\u2019UTRs of H- and L-mRNAs of the corresponding taxon. It appeared that the ratio of Obs to Exp AUG frequencies in the 5\u2019UTRs of mammalian (0.514) and plant (0.503) L-mRNAs was considerably higher than for the 5\u2019UTRs of H-mRNAs (0.326 and 0.343, respectively). Therefore, a higher AUG content in the 5\u2019UTRs in L-mRNAs is not simply due to their greater lengths. For comparison: the ratios of the observed to expected AUG frequencies in 3\u2019UTRs of mammalian H- and L-mRNAs were virtually identical (0.93 and 0.94, respectively).<\/p>\n

Thus, it was demonstrated that the difference in the expression levels of eukaryotic genes is often accompanied by a significant difference in the structural features of their mRNAs. The leaders of the H-mRNAs are poorer in G+C content, the ratio of their complementary nucleotides is much more asymmetric, their secondary structure is less stable, and they contain less AUGs. More initiation codons in the H-mRNAs are in the optimal or suboptimal contexts. In contrast, the features 5\u2019UTRs of the mRNAs encoding regulatory proteins decrease considerably the translational activity of the templates.<\/p>\n

The structural features of H-mRNAs facilitate highly efficient initiation of protein synthesis; the optimization of mRNA\u00a0<\/i>translation efficiency is perhaps essential to reach the fast and high-level protein expression. On the other hand, the regulation of transcription seems insufficient to control the expression level of low-abundant proteins, and the low translation efficiency of these mRNAs may be of general importance for the stringent control of regulatory gene expression.<\/p>\n

3.5. The computer system for prediction of the mRNA translational activity.<\/b><\/i><\/p>\n

Considerable difference between the structural characteristics of mRNA 5\u2019UTRs of the high- and low-expression genes were used to design the program for the discrimination between the high- and low-expression mRNAs of dicot plants (http:\/\/wwwmgs.bionet.nsc.ru\/Programs\/acts2\/dicot.htm<\/a>). Twelve criteria (Fi) were generated:<\/p>\n

 <\/p>\n

Expert Weights (0-10 are valid; 5 automatically)
\n<\/b>1. Translation INCREASES with DECREASING the Leader length
\n2. Translation INCREASES with DECREASING the G+C content
\n3. Translation INCREASES with INCREASING the G\/C-disbalance
\n4. Translation INCREASES with DECREASING the alt-AUG content
\n5. Translation INCREASES with DECREASING the framed AUG content
\n6. Translation INCREASES depending on the “-3 position” rule
\n7. Translation INCREASES with DECREASING the AUG inside the leader
\n8. Translation INCREASES with INCREASING the [C] content
\n9. Translation INCREASES with INCREASING the [YM] content
\n10. Translation INCREASES with INCREASING the[CnY] content
\n11. Comparison with the weight matrix for nucl. content in 5\u2019UTRs of H- mRNAs
\n12. Comparison with complex (high to low) weight matrix for nucl. content<\/i><\/p>\n

Fig.1. List of the characteristics of a mRNA 5\u2019UTR that were used for prediction of the gene expression level. Parameters 8-12 were determined for the (-35;-1) 5\u2019UTR fragments.<\/p>\n

 <\/p>\n

The program determines the 5\u2019UTR characteristics and evaluates the translational activity of the mRNA according to formula:\u00a0F(seq)<\/i> = \"\". User may change the list of criteria and their weights. The discrimination between the control samples of the high- and low-expression mRNAs of dicot plants showed a good results: 84% of the high- and 76% of the low-expression mRNAs were classified correctly (Fig.1).<\/p>\n

 <\/p>\n

\"\"<\/a><\/p>\n

Fig. 2. The control results obtained using independent data. Broken line is the predicting rule \u011eif F>0 then High expression, otherwise Low expression\u011f.<\/p>\n

 <\/p>\n

Acknowledgments<\/b><\/p>\n

A.K. was supported by the SD RAS grant for young scientists.<\/p>\n

 <\/p>\n

References<\/b><\/p>\n

    \n
  1. B.K. Ray, T.G. Brendler, S. Adya, S. Daniels-McQeen, J.K. Miller, J.W.B. Hershey, J.A. Grifo, W.C. Merrick, and R.E. Thach, \u201cRole of mRNA competition in regulating translation: further characterization of mRNA discriminatory initiation factors\u201d Proc. Natl. Acad. Sci. USA\u00a080<\/b>, 663-667 (1983).<\/li>\n
  2. M. Kozak, \u201cStructural features in eukaryotic mRNAs that modulate the initiation of translation\u201d J. Biol. Chem.\u00a0266,<\/b>\u00a019867-19870 (1991).<\/li>\n
  3. V. Pain, \u201cInitiation of protein synthesis in eukaryotic cells\u201d Eur. J. Biochem.\u00a0236,<\/b>\u00a0747-771 (1996).<\/li>\n
  4. M.E. Greenberg and J.G. Belasco, \u201cControl of the decay of labile protooncogene and cytokine mRNAs, \/\/In: Control of messenger RNA stability\u201d (Eds.: Belasco, M.E., Brawerman, G. P. 1993. Academy Press, San Diego, Calif. 199-218).<\/li>\n
  5. H.L. Pahl and P.A. Baeuerle, \u201cControl of gene expression by proteolysis\u201d Curr. Opin. Cell Biol.\u00a08,<\/b>\u00a0340-347 (1996).<\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"

    KOCHETOV A.V.+,\u00a0PONOMARENKO M.P.,\u00a0VOROBIEV D.G.,\u00a0FROLOV A.S.,\u00a0KISSELEV L.L.1,\u00a0KOLCHANOV N.A. Institute of Cytology and Genetics, Siberian Branch of the Russian Academy of Sciences, 10 Lavrentiev Ave., Novosibirsk, 630090, Russia; e-mail:ak@bionet.nsc.ru 1Engelhardt Institute of Molecular Biology, Moscow, 117984 Russia, e-mail:kissel@imb.imb.msk.ru +Corresponding author Keywords: eukarioyic … Continue reading →<\/span><\/a><\/p>\n","protected":false},"author":13,"featured_media":0,"parent":97,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":[],"_links":{"self":[{"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages\/787"}],"collection":[{"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/users\/13"}],"replies":[{"embeddable":true,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/comments?post=787"}],"version-history":[{"count":5,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages\/787\/revisions"}],"predecessor-version":[{"id":1511,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages\/787\/revisions\/1511"}],"up":[{"embeddable":true,"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/pages\/97"}],"wp:attachment":[{"href":"https:\/\/conf.icgbio.ru\/bgrs98\/wp-json\/wp\/v2\/media?parent=787"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}