Skip to main content
  • Research article
  • Open access
  • Published:

Global mRNA expression analysis in myosin II deficient strains of Saccharomyces cerevisiae reveals an impairment of cell integrity functions

Abstract

Background

The Saccharomyces cerevisiae MYO1 gene encodes the myosin II heavy chain (Myo1p), a protein required for normal cytokinesis in budding yeast. Myo1p deficiency in yeast (myo1Δ) causes a cell separation defect characterized by the formation of attached cells, yet it also causes abnormal budding patterns, formation of enlarged and elongated cells, increased osmotic sensitivity, delocalized chitin deposition, increased chitin synthesis, and hypersensitivity to the chitin synthase III inhibitor Nikkomycin Z. To determine how differential expression of genes is related to these diverse cell wall phenotypes, we analyzed the global mRNA expression profile of myo1Δ strains.

Results

Global mRNA expression profiles of myo1Δ strains and their corresponding wild type controls were obtained by hybridization to yeast oligonucleotide microarrays. Results for selected genes were confirmed by real time RT-PCR. A total of 547 differentially expressed genes (p ≤ 0.01) were identified with 263 up regulated and 284 down regulated genes in the myo1Δ strains. Gene set enrichment analysis revealed the significant over-representation of genes in the protein biosynthesis and stress response categories. The SLT2/MPK1 gene was up regulated in the microarray, and a myo1Δslt2Δ double mutant was non-viable. Overexpression of ribosomal protein genes RPL30 and RPS31 suppressed the hypersensitivity to Nikkomycin Z and increased the levels of phosphorylated Slt2p in myo1Δ strains. Increased levels of phosphorylated Slt2p were also observed in wild type strains under these conditions.

Conclusion

Following this analysis of global mRNA expression in yeast myo1Δ strains, we conclude that 547 genes were differentially regulated in myo1Δ strains and that the stress response and protein biosynthesis gene categories were coordinately regulated in this mutant. The SLT2/MPK1 gene was confirmed to be essential for myo1Δ strain viability, supporting that the up regulated stress response genes are regulated by the PKC1 cell integrity pathway. Suppression of Nikkomycin Z hypersensitivity together with Slt2p phosphorylation was caused by the overexpression of ribosomal protein genes RPL30 and RPS31. These ribosomal protein mRNAs were down regulated in the myo1Δ arrays, suggesting that down regulation of ribosomal biogenesis may affect cell integrity in myo1Δ strains.

Background

In yeast cells as in animal cells, cytokinesis requires the formation of an actomyosin contractile ring [1]. Saccharomyces cerevisiae possesses a single MYO1 gene, which encodes the myosin II heavy chain protein (Myo1p), the only conventional myosin II of the budding yeast. Myo1p is associated with F-actin during formation of the contractile apparatus, and provides the actin-activated ATPase activity driving the contraction of the actomyosin ring at cytokinesis [2, 3]. Previous studies have described the importance of Myo1p in cytokinesis in yeast by establishing that myosin II concentrates at the bud neck, contracts during the cytokinesis, and subsequently is redistributed or degraded after the completion of mitosis [2, 4].

In addition to the loss of normal cytokinesis function, myo1Δ strains present phenotypes such as the formation of attached cells, loss of axial budding pattern, cell enlargement, osmotic hypersensitivity, delocalized chitin deposition [5], increased chitin synthesis[6, 7], and hypersensitivity to Nikkomycin Z [8]. The diversity of these phenotypes, many of which can be associated with cell wall organization, suggests that myosin II may affect cellular functions other than cytokinesis. The purpose of this study was to generate a global mRNA expression profile of myo1Δ strains using yeast oligonucleotide microarrays to identify potential relationships between the multiple myo1Δ phenotypes mentioned and changes in the transcription profile of specific genes that can be functionally related to these phenotypes.

Results and Discussion

The results of independent microarray hybridization experiments conducted with two different myo1Δ knockout strains were averaged and combined for comparison with the corresponding wild-type control values. Our rationale for this experimental design was to identify deferentially regulated genes held in common between at least two different myo1Δ strain backgrounds in order to eliminate any strain specific differences and preserve strictly myo1Δ-dependent genetic phenotypes. Cells were cultured in Histidine dropout medium or CSM to maintain the culture conditions under which myo1Δ strains have been shown to require chitin synthase III expression for their survival [8] as well as to maintain equivalent growth conditions between mutant and wild-type strains, respectively.

Up regulated genes in myo1Δ

Previously published studies reported differentially expressed genes under conditions such as cell wall mutations, chemical and biochemical perturbations of the cell wall [9], diverse mutations covering 4% of the yeast genome [10], and during the cell cycle [11], employing values with p ≤ 0.05. We identified 97 genes in common with these previous studies out of a total of 2,273 genes with p ≤ 0.05 that were differentially expressed in myo1Δ strains (data not shown).

