Extending the reach of homology by using successive computational filters to find yeast pheromone genes

The mating of fungi depends on pheromones that mediate communication between two mating types. Most species use short peptides as pheromones, which are either unmodified (e.g., α-factor in Saccharomyces cerevisiae) or C-terminally farnesylated (e.g., a-factor in S. cerevisiae). Peptide pheromones have been found by genetics or biochemistry in a small number of fungi, but their short sequences and modest conservation make it impossible to detect homologous sequences in most species. To overcome this problem, we used a four-step computational pipeline to identify candidate a-factor genes in sequenced genomes of the Saccharomycotina, the fungal clade that contains most of the yeasts: we require that candidate genes have a C-terminal prenylation motif, are shorter than 100 amino acids long, and contain a proteolytic-processing motif upstream of the potential mature pheromone sequence and that closely related species contain highly conserved homologs of the potential mature pheromone sequence. Additional manual curation exploits the observation that many species carry more than one a-factor gene, encoding identical or nearly identical pheromones. From 332 Saccharomycotina genomes, we identified strong candidate pheromone genes in 241 genomes, covering 13 clades that are each separated from each other by at least 100 million years, the time required for evolution to remove detectable sequence homology among small pheromone genes. For one small clade, the Yarrowia, we demonstrated that our algorithm found the a-factor genes: deleting all four related genes in the a-mating type of Yarrowia lipolytica prevents mating.

Population Genomics of the Maize Pathogen Ustilago maydis: Demographic History and Role of Virulence Clusters in Adaptation

The tight interaction between pathogens and their hosts results in reciprocal selective forces that impact the genetic diversity of the interacting species. The footprints of this selection differ between pathosystems because of distinct life-history traits, demographic histories, or genome architectures. Here, we studied the genome-wide patterns of genetic diversity of 22 isolates of the causative agent of the corn smut disease, Ustilago maydis, originating from five locations in Mexico, the presumed center of origin of this species. In this species, many genes encoding secreted effector proteins reside in so-called virulence clusters in the genome, an arrangement that is so far not found in other filamentous plant pathogens. Using a combination of population genomic statistical analyses, we assessed the geographical, historical, and genome-wide variation of genetic diversity in this fungal pathogen. We report evidence of two partially admixed subpopulations that are only loosely associated with geographic origin. Using the multiple sequentially Markov coalescent model, we inferred the demographic history of the two pathogen subpopulations over the last 0.5 Myr. We show that both populations experienced a recent strong bottleneck starting around 10,000 years ago, coinciding with the assumed time of maize domestication. Although the genome average genetic diversity is low compared with other fungal pathogens, we estimated that the rate of nonsynonymous adaptive substitutions is three times higher in genes located within virulence clusters compared with nonclustered genes, including nonclustered effector genes. These results highlight the role that these singular genomic regions play in the evolution of this pathogen.

Design of an improved universal signal peptide based on the α-factor mating secretion signal for enzyme production in yeast

Saccharomyces cerevisiae plays an important role in the heterologous expression of an array of proteins due to its easy manipulation, low requirements and ability for protein post-translational modifications. The implementation of the preproleader secretion signal of the α-factor mating pheromone from this yeast contributes to increase the production yields by targeting the foreign protein to the extracellular environment. The use of this signal peptide combined with enzyme-directed evolution allowed us to achieve the otherwise difficult functional expression of fungal laccases in S. cerevisiae, obtaining different evolved α-factor preproleader sequences that enhance laccase secretion. However, the design of a universal signal peptide to enhance the production of heterologous proteins in S. cerevisiae is a pending challenge. We describe here the optimisation of the α-factor preproleader to improve recombinant enzyme production in S. cerevisiae through two parallel engineering strategies: a bottom-up design over the native α-factor preproleader (αnat) and a top-down design over the fittest evolved signal peptide obtained in our lab (α9H2 leader). The goal was to analyse the effect of mutations accumulated in the signal sequence throughout iterations of directed evolution, or of other reported mutations, and their possible epistatic interactions. Both approaches agreed in the positive synergism of four mutations (Aα9D, Aα20T, Lα42S, Dα83E) contained in the final optimised leader (αOPT), which notably enhanced the secretion of several fungal oxidoreductases and hydrolases. Additionally, we suggest a guideline to further drive the heterologous production of a particular enzyme based on combinatorial saturation mutagenesis of positions 86th and 87th of the αOPT leader fused to the target protein.

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.

