Results
Fresh yeast cultures were grown from a stock source for each of the mutant strains and wild type strain used in order to isolate RNA from cells which are at the same stage of growth. The yeast were allowed to grow to mid log phase by measuring the absorbance of the cell culture at an optical density of 600nm. Mid log phase for yeast is considered to be 0.4 - 1.7 absorbance at OD of 600nm and the optical density of the three strains used in this experiment are listed in Table 2 below.
Table 2. Absorbance readings for yeast cultures used to
isolate RNA samples for the reverse transcription into cDNA.
| Yeast Culture | OD600 |
| Wild Type | 0.624 |
| Δzms1 | 0.666 |
| Δzms1/Δzms2 | 0.600 |
Total RNA was isolated from the yeast cultures of all three strains and the purity and concentration of RNA samples used for the microarray slides analyzed in this experiment are listed in Table 3. The purity is a ratio of the Absorbance at 260nm to the Abosrbance at 280nm and was obtained using Nanodrop analysis. The concentration of the RNA is given in nanograms per microliter and was also obtained by Nanodrop analysis.
Table 3. Purity and Concentrations for the RNA samples of each group whose microarray data was analyzed.
| RNA Sample | Purity (A260/A280) | Concentration (ng/μl) |
| Δzms1 (Group 1) | 2.16 | 1197.0 |
| Δzms1 (Group 4) | 1.86 | 650.2 |
| Δzms1 (Group 9) | (unable to obtain) | (unable to obtain) |
| Δzms1/Δzms2 (Group 2) | 2.19 | 536.0 |
| Wild Type (Group 1) | 2.07 | 418.0 |
To be sure that the RNA samples isolated by each group were intact and not degraded, a 1.2% agarose gel was run and stained with Ethidium Bromide (Figure 2 and 3). RNA in Lane 2, 3 and 4 of Figure 2 and Lane 1, 3 and 4 of Figure 3 contain questionably intact RNA which is either a low concentration or degraded. There also appears to be DNA contaminant in all Lanes of Figure 2 and in Lane 2 of Figure 3.