To refine this analysis we used a more stringent cutoff (p ≤ 0.01) to identify 547 differentially expressed genes in myo1Δ strains. We observed a total of 263 up regulated genes in myo1Δ strains compared to the wild type control strains. The functional categories of the up regulated genes that were represented as listed in the Saccharomyces Genome Database (SGD) were: cell organization and biogenesis (17/263), metabolism (25/263), transport (9/263), carbohydrate metabolism (23/263), stress response (17/263), protein biosynthesis (2/263), cell cycle (3/263), signal transduction (7/263), transcription (1/263), sporulation (4/263), budding (4/263), RNA processing (4/263), protein transport (6/263), protein amino acid phosphorylation (2/263), protein amino acid dephosphorylation (1/263), autophagy (2/263), DNA repair (5/263), cell growth and maintenance (1/263), and unknown functions (130/263) (Additional file 1).

The results of this DNA microarrays analysis were validated by real time RT-PCR using a subset of these up regulated genes. The real time RT-PCR results indicated a fold change of 1.8, 3.5, 13.0, 3.8, 3.5 and 2.3 for SLT2, ECM4, SPI1, YHR097C, ROM1 and IRA2 respectively that were consistent with the microarray results (Table 1).

Table 1 Validation of DNA microarray results. mRNAs were selected for real time RT-PCR analysis from the list of differentially expressed genes in myo1 Δ with p < 0.01(Additional File 1). Fold-change was estimated by a semi-quantitative method as described in the Methods section.

Down regulated genes in myo1Δ

A total of 284 genes were down regulated in the myo1Δ strains. As was expected for the knockout mutation, the MYO1 gene was identified among the down regulated genes in myo1Δ strains (4.5-fold, Additional file 1). Represented functional categories established by the SGD were: cell organization and biogenesis (59/284), transport (16/284), protein biosynthesis (72/284), metabolism (39/284), DNA metabolism (12/284), carbohydrate metabolism (7/284), stress response (2/284), cell cycle (5/284), signal transduction (1/284), transcription (5/284), sporulation (2/284), budding (5/284), cell growth and maintenance (1/284), protein transport (4/284), DNA damage response (1/284), RNA localization (1/284), protein amino acid phosphorylation (1/284), RNA processing (5/284), and genes of unknown function (46/284) (Additional File 1). Real time RT-PCR results for RPL30 and RPS31 indicated negative fold changes of 6.3 and 7.7 respectively that were consistent with the microarray results (Table 1).

Gene Set Enrichment Analysis

In order to identify functionally related gene sets that were differentially expressed and behaved in a coordinated fashion in the myo1Δ strains, a Gene Set Enrichment Analysis (GSEA) [12] was performed. This kind of analysis is useful for whole genome expression assays, because it allows the visualization of the data in groups and how these groups are represented in the microarray. For this analysis a gene set was created for all the 6,256 genes contained in the array, considering the t-value for each gene, and then calculating the mean of the t-values for each category. The genes were classified in categories by biological function according to the Saccharomyces genome database (SGD) [13]. The GSEA analysis was performed using the Limma package from Bioconductor [14] as described in the materials and methods section. Of the 25 categories considered for GSEA, there were five categories with a corrected p-value below the cutoff (p-value ≤ 0.0004). These categories were: protein biosynthesis, stress response, RNA processing, carbohydrate metabolism and genes of unknown function. Histograms of density versus t-value were generated for each category to observe the distribution of genes of a specific category compared with the normal distribution of all the categories on the arrays (Figure 1A.). Plots of t-values vs. A-values were created to identify the genes in each category and observe their distribution across the array (Figure 1B). Of the five gene categories selected, the protein biosynthesis and stress response categories presented the most dramatic changes in their normal distribution. Although the other three categories had p-values below the cutoff, when we observe their respective histograms there were less dramatic differences in the distribution of each group compared to the whole microarray. The protein biosynthesis category is interesting because the histogram shows a shift from the normal distribution towards negative t-values (Figure 1A). This is also observed in the corresponding t-value vs. A-value plot where we can see a greater quantity of genes biased towards negative t-values (Figure 1B). This result was predictable due to the great quantity of genes related to protein biosynthesis that were down regulated in the microarray. Within the protein biosynthesis gene category there were an unusually high number of ribosomal protein genes that were down regulated suggesting that ribosome biogenesis may be affected in the myo1Δ strains. This interpretation would be consistent with the relatively slow growth rate reported for these cells compared to their wild-type controls [5]. Conversely, in the histogram for stress response genes, there was a shift from the normal distribution of the array towards positive t-values (Figure 1A). This shift can be explained by the up-regulation of stress response genes with high t-values, therefore causing this type of shift to the right of the normal distribution curve. Considering the t-value vs. A-value plot, the stress response genes were distributed above zero towards positives t-values. Among these stress response genes were GRE1 and GRE3 that are regulated by the High Osmolarity Glycerol response-signaling pathway (HOG) [15]. Another stress response gene that had the highest t-value was HOR7. This gene encodes a putative GPI anchored cell wall protein of unknown function that was also previously reported to be related to high osmolarity response [15]. Differential expression of these genes indicates that genes in this pathway may become activated in myo1Δ strains despite the absence of hyperosmotic conditions in our experimental design. However, a hog1Δ deletion was viable in these myo1Δ strains indicating that this pathway is not essential for myo1Δ cell viability under our culture conditions (data not shown).

