Budding Yeast Saccharomyces Cerevisiae Research Articles

A high molecular mass zinc transporter MTP12 forms a functional heteromeric complex with MTP5 in the Golgi in Arabidopsis thaliana.

Zinc (Zn) is an essential micronutrient required for plant growth and development. In Arabidopsis thaliana, several families of Zn transporters engaged in Zn import, export and intracellular compartmentalization play important roles in Zn homeostasis. We describe a novel Zn transporter, A. thaliana metal tolerance protein 12 (AtMTP12), which belongs to the cation diffusion facilitator family. AtMTP12 is predicted to consist of 798 amino acids and have 14 transmembrane segments. The expression of AtMTP12 in suspension-cultured cells was not affected by Zn deficiency or excess. Heterologous expression in a mutant of budding yeast (Saccharomyces cerevisiae) that lacks Msc2p, an orthologue of AtMTP12, revealed that AtMTP12 complements the growth phenotype of the msc2 mutant when AtMTP5t1, one of the splicing variants of AtMTP5, is coexpressed. Transient expression of AtMTP12-fused green fluorescent protein in A. thaliana mesophyll protoplasts demonstrated that AtMTP12 is localized to the Golgi apparatus. Moreover, AtMTP12 and AtMTP5t1 interact in the Golgi, as determined by a bimolecular fluorescence complementation assay. These results suggest that AtMTP12 forms a functional complex with AtMTP5t1 to transport Zn into the Golgi. Nucleotide sequence data for full-length of AtMTP12 is available in the DDBJ/EMBL/GenBank database under accession number AB986563.

Nuclear localization domains of GATA activator Gln3 are required for transcription of target genes through dephosphorylation in Saccharomyces cerevisiae

The GATA transcription activator Gln3 in the budding yeast (Saccharomyces cerevisiae) activates transcription of nitrogen catabolite repression (NCR)-sensitive genes. In cells grown in the presence of preferred nitrogen sources, Gln3 is phosphorylated in a TOR-dependent manner and localizes in the cytoplasm. In cells grown in non-preferred nitrogen medium or treated with rapamycin, Gln3 is dephosphorylated and is transported from the cytoplasm to the nucleus, thereby activating the transcription of NCR-sensitive genes. Caffeine treatment also induces dephosphorylation of Gln3 and its translocation to the nucleus and transcription of NCR-sensitive genes. However, the details of the mechanism by which phosphorylation controls Gln3 localization and transcriptional activity are unknown. Here, we focused on two regions of Gln3 with nuclear localization signal properties (NLS-K, and NLS-C) and one with nuclear export signal (NES). We constructed various mutants for our analyses: gln3 containing point mutations in all potential phosphoacceptor sites (Thr-339, Ser-344, Ser-347, Ser-355, Ser-391) in the NLS and NES regions to produce non-phosphorylatable (alanine) or mimic-phosphorylatable (aspartic acid) residues; and deletion mutants. We found that phosphorylation of Gln3 was impaired in all of these mutations and that the aspartic acid substitution mutants showed drastic reduction of Gln3-mediated transcriptional activity despite the fact that the mutations had no effect on nuclear localization of Gln3. Our observations suggest that these regions are required for transcription of target genes presumably through dephosphorylation.

Asymmetric nucleosomes flank promoters in the budding yeast genome.

Nucleosomes in active chromatin are dynamic, but whether they have distinct structural conformations is unknown. To identify nucleosomes with alternative structures genome-wide, we used H4S47C-anchored cleavage mapping, which revealed that 5% of budding yeast (Saccharomyces cerevisiae) nucleosome positions have asymmetric histone-DNA interactions. These asymmetric interactions are enriched at nucleosome positions that flank promoters. Micrococcal nuclease (MNase) sequence-based profiles of asymmetric nucleosome positions revealed a corresponding asymmetry in MNase protection near the dyad axis, suggesting that the loss of DNA contacts around H4S47 is accompanied by protection of the DNA from MNase. Chromatin immunoprecipitation mapping of selected nucleosome remodelers indicated that asymmetric nucleosomes are bound by the RSC chromatin remodeling complex, which is required for maintaining nucleosomes at asymmetric positions. These results imply that the asymmetric nucleosome-RSC complex is a metastable intermediate representing partial unwrapping and protection of nucleosomal DNA on one side of the dyad axis during chromatin remodeling.

