Thiamine-biosynthesis genes Bbpyr and Bbthi are required for conidial production and cell wall integrity of the entomopathogenic fungus Beauveria bassiana.

Beauveria bassiana is an important entomopathogenic fungus used to control a variety of insect pests. Conidia are the infective propagules of the fungus. However, some important factors that influence conidiation are still to be investigated. In this study, a mutant with decreased conidial production and hyphal growth was identified from a random T-DNA insertional library of B. bassiana. The corresponding gene (Bbthi) for this mutation encodes a putative thiazole synthase. Thiazole and pyrimidine are structural components of thiamine (vitamin B1), which is an essential nutrient for all forms of life. Disruption of Bbthi, Bbpyr, a putative pyrimidine synthetic gene, or both in B. bassiana results in a significant decrease of thiamine content. Loss of Bbthi and Bbpyr function significantly decreased the conidial production and hyphal growth, as well as disrupted the integrity of conidial cell wall. However, the defect of Bbpyr and Bbthi does not decrease the virulence of B. bassiana. Our results indicate the importance of thiamine biosynthesis in conidiation of B. bassiana, and provide useful information to produce conidia of entomopathogenic fungi for biocontrol of insect pests.

Development of a Mini-Replicon-Based Reverse-Genetics System for Rice Stripe Tenuivirus.

Negative-stranded RNA (NSR) viruses include both animal- and plant-infecting viruses that often cause serious diseases in humans and livestock and in agronomic crops. Rice stripe tenuivirus (RSV), a plant NSR virus with four negative-stranded/ambisense RNA segments, is one of the most destructive rice pathogens in many Asian countries. Due to the lack of a reliable reverse-genetics technology, molecular studies of RSV gene functions and its interaction with host plants are severely hampered. To overcome this obstacle, we developed a mini-replicon-based reverse-genetics system for RSV gene functional analysis in Nicotiana benthamiana. We first developed a mini-replicon system expressing an RSV genomic RNA3 enhanced green fluorescent protein (eGFP) reporter [MR3(-)eGFP], a nucleocapsid (NP), and a codon usage-optimized RNA-dependent RNA polymerase (RdRpopt). Using this mini-replicon system, we determined that RSV NP and RdRpopt are indispensable for the eGFP expression from MR3(-)eGFP. The expression of eGFP from MR3(-)eGFP can be significantly enhanced in the presence of four viral suppressors of RNA silencing (VSRs), NSs, and P19-HcPro-γb. In addition, NSvc4, the movement protein of RSV, facilitated eGFP trafficking between cells. We also developed an antigenomic RNA3-based replicon in N. benthamiana. However, we found that the RSV NS3 coding sequence acts as a cis element to regulate viral RNA expression. Finally, we made mini-replicons representing all four RSV genomic RNAs. This is the first mini-replicon-based reverse-genetics system for monocot-infecting tenuivirus. We believe that the mini-replicon system described here will allow studies of the RSV replication, transcription, cell-to-cell movement, and host machinery underpinning RSV infection in plants. IMPORTANCE Plant-infecting segmented negative-stranded RNA (NSR) viruses are grouped into three genera: Orthotospovirus, Tenuivirus, and Emaravirus. Reverse-genetics systems have been established for members of the genera Orthotospovirus and Emaravirus. However, there is still no reverse-genetics system available for Tenuivirus. Rice stripe virus (RSV) is a monocot-infecting tenuivirus with four negative-stranded/ambisense RNA segments. It is one of the most destructive rice pathogens and causes significant damage to the rice industry in Asian countries. Due to the lack of a reliable reverse-genetics system, molecular characterizations of RSV gene functions and the host machinery underpinning RSV infection in plants are extremely difficult. To overcome this obstacle, we developed a mini-replicon-based reverse-genetics system for RSV in Nicotiana benthamiana. This is the first mini-replicon-based reverse-genetics system for tenuivirus. We consider that this system will provide researchers a new working platform to elucidate the molecular mechanisms dictating segmented tenuivirus infections in plants.

Signatures of optimal codon usage in metabolic genes inform budding yeast ecology.

