Critique the aging paper in the upload file.

Naked mole-rat has increased translational fidelity compared with the mouse, as well as a unique 28S ribosomal RNA cleavage Jorge Azpuruaa,1, Zhonghe Kea,1, Iris X. Chena, Quanwei Zhangb, Dmitri N. Ermolenkoc, Zhengdong D. Zhangb, Vera Gorbunovaa,2, and Andrei Seluanova,2

aDepartment of Biology, University of Rochester, Rochester, NY 14627; bDepartment of Genetics, Albert Einstein College of Medicine, New York, NY 10461; and cDepartment of Biochemistry and Biophysics, University of Rochester, Rochester, NY 14642

Edited* by Eviatar Nevo, Institute of Evolution, Haifa, Israel, and approved September 3, 2013 (received for review July 18, 2013)

The naked mole-rat (Heterocephalus glaber) is a subterranean euso- cial rodent with a markedly long lifespan and resistance to tumor- igenesis. Multiple data implicate modulation of protein translation in longevity. Here we report that 28S ribosomal RNA (rRNA) of the naked mole-rat is processed into two smaller fragments of un- equal size. The two breakpoints are located in the 28S rRNA di- vergent region 6 and excise a fragment of 263 nt. The excised fragment is unique to the naked mole-rat rRNA and does not show homology to other genomic regions. Because this hidden break site could alter ribosome structure, we investigated whether trans- lation rate and amino acid incorporation fidelity were altered. We report that naked mole-rat fibroblasts have significantly increased translational fidelity despite having comparable translation rates with mouse fibroblasts. Although we cannot directly test whether the unique 28S rRNA structure contributes to the increased fidelity of translation, we speculate that it may change the folding or dynamics of the large ribosomal subunit, altering the rate of GTP hydrolysis and/or interaction of the large subunit with tRNA dur- ing accommodation, thus affecting the fidelity of protein synthe- sis. In summary, our results show that naked mole-rat cells produce fewer aberrant proteins, supporting the hypothesis that the more stable proteome of the naked mole-rat contributes to its longevity.

aging | NMR

Naked mole-rats (NMRs; Heterocephalus glaber) are extremelongevity outliers compared with other rodents of their size. Although only slightly larger than a laboratory mouse, the NMR lives an order of magnitude longer, with a maximum recorded lifespan of more than 30 y for both reproductive and non- reproductive castes (1). It was also noted that NMRs are markedly resistant to the increased frailty that accompanies ag- ing in most metazoans. Furthermore, despite observations of many aged individuals, no neoplasia has been reported in the NMR (2, 3). NMR cells have been experimentally shown to be much more

resistant to oncogenic transformation than mouse cells (4, 5). The mechanism of tumor resistance is mediated by hypersensi- tivity of NMR cells to contact inhibition triggered by hyaluronan of extremely high molecular mass (5–7). Although these studies explain why the animals are tumor-resistant, they do not explain their overall longevity and their ability to maintain somatic in- tegrity until extreme age. Previous studies have indicated that the NMR is able to tolerate high intracellular levels of oxidized proteins and oxidized lipids (8–10). Furthermore, the NMR is notably able to maintain proteomic integrity (as evidenced by lower protein unfolding after treatment with urea) until very late in its lifespan without an age-dependent increase in protein ubiquitination, such as observed in short-lived mice (11). Multiple signaling pathways implicated in aging and cancer

converge on factors that modify the rate of translation initiation and elongation (12–15). One of the main functions of the mTOR

pathway is the promotion of protein translation in response to nutrients. Modulating the mTOR pathway genetically or phar- macologically, using the inhibitor rapamycin, has been shown to extend lifespan in yeast, Caenorhabditis elegans, Drosophila, and mice (14). The exact mechanism by which modulating protein translation via inhibition of mTOR extends lifespan is unclear. Two major explanations that were proposed are global reduction in mRNA translation, leading to better maintenance of protein homeostasis, or differential translation of specific mRNAs ben- eficial for longevity and stress resistance (reviewed in ref. 13). Proper protein folding and stability are heavily implicated in

aging. Overexpression of heat shock transcription factors and chaperones that aid protein folding can extend lifespan in Dro- sophila melanogaster (16) and C. elegans (17, 18). In bacteria, aging occurs by asymmetrical segregation of protein aggregates (19), and a similar mechanism has been described in yeast (20). The formation of aggregates as a result of mistranslated poly- peptides has been implicated in several aging-related diseases such as amyotrophic lateral sclerosis, Huntington’s disease, and Alzheimer’s disease (21). A major factor determining proteome quality is translational fidelity of the ribosome. Severe disruptions of translational fidelity have been shown to cause accumulation of aggregates and pathological neurodegeneration (22). The connection between translation and aging led us to ex-

amine NMR ribosomes for longevity-promoting phenotypes such as altered translation rate and translational fidelity. We describe a unique ribosomal DNA insertion, which is processed after rRNA transcription and breaks the NMR 28S rRNA into two distinct fragments. We assayed ribosomal translation rate (com- bined initiation and elongation) by measuring bulk incorporation

Significance