During sexual reproduction, fertilization must happen between exactly two gametes to ensure genome stability. This study shows that two mechanisms – establishment of zygotic fate and re-entry to the cell cycle – combine to prevent fission yeast zygotes fusing with further gametes.

Cytoplasmic retention and degradation of a mitotic inducer enable plant infection by a pathogenic fungus

In the fungus Ustilago maydis, sexual pheromones elicit mating resulting in an infective filament able to infect corn plants. Along this process a G2 cell cycle arrest is mandatory. Such as cell cycle arrest is initiated upon the pheromone recognition in each mating partner, and sustained once cell fusion occurred until the fungus enter the plant tissue. We describe that the initial cell cycle arrest resulted from inhibition of the nuclear transport of the mitotic inducer Cdc25 by targeting its importin, Kap123. Near cell fusion to take place, the increase on pheromone signaling promotes Cdc25 degradation, which seems to be important to ensure the maintenance of the G2 cell cycle arrest to lead the formation of the infective filament. This way, premating cell cycle arrest is linked to the subsequent steps required for establishment of the infection. Disabling this connection resulted in the inability of fungal cells to infect plants.

Real-Time Genetic Compensation Defines the Dynamic Demands of Feedback Control

Biological signaling networks use feedback control to dynamically adjust their operation in real time. Traditional static genetic methods such as gene knockouts or rescue experiments can often identify the existence of feedback interactions but are unable to determine what feedback dynamics are required. Here, we implement a new strategy, closed-loop optogenetic compensation (CLOC), to address this problem. Using a custom-built hardware and software infrastructure, CLOC monitors, in real time, the output of a pathway deleted for a feedback regulator. A minimal model uses these measurements to calculate and deliver-on the fly-an optogenetically enabled transcriptional input designed to compensate for the effects of the feedback deletion. Application of CLOC to the yeast pheromone response pathway revealed surprisingly distinct dynamic requirements for three well-studied feedback regulators. CLOC, a marriage of control theory and traditional genetics, presents a broadly applicable methodology for defining the dynamic function of biological feedback regulators.

Identification and characterization of the Komagataella phaffii mating pheromone genes

The methylotrophic yeast Komagataella phaffii (Pichia pastoris) is a haploid yeast that is able to form diploid cells by mating once nitrogen becomes limiting. Activation of the mating response requires the secretion of a- and α-factor pheromones, which bind to G-protein coupled receptors on cells of opposite mating type. In K. phaffii, the genes coding for the α-factor (MFα), the pheromone surface receptors and the conserved a-factor biogenesis pathway have been annotated previously. Initial homology-based search failed to identify potential a-factor genes (MFA). By using transcriptome data of heterothallic strains under mating conditions, we found two K. phaffiia-factor genes. Deletion of both MFA genes prevented mating of a-type cells. MFA single mutants were still able to mate and activate the mating response pathway in α-type cells. A reporter assay was used to confirm the biological activity of synthetic a- and α-factor peptides. The identification of the a-factor genes enabled the first characterization of the role and regulation of the mating pheromone genes and the response of K. phaffii to synthetic pheromones and will help to gain a better understanding of the mating behavior of K. phaffii.

Membrane curvature directs the localization of Cdc42p to novel foci required for cell-cell fusion