Reverse ecology is the inference of ecological information from patterns of genomic variation. One rich, heretofore underutilized, source of ecologically relevant genomic information is codon optimality or adaptation. Bias toward codons that match the tRNA pool is robustly associated with high gene expression in diverse organisms, suggesting that codon optimization could be used in a reverse ecology framework to identify highly expressed, ecologically relevant genes. To test this hypothesis, we examined the relationship between optimal codon usage in the classic galactose metabolism (GAL) pathway and known ecological niches for 329 species of budding yeasts, a diverse subphylum of fungi. We find that optimal codon usage in the GAL pathway is positively correlated with quantitative growth on galactose, suggesting that GAL codon optimization reflects increased capacity to grow on galactose. Optimal codon usage in the GAL pathway is also positively correlated with human-associated ecological niches in yeasts of the CUG-Ser1 clade and with dairy-associated ecological niches in the family Saccharomycetaceae. For example, optimal codon usage of GAL genes is greater than 85% of all genes in the genome of the major human pathogen Candida albicans (CUG-Ser1 clade) and greater than 75% of genes in the genome of the dairy yeast Kluyveromyces lactis (family Saccharomycetaceae). We further find a correlation between optimization in the GALactose pathway genes and several genes associated with nutrient sensing and metabolism. This work suggests that codon optimization harbors information about the metabolic ecology of microbial eukaryotes. This information may be particularly useful for studying fungal dark matter-species that have yet to be cultured in the lab or have only been identified by genomic material.

Quantifying absolute gene expression profiles reveals distinct regulation of central carbon metabolism genes in yeast.

In addition to controlled expression of genes by specific regulatory circuits, the abundance of proteins and transcripts can also be influenced by physiological states of the cell such as growth rate and metabolism. Here we examine the control of gene expression by growth rate and metabolism, by analyzing a multi-omics dataset consisting of absolute-quantitative abundances of the transcriptome, proteome, and amino acids in 22 steady-state yeast cultures. We find that transcription and translation are coordinately controlled by the cell growth rate via RNA polymerase II and ribosome abundance, but they are independently controlled by nitrogen metabolism via amino acid and nucleotide availabilities. Genes in central carbon metabolism, however, are distinctly regulated and do not respond to the cell growth rate or nitrogen metabolism as all other genes. Understanding these effects allows the confounding factors of growth rate and metabolism to be accounted for in gene expression profiling studies.

Cell cycle-dependent and independent mating blocks ensure fungal zygote survival and ploidy maintenance.

To ensure genome stability, sexually reproducing organisms require that mating brings together exactly 2 haploid gametes and that meiosis occurs only in diploid zygotes. In the fission yeast Schizosaccharomyces pombe, fertilization triggers the Mei3-Pat1-Mei2 signaling cascade, which represses subsequent mating and initiates meiosis. Here, we establish a degron system to specifically degrade proteins postfusion and demonstrate that mating blocks not only safeguard zygote ploidy but also prevent lysis caused by aberrant fusion attempts. Using long-term imaging and flow-cytometry approaches, we identify previously unrecognized and independent roles for Mei3 and Mei2 in zygotes. We show that Mei3 promotes premeiotic S-phase independently of Mei2 and that cell cycle progression is both necessary and sufficient to reduce zygotic mating behaviors. Mei2 not only imposes the meiotic program and promotes the meiotic cycle, but also blocks mating behaviors independently of Mei3 and cell cycle progression. Thus, we find that fungi preserve zygote ploidy and survival by at least 2 mechanisms where the zygotic fate imposed by Mei2 and the cell cycle reentry triggered by Mei3 synergize to prevent zygotic mating.

PRD-2 directly regulates casein kinase I and counteracts nonsense-mediated decay in the Neurospora circadian clock.