Figure 2. RNA run on a 1.2% agarose gel. Lane
1 contains
Figure 3. RNA run on a 1.2%
agarose gel. Lane 1 contains
Δzms1 RNA isolated by Group 4 and Lane 2
contains
Δzms2
wild type RNA isolated by Group 3 and Lane 2 contains
Δzms1Δzms2
RNA isolated by Group 4. Lane 3 contains
Δzms2
RNA isolated
RNA isolated by Group 3. Lane 3 contains
Δzms2
RNA isolated
by Group 5 and Lane 4 contains
Δzms1Δzms2
RNA isolated by
by Group 2 and Lane 4 contains wild type RNA isolated by the same
the same group. Lane 5 and 6 contain RNA
Δzms1
and wild type
group.
RNA respectively, isolated by Group 1.
After isolating the RNA, the selected RNA samples from Table 3 were reverse transcribed to produce cDNA which was then hybridized to microarray slides containing over 6500 yeast genes. All slides were set up to compare the mutant strain expression to wild type yeast expression. Figure 3 below is a sample portion of the microarray data which was then gridded in Scanalyze software (Blocks 7 and 8), and analyzed in Magic Tool version 2.1.
Figure 4. Portion of the microarray slide #1360722 hybridized with
Δzms1 cDNA labeled with Cy5 dye and Wild Type cDNA labeled
with Cy3 dye. Image obtained from Scanalyze software used for gridding.
Data obtained from the scan at Davidson College were presented in an Excel file. Pixel intensities for both the Cy3 (CH1I) and Cy5 (CH2I) dyes for each gene were given, along with background intensities. Plots of the Cy3 vs. Cy5 intensities are expected to show a best fit line with a slope equal to one if the data are normally distributed. However, the plots in Figure 5 show that the data from group 4 (A) appears skewed in favor of Cy3 with the plot having a best fit line with slope greater than one (m = 5.0935). Similarly, the data from group 9 (D) is also skewed in favor of Cy3 (m = 1.6755). The data from groups 2 (B) and 1 (C) both show data skewed in favor of Cy5 with slopes less than one (m = 0.2378 and m = 0.01 respectively.



Figure 5. Graphs presenting the green and red intensities from the four slides analyzed before standardization. A. Group 4 Δzms1 labeled with Cy5 vs. wild type labeled with Cy3. B. Group 2 Δzms1 labeled with Cy3 vs. wild type labeled with Cy5. C. Group 1 Δzms1 labeled with Cy5 vs. wild type labeled with Cy3. D. Group 9 Δzms1Δzms2 labeled with Cy5 vs. wild type labeled with Cy3.
In order to start analysis of the data, background intensities were subtracted off of the total intensities and any resulting negative or zero values were discarded. Then ratios of mutant to wild-type expression were obtained by dividing Cy5/Cy3 intensities. Data from each of the four microarray slides were then loaded into the software program Magic Tool for further analysis. A log transformation of log to the base two was performed to have the data within a reasonable range for interpretation (from -16 to 16). Comparison between multiple slides is only reliable when the data sets have similar distributions so in order to visualize this, boxplots were created (Figure 6A). Since the distributions were not similar between the four microarray slides being analyzed, standardization to a mean of zero and a standard deviation of one of the data was required. After standardization the data were plotted again (Figure 6B) for confirmation that the distributions were approximately similar.
A
B
Figure 6. Box plots of the ratios for each microarray slide analyzed. A. Ratio values obtained after
taking the log base 2 of the ratios and before standardization. B. Ratio values obtained after standardizing
to a mean of zero and standard deviation of one. Image obtained from Magic Tool version 2.1 software.
Genes were selected as significant if they showed consistent expression levels in all three Δzms1 trials and contrasting expression levels in the Δzms1/Δzms2 trial. Table 4 shows the twenty genes selected as significant along with the specific intensity ratios observed on each microarray slide. Table 5 shows those selected genes’ molecular function and the biological processes in which they are involved.
Table 4. Genes with consistent expression patterns within three trials of ZMS1 mutants which have expression patterns in contrast to the ZMS1 and ZMS2 double mutant. The ORF name, gene name if known and chromosome number are given. Expression patterns are displayed in intensity ratios of mutant to wild type, where a negative ratio represents down-regulation and a positive ratio represents up-regulation.
|
ORF Name |
Gene |
Chromosome |
zms1 |
zms1 |
zms1 |
zms1zms2 (Group 9) |
|
YKL197C |
PEX1 |
11 |
-0.1733 |
-0.8946 |
-0.6861 |
0.8917 |
|
YAL045C |
|
1 |
-0.1253 |
-0.5438 |
-1.1493 |
0.8551 |
|
YHL008C |
|
8 |
-0.2447 |
-0.727 |
-0.5846 |
0.4618 |
|
YHR187W |
IKI1 |
8 |
-0.2355 |
-0.9531 |
-0.1768 |
0.9295 |
|
YPL213W |
LEA1 |
16 |
-0.8305 |
-0.6529 |
-0.4104 |
0.1776 |
|
YDR416W |
SYF1 |
4 |
-0.3913 |
-0.785 |
-0.6256 |
0.6892 |
|
YLL006W |
MMM1 |
12 |
-0.7703 |
-0.843 |
-0.4177 |
0.0497 |
|
YCR037C |
PHO87 |
3 |
-0.7498 |
-0.949 |
-0.825 |
1.084 |
|
YOR259C |
RPT4 |
15 |
-2.557 |
-2.2819 |
-0.3978 |
1.1816 |
|
YDL207W |
GLE1 |
4 |
-0.04549 |
-0.8132 |
-0.6312 |
0.7983 |
|
YDR062W |
LCB2 |
4 |
-1.1133 |
-0.8945 |
-1.1985 |
0.6234 |
|
YMR174C |
PAI3 |
13 |
-0.4768 |
-0.5984 |
-0.641 |
0.728 |
|
YPR189W |
SKI3 |
16 |
1.3372 |
1.125 |
0.7516 |
-0.4816 |
|
YDR261C |
EXG2 |
4 |
1.1374 |
0.7629 |
2.9679 |
-1.3738 |
|
YJL222W |
VTH2 |
10 |
0.5403 |
0.5519 |
0.5704 |
-0.8708 |
|
YKL096W |
CWP1 |
11 |
-0.05486 |
-1.6538 |
-0.7563 |
0.5739 |
|
YDR065W |
|
4 |
-1.6 |
-0.2581 |
-0.816 |
0.4641 |
|
YLR431C |
|
12 |
-0.6461 |
0.4716 |
-0.5621 |
1.2265 |
|
YML005W |
|
13 |
-0.2272 |
-0.7056 |
-0.1924 |
0.74 |
|
YLR392C |
|
12 |
-0.4203 |
-0.6674 |
-0.8158 |
0.6066 |
Table 5. Biological processes and molecular functions of genes found to be differentially expressed between the ZMS1 mutant and the ZMS1 and ZMS2 double mutant, corresponding to the genes in Table 4.
|
ORF Name |
Gene |
Chromosome |
Biological Process |
Molecular Function |
|
YKL197C |
PEX1 |
11 |
peroxisome organization and biogenesis |
ATPase activity |
|
YAL045C |
|
1 |
biological_process unknown |
molecular_function unknown |
|
YHL008C |
|
8 |
biological_process unknown |
transporter activity |
|
YHR187W |
IKI1 |
8 |
regulation of transcription from Pol II promoter |
Pol II transcription elongation factor activity |
|
YPL213W |
LEA1 |
16 |
mRNA splicing |
pre-mRNA splicing factor activity |
|
YDR416W |
SYF1 |
4 |
mRNA splicing* |
molecular_function unknown |
|
YLL006W |
MMM1 |
12 |
mitochondrion organization and biogenesis* |
molecular_function unknown |
|
YCR037C |
PHO87 |
3 |
phosphate transport |
inorganic phosphate transporter activity |
|
YOR259C |
RPT4 |
15 |
ubiquitin-dependent protein catabolism |
ATPase activity* |
|
YDL207W |
GLE1 |
4 |
poly(A)+ mRNA-nucleus export |
molecular_function unknown |
|
YDR062W |
LCB2 |
4 |
sphingolipid biosynthesis |
serine C-palmitoyltransferase activity |
|
YMR174C |
PAI3 |
13 |
vacuolar protein catabolism |
endopeptidase inhibitor activity |
|
YPR189W |
SKI3 |
16 |
mRNA catabolism |
translation repressor activity |
|
YDR261C |
EXG2 |
4 |
biological_process unknown |
glucan 1,3-beta-glucosidase activity |
|
YJL222W |
VTH2 |
10 |
Golgi to vacuole transport |
protein signal sequence binding activity |
|
YKL096W |
CWP1 |
11 |
cell wall organization and biogenesis |
structural constituent of cell wall |
|
YDR065W |
|
4 |
biological_process unknown |
molecular_function unknown |
|
YLR431C |
|
12 |
biological_process unknown |
molecular_function unknown |
|
YML005W |
|
13 |
biological_process unknown |
molecular_function unknown |
|
YLR392C |
|
12 |
biological_process unknown |
molecular_function unknown |