Molecular mechanisms responsible for differences in longevity between animal species are largely unknown. Here we show that the longest-lived rodent, the naked mole-rat, has more accurate protein translation than the mouse. Furthermore, we show that the naked mole-rat has a unique fragmented ribo- somal RNA structure. Such cleaved ribosomal RNA has been reported for only one other species of mammal. This article suggests the importance of protein translation in aging and provides insight into the mechanisms of longevity.

Author contributions: J.A., Z.K., Z.D.Z., V.G., and A.S. designed research; J.A., Z.K., I.X.C., Q.Z., and Z.D.Z. performed research; J.A., Z.K., Q.Z., D.N.E., Z.D.Z., V.G., and A.S. analyzed data; and J.A., Z.K., Z.D.Z., V.G., and A.S. wrote the paper.

The authors declare no conflict of interest.

*This Direct Submission article had a prearranged editor. 1J.A. and Z.K. contributed equally to this work. 2To whom correspondence may be addressed. E-mail: [email protected] or [email protected].

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. 1073/pnas.1313473110/-/DCSupplemental.

17350–17355 | PNAS | October 22, 2013 | vol. 110 | no. 43 www.pnas.org/cgi/doi/10.1073/pnas.1313473110

D ow

nl oa

de d

fr om

h tt

ps :/

/w w

w .p

na s.

or g

by 7

3. 20

3. 24

5. 0

on A

pr il

1 2,

2 02

2 fr

om I

P a

dd re

ss 7

3. 20

3. 24

5. 0.

of [35S]-cystein/methionine into nascent polypeptide chains. We found that NMR and mouse skin fibroblasts have roughly equal rates of translation. By using firefly luciferase reporters with a variety of mutations, including amino acid substitutions or stop codons that required misincorporation or frameshift to occur to reactivate the luciferase, we found that primary NMR fibroblasts have a fourfold lower rate of amino acid misincorporation than primary mouse fibroblasts. These results suggest that NMR cells produce fewer aberrant proteins and support the hypothesis that the NMR’s stable proteome is a contributor to its longevity.

Results NMR 28S rRNA Has an Atypical Pattern. During routine extraction of RNA, we noticed that the 28S ribosomal RNA of the NMR was processed into two fragments of unequal size. NMR RNA extracts from several tissues (liver, lung, brain, spleen, kidney, and testis) exhibited an unusual rRNA band pattern, whereas mouse RNA extracts did not (Fig. 1A). We observed this pattern after using two different methods of RNA extraction in extracts prepared from tissues and cultured cell lines (Fig. S1A). We also extracted RNA from other nonmodel rodents (Fig. S1B) and found the break to be unique for the NMR. To confirm that the break was not a result of RNase

contamination, we combined NMR cells with mouse cells and simultaneously extracted RNA after homogenizing them. When independently extracted, the fibroblast RNA showed three bands and two bands for NMR and mouse, respectively. In the mixed sample, all bands could be seen, indicating that no nuclease contamination was cleaving the mouse rRNA (Fig. 1B). This NMR 28S rDNA was undergoing specific processing in vivo.

