04.04.2013 Views

Transcriptional Characterization of Glioma Neural Stem Cells Diva ...

Transcriptional Characterization of Glioma Neural Stem Cells Diva ...

Transcriptional Characterization of Glioma Neural Stem Cells Diva ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

5.6 Differential Is<strong>of</strong>orm Expression Methods<br />

Category 4. Is there no literature implicating the gene in any cancer?<br />

If W = total weighted index and C = cancer index, then if W > 10 (internal<br />

parameter, does not appear in fig 5.6), the gene is assigned to category one; if 1<br />

< W < 10, the gene is assigned to category two; if W = 0 and C > 10 (internal<br />

parameter, does not appear in fig 5.6) then the gene is assigned to category<br />

three; finally, if W = 0 and C = 0 then the gene is assigned to category four.<br />

5.6 Differential Is<strong>of</strong>orm Expression<br />

With the aim <strong>of</strong> establishing Tag-seq as a sensitive technique for the iden-<br />

tification <strong>of</strong> 3'UTR transcript is<strong>of</strong>orms that are differentially expressed be-<br />

tween GNS and NS cell lines, we evaluated the differential expression <strong>of</strong> 4,727<br />

genes with multiple expressed tags using two alternative methods: the non-<br />

parametric χ 2 test, and the logarithmic algorithm adapted from Morrissy et<br />

al [346]. Non-parametric methods are applied when the assumptions for para-<br />

metric methods about the underlying distribution <strong>of</strong> the data are not met,<br />

namely that the data is normally distributed or that sample size is large enough<br />

to support the assessment <strong>of</strong> the normality distribution. Whilst the normal-<br />

ity assumption is satisfied when tested against the 28,351 tags mapping over<br />

16,025 genes in the six libraries altogether, sample size drops to a range <strong>of</strong> two<br />

to 14 tags per gene when assessing differential expression over a single tran-<br />

script is<strong>of</strong>orm. With a sample size smaller than 30, the use <strong>of</strong> non-parametric<br />

methods is <strong>of</strong>ten advisable. Due to the nature <strong>of</strong> the SAGE protocol, tags map<br />

at the 3'UTR <strong>of</strong> a gene and identify potential splicing is<strong>of</strong>orms. We conducted<br />

a χ 2 test with the chisq.test function in R on all tags mapping on to the same<br />

gene <strong>of</strong> the 12,794 tags expressed at higher than 10 tpm in at least one GNS<br />

and one NS library.<br />

In the second approach, following the logarithmic ratio method adapted by<br />

Morrissy et al [346], the is<strong>of</strong>orm expression between the GNS and NS cell lines<br />

was evaluated by analyzing the expression <strong>of</strong> tags that mapped either twice on<br />

the same gene or, if more than twice, for all their pair combinations. Again,<br />

only tags with a minimum expression <strong>of</strong> 10 tpm in at least one GNS and one<br />

NS library were considered (Table 5.3). The ratio <strong>of</strong> normalised expression<br />

for each tag pair was calculated in each GNS and NS library, and multiplied<br />

by the natural logarithm <strong>of</strong> the difference between the expression <strong>of</strong> each tag<br />

in the pair, ensuring that the more highly expressed tag pairs would also be<br />

more highly ranked (equations 5.4 and 5.5). For pairs <strong>of</strong> genes with positive<br />

110

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!