Important role of catalase in the cellular response of the budding yeast Saccharomyces cerevisiae exposed to ionizing radiation.

Ionizing radiation indirectly causes oxidative stress in cells via reactive oxygen species (ROS), such as hydroxyl radicals (OH(-)) generated by the radiolysis of water. We investigated how the catalase function was affected by ionizing radiation and analyzed the phenotype of mutants with a disrupted catalase gene in Saccharomyces cerevisiae exposed to radiation. The wild-type yeast strain and isogenic mutants with disrupted catalase genes were exposed to various doses of (60)Co gamma-rays. There was no difference between the wild-type strain and the cta1 disruption mutant following exposure to gamma-ray irradiation. In contrast, there was a significant decrease in the ctt1 disruption mutant, suggesting that this strain exhibited decreased survival on gamma-ray exposure compared with other strains. In all three strains, stationary phase cells were more tolerant to the exposure of gamma-rays than exponential phase cells, whereas the catalase activity in the wild-type strain and cta1 disruption mutant was higher in the stationary phase than in the exponential phase. These data suggest a correlation between catalase activity and survival following gamma-ray exposure. However, this correlation was not clear in the ctt1 disruption mutant, suggesting that other factors are involved in the tolerance to ROS induced by irradiation.

Applied usage of yeast spores as chitosan beads.

In this study, we present a nonhazardous biological method of producing chitosan beads using the budding yeast Saccharomyces cerevisiae. Yeast cells cultured under conditions of nutritional starvation cease vegetative growth and instead form spores. The spore wall has a multilaminar structure with the chitosan layer as the second outermost layer. Thus, removal of the outermost dityrosine layer by disruption of the DIT1 gene, which is required for dityrosine synthesis, leads to exposure of the chitosan layer at the spore surface. In this way, spores can be made to resemble chitosan beads. Chitosan has adsorptive features and can be used to remove heavy metals and negatively charged molecules from solution. Consistent with this practical application, we find that spores are capable of adsorbing heavy metals such as Cu(2+), Cr(3+), and Cd(2+), and removal of the dityrosine layer further improves the adsorption. Removal of the chitosan layer decreases the adsorption, indicating that chitosan works as an adsorbent in the spores. Besides heavy metals, spores can also adsorb a negatively charged cholesterol derivative, taurocholic acid. Furthermore, chitosan is amenable to chemical modifications, and, consistent with this property, dit1Δ spores can serve as a carrier for immobilization of enzymes. Given that yeast spores are a natural product, our results demonstrate that they, and especially dit1Δ mutants, can be used as chitosan beads and used for multiple purposes.

Organization and evolution of a novel cervid satellite DNA with yeast CDEI-like repeats

AbstractBackgroundIt has been proposed that pericentromeric satellite DNA arises from the progressive proximal expansion of ancient centromeric DNA. In an attempt to recover putative ancestral centromeric DNA, we microdissected the pericentromeric/centromeric DNA from the chromosome X + 3 of Indian muntjac (Muntiacus muntjak vaginalis) and constructed a microclone-library of the X + 3 centromeric DNA.ResultsA new cervid satellite DNA element, designated as satellite VI, was isolated from this library. Fluorescence in situ hybridization (FISH) studies revealed that satellite VI is predominately located on the distal pericentromeric region of the Indian muntjac chromosome X + 3 and on the pericentromeres of several Old World deer species studied. Its sequence is organized as 11-bp monomeric (ATCACGTGGGA) tandem repeats. Further sequencing on a BAC clone of Indian muntjac harboring this repeat showed that an array of this repeat stretches over approximately 5 kb followed by approximately 3 kb of interspersed repetitive sequences, such as long interspersed elements (LINEs), short interspersed elements (SINEs), and long terminal repeats (LTRs).ConclusionsBased on the chromosomal localization, genomic and sequence organization, and copy numbers of satellite VI in deer species studied, we postulate that this newly found satellite DNA could be a putative ancient cervidic centromeric DNA that is still preserved in some Old World deer. Interestingly, the first eight nucleotides of the 11-bp monomeric consensus sequences are highly conserved and identical to the CDEI element in the centromere of the budding yeastSaccharomyces cerevisiae. The centromeric/pericentromeric satellite DNA harboring abundant copies of CDEI sequences is the first found in a mammalian species. Several zipper-like d (GGGA)2motifs were also found in the (ATCACGTGGGA)nrepeat of satellite VI DNA. Whether the satellite VI is structurally and functionally correlated with the CDEI of centromere of budding yeast and whether a zipper-like structure forms in satellite VI require further studies.