Figure 1
figure 1

Histograms derived from Gene Set Enrichment Analysis for categories with a corrected p-value ≤ 0.0004. A) Density versus t-value plots for protein biogenesis, stress response, unknown, carbohydrate metabolism, and RNA processing categories. Red lines represent the distribution of genes of a specific category in the array. Black lines represent the distribution of all genes in the array. B) t-value versus A-value plots for protein biogenesis, stress response, unknown, carbohydrate metabolism, and RNA processing categories. Red dots represent the genes of a specific category in the array. Black dots represent the distribution of all genes in the array. The cutoff for a significant category is based on a p-value calculated after 10,000 permutations and then a corrected p-value was calculated using the Bonferroni correction. See Methods section for details.

Although the "unknown function" category had a greater representation of differentially expressed genes in the array, there was no shift in the density vs. t-values histogram (Figure 1A) or the t-values vs. A-values plot (Figure 1B). Due to the fact that this category also has more genes represented in the yeast genome, this group was not considered of interest for this analysis. However, we are aware that this category represents the majority population of differentially expressed genes many of which exhibit extraordinary fold-changes in the double-digit range. These could represent a pool of genes with potentially novel functions related to cell proliferation, cell division, and cell wall organization.

Carbohydrate metabolism was another functional category exhibiting a large number of genes with dramatic ranges in mRNA expression levels. The GSEA results for this category showed a bimodal distribution (Figure 1A). However, in the t-value vs. A-value plot these genes were distributed throughout the array with no apparent bias towards negative or positive t-values (Figure 1B). This bimodal distribution may be significant as it may represent the existence of two sub-populations of genes with distinct behaviors.

Differential regulation of Pkc1p/Slt2p cell integrity pathway

The SLT2/MPK1 gene that encodes a serine/threonine MAP kinase involved in regulating the maintenance of cell integrity [9, 1618] was up regulated approximately 2-fold in the myo1Δ strains (Table 1 and Additional file 1). In other studies it was reported that phospho-Slt2p was also activated in the myo1Δ mutant compared to the wild type strain [19, 20]. A previous study defined a cluster of 20 genes regulated by this PKC1-mediated signaling pathway that were strongly induced by chemical perturbation of the cell wall with Congo red dye, by enzymatic treatment with Zymolyase, or by mutations in specific cell wall genes such as KRE6, GAS1, FKS1, MNN9 and KNR4 [9]. SED1, HSP12, SLT2, FBP26 and YHR097C (Additional file 1) represent five of these 20 genes that were also up regulated in the myo1Δ strains (p ≤ 0.01) further supporting that the PKC1 signaling pathway is activated. A myo1Δslt2Δ strain expressing the MYO1 gene on a pRS316 plasmid was grown at 26°C in 5-FOA to extract the plasmid and generate the myo1Δslt2Δ double knockout. After three days, no growth was observed for the myo1Δslt2Δ double knockout (Figure 2A). This result corroborates that slt2Δ is synthetically lethal in a myo1Δ mutant and supports that the cell integrity pathway is essential for viability in these strains [20].

Figure 2
figure 2

(A) Genetic disruption of the SLT2/MPK1 gene induces lethality in a myo1Δ strain. Wild type, myo1Δ, slt2Δ and myo1Δ slt2Δ pRS316-MYO1 strains were grown in CSM (1 mg/ml) 5-FOA for three days at 26°C. (B)Suppression of Nikkomycin Z hypersensitivity in myo1Δ strains overexpressing the ribosomal protein genes RPL30 and RPS31. Wt, myo1Δ, myo1Δ pRS316-RPL30 and myo1Δ pRS316-RPS31 strains were grown CSM or CSM URA- in presence or absence of 6.25 μM Nikkomycin Z in 2% glucose or 2% galactose for 48 hours at 26°C. Percent of Nikkomycin Z resistance was calculated as: (OD600 nm treated/OD600nm untreated) × (100). Dark gray histograms represent cells where expression of the plasmid is repressed with 2% glucose and light gray histograms represent cells where expression of the plasmid gene is induced with 2% galactose. (C) Determination of the activation of the cell integrity pathway in myo1Δ strains overexpressing the ribosomal protein genes RPL30 and RPS31. Wild type, myo1Δ, wildtype pRS316-RPL30, wild type pRS316-RPS31, myo1Δ pRS316-RPL30 and myo1Δ pRS316-RPS31 strains were grown in CSM or CSM URA- in the presence of 2%glucose or 2% galactose. Equal amounts of protein (75 μg) were analyzed by Western blot as described in the Methods section. The phosphorylated levels of Slt2p were observed using a mouse monoclonal antibody against phospho-p42/p44 Slt2p (p-Slt2p). The membrane was stripped and reprobed with a monoclonal antibody against Pgk1p as a loading control.

Suppression of Nikkomycin Z hypersensitivity by overexpression of ribosomal protein genes