Circadian clocks in fungi and animals are driven by a functionally conserved transcription-translation feedback loop. In Neurospora crassa, negative feedback is executed by a complex of Frequency (FRQ), FRQ-interacting RNA helicase (FRH), and casein kinase I (CKI), which inhibits the activity of the clock's positive arm, the White Collar Complex (WCC). Here, we show that the prd-2 (period-2) gene, whose mutation is characterized by recessive inheritance of a long 26 hr period phenotype, encodes an RNA-binding protein that stabilizes the ck-1a transcript, resulting in CKI protein levels sufficient for normal rhythmicity. Moreover, by examining the molecular basis for the short circadian period of upf-1prd-6 mutants, we uncovered a strong influence of the Nonsense Mediated Decay pathway on CKI levels. The finding that circadian period defects in two classically derived Neurospora clock mutants each arise from disruption of ck-1a regulation is consistent with circadian period being exquisitely sensitive to levels of casein kinase I.

The putative transient receptor potential channel protein encoded by the orf19.4805 gene is involved in cation sensitivity, antifungal tolerance, and filamentation in Candida albicans.

Transient receptor potential (TRP) channels, an ancient family of cation channels, are highly conserved in eukaryotes and play various physiological functions, ranging from sensation of ion homeostasis to reception of pain and vision. Calcium-permeable TRP channels have been identified from the plant Arabidopsis thaliana (AtCsc1) and the budding yeast Saccharomyces cerevisiae (ScCsc1). In this study, we characterized the functions of the Csc1 homolog, orf19.4805, in Candida albicans. Orf19.4805 is a protein of 866 amino acids and 11 transmembrane domains, which shares 49% identity (69% similarity) in amino acid sequence with ScRsn1. Here, we demonstrate that deletion of the orf19.4805 gene causes C. albicans cells to be sensitive to SDS (sodium dodecyl sulfate) and antifungal drugs, and tolerance to zinc, manganese, and cadmium ions. Candida albicans cells lacking orf19.4805 show a defect in filamentation in vitro. Therefore, orf19.4805 is involved in the regulation of cation homeostasis and filamentation in C. albicans.

Genome-Wide Detection of Genes Under Positive Selection in Worldwide Populations of the Barley Scald Pathogen.

Coevolution between hosts and pathogens generates strong selection pressures to maintain resistance and infectivity, respectively. Genomes of plant pathogens often encode major effect loci for the ability to successfully infect specific host genotypes. Hence, spatial heterogeneity in host genotypes coupled with abiotic factors could lead to locally adapted pathogen populations. However, the genetic basis of local adaptation is poorly understood. Rhynchosporium commune, the pathogen causing barley scald disease, interacts at least partially in a gene-for-gene manner with its host. We analyzed global field populations of 125 R. commune isolates to identify candidate genes for local adaptation. Whole genome sequencing data showed that the pathogen is subdivided into three genetic clusters associated with distinct geographic and climatic regions. Using haplotype-based selection scans applied independently to each genetic cluster, we found strong evidence for selective sweeps throughout the genome. Comparisons of loci under selection among clusters revealed little overlap, suggesting that ecological differences associated with each cluster led to variable selection regimes. The strongest signals of selection were found predominantly in the two clusters composed of isolates from Central Europe and Ethiopia. The strongest selective sweep regions encoded protein functions related to biotic and abiotic stress responses. Selective sweep regions were enriched in genes encoding functions in cellular localization, protein transport activity, and DNA damage responses. In contrast to the prevailing view that a small number of gene-for-gene interactions govern plant pathogen evolution, our analyses suggest that the evolutionary trajectory is largely determined by spatially heterogeneous biotic and abiotic selection pressures.

Roles of phospholipid methyltransferases in pycnidia development, stress tolerance and secondary metabolism in the taxol-producing fungus Pestalotiopsis microspore.

Phosphatidylcholine (PC) is an important membrane component of the eukaryotic cell. In yeast fungi, two phospholipid methyltransferases catalyze consecutive steps of methylation in the formation of phosphatidylcholine from phosphatidylethanolamine. However, roles of phospholipid methyltransferases in filamentous fungi remains less investigated. We report here the characterization of two genes, choA and choC, that putatively encoded phospholipid methyltransferases in the taxol-producing fungus Pestalotiopsis microspora. Deletion of choC resulted in defects in PC production, vegetative growth and development of asexual structure. The mutant strains exhibited multiple morphological abnormalities, e.g. swollen hyphal tips and enhanced hyphal branching, and even mycelial autolysis. Some novel roles for the genes were also revealed, for instance, the deletion of either choC or choA impaired the development of pycnidia and conidia, the cell wall integrity. The mutant strains displayed a hypersensitivity to stress conditions, e.g. osmotic stress, cold and metal ions. The osmotic hypersensitivity indicates a crosstalk of PC pathways to other signaling pathways, such as the HOG pathway. Still more, choA, but not choC, was required for the production of secondary metabolites, e.g. pestalotiollide B, suggesting distinct roles of the two genes. This work would contribute to better understanding the function of phospholipid methyltransferases in fungi.