Precise estimates of mutation rate and spectrum in yeast.

Mutation is the ultimate source of genetic variation. The most direct and unbiased method of studying spontaneous mutations is via mutation accumulation (MA) lines. Until recently, MA experiments were limited by the cost of sequencing and thus provided us with small numbers of mutational events and therefore imprecise estimates of rates and patterns of mutation. We used whole-genome sequencing to identify nearly 1,000 spontaneous mutation events accumulated over ∼311,000 generations in 145 diploid MA lines of the budding yeast Saccharomyces cerevisiae. MA experiments are usually assumed to have negligible levels of selection, but even mild selection will remove strongly deleterious events. We take advantage of such patterns of selection and show that mutation classes such as indels and aneuploidies (especially monosomies) are proportionately much more likely to contribute mutations of large effect. We also provide conservative estimates of indel, aneuploidy, environment-dependent dominant lethal, and recessive lethal mutation rates. To our knowledge, for the first time in yeast MA data, we identified a sufficiently large number of single-nucleotide mutations to measure context-dependent mutation rates and were able to (i) confirm strong AT bias of mutation in yeast driven by high rate of mutations from C/G to T/A and (ii) detect a higher rate of mutation at C/G nucleotides in two specific contexts consistent with cytosine methylation in S. cerevisiae.

Phosphorylation-Regulated Transitions in an Oligomeric State Control the Activity of the Sae2 DNA Repair Enzyme

In the DNA damage response, many repair and signaling molecules mobilize rapidly at the sites of DNA double-strand breaks. This network of immediate responses is regulated at the level of posttranslational modifications that control the activation of DNA processing enzymes, protein kinases, and scaffold proteins to coordinate DNA repair and checkpoint signaling. Here we investigated the DNA damage-induced oligomeric transitions of the Sae2 protein, an important enzyme in the initiation of DNA double-strand break repair. Sae2 is a target of multiple phosphorylation events, which we identified and characterized in vivo in the budding yeast Saccharomyces cerevisiae. Both cell cycle-dependent and DNA damage-dependent phosphorylation sites in Sae2 are important for the survival of DNA damage, and the cell cycle-regulated modifications are required to prime the damage-dependent events. We found that Sae2 exists in the form of inactive oligomers that are transiently released into smaller active units by this series of phosphorylations. DNA damage also triggers removal of Sae2 through autophagy and proteasomal degradation, ensuring that active Sae2 is present only transiently in cells. Overall, this analysis provides evidence for a novel type of protein regulation where the activity of an enzyme is controlled dynamically by posttranslational modifications that regulate its solubility and oligomeric state.

Electrostatics-Dependent Shape of the Intrinsically-Disordered Protein Sic1

For intrinsically disordered proteins (IDPs) containing a high density of charged residues, a Gaussian chain is an inadequate description of the protein's shape. One such IDP is Sic1, a highly positively charged IDP in the budding yeast Saccharomyces Cerevisiae which prevents the cell cycle from entering the S-phase from the G1-phase. Sic1's binding affinity for Cdc4 is highly phosphorylation state dependent, although the corresponding physical basis is not fully understood. NMR data supports the presence of a dynamic complex of Sic1:Cdc4 and a poly-electrostatic model has been proposed by Forman-Kay and coworkers.We studied the Sic1 N-terminal targeting region (1-90) to better understand the role of intrachain electrostatics, and to compare with the structural ensembles calculated from NMR and SAXS data. Sic1 is characterized using time-resolved fluorescence anisotropy (TRFA), which is sensitive to rotational and conformational dynamics. Sic1 exists in a dynamic ensemble of conformations and ensemble-averaged experiments have limitations in identifying and characterizing static and/or dynamic inhomogeneity in the motional dynamics. Therefore, we also performed burst spectroscopy and single-molecule TRFA on freely diffusing Sic1. Shape and flexibility were probed by varying the degree of intrachain repulsion screening through adjusting salt concentrations and described within a polymer physics framework, the polyelectrolyte model. To investigate the relative contributions of “global rotation” and “conformational flexibility”, Sic1 was labelled at three different sites. Sic1 is found to be well modelled as a prolate ellipsoid with internal flexibility.The single-molecule derived TRFA distribution data may be valuable as a constraint incorporated in future conformational ensemble calculations, complementary to SAXS and NMR data. Additionally, these measurements raise questions about the accuracy of highly charged IDP's radii of gyration when calculated from sm-FRET derived end-to-end distances, which often assume a Gaussian chain model for the end-to-end distance probability distribution.