Gene Set Enrichment Analysis showed that protein biosynthesis genes were coordinately down regulated in these myo1Δ strains (Figure 1B). RPL30 and RPS31 ribosomal protein genes were among those differentially expressed in the myo1Δ microarrays. The RPL30 gene is essential in budding yeast [21] and its protein product is a component of the large ribosomal (60S) subunit. The RPS31 gene product constitutes a fusion of a non-essential small ribosomal (40S) subunit protein with ubiquitin, which facilitates the assembly of ribosomal proteins into ribosomes and is required for the proper translation of GCN4, a transcription factor that regulates amino acid biosynthesis [22]. A previous screen for suppressors of Nikkomycin Z (NZ) hypersensitivity in myo1Δ strains identified the cDNAs encoding RPL30 and RPS31 among others [24]. Nikkomycin Z is a competitive inhibitor of Chs3p and our previous studies established the hypersensitivity of myo1Δ strains to this antibiotic which inhibits cell growth and affects cell integrity by blocking chitin synthesis [8]. Recombinant RPL30 and RPS31 conferred 67% and 69% resistance to Nikkomycin Z respectively, when overexpressed by induction with 2% galactose (Figure 2B). We confirmed increased transcription of both ribosomal protein cDNAs in a wild-type strain by real time RT-PCR (9.13 and 8.03 respectively). These results support that overexpression of ribosomal protein genes may suppresses NZ hypersensitivity in myo1Δ strains by contributing to restore cell integrity. Though the mechanism is unclear, other RP genes identified in our previous study had a similar effect.

Western blot analysis of the level of phosphorylated Slt2p in myo1Δ strains overexpressing RPL30 or RPS31

To relate suppression of NZ hypersensitivity with the restoration of cell integrity in myo1Δ strains, we assayed the level of phosphorylated Slt2p (p-Slt2p) by Western blot. It has been previously established that the cell integrity pathway is activated in myo1Δ mutants [19, 20]. Our results show increases in p-Slt2p levels for myo1Δ strains overexpressing RPL30 (myo1Δ pRS316-RPL30) and RPS31 (myo1Δ pRS316-RPS31) (Figure 2C). Moreover, we also detected increases in p-Slt2p levels for wild type strains overexpressing RPL30 (wild type pRS316-RPL30) and RPS31 (wild type pRS316-RPS31) genes in repressing and inducing conditions (Figure 2C). The former observation can be explained by leakage from the GAL1 promoter under repressing conditions (data not shown) and suggested that the activation of the cell integrity pathway can respond to very minor changes in RP gene transcription. When SLT2 mRNA levels were measured, the fold changes for strains wild type pRS316-RPL30, wild type pRS316-RPS31, myo1Δ pRS316-RPL30, and myo1Δ pRS316-RPS31 were 1.21, 2.73, 1.6 and 6.5, respectively (see Table 2). While both RP cDNAs had similar effects on p-Slt2p levels, fold changes in SLT2 mRNA levels induced by RP overexpression were consistently lower for RPL30 than for RPS31.

Table 2 Determination of the SLT2 gene mRNA fold change for wild type pRS316-RPL30, wild type pRS316-RPS31, myo1Δ pRS316-RPL30, and myo1Δ pRS316-RPS31 strains grown under inducing (2% galactose) and repressing (2% glucose) conditions.

Conclusion

Following our analysis of the global mRNA expression profile in myo1Δ strains of Saccharomyces cerevisiae we found that there were 547 differentially regulated genes. Gene Set Enrichment Analysis indicated that stress response and protein biosynthesis gene categories are inversely related in this mutant. These results support that cell wall morphogenetic processes were affected by myosin II deficiency. The up-regulation of stress response genes in the myo1Δ strains reflects the essentiality of the cell integrity pathway and the importance of the PKC1-mediated signaling pathway for survival in these strains. The coordinated down-regulation of protein biosynthesis genes can be related to a loss of cell wall integrity in this mutant because overexpression of RPL30 and RPS31, two differentially regulated ribosomal protein (RP) genes, conferred Nikkomycin Z resistance to the hypersensitive myo1Δ strains. Overexpression of RPL30 and RPS31 did not shut down the cell integrity pathway supporting that their overexpression promotes activation of the PKC1 pathway.

The high frequency of ribosomal protein genes in the down regulated gene category implies that ribosome biogenesis is greatly affected by myosin type II deficiency. It has been proposed that yeast cells regulate ribosome biogenesis primarily by responding to growth limiting environmental signals [23]. Cell wall and plasma membrane stress may therefore represent such signals that can activate signaling pathways for regulating growth rate in the absence of myosin II. We propose that this alteration in the normal regulation of these processes may be necessary for the completion of cytokinesis in the absence of a contractile ring that must involve increased chitin synthesis and extensive cell wall remodeling at the bud neck.

