Main

The NMD pathway is conserved in all eukaryotic organisms that have been studied1,2. Strict evolutionary maintenance of this function may result from its ability to protect cells from deleterious truncated proteins that would be produced from stable nonsense transcripts. There are many well-studied examples of human phenotypes resulting from nonsense or frameshift mutations that are modulated by NMD6,7. It seems unlikely, however, that the rare frequency of de novo nonsense mutations could provide sufficient selective pressure for complete evolutionary maintenance of NMD. Another view posits that physiologic transcripts that structurally mimic nonsense transcripts are the predominant substrates for NMD. Observations in model systems support this hypothesis. For example, expression profiling of NMD-deficient yeast strains found altered expression of 10% of the transcriptome, allowing identification of specific classes of physiologic NMD substrates8,9. In mammals, a distinct mechanism for recognizing nonsense codons that relies on splicing is used10. As yeast NMD occurs independently of splicing11, unique classes of endogenous NMD substrates probably exist in higher eukaryotes.

We used microarrays to examine transcript expression profiles after short interfering RNA (siRNA)-mediated depletion of Rent1 in HeLa cells3. We expressed 4,000 transcripts, representing 33% of the probe sets on the arrays, in HeLa cells and then assayed them in these experiments. We found that 197 transcripts (4.9%) were consistently upregulated (Supplementary Table 1 online) and 176 transcripts (4.4%) were consistently downregulated (Supplementary Table 2 online) by a factor of at least 1.9 in duplicate experiments.

An intron located at least 50–55 nucleotides downstream of a termination codon is sufficient to initiate mammalian NMD5. More than half of the upregulated transcripts (104 of 197) had identifiable features that satisfied this constraint (Table 1 and Supplementary Table 1 online). Putative NMD-inducing features included upstream open reading frames (uORFs; 70 transcripts), alternative splicing that introduces nonsense codons or frameshifts (21 transcripts, some of which undergo alternative splicing specifically in HeLa cells12,13) and introns in the 3′ untranslated region (UTR; 9 transcripts). Two transcripts encoding proteins containing selenocysteine were also upregulated, consistent with the alternative recognition of the UGA selenocysteine codon as a signal for translation termination14. Additionally, we observed upregulation of nonfunctional transcripts derived from mariner 2 transposon (HSMAR2) remnants and human endogenous retrovirus H (HERV-H) sequences that have acquired premature termination codons (PTCs) during evolution. Previous studies15,16,17 and our analyses of expressed-sequence tags showed that at least some HSMAR2 and HERV-H transcripts undergo splicing. Accordingly, we observed 1.0-kb transcripts corresponding to the size of the spliced HERV-H message (Fig. 1c).

Table 1 Classes and selected examples of putative NMD-regulated transcripts
Figure 1: Prolonged decay rates and increased steady-state abundance of transcripts identified through expression profiling of Rent1-depleted cells.
figure 1

(a) Half-lives of selected upregulated transcripts. Seventy-two hours after treatment with RNA interference directed against firefly luciferase (Luc) or RENT1, transcription was inhibited with 5,6-dichloro-1-β-D-ribofuranosylbenzimidazole. RNA was collected at the indicated time points and transcript levels were monitored by northern blotting. Levels of GAPD, a stable housekeeping transcript, were used to correct for differences in RNA loading. (b) Northern blots shown in a were quantified, and values are represented as semi-log plots. All transcripts showed exponential decay except SLC3A2, which showed a biphasic decay pattern consistent with the behavior of other known nonsense transcripts30. (c) Northern-blot analysis of transcripts derived from endogenous retroviral sequences (HERV-H) in Rent1-depleted cells. (d) Steady-state abundance of selected upregulated transcripts in Rent1-depleted and Rent2-depleted cells. Two independent experiments are shown for each experimental condition (1 and 2). To quantify each transcript (normalized to GAPD), one luciferase siRNA–treated sample was arbitrarily set to 100 and other signals were adjusted relative to this value.

