BACKGROUND
nonsense-mediated decay  is a major eukaryotic surveillance mechanism that targets for degradation transcripts that harbor a premature termination codon .
a primary role of this mechanism is in posttranscriptional quality control, preventing the formation of truncated proteins that are potentially detrimental to the cell .
in addition, nmd is often combined with alternative splicing to form an important regulatory program of gene expression .
altogether, nmd is central to maintaining normal cellular activity, and disruption of its proper function is estimated to be associated with about one third of inherited genetic disorders, as well as with many forms of cancer .
however, the decision-making process at the basis of nmd, determining whether a termination codon  is premature or not, is still not fully understood.
it is generally believed that the context of the tc within the transcript determines whether it is normal or premature, but the nature of the contextual signals is unclear.
the leading hypothesis in mammals is that a tc is recognized as premature if the pre-mrna harbors an intron more than 50– <dig> bases downstream of the tc  <cit> .
the splicing reaction usually leaves traces in the form of a protein complex called the exon junction complex , which is deposited upon the mrna ~20– <dig> bases upstream of the exon-exon junction  <cit> .
it is believed that the nmd-triggering feature is the presence of ejcs in the 3’ untranslated region  at the time of translation termination.
the underlying mechanistic explanation is that the first ribosome that scans the transcript clears ejcs from the coding sequence , but is unable to clear ejcs that are downstream enough to the tc  <cit>  .
this model explains why the presence of a translated upstream open reading frame  is also known to trigger nmd , as the downstream cds likely harbors ejcs which behave, in the context of the tuorf, as 3’utr ejcs.
in fact, it was suggested that many of the nmd-targets result from short tuorfs, leading to an enrichment of transcripts with short cds among nmd-targets  <cit> , explaining why short cds was also reported as nmd-triggering feature in yeast  <cit>  and drosophila  <cit> .
despite overall support, this model of ptc recognition is incomplete, as in various cases it was reported to be violated: some transcripts trigger nmd despite of not harboring 3’utr ejcs _enref_ <dig>  while others are known to evade nmd despite of presenting 3’utr ejcs  or tuorfs  <cit> .
these exceptions brought about the idea that nmd works through several different pathways, and led to a distinction between ejc-dependent nmd – which is the nmd that degrades transcripts that harbor 3’utr ejcs – and ejc-independent nmd, or failsafe nmd – which is the nmd that works on transcripts that lack 3’utr ejcs  <cit> .
an alternative model , known as the faux-utr model, suggests that a normal transcript is characterized by certain signals in its 3’utr, and that disruption of these signals triggers nmd  <cit> _enref_ <dig>  substantial evidence supports the notion that one of the central signals that marks a transcript as normal is physical proximity between the tc and the poly binding protein pabpc <dig> , and that a tc is tagged as premature by the lack of such physical proximity.
presumably, in a normal transcript the translation termination factor erf <dig> binds pabpc <dig>  whereas anything that promotes nmd, like an ejc situated downstream of the tc, competes with the pabpc <dig> and allows erf3-bound up-frameshift protein  <dig>  to trigger nmd.
this model not only clarifies the part that 3’utr ejcs play in nmd, but also explains why very long 3’utrs are also known to trigger nmd  <cit> , as the tc and pabpc <dig> are separated by a very long physical distance.
however, as for other nmd-triggering features, a long 3’utr does not always result in nmd  <cit> .
these models assume that upf <dig>  on which nmd critically depends, works through specific binding with the translation termination complex  <cit> _enref_ <dig>  further evidence, however, points to the possibility that upf <dig> is not only recruited by the terminating ribosome, but is associated directly with the mrna, and is even thought to have helicase activity allowing it to slide along it  <cit> .
it is still debated whether this association is translation-dependent  or independent  <cit> , but all agree that upf <dig> is displaced from the cds during translation, leading to its enrichment along the 3’utr  <cit> .
opinions also diverge regarding the question whether nmd-targets have excess of upf <dig> molecules bound to their 3’utr.
some report that upf <dig> density along the 3’utr is higher in nmd-targets  <cit> _enref_ <dig>  others do not see significant difference in upf <dig> densities between nmd-targets and nmd-insensitive transcripts, claiming that it is not the upf <dig> density that matters for nmd, but rather the density of the activated  upf <dig>  <cit> .
nmd may be directly linked to the efficiency of translation termination.
it was proposed that unfavorable context of the tc may affect the recruitment of terminating factors, leading to changes in the kinetics of translation termination, and eventually to ribosome stalling at the tc that is thought to be the event triggering nmd  <cit> .
combined with the presence of 3’utr-bound upf <dig>  this led to the suggestion that prolonged translation termination allows more upf <dig> molecules to bind to the 3’utr and become activated via regulated phosphorylation, resulting in activation of nmd  <cit> _enref_ <dig>
one of the biggest challenges in understanding ptc recognition stems from the observation that neither of the proposed nmd-triggering features characterizes a large fraction of nmd-targets.
even when all features are considered together, it was estimated that they characterize roughly 30% of nmd-targets  <cit> .
this suggests that there may be many more yet unidentified features that are important for ptc recognition.
surprisingly, the nucleotide composition around the tc has not received much attention.
hurt et al.
found that upf <dig> tends to bind 3’utrs in g-rich regions, possibly because of pausing of upf <dig> scanning in these regions  <cit> .
the same authors also reported that nmd-insensitive transcripts with long 3’utrs are enriched in poly stretches, perhaps because those stretches recruit pabpc <dig>  thus marking the tc as normal  <cit> .
studies on codon usage bias show that codons near the tc tend to be at-rich , conceivably as a result of selection against the formation of rna secondary structures near the tc that may interfere with the recruitments of release factors and lead to improper translation termination.
yet, the nucleotide composition of the sequence upstream the stop codon in the coding region was never regarded as a potential nmd triggering feature.
here, we wished to identify additional features that are involved in the decision-making process of nmd.
for this, we have used available genome-wide data on rna stability following upf <dig> knockdown  <cit> , and identified transcripts that are targets for degradation.
specifically, we focused on transcripts that lack 3’utr ejcs, and therefore our targets degrade through the ejc-independent nmd.
for each transcript we have computed a long list of features, and tested their relevance to degradation by evaluating whether any of them shows a unique distribution within nmd-targets, when compared to non-targets.
the features we have tested include all those that had been previously described as associated with nmd, as well as many new ones, mainly measuring nucleotide composition around the tc.
we found that nmd-targets are characterized by high g content around the tc and throughout the 3’utr.
however, we noticed a group of transcripts that are not nmd-targets and yet demonstrate the same, albeit weaker, pattern of enrichment.
interestingly, these transcripts are not nmd-insensitive either, but show decreased mrna stability following upf <dig> knockdown, and are thus called ‘destabilized’.
these transcripts might show instability due to factors other than nmd.
surprisingly, the 3’utr length is similar for nmd-targets and nmd-insensitive transcripts, but it is substantially shorter for the destabilized transcripts.
based on these observations we suggest that high g content around the tc and throughout the 3’utr leads to aberrant translation termination that triggers nmd.
we show that this model has a far greater discriminative power than any other ejc-independent nmd-triggering feature known today.
RESULTS
identifying nmd-targets
the ability to identify transcripts that are nmd-targets is a critical starting point in any analysis of nmd-triggering features.
nmd-targets are traditionally identified as transcripts whose characteristic behavior is affected by knocking-down a core nmd component.
the most common approach is to mark as nmd-targets transcripts whose expression level is increased following upf1-knockdown.
however, this approach is expected to have a high rate of false positives due to transcripts that are secondarily up-regulated.
moreover, the false negative rate may also be high due to genuine nmd-targets whose expression level is stabilized due to feedback regulation  <cit> .
recently, a different approach was introduced by tani et al.
<cit> , who used 5’-bromo-uridine immunoprecipitation chase-deep sequencing  at four time points to measure genome-wide mrna stability.
they compared control to upf1-knockdown samples, and marked as nmd-targets transcripts whose half-life was significantly increased following upf1-knockdown.
while this approach may still suffer from secondary effects, it is nonetheless a more direct measure of rna stability.
the current analysis starts with the raw data of tani et al.
<cit> .
we have followed parts of their analysis pipeline, but developed a novel method to detect nmd-targets.
tani et al.
estimated the decay rate of each transcript by modeling its mrna abundance decay curve.
however, decay curves of nmd-targets often deviate from a simple exponential model, reflecting a mixtures of currently degrading transcripts with older transcripts that had escaped degradation  <cit> .
given that tani et al.
only used four time points, their half-life estimates – especially when the decay is not a simple exponent– may have large standard errors.
to circumvent this, we did not attempt to estimate decay rates, and instead developed a test that compares the shapes of the mrna abundance decay curves and marks as nmd-targets transcripts whose decay curve following upf1-knockdown is significantly above their control decay curve .
this group of transcripts is denoted ‘stabilized’.fig.
<dig> schematic description of the transcript classification procedure.
red denotes mrna abundance decay curve in the control experiments, blue denotes mrna abundance decay curve following upf <dig> knockdown.
a stabilized transcripts .
b destabilized transcripts.
c nmd-insensitive  transcripts
theoretically, knocking down upf <dig> should result in stability increase of nmd-targets, and stability stasis for nmd-insensitive transcripts.
however, by using the technique above in the reverse direction we found transcripts whose stability decreases upon upf1-knockdown, namely, their mrna abundance decay curves following upf1-knockdown are significantly below their control decay curves .
such transcripts are unexpected by current theory, and are denoted here ‘destabilized’.
transcripts whose decay curves did not change significantly following upf1-knockdown were marked as ‘unaffected’, and represent transcripts that are not targeted for degradation by nmd .
finally, transcripts whose decay curves showed inconsistent behavior were left unclassified and were excluded from further analysis.
we have re-analyzed the data of tani et al.
and obtained mrna abundance information for a total of  <dig> transcripts.
from these, we removed  <dig> whose mrna abundance seemed to increase with time.
of the remaining  <dig> transcripts,  <dig>  were labeled as stabilized ,  <dig>  were labeled as destabilized,  <dig>  were labeled as unaffected, and  <dig>  were left unclassified.
in total, we could classify  <dig> transcripts for which we carried out all subsequent analyses.
to test the robustness of our results, we have repeated the process using a different, stricter, classification criterion totaling in  <dig> classifiable transcripts .
generating a set of nmd-triggering features
for each of the  <dig> classified transcripts we computed the following features, which were all previously suggested as characterizing nmd-targets in mammals:
tuorf.
a binary feature, indicating the presence/absence of translated uorfs strictly within the 5’utr.
3’utr length, in bases.
3’utr a and g content.
a and g density at the entire sequence of the 3’utr.
cds length, in bases.
alu.
a binary feature, indicating the presence/absence of alu elements in the 3’utr.
these elements are not directly associated with nmd, but are nevertheless thought to be involved in mrna degradation through staufen-mediated decay  <cit>  .
alu density.
the density of alu elements within the 3’utr.
we did not consider 3’utr ejc presence/absence as a feature, because the transcript list used by tani et al.
was based on the refseq database, which normally excludes transcripts that harbor 3’utr ejcs  <cit> .
indeed, our data include only  <dig> transcripts with 3’utr ejcs.
ten of them were classified as unaffected, whereas the other four were left unclassified.
this may reflect a bias in refseq annotations towards the inclusion of 3’utr ejc-harboring transcripts only if confirmed as resistant to nmd.
therefore, the transcripts we have identified in this work as nmd-target are targets of ejc-independent nmd.
these features explain just part of the differences between nmd-targets and transcripts that are not marked for degradation  <cit> .
in order to better characterize nmd-targets, we have examined many more features, none had been previously reported to be associated with nmd.
most of these features capture various aspects of the nucleotide composition around the tc.
these features are:
sejc.
a binary feature, indicating the presence/absence of a ‘shallow-ejc’, meaning an ejc that comes from an intron that resides in a 3’utr, but less than  <dig> bases downstream of the tc.
tovorf.
a binary feature, indicating whether a tuorf overlaps the main orf.
mononucleotide content.
all four nucleotide densities at the entire sequence of the 3’utr, as well as at the last  <dig> bases of the cds.
dinucleotide content.
all  <dig> dinucleotide densities at the entire sequence of the 3’utr, as well as at the last  <dig> bases of the cds.
nucleotide runs: the number of nucleotide runs  at the entire sequence of the 3’utr, as well as at the last  <dig> bases of the cds.
this number was compared to the number expected by chance based on the nucleotide densities, and the resulting χ
<dig> -statistics were also used as features.
we ended up with a list of  <dig> different features.
it is expected that this set of features contains many inter-dependencies, and we have therefore removed features that are highly correlated with other features.
this had left us with a set of  <dig> non-redundant features, of which four are binaries .
instable transcripts are characterized by high gc content around the termination codon
we next wished to find which features display a different distribution in stabilized transcripts  when compared to unaffected ones.
for the non-binary features we used the mann–whitney u-test, whereas for the binary features we used the fisher exact test.
in both cases the p-values were fdr-corrected.
none of the binary features was found to be significantly different between the two groups of transcripts , but many non-binary features did.
however, we noticed that the actual difference in the distributions between the two groups is typically very small.
in order to account for this, we computed for each non-binary feature the non-parametric common language effect size  <cit> .
for a feature measured for samples coming from two groups c
<dig> and c
<dig>  the effect size a is defined as the probability that the feature value for a random sample from c
<dig> would be higher than its value for a random sample from c
<dig>  therefore, the closer a is to  <dig> , the weaker is the effect.
in total,  <dig> features were found to be both significantly different between stabilized and unaffected transcripts  and to have high effect size  .
interestingly, all of them describe the nucleotide composition around the tc.
at the 3’utr side they include g content, several dinucleotide contents , the number of g-runs, and the χ
<dig> g-run statistic.
at the cds side they include c content, and several dinucleotide contents .
all these features point at higher gc content upstream to a ptc as well along the entire 3’utr .
differences in 3’utr g content between stabilized and unaffected transcripts are compatible with a recent study that found a preference for g in upf <dig> binding sites at the 3’utr of putative nmd-targets  <cit> .
however, binding of upf <dig> is unlikely to be the reason for this enrichment, as high g content extends also to the cds, and significant enrichment is also observed for c’s.fig.
<dig> the magnitude a −  <dig> , where a is the effect size of each non-binary feature in the non-redundant feature set.
red depicts features with effect size a −  <dig>  ≥  <dig> , blue depicts features with effect size a −  <dig>  ≤ −  <dig> , and green depicts features with effect size −  <dig>  < a −  <dig>  <  <dig> .
a nmd-targets are compared to unaffected transcripts.
b destabilized transcripts are compared to unaffected ones.
c nmd-targets are compared to destabilized transcripts
fig.
<dig> mononucleotide densities around the tc in stabilized and unaffected transcripts
we have repeated the above analysis to compare destabilized transcripts with unaffected ones, and found a total of six non-binary features that are significantly different and with high effect size .
two of them  are also the features that are most discriminative between stabilized transcripts and unaffected ones .
more generally, all these non-binary features show a similar trend of high gc content around the tc for stabilized and destabilized transcripts.
this suggests that this characterizes transcripts that are more prone to stability modifications, regardless of the nature of the specific degradation pathway.
nmd-targets are characterized by high g content around the termination codon
although stabilized and destabilized transcripts show similar nucleotide composition around the tc, we noticed that some features behave differently in the two groups .
in order to test whether any of the features discriminates directly between stabilized and destabilized transcripts, we have carried out the above analysis for a third time, now comparing destabilized transcripts with stabilized ones.
we found that four such features, all related to 3’utr g content .
this suggests that nmd-targets have higher g content around the tc, particularly along the 3’utr.
supporting our observation that g content is elevated on both sides of the tc, we note that the dinucleotide gg is significantly enriched upstream to the tc of stabilized transcripts, but not in destabilized ones.
among the binary features, the only one that significantly separates destabilized transcripts from unaffected ones is the presence of alu elements in the 3’utr , suggesting a strong depletion in destabilized transcripts.
however, this effect is attributed to the shorter 3’utr of destabilized transcripts, as the alu density along their 3’utr shows low effect size when compared to stabilized transcripts  and to unaffected ones .
our results suggest that true nmd-targets are characterized by elevated g content around the tc, especially along the 3’utr.
to test the robustness of this conclusion, we have repeated the above analyses using the strict classification criterion.
we observe qualitatively similar results , although the smaller size of the stabilized and destabilized groups reduce the statistical power of the analysis.
in addition, we wanted to account for the possible scenario that a transcript in our data set appears as lacking 3’utr ejc in refseq, and that another transcript that is indistinguishable by the bric-seq setup does harbor 3’utr ejc but is absent in refseq.
to this end, we have examined all transcripts against the more permissive ensembl database, and found that  <dig> out of the  <dig> nmd-targets  do not have any isoform with 3’utr ejcs in ensembl.
repeating the analysis above for the reduced set of these  <dig> nmd-targets led to qualitatively similar results .
more rna secondary structures near the tc of nmd-targets
why high g content around the tc promotes nmd?
we suggest that high g content around the tc increases the likelihood of the formation of secondary structures there, which hinders normal translation termination.
several examinations of codon usage bias towards the gene end had found a decrease in g in normally translated transcripts.
early works on e.coli genes revealed that they tend to end with at-rich codons  <cit> .
this was shortly followed with identical observations in b. subtilis  <cit>  and yeast  <cit> .
it has been previously shown that the presence of mrna secondary structures can stall the ribosome during the elongation phase  <cit> , or even lead to translation abortion  <cit> .
impeding the kinetics of translation termination may lead to a near cognate trna recognizing the termination codon as a sense codon and reads through it  <cit> .
by this hypothesis, high g content near the gene end promotes the formation of rna secondary structures, which interfere with normal translation termination and activates nmd.
enrichment in g along the 3’utr even in regions that are far from the tc may have other reasons, which will be detailed in the discussion.
we therefore wanted to test whether stabilized transcripts are associated with increased stability of secondary structures around the tc.
for this, we used rnafold  <cit>  to compute the minimum free energy  around the tc.
since mfe is negative, we used its absolute value so that higher values are indicative of more stable secondary structures.
we found that mfe absolute values are significantly different between stabilized, unaffected, and destabilized transcripts upstream to the tc  and along the 3’utr .
we then used dunn’s test to find which transcript groups significantly differ from each other , and computed their corresponding nonparametric common language effect size.
we found that stabilized transcripts have significantly higher mfe absolute value than unaffected transcripts both upstream to the tc and along the 3’utr, whereas destabilized transcripts only differ from unaffected ones upstream to the tc.
this result fully agrees with our inference from nucleotide composition, and suggests that nmd-targets harbor a higher density of stable rna secondary structures around the tc and, especially, along the 3’utr.table  <dig> comparison of mfe absolute values between stabilized, unaffected, and destabilized transcripts, both upstream to the tc and along the 3’utr
stabilized vs.
0
0
stabilized vs.
0
0
destabilized vs.
0
0
h
<dig> represents the null hypothesis of no difference between the two transcript groups.
effect size is defined as |a −  <dig> |, where a is the common language effect size
in order to make the above analysis more quantitative, we wished to use a measure for the level of stability change of a transcript upon upf <dig> knockdown, and to test whether it is correlated with the absolute value of mfe.
to this end, we defined a measure for the stability change of a transcript , and used it to divide the stabilized transcripts into three groups .
we used kruskal-wallis test to check whether these groups are characterized by different absolute values of mfe, and found that while mfe absolute values upstream to the tc are not related to the strength of the stability change, those along the 3’utr certainly do .
moreover, transcripts with large stability change are characterized by significantly more stable rna secondary structures along their 3’utr than transcripts with small and moderate stability change.
although transcripts with moderate stability change show statistically indistinguishable mfe absolute values than those of small stability change, they are nevertheless slightly higher .
similarly, while the three groups of transcripts show statistically identical mfe absolute values upstream of the tc, there is still a clear trend of higher values related to larger stability change.
these results suggest that higher densities of stable rna secondary structures around the tc and particularly along the 3’utr are associated with higher stability change of the transcript following upf <dig> knockdown.table  <dig> comparison of mfe absolute values between transcripts with small, moderate, and large stability change, both upstream to the tc and along the 3’utr
large vs. small
0
0
large vs. moderate
0
moderate vs. small
0
0
h
<dig> represents the null hypothesis of no difference between the two transcript groups.
effect size is defined as |a −  <dig> |, where a is the common language effect size
validating the model on independent nmd-target set
to provide further support to our model, we had examined an independent set of experimentally validated ejc-independent nmd-targets.
for this, we used a recent experiment by schmidt et al.
<cit> , who identified endogenous transcripts that are cleaved by smg <dig> using parallel analysis of rna ends .
these transcripts are potential nmd-targets, as cleavage by smg <dig> is the initial step in nmd  <cit> .
overall, schmidt et al.
identified  <dig> transcripts that are potential nmd-targets.
of these, we identified  <dig> transcripts that lack 3’utr ejcs and are therefore potential ejc-independent nmd-target .
we have used this set of transcripts as our validation set, as well as a subset of  <dig> transcripts that are also up-regulated upon knocking down upf <dig>  smg <dig>  or both.
we compared these ejc-independent nmd-targets to a list of nmd-insensitive transcripts that we prepared by taking the canonical transcripts  of all human genes , excluding genes that are reported as nmd-targets by tani et al.
or by schmidt et al.
.
analyzing this set of ejc-independent nmd-targets, we noticed that they tend to be short, both in the cds and in the 3’utr , possibly suggesting previously unnoticed bias in the experiment.
due to their short 3’utr, they are depleted in nucleotide runs, reflected by low values of the χ
2-statistics.
reassuringly, this set of nmd-targets show very similar nucleotide composition to the one we observed for our set of nmd-targets, namely higher c and particularly g contents upstream to the tc and throughout the 3’utr .fig.
<dig> the magnitude a −  <dig> , where a is the effect size of all non-binary features in the non-redundant feature set.
red depicts features with effect size a −  <dig>  ≥  <dig> , blue depicts features with effect size a −  <dig>  ≤ −  <dig> , and green depicts features with effect size −  <dig>  < a −  <dig>  <  <dig> .
a full validation set of  <dig> ecj-independent nmd-targets compared to nmd-insensitive transcripts.
b validation set of  <dig> ejc-independent nmd-targets that are also up-regulated upon smg6/upf <dig> knockout compared to nmd-insensitive transcripts
discussion
one of the challenging tasks in studying nmd is the genome-wide identification of its targets.
most studies equate nmd-targets with transcripts whose expression level increases following upf <dig> knockdown .
however, the complex networks of protein and gene interactions lead to many indirect dependencies on upf <dig> levels.
as a result, transcripts that increase their expression level following upf <dig> knockdown are not necessarily nmd-targets, and transcripts that are nmd-targets may seem insensitive to upf <dig> knockdown  <cit> .
in addition, expression levels show high variability among cell populations and across cell types, and the efficiency of nmd varies between cells and under different physiological conditions .
as a result, different studies identified only partially overlapping sets of nmd-targets.
in order to reduce the fraction of misidentified transcripts, we have defined nmd-targets as transcripts whose mrna half-life increases following upf <dig> knockdown.
while it is still plausible that some nmd-targets would be misidentified by this criterion, it is nonetheless a more direct approach, and less prone to secondary effects  <cit> .
we have used tani et al.
data of mrna half-life alteration following upf <dig> knockdown  <cit> , and developed our own technique to identify nmd-targets.
it is important to reiterate what these targets really are.
first, nmd is known to operate through ejc-dependent and ejc-independent pathways  <cit> .
in order to be compatible with tani et al.’s transcript nomenclature, we have used a transcript database  that is heavily depleted with 3’utr ejc-bearing transcripts  <cit> .
in total, among the  <dig> transcripts that could be classified, only  <dig> harbor 3’utr ejcs.
hence, the current study focuses on 3’utr ejc-independent nmd-targets.
second, upf <dig> is not critical only to nmd, but also to staufen-mediated mrna decay   <cit> .
so what we call here nmd-targets are in fact transcripts that are degraded by either smd or 3’utr ejc-independent nmd.
we assume, however, that nmd-targets form the majority of these degraded transcripts.
in addition to stabilized  and unaffected  transcripts, we have identified a third group of transcripts that are destabilized upon upf <dig> knockdown.
such transcripts are instable in the absence of nmd, possibly suggesting that their stability is maintained through the suppression of certain nmd-targets.
it remains to be studied what is the mechanism that regulates the stability of this unique group of transcripts.
we found that ejc-independent nmd-targets have high g content upstream to their tc and throughout their 3’utr.
the most discriminative feature of nmd-targets is 3’utr g content .
while destabilized transcripts also display high 3’utr g content , it is still significantly lower than that of true nmd-targets .
this enrichment in g is also observed upstream to the tc of nmd-targets, although to a lesser extent than the enrichment along the 3’utr.
for example, the dinucleotide gg is enriched upstream of the tc in nmd-targets , but not in destabilized transcripts .
this result was validated using an independent set of nmd-targets .
we showed that the enrichment of g nucleotides in nmd-targets results in higher density of stable rna secondary structures around the tc , and that nmd efficiency increases with the density of these secondary structures .
while the power of 3’utr mfe to discriminate stabilized transcripts from unaffected ones  is lower than that of 3’utr g content , the density of 3’utr mfe values  has a comparable discriminative power as the 3’utr g content .
the effect of 3’utr g content on nmd is probably not fully explained by the formation of secondary structures, as they may form quite far from the tc.
there are other arguments suggesting that upf <dig> activity is directly related to the g content along the 3’utr.
first, it was observed that upf <dig> binding sites have high g content  <cit> .
second, upf <dig> exhibits a helicase activity, and helicases are known to pauses at g-rich regions  <cit> .
therefore, hurt et al.
asserted that g-runs along the 3’utr of nmd-targets block upf <dig> scanning and lead to its accumulation at the 3’utr  <cit> .
however, this explanation may be challenged.
first, it was shown that upf <dig> binding is not an indicator of nmd activity, as both hyper-phosphorylated and hypo-phosphorylated upf <dig> bind to the same sites, but it is only the hyper-phosphorylated upf <dig> that is relevant to nmd  <cit> .
second, the claim that upf <dig> binding sites are g-rich was not replicated in a recent study by zund et al.
<cit> .
we believe that an alternative model is also plausible, by which high g content along the 3’utr increases the propensity of forming secondary structures, which promotes the recruitment of various rna-binding proteins that may interfere in different ways with the translation termination complex.
we believe that further studies should be conducted in order to determine the correct explanation.
long 3’utrs are believed to be enriched in nmd-targets  <cit> , presumably because they extend the physical distance between the tc and the poly tail.
it is therefore surprising that in our data 3’utr length is similar between nmd-targets and nmd-insensitive transcripts .
even more striking is our analysis of schmidt et al.’s ejc-independent nmd-targets, which showed that they have significantly shorter 3’utrs .
these observations suggest that ejc-independent nmd is independent of the 3’utr length.
similarly, none of the other known nmd-triggering features  are enriched in our set of nmd-targets .
we therefore conclude that it is the nucleotide composition around the tc that mainly drives ejc-independent nmd.
while several other studies identified sets of putative nmd-targets, our focus on reliable ejc-independent nmd-targets made most of them irrelevant for the current study.
for example, andersen et al.
<cit>  identified nmd-targets by knocking down smg <dig>  but all his identified nmd-targets harbor 3’utr ejcs.
also, we avoided nmd-target sets that relied on expression data.
the data produced by schmidt et al.
best suits our criteria.
first, it avoids the use of expression to define nmd-targets.
instead, they were defined by examining smg <dig> cleavage.
second, it is not exclusively made up of 3’utr ejc-containing transcripts, and thus includes sufficient number of targets of the ejc-independent nmd pathway.
CONCLUSIONS
we found that transcript features traditionally associated with nmd, such as the length of the 3’utr and the presence of tuorf, do not characterize ejc-independent nmd targets.
rather, ejc-independent nmd is likely triggered following certain nucleotide composition features that lead to elevated levels of rna secondary structures.
the main feature is high density of g nucleotides upstream of the tc, and even more so along the 3’utr.
we propose that stable rna secondary structures that are formed around the tc or along the 3’utr interfere with normal translation termination leading to transcript instability and to the activation of nmd.

 nonsense-mediated decay is a major eukaryotic surveillance mechanism that targets for degradation transcripts that harbor a premature termination codon. a primary role of this mechanism is in posttranscriptional quality control, preventing the formation of truncated proteins that are potentially detrimental to the cell. in addition, nmd is often combined with alternative splicing to form an important regulatory program of gene expression. altogether, nmd is central to maintaining normal cellular activity, and disruption of its proper function is estimated to be associated with about one third of inherited genetic disorders, as well as with many forms of cancer. however, the decision-making process at the basis of nmd, determining whether a termination codon is premature or not, is still not fully understood. it is generally believed that the context of the tc within the transcript determines whether it is normal or premature, but the nature of the contextual signals is unclear. the leading hypothesis in mammals is that a tc is recognized as premature if the pre-mrna harbors an intron more than 50– bases downstream of the tc. the splicing reaction usually leaves traces in the form of a protein complex called the exon junction complex, which is deposited upon the mrna ~20– bases upstream of the exon-exon junction. it is believed that the nmd-triggering feature is the presence of ejcs in the 3’ untranslated region at the time of translation termination. the underlying mechanistic explanation is that the first ribosome that scans the transcript clears ejcs from the coding sequence, but is unable to clear ejcs that are downstream enough to the tc. this model explains why the presence of a translated upstream open reading frame is also known to trigger nmd, as the downstream cds likely harbors ejcs which behave, in the context of the tuorf, as 3’utr ejcs. in fact, it was suggested that many of the nmd-targets result from short tuorfs, leading to an enrichment of transcripts with short cds among nmd-targets, explaining why short cds was also reported as nmd-triggering feature in yeast and drosophila. despite overall support, this model of ptc recognition is incomplete, as in various cases it was reported to be violated: some transcripts trigger nmd despite of not harboring 3’utr ejcs _enref_ while others are known to evade nmd despite of presenting 3’utr ejcs or tuorfs. these exceptions brought about the idea that nmd works through several different pathways, and led to a distinction between ejc-dependent nmd – which is the nmd that degrades transcripts that harbor 3’utr ejcs – and ejc-independent nmd, or failsafe nmd – which is the nmd that works on transcripts that lack 3’utr ejcs. an alternative model, known as the faux-utr model, suggests that a normal transcript is characterized by certain signals in its 3’utr, and that disruption of these signals triggers nmd _enref_ substantial evidence supports the notion that one of the central signals that marks a transcript as normal is physical proximity between the tc and the poly binding protein pabpc, and that a tc is tagged as premature by the lack of such physical proximity. presumably, in a normal transcript the translation termination factor erf binds pabpc whereas anything that promotes nmd, like an ejc situated downstream of the tc, competes with the pabpc and allows erf3-bound up-frameshift protein to trigger nmd. this model not only clarifies the part that 3’utr ejcs play in nmd, but also explains why very long 3’utrs are also known to trigger nmd, as the tc and pabpc are separated by a very long physical distance. however, as for other nmd-triggering features, a long 3’utr does not always result in nmd. these models assume that upf on which nmd critically depends, works through specific binding with the translation termination complex _enref_ further evidence, however, points to the possibility that upf is not only recruited by the terminating ribosome, but is associated directly with the mrna, and is even thought to have helicase activity allowing it to slide along it. it is still debated whether this association is translation-dependent or independent, but all agree that upf is displaced from the cds during translation, leading to its enrichment along the 3’utr. opinions also diverge regarding the question whether nmd-targets have excess of upf molecules bound to their 3’utr. some report that upf density along the 3’utr is higher in nmd-targets _enref_ others do not see significant difference in upf densities between nmd-targets and nmd-insensitive transcripts, claiming that it is not the upf density that matters for nmd, but rather the density of the activated upf. nmd may be directly linked to the efficiency of translation termination. it was proposed that unfavorable context of the tc may affect the recruitment of terminating factors, leading to changes in the kinetics of translation termination, and eventually to ribosome stalling at the tc that is thought to be the event triggering nmd. combined with the presence of 3’utr-bound upf this led to the suggestion that prolonged translation termination allows more upf molecules to bind to the 3’utr and become activated via regulated phosphorylation, resulting in activation of nmd _enref_ one of the biggest challenges in understanding ptc recognition stems from the observation that neither of the proposed nmd-triggering features characterizes a large fraction of nmd-targets. even when all features are considered together, it was estimated that they characterize roughly 30% of nmd-targets. this suggests that there may be many more yet unidentified features that are important for ptc recognition. surprisingly, the nucleotide composition around the tc has not received much attention. hurt et al. found that upf tends to bind 3’utrs in g-rich regions, possibly because of pausing of upf scanning in these regions. the same authors also reported that nmd-insensitive transcripts with long 3’utrs are enriched in poly stretches, perhaps because those stretches recruit pabpc thus marking the tc as normal. studies on codon usage bias show that codons near the tc tend to be at-rich, conceivably as a result of selection against the formation of rna secondary structures near the tc that may interfere with the recruitments of release factors and lead to improper translation termination. yet, the nucleotide composition of the sequence upstream the stop codon in the coding region was never regarded as a potential nmd triggering feature. here, we wished to identify additional features that are involved in the decision-making process of nmd. for this, we have used available genome-wide data on rna stability following upf knockdown, and identified transcripts that are targets for degradation. specifically, we focused on transcripts that lack 3’utr ejcs, and therefore our targets degrade through the ejc-independent nmd. for each transcript we have computed a long list of features, and tested their relevance to degradation by evaluating whether any of them shows a unique distribution within nmd-targets, when compared to non-targets. the features we have tested include all those that had been previously described as associated with nmd, as well as many new ones, mainly measuring nucleotide composition around the tc. we found that nmd-targets are characterized by high g content around the tc and throughout the 3’utr. however, we noticed a group of transcripts that are not nmd-targets and yet demonstrate the same, albeit weaker, pattern of enrichment. interestingly, these transcripts are not nmd-insensitive either, but show decreased mrna stability following upf knockdown, and are thus called ‘destabilized’. these transcripts might show instability due to factors other than nmd. surprisingly, the 3’utr length is similar for nmd-targets and nmd-insensitive transcripts, but it is substantially shorter for the destabilized transcripts. based on these observations we suggest that high g content around the tc and throughout the 3’utr leads to aberrant translation termination that triggers nmd. we show that this model has a far greater discriminative power than any other ejc-independent nmd-triggering feature known today.  identifying nmd-targets the ability to identify transcripts that are nmd-targets is a critical starting point in any analysis of nmd-triggering features. nmd-targets are traditionally identified as transcripts whose characteristic behavior is affected by knocking-down a core nmd component. the most common approach is to mark as nmd-targets transcripts whose expression level is increased following upf1-knockdown. however, this approach is expected to have a high rate of false positives due to transcripts that are secondarily up-regulated. moreover, the false negative rate may also be high due to genuine nmd-targets whose expression level is stabilized due to feedback regulation. recently, a different approach was introduced by tani et al., who used 5’-bromo-uridine immunoprecipitation chase-deep sequencing at four time points to measure genome-wide mrna stability. they compared control to upf1-knockdown samples, and marked as nmd-targets transcripts whose half-life was significantly increased following upf1-knockdown. while this approach may still suffer from secondary effects, it is nonetheless a more direct measure of rna stability. the current analysis starts with the raw data of tani et al.. we have followed parts of their analysis pipeline, but developed a novel method to detect nmd-targets. tani et al. estimated the decay rate of each transcript by modeling its mrna abundance decay curve. however, decay curves of nmd-targets often deviate from a simple exponential model, reflecting a mixtures of currently degrading transcripts with older transcripts that had escaped degradation. given that tani et al. only used four time points, their half-life estimates – especially when the decay is not a simple exponent– may have large standard errors. to circumvent this, we did not attempt to estimate decay rates, and instead developed a test that compares the shapes of the mrna abundance decay curves and marks as nmd-targets transcripts whose decay curve following upf1-knockdown is significantly above their control decay curve. this group of transcripts is denoted ‘stabilized’.fig. schematic description of the transcript classification procedure. red denotes mrna abundance decay curve in the control experiments, blue denotes mrna abundance decay curve following upf knockdown. a stabilized transcripts. b destabilized transcripts. c nmd-insensitive transcripts theoretically, knocking down upf should result in stability increase of nmd-targets, and stability stasis for nmd-insensitive transcripts. however, by using the technique above in the reverse direction we found transcripts whose stability decreases upon upf1-knockdown, namely, their mrna abundance decay curves following upf1-knockdown are significantly below their control decay curves. such transcripts are unexpected by current theory, and are denoted here ‘destabilized’. transcripts whose decay curves did not change significantly following upf1-knockdown were marked as ‘unaffected’, and represent transcripts that are not targeted for degradation by nmd. finally, transcripts whose decay curves showed inconsistent behavior were left unclassified and were excluded from further analysis. we have re-analyzed the data of tani et al. and obtained mrna abundance information for a total of transcripts. from these, we removed whose mrna abundance seemed to increase with time. of the remaining transcripts, were labeled as stabilized, were labeled as destabilized, were labeled as unaffected, and were left unclassified. in total, we could classify transcripts for which we carried out all subsequent analyses. to test the robustness of our results, we have repeated the process using a different, stricter, classification criterion totaling in classifiable transcripts. generating a set of nmd-triggering features for each of the classified transcripts we computed the following features, which were all previously suggested as characterizing nmd-targets in mammals: tuorf. a binary feature, indicating the presence/absence of translated uorfs strictly within the 5’utr. 3’utr length, in bases. 3’utr a and g content. a and g density at the entire sequence of the 3’utr. cds length, in bases. alu. a binary feature, indicating the presence/absence of alu elements in the 3’utr. these elements are not directly associated with nmd, but are nevertheless thought to be involved in mrna degradation through staufen-mediated decay. alu density. the density of alu elements within the 3’utr. we did not consider 3’utr ejc presence/absence as a feature, because the transcript list used by tani et al. was based on the refseq database, which normally excludes transcripts that harbor 3’utr ejcs. indeed, our data include only transcripts with 3’utr ejcs. ten of them were classified as unaffected, whereas the other four were left unclassified. this may reflect a bias in refseq annotations towards the inclusion of 3’utr ejc-harboring transcripts only if confirmed as resistant to nmd. therefore, the transcripts we have identified in this work as nmd-target are targets of ejc-independent nmd. these features explain just part of the differences between nmd-targets and transcripts that are not marked for degradation. in order to better characterize nmd-targets, we have examined many more features, none had been previously reported to be associated with nmd. most of these features capture various aspects of the nucleotide composition around the tc. these features are: sejc. a binary feature, indicating the presence/absence of a ‘shallow-ejc’, meaning an ejc that comes from an intron that resides in a 3’utr, but less than bases downstream of the tc. tovorf. a binary feature, indicating whether a tuorf overlaps the main orf. mononucleotide content. all four nucleotide densities at the entire sequence of the 3’utr, as well as at the last bases of the cds. dinucleotide content. all dinucleotide densities at the entire sequence of the 3’utr, as well as at the last bases of the cds. nucleotide runs: the number of nucleotide runs at the entire sequence of the 3’utr, as well as at the last bases of the cds. this number was compared to the number expected by chance based on the nucleotide densities, and the resulting χ -statistics were also used as features. we ended up with a list of different features. it is expected that this set of features contains many inter-dependencies, and we have therefore removed features that are highly correlated with other features. this had left us with a set of non-redundant features, of which four are binaries. instable transcripts are characterized by high gc content around the termination codon we next wished to find which features display a different distribution in stabilized transcripts when compared to unaffected ones. for the non-binary features we used the mann–whitney u-test, whereas for the binary features we used the fisher exact test. in both cases the p-values were fdr-corrected. none of the binary features was found to be significantly different between the two groups of transcripts, but many non-binary features did. however, we noticed that the actual difference in the distributions between the two groups is typically very small. in order to account for this, we computed for each non-binary feature the non-parametric common language effect size. for a feature measured for samples coming from two groups c and c the effect size a is defined as the probability that the feature value for a random sample from c would be higher than its value for a random sample from c therefore, the closer a is to, the weaker is the effect. in total, features were found to be both significantly different between stabilized and unaffected transcripts and to have high effect size. interestingly, all of them describe the nucleotide composition around the tc. at the 3’utr side they include g content, several dinucleotide contents, the number of g-runs, and the χ g-run statistic. at the cds side they include c content, and several dinucleotide contents. all these features point at higher gc content upstream to a ptc as well along the entire 3’utr. differences in 3’utr g content between stabilized and unaffected transcripts are compatible with a recent study that found a preference for g in upf binding sites at the 3’utr of putative nmd-targets. however, binding of upf is unlikely to be the reason for this enrichment, as high g content extends also to the cds, and significant enrichment is also observed for c’s.fig. the magnitude a − , where a is the effect size of each non-binary feature in the non-redundant feature set. red depicts features with effect size a −   ≥ , blue depicts features with effect size a −   ≤ − , and green depicts features with effect size −   < a −   < . a nmd-targets are compared to unaffected transcripts. b destabilized transcripts are compared to unaffected ones. c nmd-targets are compared to destabilized transcripts fig. mononucleotide densities around the tc in stabilized and unaffected transcripts we have repeated the above analysis to compare destabilized transcripts with unaffected ones, and found a total of six non-binary features that are significantly different and with high effect size. two of them are also the features that are most discriminative between stabilized transcripts and unaffected ones. more generally, all these non-binary features show a similar trend of high gc content around the tc for stabilized and destabilized transcripts. this suggests that this characterizes transcripts that are more prone to stability modifications, regardless of the nature of the specific degradation pathway. nmd-targets are characterized by high g content around the termination codon although stabilized and destabilized transcripts show similar nucleotide composition around the tc, we noticed that some features behave differently in the two groups. in order to test whether any of the features discriminates directly between stabilized and destabilized transcripts, we have carried out the above analysis for a third time, now comparing destabilized transcripts with stabilized ones. we found that four such features, all related to 3’utr g content. this suggests that nmd-targets have higher g content around the tc, particularly along the 3’utr. supporting our observation that g content is elevated on both sides of the tc, we note that the dinucleotide gg is significantly enriched upstream to the tc of stabilized transcripts, but not in destabilized ones. among the binary features, the only one that significantly separates destabilized transcripts from unaffected ones is the presence of alu elements in the 3’utr, suggesting a strong depletion in destabilized transcripts. however, this effect is attributed to the shorter 3’utr of destabilized transcripts, as the alu density along their 3’utr shows low effect size when compared to stabilized transcripts and to unaffected ones. our results suggest that true nmd-targets are characterized by elevated g content around the tc, especially along the 3’utr. to test the robustness of this conclusion, we have repeated the above analyses using the strict classification criterion. we observe qualitatively similar results, although the smaller size of the stabilized and destabilized groups reduce the statistical power of the analysis. in addition, we wanted to account for the possible scenario that a transcript in our data set appears as lacking 3’utr ejc in refseq, and that another transcript that is indistinguishable by the bric-seq setup does harbor 3’utr ejc but is absent in refseq. to this end, we have examined all transcripts against the more permissive ensembl database, and found that out of the nmd-targets do not have any isoform with 3’utr ejcs in ensembl. repeating the analysis above for the reduced set of these nmd-targets led to qualitatively similar results. more rna secondary structures near the tc of nmd-targets why high g content around the tc promotes nmd? we suggest that high g content around the tc increases the likelihood of the formation of secondary structures there, which hinders normal translation termination. several examinations of codon usage bias towards the gene end had found a decrease in g in normally translated transcripts. early works on e.coli genes revealed that they tend to end with at-rich codons. this was shortly followed with identical observations in b. subtilis and yeast. it has been previously shown that the presence of mrna secondary structures can stall the ribosome during the elongation phase, or even lead to translation abortion. impeding the kinetics of translation termination may lead to a near cognate trna recognizing the termination codon as a sense codon and reads through it. by this hypothesis, high g content near the gene end promotes the formation of rna secondary structures, which interfere with normal translation termination and activates nmd. enrichment in g along the 3’utr even in regions that are far from the tc may have other reasons, which will be detailed in the discussion. we therefore wanted to test whether stabilized transcripts are associated with increased stability of secondary structures around the tc. for this, we used rnafold to compute the minimum free energy around the tc. since mfe is negative, we used its absolute value so that higher values are indicative of more stable secondary structures. we found that mfe absolute values are significantly different between stabilized, unaffected, and destabilized transcripts upstream to the tc and along the 3’utr. we then used dunn’s test to find which transcript groups significantly differ from each other, and computed their corresponding nonparametric common language effect size. we found that stabilized transcripts have significantly higher mfe absolute value than unaffected transcripts both upstream to the tc and along the 3’utr, whereas destabilized transcripts only differ from unaffected ones upstream to the tc. this result fully agrees with our inference from nucleotide composition, and suggests that nmd-targets harbor a higher density of stable rna secondary structures around the tc and, especially, along the 3’utr.table comparison of mfe absolute values between stabilized, unaffected, and destabilized transcripts, both upstream to the tc and along the 3’utr stabilized vs. 0 0 stabilized vs. 0 0 destabilized vs. 0 0 h represents the null hypothesis of no difference between the two transcript groups. effect size is defined as |a −  |, where a is the common language effect size in order to make the above analysis more quantitative, we wished to use a measure for the level of stability change of a transcript upon upf knockdown, and to test whether it is correlated with the absolute value of mfe. to this end, we defined a measure for the stability change of a transcript, and used it to divide the stabilized transcripts into three groups. we used kruskal-wallis test to check whether these groups are characterized by different absolute values of mfe, and found that while mfe absolute values upstream to the tc are not related to the strength of the stability change, those along the 3’utr certainly do. moreover, transcripts with large stability change are characterized by significantly more stable rna secondary structures along their 3’utr than transcripts with small and moderate stability change. although transcripts with moderate stability change show statistically indistinguishable mfe absolute values than those of small stability change, they are nevertheless slightly higher. similarly, while the three groups of transcripts show statistically identical mfe absolute values upstream of the tc, there is still a clear trend of higher values related to larger stability change. these results suggest that higher densities of stable rna secondary structures around the tc and particularly along the 3’utr are associated with higher stability change of the transcript following upf knockdown.table comparison of mfe absolute values between transcripts with small, moderate, and large stability change, both upstream to the tc and along the 3’utr large vs. small 0 0 large vs. moderate 0 moderate vs. small 0 0 h represents the null hypothesis of no difference between the two transcript groups. effect size is defined as |a −  |, where a is the common language effect size validating the model on independent nmd-target set to provide further support to our model, we had examined an independent set of experimentally validated ejc-independent nmd-targets. for this, we used a recent experiment by schmidt et al., who identified endogenous transcripts that are cleaved by smg using parallel analysis of rna ends. these transcripts are potential nmd-targets, as cleavage by smg is the initial step in nmd. overall, schmidt et al. identified transcripts that are potential nmd-targets. of these, we identified transcripts that lack 3’utr ejcs and are therefore potential ejc-independent nmd-target. we have used this set of transcripts as our validation set, as well as a subset of transcripts that are also up-regulated upon knocking down upf smg or both. we compared these ejc-independent nmd-targets to a list of nmd-insensitive transcripts that we prepared by taking the canonical transcripts of all human genes, excluding genes that are reported as nmd-targets by tani et al. or by schmidt et al.. analyzing this set of ejc-independent nmd-targets, we noticed that they tend to be short, both in the cds and in the 3’utr, possibly suggesting previously unnoticed bias in the experiment. due to their short 3’utr, they are depleted in nucleotide runs, reflected by low values of the χ 2-statistics. reassuringly, this set of nmd-targets show very similar nucleotide composition to the one we observed for our set of nmd-targets, namely higher c and particularly g contents upstream to the tc and throughout the 3’utr.fig. the magnitude a − , where a is the effect size of all non-binary features in the non-redundant feature set. red depicts features with effect size a −   ≥ , blue depicts features with effect size a −   ≤ − , and green depicts features with effect size −   < a −   < . a full validation set of ecj-independent nmd-targets compared to nmd-insensitive transcripts. b validation set of ejc-independent nmd-targets that are also up-regulated upon smg6/upf knockout compared to nmd-insensitive transcripts discussion one of the challenging tasks in studying nmd is the genome-wide identification of its targets. most studies equate nmd-targets with transcripts whose expression level increases following upf knockdown. however, the complex networks of protein and gene interactions lead to many indirect dependencies on upf levels. as a result, transcripts that increase their expression level following upf knockdown are not necessarily nmd-targets, and transcripts that are nmd-targets may seem insensitive to upf knockdown. in addition, expression levels show high variability among cell populations and across cell types, and the efficiency of nmd varies between cells and under different physiological conditions. as a result, different studies identified only partially overlapping sets of nmd-targets. in order to reduce the fraction of misidentified transcripts, we have defined nmd-targets as transcripts whose mrna half-life increases following upf knockdown. while it is still plausible that some nmd-targets would be misidentified by this criterion, it is nonetheless a more direct approach, and less prone to secondary effects. we have used tani et al. data of mrna half-life alteration following upf knockdown, and developed our own technique to identify nmd-targets. it is important to reiterate what these targets really are. first, nmd is known to operate through ejc-dependent and ejc-independent pathways. in order to be compatible with tani et al.’s transcript nomenclature, we have used a transcript database that is heavily depleted with 3’utr ejc-bearing transcripts. in total, among the transcripts that could be classified, only harbor 3’utr ejcs. hence, the current study focuses on 3’utr ejc-independent nmd-targets. second, upf is not critical only to nmd, but also to staufen-mediated mrna decay. so what we call here nmd-targets are in fact transcripts that are degraded by either smd or 3’utr ejc-independent nmd. we assume, however, that nmd-targets form the majority of these degraded transcripts. in addition to stabilized and unaffected transcripts, we have identified a third group of transcripts that are destabilized upon upf knockdown. such transcripts are instable in the absence of nmd, possibly suggesting that their stability is maintained through the suppression of certain nmd-targets. it remains to be studied what is the mechanism that regulates the stability of this unique group of transcripts. we found that ejc-independent nmd-targets have high g content upstream to their tc and throughout their 3’utr. the most discriminative feature of nmd-targets is 3’utr g content. while destabilized transcripts also display high 3’utr g content, it is still significantly lower than that of true nmd-targets. this enrichment in g is also observed upstream to the tc of nmd-targets, although to a lesser extent than the enrichment along the 3’utr. for example, the dinucleotide gg is enriched upstream of the tc in nmd-targets, but not in destabilized transcripts. this result was validated using an independent set of nmd-targets. we showed that the enrichment of g nucleotides in nmd-targets results in higher density of stable rna secondary structures around the tc, and that nmd efficiency increases with the density of these secondary structures. while the power of 3’utr mfe to discriminate stabilized transcripts from unaffected ones is lower than that of 3’utr g content, the density of 3’utr mfe values has a comparable discriminative power as the 3’utr g content. the effect of 3’utr g content on nmd is probably not fully explained by the formation of secondary structures, as they may form quite far from the tc. there are other arguments suggesting that upf activity is directly related to the g content along the 3’utr. first, it was observed that upf binding sites have high g content. second, upf exhibits a helicase activity, and helicases are known to pauses at g-rich regions. therefore, hurt et al. asserted that g-runs along the 3’utr of nmd-targets block upf scanning and lead to its accumulation at the 3’utr. however, this explanation may be challenged. first, it was shown that upf binding is not an indicator of nmd activity, as both hyper-phosphorylated and hypo-phosphorylated upf bind to the same sites, but it is only the hyper-phosphorylated upf that is relevant to nmd. second, the claim that upf binding sites are g-rich was not replicated in a recent study by zund et al.. we believe that an alternative model is also plausible, by which high g content along the 3’utr increases the propensity of forming secondary structures, which promotes the recruitment of various rna-binding proteins that may interfere in different ways with the translation termination complex. we believe that further studies should be conducted in order to determine the correct explanation. long 3’utrs are believed to be enriched in nmd-targets, presumably because they extend the physical distance between the tc and the poly tail. it is therefore surprising that in our data 3’utr length is similar between nmd-targets and nmd-insensitive transcripts. even more striking is our analysis of schmidt et al.’s ejc-independent nmd-targets, which showed that they have significantly shorter 3’utrs. these observations suggest that ejc-independent nmd is independent of the 3’utr length. similarly, none of the other known nmd-triggering features are enriched in our set of nmd-targets. we therefore conclude that it is the nucleotide composition around the tc that mainly drives ejc-independent nmd. while several other studies identified sets of putative nmd-targets, our focus on reliable ejc-independent nmd-targets made most of them irrelevant for the current study. for example, andersen et al. identified nmd-targets by knocking down smg but all his identified nmd-targets harbor 3’utr ejcs. also, we avoided nmd-target sets that relied on expression data. the data produced by schmidt et al. best suits our criteria. first, it avoids the use of expression to define nmd-targets. instead, they were defined by examining smg cleavage. second, it is not exclusively made up of 3’utr ejc-containing transcripts, and thus includes sufficient number of targets of the ejc-independent nmd pathway.  we found that transcript features traditionally associated with nmd, such as the length of the 3’utr and the presence of tuorf, do not characterize ejc-independent nmd targets. rather, ejc-independent nmd is likely triggered following certain nucleotide composition features that lead to elevated levels of rna secondary structures. the main feature is high density of g nucleotides upstream of the tc, and even more so along the 3’utr. we propose that stable rna secondary structures that are formed around the tc or along the 3’utr interfere with normal translation termination leading to transcript instability and to the activation of nmd. 