We previously demonstrated that overexpression of the cDNA for ubiquitin conjugating enzyme UBC4 lowered the p-Slt2p levels and conferred resistance to NZ in myo1Δ strains, suggesting that cell integrity was restored [24]. We have now discovered that overexpression of RPL30 and RPS31 can also up regulate the cell integrity pathway to achieve a similar effect. It is unclear how such an alteration in the level of transcription of individual RP genes can affect this pathway but it is known that ribosome biosynthesis is a tightly regulated process in yeast that can be affected by altering the stoichiometry of the individual ribosomal subunits [25]. This observation could therefore represent an effect of up regulation of ribosome biosynthesis on cell integrity. Alternatively, the response of the PKC1 pathway to changes in RP gene transcription may represent one aspect of a novel feedback regulatory loop. Others have reported that inhibition of the TOR pathway by rapamycin treatment activated the PKC1 cell integrity pathway in yeast cells [26]. It was acknowledged that inhibition of the TOR pathway also results in repression of the ribosomal biogenesis process [26]. The similarity of the described phenotypes with those of myo1Δ strains merits further investigation to determine if down regulation of the TOR pathway is a contributing factor to the phenotypes in these strains.

Methods

Strains and culture conditions

All the experiments were performed using Saccharomyces cerevisiae wild type and myo1Δ strains listed in Table 3. YJR6 is a myo1Δ ::HIS5 strain generated by homologous recombination in the parental haploid wild type strain, MGD-353-46D using a PCR based method. YJR12 (wild type) and YJR13 (myo1Δ) strains were obtained as haploid segregants from a cross between YJR6 and BY4741 (obtained from ATCC). Cultures were grown overnight at 26°C to an optical density between 0.5–0.8 (OD600) in complete synthetic media (CSM, 2% glucose, 1× Nitrogen base) or Histidine dropout media (CSM-HIS-) with continuous shaking at 200 rpm. Where indicated in Table 3, myo1Δ strains were transformed with plasmids containing cDNAs for RPL30 (pRS316-RPL30) or RPS31 (pRS316-RPS31) regulated by the GAL1 promoter, or MYO1 (pRS316-MYO1) regulated by its natural promoter.

Table 3 Strains used in this study.

RNA extraction procedure

Total RNA was extracted from 4 × 107 cells derived from triplicate biological replicate cultures of strains MGD353-46D, YJR6, YJR12, and YJR13 using the RNeasy Mini Kit for isolation of total RNA (Qiagen, Hilden, Germany) following manufacturer's instructions. RNA concentrations were determined by measuring absorbance at 260 nm using a Nanodrop spectrophotometer (Nanodrop Technologies, Wilmington, DE). The purity and integrity of the RNA was monitored using an Agilent Bioanalyzer (Agilent Technologies, Palo Alto, CA) following manufacturer's instructions.

DNA microarrays hybridization and analysis

1.0 μg of total RNA from each sample was amplified using the Low RNA Input Fluorescent Linear Amplification kit (Agilent Technologies, Palo Alto, CA). The amplified cRNA was labeled with 10 mM Cyanine 5-CTP (Cy5) or Cyanine 3-CTP (Cy3) (Perkin Elmer Life Sciences, Boston, MA). Labeled cRNA's were purified with Qiagen RNeasy mini spin columns and dye incorporation was monitored on an Agilent Bioanalyzer. Hybridization of Cy5 and Cy3 labeled cRNA's were performed using Yeast Oligo Microarray slides and hybridization kit from Agilent Technologies in 1012AG hybridization chambers (Sheldon Manufacturing, Cornelius, Oregon) at 60°C overnight. Slides were washed at high stringency and scanned with a VersArray Chip Reader system (BioRad, Hercules, CA) at a resolution of 5 μm with detector sensitivity values between 704–800 and laser power at 85%. Scanned images were transferred to the Imagene 3.0 software program (Biodiscovery, El Segundo, CA) for further analysis to locate spots, adjust the appropriate grid, and obtain the Cy3 and Cy5 TIFF files. The microarrays raw data generated with Imagene 3.0 were analyzed using Limma software (Bioconductor Package 1.7). The data was prepared for analysis by correcting for background intensity. The individual data sets were normalized using the locally weighted linear regression (Lowess) within each array. After normalization, the difference between the experimental and control signal was calculated, replicates were combined, and their averages were calculated. The fold change in gene expression was calculated by 2(M), where M is the log2-fold change after background correction and normalization. An Empirical Bayes Statistics for differential expression analysis [27] (eBayes statistics) was performed by Limma. Genes with a p-value ≤ 0.01 were established as a cutoff for differential expression. In addition, a false discovery rate (FDR) test [28] was performed by Limma [14]. The microarray raw data and the processed data are available at Gene Expression Omnibus of NCBI (GSE5931) [29].

Gene Set Enrichment Analysis

Gene Set Enrichment Analysis was performed using the Limma package of Bioconductor [14]. For this purpose a gene set file was created classifying all the genes included in the microarray into groups according to their specific involvement in a biological process. From this gene set, a total of 25 categories were represented according to the Osprey network visualization software [30], matching the category with the t-value. Then, the Limma software calculated the average of the t-value for each category or biological process. Significance of differential expression as determined by the enrichment analysis was recalculated 10,000 times. A corrected p-value was obtained from the analysis using the Bonferroni correction. The initial cutoff to determine that the gene set was differentially expressed was p-value ≤ 0.01. Based on this correction, the cutoff for significance was established at a p-value ≤ 0.0004.