Cdk1 Regulates the Temporal Recruitment of Telomerase and Cdc13-Stn1-Ten1 Complex for Telomere Replication

In budding yeast (Saccharomyces cerevisiae), the cell cycle-dependent telomere elongation by telomerase is controlled by the cyclin-dependent kinase 1 (Cdk1). The telomere length homeostasis is balanced between telomerase-unextendable and telomerase-extendable states that both require Cdc13. The recruitment of telomerase complex by Cdc13 promotes telomere elongation, while the formation of Cdc13-Stn1-Ten1 (CST) complex at the telomere blocks telomere elongation by telomerase. However, the cellular signaling that regulates the timing of the telomerase-extendable and telomerase-unextendable states is largely unknown. Phosphorylation of Cdc13 by Cdk1 promotes the interaction between Cdc13 and Est1 and hence telomere elongation. Here, we show that Cdk1 also phosphorylates Stn1 at threonine 223 and serine 250 both in vitro and in vivo, and these phosphorylation events are essential for the stability of the CST complexes at the telomeres. By controlling the timing of Cdc13 and Stn1 phosphorylations during cell cycle progression, Cdk1 regulates the temporal recruitment of telomerase complexes and CST complexes to the telomeres to facilitate telomere maintenance.

Structure based Functional Distinction between Cln1 and Cln2 Depends on the Ubiquitin-Proteasome Pathway

Cln1 and Cln2, G1/S cyclins of the ascomycetous budding yeast Saccharomyces cerevisiae (S. cerevisiae), oscillate during the cell cycle, rising in late G1 and falling in early S phase. We have been tried to elucidate the structure basis of the functional distinction between Cln1 and Cln2. Here we performed in silico simulations: construction and evaluation of three dimensional structures of Cln1-Cdc28 and Cln2-Cdc28 complexes. Our in silico simulations suggested that the interaction of Cln1 and Cln2 with Cdc28 were in the two distinct situations, designated as flip and flop conformation, at the extra amino acid region in the cyclin box of Cln1 and Cln2. We speculated the trigger of this flip-flop conversion of the extra amino acid region in the cyclin box of Cln1-Cdc28 and Cln2-Cdc28 might be regulated by the ubiquitination of the sequences rich in Pro (P), Glu (E), Ser (S) and Thr (T), so-called PEST motifs, in Cln1 and Cln2. Furthermore, we presumed that the functional superiority between Cln1 and Cln2 in the G1/S phase of S. cerevisiae might be controlled by flip-flop conversion and ubiquitin-proteasome pathway.

Cap Completion and C-Terminal Repeat Domain Kinase Recruitment Underlie the Initiation-Elongation Transition of RNA Polymerase II

After transcription initiation, RNA polymerase (Pol) II escapes from the promoter and recruits elongation factors. The molecular basis for the initiation-elongation factor exchange during this transition remains poorly understood. Here, we used chromatin immunoprecipitation (ChIP) to elucidate the initiation-elongation transition of Pol II in the budding yeast Saccharomyces cerevisiae. We show that the early Pol II elongation factor Spt5 contributes to stable recruitment of the mRNA capping enzymes Cet1, Ceg1, and Abd1. Genome-wide occupancy for Cet1 and Ceg1 is restricted to the transcription start site (TSS), whereas occupancy for Abd1 peaks at ~110 nucleotides downstream, and occupancy for the cap-binding complex (CBC) rises subsequently. Abd1 and CBC are important for recruitment of the kinases Ctk1 and Bur1, which promote elongation and capping enzyme release. These results suggest that cap completion stimulates productive Pol II elongation.

Antimicrobial activities of ant Ponericin W1 against plant pathogens in vitro and the disease resistance in its transgenic Arabidopsis