instead, we saw that nmd-targets are strongly enriched with g nucleotides upstream of the termination codon, and even more so along their 3’utr. this is particularly true for transcripts lacking an exon junction complex along their 3’ untranslated region, and thus degrade through the ejc-independent nmd pathway. we examined many features characterizing these transcripts, and compared them to nmd-insensitive transcripts, as well as to a group of transcripts that are destabilized following nmd repression. nonsense-mediated decay ejc-independent nmdnmd-triggering featuresstop codon gc contentstop codon nucleotide compositionrna secondary structureexon junction complex transcription terminationhttp://dx.doi.org/ /501100003977israel science foundation1431/13carmel liran issue-copyright-statement© the author 2016 it is not fully understood how a termination codon is recognized as premature by the nonsense-mediated decay machinery. we suggest that high g content around the termination codon impedes translation termination as a result of mrna folding, thus triggering nmd. we also suggest that high g content in the 3’utr helps to activate nmd by allowing for the accumulation of upf or other nmd-promoting proteins, along the 3’utr.
 it is not fully understood how a termination codon is recognized as premature by the nonsense-mediated decay machinery. this is particularly true for transcripts lacking an exon junction complex along their 3’ untranslated region, and thus degrade through the ejc-independent nmd pathway.  here, we analyzed data of transcript stability change following nmd repression and identified over ejc-independent nmd-targets. we examined many features characterizing these transcripts, and compared them to nmd-insensitive transcripts, as well as to a group of transcripts that are destabilized following nmd repression.  we found that none of the known nmd-triggering features, such as the presence of upstream open reading frames, significantly characterizes ejc-independent nmd-targets. instead, we saw that nmd-targets are strongly enriched with g nucleotides upstream of the termination codon, and even more so along their 3’utr. we suggest that high g content around the termination codon impedes translation termination as a result of mrna folding, thus triggering nmd. we also suggest that high g content in the 3’utr helps to activate nmd by allowing for the accumulation of upf or other nmd-promoting proteins, along the 3’utr. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords nonsense-mediated decay ejc-independent nmdnmd-triggering featuresstop codon gc contentstop codon nucleotide compositionrna secondary structureexon junction complex transcription terminationhttp://dx.doi.org/ /501100003977israel science foundation1431/13carmel liran issue-copyright-statement© the author 2016 