Multispecies reconstructions uncover widespread conservation, and lineage-specific elaborations in eukaryotic mRNA metabolism.

The degree of conservation and evolution of cytoplasmic mRNA metabolism pathways across the eukaryotes remains incompletely resolved. In this study, we describe a comprehensive genome and transcriptome-wide analysis of proteins involved in mRNA maturation, translation, and mRNA decay across representative organisms from the six eukaryotic super-groups. We demonstrate that eukaryotes share common pathways for mRNA metabolism that were almost certainly present in the last eukaryotic common ancestor, and show for the first time a correlation between intron density and a selective absence of some Exon Junction Complex (EJC) components in eukaryotes. In addition, we identify pathways that have diversified in individual lineages, with a specific focus on the unique gene gains and losses in members of the Excavata and SAR groups that contribute to their unique gene expression pathways compared to other organisms.

Relative Abundances of Candida albicans and Candida glabrata in In Vitro Coculture Biofilms Impact Biofilm Structure and Formation.

Candida is a member of the normal human microbiota and often resides on mucosal surfaces such as the oral cavity or the gastrointestinal tract. In addition to their commensality, Candida species can opportunistically become pathogenic if the host microbiota is disrupted or if the host immune system becomes compromised. An important factor for Candida pathogenesis is its ability to form biofilm communities. The two most medically important species-Candida albicans and Candida glabrata-are often coisolated from infection sites, suggesting the importance of Candida coculture biofilms. In this work, we report that biofilm formation of the coculture population depends on the relative ratio of starting cell concentrations of C. albicans and C. glabrata When using a starting ratio of C. albicans to C. glabrata of 1:3, ∼6.5- and ∼2.5-fold increases in biofilm biomass were observed relative to those of a C. albicans monoculture and a C. albicans/C. glabrata ratio of 1:1, respectively. Confocal microscopy analysis revealed the heterogeneity and complex structures composed of long C. albicans hyphae and C. glabrata cell clusters in the coculture biofilms, and reverse transcription-quantitative PCR (qRT-PCR) studies showed increases in the relative expression of the HWP1 and ALS3 adhesion genes in the C. albicans/C. glabrata 1:3 biofilm compared to that in the C. albicans monoculture biofilm. Additionally, only the 1:3 C. albicans/C. glabrata biofilm demonstrated an increased resistance to the antifungal drug caspofungin. Overall, the results suggest that interspecific interactions between these two fungal pathogens increase biofilm formation and virulence-related gene expression in a coculture composition-dependent manner.IMPORTANCECandida albicans and Candida glabrata are often coisolated during infection, and the occurrence of coisolation increases with increasing inflammation, suggesting possible synergistic interactions between the two Candida species in pathogenesis. During the course of an infection, the prevalence of each Candida species may change over time due to differences in metabolism and in the resistance of each species to antifungal therapies. Therefore, it is necessary to understand the dynamics between C. albicans and C. glabrata in coculture to develop better therapeutic strategies against Candida infections. Existing in vitro work has focused on understanding how an equal-part culture of C. albicans and C. glabrata impacts biofilm formation and pathogenesis. What is not understood, and what is investigated in this work, is how the composition of Candida species in coculture impacts overall biofilm formation, virulence gene expression, and the therapeutic treatment of biofilms.

Cross-pathway control gene CPC1/GCN4 coordinates with histone acetyltransferase GCN5 to regulate catalase-3 expression under oxidative stress in Neurospora crassa.