Cell fusion is ubiquitous in eukaryotic fertilization and development. The highly conserved Rho-GTPase Cdc42p promotes yeast fusion through interaction with Fus2p, a pheromone-induced amphiphysin-like protein. We show that in prezygotes, Cdc42p forms a novel Fus2p-dependent focus at the center of the zone of cell fusion (ZCF) and remains associated with remnant cell walls after initial fusion. At the ZCF and during fusion, Cdc42p and Fus2p colocalized. In contrast, in shmoos, both proteins were near the cortex but spatially separate. Cdc42p focus formation depends on ZCF membrane curvature: mutant analysis showed that Cdc42p localization is negatively affected by shmoo-like positive ZCF curvature, consistent with the flattening of the ZCF during fusion. BAR-domain proteins such as the fusion proteins Fus2p and Rvs161p are known to recognize membrane curvature. We find that mutations that disrupt binding of the Fus2p/Rvs161p heterodimer to membranes affect Cdc42p ZCF localization. We propose that Fus2p localizes Cdc42p to the flat ZCF to promote cell wall degradation.

Pheromone independent unisexual development in Cryptococcus neoformans

The fungus Cryptococcus neoformans can undergo a-α bisexual and unisexual reproduction. Completion of both sexual reproduction modes requires similar cellular differentiation processes and meiosis. Although bisexual reproduction generates equal number of a and α progeny and is far more efficient than unisexual reproduction under mating-inducing laboratory conditions, the α mating type dominates in nature. Population genetic studies suggest that unisexual reproduction by α isolates might have contributed to this sharply skewed distribution of the mating types. However, the predominance of the α mating type and the seemingly inefficient unisexual reproduction observed under laboratory conditions present a conundrum. Here, we discovered a previously unrecognized condition that promotes unisexual reproduction while suppressing bisexual reproduction. Pheromone is the principal stimulus for bisexual development in Cryptococcus. Interestingly, pheromone and other components of the pheromone pathway, including the key transcription factor Mat2, are not necessary but rather inhibitory for Cryptococcus to complete its unisexual cycle under this condition. The inactivation of the pheromone pathway promotes unisexual reproduction despite the essential role of this pathway in non-self-recognition during bisexual reproduction. Nonetheless, the requirement for the known filamentation regulator Znf2 and the expression of hyphal or basidium specific proteins remain the same for pheromone-dependent or independent sexual reproduction. Transcriptome analyses and an insertional mutagenesis screen in mat2Δ identified calcineurin being essential for this process. We further found that Znf2 and calcineurin work cooperatively in controlling unisexual development in this fungus. These findings indicate that Mat2 acts as a repressor of pheromone-independent unisexual development while serving as an activator for a-α bisexual development. The bi-functionality of Mat2 might have allowed it to act as a toggle switch for the mode of sexual development in this ubiquitous eukaryotic microbe.

Barcode Sequencing Screen Identifies SUB1 as a Regulator of Yeast Pheromone Inducible Genes

The yeast pheromone response pathway serves as a valuable model of eukaryotic mitogen-activated protein kinase (MAPK) pathways, and transcription of their downstream targets. Here, we describe application of a screening method combining two technologies: fluorescence-activated cell sorting (FACS), and barcode analysis by sequencing (Bar-Seq). Using this screening method, and pFUS1-GFP as a reporter for MAPK pathway activation, we readily identified mutants in known mating pathway components. In this study, we also include a comprehensive analysis of the FUS1 induction properties of known mating pathway mutants by flow cytometry, featuring single cell analysis of each mutant population. We also characterized a new source of false positives resulting from the design of this screen. Additionally, we identified a deletion mutant, sub1Δ, with increased basal expression of pFUS1-GFP. Here, in the first ChIP-Seq of Sub1, our data shows that Sub1 binds to the promoters of about half the genes in the genome (tripling the 991 loci previously reported), including the promoters of several pheromone-inducible genes, some of which show an increase upon pheromone induction. Here, we also present the first RNA-Seq of a sub1Δ mutant; the majority of genes have no change in RNA, but, of the small subset that do, most show decreased expression, consistent with biochemical studies implicating Sub1 as a positive transcriptional regulator. The RNA-Seq data also show that certain pheromone-inducible genes are induced less in the sub1Δ mutant relative to the wild type, supporting a role for Sub1 in regulation of mating pathway genes. The sub1Δ mutant has increased basal levels of a small subset of other genes besides FUS1, including IMD2 and FIG1, a gene encoding an integral membrane protein necessary for efficient mating.