The antimicrobial peptides (AMPs) exhibit a broad antimicrobial spectrum. The application of AMPs from non-plant organisms attracts considerable attention in plant disease resistance engineering. Ponericin W1, isolated from the venom of ant (Pachycondyla goeldii), shows antimicrobial activities against Gram-positive bacteria, Gram-negative bacteria and the budding yeast (Saccharomyces cerevisiae); however, it is not clear whether Ponericin W1 is effective against plant pathogens. The results of this study indicated synthesized Ponericin W1 inhibited mycelial growth of Magnaporthe oryzae and Botrytis cinerea, as well as hyphal growth and spore production of Fusarium graminearum. Besides, Ponericin W1 exhibited antibacterial activities against Pseudomonas syringae pv. tomato and Xanthomonas oryzae pv. oryzae. After codon optimization, Ponericin W1 gene was constructed into plant expression vector, and transformed into Arabidopsis thaliana by floral dip method. The Ponericin W1 was located in intercellular space of the transgenic plants as expected. Compared with the wild-type plants, there were ungerminated spores and less hyphal, conidia on the leaves of transgenic plants after innoculation with the powdery mildew fungus Golovinomyces cichoracearum. After innoculation with the pathogenic bac-terium Pseudomonas syringae pv. tomato, the baceria in the leaves of transgenic plants was significantly less than the wild-type plants, indicating that the transgenic plants displayed enhanced disease resistance to pathogens. These results demonstrate a potential use of Ponericin W1 in genetic engineering for broad-spectrum plant disease resistance.

Genetic and Nongenetic Determinants of Cell Growth Variation Assessed by High-Throughput Microscopy

In microbial populations, growth initiation and proliferation rates are major components of fitness and therefore likely targets of selection. We used a high-throughput microscopy assay, which enables simultaneous analysis of tens of thousands of microcolonies, to determine the sources and extent of growth rate variation in the budding yeast (Saccharomyces cerevisiae) in different glucose environments. We find that cell growth rates are regulated by the extracellular concentration of glucose as proposed by Monod (1949), but that significant heterogeneity in growth rates is observed among genetically identical individuals within an environment. Yeast strains isolated from different geographic locations and habitats differ in their growth rate responses to different glucose concentrations. Inheritance patterns suggest that the genetic determinants of growth rates in different glucose concentrations are distinct. In addition, we identified genotypes that differ in the extent of variation in growth rate within an environment despite nearly identical mean growth rates, providing evidence that alleles controlling phenotypic variability segregate in yeast populations. We find that the time to reinitiation of growth (lag) is negatively correlated with growth rate, yet this relationship is strain-dependent. Between environments, the respirative activity of individual cells negatively correlates with glucose abundance and growth rate, but within an environment respirative activity and growth rate show a positive correlation, which we propose reflects differences in protein expression capacity. Our study quantifies the sources of genetic and nongenetic variation in cell growth rates in different glucose environments with unprecedented precision, facilitating their molecular genetic dissection.

Phylogenetic Portrait of the Saccharomyces cerevisiae Functional Genome

The genome of budding yeast (Saccharomyces cerevisiae) contains approximately 5800 protein-encoding genes, the majority of which are associated with some known biological function. Yet the extent of amino acid sequence conservation of these genes over all phyla has only been partially examined. Here we provide a more comprehensive overview and visualization of the conservation of yeast genes and a means for browsing and exploring the data in detail, down to the individual yeast gene, at http://yeast-phylogroups.princeton.edu. We used data from the OrthoMCL database, which has defined orthologs from approximately 150 completely sequenced genomes, including diverse representatives of the archeal, bacterial, and eukaryotic domains. By clustering genes based on similar patterns of conservation, we organized and visualized all the protein-encoding genes in yeast as a single heat map. Most genes fall into one of eight major clusters, called “phylogroups.” Gene ontology analysis of the phylogroups revealed that they were associated with specific, distinct trends in gene function, generalizations likely to be of interest to a wide range of biologists.

Cks1 Enhances Transcription Efficiency at the GAL1 Locus by Linking the Paf1 Complex to the 19S Proteasome

Cks1 was originally identified based on genetic interactions with CDC28, the gene that encodes Cdk1 in the budding yeast Saccharomyces cerevisiae. Subsequent work has shown that Cks1 binds Cdc28 and modulates its activity against certain substrates. However, the Cks1/Cdc28 complex also has a role in transcriptional chromatin remodeling not related to kinase activity. In order to elucidate protein networks associated with Cks1 transcriptional functions, proteomic analysis was performed on immunoaffinity-purified Cks1, identifying a physical interaction with the Paf1 complex. Specifically, we found that the Paf1 complex component Rtf1 interacts directly with Cks1 and that this interaction is essential for efficient recruitment of Cks1 to chromatin in the context of GAL1 gene induction. We further found that Cks1 in this capacity serves as an adaptor allowing Rtf1 to recruit 19S proteasome particles, shown to be required for efficient RNA production from some rapidly inducible genes such as GAL1.