To confirm that transcripts that were upregulated in cells depleted of Rent1 were primary substrates of the NMD pathway, we determined representative mRNA decay rates (Fig. 1a,b and Table 2). All 14 transcripts that we examined had prolonged half-lives in Rent1-depleted cells, consistent with direct regulation by nonsense surveillance. Reliable half-lives of transposon- or retrovirus-derived transcripts could not be measured as they are encoded by thousands of loci15,16 that cannot be individually discriminated by northern blotting and probably have widely variable decay rates. We observed increased steady-state abundance of these RNAs, validating the microarray results (Fig. 1c,d). To confirm that the upregulated transcripts represent bona fide NMD substrates, we examined the steady-state abundance of a representative set in cells depleted of a second essential trans-effector of nonsense surveillance, Rent2 (also called hUpf2). We observed concordant upregulation of all five transcripts examined (Fig. 1d). We also observed increased expression of selected transcription factors (e.g., ATF3, ATF4, TCFL4 and CEBPG), indicating that there are probably also some indirect effects on transcript abundance.

Table 2 Half-lives of putative NMD-regulated transcripts

In contrast to the upregulated transcripts, the seven downregulated transcripts that we tested either showed no change in steady-state abundance or stability on depletion of Rent1 or showed no concordant changes on depletion of Rent2, as assessed by northern blotting (data not shown). Thus, as in yeast8, we conclude that few, if any, mammalian transcripts that showed decreased abundance by microarray analysis are regulated by NMD.

As in yeast8, we observed upregulation of numerous transposon-derived RNAs and transcripts containing uORFs in NMD-deficient mammalian cells, indicating that these are conserved classes of NMD-regulated transcripts. The other classes of NMD substrates in yeast were not identified in our analyses (e.g., pseudogene and polycistronic transcripts; and transcripts with inefficient pre-mRNA splicing, leaky scanning or +1 frameshifting); this absence probably reflects the limited composition of current mammalian microarrays and incomplete annotation of the human genome.

The requirement for splicing in mammalian nonsense-codon recognition has given rise to new classes of NMD substrates, including alternatively spliced transcripts and transcripts with introns in their 3′ UTRs. Bioinformatic analyses have suggested that introduction of PTCs through alternative splicing is a relatively common feature of the human transcriptome18. When coupled with NMD, this provides a mechanism to achieve regulated titration of gene expression using the complex alternative splicing machinery, while avoiding the potentially deleterious consequences of expression of truncated proteins.

Among physiologic NMD substrates, we observed increased representation of genes involved in amino acid metabolism (Table 3). Of the 80 genes on the array assigned to gene ontology terms 'amino acid transport', 'amino acid biosynthesis' and 'amino acid activation', 12 (15%) were upregulated after abrogation of NMD. Given that only 4.9% of assayed transcripts were upregulated overall, this is a significant enrichment (χ2 = 16.3, P < 10−4). We also observed upregulation of two transcription factors that coordinate cellular responses to amino acid starvation (ATF3 and ATF4; refs. 19,20). Depletion of amino acids inhibits translation21, and as NMD requires ongoing translation, regulation of these transcripts by nonsense surveillance probably couples their expression level to translational efficiency. Thus, under conditions of amino acid starvation, inhibition of translation and NMD increases expression of transcripts that promote restoration of amino acid homeostasis. In keeping with this hypothesis, we observed attenuation of NMD under conditions of amino acid starvation, as assessed by the increased abundance of nonsense transcripts derived from a β-globin minigene and the increased abundance of representative NMD-regulated transcripts both with and without functions related to amino acid homeostasis (Fig. 2). As expected, some transcripts related to amino acid homeostasis that were not upregulated by NMD inhibition showed increased abundance in response to starvation (e.g., MARS), whereas others did not (e.g., KARS). These data attest to the complexity of the response to nutrient deprivation. Our observations support the view that NMD contributes to the regulation of transcripts that encode key regulators of the starvation response. Moreover, there is adequate precedent that starvation-induced translational inhibition can be associated with an increased translational yield for selected transcripts, including ATF3 and ATF4 (refs. 19,22).