Catalase is an important enzyme found in nearly all aerobic organisms and plays an essential role in protecting cells from oxidative damage by catalyzing the degradation of hydrogen peroxide into water and oxygen. In filamentous fungus Neurospora crassa, the expression levels of catalases are rigorously regulated by morphogenetic transition during growth and development in cells. Our study revealed that catalase-3 transcription is positively regulated by histone acetyltransferase GCN5 and the cross-pathway control gene cpc-1, as the cat-3 expression level is significantly decreased in gcn5KO and cpc-1 (j-5) mutants. Moreover, gcn5KO and cpc-1 (j-5) mutants could not respond to H2O2 treatment due to the inadequate cat-3 transcription, while wild-type strains showed high expression levels of catalase upon H2O2 treatment. The global H3 acetylation and the acetylation of H3 at cat-3 locus dramatically decreased in gcn5KO under normal or oxidative stress conditions. Meanwhile, the expression of CAT-3 is reduced in gcn5E146Q, the catalytically dead mutant, suggesting that the catalytic activity of GCN5 functions in regulation of cat-3 transcription. In addition, GCN5 cannot acetylate histone H3 efficiently at cat-3 locus in cpc-1 (j-5) mutant strains under normal or oxidative stress conditions. Furthermore, ChIP assays data revealed that the CPC1/GCN4 can directly target the cat-3 promoter region, which may recruit GCN5 to modify the histone acetylation of this region. These results disclosed a distinctive function of CPC1/GCN4 in the regulatory pathway of cat-3 transcription, which is mediated by GCN5-dependent acetylation.

Comparative transcriptome and gene co-expression network analysis reveal genes and signaling pathways adaptively responsive to varied adverse stresses in the insect fungal pathogen, Beauveria bassiana.

Sensing, responding, and adapting to the surrounding environment are crucial for all living organisms to survive, proliferate, and differentiate in their biological niches. Beauveria bassiana is an economically important insect-pathogenic fungus which is widely used as a biocontrol agent to control a variety of insect pests. The fungal pathogen unavoidably encounters a variety of adverse environmental stresses and defense response from the host insects during application of the fungal agents. However, few are known about the transcription response of the fungus to respond or adapt varied adverse stresses. Here, we comparatively analyzed the transcriptome of B. bassiana in globe genome under the varied stationary-phase stresses including osmotic agent (0.8 M NaCl), high temperature (32 °C), cell wall-perturbing agent (Congo red), and oxidative agents (H2O2 or menadione). Total of 12,412 reads were obtained, and mapped to the 6767 genes of the B. bassiana. All of these stresses caused transcription responses involved in basal metabolism, cell wall construction, stress response or cell rescue/detoxification, signaling transduction and gene transcription regulation, and likely other cellular processes. An array of genes displayed similar transcription patterns in response to at least two of the five stresses, suggesting a shared transcription response to varied adverse stresses. Gene co-expression network analysis revealed that mTOR signaling pathway, but not HOG1 MAP kinase pathway, played a central role in regulation the varied adverse stress responses, which was verified by RNAi-mediated knockdown of TOR1. Our findings provided an insight of transcription response and gene co-expression network of B. bassiana in adaptation to varied environments.

Discovering the role of the apolipoprotein gene and the genes in the putative pullulan biosynthesis pathway on the synthesis of pullulan, heavy oil and melanin in Aureobasidium pullulans.

Pullulan produced by Aureobasidium pullulans presents various applications in food manufacturing and pharmaceutical industry. However, the pullulan biosynthesis mechanism remains unclear. This work proposed a pathway suggesting that heavy oil and melanin may correlate with pullulan production. The effects of overexpression or deletion of genes encoding apolipoprotein, UDPG-pyrophosphorylase, glucosyltransferase, and α-phosphoglucose mutase on the production of pullulan, heavy oil, and melanin were examined. Pullulan production increased by 16.93 and 8.52% with the overexpression of UDPG-pyrophosphorylase and apolipoprotein genes, respectively. Nevertheless, the overexpression or deletion of other genes exerted little effect on pullulan biosynthesis. Heavy oil production increased by 146.30, 64.81, and 33.33% with the overexpression of UDPG-pyrophosphorylase, α-phosphoglucose mutase, and apolipoprotein genes, respectively. Furthermore, the syntheses of pullulan, heavy oil, and melanin can compete with one another. This work may provide new guidance to improve the production of pullulan, heavy oil, and melanin through genetic approach.