Suppression of Chromosome Healing and Anticheckpoint Pathways in Yeast Postsenescence Survivors

Telomere repeat-like sequences at DNA double-strand breaks (DSBs) inhibit DNA damage signaling and serve as seeds to convert DSBs to new telomeres in mutagenic chromosome healing pathways. We find here that the response to seed-containing DSBs differs fundamentally between budding yeast (Saccharomyces cerevisiae) cells that maintain their telomeres via telomerase and so-called postsenescence survivors that use recombination-based alternative lengthening of telomere (ALT) mechanisms. Whereas telomere seeds are efficiently elongated by telomerase, they remain remarkably stable without de novo telomerization or extensive end resection in telomerase-deficient (est2Δ, tlc1Δ) postsenescence survivors. This telomere seed hyper-stability in ALT cells is associated with, but not caused by, prolonged DNA damage checkpoint activity (RAD9, RAD53) compared to telomerase-positive cells or presenescent telomerase-negative cells. The results indicate that both chromosome healing and anticheckpoint activity of telomere seeds are suppressed in yeast models of ALT pathways.

Pdc1 Functions in the Assembly of P Bodies in Schizosaccharomyces pombe

P bodies are cytoplasmic RNA granules containing the Dcp1-Dcp2 decapping enzymes where mRNA decay can occur. Here, we describe the characterization of P bodies in the fission yeast Schizosaccharomyces pombe. Most information on the property and function of P bodies stems from studies in the distantly related budding yeast Saccharomyces cerevisiae, and Edc3 was identified as a scaffold protein required for P-body assembly. However, we found that, unlike in S. cerevisiae, fission yeast Edc3 was dispensable for P-body formation. Pdc1, a novel partner of the fission yeast decapping enzyme, with a limited similarity to plant Edc4/Varicose that is required for the assembly of P bodies, was identified (tandem affinity purification-matrix-assisted laser desorption ionization tandem mass spectrometry [TAP-MALDI MS/MS]). Pdc1 interacts with Dcp2 through its C terminus and contains a coiled-coil region for self-interaction to mediate P-body formation. In line with the model that Pdc1 cross-bridges different proteins, additional interactions can be demonstrated with components such as Edc3 and Ste13. Although Pdc1 is not required for the interaction between Dcp1 and Dcp2, our data suggest that Pdc1 acts as a functional homologue of Edc4, a third component of the decapping enzymes that is thought to be absent from fungi. Together, these results highlight the diverse P-body protein compositions between different species and might help to provide insight into their evolutionary paths.

Mutational effects depend on ploidy level: all else is not equal

Ploidy is predicted to influence adaptation directly, yet whether single mutations behave the same in different ploidy backgrounds has not been well studied. It has often been assumed theoretically that aside from dominance, selective parameters do not differ between cells of varying ploidy. Using the budding yeast Saccharomyces cerevisiae, I compared the effect size of 20 adaptive mutations in haploids and homozygous diploids and found, surprisingly, that the same mutations often had a much larger effect in haploids than homozygous diploids. This empirical result demonstrates that it cannot be assumed that mutations will have the same effect in haploids and homozygous diploids.

Massive-Scale Gene Co-Expression Network Construction and Robustness Testing Using Random Matrix Theory

The study of gene relationships and their effect on biological function and phenotype is a focal point in systems biology. Gene co-expression networks built using microarray expression profiles are one technique for discovering and interpreting gene relationships. A knowledge-independent thresholding technique, such as Random Matrix Theory (RMT), is useful for identifying meaningful relationships. Highly connected genes in the thresholded network are then grouped into modules that provide insight into their collective functionality. While it has been shown that co-expression networks are biologically relevant, it has not been determined to what extent any given network is functionally robust given perturbations in the input sample set. For such a test, hundreds of networks are needed and hence a tool to rapidly construct these networks. To examine functional robustness of networks with varying input, we enhanced an existing RMT implementation for improved scalability and tested functional robustness of human (Homo sapiens), rice (Oryza sativa) and budding yeast (Saccharomyces cerevisiae). We demonstrate dramatic decrease in network construction time and computational requirements and show that despite some variation in global properties between networks, functional similarity remains high. Moreover, the biological function captured by co-expression networks thresholded by RMT is highly robust.