Table 3 Upregulated human transcripts with amino acid metabolism and transport functions
Figure 2: Amino acid starvation inhibits NMD and upregulates transcripts required for amino acid homeostasis.
figure 2

(a) Cells were transfected with a wild-type (WT) or nonsense-containing (PTC) β-globin minigene and grown in the presence (Fed) or absence (Starved) of amino acids. The steady-state abundance of β-globin messages was determined by northern blotting. Levels of neomycin phosphotransferase II (Neo), also encoded on the minigene plasmid, were used to control for differences in transfection efficiency and loading. (b) Northern-blot analysis of physiologic NMD substrates in the presence and absence of amino acids. The six transcripts in the upper portion of the figure (group 1) have functions related to amino acid homeostasis, whereas UPP1 and HSMAR2 (group 2) do not. MARS and KARS encode transcription factors related to amino acid metabolism but were not upregulated by NMD inhibition (group 3). β-actin serves as a loading control.

Notably, analysis of a previously determined set of yeast genes that have increased expression in response to NMD inhibition8 uncovered evidence that this mechanism to preserve amino acid homeostasis is evolutionarily conserved (Table 4). Upregulated transcripts include SSY1, encoding a sensor of external amino acid concentration that also coordinates transcription of amino acid permease genes in response to nutrient deprivation; LST7, encoding a positive regulator of transport of amino acid permeases from the Golgi to the cell surface; and AGP3, encoding a low-affinity permease that might have a prominent role in amino acid transport induced by nitrogen starvation23,24. Other upregulated transcripts encode factors involved in the transport and catabolism of secondary nitrogen sources25, including the only known allantoate permeases, two putative allantoate permeases, five of six factors known to mediate allantoin degradation, a urea transport and degradation enzyme, and a sensor and transporter of ammonia. NMD inhibition was also associated with increased expression of crucial mediators of autophagy (bulk vacuolar degradation of cytosolic proteins to recycle amino acids in response to starvation) and sporulation. Multiple yeast transcripts that show increased expression after NMD inhibition (e.g., TRP1, ATG12 (also called APG12) and DAL5) have been shown experimentally to lack responsiveness to Gcn4p, a master transcriptional activator of amino acid biosynthetic genes in multiple pathways26, suggestive of complementarity between transcriptional and post-transcriptional responses to amino acid deprivation. Consistent with the paradigm that NMD has been functionally incorporated into homeostatic mechanisms, we also observed upregulation of the human homolog of Caenorhabditis elegans smg-5, a gene that is essential for mammalian NMD27, in NMD-deficient HeLa cells. This may be a mechanism by which NMD regulates its own efficiency.

Table 4 Starvation-response yeast transcripts upregulated upon NMD inhibition

Our results suggest that the predominant role of mammalian NMD is to regulate the expression of thousands of physiologic transcripts. In this view, its ability to modify disease caused by nonsense mutations, though medically important, is inconsequential for evolutionary selection. Given the number and diversity of transcripts regulated by NMD, the process has probably been functionally integrated into additional developmental and homeostatic mechanisms. This may explain the unique dependence of mammals on NMD for viability28 and limit therapeutic strategies based on modulation of the pathway. NMD also mutes a noisy genome and may thus provide tolerance for nonproductive events that occur through genome evolution. This may serve to buffer and hence maintain variation that proves beneficial only under selected physiologic or environmental conditions.

Note: Supplementary information is available on the Nature Genetics website.

Methods

Cell culture, siRNA treatment and RNA analysis.