The functions of PEX genes in peroxisome biogenesis and pathogenicity in phytopathogenic fungi.

Peroxisomes are cellular organelles present ubiquitously in eukaryotic cells and are involved in ß-oxidation, glyoxylate cycle and a variety of biochemical metabolisms. Recently peroxisomes have been demonstrated to play vital roles in the host infection processes by plant fungal pathogens. The biogenesis of peroxisomes requires a category of proteins named peroxins, which are encoded by the PEX genes. So far, more than 10 PEX genes were isolated in phytopathogenic fungi, and significant research efforts are focused on the mechanism of peroxisome formation and the roles of peroxisome in the development and pathogenicity of fungal pathogens. In this review, we summarize the latest advances in peroxisome biogenesis and functions in pathogenic fungi, including the roles of PEXs in life cycle of peroxisome, peroxisome related metabolisms, and fungal development, infection and pathogenicity, in order to provide references for future studies in plant pathogenic fungi and the control of disease.

Multiple Taf subunits of TFIID interact with Ino2 activation domains and contribute to expression of genes required for yeast phospholipid biosynthesis.

Expression of phospholipid biosynthetic genes in yeast requires activator protein Ino2 which can bind to the UAS element inositol/choline-responsive element (ICRE) and trigger activation of target genes, using two separate transcriptional activation domains, TAD1 and TAD2. However, it is still unknown which cofactors mediate activation by TADs of Ino2. Here, we show that multiple subunits of basal transcription factor TFIID (TBP-associated factors Taf1, Taf4, Taf6, Taf10 and Taf12) are able to interact in vitro with activation domains of Ino2. Interaction was no longer observed with activation-defective variants of TAD1. We were able to identify two nonoverlapping regions in the N-terminus of Taf1 (aa 1-100 and aa 182-250) each of which could interact with TAD1 of Ino2 as well as with TAD4 of activator Adr1. Specific missense mutations within Taf1 domain aa 182-250 affecting basic and hydrophobic residues prevented interaction with wild-type TAD1 and caused reduced expression of INO1. Using chromatin immunoprecipitation we demonstrated Ino2-dependent recruitment of Taf1 and Taf6 to ICRE-containing promoters INO1 and CHO2. Transcriptional derepression of INO1 was no longer possible with temperature-sensitive taf1 and taf6 mutants cultivated under nonpermissive conditions. This result supports the hypothesis of Taf-dependent expression of structural genes activated by Ino2.

A tree of life based on ninety-eight expressed genes conserved across diverse eukaryotic species.

Rapid advances in DNA sequencing technologies have resulted in the accumulation of large data sets in the public domain, facilitating comparative studies to provide novel insights into the evolution of life. Phylogenetic studies across the eukaryotic taxa have been reported but on the basis of a limited number of genes. Here we present a genome-wide analysis across different plant, fungal, protist, and animal species, with reference to the 36,002 expressed genes of the rice genome. Our analysis revealed 9831 genes unique to rice and 98 genes conserved across all 49 eukaryotic species analysed. The 98 genes conserved across diverse eukaryotes mostly exhibited binding and catalytic activities and shared common sequence motifs; and hence appeared to have a common origin. The 98 conserved genes belonged to 22 functional gene families including 26S protease, actin, ADP-ribosylation factor, ATP synthase, casein kinase, DEAD-box protein, DnaK, elongation factor 2, glyceraldehyde 3-phosphate, phosphatase 2A, ras-related protein, Ser/Thr protein phosphatase family protein, tubulin, ubiquitin and others. The consensus Bayesian eukaryotic tree of life developed in this study demonstrated widely separated clades of plants, fungi, and animals. Musa acuminata provided an evolutionary link between monocotyledons and dicotyledons, and Salpingoeca rosetta provided an evolutionary link between fungi and animals, which indicating that protozoan species are close relatives of fungi and animals. The divergence times for 1176 species pairs were estimated accurately by integrating fossil information with synonymous substitution rates in the comprehensive set of 98 genes. The present study provides valuable insight into the evolution of eukaryotes.