Confirmation of microarray data by real time RT-PCR

Real time RT-PCR assays were performed with 30 ng of total RNA using the Quantitec SYBR Green RT-PCR kit (Qiagen, Valencia, CA) with primers at 0.5 μM (Table 4) and a final reaction volume of 25.0 μL, following manufacturer recommendations. The reactions were performed following the conditions recommended by the manufacturer using the iCycler iQ Multicolor real time PCR detection system (BioRad, Hercules, CA). The PCR quantification and melting curves were generated using the iCycler new version software (BioRad, Hercules, CA). The sequences of the forward and reverse primers for the selected mRNAs are listed in Table 4. The fold change was determined by the 2ΔΔCt method [31]. The ΔCt of the control and experimental samples was calculated from the threshold cycle of the target gene minus the threshold cycle of the reference gene (ACT1). The ΔΔCt was calculated by subtracting the ΔCt of the control sample minus ΔCt of the experimental sample.

Table 4 Primers used for real time RT-PCR in this study.

Nikkomycin Z sensitivity assay in myo1Δ strains overexpressing ribosomal protein genes

The screening for Nikkomycin Z sensitivity suppression in myo1Δ was performed as described in [24]. Briefly, 1 × 105 cells/ml of WT, myo1Δ, myo1Δ p RPL30 and myo1Δ p RPS31 cells were grown in ura-glu- or ura-gal- in the presence or absence of 6.25 μM Nikkomycin Z at 26°C and 225 rpm. The OD600 nm was determined after 48 hrs of growth. The % Nikkomycin Z resistance represents: (OD600 nm treated/OD600 nm untreated) × 100.

SLT2 gene disruption in myo1Δ a background

The SLT2 gene disruption in the YJR13 pRS316-MYO1 strain was performed by replacing the SLT2 gene with a KanMX4 module by homologous recombination using a PCR based method [32]. The myo1Δ slt2Δ pRS316-MYO1 strain was grown in CSM 1 mg/ml 5-Fluoroorotic acid (CSM 5-FOA) to extract the plasmid and obtain a myo1Δslt2Δ double mutant.

Western blot analysis of phosphorylated Slt2p levels

Wild type, myo1Δ, wild type pRS316-RPL30, wild type pRS316-RPS31, myo1Δ pRS316-RPL30, and myo1Δ pRS316-RPS31 yeast strains were grown in selective medium in 2% glucose or 2% galactose to an OD600 0.5–0.8 at 26°C, overnight. Cells were centrifuged for 5 minutes, washed in ice cold medium, and resuspended in lysis buffer as described in [24]. Protein extract were centrifuged at 13,000 for 10 minutes at 4°C, the supernatant was removed and quantified using the DC Protein Assay method (BioRad). Total protein extracts (75 μg) were separated in a 10% SDS-PAGE gel and transferred to a nitrocellulose membrane at 30 V for 16 h at 4°C. Membrane was incubated with anti-phospho-p42/44 MAP kinase monoclonal antibody (1:2000) (Cell Signaling Technology) to detect phosphorylated Slt2p levels. Membrane was stripped and reprobed with a mouse monoclonal antibody against Pgk1p (1:125 dilution) (Molecular Probes, Invitrogen) as loading control.