We cultured HeLa cells in Dulbecco's modified Eagle medium supplemented with 10% fetal bovine serum and carried out siRNA-mediated silencing of Rent1 and Rent2 as described3. We isolated RNA using Trizol (Invitrogen) 72 h after siRNA treatment and purified it using the RNeasy Mini Kit (Qiagen). As a positive control for the efficacy of siRNA treatment, we transfected an aliquot of the cells treated by RNA interference with a PTC-containing T-cell receptor β minigene (a gift from M. Wilkinson, University of Texas, MD Anderson Cancer Center) as described3 and determined the efficiency of NMD by northern blotting. We observed potent inhibition of NMD in the cells treated with siRNA (data not shown). Successful siRNA treatment was also indicated by decreased RENT1 mRNA levels, as assessed by microarray analysis (Supplementary Table 2 online). For mRNA half-life experiments, we treated cells with 100 μg ml of 5,6-dichloro-1-β-D-ribofuranosylbenzimidazole and collected RNA. For amino acid starvation experiments, we starved cells for 20 h in Krebs-Ringer bicarbonate buffer supplemented with 10% dialyzed fetal bovine serum. We grew fed cells in parallel in Dulbecco's modified Eagle medium supplemented with 10% dialyzed fetal bovine serum. We obtained the β-globin minigene from J. Lykke-Andersen (University of Colorado) and J. Steitz (Yale University). We generated northern-blot probes by RT-PCR (primer sequences available on request) and carried out hybridizations with UltraHyb (Ambion). We quantified radioactive signals using an Instant Imager (Packard) and calculated half-lives as described29.

Microarray hybridization and data analysis.

We processed total RNA from control and experimental cell preparations (two independent replicates for each) using single-round RNA amplification protocols as described in the Affymetrix GeneChip Expression Analysis Technical Manual. We used the SuperScript Choice System (Invitrogen) to synthesize first-strand cDNA from 5 μg of total RNA using a primer with oligo dT and T7 promoter sequences (Proligo LLC). After double-stranded cDNA synthesis, we purified the product by phenol-chloroform extraction and generated biotinylated antisense cRNA through in vitro transcription using the BioArray RNA High Yield Transcript Labelling kit (ENZO Life Sciences). After fragmenting cRNA at 94 °C for 35 min in 100 mM Tris-acetate (pH 8.2), 500 mM potassium acetate and 150 mM magnesium acetate, we hybridized 10 μg of cRNA to the Affymetrix human genome GeneChip array U95Av2 for 16 h at 45 °C with constant rotation. We used a Fluidics Station 400 (Affymetrix) to wash and stain the chips. Fluorescence was detected using a G2500 GeneArray Scanner (Hewlett-Packard) and image analysis of each GeneChip was done through the Microarray Suite 5.0 software from Affymetrix (MAS 5.0). For comparison between different chips, we used global scaling with a user-defined target intensity of 150.

We carried out quality control analysis of RNA samples using the Agilent Bioanalyzer Lab on a Chip (Agilent Technologies), confirming that all samples had optimal rRNA ratios (1:2 for 18S:28S). For quality control of hybridizations, GeneChip images and comparisons between chips, we studied the following parameters: scaling factor (all resulting values were within comparable range), background (all values were between 36 and 48), percentage of present calls (all chips had 31–33% present), housekeeping genes (3′:5′ ratios of GAPD were consistently between 1.1 and 1.3) and presence or absence of internal spike controls (Bio B was present in all cases).

The initial analysis of the expression results was based on pairwise comparisons among the different experimental conditions. Any transcript whose expression level changed by a factor of at least 1.9 between experimental sample and control sample was considered to be significantly differentially expressed. We chose 1.9 as the cut-off value because at or above this value, transcripts showed full reproducibility in our replicate analyses. The results obtained from the duplicate samples were filtered independently for significance on each of the four iterative pairwise comparisons. We selected transcripts that were significant in at least two of the four comparisons for the final gene candidate list. The relative change values were converted from the average signal log ratio values obtained from MAS 5.0.

Analysis of genomic organization and alternative splicing of upregulated transcripts.

We used publicly available databases, including the National Center for Biotechnology Information and the University of California Santa Cruz genome browser, to determine transcript structures.

Accession numbers.

GenBank: spliced HSMAR2 expressed-sequence tags, T98067, T40030, R96154, U92014 and U92019. GEO: raw microarray data: luciferase control 1, GSM29530; luciferase control 2, GSM29531; RENT1-1, GSM29532; RENT1-2, GSM29534; series 'determination of mRNA transcripts in HeLa cells that are regulated by RENT1', GSE1703.