Funneled potential and flux landscapes dictate the stabilities of both the states and the flow: Fission yeast cell cycle.

Using fission yeast cell cycle as an example, we uncovered that the non-equilibrium network dynamics and global properties are determined by two essential features: the potential landscape and the flux landscape. These two landscapes can be quantified through the decomposition of the dynamics into the detailed balance preserving part and detailed balance breaking non-equilibrium part. While the funneled potential landscape is often crucial for the stability of the single attractor networks, we have uncovered that the funneled flux landscape is crucial for the emergence and maintenance of the stable limit cycle oscillation flow. This provides a new interpretation of the origin for the limit cycle oscillations: There are many cycles and loops existed flowing through the state space and forming the flux landscapes, each cycle with a probability flux going through the loop. The limit cycle emerges when a loop stands out and carries significantly more probability flux than other loops. We explore how robustness ratio (RR) as the gap or steepness versus averaged variations or roughness of the landscape, quantifying the degrees of the funneling of the underlying potential and flux landscapes. We state that these two landscapes complement each other with one crucial for stabilities of states on the cycle and the other crucial for the stability of the flow along the cycle. The flux is directly related to the speed of the cell cycle. This allows us to identify the key factors and structure elements of the networks in determining the stability, speed and robustness of the fission yeast cell cycle oscillations. We see that the non-equilibriumness characterized by the degree of detailed balance breaking from the energy pump quantified by the flux is the cause of the energy dissipation for initiating and sustaining the replications essential for the origin and evolution of life. Regulating the cell cycle speed is crucial for designing the prevention and curing strategy of cancer.

Organization and evolution of mating-type genes in three Stagonosporopsis species causing gummy stem blight of cucurbits and leaf spot and dry rot of papaya.

Population divergence and speciation of closely related lineages can result from reproductive differences leading to genetic isolation. An increasing number of fungal diseases of plants and animals have been determined to be caused by morphologically indistinguishable species that are genetically distinct, thereby representing cryptic species. We were interested in identifying if mating systems among three Stagonosporopsis species (S. citrulli, S. cucurbitacearum, and S. caricae) causing gummy stem blight (GSB) of cucurbits or leaf spot and dry rot of papaya differed, possibly underlying species divergence. Additionally, we were interested in identifying evolutionary pressures acting on the genes controlling mating in these fungi. The mating-type loci (MAT1) of three isolates from each of the three species were identified in draft genome sequences. For the three species, MAT1 was structurally identical and contained both mating-type genes necessary for sexual reproduction, which suggests that all three species are homothallic. However, both MAT1-1-1 and MAT1-2-1 were divergent among species showing rapid evolution with a much greater number of amino acid-changing substitutions detected for the reproductive genes compared with genes flanking MAT1. Positive selection was detected in MAT1-2-1, especially in the highly conserved high mobility group (MATA_HMG-box) domain. Thus, the mating-type genes are rapidly evolving in GSB fungi, but a difference in mating systems among the three species does not underlie their divergence.

Finding optimum width of discretization for gene expressions using functional annotations.

Discretizing gene expression values is an important step in data preprocessing as it helps in reducing noise and experimental errors. This in turn provides better results in various tasks such as gene regulatory network analysis and disease prediction. A supervised discretization method for gene expressions using gene annotation is developed. The method is called "Gene Annotation Based Discretization" (GABD) where the discretization width is determined by maximizing the positive predictive value (PPV), computed using gene annotations, for top 20,000 gene pairs. The method can capture the gene similarity better than those obtained using original expressions. The performance of GABD is compared with some existing discretization methods like equal width discretization, equal frequency discretization and k-means discretization in terms of positive predictive value (PPV). The utility of GABD is also shown by clustering genes using k-medoid algorithm and thereby predicting the function of 23 unclassified Saccharomyces cerevisiae genes using p-value cut off 10-10. The source code for GABD is available at http://www.sampa.droppages.com/GABD.html.