Mapping of the 28S rRNA Cleavage Sites. To determine the site and sequence of the NMR 28S rRNA breakage point, we first designed primers to conserved 28S rDNA sequence from the mouse genome. We estimated the location of the breakpoint region by noting that the sizes of the two fragments (2.5 and

3 kb) were only slightly different, meaning the break would be within ∼1 kb of the center. This allowed us to amplify a central 1-kb fragment of genomic NMR 28S rDNA. This fragment included 28S divergent region 6 (D6) that exhibited a 118-nt unique insertion compared with the mouse sequence flanked by regions of weak homology with the mouse D6 (Fig. 2). Using the 1-kb sequence and a publically available NMR coding DNA sequence, we designed primers to perform 5′ RACE and 3′ RACE. After aligning with the NMR 28S rDNA sequence, we identified the two cleavage sites and found that a 263-nt fragment flanked by 5-nt direct repeats 5′-CGGAC-3′ was cleaved out in the 28S rRNA, leaving behind the fragmented 28S rRNA. The 263-nt fragment included the 118-nt insertion that was not homologous to the mouse sequence (Fig. 2). The ribosomal DNA loci are repeated thousands of times, and from our RNA extractions, most of the products seemed to be cleaved. We used a bioinformatics approach to test whether the identified cut sites were conserved among all of the 28S rDNA loci in the NMR. From whole-genome shotgun reads from two different NMR genome sequences [Broad Institute (http://www.ncbi.nlm. nih.gov/nuccore/AHKG00000000) and Kim et al. (23)], the cut sites were found to be conserved among all of the rDNA loci (Fig. S2 and Table S1), which confirmed that the majority of 28S rRNAs are cleaved in the NMR.

NMR Translation Rate Does Not Differ from the Mouse. We hy- pothesized that the fragmented 28S rRNA might cause alternate assembly or structure of the NMR ribosome, with detectable consequences for protein synthesis rate or fidelity. To investigate whether the NMR ribosome had a different translation rate than

A

Mouse NMR Mix B

Naked Mole-Rat

28S

18S

Li ve

r Lu

ng

B ra

in

Sp le

en

K id

ne y

Te st

es

Li ve

r

Mouse

18S

28S cleaved

Fig. 1. NMR rRNA shows an unusual fragmented pattern. (A) 28R rRNA from multiple organs of the NMR appears as two fragments. Arrows indicate 28S fragments and an 18S band. (B) The fragmentation of NMR 28S rRNA is not caused by degradation during RNA extraction. Either NMR and mouse RNA were extracted from fibroblasts independently or the cells from both species were mixed together before extraction. RNA was separated on 1% denaturing agarose gel.

D6 region

Mouse 28S rRNA

Naked mole-rat 28S rRNA

Ctenomys 28S rRNA

A

B

154 bp

No break

Mouse

NMR

Mouse

NMR

Mouse

NMR

Mouse

NMR

296 bp

263 bp

255 bp

1864 1931

1965

2036

1932

1967

Fig. 2. Structure of the NMR 28S rRNA and localization of break sites. (A) Overall structures of the 28S rRNA in mouse, NMR, and Ctenomys. The mouse 28S rRNA has no cleavage sites. The D6 region of the NMR 28S rRNA contains a unique 118-nt insertion (shaded box). The two cleavage sites then excise a 263-nt fragment. Ctenomys is the only other known mammal with fragmented rRNA. The D6 region of the Ctenomys 28S rRNA contains a unique 106-bp insertion (shaded box) with the cleavage site located within this unique sequence. Thick black lines represent the 28S rRNA, the white box represents the D6 region, the arrows denote the cleavage sites, and the size of the D6 region and the distance between the two cleavage sites are indicated. (B) The alignment of the NMR and mouse 28S rDNA sequences. The box highlights the D6 region, the arrows indicate the cleavage sites, and hyphens indicate the gap.

Azpurua et al. PNAS | October 22, 2013 | vol. 110 | no. 43 | 17351

C EL L B IO LO

G Y

D ow

nl oa

de d

fr om

h tt

ps :/

/w w

w .p

na s.

or g

by 7

3. 20

3. 24

5. 0

on A

pr il

1 2,

2 02

2 fr

om I

P a

dd re

ss 7

3. 20

3. 24

5. 0.

the mouse ribosome, we implemented an assay based on mea- suring [35S]-cystein/methionine incorporation into nascent pro- tein. As the physiological temperature of the NMR is 32 °C, rather than 37 °C, and NMR fibroblasts are usually cultured at 32 °C, we tested the translation rate at both temperatures for both mouse and NMR cells. Mouse and NMR fibroblasts were given a pulse of [35S]-cystein/methionine for 30 min, cultured at 32 °C and 37 °C, total cell proteins were extracted, and the amount of radioactive isotope incorporation was measured by a scintillation counter. The translation rate was found to be altered twofold by shifting

the temperature by 5° in both species (Student t test, NMR, P = 0.02; mouse, P = 0.0027), but no significant differences were seen between the NMR and mouse cells when they were temperature- matched (Fig. 3A). To confirm that the protein quantification was accurate and that mouse and NMR fibroblasts have similar amounts of ribosomes, we performed Western blots on ribosomal protein S6 (RPS6). No large differences were seen in the overall amount of ribosomal protein or loading controls from mouse and NMR fibroblasts (Fig. 3B). Taken together, these results show that the NMR translation rate is equivalent to that of the mouse when they are temperature-matched, and that the cells express similar levels of a ribosomal protein.

NMR Has Higher Translational Fidelity than the Mouse. We next ex- amined whether altered ribosomal structure might affect the overall translation fidelity of proteins. To assay mouse and NMR cells for translational fidelity, we constructed a set of firefly lu- ciferase mutant reporters. In these reporters, a key amino acid in the luciferase protein (lysine 529) was mutated to completely destroy luciferase activity (24), and only if the wrong amino acid were incorporated into the nascent polypeptide would enzymatic activity be restored. The reporters included all triplet positions,

premature stop (third codon mutation at the beginning of the protein), and both positive and negative frameshifts by nucleo- tide insertion and deletion, respectively. The luciferase assay indicated that NMR cells had higher

protein translation fidelity than mouse cells. This was particularly noticeable when mouse and NMR cells were grown at 32 °C, as the NMR cells had ∼10-fold lower misincorporation frequency for the lysine triplet mutants, and these differences were highly statistically significant for all reporters (Student t test, P < 0.005) (Fig. 4A). At 37 °C, NMR cells had approximately a fourfold lower misincorporation; this was also statistically significant for all reporters (Student t test, P < 0.05) (Fig. 4B). Notably, the NMR cells were able to maintain a high trans-

lational fidelity even when forced to grow more quickly. We in- tegrated simian virus 40 (SV40) Large T-antigen (LargeT) and human Ras Glycine 12 to Valine (H-RasV12) mutant protein into wild-type NMR fibroblasts to induce more rapid growth via constant Ras signaling and abrogation of the p53/retinoblastoma protein pRB tumor suppressor pathways to prevent oncogene- induced senescence. The fidelity assays were then repeated on these cells grown at 32 °C and 37 °C. The translation error rate was not significantly different from in the primary NMR fibroblasts at either temperature (Fig. 4 A and B), indicating that even when NMR cells proliferate at a similar rate as mouse cells, they still maintained higher translational fidelity. This result suggests the NMR ribosomes are intrinsically more accurate than the mouse ribosomes.

Discussion Our study shows an unexpected processing of NMR 28S rRNA into two molecules of 2.5 and 3 kb, corresponding to 3′ and 5′ segments of the rRNA gene. The break occurs within the D6 region and results from two cuts that excise a 263-nt fragment containing a unique 118-nt sequence. This 118-nt sequence bears no significant homology to other known sequences. The cleavage occurs in all NMR organs that we analyzed, as well as in the primary fibroblast cultures. The excised fragment is flanked by 5-nt direct repeats, which

supports that the insert originated from a transposable element. The sequences distal to the cut sites share weak homology with the mouse D6, suggesting that the insert could have also been introduced by homologous recombination. We do not know the mechanism responsible for the cleavage of the NMR 28S rRNA. The cuts could be generated by a ribozyme activity residing within the D6 sequence itself, as was reported for R2 trans- posons residing within the rDNA locus (25). We tested this possibility by in vitro transcribing 28S rRNA; however, the in vitro transcribed 28S rRNA molecules remained uncleaved. Fi- nally, the cleavage may be mediated by a ribonuclease that rec- ognizes the D6 region. Fragmented rRNA has been found in Salmonella (26), pro-

tozoa (27, 28), a worm (29), and several arthropods (ref. 30 and references therein). However, in vertebrates there is only one other example of a fragmented rRNA, found in rodents of the genus Ctenomys, also known as tuco-tuco (31). In Ctenomys, 28S rRNA is cleaved within the D6 region, but the cleavage occurs with only a single cut site, and sequences surrounding the break share no homology with the NMR sequence (Fig. 2A). Interestingly, Cten- omys is a South American rodent that is phylogenetically distant from the NMR but shares similar ecology. Ctenomys is a sub- terranean rodent, with some species of the genus leading a social lifestyle. The maximum lifespan in captivity is not determined for Ctenomys, but in the wild, Ctenomys is short-lived, with an average lifespan of 3 y (32). We found that NMR cells had equivalent translation rates

(adjusted for physiological temperature). However, by construct- ing a set of unique luciferase reporters, we observed that NMR cells have substantially higher translation fidelity than mouse cells,

RPS6

-actin

NMR Mouse

p = 0.02

p = 0.0027

B

0

5

10

15

20

25

Tr an

sl at

io n

ra te

, C

P M

/ μg

p ro

te in

1,

00 0

NMR SF 32oC

NMR SF 37oC

MSF 32oC

MSF 37oC

A

Fig. 3. NMR and mouse have similar translation rates. (A) [35S]-cystein/ methionine incorporation assay performed on growing fibroblasts showing NMR translation occurs at roughly equal speed to that of mouse fibroblasts. The assays were performed at 32 °C and 37 °C for both mouse and NMR cells. NMR SF are primary NMR skin fibroblast and MSF are primary mouse skin fibroblast. (B) Western blot against RPS6, a protein component of the ma- ture ribosome (Upper) and β-actin as loading control (Lower).

17352 | www.pnas.org/cgi/doi/10.1073/pnas.1313473110 Azpurua et al.

D ow

nl oa

de d

fr om

h tt

ps :/

/w w

w .p

na s.

or g

by 7

3. 20

3. 24

5. 0

on A

pr il

1 2,

2 02

2 fr

om I

P a

dd re

ss 7

3. 20

3. 24

5. 0.

even when transformed with oncogenes to increase their cell cycle speed. NMR cells also displayed higher fidelity of translation for all errors we tested, including misincorporation, stop codon skip- ping, and frameshifts. A potential caveat of our translation fidelity assay is that active firefly luciferase could be generated not only by incorrect translation but also by incorrect transcription or random mutations of the plasmid DNA during cell growth. This seems unlikely given that if replication or transcription produced active luciferase, we would not see the expected nucleotide position bias that the assay clearly shows (with the third position being the least mistranslated). Furthermore, we would not expect to see such a large bias in favor of insertions over deletions. The mis- incorporation rate we detect in our assay is also between 1:1,000 and 1:10,000, which is consistent with reports in the literature (24). Thus, we conclude that most of the luminescence from our assay comes from mistranslation, rather than from mistakes at the nucleotide level. Pérez et al. reported that the NMR proteome has higher levels

of cysteine than the mouse proteome and that the proteins are resistant to unfolding stressors such as 1 M urea. These results are consistent with high levels of translation fidelity (11). According to our results, the mouse proteome would have two- to 10-fold higher levels of proteins with misincorporated amino acids relative to the NMR. These misincorporated amino acids could increase the number of proteins that are unstable under challenging conditions, such as urea treatment. Pérez et al. also report that the NMR does not show an age-dependent increase in protein ubiquitination. This is consistent with the high fidelity of protein translation in the NMR. More recently, it has been reported that NMR proteins maintain a higher level of activity after carbonylation and oxidation than the equivalent mouse proteins (33). Again, it is possible that fewer amino acid mis- incorporations allow proteins to tolerate more stochastic improper modifications and oxidative damage, and thus retain tertiary struc- ture and activity. Given the remarkable longevity of the NMR, it is possible that

increased translational fidelity plays a role in its longevity. The “error catastrophe theory of aging” first proposed by Orgel in 1963 (34) and explored in Escherichia coli by Edelmann and Gallant in 1977 posits that early defects in translational fidelity lead to a vicious cycle of progressively more severe defects in protein homeostasis (35). Differences in translational fidelity may play an important role in determining lifespan. Many age- related disorders are caused by protein aggregation (21), and problems with translational fidelity can lead to the accumulation of misfolded protein and the formation of inclusion bodies (22). The proteome is maintained by various mechanisms, all of which show an age-dependent decline (36). During heat stress, chaper- ones are up-regulated, but this mechanism fails with age (37, 38). Proteasome activity also declines with age (39). Aging also leads to a decline in macroautophagy, which eliminates misfolded or in- correctly glycosylated proteins from the cell (40). It is possible that when proteome maintenance mechanisms fail, it allows the aberrant mistranslated proteins to accumulate. Species, such as the NMR, with a lower overall level of mistranslated proteins would thus be buffered against some of the decline of old age, given that they have fewer misfolded proteins to begin with. The main source of mistranslated proteins in E. coli is discrimi-

nation of noncognate or near-cognate tRNAs during translation (24). The ribosome structure itself confers at least some of the fidelity of the process, as shown by hyperaccurate ribosomes (41), and its sequence and structure (42) may be altered by natural se- lection in species that evolve long lifespans. Indeed, in E. coli, there is an observed evolutionary tradeoff between ribosomal speed and accuracy (41). Because the NMR in vitro translation system cur- rently cannot be reconstituted from purified components, we are unable to directly test whether 28S cleavage is responsible for im- proved translation fidelity in the NMR. The cleavage leads to

A

WT Control 1st 2nd 3rd Stop FS+ FS-

Firefly Luciferase Reporter

Fi re

fly /

R en

ill a

Lu m

in es

ce nc

e

0

0.01

0.02

0.03

0.04

5

10

15

20

* *

* * *

*

B

WT Control 1st 2nd 3rd Stop FS+ FS-

Firefly Luciferase Reporter

Fi re

fly /

R en

ill a

Lu m

in es

ce nc

e

0

0.01

0.02

0.03

0.04

5

10

15

20

*

*

* * * *

25

30

Fig. 4. NMR has higher fidelity of translation than the mouse. Translation fidelity of NMR and mouse skin fibroblasts was quantified as a frequency of amino acid misincorporation into a set of firefly luciferase reporters with a mutated critical lysine 529 (Fig. S3). Misincorporation events can rescue luciferase activity. Luminescence of mutated firefly luciferase was measured and normalized to the luminescence of Renilla luciferase, both under control of the CMV promoter. The positive control (wild-type luciferase) ratio was several orders of magnitude higher and is cropped for visibility. The mutant reporters E - K529E; I - K529I; N - K529N measured misincorporation into the first, second, and third codon positions, respectively; the STOP reporter contained a premature stop codon mutation at amino acid 81; FS+ con- tained a positive frameshift mutation by nucleotide insertion at position 81; and FS contained a negative frameshift mutation by nucleotide deletion at position 81. Results for the assay are shown for NMR and mouse cells grown at 32 °C (A) and 37 °C (B). Error bars show standard deviation. Asterisks in- dicate values significantly different between NMR and mouse (P < 0.005) Student t test.

Azpurua et al. PNAS | October 22, 2013 | vol. 110 | no. 43 | 17353

C EL L B IO LO

G Y

D ow

nl oa

de d

fr om

h tt

ps :/

/w w

w .p

na s.

or g

by 7

3. 20

3. 24

5. 0

on A

pr il

1 2,

2 02

2 fr

om I

P a

dd re

ss 7

3. 20

3. 24

5. 0.

deletion of the major part of the D6 region corresponding to a helix 45 of 23S rRNA in E. coli and leaves two fragments of 28S rRNA disconnected. This may change the folding or dy- namics of the large ribosomal subunit, altering the rate of GTP hydrolysis by elongation factor 1A and/or interaction of the large subunit with tRNA during accommodation, thus affect- ing the fidelity of protein synthesis. In particular, in bacterial ribosome, helix 45 of 23S rRNA is connected to the L11-stalk (L12 stalk in mammals), which is formed by helices 42, 43, and 44 through coaxial stacking with helices 40 and 41. Thus, the deletion of helix 45 may affect dynamics of the L11 stalk, which directly interacts with the elbow of the aminoacyl-tRNA in the intermediate (A/T) state of tRNA accommodation (43). In summary, we have shown that NMR has a unique form of 28S

rRNA that is processed in two fragments. The processing occurs by two cuts within the D6 region. Furthermore, we demonstrate that NMR has higher fidelity of translation than the mouse. We hy- pothesize that the unique structure of the NMR 28S rRNA may contribute to the higher fidelity of translation, which in turn may contribute to exceptional longevity of the NMR.

Materials and Methods RNA Extraction and Denaturing Gel Electrophoresis. RNA from rodent tissues and cells was extracted either by RNEasy Kit (Qiagen), according to the manufacturer’s instructions, or by TriReagent extraction. TriReagent ex- traction was performed by mixing cells or tissues with TriReagent, rupturing in a glass homogenizer, and centrifuging at high speed (>10,000 × g) at 4 °C on a tabletop centrifuge. The aqueous phase was then washed with iso- propanol to precipitate RNA and then washed with ethanol, dried, and resuspended in RNase-free water.

Denaturing agarose gels were prepared in double-distilled H2O and 20× Mops buffer diluted to 1× concentration. Formaldehyde was added to a final concentration of 2% (vol/vol). RNA samples were denatured in formaldehyde/Mops buffer at 80 °C for 5 min and then immediately placed on ice for at least 2 min. Samples were loaded and run at ∼100 V for 2–3 h.

rRNA Break Mapping. The unbroken 28S rRNA was estimated to be ∼5 kbp, and because the two fragments were close in size, we estimated that the break probably happened within 1,000 bp of the middle. 5′ RACE primers specific to highly conserved regions downstream of the putative 28S mid- point were designed. A 5′/3′ second-generation RACE kit (Roche) was used to obtain a fragment that mapped from the conserved region to the 5′ break point. Primers were designed that mapped upstream of the putative breakpoint, and 3′ RACE was performed by ligating an adenylated linker to the 3′ end of the ribosomal RNA. PCR was then used to extend from the upstream conserved site to the start point of the break. Primer sequences are shown in Table S2.

Tissue Culture. NMR fibroblasts were grown at 32 °C (in vivo body temperature, unless otherwise indicated) 5% CO2, 3% O2 on treated polystyrene culture dishes (Corning) in eagle’s minimum essential medium (EMEM) media (ATCC) supplemented with 15% FBS (Gibco), nonessential amino acids, sodium pyru- vate, 100 U/mL penicillin, and 100 μg/mL streptomycin (Gibco). Mouse cells were grown in identical conditions except at 37 °C, unless otherwise indicated.

Translation Rate Assay. Cells were grown to ∼75% confluence, and media was removed and incubated for 15 min with cysteine/methionine-free media to deplete intracellular pools. The media was then removed, and cells were in- cubated for 30 min in media containing [35S]-methionine/cystein. Radioactive media was then removed, and cells were incubated with 1.5 mL ice-cold lysis buffer and scraped from the plate, using a disposable plastic scraper. The cells were then stored at −80 °C overnight. Trichloroacetic acid (TCA) precipitation was used the following day to pellet proteins, and the pellet was washed with

acetone. Lowry assay was performed to determine protein concentration. Radioactivity was evaluated by scintillation counter (Beckman).

Western Blot. Proteins samples were run through 15% SDS/PAGE and transferred to nitrocellulose via the TransBlot Turbo system (BioRad). The membrane was cut and blotted separately for ribosomal protein RPS6 and beta-actin (loading control). The membrane was incubated with primary antibody overnight at 4 °C and with secondary antibody for 1 h at room temperature. Antibodies used were Santa Cruz biotechnology SC-47778 (beta-actin) and Cell Signaling 2317 (RPS6).

Plasmid Construction and Transfection. pGL3 wild-type luciferase plasmid (Promega) was mutated as indicated by site-directed mutagenesis, using a QuikChange kit (Agilent). Several clones were picked, miniprepped, and sent for sequencing to confirm proper mutagenesis. Correct plasmids were prepped using EndoFree Maxi prep kit (Qiagen). Cells were thawed and grown for two passages before being harvested. Cells were counted, and 1 × 106 cells were transfected using Nuclefector transfection (Lonza) with program U-020 in human dermal fibroblast (NHDF) solution with 5 μg reporter plasmid [mutant, wild-type-positive control, or hypoxanthine-guanine phosphoribosyltransferase (HPRT) negative control] and 0.1 μg Renilla luciferase plasmid (transfection/promoter control). Cells were then pipetted onto 10 cm Corning treated plates and allowed to grow for 24 h at the indicated temperature. Mutation sites within Luciferase gene are shown in Fig. S3.

Luciferase Assay. A Promega Dual-Luciferase assay kit was used to perform analysis on transfected cells. Cells were harvested, counted, and lysed, using the provided buffer. Extracts were thoroughly mixed with LAR-II reagent, and luminescence was measured. The sample was then mixed with STOP&GLO solution (halting firefly luciferase and starting Renilla luciferase), and lumi- nescence was measured. Ratio of firefly to Renilla was used as an indicator of translational fidelity.

Assembly of 28S rRNA Sequence from Illumina Raw Reads. First, we mapped the paired-end reads to the mouse 28S DNA sequence, using Bowtie2 (44), which allows mismatches and small indels/gaps. Then Velvet (45) was used to as- semble the reads that are mapped to the mouse 28S DNA sequence. The set of contigs with the largest N50 were selected. However, it is hard to obtain the reads from divergent regions by this method. To find scaffolds that the cover NMR 28S gene, we mapped our contigs to all of the scaffolds obtained by genome-wide assembly. Finally, we find one scaffold that matched our contigs the most.

Divergence Analysis of the Cleavage Sites in the D6 Region. To assess the di- vergence/conservation at the two cutting sites in D6 among all of the copies of 28S in NMR, we used the following strategy: Given a single sequenced copy of NMR D6 and its flanking regions, the conserved 50-bp sequences immediately upstream to the start of D6 and downstream to the end of D6 were used as the target for the 5′ and 3′ end cutting sites, separately. We mapped the target to the reads with the maximum 3 mismatch (i.e., we take the sequencing short reads as the references and do alignment for the target). After we get the matched reads, if the read covers the cutting site, we save the subsequence around the cutting site. We kept ±15 bp from the 5′ end cutting site and ±18 bp from the 3′ end cutting site, by which we centered the cutting site in the subsequences around the cutting sites. Finally, the multiple sequence alignment program Clustal Omega (46) was used to align the subsequences and the nucleotide frequency was counted based on the multiple sequence alignment result. Fig. S2 shows the divergence analysis strategy. Table S1 shows that nucleotide frequency at the two cutting sites. The cutting sites were highly conserved among the 28S rRNA copies.

ACKNOWLEDGMENTS. We thank members of the T. Eickbush and G. Culver laboratories at the University of Rochester for useful discussion and sugges- tions. This work was supported by the National Institutes of Health and Ellison Medical Foundation grants (to V.G.) and by Life Extension Foundation grants (to V.G. and A.S).

1. Buffenstein R (2005) The naked mole-rat: A new long-living model for human aging

research. J Gerontol A Biol Sci Med Sci 60(11):1369–1377. 2. Lewis KN, Andziak B, Yang T, Buffenstein R (2012) The Naked Mole-Rat Response to

Oxidative Stress: Just Deal with It. Antioxid Redox Signal, 10.1089/ars.2012.4911. 3. Delaney MA, Nagy L, Kinsel MJ, Treuting PM (2013) Spontaneous histologic lesions of

the adult naked mole rat (Heterocephalus glaber): A retrospective survey of lesions in

a zoo population. Vet Pathol 50(4):607–621.

4. Liang S, Mele J, Wu Y, Buffenstein R, Hornsby PJ (2010) Resistance to experimental

tumorigenesis in cells of a long-lived mammal, the naked mole-rat (Heterocephalus

glaber). Aging Cell 9(4):626–635. 5. Seluanov A, et al. (2009) Hypersensitivity to contact inhibition provides a clue to cancer

resistance of naked mole-rat. Proc Natl Acad Sci USA 106(46):19352–19357. 6. Tian X, et al. (2013) High-molecular-mass hyaluronan mediates the cancer resistance

of the naked mole rat. Nature 499(7458):346–349.

17354 | www.pnas.org/cgi/doi/10.1073/pnas.1313473110 Azpurua et al.

D ow

nl oa

de d

fr om

h tt

ps :/

/w w

w .p

na s.

or g

by 7

3. 20

3. 24

5. 0

on A

pr il

1 2,

2 02

2 fr

om I

P a

dd re

ss 7

3. 20

3. 24

5. 0.

7. Azpurua J, Seluanov A (2013) Long-lived cancer-resistant rodents as new model species for cancer research. Frontiers Genet 3:319.

8. Hulbert AJ, Faulks SC, Buffenstein R (2006) Oxidation-resistant membrane phospho- lipids can explain longevity differences among the longest-living rodents and simi- larly-sized mice. J Gerontol A Biol Sci Med Sci 61(10):1009–1018.

9. Andziak B, et al. (2006) High oxidative damage levels in the longest-living rodent, the naked mole-rat. Aging Cell 5(6):463–471.

10. Andziak B, Buffenstein R (2006) Disparate patterns of age-related changes in lipid peroxidation in long-lived naked mole-rats and shorter-lived mice. Aging Cell 5(6): 525–532.

11. Pérez VI, et al. (2009) Protein stability and resistance to oxidative stress are deter- minants of longevity in the longest-living rodent, the naked mole-rat. Proc Natl Acad Sci USA 106(9):3059–3064.

12. Mamane Y, Petroulakis E, LeBacquer O, Sonenberg N (2006) mTOR, translation ini- tiation and cancer. Oncogene 25(48):6416–6422.

13. Johnson SC, Rabinovitch PS, Kaeberlein M (2013) mTOR is a key modulator of ageing and age-related disease. Nature 493(7432):338–345.

14. Kaeberlein M, Kennedy BK (2011) Hot topics in aging research: Protein translation and TOR signaling, 2010. Aging Cell 10(2):185–190.

15. Zoncu R, Efeyan A, Sabatini DM (2011) mTOR: From growth signal integration to cancer, diabetes and ageing. Nat Rev Mol Cell Biol 12(1):21–35.

16. Morrow G, Samson M, Michaud S, Tanguay RM (2004) Overexpression of the small mitochondrial Hsp22 extends Drosophila life span and increases resistance to oxida- tive stress. FASEB J 18(3):598–599.

17. Morley JF, Morimoto RI (2004) Regulation of longevity in Caenorhabditis elegans by heat shock factor and molecular chaperones. Mol Biol Cell 15(2):657–664.

18. Hsu AL, Murphy CT, Kenyon C (2003) Regulation of aging and age-related disease by DAF-16 and heat-shock factor. Science 300(5622):1142–1145.

19. Lindner AB, Madden R, Demarez A, Stewart EJ, Taddei F (2008) Asymmetric segre- gation of protein aggregates is associated with cellular aging and rejuvenation. Proc Natl Acad Sci USA 105(8):3076–3081.

20. Budovsky A, Fraifeld VE, Aronov S (2011) Linking cell polarity, aging and re- juvenation. Biogerontology 12(2):167–175.

21. Grune T, Jung T, Merker K, Davies KJ (2004) Decreased proteolysis caused by protein aggregates, inclusion bodies, plaques, lipofuscin, ceroid, and ‘aggresomes’ during oxidative stress, aging, and disease. Int J Biochem Cell Biol 36(12):2519–2530.

22. Lee JW, et al. (2006) Editing-defective tRNA synthetase causes protein misfolding and neurodegeneration. Nature 443(7107):50–55.

23. Kim EB, et al. (2011) Genome sequencing reveals insights into physiology and lon- gevity of the naked mole rat. Nature 479(7372):223–227.

24. Kramer EB, Farabaugh PJ (2007) The frequency of translational misreading errors in E. coli is largely determined by tRNA competition. RNA 13(1):87–96.

25. Eickbush DG, Eickbush TH (2010) R2 retrotransposons encode a self-cleaving ribozyme for processing from an rRNA cotranscript. Mol Cell Biol 30(13):3142–3150.

26. Burgin AB, Parodos K, Lane DJ, Pace NR (1990) The excision of intervening sequences from Salmonella 23S ribosomal RNA. Cell 60(3):405–414.

27. Schnare MN, Gray MW (1990) Sixteen discrete RNA components in the cytoplasmic ribosome of Euglena gracilis. J Mol Biol 215(1):73–83.

28. Boer PH, Gray MW (1988) Scrambled ribosomal RNA gene pieces in Chlamydomonas reinhardtii mitochondrial DNA. Cell 55(3):399–411.

29. van Keulen H, Mertz PM, LoVerde PT, Shi H, Rekosh DM (1991) Characterization of a 54-nucleotide gap region in the 28S rRNA gene of Schistosoma mansoni. Mol Bio- chem Parasitol 45(2):205–214.

30. Fujiwara H, Ishikawa H (1986) Molecular mechanism of introduction of the hidden break into the 28S rRNA of insects: Implication based on structural studies. Nucleic Acids Res 14(16):6393–6401.

31. Melen GJ, Pesce CG, Rossi MS, Kornblihtt AR (1999) Novel processing in a mammalian nuclear 28S pre-rRNA: Tissue-specific elimination of an ‘intron’ bearing a hidden break site. EMBO J 18(11):3107–3118.

32. Nowak R (1999) Walker’s Mammals of the World (John Hopkins Univ Press, Baltimore). 33. De Waal EM, et al. (2013) Elevated protein carbonylation and oxidative stress do not

affect protein structure and function in the long-living naked-mole rat: A proteomic approach. Biochem Biophys Res Commun 434(4):815–819.

34. Orgel LE (1963) The maintenance of the accuracy of protein synthesis and its rele- vance to ageing. Proc Natl Acad Sci USA 49:517–521.

35. Edelmann P, Gallant J (1977) On the translational error theory of aging. Proc Natl Acad Sci USA 74(8):3396–3398.

36. Koga H, Kaushik S, Cuervo AM (2011) Protein homeostasis and aging: The importance of exquisite quality control. Ageing Res Rev 10(2):205–215.

37. Arslan MA, Csermely P, Soti C (2006) Protein homeostasis and molecular chaperones in aging. Biogerontology 7(5-6):383–389.

38. Shpund S, Gershon D (1997) Alterations in the chaperone activity of HSP70 in aging organisms. Arch Gerontol Geriatr 24(2):125–131.

39. Baraibar MA, Friguet B (2012) Changes of the proteasomal system during the aging process. Prog Mol Biol Transl Sci 109:249–275.

40. Cuervo AM, et al. (2005) Autophagy and aging: The importance of maintaining “clean” cells. Autophagy 1(3):131–140.

41. von Ahsen U (1998) Translational fidelity: Error-prone versus hyper-accurate ribo- somes. Chem Biol 5(1):R3–R6.

42. Vila-Sanjurjo A, et al. (2003) X-ray crystal structures of the WT and a hyper-accurate ribosome from Escherichia coli. Proc Natl Acad Sci USA 100(15):8682–8687.

43. Schmeing TM, et al. (2009) The crystal structure of the ribosome bound to EF-Tu and aminoacyl-tRNA. Science 326(5953):688–694.

44. Langmead B, Salzberg SL (2012) Fast gapped-read alignment with Bowtie 2. Nat Methods 9(4):357–359.

45. Zerbino DR, Birney E (2012) Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 18(5):821–829.

46. Sievers F, et al. (2011) Fast, scalable generation of high-quality protein multiple se- quence alignments using Clustal Omega. Mol Syst Biol 7:539.

Azpurua et al. PNAS | October 22, 2013 | vol. 110 | no. 43 | 17355

C EL L B IO LO

G Y

D ow

nl oa

de d

fr om

h tt

ps :/

/w w

w .p

na s.

or g

by 7

3. 20

3. 24

5. 0

on A

pr il

1 2,

2 02

2 fr

om I

P a

dd re

ss 7

3. 20

3. 24

5. 0.