References

  1. VerPlank L, Li R: Cell cycle-regulated trafficking of Chs2 controls actomyosin ring stability during cytokinesis. Mol Biol Cell. 2005, 16: 2529-2543. 10.1091/mbc.E04-12-1090.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  2. Lippincott J, Li R: Sequential assembly of myosin II, an IQGAP-like protein, and filamentous actin to a ring structure involved in budding yeast cytokinesis. J Cell Biol. 1998, 140: 355-366. 10.1083/jcb.140.2.355.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  3. Bi E, Maddox P, Lew DJ, Salmon ED, McMillan JN, Yeh E, Pringle JR: Involvement of an actomyosin contractile ring in Saccharomyces cerevisiae cytokinesis. J Cell Biol. 1998, 142: 1301-1312. 10.1083/jcb.142.5.1301.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Vallen EA, Caviston J, Bi E: Roles of Hof1p, Bni1p, Bnr1p, and myo1p in cytokinesis in Saccharomyces cerevisiae. Mol Biol Cell. 2000, 11: 593-611.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Rodriguez JR, Paterson BM: Yeast myosin heavy chain mutant: Maintenance of the cell type specific budding pattern and the normal deposition of chitin and cell wall components requieres an intact myosin heavy chain gene. Cell MotilCytoskeleton. 1990, 17: 301-308. 10.1002/cm.970170405.

    Article  CAS  Google Scholar 

  6. Cruz JA, Garcia R, Rodriguez-Orengo JF, Rodriguez-Medina JR: Increased chitin synthesis in response to type II myosin deficiency in Saccharomyces cerevisiae. Mol Cell Biol Res Commun. 2000, 3: 20-25. 10.1006/mcbr.2000.0180.

    Article  PubMed  CAS  Google Scholar 

  7. Rodriguez-Medina JR, Cruz JA, Robbins PW, Bi E, Pringle JR: Elevated expression of chitinase 1 and chitin synthesis in myosin II-deficient Saccharomyces cerevisiae. Cell Mol Biol (Noisy-le-grand). 1998, 44: 919-925.

    CAS  Google Scholar 

  8. Rivera-Molina FE, González-Crespo S, Maldonado-De la Cruz Y, Ortiz-Betancourt JM, Rodríguez-Medina JR: 2,3-Butanedione monoxime increases sensitivity to Nikkomycin Z in the budding yeast Saccharomyces cerevisiae. World J Microbio Biotechnol. 2006, 22: 255-260. 10.1007/s11274-005-9028-x.

    Article  Google Scholar 

  9. Garcia R, Bermejo C, Grau C, Perez R, Rodriguez-Pena JM, Francois J, Nombela C, Arroyo J: The global transcriptional response to transient cell wall damage in Saccharomyces cerevisiae and its regulation by the cell integrity signaling pathway. J Biol Chem. 2004, 279: 15183-15195. 10.1074/jbc.M312954200.

    Article  PubMed  CAS  Google Scholar 

  10. Hughes TR, Marton MJ, Jones AR, Roberts CJ, Stoughton R, Armour CD, Bennett HA, Coffey E, Dai H, He YD, Kidd MJ, King AM, Meyer MR, Slade D, Lum PY, Stepaniants SB, Shoemaker DD, Gachotte D, Chakraburtty K, Simon J, Bard M, Friend SH: Functional discovery via a compendium of expression profiles. Cell. 2000, 102: 109-126. 10.1016/S0092-8674(00)00015-5.

    Article  PubMed  CAS  Google Scholar 

  11. Cho RJ, Campbell MJ, Winzeler EA, Steinmetz L, Conway A, Wodicka L, Wolfsberg TG, Gabrielian AE, Landsman D, Lockhart DJ, Davis RW: A genome-wide transcriptional analysis of the mitotic cell cycle. Mol Cell. 1998, 2: 65-73. 10.1016/S1097-2765(00)80114-8.

    Article  PubMed  CAS  Google Scholar 

  12. Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, Houstis N, Daly MJ, Patterson N, Mesirov JP, Golub TR, Tamayo P, Spiegelman B, Lander ES, Hirschhorn JN, Altshuler D, Groop LC: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34: 267-273. 10.1038/ng1180.

    Article  PubMed  CAS  Google Scholar 

  13. [http://www.yeastgenome.org]

  14. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Kapteyn JC, ter Riet B, Vink E, Blad S, De Nobel H, Van Den Ende H, Klis FM: Low external pH induces HOG1-dependent changes in the organization of the Saccharomyces cerevisiae cell wall. Mol Microbiol. 2001, 39: 469-479. 10.1046/j.1365-2958.2001.02242.x.

    Article  PubMed  CAS  Google Scholar 

  16. Hohmann S: Osmotic stress signaling and osmoadaptation in yeasts. Microbiol Mol Biol Rev. 2002, 66: 300-372. 10.1128/MMBR.66.2.300-372.2002.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. de Nobel H, Ruiz C, Martin H, Morris W, Brul S, Molina M, Klis FM: Cell wall perturbation in yeast results in dual phosphorylation of the Slt2/Mpk1 MAP kinase and in an Slt2-mediated increase in FKS2-lacZ expression, glucanase resistance and thermotolerance. Microbiology. 2000, 146 ( Pt 9): 2121-2132.

    Article  CAS  Google Scholar 

  18. Roberts CJ, Nelson B, Marton MJ, Stoughton R, Meyer MR, Bennett HA, He YD, Dai H, Walker WL, Hughes TR: Signaling and circuitry of multiple MAPK pathways revealed by a matrix of global gene expression profiles. Science. 2000, 287: 873-880. 10.1126/science.287.5454.873.

    Article  PubMed  CAS  Google Scholar 

  19. Rivera-Molina F: Biochemical and Genetic Analysis of Myosin II Function in the Trafficking of the Chitin Synthase III Catalytic Subunit, Ph.D. thesis. Department of Biochemistry. 2005, , University of Puerto Rico, Medical Sciences Campus, Department of Biochemistry

    Google Scholar 

  20. Nishihama R: Septin Dependent, Actomyosin Ring Independent Cytokinesis in Yeast [Abstract]. Amer Soc Cell Biol. 2003, 43:L81.:

    Google Scholar 

  21. Giaever G, Chu AM, Ni L, Connelly C, Riles L, Veronneau S, Dow S, Lucau-Danila A, Anderson K, Andre B, Arkin AP, Astromoff A, El-Bakkoury M, Bangham R, Benito R, Brachat S, Campanaro S, Curtiss M, Davis K, Deutschbauer A, Entian KD, Flaherty P, Foury F, Garfinkel DJ, Gerstein M, Gotte D, Guldener U, Hegemann JH, Hempel S, Herman Z, Jaramillo DF, Kelly DE, Kelly SL, Kotter P, LaBonte D, Lamb DC, Lan N, Liang H, Liao H, Liu L, Luo C, Lussier M, Mao R, Menard P, Ooi SL, Revuelta JL, Roberts CJ, Rose M, Ross-Macdonald P, Scherens B, Schimmack G, Shafer B, Shoemaker DD, Sookhai-Mahadeo S, Storms RK, Strathern JN, Valle G, Voet M, Volckaert G, Wang CY, Ward TR, Wilhelmy J, Winzeler EA, Yang Y, Yen G, Youngman E, Yu K, Bussey H, Boeke JD, Snyder M, Philippsen P, Davis RW, Johnston M: Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002, 418: 387-391. 10.1038/nature00935.

    Article  PubMed  CAS  Google Scholar 

  22. Mueller PP, Grueter P, Hinnebusch AG, Trachsel H: A ribosomal protein is required for translational regulation of GCN4 mRNA. Evidence for involvement of the ribosome in eIF2 recycling. J Biol Chem. 1998, 273: 32870-32877. 10.1074/jbc.273.49.32870.

    Article  PubMed  CAS  Google Scholar 

  23. Levy S, Ihmels J, Carmi M, Weinberger A, Friedlander G, Barkai N: Strategy of transcription regulation in the budding yeast. PLoS ONE. 2007, 2: e250-10.1371/journal.pone.0000250.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Diaz-Blanco NL, Rodriguez-Medina JR: Dosage rescue by UBC4 restores cell wall integrity in Saccharomyces cerevisiae lacking the myosin type II gene MYO1. Yeast. 2007, 24: 343-355. 10.1002/yea.1481.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  25. Zhao Y, Sohn JH, Warner JR: Autoregulation in the biosynthesis of ribosomes. Mol Cell Biol. 2003, 23: 699-707. 10.1128/MCB.23.2.699-707.2003.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Torres J, Di Como CJ, Herrero E, De La Torre-Ruiz MA: Regulation of the cell integrity pathway by rapamycin-sensitive TOR function in budding yeast. J Biol Chem. 2002, 277: 43495-43504. 10.1074/jbc.M205408200.

    Article  PubMed  CAS  Google Scholar 

  27. Smyth G: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: 1-25.

    Google Scholar 

  28. Benjamini Y, Hochberg Y: Controlling the False Discovery Rate:a practical and powerful approach to multiple testing. JR Statist Soc B. 1995, 57: 289-300.

    Google Scholar 

  29. [http://www.ncbi.nlm.nih.gov/geo/]

  30. Breitkreutz BJ, Stark C, M. T: Osprey: "a network visualization system". Genome Biology. 2003, 4: R22-10.1186/gb-2003-4-3-r22.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

    Article  PubMed  CAS  Google Scholar 

  32. Wach A, Brachat A, Pohlmann R, Philippsen P: New heterologous modules for classical or PCR-based gene disruptions in Saccharomyces cerevisiae. Yeast. 1994, 10: 1793-1808. 10.1002/yea.320101310.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We acknowledge outstanding technical support by Lilliam Villanueva-Alicea and Sahily Gonzalez-Crespo. Support for this work was provided by United States Public Health Service grants to the University of Puerto Rico from NIGMS-SCORE (S06-GM008224) and NCRR-RCMI (G12-RR03051). JFRQ and FRM were supported by NIGMS-RISE (R25-GM61838). FRM received partial support from NIH-MIRT (T37-TW0096-01A1). NDB was supported by NIGMS-MARC (F34-GM69277).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to José R Rodríguez-Medina.

Additional information

Authors' contributions

JFRQ participated in the data analysis, real time RT-PCR experiments, genetic knockout experiments, Western blot experiments, and writing of the manuscript. DGG and FERM participated in the data analysis. NDB performed Nikkomycin Z dosage suppressor screening, identification and characterization. RAI participated in the design of experiments, supervised data analysis, and provided expertise in DNA microarrays data processing and statistical analysis software. JRRM participated in the design of experiments, data analysis, interpretation of results, and writing of the manuscript.

Electronic supplementary material

12864_2007_1228_MOESM1_ESM.pdf

Additional File 1: A text table of differentially expressed genes (p-value ≤ 0.01) in myo1Δ strains divided according to biological process. All differentially expressed genes are included except those with unknown biological process. The list of genes with unknown biological process is available at Gene Expression Omnibus on the NCBI website (GSE5931). (PDF 136 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Rodríguez-Quiñones, J.F., Irizarry, R.A., Díaz-Blanco, N.L. et al. Global mRNA expression analysis in myosin II deficient strains of Saccharomyces cerevisiae reveals an impairment of cell integrity functions. BMC Genomics 9, 34 (2008). https://doi.org/10.1186/1471-2164-9-34

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-9-34

Keywords